24小时热门版块排行榜    

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

xiashaoyan

新虫 (正式写手)

[求助] ifft的问题

现在遇到一个问题就是有关ifft的使用,现在问题如下文所述:
注意:A中有负数,不都是正数
A =
    -8     1     6
    -3     5     7
     4     9     2
对Afft2得到A的谱pu

pu =

  23.0000 + 0.0000i -22.0000 + 0.0000i -22.0000 + 0.0000i
-13.0000 + 5.1962i   0.5000 +12.9904i -13.0000 + 0.0000i
-13.0000 - 5.1962i -13.0000 - 0.0000i   0.5000 -12.9904i

然后对A的pu乘以一个常数exp(i*0.5),即:new_pu = pu.*exp(i*0.5),得到新的一个谱 new_pu

new_pu =

  20.1844 +11.0268i -19.3068 -10.5474i -19.3068 -10.5474i
-13.8997 - 1.6725i  -5.7891 +11.6398i -11.4086 - 6.2325i
  -8.9174 -10.7926i -11.4086 - 6.2325i   6.6667 -11.1604i

然后对new_pu进行逆变换,得到类似于A的数据 new_A  = iff2(new_pu)
new_A =

  -7.0207 - 3.8354i   0.8776 + 0.4794i   5.2655 + 2.8766i
  -2.6327 - 1.4383i   4.3879 + 2.3971i   6.1431 + 3.3560i
   3.5103 + 1.9177i   7.8982 + 4.3148i   1.7552 + 0.9589i

现在很奇怪的是,对new_pu(复数)进行iff2得到的数据应该是double数据,不应该是复数域内的数值呀,但是如果对A的谱pu直接ifft2就得到实数,不是复数,现在我想对iff2(new_pu)后也得到一个实数,但是不知道怎么办?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bamaxingfu12

银虫 (小有名气)

反变换回去不应该是a.*exp(i*0.5)吗?
3楼2015-12-10 17:31:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

feixiaolin

荣誉版主 (文坛精英)

优秀版主

★ ★ ★ ★ ★
xiashaoyan(conanwj代发): 金币+5, 感谢应助 2015-12-19 21:25:21
设想
A =
    -8     1     6
    -3     5     7
     4     9     2
B=
    -8     1     6
    -3     5     7
     4     9     2
构建
C=A+iB
    对C变换,只取实部,如何?
2楼2015-12-07 07:28:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +6 Lanmanbaby 2026-08-09 8/400 2026-08-09 16:24 by zhjj370
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见