【文章內容簡介】 EN n=b*2 ELSE nu=int(log(float(b))/+factor_m) n=2**nu END IF n1=1+(nb)/2 n2=n1+b1 write(*,*)m1,m2,n1,n2,m,n pause end subroutine calculate_mn !************************************************************************* ! 子程序： INPUT_GRD ! 功能：讀取 grd文件中的數(shù)據(jù) ! 輸入?yún)?shù)說明： ! filename_obser：輸入文件名 ! m1,m2: x方向實際數(shù)據(jù)起點位置和終點位置點號 ! m:擴邊后數(shù)據(jù)終點位置點號 (起點位置點號為 1) ! n1,n2: y方向實際數(shù)據(jù)起點位置和終點位置點號 ! n:擴邊后數(shù)據(jù)終點位置點號 (起點位置點號為 1) ! 輸出參數(shù)說明： ! A：存放輸出數(shù)據(jù)的二維數(shù)組名 !************************************************************************* SUBROUTINE INPUT_GRD(A,input_file,m1,m2,n1,n2,m,n) character*(*)input_file integer m1,m2,n1,n2,m,n real xmin,xmax,ymin,ymax real A(1:m,1:n) real i,j,k do j=1,n,1 do i=1,m,1 A(i,j)=*1e10 enddo enddo Open(20,file=input_file,status=39。old39。) read(20,*) read(20,*) read(20,*)xmin,xmax read(20,*)ymin,ymax read(20,*) read(20,*) ((A(i,j),i=m1,m2),j=n1,n2) Close(20) END SUBROUTINE INPUT_GRD !************************************************************************* ! 子程序： expand_cos_2D ! 功能：二維擴邊子程序并為信號虛部賦值 ! 輸入?yún)?shù)說明： ! m1,m2: x方向實際數(shù)據(jù)起點位置和終點位置點號 ! m:擴邊后數(shù)據(jù)終點位置點號 (起點位置點號為 1) ! n1,n2: y方向實際數(shù)據(jù)起點位置和終點位置點號 ! n:擴邊后數(shù)據(jù)終點位置點號 (起點位置點號為 1) ! Ur:初始觀測面信號的實部 ! Ui：初始觀測面信號的虛部 ! 輸出參數(shù)說明： ! Ur:初始觀測面信號的實部 ! Ui：初始觀測面信號的虛部 !************************************************************************* subroutine expand_cos_2D(m1,m2,m,n1,n2,n,Ur,Ui) implicit none integer m,n,m1,m2,n1,n2 real Ur(1:m,1:n),Ui(1:m,1:n) real,allocatable::u(:),r(:) integer j,i,k allocate(u(1:m)) do j=n1,n2,1 do i=1,m,1 u(i)=Ur(i,j) enddo call expand_cos_1d(1,m1,m2,m,u(1)) do i=1,m,1 Ur(i,j)=u(i) enddo enddo deallocate(u) allocate(r(1:n)) do i=1,m,1 do j=1,n,1 r(j)=Ur(i,j) enddo call expand_cos_1d(1,n1,n2,n,r(1)) do j=1,n,1 Ur(i,j)=r(j) enddo enddo deallocate(r) do i=1,m do j=1,n Ui(i,j)=0 enddo enddo end subroutine expand_cos_2D !************************************************************************** ! 子程序： expand_cos_1d ! 功能：一維擴邊子程序 ! 輸入?yún)?shù)說明： ! n0,n3：擴邊后數(shù)據(jù)起點位置和終點位置 ! n1,n2：實際數(shù)據(jù)起點位置和終點位置 ! feild(i),(i=n1,n1+1,...,n2)：實際數(shù)據(jù) ! 輸出參數(shù)說明： ! field(i),(i=n0,...,n11)：擴邊后左邊的數(shù)據(jù) ! field(i),(i=n2+1,...,n3)：擴邊后右邊的數(shù)據(jù) !************************************************************************** Subroutine expand_cos_1d(n0,n1,n2,n3,Field) Real Field(n0:n3) pi= Field (n0)=(Field (n1)+Field (n2))/ Field (n3)=Field (n0) i1=n0 i2=n1 DO i=i1,i21,1 Field(i)=Field(i1)+cos(pi/*(i2i)/(i2i1))*(Field(i2)Field(i1)) End do i1=n2 i2=n3 DO i=i1+1,i2,1 Field(i)=Field(i1)+cos(pi/*(i2i)/(i2i1))*(Field(i2)Field(i1)) End do End subroutine expand_cos_1d !************************************************************************ ! 功能： FFT2 ! 功能：復數(shù)組 2D快速 Fourier變換 ! 輸入?yún)?shù)說明： ! m0,m3： x方向的起點和終點 ! n0,n3： y方向的起點和終點 ! field：輸入信號 (需要賦值給 Freal，實部 ) ! m,n： x,y方向擴邊后數(shù)據(jù)終點點號位置（起始點號為 1） ! NF: 正、反變換標志量 (1:反變換 。2:正變換 ) ! 輸出參數(shù)說明： ! Freal：信號的實部 ! Fimage：信號的虛部 (對于實信號而言，賦值為 0) ! 對應頻率分布說明： ! 數(shù)據(jù) Freal(m,n)和 Fimage(m,n)對應的頻率分布位置為 : ! m 方向： 0,1,.......,m/21,m/2,(m/21),......,1 ! n 方向： 0,1,.......,n/21,n/2,(n/21),......,1 !************************************************************************ SUBROUTINE FFT2(Freal,Fimage,m,n,nf) implicit none INTEGER m,n,nf REAL Freal(1:m,1:n),Fimage(1:m,1:n) real,ALLOCATABLE::Treal(:),Timage(:) integer nmmax,ierr,i,j nmmax=max(m,n) allocate(Treal(1:nmmax),Timage(1:nmmax),STAT=ierr) if() STOP DO i=1,m,1 IF () THEN do j=1