天线计算＿.pdf

麻烦哪位高手帮忙将天线计算的Fortran语言翻译成MATLAB语言，急需！！！谢谢

我帮你把文字抠出来了，你自己转换吧，仔细检查，文字在转换过过程中有一些错误，没办法
C************************************ C
C THIS PROGRAM USES HALLEN'S EQUATION ON SYMMETRICAL DIPOLES
C [EQU.(7-46)] AND COMPUTES THE CURRENT DISTRIBUTION,INPUT
C IMPEDANCE,AND RADIATION PATTERN OP A FINITE DIAMETER
C DIPOLE.EQUATION (7-46) IS SOLVED POR THE DIPOLE CURRC
ENT DISTRIBUTION USING MOMENT METHOD TECHNIQUES,WITH
C ENTIRE DOMAIN BASIS FUNCTIONS,AND POINT MATCHING,AS
C OUTLINED IN SECTION 7.S.2.C
C************************************ C
EXTERNAL GE,BMODE,POWER COMPLEX OP(5O,5O),EFRIL(50),ZIN
DIMENSION CRL(100),CIM(100),ZE(100),PWR(181),PHI(181)
1,PHID(181) COMMON RA,POINT,EL,N,Nl,NMODE,CRL,CIM,ZE DIMENSION
PIVOT(50),NPT(1),ICR(1),ICI(1)
DATA NPT/100/,ICR/'RRRR'/,ICI/'IIII'/ PI=3.1415927
C EACH DATA INPUT CARD HAS: NUMBER OF MODES TO BE USED (12),
C RADIUS (E10.3) AND LENGTH OF WIRE (F5.2) IN WAVELENGTHS.5 READ
(5,100,END=500) NMODE,RA,EL
100 FORMAT(12,E10.3,F5.2)


