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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 其他交流 » 复数数组的傅里叶变换，求一个fortran的例子
北京石油化工学院2026年研究生招生接收调剂公告
192/2返回列表上一页12
查看: 2504  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602 化工专硕 338分 求调剂 +9 路痴小琪 2026-03-27 9/450 2026-03-28 01:11 by 飞行日记西
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 考研调剂 +4 Sanmu-124 2026-03-26 4/200 2026-03-27 17:49 by kiokin
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 求调剂323材料与化工 +7 1124361 2026-03-24 7/350 2026-03-27 10:22 by wangjy2002
[考研] 材料调剂 +8 匹克i 2026-03-23 8/400 2026-03-27 08:11 by hypershenger
[考研] 343求调剂 +4 赠我一本书 2026-03-23 4/200 2026-03-27 00:40 by wxiongid
[考研] 【双一流院校新能源、环境材料，材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 8/400 2026-03-26 20:34 by Higraduate
[考研] 081700 调剂 267分 +11 迷人的哈哈 2026-03-23 11/550 2026-03-26 15:41 by zzll406
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 285求调剂 +3 AZMK 2026-03-24 3/150 2026-03-25 12:23 by userper
[考研] 化学调剂 +6 yzysaa 2026-03-21 6/300 2026-03-25 09:27 by aa331100
[考研] 0703化学调剂，求导师收 +7 天天好运来上岸� 2026-03-24 7/350 2026-03-24 20:26 by peike
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 一志愿国科过程所081700，274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
信息提示
关闭
请填处理意见
关闭