| ²é¿´: 1491 | »Ø¸´: 25 | |||
| µ±Ç°Ö÷ÌâÒѾ´æµµ¡£ | |||
| µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû | |||
LIHY3763½ð³æ (СÓÐÃûÆø)
|
[½»Á÷]
¡¾ÇóÖú¡¿2άFFTµÄFortran»òC³ÌÐò
|
||
| Çë´ó¼Ò°ï°³ÕÒ¸ö2άFFTµÄ³ÌÐò£¬¼±Ó𡣡Íò·Ö¸Ðл£¡ |
» ²ÂÄãϲ»¶
0856Çóµ÷¼Á
ÒѾÓÐ13È˻ظ´
0703»¯Ñ§
ÒѾÓÐ11È˻ظ´
329Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
321Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
0856²ÄÁÏ»¯¹¤µ÷¼Á ×Ü·Ö330
ÒѾÓÐ12È˻ظ´
309Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
070300»¯Ñ§354Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
330·ÖÇóµ÷¼Á
ÒѾÓÐ3È˻ظ´
299Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
Ò»Ö¾Ô¸±±¾©Àí¹¤´óѧ±¾¿Æ211²ÄÁϹ¤³Ì294Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
wangen994
ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)
½«¾ü
- ³ÌÐòÇ¿Ìû: 1
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ¹ó±ö: 3.348
- ½ð±Ò: 7797.4
- É¢½ð: 404
- ºì»¨: 5
- Ìû×Ó: 2007
- ÔÚÏß: 242Сʱ
- ³æºÅ: 649482
- ×¢²á: 2008-11-08
- ÐÔ±ð: GG
- רҵ: ÎïÀí
- ¹ÜϽ: ³ÌÐòÓïÑÔ
3Â¥2009-06-25 10:41:37
woshilsh
ÈÙÓþ°æÖ÷ (Ö°Òµ×÷¼Ò)
- Ó¦Öú: 43 (СѧÉú)
- ¹ó±ö: 10.33
- ½ð±Ò: 21115.1
- É¢½ð: 737
- ºì»¨: 22
- Ìû×Ó: 3777
- ÔÚÏß: 639.7Сʱ
- ³æºÅ: 293038
- ×¢²á: 2006-11-04
- רҵ: ¼ÆËã»úÓ¦Óü¼Êõ
- ¹ÜϽ: ·ÂտģÄâ

2Â¥2009-06-25 00:19:51
calos818
ľ³æ (ÖøÃûдÊÖ)
ÐÂÈ˲ËÄñ
- Ó¦Öú: 11 (СѧÉú)
- ½ð±Ò: 3125.7
- É¢½ð: 77
- ºì»¨: 3
- Ìû×Ó: 1201
- ÔÚÏß: 730.8Сʱ
- ³æºÅ: 698507
- ×¢²á: 2009-02-08
- רҵ: Äý¾Û̬ÎïÐÔ II £ºµç×ӽṹ
4Â¥2009-06-25 22:39:28
calos818
ľ³æ (ÖøÃûдÊÖ)
ÐÂÈ˲ËÄñ
- Ó¦Öú: 11 (СѧÉú)
- ½ð±Ò: 3125.7
- É¢½ð: 77
- ºì»¨: 3
- Ìû×Ó: 1201
- ÔÚÏß: 730.8Сʱ
- ³æºÅ: 698507
- ×¢²á: 2009-02-08
- רҵ: Äý¾Û̬ÎïÐÔ II £ºµç×ӽṹ
¡ï ¡ï ¡ï
sunxiao(½ð±Ò+2,VIP+0):thanks a lot 6-26 03:08
jove1782(½ð±Ò+1,VIP+0):½±Àø·¢´úÂë 9-27 10:45
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














»Ø¸´´ËÂ¥
10