write(6,101)
NMODE.RA.EL
101 FORMAT(1H1,13,* MODES,RADIUS-1,E10.3,1 LENGTH-1,FS.2) C
A AND B ARE LIMITS OF NUMERICAL INTEGRATION
A=EL/2.B=EL/2.
DEL-EL*.05 KK-NMODE+1
DO 20 M-1,KK
POINT-EL*FLOAT(M-l)/FLOAT(NMODE)/2.C OP IS THE N*
H OPERATOR MATRIX OF "MUTUAL IMPEDANCES"
OP(M,KK)-(0.,1.)*(COS(2.*PI*POINT))/120./PI C
EFRIL IS THE EXCITATION MATRIX (CURRENT MATRIX AFTER
SOLNI
EFRIHM)-(0,,-1.)# (SIN(2.* PI* ABS (POINT)))/240 ./PI
C LIMITS ARE DIVIDED INTO 3 PARTS FOR SEPARATE
INTEGRATION C-POINT-DEL D-POINT+DEL IF(C.LT.A) C*A
IF(D.GT.B) D-B DO 10 N-l,NMODE C SEPARATE INTEGRATIONS
DONE FOR REAL AND IMAGINARY PARTS Nl-1
CALL QG32A (A,C,GE,ANSI) CALL QG32A (C,D,GE,ANS2) CALL
QG32A (D,B,GE,ANS3) ANSR-ANS1+ANS2+ANS3 Nl-2
CALL QG32A (A,C,GE,ANS1) CALL
QG32A (C,D,GE,ANS2) CALL QG32A
(D,B,GE,ANS3)
ANSI"AMS1+ANS2+ANS3 10
OP(M,N)-CMPLX(ANSR,ANSI) 20
CONTINUE
C LEQT1C IS AH IMSL ROUTINE FOR COMPLEX MATRIX EQN SOLN
CALL LEQT1C (OP,KK,50,EFRIL,1,50,0,PIVOT,IER)
C SET UP TWO (100 ELMT) ARRAYS HAVING REAL AND IMAG C
CURRENTS
DO 3 0 J=l,100
CRL(J) =0.
CIM(J) =0.
ZE(J)«FLOAT(J-l>*EL/198.
DO 30 I=l,NMODE N«I
CRL(J)iCRL(J)+REAL(EFRIL(I))*BMODE(ZE(J)) 30
CIM(J)=CIMfJJ+AIMAG(EFRIL(I))*BMODE(2E(J))
C CALCULATE THE INPUT IMPEDANCE
C 2IN IS GAP VOLTAGE (1 V) DIVIDED BY GAP CURRENT
ZIN-(1.,0.)/CMPLX(CRL(1),CIM(1))

write(6,104) ZIN 104
FORMAT (/,' THE INPUT IMPEDANCE IS (COMPLEX OHMS)' 1,2F7.2)
C
C PRINT OUT CURRENTS ALONG ONE HALF OF DIPOLE

write(6,107)
DO 40 J=l,100,3
CMAG-SQRT(CRL(J)** 2+(IM (J)** 2) 40

write   (6,108)
ZE(J),CRL(J),CIM(J),CMAG
C
C CALCULATE THE POWER PATTERN OF THE ANTENNA

write(6,110)
PMAX-0.DO 35 i=1,181
PHID(I)-l.0*FLOAT(I)-0.9999 PffI(Ii«PHID(I)*PI/180.
PWR(I)&#187;POWER(PHI< U)*60.*PI/SIN(PHI(I))
35 IF (PWR(I).GT.PMAX) PMAX-PWR(I) C PUT PATTERN IN
DB
DO 36 1-1,181
PWR(I)&#187;10.*ALOG10(PWR(I)/PMAX) IP
(PWR(I) .LT.-40.) PWR(I)-40.
36

write(6,111) PHID(I),PWR(I)
107 FORMAT(//,' CURRENT: Z',6X,'REAL',6X,'IMAGINARY',
44X,'MAGNITUDE')
108 FORMAT (Fl 2.,3,E11.3,2E13.3)
110 FORMAT(//,* POWER PATTERN: THETA POWER(DB)')
111 FORMAT(T17,F6.1,F11.2)
C REPEAT PROGRAM AS LONG AS DATA CARDS
REMAIN GO TO 5 500 STOP END
C FUNCTION SUBROUTINE NO.l FUNCTION
BMODE(Z)
C BHODE EVALUATES THE NTH MODE AT DISTANCE Z FROM
CENTER COMMON RA,POINT,￡L,N,Nl,NMODE
BMODE-SIN((FLOAT(N)* 3.141592 7/EL)*(EL/2.-ABS(Z)))
RETURN
END
C FUNCTION SUBROUTINE NO.2 FUNCTION GEIZ)
C GE PROVIDES THE KERNEL OF H&#212;LLENS EQN FOR INTEGRATION
C REAL OR IMAG PART CHOSEN BY PARAMETER N1 EQUAL TO 1 OR
2 EXTERNAL BMODE
COMMON RA,POINT,EL,N,N1,NMODE COMPLEX' G
PI-3.1415927
R"SQRT(RA**2+(POINT-Z)**2) K-2.*PI
G-CEXP((0.,-1..)*K*R) &#187;BMODE(Z)/4 ./PI/R IF (N1.EQ.2) GO
TO 10 GE-REAL(G) GO TO 20 10 GE*AIMAG(G) 20 RETURN
END
C FUNCTION SUBROUTINE Nò.3
FUNCTION POWER(THETA) C POWER CALCULATES THE
RADIATED POWER LEVEL AT ANGLE THETA C RADIANS TO THE
DIPOLE AXIS.CURRENT DISTRIBUTION GIVEN BY C ARRAYS CRL,
CIM,AT POSITIONS GIVEN BY ARRAY ZE DIMENSION
CRL(100),CIM(100),Zé(100) COMMON
RA,POINT,EL,N,N1,NMODE,CRL rCIM,ZE PI=3.1415927
PRR-0.0 ?RI*0.0 DO 10 1-1,99 K-I+l
PHAS&#171;PI*COS(THETA)*(ZE(I)+ZE(K))
PRR-PRR+(CRL(I)+CRL(K))*COS(PHAS)
10 PRI-PRI+(CIM(I)+CIM(K))*COS(PHAS)
POWER"((SIN(THETA))**3)*((EL/396.)** 2)*(PRR**2+PRI**2)
RETURN .END
C NUMERICAL INTEGRATION
SUBROUTINE
SUBROUTINE 0GJ 2 A (XL,XU,FCT,Y) C QG32A DOES
32 POINT GAUSSIAN QUADRATURE INTEGRATION OF C
FCT OVER THE LIMITS XL TO XU,AND RETURNS THE
ANSWE R C I N Y ,
A=-S*(XL+XU)
B=XU-XL
C=.4986319*B
Y=.003509305*(FCT(A+C)+FCT(A-C))C=.4928058*B
Y=Y+.008137197*(FCT(A+C)+FCT{A-C))
C=.4823811*B
Y=Y+.01269603*(FCT{A+C)+FCT(A-C))
C=.4674530*B
Y=Y+01713693*(FCT(A+C)+FCT(A-C))C=.4481606*B
Y=Y+.02141795*(FCT(A+C)+FCT(A-0)C=.4246838*B
Y=Y+.02549903*{FCT(A+C)+FCT{A-C))C=.3972419*B
Y*Y+.02934205*(FCT(A+C)+FCT(A-C))C=.3660911*B
Y=Y+.03291111*(FCT(A+C)+FCT(A-C))C=-.3315221*B
Y*Y+.03617290*(FCT(A+C)+FCT(A-C))C=.2938579*B
Y-Y+.03909695*(FCT(A+C)+FCT(A-C))C=.2534500*B
Y=IV.04165596*(FCT(A+C)+FCT(A-C))C=.2106756*B
Y=Y+.04382605*{FCT(A+C)+FCT(A-O)C=.1659343*B
V=Y+.04558694*(FCT(A+C)+FCT(A-C))C=.1196437*B
Y=Y+.04692220*(FCT(A+C)+FCT(A-C))
C=.07223598*9
Y=Y+.04781936*(FCT(A+C)+FCT(A-C))
C=.02415383*B
Y=B*(Y+.04837004*(FCT(A+C)+FCT(A-C)))
RETURN
END


write(rem,110) Phi*180/Pi
write(1,*)rem
write(rem,111)TimeStop-Wl/2
write(1,*) rem
write(1,*)'" "'
write(1,*)'" "'
write(1,*)MaxX,MinX,MaxY,MinY
write(1,*)XNum,YNum,Step close(1)
open(1,file='FDTD.El') close(1,status=delete')
open(2,file='FDTD.ER') closc(2,statu$=delete')
open(3 ,file='FDTD.HXI') closc(3,status='delete') open(4,file='FDTD.HXR')
close(4,status='delete') open(5,file='FDTD.HYI') close(5,status='delete')
open(6,file='FDTD.HYR') close(6,status='delete') open(7,file='FDTD.BI')
close(7,status='delete') open(8,file='FDTD.BR') close(8,status='delete')
end
*************************************:*
* Extrapolation of Radar Cross Section of Far Field**
* Based on FDTD Calculations**
* TYPE : main program**
* CALLS &#8226; Bessel Functions**
* LANGUAGE : Fortran77**** INPUT FILE NAME
*j ???????.bnd**
* OUTPUT FILE NAME : ???????.dat**
**************************************
*****************
complex Jz
complex Mx
complex My
complex Ctemp,ctempl
cotnplex E(-1250:1250,4)
complex H(-1250: 1250,4)
real WaveLength
integer WL
integer TEMFlag
real IncidentAngle
real Phi
real
StartAngle,EndAngic,Degree
Step character* 30 FileName
character* 1 Blank,Respond
logical FileExist
integer FileNameLength,Record
10000 write(*,*) char(7)
write(*,'(/////////////////////////)')
write(*,*)' '
write(*,*) Extrapolate Base on FD-TD Boundary to
Far Zone'
write(*,*)' to Caculate Radar Cross Section
write(*,*)' (Two-Dimension TM Mode >'
write(*,*)' '
write(*,*) Blank=' '
write(*,'(llX,34hThe scatterer data file name is,,$)') !
read(*,'(a30)')FileName Fi!eNameLength=index(FileName,
Blank)-1
inquire(file=FileNamed &#161;FileNameLength)//'.BND &#187;exist=FileExist)
if (FileExist) then
OPEN(1,file=FileNameC 11FileNameLength)//'.BND',+
access='direct',reel=8)
else
write(*,'(lx,15x,15hFilc not found.)')
write(*,(lx,15x,33hPress [Enter] to redo from start.)')
rcad(*,'(a1)') Respond
goto 10000
endif
read(1,rec=1) WaveLength
read(1,rec=2) WL
read(1,rec=3) TEMFlag
read(1,rec=4) Incident Angle
read C1,rec=5) Imin
read(1,rec=6) Imax
read(1,rec=7) jmin
read(1,rec=8) Jmax
read(1,rec=9) Itmin
read(1 frec=10)Itmax
read(1,rec=ll) jtmin
read(1 frec=12)Jtmax
read(1,rec=13)Iomin
read(1,rcc=14)Iomax
read(1,rec=15)Jomin
read(1,rcc=16)Jomax
read(1,rec=17)TimeStep
write(*,*)''
write(*,'(HXr15hKnown message :)')
write(*,'(13x,llhWaveLength=,el2.4,4H m=ti4,5H Grid)')-r
WaveLength,WL
if (TEMFlag.eq.1) then
write(*,'(13x,12hWaveMode: TM,5x,17hlncident angle=,
+f7.2,8h degrees)') Incident Angle
else
write(*,'(13x,12hWaveMode; TE,5x,17hIncident angle
=,+{7.2,8h degrees)') Incident Angle
end if
write(*/(13X,25hAbsorbing boundary is ;,4i6)')+
Iminrlrnax.jmin.Jmax

write(*() 3Xt25hConnective boundary is : ,4i6)') f-
Itmin,Itmax,Jtmin,Jtmax

write(*,'(13X,25hThe output boundary is ; ,4i6)')+
Iomin,Iomax,Jomin,Jomax
write(*,'(13X,14hTimeStop is : ,i6)')TimeStep
wriie(*,*)''
write(*,'CllX,19hMake choice below .)')
write(*,'(13x,34hStartDegree,EndDegree,DegreeStep:,$)')
read(*,*)StartAngle,EndAngle,DegreeStep
! begin
pi=3.1415927
rk=2*pi/iloat(WL) ! contant of propogation
K !
2=sqrt(pi*4.7/8.854e-12) ! wave inpedence of free space !
Rceord=17 do i=Iomin,Iomax Record=Record+1
read(1,rec=Record)E(i,l) ! MzO=Esz,JxO=Hsx
Record=Record+1
read (1,rec=Record) H (i,1)
Record=Record+1
read(l i rec-Record)Ctemp
E(i,3)"=-Ctemp ! MzO=-Esz
Record=Record+1
read(l,rec=Record) Ctempl
H(i,3)=-Ctempl end do
do i=jomin,Jomitx Record=Record+l
read(1,ree=Record)E(i,2) Record=Record 4"1
read (ì,rec=Record) H(i,2) Record=Reco rd
+1
read(1,rec=Record)Ctemp E(i,4)=-Ctemp
Record=Record+1
read(l,ree=Record) Ctempl H(i,4)=-Ctempl end
do close(l)
if CTEMFlag.eq.-1) then do L=Lomintiomax
E(i,l)=E(i,l) K(i,3)=E(i,3) HCi,1)=-
H(i,1) HCi,3)=-HO,3) end do
do i=jomin,jomax E(it2)=E(i,2)
E<i,4)=E(i,4) H(i,2)=-H(i,2)
H(l,4)=-H(i,4)
end do
Z=l,/Z
end if
OPEN(2 &#187;file=FiteName(1 (FileNameLength) //',dat') do dg=Start Angle)
End Angle ,DegreeStep
Jz=(0.,0,) Mx=(0.,0,) My=(0.,0,)

write(*,1) dg
formatClH+,10x,'Now calculate degree \f6.1)
Phi=dg*pi/180.
do I=Iomin,Iomax
Ctemp=cmplxC0.,(float(I)*cosCPhi)+Cjomin-0.5)*sirtCPhi))*rk)
Ctemp=cexp(Ctemp) ! exp(j*k*r);x=i ,j=jomin-.5
Mx=Mx+E
Citl)*
Ctemp
Jz=Jz+H(i,l)*Ctemp
Ctemp=cmplx(0.,Cfloat(I)*cosCPhi)+Cjomax+.5)*sin(Phi))*rk)
Cternp=cexp(Ctemp) ! exp(j*k*
r) ;x=i,j=jomin-.5
Mx=Mx+E(i,3)*Ctemp
Jz=Jz+H(i,3)*Ctemp end
do
do I=Jomin,Jomax
Ctemp=cmplx(0.,(fioati I)*sin(Phi)+Gomax+.5)*cosCPhi))*rk)
Ctemp=cexp(Ctemp) ! exp(j*k*r);
x=iomax+.5&#187;y=i
My=My+E(
i,2)*
Ctemp
Ji=Jz+H(i,2)*Ctemp
Ctemp=cmplx(0.,CfloatC 1)*sin(Phi)4"Ciomin-.5)*cosCPhi))*rk)
Ctemp=cexp(Ctemp) ! exp(j*k*r) ;x=iomin&#8482;.5,y=i
My=My+E(i,4)*
Ctemp
Jz=Jzl
H(i,4)*Ctemp
end
do
RCS=abs(Mx* sinCPhi)-My*cos(Fhi)+Z*Jz) RCS=.25*rk* res*
res*Wavelength/floatCWL)/Wavelength

write(1,*)dg,rcs RCS=10.*
loglO(RCS)

write(2,*)dg,RCS end do close(2)


write(*,*) CHAR(7) ! beep
end

本帖最后由 fong33 于 2010-12-10 06:15 编辑

kerbcurb 的Fortran Translation  {:7_1270:}

回复 kerbcurb 的帖子

谢谢你帮我那么多，但是我还是不会，我没学过fortran，看不懂，你能否帮我帮我把前几行和最后几行翻译一下，我要交个作业，麻烦你了！！！

我不会

这种行为，不知道该怎么说
别人交作业难道都是懂fortran的吗？不懂就去学习哈！

行为不好，自己动手，走捷径？

楼主的行为不厚道，难道高手是枪手吗

本帖最后由 fong33 于 2010-12-10 06:20 编辑

自己不动手，滥用高手当枪手，行为不可取啊！ 及早醒悟，翻开Fortran的书才是正道。

楼上几位说的好

恩，不会的可以自己学，不能过度依赖他人