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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 其他交流 » 复数数组的傅里叶变换，求一个fortran的例子
北京石油化工学院2026年研究生招生接收调剂公告
52/1返回列表
查看: 2514  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chenxiankai

禁虫 (正式写手)
本帖内容被屏蔽

» 猜你喜欢
1楼 2014-09-02 08:06:16
已阅   关注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的回帖
查看全部 19 个回答

chenxiankai

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

chenxiankai

禁虫 (正式写手)
本帖内容被屏蔽
4楼2014-09-02 10:59:30
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
4楼: Originally posted by chenxiankai at 2014-09-02 10:59:30
请问，我的复数数组，是一个与时间相关的数组，从-t-t之间的，所以数组个数是奇数，因为包括了t=0的时候，请问我该怎么来写这个呢...

可以补一个数，取最末三个数的平均。
赞一下 回复此楼
5楼2014-09-02 11:17:06
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 083000学硕274求调剂 +7 Li李鱼 2026-03-26 7/350 2026-03-28 08:01 by baoball
[考研] 压国家一区线，求导师收留，有恩必谢！ +3 迷人的哈哈 2026-03-28 3/150 2026-03-28 07:34 by ACS Nano——
[考研] 304求调剂 +4 曼殊2266 2026-03-27 4/200 2026-03-28 07:22 by Iveryant
[考研] 085701环境工程，267求调剂 +13 minht 2026-03-26 13/650 2026-03-28 07:14 by mmm just
[考研] 求调剂 +8 张zz111 2026-03-27 9/450 2026-03-28 03:41 by fmesaito
[考研] 085600，材料与化工321分调剂 +4 大馋小子 2026-03-27 6/300 2026-03-27 14:11 by 松花缸1201
[考研] 0703一志愿9，初试成绩：338，四六级已过，有科研经历，求调剂！ +3 Zuhui0306 2026-03-25 3/150 2026-03-27 14:09 by shangxh
[考研] 085601 材料工程 313分 求调剂 +5 Ong3 2026-03-27 5/250 2026-03-27 12:24 by goldfish51
[考研] 22408 359分调剂 +3 Qshers 2026-03-27 3/150 2026-03-27 12:22 by wxiongid
[论文投稿] Journal of Mechanical Science and Technology +3 Russ_ss 2026-03-25 5/250 2026-03-27 10:49 by 陆小果画大饼
[考研] 求调剂323材料与化工 +7 1124361 2026-03-24 7/350 2026-03-27 10:22 by wangjy2002
[硕博家园] 北京林业大学硕导招生广告 +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
[考研] 316求调剂 +5 江辞666 2026-03-26 5/250 2026-03-27 08:08 by hypershenger
[考研] 329求调剂 +7 钮恩雪 2026-03-25 7/350 2026-03-27 04:28 by wxiongid
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
[考研] 324求调剂 +4 wysyjs25 2026-03-21 4/200 2026-03-26 20:38 by fmesaito
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-26 3/150 2026-03-26 19:57 by nihaoar
[考研] 0854电子信息求调剂 +7 α____ 2026-03-22 9/450 2026-03-25 13:37 by α____
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
信息提示
关闭
请填处理意见
关闭