应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 其他交流 » 复数数组的傅里叶变换，求一个fortran的例子
192/2返回列表上一页12
查看: 2750  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

feixiaolin

荣誉版主 (文坛精英)
引用回帖:
10楼: Originally posted by chenxiankai at 2014-09-02 14:30:14
那这样就不用吧中间那个时间t=0的那个数据给删除了，是吧？...

是
11楼2014-09-02 14:39:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxiankai

禁虫 (正式写手)
本帖内容被屏蔽
12楼2014-09-02 15:19:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
12楼: Originally posted by chenxiankai at 2014-09-02 15:19:52
请问那两个程序是不是计算速度很慢呢？...

不会的。
回复此楼
13楼2014-09-02 15:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxiankai

禁虫 (正式写手)
本帖内容被屏蔽
14楼2014-09-02 15:41:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenxiankai

禁虫 (正式写手)
本帖内容被屏蔽
15楼2014-09-02 15:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
再看一下这个：

http://blog.renren.com/blog/247544700/703436782?bfrom=01020100300
!!时域分解FFT算法
!!离散时域信号的序号对应的时刻为nΔ，离散频域信号的序号对应的频率为m/（nΔ）
!!离散时域信号保存在文件A1.TXT中
!! 变换后的离散频域信号保存在文件 FFT_RESULT.TXT
!!编译器为cvf6.6

      PROGRAM FFT
      PARAMETER(K=4)
      !COMPLEX用结构体定义复数
      TYPE CMPLX0
      REAL RE0,IM0
      END TYPE

      REAL M
      TYPE(CMPLX0)A(0:2**K-1),A0(0:2**K-1)

      OPEN(1,FILE='A1.TXT')
      DO I=0,2**K-1
           READ(1,*) A(I).RE0,A(I).IM0
      ENDDO
!     PRINT*,A
      CLOSE(1)
      DO I=0,2**K-1
           A0(I)=A(IFIX(CTY_I(I,K)))
      ENDDO
!     PRINT*,A0
!     PAUSE

      DO I=1,K                                                      !大循环
           DO J0=1,2**(K-I)                                !单元个数
                 M=0
                 DO J1=(J0-1)*2**I,(J0-1)*2**I+2**(I-1)-1 !每个单元的运算次数
                      A(J1)=PLUS(A0(J1),MULT(A0(J1+2**(I-1)),W(M,2.**K)))
                      A(J1+2**(I-1))=SUBT
     &(A0(J1),MULT(A0(J1+2**(I-1)),W(M,2.**K)))
                      M=M+2**(K-I)
                      IF(M>=2**(K-1)) M=0
                 ENDDO
           ENDDO
      A0=A
      ENDDO
      OPEN(2,FILE='FFT_RESULT.TXT',ACTION='WRITE')
      DO I=0,2**K-1
           IF(A(I).RE0==0) THEN
                 WRITE(2,100) A(I).IM0
           ELSEIF(A(I).IM0==0)THEN
                 WRITE(2,200) A(I).RE0
           ELSE
                 WRITE(2,300)A(I).RE0,A(I).IM0
           ENDIF
      ENDDO
100 FORMAT(F18.3,'i')
200 FORMAT(F18.3)
300 FORMAT(F18.3,'           +',F18.3,'i')

      CONTAINS
      !PLUS 定义复数的加法
      FUNCTION PLUS(Z1,Z2)
      TYPE(CMPLX0) Z1,Z2,PLUS
      PLUS.RE0=Z1.RE0+Z2.RE0
      PLUS.IM0=Z1.IM0+Z2.IM0
      ENDFUNCTION
      !SUBTRACTION 定义复数的减法
      FUNCTION SUBT(Z1,Z2)
      TYPE(CMPLX0) Z1,Z2,SUBT
      SUBT.RE0=Z1.RE0-Z2.RE0
      SUBT.IM0=Z1.IM0-Z2.IM0
      ENDFUNCTION
      !MULTIPLACATION 定义复数的乘法
      FUNCTION MULT(Z1,Z2)
      TYPE(CMPLX0) Z1,Z2,MULT
      MULT.RE0=Z1.RE0*Z2.RE0-Z1.IM0*Z2.IM0
      MULT.IM0=Z1.RE0*Z2.IM0+Z1.IM0*Z2.RE0
      ENDFUNCTION
      !W(N)**M 定义WN的M次方
      FUNCTION W(M,N)
      REAL M,N
      TYPE(CMPLX0) W
      W.RE0=COS(2*3.1415926*M/N)
      W.IM0=-SIN(2*3.1415926*M/N)
      ENDFUNCTION
      !CONTRARY I 求I的二进制逆序数
      FUNCTION CTY_I(I0,N)
      INTEGER I,A(N),N,IO
      REAL CTY_I
      A=0
      I=I0
      CTY_I=0
      DO K0=1,N
      IF(I-2**(N-K0)<0)THEN
           A(K0)=0
      ELSE
           A(K0)=1
           I=I-2**(N-K0)
      ENDIF
      CTY_I=CTY_I+A(K0)*2**(K0-1)
      ENDDO
      ENDFUNCTION
      END
赞一下 回复此楼
16楼2014-09-02 16:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guangcun

禁虫 (初入文坛)
本帖内容被屏蔽
17楼2015-02-04 19:23:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
17楼: Originally posted by guangcun at 2015-02-04 19:23:05
请问二维傅里叶变换可以用这个吗，还是只是一维可以 ？...

可以作为二维傅里叶变换调用的子函数
赞一下 回复此楼
18楼2015-02-04 20:23:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

575406210
就
19楼2015-10-13 15:53:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenxiankai 的主题更新
192/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 16/800 2026-08-20 19:26 by wuchongjun
[基金申请] 科研孤儿太难了 +9 我4大白菜 2026-08-20 9/450 2026-08-20 18:23 by 2549091707
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天系统多次维护，明天很可能放榜！ +9 zju2000 2026-08-16 10/500 2026-08-20 15:56 by Leogzhya
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 言之凿凿，明天官网会通知放榜？ +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗？20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