Main Archive Specials Wiki | FAQ Links Submit Forum

 Fast sine and cosine calculationType : waveform generationReferences : Lot's or references... Check Julius O. SMith mainlyCode : init: float a = 2.f*(float)sin(Pi*frequency/samplerate); float s[2]; s[0] = 0.5f; s[1] = 0.f; loop: s[0] = s[0] - a*s[1]; s[1] = s[1] + a*s[0]; output_sine = s[0]; output_cosine = s[1]

 CommentsAdded on : 05/04/03 by DFLComment : Yeah, this is a cool trick! :) FYI you can set s[0] to whatever amplitude of sinewave you desire. With 0.5, you will get +/- 0.5 Added on : 05/04/03 by bigtick[ AT ]pastnotecut[ DOT ]orgComment : After a while it may drift, so you should resync it as follows: const float tmp=1.5f-0.5f*(s[1]*s[1]+s[0]*s[0]); s[0]*=tmp; s[1]*=tmp;     This assumes you set s[0] to 1.0 initially. 'Tick        Added on : 08/04/03 by DFLComment : Just to expalin the above "resync" equation (3-x)/2 is an approximation of 1/sqrt(x) So the above is actually renormalizing the complex magnitude. [ sin^2 (x) + cos^2(x) = 1 ] Added on : 15/05/03 by nigelComment : This is the Chamberlin state variable filter specialized for infinite Q oscillation. A few things to note: Like the state variable filter, the upper frequency limit for stability is around one-sixth the sample rate. The waveform symmetry is very pure at low frequencies, but gets skewed as you get near the upper limit. For low frequencies, sin(n) is very close to n, so the calculation for "a" can be reduced to a = 2*Pi*frequency/samplerate. You shouldn't need to resync the oscillator--for fixed point and IEEE floating point, errors cancel exactly, so the osciallator runs forever without drifting in amplitude or frequency.Added on : 02/11/03 by antiprosynthesis[ AT ]hotmail[ DOT ]comComment : I made a nice little console 'game' using your cordic sinewave approximation. Download it at http://users.pandora.be/antipro/Other/Ascillator.zip (includes source code). Just for oldschool fun :).Added on : 22/12/04 by hplusComment : Note that the peaks of the waveforms will actually be between samples, and the functions will be phase offset by one half sample's worth. If you need exact phase, you can compensate by interpolating using cubic hermite interpolation.Added on : 24/07/07 by more on that topic[ DOT ][ DOT ][ DOT ]Comment : ... can be found in Jon Datorro, Effect Design, Part 3, a paper that can be easily found in the web. Funny, this is just a complex multiply that is optimized for small angles (low frequencies) When the CPU rounding mode is set to nearest, it should be stable, at least for small frequencies.Added on : 24/07/07 by mail[ AT ]mroc[ DOT ]deComment : More on that can be found in Jon Datorro, Effect Design, Part 3, a paper that can be easily found in the web. Funny, this is just a complex multiply that is optimized for small angles (low frequencies) When the CPU rounding mode is set to nearest, it should be stable, at least for small frequencies. Added on : 21/09/08 by bob[ AT ]yahoof[ DOT ]comComment : How do I set a particular phase for this? I've tried setting s[0] = cos(phase) and s[1] = sin(phase), but that didn't seem to be accurate enough. Thanks

 Add your own comment Comments are displayed in fixed width, no HTML code allowed! Email: Comment: Are you human?

Site created and maintained by Bram
Graphic design by line.out | Server sponsered by fxpansion