Question

fortran语言翻译成matlab语言

005DIMENSIONV(60,40)
010DATA VC,ER,E0/3.0E+08,9.6,8.854E-12/
016Z1=44.0
020ALFA=1.75
025I1=31
035I2=51
045DO 10 I=1,I2
050DO 10 J=1,31
06510V(I,J)=0
075DO 20 I=21,I1
08520V(I,11)=1.0
105I2M1=I2-1
10718K=0
110ER=9.6
120DO 211 ITE=1,2
135IF(ITE.EQ.2)ER=1.0
1408K=K+1
145DO 100 I=2,I2M1
155DO 100 J=2,30
165IF (V(I,J).EQ.1.0)GOTO 100
175IF (J-11) 60,40,60
18540VNEW=((V(I,J+1)+ER*V(I,J-1))*2.0/(ER+1.0)+V(I-1,J)+V(I+1,J))/4.0
205GOTO 70
21560VNEW=(V(I+1,J)+V(I-1,J)+V(I,J+1)+V(I,J-1))*0.25
22570VNEW=ALFA*VNEW+(1.0-ALFA)*V(I,J)
255V(I,J)=VNEW
265100 CONTINUE
267IF(K.LT.5)GOTO 8
305WE=0.
315DO 150 I=1,I2M1
325DO 150 J=1,30
335ENERGY=(V(I,J)-V(I+1,J+1))**2+(V(I+1,J)-V(I,J+1))**2
345IF(J.LT.11) ENERGY=ENERGY*ER
350ENERGY=ENERGY*E0/4.0
355150 WE=WE+ENERGY
365C0=WE*2.0
375IF(ITE.EQ.1)GOTO 140
385Z0=1.0/(VC*SQRT(C0*CER))
395GOTO 211
405140 CER=C0
415211 CONTINUE
417DEL=ABS(Z0-Z1)/Z0
420Z1=Z0
422IF(DEL.GT.1.0E-3)GOTO 18
425WRITE(6,17)C0,CER,Z0
43517FORMAT(1X,3HC0=,E14.3,1X,4HCER=,E14.3,1X,3HZ0=,F12.3)
445STOP
455END
程序解释：
句45~85给导体边界赋值和给内部网格节点赋初值。
句120~135当时是计算介质存在时的电容，如则是求介质不存在时的电容。
句145~267计算线内部节点的电位，其中185句是根据式（1-29）求介质交界面上的电位，每迭代5次计算电位后，跟着求线的电容和特性阻抗。
句305~355求场区电场的总储能。
句385求特性阻抗。
句422检查特性阻抗的计算是否已收敛，未收敛则重复以上的计算；否则打印输出。

Answer 1

不知道