弱弱地问个问题,为什么我的频谱图会有额外频率成分?

banforks

程序的代码如下:  不明白900HZ处为什么会有频谱,指教一下:

%
clear
fs=1000;
Ts=1/fs;
n=1024;
for i=1:n
x(i)=(1+cos(16*pi*i*Ts)+cos(40*pi*i*Ts))*cos(240*pi*i*Ts);
end
figure(1);
plot(x);
title('原始信号');
f=0:length(x)-1;
figure(2);
plot(fs*f/length(x),abs(fft(x))/fs);
%stem(abs(fft(x)));
title('频谱图');

破凰

DFT结果是共轭对称的
本版有单边谱和双边谱完整的程序，搜索一下便知。

[ 本帖最后由 破凰 于 2007-8-12 17:58 编辑 ]

banforks

谢谢楼上的仁兄给予解答!

HolySaint

可以用滤波器滤掉

楼主可以试试多个正弦波叠加

那样的频谱看着更明显

滤波效果也很明显

songzy41

因为信号是实数序列，它的DFT是共轭对称的，采样频率是500，在谱图中0~500反映正频率的分量，而500~1000是反映负频率的分量。正如破凰所讲是单边谱和双边谱的问题，而不是能用滤波器滤掉的。