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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 帮忙看一下这个程序啊！！
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 830  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

破碎东瓜

金虫 (初入文坛)

[求助] 帮忙看一下这个程序啊！！ 已有1人参与

SUBROUTINE KKFFT(PR,PI,N,K,FR,FI,L,IL)
        DIMENSION PR(N),PI(N),FR(N),FI(N)
        DOUBLE PRECISION PR,PI,FR,FI,P,Q,S,VR,VI,PODDR,PODDI
        DO 20 IT=0,N-1
          M=IT
          IS=0
          DO 10 I=0,K-1
            J=M/2
            IS=2*IS+(M-2*J)
            M=J
10          CONTINUE
          FR(IT+1)=PR(IS+1)
          FI(IT+1)=PI(IS+1)
20        CONTINUE
        PR(1)=1.0
        PI(1)=0.0
        PR(2)=COS(6.283185306/N)
        PI(2)=-SIN(6.283185306/N)
        IF (L.NE.0) PI(2)=-PI(2)
        DO 30 I=3,N
          P=PR(I-1)*PR(2)
          Q=PI(I-1)*PI(2)
          S=(PR(I-1)+PI(I-1))*(PR(2)+PI(2))
          PR(I)=P-Q
          PI(I)=S-P-Q
30        CONTINUE
        DO 40 IT=0,N-2,2
          VR=FR(IT+1)
          VI=FI(IT+1)
          FR(IT+1)=VR+FR(IT+2)
          FI(IT+1)=VI+FI(IT+2)
          FR(IT+2)=VR-FR(IT+2)
          FI(IT+2)=VI-FI(IT+2)
40        CONTINUE
        M=N/2
        NV=2
        DO 70 L0=K-2,0,-1
          M=M/2
          NV=2*NV
          DO 60 IT=0,(M-1)*NV,NV
          DO 60 J=0,(NV/2)-1
            P=PR(M*J+1)*FR(IT+J+1+NV/2)
            Q=PI(M*J+1)*FI(IT+J+1+NV/2)
            S=PR(M*J+1)+PI(M*J+1)
            S=S*(FR(IT+J+1+NV/2)+FI(IT+J+1+NV/2))
            PODDR=P-Q
            PODDI=S-P-Q
            FR(IT+J+1+NV/2)=FR(IT+J+1)-PODDR
            FI(IT+J+1+NV/2)=FI(IT+J+1)-PODDI
            FR(IT+J+1)=FR(IT+J+1)+PODDR
            FI(IT+J+1)=FI(IT+J+1)+PODDI
60          CONTINUE
70        CONTINUE


        IF (L.NE.0) THEN
          DO 80 I=1,N
            FR(I)=FR(I)/N
            FI(I)=FI(I)/N
80          CONTINUE
        END IF
        IF (IL.NE.0) THEN
          DO 90 I=1,N
            PR(I)=SQRT(FR(I)*FR(I)+FI(I)*FI(I))
            PI(I)=ATAN(FI(I)/FR(I))*360.0/6.283185306
90          CONTINUE
        END IF
        RETURN
        END
这个是徐士良的快速傅里叶变换变换算法的源代码，能不能帮忙看一下划线部分到底是什么意思
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
如果你想拥有你从未有过的东西，那么你就要做你从未做过的事情。
1楼 2015-09-15 21:06:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cob

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
do Number xxx
xxxxx
Number continue

循环语句
赞一下(1人) 回复此楼
2楼2015-09-16 20:15:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

破碎东瓜

金虫 (初入文坛)
引用回帖:
2楼: Originally posted by cob at 2015-09-16 20:15:17
do Number xxx
xxxxx
Number continue

循环语句

我想问的是这个部分在整个FFT中起什么作用，那个代码我肯定是懂的，就是不懂为什么要这样编写
赞一下 回复此楼
如果你想拥有你从未有过的东西，那么你就要做你从未做过的事情。
3楼2015-09-17 08:11:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cob

木虫 (正式写手)
引用回帖:
3楼: Originally posted by 破碎东瓜 at 2015-09-17 08:11:06
我想问的是这个部分在整个FFT中起什么作用，那个代码我肯定是懂的，就是不懂为什么要这样编写...

作用应该是找出对应  PR(IS+1) / PI(IS+1) 赋值给 对应 的 FR(IT+1)/FI(IT+1)。
至于为啥这么写，书的公式应该有对应的解释。我没有这本书，所以没办法解释。
赞一下 回复此楼
4楼2015-09-17 19:48:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cob

木虫 (正式写手)
引用回帖:
4楼: Originally posted by cob at 2015-09-17 19:48:03
作用应该是找出对应  PR(IS+1) / PI(IS+1) 赋值给 对应 的 FR(IT+1)/FI(IT+1)。
至于为啥这么写，书的公式应该有对应的解释。我没有这本书，所以没办法解释。...

补充一点，
J=M/2 ， J 是等于 M/2 的整数部分。         
IS=2*IS+(M-2*J)， 所以 M-2*J ，M是奇数等于1，偶数等于0.
赞一下 回复此楼
5楼2015-09-17 20:02:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 破碎东瓜 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 8/400 2026-04-02 16:26 by psa1234
[考研] 一志愿同济大学323分（080500）求调剂 +6 yikeniu 2026-04-01 6/300 2026-04-02 14:19 by smileboy2006
[考研] 化学308分调剂 +14 你好明天你好 2026-03-30 15/750 2026-04-02 13:39 by 嘻嘻笑22
[考研] 085601 材料工程 313分 求调剂 +7 Ong3 2026-03-27 7/350 2026-04-02 11:29 by 3041
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 303求调剂 +8 DLkz1314. 2026-03-30 8/400 2026-04-02 09:08 by sanrepian
[考研] 求调剂，一志愿南京师范大学计算机专硕，初试373，六级通过， +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 296求调剂 +4 汪！？！ 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
[考研] 建环，能源，土木老师路过看一看！！！ +4 嘿嘿uu 2026-04-01 4/200 2026-04-01 20:42 by 无懈可击111
[考研] 286求调剂 +5 lim0922 2026-03-26 5/250 2026-04-01 19:08 by 客尔美德
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 322求调剂 +8 三水sss 2026-04-01 8/400 2026-04-01 10:19 by 唐沐儿
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 375求调剂 +7 雨夏整夜 2026-03-29 7/350 2026-03-31 18:52 by xhai2011
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 求化学调剂 +12 wulanna 2026-03-28 12/600 2026-03-31 16:38 by 690616278
[考研] 328求调剂 +8 嗯滴的基本都 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 298求调剂 +3 种圣赐 2026-03-29 3/150 2026-03-29 12:06 by longlotian
信息提示
关闭
请填处理意见
关闭