MATLAB的答案

奇怪的傅里叶变换的结果

2视图(30天)
大卫
大卫 2013年6月19日
你好,所有人
Matlab我很新,我想用它来发现一个信号的频率。
在30赫兹信号采样,长约3秒(88)样品。看着似乎有一个周期振荡信号在6-8Hz,这正是我希望看到的傅里叶变换。但是看起来不像我预期的傅里叶变换(我曾希望在别处在6-8Hz飙升和噪音)。
我也是新傅里叶变换,所以我用一个例子从Matlab turorials这里 http://www.mathworks.co.uk/help/matlab/ref/fft.html
信号是这样的
http://img801.imageshack.us/img801/6610/vud.png
和英国《金融时报》是这样的
http://img405.imageshack.us/img405/531/c12l.png
这是一个最小的工作示例代码
y = (207 206 177 203 204 203 183 201 173 208 207 194 205 207 176 206 198…
205 208 207 191 184 201 206 208 185 193 205 204 202 189 205 207 198…
185 194 209 210 197 201 208 208 191 204 209 209 199 205 210 211 203…
194 196 206 194 194 208 208 212 196 203 210 204 192 193 209 210 207…
202 209 206 194 194 197 205 197 197 207 201 206 192 192 196 196 196…
199 205 205];
samplingRate = 30;
样品= 88;
timeSeries = (0: samples-1)。/ samplingRate;
这个信号%的阴谋
情节(timeSeries y);
包含(“时间”)
ylabel (的信号值)
%找到英国《金融时报》
NFFT = 2 ^ nextpow2(样本);%的下一个2的幂的长度的信号
谱= fft (y, NFFT) /样品;
频率= samplingRate / 2 * linspace (0, 1, NFFT / 2 + 1);
%绘制英国《金融时报》
情节(频率、2 * abs(光谱(1:NFFT / 2 + 1)))
标题(单面振幅谱的pixelSeries (t)”)
包含(的频率(赫兹))
ylabel (“FFT (pixelSeries)”)
我不确定有问题的代码,或者信号是否太短,噪音对英国《金融时报》寻找什么吗?
非常感谢你的帮助
戴夫

在回答这个问题。

接受的答案

韦恩王
韦恩王 2013年6月19日
嗨,你是正确的。有一个峰值7赫兹左右。几个事情,信号有一个非零的意思,导致直流(零频率分量)掩盖了情节。
信号是短暂的,所以我填补它与零长度的两倍,帮助可视化高峰。
y = (207 206 177 203 204 203 183 201 173 208 207 194 205 207 176 206 198…
205 208 207 191 184 201 206 208 185 193 205 204 202 189 205 207 198…
185 194 209 210 197 201 208 208 191 204 209 209 199 205 210 211 203…
194 196 206 194 194 208 208 212 196 203 210 204 192 193 209 210 207…
202 209 206 194 194 197 205 197 197 207 201 206 192 192 196 196 196…
199 205 205];
y, y =去趋势(0);
npad = 2 *长度(y);
ydft = fft (y, npad);
频率=破发/ npad: 15;
情节(频率、abs (ydft (1: npad / 2 + 1)))
包含(“赫兹”);ylabel (“级”);
1评论
大卫
大卫 2013年6月19日
它完美地解决,多谢韦恩!

登录置评。

更多的答案(0)

类别

找到更多的在离散傅里叶余弦变换在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!