Additive synthesis

Explanation

The sounds that are heard in everyday life are not characterized by a single frequency. Instead, they consist of a sum of pure sine frequencies, each one at a different amplitude. When humans hear these frequencies simultaneously, we can recognize the sound. This is true for both "non-musical" sounds (e.g. water splashing, leaves rustling, etc.) and for "musical sounds" (e.g. a piano note, a bird's tweet, etc.). This set of parameters (frequencies, their relative amplitudes, and how the relative amplitudes change over time) are encapsulated by the timbre of the sound. Fourier analysis is the technique that is used to determine these exact timbre parameters from an overall sound signal; conversely, the resulting set of frequencies and amplitudes is called the Fourier series of the original sound signal.

In the case of a musical note, the lowest frequency of its timbre is designated as the sound's fundamental frequency. For simplicity, we often say that the note is playing at that fundamental frequency (e.g. "middle C is 261.6 Hz"), even though the sound of that note consists of many other frequencies as well. The set of the remaining frequencies is called the overtones (or the harmonics, if their frequencies are integer multiples of the fundamental frequency) of the sound. In other words, the fundamental frequency alone is responsible for the pitch of the note, while the overtones define the timbre of the sound. The overtones of a piano playing middle C will be quite different from the overtones of a violin playing the same note; that's what allows us to differentiate the sounds of the two instruments. There are even subtle differences in timbre between different versions of the same instrument (for example, an upright piano vs. a grand piano).

Additive synthesis aims to exploit this property of sound in order to construct timbre from the ground up. By adding together pure frequencies (sine waves) of varying frequencies and amplitudes, we can precisely define the timbre of the sound that we want to create.

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:

: \begin{align} y(t) &= \frac{a_0}{2} + \sum_{k=1}^{\infty} \left[ a_k \cos(2 \pi k f_0 t ) - b_k \sin(2 \pi k f_0 t ) \right] \ &= \frac{a_0}{2} + \sum_{k=1}^{\infty} r_k \cos\left(2 \pi k f_0 t + \phi_k \right) \ \end{align}

where

• f_0 = 1/P is the fundamental frequency of the waveform and is equal to the reciprocal of the period,
• a_k = r_k \cos(\phi_k) = 2 f_0 \int_{0}^P y(t) \cos(2 \pi k f_0 t), dt, \quad k \ge 0
• b_k = r_k \sin(\phi_k) = -2 f_0 \int_{0}^P y(t) \sin(2 \pi k f_0 t), dt, \quad k \ge 1
• r_k = \sqrt{a_k^2 + b_k^2} is the amplitude of the kth harmonic,
• \phi_k = \operatorname{atan2}(b_k, a_k) is the phase offset of the kth harmonic. atan2 is the four-quadrant arctangent function,

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 \phi_k are the amplitude, frequency, and the phase offset, respectively, of the kth 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. | chapter-url = https://ccrma.stanford.edu/~jos/parshl/Additive_Synthesis.html | access-date = 11 January 2015 | chapter-url = https://ccrma.stanford.edu/~jos/sasp/Additive_Synthesis_Early_Sinusoidal.html | access-date = 9 January 2012

: y(t) = \sum_{k=1}^{K} r_k(t) \cos\left(2 \pi f_k t + \phi_k \right),

where f_k is the constant frequency of kth 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 \pi. 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. | author-link = Curtis Roads | url-access=limited 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. | author-link = Curtis Roads | url-access=limited

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. | access-date = 14 January 2012

Wavetable synthesis

Main article: 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. |access-date = 21 May 2005 |archive-url = https://web.archive.org/web/20130615202748/http://musicdsp.org/files/Wavetable-101.pdf |archive-date = 15 June 2013

Group additive synthesis

Group additive synthesis | access-date = 12 May 2011 | archive-url = https://web.archive.org/web/20110606200135/https://ccrma.stanford.edu/~jos/sasp/Group_Additive_Synthesis.html | archive-date = 6 June 2011 | url-status= live}} | chapter-url = http://quod.lib.umich.edu/cgi/p/pod/dod-idx?c=icmc;idno=bbp2372.1990.015

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

|author-link2 = Thomas F. Quatieri |access-date = 9 December 2013 |archive-url = https://web.archive.org/web/20120521071601/http://www.ll.mit.edu/publications/journal/pdf/vol01_no2/1.2.3.speechprocessing.pdf |archive-date = 21 May 2012

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. | access-date = 13 January 2012

Additive analysis/resynthesis has been employed in a number of techniques including Sinusoidal Modelling, | access-date = 9 January 2012

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 See "Excitation plus resonances voice model" (p. 51) 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).

