FFT Analysis of two signals

5 次查看（过去 30 天）
Christoph Bauer
Christoph Bauer 2012-3-14
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
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
David Young 2012-3-14
I think you may need to show the complete code that compares the two signals.

请先登录，再进行评论。

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

采纳的回答

Wayne King
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 个）

Christoph Bauer
Christoph Bauer 2012-3-15
I want to compare the same signal before and after filtering.
There is no other code to compare the signals yet. The signal was measured and later filtered in another software, using a butterworth 4th order filter. After that I imported the signal into matlab, first the raw signal and then the filtered signal. Then I ploted both and used the above code to compare the signal before and after filtering.

类别

在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by