Additive synthesis

Definitions

Harmonic additive synthesis is closely related to the concept of a Fourier series which is a way of expressing a periodic function as the sum of sinusoidal functions with frequencies equal to integer multiples of a common fundamental frequency. These sinusoids are called harmonics, overtones, or generally, partials. In general, a Fourier series contains an infinite number of sinusoidal components, with no upper limit to the frequency of the sinusoidal functions and includes a DC component (one with frequency of 0 Hz). Frequencies outside of the human audible range can be omitted in additive synthesis. As a result, only a finite number of sinusoidal terms with frequencies that lie within the audible range are modeled in additive synthesis.

A waveform or function is said to be periodic if

y ( t ) = y ( t + P )  

for all t and for some period P .

The Fourier series of a periodic function is mathematically expressed as:

y ( t ) = a 0 2 + ∑ k = 1 ∞ [ a k cos ⁡ ( 2 π k f 0 t ) − b k sin ⁡ ( 2 π k f 0 t ) ] = a 0 2 + ∑ k = 1 ∞ r k cos ⁡ ( 2 π k f 0 t + ϕ k )

where

Being inaudible, the DC component, a 0 / 2 , and all components with frequencies higher than some finite limit, K f 0 , are omitted in the following expressions of additive synthesis.

Harmonic form

The simplest harmonic additive synthesis can be mathematically expressed as:

where y ( t ) is the synthesis output, r k , k f 0 , and ϕ k are the amplitude, frequency, and the phase offset, respectively, of the k th harmonic partial of a total of K harmonic partials, and f 0 is the fundamental frequency of the waveform and the frequency of the musical note.

Time-dependent amplitudes

More generally, the amplitude of each harmonic can be prescribed as a function of time, r k ( t ) , in which case the synthesis output is

Each envelope r k ( t ) should vary slowly relative to the frequency spacing between adjacent sinusoids. The bandwidth of r k ( t ) should be significantly less than f 0 .

Inharmonic form

Additive synthesis can also produce inharmonic sounds (which are aperiodic waveforms) in which the individual overtones need not have frequencies that are integer multiples of some common fundamental frequency. While many conventional musical instruments have harmonic partials (e.g. an oboe), some have inharmonic partials (e.g. bells). Inharmonic additive synthesis can be described as

y ( t ) = ∑ k = 1 K r k ( t ) cos ⁡ ( 2 π f k t + ϕ k ) ,

where f k is the constant frequency of k th partial.

Time-dependent frequencies

In the general case, the instantaneous frequency of a sinusoid is the derivative (with respect to time) of the argument of the sine or cosine function. If this frequency is represented in hertz, rather than in angular frequency form, then this derivative is divided by 2 π . This is the case whether the partial is harmonic or inharmonic and whether its frequency is constant or time-varying.

In the most general form, the frequency of each non-harmonic partial is a non-negative function of time, f k ( t ) , yielding

Broader definitions

Additive synthesis more broadly may mean sound synthesis techniques that sum simple elements to create more complex timbres, even when the elements are not sine waves. For example, F. Richard Moore listed additive synthesis as one of the "four basic categories" of sound synthesis alongside subtractive synthesis, nonlinear synthesis, and physical modeling. In this broad sense, pipe organs, which also have pipes producing non-sinusoidal waveforms, can be considered as a variant form of additive synthesizers. Summation of principal components and Walsh functions have also been classified as additive synthesis.

Implementation methods

Modern-day implementations of additive synthesis are mainly digital. (See section Discrete-time equations for the underlying discrete-time theory)

Oscillator bank synthesis

Additive synthesis can be implemented using a bank of sinusoidal oscillators, one for each partial.

Wavetable synthesis

In the case of harmonic, quasi-periodic musical tones, wavetable synthesis can be as general as time-varying additive synthesis, but requires less computation during synthesis. As a result, an efficient implementation of time-varying additive synthesis of harmonic tones can be accomplished by use of wavetable synthesis.

Group additive synthesis

Group additive synthesis is a method to group partials into harmonic groups (having different fundamental frequencies) and synthesize each group separately with wavetable synthesis before mixing the results.

