24小时热门版块排行榜    

查看: 2294  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 阿乐斗士 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 【长期招聘】复旦大学 AI 制药实验室 博士后招聘 +3 第三种态度 2026-10-28 9/450 2026-11-02 13:06 by Grace_PepLib
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-10-31 6/300 2026-11-02 12:44 by n0P4fnsXKbBK
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 YXM1FDCSE8Bn 2026-10-31 8/400 2026-11-02 12:31 by n0P4fnsXKbBK
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 12:25 by n0P4fnsXKbBK
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 BCdvBAEpt3zm 2026-11-01 4/200 2026-11-02 10:23 by dcc2yOIRGGTT
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 09:57 by dcc2yOIRGGTT
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 06:29 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 06:13 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 05:54 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 X5bHxOULxFqb 2026-10-30 27/1350 2026-11-02 05:28 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:05 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
信息提示
请填处理意见