 |
 |
 |
 |
|
Ãëàâíàÿ -> Ìíîãîñâÿçíûå ïîëîñêîâûå ñòðóêòóðû 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 [ 25 ] 26 27 28 KE=(VE*VE-ZY1)/ZY2 KO= (VO*VO-ZYl) /ZY2 Y1E=(YII+KE*Y12) /VE YIO=(YII + KO*Y12)/VO Y2E= (YI2+KE*Y22)/VE Y20= (Y12--KO*Y22)VO Ñ Ñ ÔÎÐÌÈÐÎÂÀÍÈÅ ÌÀÒÐÈÖÛ [Â] ÄËß ÐÀÑ×ÅÒÀ ÏÎÑÒÎßÍÍÛÕ ÈÍÒÅÃÐÈÐÎÂÀÍÈß Ñ â (1,1) = ED+Y1E*Z1 B(1,2) = ED--Y10*Z1 B(1,3) = ED-Y1E*Z1 B(1,4) = ED-Y10*Z1 B(2,1) = KE+Y2E*Z2 B(2,2) = K0+Y20*Z2 B(2,3) = KE-Y2E*Z2 B(2,4) = KO-Y20*Z2 VE=VE*L VO= VO L  (3,1) = (ED-Y1E*Z3)*CEXP (-VE) B(3,2)= (ED-Y10*Z3)*CEXP (-VO)  (3,3) = (ED+Y1E*Z3)*CEXP (VE) B(3.4)=(ED+YIO*Z3)*CEXP (VO) B(4,1)=(KE-Y2E*Z4)*CEXP (-VE) B(4,2)= (ÊÎ-Y20*Z4)*CEXP (-VO) B(4,3)=(KE+Y2E*Z4)*CEXP (VE) B(4,4)= (K04-Y20*Z4)*CEXP (VO) VE=VE/L VO= VO/L Ñ Ñ ÐÀÑ×ÅÒ ÏÎÑÒÎßÍÍÛÕ ÈÍÒÅÃÐÈÐÎÂÀÍÈß ÌÅÒÎÄÎÌ ÃÀÓÑÑÀ Ñ N=1 ! ×ÈÑËÎ ÏÐÀÂÛÕ ×ÀÑÒÅÉ (ÂÅÊÒÎÐ [Å1,Å2,ÅÇ,Å41) Ì=4 ! ×ÈÑËÎ ÓÐÀÂÍÅÍÈÉ ILL=0 ! ÏÀÐÀÌÅÒÐ CALL GAUSSC NL=2 ! ×ÈÑËÎ ÑÅ×ÅÍÈÉ -N-I DL=L/NL X=-DL DO 21= 1,3 ! ÖÈÊË ÏÎ ÄËÈÍÅ ÎÒ X=O...L, Ñ ØÀÃÎÌ DL X=X + DL è (1) =  (1,5)*ÑÅÕÐ (-VE*X) ! ÏÀÄÀÞÙÀß ×ÅÒÍÀß ÂÎËÍÀ è (2) =  (2,5)*ÑÅÕÐ (-VO X) ÏÀÄÀÞÙÀß ÍÅ×ÅÒÍÀß ÂÎËÍÀ è (3) =  ÃÇ,5)*ÑÅÕÐ (VE*X) ! ÎÒÐÀÆÅÍÍÀß ×ÅÒÍÀß ÂÎËÍÀ è (4) =  (4,5)*ÑÅÕÐ (VO*X) ! ÎÒÐÀÆÅÍÍÀß ÍÅ×ÅÒÍÀß ÂÎËÍÀ è (5) = è (1) + è (2) ] ÏÀÄÀÞÙÅÅ ÍÀÏÐßÆ.  ËÈÍÈÈ1 è (6) = è (3) + è (4) ! ÎÒÐÀÆÅÍÍÎÅ ÍÀÏÐßÆ.  ËÈÍÈÈ 1 è (7) = ÊÅ è (1) + ÊÎ*è (2) ! ÏÀÄÀÞÙÅÅ ÍÀÏÐßÆ.  ËÈÍÈÈ 2 è (8)=ÊÅ è (3)+ÊÎ*è (4) ! ÎÒÐÀÆÅÍÍÎÅ ÍÀÏÐßÆ.  ËÈÍÈÈ 2 è (9) = è (5)+è (6) ! ÍÀÏÐßÆ.  ÑÅ×ÅÍÈÈ X ËÈÍÈÈ 1 è 10) = è (7)+è (8) ! ÍÀÏÐßÆ.  ÑÅ×ÅÍÈÈ X ËÈÍÈÈ 2 U (11) = Y1E*U(1) 4-YlO*U (2) ! ÒÎÊ Â ÑÅ×ÅÍÈÈ X ËÈÍÈÈ 1 ! ÒÎÊ Â ÑÅ×ÅÍÈÈ X ËÈÍÈÈ 2 U(12) = -Y1E*U(3)-Y10*U(4) U(13) = Y2E*U(1) + Y20*U(2) U(14)=-Y2E*U(3)-Y20*U(4) è (15) = U (11)+U (12) è (16) = U (13) + U (14) ÀÏÏ = ÑÀÂÌÜ (II)) ! ÌÎÄÓËÜ ÍÀÏÐßÆÅÍÈß, ÒÎÊÀ PI (II)=(ATAN2 (AIMAG (U (II)), REAL (U (II))) ) 57.2958 ! ÔÀÇÛ PRINT 993,X PRINT 994 PRINT 995,A (1),A (2),A (3),A (4) PRINT 996,FI(1),FI(2),FI(3),FI(4) PRINT 999 PRINT 998,A(5),A(6),A(9) PRINT 996,FI(5),FI(6),FI(9) PRINT 997,A(7),A(8),A(10) PRINT 996,F1(7),FI(8),FI(10) PRINT 899 PRINT 998,A(11),A(12).A(15) PRINT 996,FI(11),FI(12),FI(15) PRINT 997,A(13),A(14),A(16) PRINT 996,F1(13),F1(14),F1(I6) 2 CONTINUE 1 CONTINUE U1E=REAL(Y1E) UIO=REAL (YIO) U2E = REAL (Y2E) U20=REAL (Y20) GE=CABS(VE) GO=CABS(VO) UE=REAL(KE) UO=REAL(KO) PRINT 600,GE,GO,UE,UO PRINT 700,U IE,U 10,U2E,U20 GO TO 1000 100 FORMAT(lX,4(2X,G125)/) 200 FORMAT(lX,F=,G125,MrU) 201 FORMAT(lX,----------, 9 ---------- ) 600 FORMAT(/, VE=,G12.5, VO=,G12.5, KE=,G12.5, 1 KO=,G125) 700 FORMAT (1X,Y1E=,G12.5,2X,Y10=,GI2.5,2X,Y2E=,GI2.5,2X, 2 Y20=,G12.5) 800 FORMATC Z1 = ,2G125, Z2=,2G12.5,/, 3 Z3=,2G12.5, Z4=,2G12.5/) 999 ÐÎÐÌÀÒ(1Õ,ÍÀÏÐßÆÅÍÈß-!.....ÏÀÄÀÞÙÈÅ.....!, 4 ...ÎÒÐÀÆÅÍÍÛÅ.....!.....ÏÎËÍÛÅ.....!) 899 FORMAT(lX,.....ÒÎÊÈ.....-!.....ÏÀÄÀÞÙÈÅ.....!, 4 .....ÎÒÐÀÆÅÍÍÛÅ.....!.....ÏÎËÍÛÅ.....!) 998 FORMAT(lX,  1 ËÈÍÈÈ:4(1Õ,0125)) 997 FORMAT(lX,BO 2 ËÈÍÈÈ:4 (1X,G125)) 996 FORMAT (IX, ÔÀÇÛ, ÃÐÀÄ.4 (1X,G12.5)) 995 FORMAT (1Õ,ÌÎÄÓËÜ,  4 (IX, G12.5)) 994 FORMAT (13Õ,!...ÏÀÄ. ×ÅÒÍÀß ...ÏÀÄ. ÍÅ×ÅÒ.. 5 I...OTP. ×ÅÒÍÀß! ...OTP. ÍÅ×ÅÒ.) 993 F0RMAT(1X,X=.G12.5.M) 99 STOP ÏÐÈÌÅÐ ÐÀÑ×ÅÒÀ ÏÎ ÏÐÎÃÐÀÌÌÅ CANON Ñ11= 429.17 L11 = 0.25840 Z 1 = 0.00000 Z3= 0.00000 El = 0.00000 Å3= 0.00000 Ñ12= 389.01 L12= 0.11620 0.00000 0.00000 0.00000 0.00000 Ñ22= 468.65 L22= 0.1 7850 Z2= 0.00000 Z4= 0.00000 E2= 0.00000 E4= 0.00000 - ÏÔ/Ì ÌÊÃÍ/Ì 0.00000 0.00000 0.00000 0.00000 F= 100.00 X= 0.00000 1...ÏÀÄ. ×ÅÒÍÀß!..ÏÀß. ÍÅ×ÅÒ!..OTP. ×ÅÒÍÀß..OTP. ÍÅ×ÅÒ! ÌÎÄÓËÜ,  0.00000 0.00000 0.00000 î 00000 ÔÀÇÛ, ÃÐÀÄ. 90.000 90.000 90.000 90.000 ÍÀÏÐßÆÅÍÈß-!..ÏÀÄÀÞÙÈÅ...!...ÎÒÐÀÆÅÍÍÛÅ. ...ÏÎËÍÛÅ...!  1 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ. ÂÎ 2 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 ..ÒÎÊÈ.. - ..ÏÀÄÀÞÙÈÅ ..!.ÎÒÐÀÆÅÍÍÛÅ...!...ÏÎËÍÛÅ...!  1 ËÈÍÈÈ: 0,00000 0.00000 0 00000 ÔÀÇÛ, ÃÐÀÄ. 90.000 90.000 90.000 ÂÎ 2 ËÈÍÈÈ: 0.00000 0.00000 0.00000 ÔÀÇÛ, ÃÐÀÄ. 90.000 90.000 90.000 Õ = 0.750Î0Å-01Ì !...ÏÀÄ. ×ÅÒÍÀß!. ÏÀÄ. ÍÅ×ÅÒ.!...OTP. ×ÅÒÍÀß!...OTP. ÍÅ×ÅÒ! ÌÎÄÓËÜ,  ÔÀÇÛ, ÃÐÀÄ. 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 ÍÀÏÐßÆÅÍÈß-!.-ÏÀÄÀÞÙÈÅ...!...ÎÒÐÀÆÅÍÍÛÅ...!...ÏÎËÍÛÅ...  1 ËÈÍÈÈ: 0.00000 0.00000 0.00000 ÔÀÇÛ, ÃÐÀÄ. ÂÎ 2 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ. 90.000 0.00000 90.000 90.000 0.00000 90.000 90.000 0.00000 90.000 ÒÎÊÈ. -!. ÏÀÄÀÞÙÈÅ..!...ÎÒÐÀÆÅÍÍÛÅ...!;..ÏÎËÍÛÅ. !  1 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ. ÂÎ 2 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ. Õ= 0.15000 Ì 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 ÍÅ×ÅÒ!..OTP. ×ÅÒÍÀß!..ÎÒÐ. ÍÅ×ÅÒ! 0.00000 90.000 0.00000 90.000 1...ÏÀÄ. ×ÅÒÍÀß!...ÏÀÄ ÌÎÄÓËÜ Â Î.ÎÎÎÎÎ 0.00000 ÔÀÇÛ, ÃÐÀÄ. 90.000 90.000 ÍÀÏÐßÆÅÍÈß-!...ÏÀÄÀÞÙÈÅ...!...ÎÒÐÀÆÅÍÍÛÅ...!...ÏÎËÍÛÅ...! Î 00000 0.00000 0.00000 90.0(É 90.000 90.000 0.00000 0.00000 0.00000 90.000 90.000 90.000 ...ÒÎÊÈ...-!...ÏÀÄÀÞÙÈÅ..!.ÎÒÐÀÆÅÍÍÛÅ...!..ÏÎËÍÛÅ...! 0.00000 90.000 0.00000  1 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ. ÂÎ 2 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ.  1 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ. ÂÎ 2 ËÈÍÈÈ: ÔÀÇÛ, ÃÐÀÄ 0.00000 90.000 0.00000 90.000 0.00000 90.000 0.00000 90.000 90.000 VE= 2.7463 YlE=0.81689E-02 VO 5.7942 Y10= 0,64256E-01 KE= 1.0115 K0= 0,42001 Y2E = 0 19449E-0I Y20= -0.63529E-01 ÏÐÎÃÐÀÌÌÀ ÐÀÑ×ÅÒÀ ÓÑÒÐÎÉÑÒ ÑÎ ÑÒÓÏÅÍ×ÀÒÎ-ÍÅÎÄÍÎÐÎÄÍÎÉ ÑÒÐÓÊÒÓÞÉ ANALIZ LOGICAL*I TIMER(8) COMPLEX XM,XD,XO,S(2,2),ZXR,ZX INTEGER KI(20),J,JI,J2,J3,IW,NNNN, I TIMOLD(2),TIMNEW(2),TIM(2) REAL F,FB,CD,LPD,LPA,KCB,P,FZ,DFZ, 1 C(I2),DLI(I2),LO(I2),WD(I0), 2 RO(I6),LOD(I6),ZG,ZH,K2(I5), 3 BP(4),PHI(4),SL(4),FF,RD(2), 4 CW(5),YCI (5),yCl2(5),YC2(5), 5 YLI(5),YL12(5),YL2(5) COMMON l/FP /F.FB.IW 2/UM /KI 3/COMP/XM,XD,XO 4/ZZ /ZG.ZH 5/SK /DLI 6/PS /J,JLJ2,J3 7/CH /KCB,P,FZ,DFZ 8/CK /WD 9/PD /RD,CD,LPD,LPA 1/ED /C,LO 2/DLN /RO.LOD 3/KK /Ê2 4/PP /CW,YCI,YCI2,YC2,YLI,YL12,YL2 5/BPC /BP ã 9999 NNNN=0 CALL READOP(AHAËÈÇ VI.4 CM) CALL CLOCK(l)-KII=I KI2=I K13=l I1 = KI (I) 12= ÊI (2) 13=K1 (3) IF I IF (NNNN.EQ.O) TYPE 44 (NNNN.EQ. 1) PRINT 44 DO 12 IW=I1,13,12 F=1\V/K2(10) FB=F/1000. IF (Kl(4).EQ.l) GO TO 3 K10=KI(10) DO 2 L=1,K10 DO 2 J=l,2 CALL CEK(S,L,J,K1 (5),K1(18)) CALL RP(S) ZX=ZXR(S(1,1),S(1,2),S(2.1),S(2,2)) IF (NNNN.EQ.O) TYPE 77,L,F,KCB,P,FZ,DFZ,ZX 16 15 14 CONTINUE GO TO 12 IF (K1(6).LT.2) GO TO 17 IF (Kl(6)-3) 14, 15, 16 K13=2 KI2=2 K1I = 2 DO II J= 1,2 DO II DO 11 DO 11 CALL CALL J1=1,K1I J2-=1,KI2 J3=1,K13 FSH(S.K1(6)) RP(S) TYPE 33,F,KCB,P,FZ,DFZ PRINT 33,F,KCB,P,FZ,DFZ (<CR>-HET, 1>-ÄÀ)) IF (NNNN.EQ.O) IF (NNNN.EQ.l) 11 CONTINUE 12 CONTINUE IF (NNNN.EQ.l) GO TO 9998 CALL CLOCK(2) TYPE 55 ACCEPT 66,NNNN IF (NNNN.EQ.l)GO TO 1 9998 TYPE 88 ACCEPT 66,NNNN IF (NNNN,EQ.O) GO TO 9999 STOP ÍÎÐÌÀËÜÍÎÅ ÇÀÂÅÐØÅÍÈÅ 33 FORMAT (IX,F5,2,4X,F5.2,3X,F5.2,4X,F6.I, 1 4X,F6.1,5X,F6.1,5X,F6.1) 44 FORMAT (F,9X,KCB,7X,P,8X,FZ, 1 9X,DFZ,9X,RZXR,9X,IZXR) 55 FORMAT ( ÍÀÏÅ×ÀÒÀÒÜ ÐÅÇÓËÜÒÀÒ? 66 FORMAT (II) 77 FORMAT (L=,I2,3X,F5,2,4X,F5.2,3X,F5.2,4X,F6.1, 1 4X,F6.1,5X,F6.1,5X,F6.1) 88 FORMAT ( ÏÐÎÄÎËÆÈÌ (<CR>) ÈËÈ ÇÀÊÎÍ×ÈÌ (ÂÂÅÄÈ U)?) FUNCTION ZXR(A,B,C,D) COMPLEX ZXR,A,B,C,D ZXR CSQRT((A B)/(C*D)) RETURN END SUBROUTINE RAW (A,B,N) COMPLEX A(N,N),B(N,N) DO 1 1=1,N DO 1 J=1,N 1 B(I,J) = A(I,J) RETURN END SUBROUTINE PR (A,B,C,N) COMPLEX A(N,N),B(N,N),C(N,N] DO 1 I=I,N DO 1 J=1,N C(l,J)=(0,0.) |
|
© 2024 Constanta-Kazan.ru
Òåë: 8(843)265-47-53, 8(843)265-47-52, Ôàêñ: 8(843)211-02-95 |