Help with Chirp FFT

Question

Randa Qashoa 2019-5-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/463554-help-with-chirp-fft

评论： Daniel M 2019-11-21

I am trying to make a fft for a chirp signal with a frequency range of 15-25 MHz over a period of 10 micro seconds. When I attempt to plot the power vs. the frequency, I end up with a "square"-like wave from 15 MHz to 25 MHz and another one from 25 MHz to 35 MHz. I only need to generate the one wave that is from 15-25 MHz. Please find attached the plot and I have included the code below:

SPF=50e6;
frames = 16;
nfft = SPF/frames;
f_max = 25e6;
f_min = 15e6;
start_time=2e-6;
end_time=12e-6;
d_chirp = 10e-6;
Fs=2*f_max;
hchirp = dsp.Chirp('InitialFrequency',f_min,'TargetFrequency',f_max,... 
 'TargetTime', d_chirp,'SweepTime', d_chirp, 'SampleRate', SPF, 'SamplesPerFrame', nfft);
 
chirpData = (step(hchirp))';
freq_chirp=fft(chirpData,nfft);
power=abs(freq_chirp);
freq =(0:nfft-1)*Fs/nfft;
figure();
plot(freq,power);
xlabel('Frequency (Hz)');
ylabel('Power');
xlim([0 50e6]);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Hieu Nguyen 2019-11-21

I have the same problem; the bandwidth seems to double its length. Also, sweeping up and down have different results

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Daniel M 2019-11-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/463554-help-with-chirp-fft#answer_402671

在 MATLAB Online 中打开

You are plotting the incorrect frequency vector. fft() does not return the frequencies in order from most negative to most positive. Follow the example in the documentation for more information.

https://www.mathworks.com/help/matlab/ref/fft.html

Change your freq and power vectors to:

freq2 = (0:(nfft-1)/2)*Fs/nfft;
power2 = power(1:(nfft/2));

Also note that this isn't typically what is considered the power. This is just the energy. Power is energy squared. You would also have to divide by nfft to scale properly.