Applications

Musical instruments

Main article: Synthesizer, Electronic musical instrument, Software synthesizer

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

Speech synthesis

Main article: Speech synthesis

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

Later, in the 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. |script-title=ja:複合正弦波による音声合成 | trans-title = Speech Synthesis by Composite Sinusoidal Wave |script-chapter=ja:複合正弦波による簡易な音声合成法 |trans-chapter=Simple Speech Synthesis method by Composite Sinusoidal Wave used on a singing speech synthesis feature on the Yamaha CX5M (1984), is known to use a similar approach which was independently developed during 1966–1979. | author-link2 = Fumitada Itakura |author-link = Fumitada Itakura |publication-date = April 2004 |archive-date = 24 May 2022 |access-date = 24 October 2014 |archive-url = https://web.archive.org/web/20220524112852/https://www.icacommission.org/Proceedings/ICA2004Kyoto/pdf/We3.D.pdf |url-status = dead |chapter-url = http://dl.acm.org/citation.cfm?id=131158 |url-access = registration | chapter-url= https://ccrma.stanford.edu/~bilbao/booktop/node14.html | editor-last = Greenebaum | editor-first = K.

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, William Thomson (later ennobled as 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_\mathrm{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_\mathrm{s},.

Beginning with (),

: y(t) = \sum_{k=1}^{K} r_k(t) \cos\left(2 \pi \int_0^t f_k(u)\ du + \phi_k \right)

and sampling at discrete times t = nT = n/f_\mathrm{s} , results in

: \begin{align} y[n] & = y(nT) = \sum_{k=1}^{K} r_k(nT) \cos\left(2 \pi \int_0^{nT} f_k(u)\ du + \phi_k \right) \ & = \sum_{k=1}^{K} r_k(nT) \cos\left(2 \pi \sum_{i=1}^{n} \int_{(i-1)T}^{iT} f_k(u)\ du + \phi_k \right) \ & = \sum_{k=1}^{K} r_k(nT) \cos\left(2 \pi \sum_{i=1}^{n} (T f_k[i]) + \phi_k \right) \ & = \sum_{k=1}^{K} r_k[n] \cos\left(\frac{2 \pi}{f_\mathrm{s}} \sum_{i=1}^{n} f_k[i] + \phi_k \right) \ \end{align} where

: r_k[n] = r_k(nT) , is the discrete-time varying amplitude envelope : f_k[n] = \frac{1}{T} \int_{(n-1)T}^{nT} f_k(t)\ dt , is the discrete-time backward difference instantaneous frequency.

This is equivalent to

: y[n] = \sum_{k=1}^{K} r_k[n] \cos\left( \theta_k[n] \right)

where

: \begin{align} \theta_k[n] &= \frac{2 \pi}{f_\mathrm{s}} \sum_{i=1}^{n} f_k[i] + \phi_k \ &= \theta_k[n-1] + \frac{2 \pi}{f_\mathrm{s}} f_k[n] \ \end{align} for all n0, and : \theta_k[0] = \phi_k. ,

References

References

1. Mottola, Liutaio. (31 May 2017). "Table of Musical Notes and Their Frequencies and Wavelengths".
2. "Fundamental Frequency and Harmonics".
3. [http://www.klingbeil.com/spear/ SPEAR Sinusoidal Partial Editing Analysis and Resynthesis for Mac OS X, MacOS 9 and Windows]
4. "Loris Software for Sound Modeling, Morphing, and Manipulation".
5. [http://mtg.upf.edu/technologies/sms SMSTools application for Windows]
6. [https://arss.sourceforge.net/ ARSS: The Analysis & Resynthesis Sound Spectrograph]
7. {{harvnb. Loscos. 2007
8. {{harvnb. Loscos. 2007
9. Rubin. P.E.. (1980). "Sinewave Synthesis Instruction Manual (VAX)". link. Internal Memorandum. Haskins Laboratories, New Haven, CT
10. Prestini, Elena. (2004). "The Evolution of Applied Harmonic Analysis: Models of the Real World". Birkhäuser Boston.
11. Fourier, Jean Baptiste Joseph. (1822). "Théorie analytique de la chaleur". Chez Firmin Didot, père et fils.
12. (1875). The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science. Taylor & Francis. (July 2024)
13. Thomson, Sir W.. (1878). "Harmonic analyzer". Taylor and Francis.
14. Helmholtz, von, Hermann. (1863). "Die Lehre von den Tonempfindungen als physiologische Grundlage für die Theorie der Musik". Leopold Voss.
15. Christensen, Thomas Street. (2002). "The Cambridge History of Western Music". Cambridge University Press.
16. Cahan, David. (1993). "Hermann von Helmholtz and the foundations of nineteenth-century science". University of California Press.
17. von Helmholtz, Hermann. (1875). "On the sensations of tone as a physiological basis for the theory of music". Longmans, Green, and co..
18. Miller, Dayton Clarence. (1926). "The Science of Musical Sounds". The Macmillan Company.
19. Russell, George Oscar. (1936). "Year book - Carnegie Institution of Washington (1936)". Carnegie Institution of Washington.
20. Lodge, John E.. (April 1938). "Odd Laboratory Tests Show Us How We Speak: Using X Rays, Fast Movie Cameras, and Cathode-Ray Tubes, Scientists Are Learning New Facts About the Human Voice and Developing Teaching Methods To Make Us Better Talkers". Popular Science Publishing.
21. Comerford, P.. (1993). "Simulating an Organ with Additive Synthesis". Computer Music Journal.
22. (1940). "Institute News and Radio Notes". Proceedings of the IRE.
23. Douglas, A.. (1948). "Electrotonic Music". Proceedings of the Royal Musical Association.
24. Douglas, Alan Lockhart Monteith. (1957). "The Electrical Production of Music". Macdonald.
25. Pejrolo, Andrea. (2007). "Acoustic and MIDI orchestration for the contemporary composer". Elsevier.
26. Weidenaar, Reynold. (1995). "Magic Music from the Telharmonium". Scarecrow Press.
27. Moog, Robert A.. (October–November 1977). "Electronic Music". Journal of the Audio Engineering Society.
28. Olsen, Harvey. (14 December 2011). "Leslie Speakers and Hammond organs: Rumors, Myths, Facts, and Lore". Hammond Organ in the U.K..
29. Cooper, F. S.. (May 1951). "The interconversion of audible and visible patterns as a basis for research in the perception of speech". Proc. Natl. Acad. Sci. U.S.A..
30. Holzer, Derek. (22 February 2010). "A brief history of optical synthesis".
31. Vail, Mark. (1 November 2002). "Eugeniy Murzin's ANS – Additive Russian synthesizer".
32. Young, Gayle. "Oscillator Bank (1959)".
33. Young, Gayle. "Spectrogram (1959)".
34. Luce, David Alan. (1963). "Physical correlates of nonpercussive musical instrument tones". Massachusetts Institute of Technology.
35. Beauchamp, James. (17 November 2009). "The Harmonic Tone Generator: One of the First Analog Voltage-Controlled Synthesizers".
36. Beauchamp, James W.. (October 1966). "Additive Synthesis of Harmonic Musical Tones". Journal of the Audio Engineering Society.
37. "RMI Harmonic Synthesizer". Synthmuseum.com.
38. "Allen Organ Company".
39. Cosimi, Enrico. (20 May 2009). "EMS Story - Prima Parte". Audio Accordo.it.
40. Hinton, Graham. (2002). "EMS: The Inside Story". Electronic Music Studios (Cornwall).
41. (1979). "The New Sound of Music". BBC.
42. Twyman, John. (1 November 2004). "(inter)facing the music: The history of the Fairlight Computer Musical Instrument". Unit for the History and Philosophy of Science, University of Sydney.
43. Street, Rita. (8 November 2000). "Fairlight: A 25-year long fairytale". IMAS Publishing UK.
44. (1978). "Computer Music Journal".
45. Leider, Colby. (2004). "Digital Audio Workstation". [[McGraw-Hill]].
46. Joel, Chadabe. (1997). "Electric Sound". Prentice Hall.
47. "Kawai K5000 {{!}} Vintage Synth Explorer".
48. "Kawai K5000R & K5000S".