Inverse FFT synthesis

An inverse Fast Fourier transform can be used to efficiently synthesize frequencies that evenly divide the transform period or "frame". By careful consideration of the DFT frequency-domain representation it is also possible to efficiently synthesize sinusoids of arbitrary frequencies using a series of overlapping frames and the inverse Fast Fourier transform.

Additive analysis/resynthesis

It is possible to analyze the frequency components of a recorded sound giving a "sum of sinusoids" representation. This representation can be re-synthesized using additive synthesis. One method of decomposing a sound into time varying sinusoidal partials is short-time Fourier transform (STFT)-based McAulay-Quatieri Analysis.

By modifying the sum of sinusoids representation, timbral alterations can be made prior to resynthesis. For example, a harmonic sound could be restructured to sound inharmonic, and vice versa. Sound hybridisation or "morphing" has been implemented by additive resynthesis.

Additive analysis/resynthesis has been employed in a number of techniques including Sinusoidal Modelling, Spectral Modelling Synthesis (SMS), and the Reassigned Bandwidth-Enhanced Additive Sound Model. Software that implements additive analysis/resynthesis includes: SPEAR, LEMUR, LORIS, SMSTools, ARSS.

Products

New England Digital Synclavier had a resynthesis feature where samples could be analyzed and converted into ”timbre frames” which were part of its additive synthesis engine. Technos acxel, launched in 1987, utilized the additive analysis/resynthesis model, in an FFT implementation.

Also a vocal synthesizer, Vocaloid have been implemented on the basis of additive analysis/resynthesis: its spectral voice model called Excitation plus Resonances (EpR) model is extended based on Spectral Modeling Synthesis (SMS), and its diphone concatenative synthesis is processed using spectral peak processing (SPP) technique similar to modified phase-locked vocoder (an improved phase vocoder for formant processing). Using these techniques, spectral components (formants) consisting of purely harmonic partials can be appropriately transformed into desired form for sound modeling, and sequence of short samples (diphones or phonemes) constituting desired phrase, can be smoothly connected by interpolating matched partials and formant peaks, respectively, in the inserted transition region between different samples. (See also Dynamic timbres)

Musical instruments

Additive synthesis is used in electronic musical instruments. It is the principal sound generation technique used by Eminent organs.

Speech synthesis

In linguistics research, harmonic additive synthesis was used in 1950s to play back modified and synthetic speech spectrograms.

Later, in early 1980s, listening tests were carried out on synthetic speech stripped of acoustic cues to assess their significance. Time-varying formant frequencies and amplitudes derived by linear predictive coding were synthesized additively as pure tone whistles. This method is called sinewave synthesis. Also the composite sinusoidal modeling (CSM) used on a singing speech synthesis feature on Yamaha CX5M (1984), is known to use a similar approach which was independently developed during 1966–1979. These methods are characterized by extraction and recomposition of a set of significant spectral peaks corresponding to the several resonance modes occurred in the oral cavity and nasal cavity, in a viewpoint of acoustics. This principle was also utilized on a physical modeling synthesis method, called modal synthesis.

History

Harmonic analysis was discovered by Joseph Fourier, who published an extensive treatise of his research in the context of heat transfer in 1822. The theory found an early application in prediction of tides. Around 1876, Lord Kelvin constructed a mechanical tide predictor. It consisted of a harmonic analyzer and a harmonic synthesizer, as they were called already in the 19th century. The analysis of tide measurements was done using James Thomson's integrating machine. The resulting Fourier coefficients were input into the synthesizer, which then used a system of cords and pulleys to generate and sum harmonic sinusoidal partials for prediction of future tides. In 1910, a similar machine was built for the analysis of periodic waveforms of sound. The synthesizer drew a graph of the combination waveform, which was used chiefly for visual validation of the analysis.

Georg Ohm applied Fourier's theory to sound in 1843. The line of work was greatly advanced by Hermann von Helmholtz, who published his eight years worth of research in 1863. Helmholtz believed that the psychological perception of tone color is subject to learning, while hearing in the sensory sense is purely physiological. He supported the idea that perception of sound derives from signals from nerve cells of the basilar membrane and that the elastic appendages of these cells are sympathetically vibrated by pure sinusoidal tones of appropriate frequencies. Helmholtz agreed with the finding of Ernst Chladni from 1787 that certain sound sources have inharmonic vibration modes.

