24小时热门版块排行榜    

查看: 1941  |  回复: 14

sxf2012

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
阿乐斗士: 金币+4, ★★★很有帮助 2014-10-25 20:55:37
CODE:
dx = x(2)-x(1);
w = linspace(-2*pi/dx,2*pi/dx,length(x));
Fy = fftshift(fft(y));
plot(w,abs(Fy))

11楼2014-10-24 18:16:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
阿乐斗士: 金币+2, 有帮助 2014-10-25 20:56:17
采样点没有零点,要做相位校正,如果你只关心频谱实部的话。快速傅里叶变换采样点的间隔决定谱宽,具体采样点不用给出,只要知道第一个采样点,采样间隔和采样点数三者即可。如果随机采样做离散傅里叶变换或压缩感知重建,才需要具体的采样点信息。

[ 发自手机版 http://muchong.com/3g ]
12楼2014-10-24 19:50:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿乐斗士

木虫 (著名写手)

引用回帖:
11楼: Originally posted by sxf2012 at 2014-10-24 18:16:25
dx = x(2)-x(1);
w = linspace(-2*pi/dx,2*pi/dx,length(x));
Fy = fftshift(fft(y));
plot(w,abs(Fy))

求大神解释一下,不太懂,愿意增加金币给你
傅里叶变换
GBC8I01PF{U5LB6U9JDJ_]S.jpg

且行且珍惜
13楼2014-10-25 20:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿乐斗士

木虫 (著名写手)

引用回帖:
12楼: Originally posted by 戴钢盔的猪头 at 2014-10-24 19:50:50
采样点没有零点,要做相位校正,如果你只关心频谱实部的话。快速傅里叶变换采样点的间隔决定谱宽,具体采样点不用给出,只要知道第一个采样点,采样间隔和采样点数三者即可。如果随机采样做离散傅里叶变换或压缩感知 ...

该怎么做相位校正呢?
且行且珍惜
14楼2014-10-25 20:56:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

引用回帖:
14楼: Originally posted by 阿乐斗士 at 2014-10-25 20:56:09
该怎么做相位校正呢?...

这个一两句说不清楚。你采样缺失零点,实际上每个频率成分在第一个采样点相对零点有相移,重建出来的频谱实部不是严格的洛伦兹线性,有一定的扭曲,零频除外。频率离零频越远,产生的相移越大。如果只关心模值频谱,则没有关系。如果关心实部频谱,你会发现频率越高,线型扭曲越严重。简单的相位校正即零级相位校正,是令频谱每一点同乘一个exp(i*φ),对损失的相位进行弥补,其中φ是相位角。实际上这样只能令某个关注的谱峰及其临近的谱峰的相位得到校正。因为每个频率成分相移不同,故实际上频谱上每一点应当用不同的相位角进行相位校正,就是一级相位校正。具体请自行学习。

[ 发自手机版 http://muchong.com/3g ]
15楼2014-10-26 00:44:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿乐斗士 的主题更新
信息提示
请填处理意见