24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1183  |  回复: 25
当前主题已经存档。

LIHY3763

金虫 (小有名气)

[交流] 【求助】2维FFT的Fortran或C程序

请大家帮俺找个2维FFT的程序,急用啊!万分感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主


LIHY3763(金币+1):谢谢参与
偶偶,你网上搜搜,一大堆哦
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
2楼2009-06-25 00:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangen994

荣誉版主 (著名写手)

将军


LIHY3763(金币+1):谢谢参与
你们应该是一个实验室的吧,ip都一样,哈哈
3楼2009-06-25 10:41:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟


LIHY3763(金币+1):谢谢参与
我也想找这个源程序,fortran
4楼2009-06-25 22:39:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟

★ ★ ★
sunxiao(金币+2,VIP+0):thanks a lot 6-26 03:08
jove1782(金币+1,VIP+0):奖励发代码 9-27 10:45
在Numerical Recipes 上找的,看看是否对你有用      
SUBROUTINE four1(data,nn,isign)
      INTEGER isign,nn
      REAL data(2*nn)
      INTEGER i,istep,j,m,mmax,n
      REAL tempi,tempr
      DOUBLE PRECISION theta,wi,wpi,wpr,wr,wtemp
      n=2*nn
      j=1
      do 11 i=1,n,2
        if(j.gt.i)then
          tempr=data(j)
          tempi=data(j+1)
          data(j)=data(i)
          data(j+1)=data(i+1)
          data(i)=tempr
          data(i+1)=tempi
        endif
        m=n/2
1       if ((m.ge.2).and.(j.gt.m)) then
          j=j-m
          m=m/2
        goto 1
        endif
        j=j+m
11    continue
      mmax=2
2     if (n.gt.mmax) then
        istep=2*mmax
        theta=6.28318530717959d0/(isign*mmax)
        wpr=-2.d0*sin(0.5d0*theta)**2
        wpi=sin(theta)
        wr=1.d0
        wi=0.d0
        do 13 m=1,mmax,2
          do 12 i=m,n,istep
            j=i+mmax
            tempr=sngl(wr)*data(j)-sngl(wi)*data(j+1)
            tempi=sngl(wr)*data(j+1)+sngl(wi)*data(j)
            data(j)=data(i)-tempr
            data(j+1)=data(i+1)-tempi
            data(i)=data(i)+tempr
            data(i+1)=data(i+1)+tempi
12        continue
          wtemp=wr
          wr=wr*wpr-wi*wpi+wr
          wi=wi*wpr+wtemp*wpi+wi
13      continue
        mmax=istep
      goto 2
      endif
      return
      END
C  (C) Copr. 1986-92 Numerical Recipes Software ,4-#.
5楼2009-06-25 22:41:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfq0730

木虫 (文坛精英)

顺天府府尹


LIHY3763(金币+1):谢谢参与
wangen994(金币+0,VIP+0):斑竹能看见,哈哈 7-5 08:28
在哪能看见IP?
...
6楼2009-07-04 17:40:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主


wangen994(金币+1,VIP+0):感谢你的讨论 7-5 08:28
matlab下有现成的源代码,有兴趣自己 参考下
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
7楼2009-07-04 17:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主


wangen994(金币+1,VIP+0):感谢你的讨论 7-5 08:28
调用格式y = fft2(A)
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
8楼2009-07-04 17:55:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianlangxingaa

铁杆木虫 (著名写手)

★ ★ ★
LIHY3763(金币+1):谢谢参与
wangen994(金币+2,VIP+0):感谢你的讨论 7-5 08:28
去看看fftw的主页,上面有各种fft的代码,但大部分都是C的。

Cheers
9楼2009-07-04 19:07:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

song_emu

金虫 (正式写手)


LIHY3763(金币+1):谢谢参与
sunxiao(金币+0,VIP+0):版主才能看到滴,呵呵 7-24 00:45
IP怎么看到的?
偶也想知道
10楼2009-07-05 23:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 LIHY3763 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见