In Helmholtz's time, electronic amplification was unavailable. For synthesis of tones with harmonic partials, Helmholtz built an electrically excited array of tuning forks and acoustic resonance chambers that allowed adjustment of the amplitudes of the partials. Built at least as early as in 1862, these were in turn refined by Rudolph Koenig, who demonstrated his own setup in 1872. For harmonic synthesis, Koenig also built a large apparatus based on his wave siren. It was pneumatic and utilized cut-out tonewheels, and was criticized for low purity of its partial tones. Also tibia pipes of pipe organs have nearly sinusoidal waveforms and can be combined in the manner of additive synthesis.

In 1938, with significant new supporting evidence, it was reported on the pages of Popular Science Monthly that the human vocal cords function like a fire siren to produce a harmonic-rich tone, which is then filtered by the vocal tract to produce different vowel tones. By the time, the additive Hammond organ was already on market. Most early electronic organ makers thought it too expensive to manufacture the plurality of oscillators required by additive organs, and began instead to build subtractive ones. In a 1940 Institute of Radio Engineers meeting, the head field engineer of Hammond elaborated on the company's new Novachord as having a “subtractive system” in contrast to the original Hammond organ in which “the final tones were built up by combining sound waves”. Alan Douglas used the qualifiers additive and subtractive to describe different types of electronic organs in a 1948 paper presented to the Royal Musical Association. The contemporary wording additive synthesis and subtractive synthesis can be found in his 1957 book The electrical production of music, in which he categorically lists three methods of forming of musical tone-colours, in sections titled Additive synthesis, Subtractive synthesis, and Other forms of combinations.

A typical modern additive synthesizer produces its output as an electrical, analog signal, or as digital audio, such as in the case of software synthesizers, which became popular around year 2000.

Timeline

The following is a timeline of historically and technologically notable analog and digital synthesizers and devices implementing additive synthesis.

Discrete-time equations

In digital implementations of additive synthesis, discrete-time equations are used in place of the continuous-time synthesis equations. A notational convention for discrete-time signals uses brackets i.e. y [ n ] and the argument n can only be integer values. If the continuous-time synthesis output y ( t ) is expected to be sufficiently bandlimited; below half the sampling rate or f s / 2 , it suffices to directly sample the continuous-time expression to get the discrete synthesis equation. The continuous synthesis output can later be reconstructed from the samples using a digital-to-analog converter. The sampling period is T = 1 / f s .

Beginning with (3),

y ( t ) = ∑ k = 1 K r k ( t ) cos ⁡ ( 2 π ∫ 0 t f k ( u )   d u + ϕ k )

and sampling at discrete times t = n T = n / f s results in

y [ n ] = y ( n T ) = ∑ k = 1 K r k ( n T ) cos ⁡ ( 2 π ∫ 0 n T f k ( u )   d u + ϕ k ) = ∑ k = 1 K r k ( n T ) cos ⁡ ( 2 π ∑ i = 1 n ∫ ( i − 1 ) T i T f k ( u )   d u + ϕ k ) = ∑ k = 1 K r k ( n T ) cos ⁡ ( 2 π ∑ i = 1 n ( T f k [ i ] ) + ϕ k ) = ∑ k = 1 K r k [ n ] cos ⁡ ( 2 π f s ∑ i = 1 n f k [ i ] + ϕ k )

where

r k [ n ] = r k ( n T ) is the discrete-time varying amplitude envelope f k [ n ] = 1 T ∫ ( n − 1 ) T n T f k ( t )   d t is the discrete-time backward difference instantaneous frequency.

This is equivalent to

y [ n ] = ∑ k = 1 K r k [ n ] cos ⁡ ( θ k [ n ] )

where

θ k [ n ] = 2 π f s ∑ i = 1 n f k [ i ] + ϕ k = θ k [ n − 1 ] + 2 π f s f k [ n ] for all n > 0

and

θ k [ 0 ] = ϕ k .