24小时热门版块排行榜    

查看: 1028  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

个男子

新虫 (初入文坛)

[求助] 快速傅里叶变换

大家早上好,请教一个问题哈。小弟用快速傅里叶变换得到对应信号的实部和虚部,现在想利用实部和虚部绘制曲线,利用这个公式:快速傅里叶变换,其中d0是直流分量,dn是第n点振幅,w1是2PI/T,最后那个是虚部除以实部的反正切(0~PI)。可是求出来的曲线为啥不对那
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
3楼: Originally posted by 个男子 at 2015-01-31 17:51:44
theta是不是只能是PI/2的倍数,还是由计算得来虚部比实部的反正切的值。我现在用的模型和理论上的有些不同,为:y=d0+(-1)^n *dn* sin(2PI/T *n + theta).....其中theta取值为-PI到PI。。。但是拟合出来的曲线n=1, ...

正负号对应象限

[ 发自小木虫客户端 ]
4楼2015-01-31 18:09:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

feixiaolin

荣誉版主 (文坛精英)

优秀版主

dn=sqrt(re.FFT(n)^2+im.FFT(n)^2)
theta=反正切要讨论【根据正负号组合】
2楼2015-01-31 15:59:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

个男子

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by feixiaolin at 2015-01-31 15:59:26
dn=sqrt(re.FFT(n)^2+im.FFT(n)^2)
theta=反正切要讨论【根据正负号组合】

theta是不是只能是PI/2的倍数,还是由计算得来虚部比实部的反正切的值。我现在用的模型和理论上的有些不同,为:y=d0+(-1)^n *dn* sin(2PI/T *n + theta).....其中theta取值为-PI到PI。。。但是拟合出来的曲线n=1,2的时候还行,越大差的越大,不知道问题出在哪了
3楼2015-01-31 17:51:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 YXM1FDCSE8Bn 2026-10-31 8/400 2026-11-02 12:31 by n0P4fnsXKbBK
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 YXM1FDCSE8Bn 2026-11-01 6/300 2026-11-02 10:02 by dcc2yOIRGGTT
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-10-31 5/250 2026-11-02 09:42 by dcc2yOIRGGTT
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 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一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:34 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 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.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见