feilei16 发表于 2010-10-30 17:25:10

离散傅里叶变换出现问题

进行一个简单的信号从时域到频域变换,为什么会出现这样的问题,请高手指点一下。谢谢。附带图4为时域图,3转换后的。程序如下
load t.txt
t=t;
load p.txt
x=p;
plot(t,x)
Y = fft(x,512);
Pyy = Y.* conj(Y) / 512;
f = (0:256)/512;
plot(f,Pyy(1:257))

allenfieldin 发表于 2010-10-30 18:35:21

1# feilei16
从图像上看其实挺正常,因为你的时间域上的离散点的频率较为单一(图像从2s以后貌似都是按一个频率如omiga=2波动,而2s前的频率成分也较为单一),因此谱分析以后的结果会在那个单一频率omiga=2附近谱强度无穷大,而在其他频率处为零,就是这个意思

feilei16 发表于 2010-10-31 10:04:51

时域上除去2s之前的不说,2s后的也应该有个固定的频率值,频域上也不该在0的地方有个幅值。而且2s前的怎么说也应该有个小毛刺吧,怎么会那么平整呢?

allenfieldin 发表于 2010-10-31 14:18:43

3# feilei16
这与能量成分有关吧,比如说你的信号绝大部分都是一个频率omiga为主导,那么你分析出来的结果肯定在omiga=2这点很大,至于2s之前,对能量贡献很小,所得到的谱强度可能也会很小了,你觉得呢?你再把图像低频段放大再看看吧

ChaChing 发表于 2010-10-31 23:32:10

两个疑问:)
1.LZ的取样率多少? 10秒512点吗?
2.f为1*257, 怎图3大於450?

ZanDZ 发表于 2010-11-21 19:29:10

接楼上的两个疑问:
1. x是一个512点的信号吗?就是说图4画的信号,都进行了FFT,还是取一段进行FFT?
2. 和楼上一样,图3的横坐标是不是有问题。
还有,因为你信号的值都比较大,所以在频域,频率为0的地方就会有一个很大的值,所以说,是不是有些小的波动(图上在非零的地方有个小波动)就看不清了(如果把信号的均值控制在0,就不会有这个问题)。
页: [1]
查看完整版本: 离散傅里叶变换出现问题