FFT Analysis of two signals
显示 更早的评论
Hello everyone,
I have a beginners question about FFT analysis. I have two signals and I want to see the amplitudespectrum of them. I used the following code which was derived from the program desribed under mathworks.
Code:
Fs = 200; % Sampling frequency
T = 1/Fs; % Sample time
L = length (marker(1,2).C7); % Length of signal
t = (0:L-1)*T; % Time vector
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
y= marker(1,2).C7;
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
figure(1)
plot(f,2*abs(Y(1:NFFT/2+1)))
title('....')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
marker(1,2).C7 is the signal I want to analyse.
The following problem occured: If I run this program on two different signals I get the same amplitudespectrum for both although the signals are different. This cannnot be true as I wanted to compare the amplitudespectrum of the signal before and after it was filtered. Although they were filtered and they looked different when I ploted them the amplitudespectrum plotted with the above code showed no differences at all. Can anybody point out where I made a mistake?
Thank you for your help,
Christoph
2 个评论
Honglei Chen
2012-3-14
It's not clear to me what you are comparing. Are you comparing two signals before filtering, or you are comparing one signal before and after filtering?
David Young
2012-3-14
I think you may need to show the complete code that compares the two signals.
采纳的回答
Wayne King
2012-3-15
Let's assume that x is one signal and y is the other signal. If you have the Signal Processing Toolbox, you can quickly compare their power spectral densities with the following. fs below is your sampling frequency, so you need to put the correct value in for that.
psdestx = psd(spectrum.periodogram,x,'Fs',fs,'NFFT',length(x));
psdesty = psd(spectrum.periodogram,y,'Fs',fs,'NFFT',length(y));
plot(psdestx.Frequencies,10*log10(psdestx.Data));
title('PSD of x'); xlabel('Hz'); ylabel('dB/Hz');
figure;
plot(psdesty.Frequencies,10*log10(psdesty.Data));
title('PSD of y'); xlabel('Hz'); ylabel('dB/Hz');
更多回答（1 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Spectral Measurements 的更多信息
2012-3-14
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
