實(shí)驗(yàn)四-位場邊緣識別程序設(shè)計(jì)實(shí)驗(yàn)-資料下載頁
2025-06-05 17:36本頁面
【正文】 xm u(i)=delx*(xx1) end do do j=1,n,1 do i=1,m,1 w(i,j)=sqrt(u(i)*u(i)+v(j)*v(j)) end do end do END SUBROUTINE WAVE2D !*********************************************************************************** ! 子程序： factor ! 功能：計(jì)算 x,y,z方向?qū)?shù)的實(shí)部和虛部 ! ! 輸入?yún)?shù)說明： ! m: 點(diǎn)數(shù)（ M 必須為 2 的冪次方） ! n: 線數(shù)（ N 必須為 2 的冪次方） ! field_re:初始觀測面信號的實(shí)部 ! field_im：初始觀測面信號的虛部 ! W(m,n)：徑向圓頻率 ! 輸出參數(shù)說明： ! Px_re:x方向?qū)?shù)信號的實(shí)部 ! Px_im:x方向?qū)?shù)信號的虛部 ! Py_re:y方向?qū)?shù)信號的實(shí)部 ! Py_im:y方向 導(dǎo)數(shù)信號的虛部 ! Pz_re:z方向?qū)?shù)信號的實(shí)部 ! Pz_im:z方向?qū)?shù)信號的虛部 !*********************************************************************************** subroutine factor(m,n,field_re,field_im,u,v,w,px_re,px_im,py_re,py_im,pz_re,pz_im) implicit none integer m,n real field_re(1:m,1:n),field_im(1:m,1:n) real px_re(1:m,1:n),px_im(1:m,1:n),py_re(1:m,1:n),py_im(1:m,1:n),pz_re(1:m,1:n),pz_im(1:m,1:n) real U(1:m),V(1:n),W(1:m,1:n) integer i,j real pi pi= do i=1,m,1 do j=1,n,1 px_re(i,j)=field_im(i,j)*u(i)*(1)*pi* px_im(i,j)=field_re(i,j)*u(i)*pi* py_re(i,j)=field_im(i,j)*v(j)*(1)*pi* py_im(i,j)=field_re(i,j)*v(j)*pi* pz_re(i,j)=field_re(i,j)*w(i,j)*pi* pz_im(i,j)=field_im(i,j)*w(i,j)*pi* end do end do end subroutine factor !*************************************************************************** ! 子程序： deviration ! 功能：計(jì)算異常的水平導(dǎo)數(shù)、垂向?qū)?shù)及總導(dǎo)數(shù) ! ! 輸入?yún)?shù)說明： ! m: 點(diǎn)數(shù)（ M 必須為 2 的冪次方） ! n: 線數(shù)（ N 必須為 2 的冪次方 ） ! Px_re:x方向?qū)?shù)信號的實(shí)部 ! Px_im:x方向?qū)?shù)信號的虛部 ! Py_re:y方向?qū)?shù)信號的實(shí)部 ! Py_im:y方向?qū)?shù)信號的虛部 ! Pz_re:z方向?qū)?shù)信號的實(shí)部 ! Pz_im:z方向?qū)?shù)信號的虛部 ! 輸出參數(shù)說明： ! field_vdr:對場值作水平導(dǎo)數(shù)的結(jié)果 ! field_thdr:對場值作垂向?qū)?shù)的結(jié)果 ! field_asm:對場值作總導(dǎo)數(shù)的結(jié)果 !************************************************************************** subroutine deviration(m,n,px_re,py_re,pz_re,field_vdr,field_thdr,field_asm) implicit none integer m,n real px_re(1:m,1:n),py_re(1:m,1:n),pz_re(1:m,1:n) real field_vdr(1:m,1:n),field_thdr(1:m,1:n),field_asm(1:m,1:n) integer i,j do i=1,m,1 do j=1,n,1 field_vdr(i,j)=pz_re(i,j) field_thdr(i,j)=sqrt(px_re(i,j)**2+py_re(i,j)**2) field_asm(i,j)=sqrt(field_vdr(i,j)**2+field_thdr(i,j)**2) end do end do end subroutine deviration !*************************************************************************** ! 子程序： OUTPUT_GRD ! 功能：輸出結(jié)果 ! ! 輸入?yún)?shù)說明： ! A：存放輸出數(shù)據(jù)的二維數(shù)組 ! m: 點(diǎn)數(shù)（ M 必須為 2 的冪次方） ! n: 線數(shù)（ N 必須為 2 的冪次方） ! m1,m2：實(shí)際數(shù)據(jù) x方向起點(diǎn)位置和終點(diǎn)位置 ! n1,n2：實(shí)際數(shù)據(jù) y方向起點(diǎn)位置和終點(diǎn)位置 ! 輸出參數(shù)說明： ! x_min,x_max： x的最小值和最大值 ! y_min,y_max： y的最小值和最大值 ! filename:存放輸出結(jié)果的文件名 ! field_vdr:對場值作水平導(dǎo)數(shù)的結(jié)果 ! field_thdr:對場值作垂向?qū)?shù)的結(jié)果 ! field_asm:場值的解析信號振幅的結(jié)果 !************************************************************************** SUBROUTINE OUTPUT_GRD(A,filename,m1,m2,n1,n2,m,n,eigval,xmin,xmax,ymin,ymax) real eigval integer m1,m2,n1,n2,m,n real A(1:m,1:n) character*(*) filename real amin,amax amin=HUGE(amin) amax=HUGE(amax) write(*,*)m1,m2,n1,n2,m,n pause eigval=*1e5 DO j=n1,n2,1 do i=m1,m2,1 if(ABS(A(i,j)).) then amin=MIN(amin,A(I,j)) amax=MAX(amax,A(I,j)) end if end do END DO OPEN(20,file=filename,status=39。unknown39。,form=39。formatted39。) write (20,39。(A)39。)39。DSAA39。 write (20,*)m2m1+1,n2n1+1 write (20,*)xmin,xmax write (20,*)ymin,ymax write(20,*)amin,amax Do j=n1,n2,1 write(20,*) (A(i,j),i=m1,m2) End do Close(20) END subroutine OUTPUT_GRD
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1