how to generate a smooth sine wave with time varying frequency in matlab

Question

I'm currently designing a knock sensor for my FSAE team's car and I'd like to be able to test various software strategies with a fairly realistic signal (at least more realistic than a couple sine waves added together with some white noise).

From some reading, the frequencies I'll be getting will all be related to a fundamental frequency directly related to engine speed. As engine speed varies continuously, I'd like to be able to take a vector of times and transform it into a vector of frequencies that smoothly and randomly transition within the range 20-170 Hz, ideally with a uniform distribution. I can then create a couple sine waves based on these frequencies and a couple multiples, add them, add some white noise, add a knock at random points, then feed that into my knock detection code.

It's generating a smooth, uniformly distributed frequency vector that is my current issue, but I'm also open to other solutions to the problem.

Answer 1

rpatel: That's an interesting Tex. Instruments Application Note. On pages 8 & 9 they mention the sliding-mode DFT (known as the "sliding DFT") and the Goertzel algorithm. They wrote:

"The disadvantages of these techniques are that the computational overhead is slightly higher and the data memory requirements increase, because a sample-block-length number of samples (N) must be retained for the calculations."

I don't believe that statement is true. The sliding DFT computes updated spectral results upon the arrival of each new input sample, and past input samples need NOT be saved (retained). The Goertzel algorithm computes updated temporary results upon the arrival of each new input sample, and past input samples need NOT be saved (retained). Upon the arrival of the Nth input sample the Goertzel algorithm computes its final spectral result. (I'm merely telling you this for your information.) Maybe the sliding DFT or the Goertzel algorithm would be useful to you. But if the Tex. Instr. Application Note's method does what you want then stick with its method.

To start, try the following code to generate an increasing-freq chirp signal plus its first harmonic, followed by a decreasing-freq chirp signal plus its first harmonic:

Fs = 1000;  % Sample rate in Hz

Dur = 1;    % One second time duration

Time = 0:1/Fs:Dur;  

% Generate the fundamental increasing-freq chirp

Fstart = 100;   % Start at 100 Hz

Fstop = 150;    % Final freq of chirp

Sig_1 = chirp(Time,Fstart,Dur,Fstop); % Sweep from 100 –to- 150 Hz 

% Generate 1st harmonic of the increasing-freq fundamental chirp

Fstart = 200;    % Start at 200 Hz

Fstop = 300;     % Final freq of chirp

Sig_2 = chirp(Time,Fstart,Dur,Fstop); % Sweep from 200 –to- 300 Hz

Sig = Sig_1 + Sig_2;     % Increasing-freq chirp plus 1st harmonic

Sig = [Sig, fliplr(Sig)]; % Increasing-freq chirp followed by decreasing-freq chirp

Change the above variables to suit your needs, add your noise and ping signal, and start experimenting with your 'ping detection' method. (Add higher-freq harmonics if necessary.)

Now if you need a fundamental sine wave signal that randomly changes in frequency (plus its harmonics) then that's a whole different ballgame. To satisfy that need I'd have to think about ways of using 'continuous multitone FSK (freq shift keying) signal generation' techniques. Let's keep in touch on this.

Answer 2

The wording of your question (your terminology) makes it difficult to understand exactly what you desire. Is your question, "Using Matlab, how do I generate a time-domain sequence representing a sinusoidal wave whose instantaneous frequency increases as a linear function of time so that its spectrum contains a broad band of spectral energy?" If so, then try the following example Matlab code:

Fs = 1000;       % Signal sample rate in Hz

Dur = 1;         % Signal's 1-second time duration

Time = 0:1/Fs:Dur;  

Fstart = 100;    % Start freq is 100 Hz

Fstop = 400;     % Final freq is 400 Hz

Sig = chirp(Time,Fstart,Dur,Fstop);   % Generate a "chirp" sequence

rpatel, generate the above 'Sig' signal, compute its DFT samples, and plot the DFT's spectral magnitude values. Maybe that's what you want.