FAZAMI=AMIN1 (FAZAMI, FAZAM-KRFAZF) KCB1MI=AMIN1 (KCBIMI, KRKCB-KCBIM) KCB2MI=AMIN1 (KCB2MI, KRKCB-KCB2M)

IF (KI-KRK.GE.VK1.AND.K2-KRK.GE.VK2.AND.KRKCB-KCBl.GE. *VKCB LAND.KRKCB-KCB2.GE.VKCB2.AND.FAZA-KRFAZF.QE. VFAZA) GO T02

PRINT3,F,K1,K2,FAZA,KCB1,KCB2,KRK,KRFAZF,KRKCB,K1M,K2M,

FAZAM,KCB1M,KCB2M,PR1,PR2,PR3,PR4,PR5,PR6,PR7,PR8,PR9,PR10 PRINT 9,PR11,PR12,PR13.PR14,PR15,PR16,PR17,PR18,PR19,PR20,PR21,

PR22,PR23,PR24 PR25

3 FORMAT (2X,-ДОПУСКИ ПРЕВЫШЕНЫ72Х,Т=.Е5.0.МГЦ'./2Х,

К1 = ,Е7.2,ДБ,ЗХ,К2= ,F7.2,flB.3X,FAZA=.El 1.4,ТРАД,5Х, KCBl = ,F5.2,3X,KCB2=,F5.2/2X,KRK=,F7.2,fl5M5X,KRFAZA==-, *F5. l,TPAД ,ЗX,KRKCB=F5.1,/2X, MAKCИMAЛЬHO ОПАСНЫЕ *ОТКЛОНЕНИЯ72Х,К1М= .Е7.2,ДБ,2Х,К2М= ,F7.2,flB,FAZAM,F5.1, ГРАД,4Х;КСВ1М=,Е5.2,2Х,КСВ2М=,Е5.2/2Х,ЧУВСТВИТЕЛЬ-

НОСТИ KI ПО E2,H1.H2,H3.X72X,PRK1E2=,E11.4 ДБ,2X,PRK1H1 =

Е11.4,ДБ/ММ,2Х,РРК1Н2=,Е11.4,ДБ/ММ72Х,РРК1НЗ=,Е11.4,

flB/MM,2X,PRK 1Х= ,Е 11.4,ДБ/ММ72Х;ЧУВСТВИТЕЛЬНОСТИ

K2nOE2,Hl,H2,E11.4,flB/MM72X,PRK2H3=,E11.4.flB/MM2X -

PRK2X=,E11.4,flB/MM)

9 FORMAT (2Х,ЧУВСТВИТЕЛЬНОСТИ KCBl ПО Е2,Н1 ,Н2,НЗ,Х72Х,

PRKCB1E2=:,E11.4,2X,PRKC1H1 = ,E11.4,1/MM,2X,PRKCB1H2=.E11.4,

1/mm72X,PRKCB1H3=,E11.4,1/MM,2X,PRKCB1X=,E11.4,1/MM72X,

ЧУВСТВИТЕЛЬНОСТИ KCB2nOE2,Hl.H2,X/2X,PRKCB2E2=, * E11.4,2X,PRKCB2H1 = ,E11.4,1/MM,2X,PRKCB2H2= E11.4

1/MM/2X,PRKCB2H3=.E11 .4,l/MM,2X,PRKCB2X =.El 1.4, 1/MM,

/2Х,ЧУВСТВИТЕЛЬНОСТИ FAZA ПО E2,H1 ,H2,H3,X72X,PRFAZAE2=

El 1.4,ГPAД ,2X,PRFAZAHЗ= .El 1.4,ГPAД/MM,2X,PRFAZAX= El 1.4,

ГРАД/ММ) GO TO 4

2PRINT 5,F K1,K2,FAZA.KCB1,KCB2

5 FORMAT (2X,F= .FS.O.SX.K 1 = ,F7.2,ДБ,ЗХ,К2= ,F7.2,flB,3X,FAZA,

Е5.1,ГРАД,5Х,КСВ1 = ,F5.23X,KCB2 = ,F5.2) 4 F=F + DF

IF (F.LE.FKl GO TO 6

PRINT 7,K1M1,K2M1,FAZAMI,KCB1M1,KCB2MI 7 ФОРМАТ (2Х,МИНИМАЛЬНЫЕ ЗАПАСЫ ПО ВСЕМ ЧАСТОТАМ

/2X,KlMI=,F7.2,flB,/2X,K2MI=,F7.2,flB/2X,FAZAMI=,F5.1,

ГРАД,/2Х,КСВ 1М1= ,F5.3,/2X,KCB2MI= \F5.3) STOP END

SUBROUTINE ANALIZ (Cl 1 ,C 12,C22,L1 l,L 12.L22,R11,R12,R22,G11,G12,

G22,F,L,CN,CK,R,LV,LD,CP,D1,D2,D3,D4,D5.D6,NS) REAL L11,L12,L22,L,KCB(2),K(2),FAZA(2),LV,LD COMPLEX AP (4,4),B(4,4),AA (4,4),A1 (2,2),A(2,2),Y1,T11

CALL BAZA (Cl 1,C12,C22,L11,L12,L22,R11,R12,R22,G11,G12,G22,F,L.B,0) DO 1 1=1,4 I

DO 1 J= 1,4 AP (I,J)= (0.,0.) 1 IF (I.EQ.J.)AP (I,J)= (l.,0.> PI= 3,14593

DO 2 1=1,2

IF (I.EQ.I) C=CN

IF (I.EQ.2) C=CK

Yl=l./ (R4. (0.,1.) (2. PI F* LV -1./ (2.E-6* PI F* C* ))) +

(0.,1) (2.E-6* PI* F* CP-1./ (2. F* LD* PI))) AP (4,2) = Y1

CALL UM3 (B.AP.B.AA )

Al (1,1) = AA (1,1)-AA (1,2) * AA (4,1)/AA (4,2) Al (1,2) = AA (1,3)-AA (1,2) * AA (4,3)/AA (4,2) Al (2,1) = AA (3,1)-AA(3,2) * AA (4,1) /АА (4,2) Al (2,2) = AA (3,3)-AA (3,2) * AA (4,3) /АА (4,2) CALL US (Al.A.NS)

Tll= (A (1,1)+A (1,2) /50. + А (2,1) 50. + A (2,2)) /2 К (I) = -20.* ALOGIO (CABS(Tll)) FAZA (I) = ATAN2 (AIMAG(T11 ),REAL (Til)) * 57.2958 S11 = CABS ((Tll-A (2.1)*50.-A (2.2))/Tll) 2 KCB (I)=(1.+ S11) / (l.-Sll) D1=K(1) D2=K (2) D3=FAZA (1) D4=FAZA (2) D5= KCB (1) D6= KCB (2) RETURN END

SUBROUTINE BAZA (Cl 1,C12,C22,L11,L12.L22,R1,R12,R2,G1,G12,G2,F, L,B,IG;

REAL L, L11,L12,L22,L1,L2,L3

COMPLEX B(4,4).ZYl,ZY2,ZY3,ZY4.GAMMAO, GAMMAE.BD.TE TO.SE.SO

* CE,C0,KE,K0,Y1E,Y10,Y2E,Y20,W,U,Z11,Z12 Z22 Yl 1 Y12 Y22 WL=6.28319 *F

WC=WL* l.E-6

Z11 = CMPLX (Rl.WL* Lll)

Z12=CMPLX (R12,WL*L12)

Z22=CMPLX (R2,WL* L22)

Y11 = CMPLX (G1,WC*C11)

Y12=CMPLX (G12,-WC * C12)

Y22=CMPLX(G2,WC * C22)

ZY1 = Z11 Yl 1 + Z12 * Y12

ZY2=Z11 Y12+Z12* Y22

ZY3=Z12* Y11 + Z22* Y12

ZY4=Z12 Y12+Z22 * Y22

BD=CSQRT ((ZYl-ZY4)**2+4* ZY2* ZY3)

GAMMAE=CSQRT ((ZY1 + ZY4 + BD) /2.)

GAMMAO=CSQRT ((ZY1 + ZY4-BD) /2.)

TE=GAMMAE*L

TO=GAMMAO *L

U=CEXP IE)

W=CEXP (-ТЕ)

SE= (U-W) / (2.,0.)

CE=(U+W) / (2.,0.)

U=CEXP (TO)

W=CEXP (-TO)

SO= (U-W) / (2.,0.)

СО= (lJ + W)/(2..0.) KE=(GAMMAE**2-ZY1) /ZV2 КО= (OAMMAO ** 2-ZYl) /ZY2 Y1E= (Yll+KE*Y12)/GAMMAE Y10= (Yll + KO*Y12)/GAMMAO Y2E= Y12+KE*Y22)/GAMMAE Y20= (YI2+KO*Y22)/GAMMAO и=\\ОЛ:\: Y1E*Y20 W=KE KO

IF (IG.Ey. I) PRINT IF, ТЕ, TO

1 FORMAT (IX Т=,Р6.1,ЗХ,ТЕ=,2Е10.3,ЗХ,ТО=,2Е1О.З) В (l,2)=(CE-CO) /W

В (4,3)=B (1,2)

В (l,l)=(-KO*CE+KE*CO) /W В (3,3)=B (1,1)

B(l,3)= (-Y20*SE+Y2E*SO)/U B(l,4)= (Y10*SE-YlE*SO) /U В (2,3) = B (1,4)

В (2,1)=( -CE+CO) /W* KE* KO В (3,4) = B (2.1)

В (2,2) = (KE * CE-KO* CO) /W В (4,4) = B (2.2)

B(2,4)= (KE*YlO*SE-KO*YlE*SO)/U B(3,I)=(-K0*Y1E*SE + KE*Y10*S0)/W B(3,2)= (YlE*SE~Y10*S0)/W B(4,1) = B(3,2)

B(4,2)= (Y2E*SE-Y20*S0)/W

RETURN

SUBROUTINE UM3 (A,B,C,D)

COMPLEX A (4,4),B(4,4) E(4,4),C(4,4),D(4,4),S

DO 1 1=1.N DO 1 J=1,N S= (0.,0.) DO 2 K= 1,N

2 S=S+A (LK) * В (K.J) 1 E (I,J) = S

DO 3 1=1,N DO 3 J=1,N S= (0.,0.) DO 4 K= l.N 4 S=S + E (I,K)*C(K,J)

3 D (LJ) = S RETURN END

SUBROUTINE US (A.B.N )

COMPLEX E,F.G.H.D,L1,M1.M2.T.U,S,A (2,2),B(2,2) E=A (1,1) F=A (1,2) G=A (2,1) H=A (2,2)

D=CSQRT ((E-H) **2+ (4.,0.) * F * G)

Ll= (E + H + D)/(2.,0.)

Ml= LI * N

M2= (Ll-D) ** N

T= (M1 + M2) / (2..0.)

U= (Ml -M2) /D

S=(E-H) * и/ (2.,0.)

В (1,1) = T+S

B(1,2) = F*U

B(2,1) = G*U

В (2,2) = T-S

RETURN

SUBROUTINE CONl (E2,H1,H2,H3,W,S,WP,LM,N,NS,CN,CK,R,LV,LD,

* CP,TGD,D1,D2,D3,D4,D5,D6,F) REAL L11,L12,L22,LM,L,LV,LD

DOUBLE PRECISION ZZ,ZM,ZZH1,ZMH1,ZZH2,ZMH2,ZZH3,ZMH3,ZZW,

* ZMW,ZZS,ZMS,ZZWP,ZMWP,ZZLM,ZMLM,ZZT,ZMT,ZZ1,ZM1

4 FORMAT (1X,7F5.2.1X,7F5.2,1XI5,1XF5.0) 1 FORMAT (5F5.2/9F5.2/2F5.2,I5,F5.2)

El=l. E3= 1. E4= 1. E5=E1 E6=E2 E7= E3 E8= E4

CALL CLl (E1,E2.E3,E4.H1.H2,H3.W,S,WP.LM,N,T,ZZ.ZM,C11,C12.C22. * L11,L12,L22,L,CX12.CY12) El= 1. E2= 1. E3=l. E4= 1.

CALL CLl (E1,E2,E3,E4,H1,H2.H3,W,S,WP,LM,N,T,ZZ1,ZM1,F1,F2,F3,F4,

* F5,F6,F7,0,F8,F9)

Н1=Н1 + ГЕ-Э

CALL CLl (E1,E2,E3,E4,H1,H2,H3,W,S,WP,LM,N,T,ZZH1,ZMH1,F1,F2,F3

* F4 F5,F6 F7,0,F8,F9) Hl = Hl-l.E-3 H2= H2+1.E~3

CALL CLl (E1,E2,E3.E4,H1,H2,H3,W,S,WP,LM,N,TZZH2,ZMH2.F1,F2F3

* F4,F5 F6 F7,0 F8,F9) H2=H2-l.E-3 H3= H3+1.E-3

CALL CLl (E1,E2,E3,E4,H1,H2,H3,W,S,WP,LM,N,T,ZZH3,ZMH3 F1.F2.F3.

* F4,F5,F6,F7,0,F8,F9) H3= H3-IE-3 W=W+l.E-3

CALL CLl (E1,E2,E,3,E4,H1,H2,H3,W,S,WP,LM,N,T,ZZW,ZMW Fl F2 F3

* F4,F5,F6,F7,0,F8,F9) W=W-l.E-3

S=S+1.E-3

CALL CLI (ELE2,E3,E4,HLH2,H3,W,S,WP,LM,N,T,ZZS,ZMS,F1,F2,F3,F4,

F5,F6,F7,0,F8,F9) S= S-LE-3 WP=WP+LE-3

CALL CLI (ELE2,E3,E4,H1,H2,H3,W,S,WP,LM.N,T,ZZWP,ZMWP,FLF2,

F3,F4,F5,F6,F7,0,F8,F9) WP=WP-LE-3

CALLMil,E2,E3,E4,Hl,H2.H3,W,S,WP,LM,N,T,ZZLM,ZMLM,Fl,F2,

F3,F4,F5,F6.F7,0,F8,F9) LM=LM-LE-3 T=T+LE-3

CALL CLI (E1,E2,E3,E4,HLH2,H3,W,S,WP,LM,N,T,ZZT,ZMT,FLF2,F3,F4,

F5,F6,F7,0,F8,F9) T=T-LE-3 RS=2.607E-4*SQRT (F)

ALFAZ=L153E4*RS*(2.*ZZH2-2. ZZW+2.*ZZS-2.*ZZWP+ 2. ZZLM-2.*ZZT)/ZZ

ALFAM= L153E4*RS*(2.*ZMH2-2.*ZMW+2.*ZMS-2.*ZMWP+ 2.

ZMLM-2*ZMT)/ZM PRZH1=LE6 (ZZH1-ZZ1)*ZZ/ZZ1 PRMH1 = LE6*(ZMH1-ZM1)*ZM/ZM1 PRZH2= LE6*(ZZH2-ZZ1)*ZZ/ZZ1 PRMH2= LE6*(ZMH2-ZM1) ZM/ZM1 PRZH3= 1 .E6* (ZZH3-ZZI) ZZ/ZZ1 PRMH3= 1 .E6* (ZMH3-ZM 1) ZM/ZM1 PRZW= LE6*(ZZW-ZZ1)*ZZ/ZZ1 PRMW= LE6*(ZMW-ZM1)*ZM/ZM1 PRZS=LE6*(ZZS-ZZ1)*ZZ/ZZ1 PRMS= LE6*(ZMS-ZM1)*ZM/ZM1 PRZWP= LE6 (ZZWP-ZZ1 ) ZZ/ZZ1 PRMWP= LE6*(ZMWP-ZM1)*ZM/ZM1 PRZLM= LE6* (ZZLM-ZZ1)ZZ/ZZ1 PRMLM= LE6* (ZMLM-ZM 1)ZM/ZM1 PRZT=LE6*(ZZT-ZZ1)*ZZ/ZZI PRMT= LE6* (ZMT-ZM 1) ZM/ZM1

E1 = E5 E2=E6 E3= E7 E4= E8

3 CALL MODA (Cl 1,C12,C22,L11,L12,L22,ZZ,ZM,ALFAZ,ALFAM,0,R11, R12 R22)

2 FORMAT (1Х,ИСХ0ДНЫЕ ДАННЫЕ71Х,Е1 = ,Р5.2,2Х,Н1 = ,Е5.2,2Х,

W= Т5.2,2Х,Т= ,F5.2/1X,E2= ,F5.2,2X,H2= ,F5.2,2X,S= ,F5.2,2X,N=,

12/ 1X,E3= ,F5.2,2X.H3= ,F5.2,2X,WP= ,F5.2,2X.F= ,F7.2/lX,E4=

F5.2,2X,LM=,F5.2,2X,T=-,F5.2/lX,norOHHbIE ПАРАМЕТРЫ'/1Х. C11 = ,F12.8,2X,LI1 = ,F12.8,/1X,C12=,F12.8,2X,LI2= ,F12.8./1X, C22=,F12.8,2X,L22=,F12.8/lX,3ATyXAHHE7lX,ALFAE=,F11.8,

ДБ/M,2X,ALFAO=.F11.8,ДБ/M71X,BOЛHOBЫE СОПРОТИВЛЕ-

НИЯ71Х,гЕ=,Р8.4,2Х,гО=,Р8.4/1Х.ПРОИЗВОДНЫЕ ВОЛНОВЫХ СОПРОТИВЛЕНИЙ В OM/M71X,PREHl = ,E11.4,2X,PROHl = ,

E11.4/lX,PREH2=,En.4,2X,PROH2=,EII.4/lX/PREH3=,Ell.4,2X,

PROH3=,En.4/IX,PREW=,EII.4,3X,PROW=,EII.4/lXPRES=-

ElI.4.3X,PROS=,EII.4/lX,PREWP=,EII.4,2X,PROWP= Ell 4/IX

PRELM=,EI1.4,2X,PROLM=,EII.4/IX,PRET=,ElI.4,3X,=PROT=.

El 1.4)

L= L/2.E3 PI=3.141593

011= (CX12+CYl2)*2.*PI*F TGD*l.E-6

G12=C12*2.*PI F*TGD*1.E-6 022=011

ANALIZ (C11,C12,C22,L11,L12,L22,R11,R12,R22,G11,G12,G22,F,L

CN,CK,R,LV,LD,CP.D1.D2,D3,D4,D5 D6,NS) RETURN

JHPAV (E.E2,E3,E4,H1.H2,H3,W,S,WP.LM,N,T,ZZ,ZMC11

C12,C22,L11,L12,L22,L,I,CX12,CY12) REAL LM,MO

DOUBLE PRECISION Q,Ql 1.Q21,Q31,Q12,Q22,Q32,Q42,Q13,Q23,Q33 . 9il Я'9 .C12,C22,L11,L12,L22,L,LXI I LY l.CXl l.CYl 1 CXlb CY12 CX22,CY22,Q,C1,C2,C3,C4,C5,C6C7,C8

EO= 8.85419

M0= 1.25663

Q11 = C1 (WP,H1-T)

Q21 = C2 (WP,LM,H2)

Q31 = C1(LM,H3-T)

g=(Qll+Q2l)*(Q31+Q21)-Q21*Q21

L12=MO*Q2l/Q

L22=L12+MO*Qll/Q

LXll = L12+MO*Q31/Q

CXl2=EO*E2*Q21

CX22=CX12+EO*E3*Q31

CX1 = CXI2+E0*E1*Q1I

L=N* (W+S)-S

G=(LM-WP)/L

N1=N-1

N2= N-2

Q12=C3 (W.S.Hl-T)

Q22=C3 (W,S,H2)

Q32=C5 (W,S,Hl-T)

Q42=C7 (W,S,H2)

CYl2=EO*E2*G* (N2*Q22+2*Q42)

CYl l = CY12 + EO*El*G* (N2*Q12+2*Q32)

CY22=CY12 + EO*G*Nl*S/ ((Hl-T)/El + T/E4+H2/E2)

Q13=C4 (W.S,H1-T)+4*T/S

Q23=C4 (W,S,H2)

Q33=C6 (W,S,Hl-T)+2*T/S

Q43=C8 (W,S,H2)

LYl 1 = MOG* (N2/(Q13+Q23) +2/(Q33+Q43))

C11 = CX11+CY11

C12=CX12+CY12

C22=CX22+CY22

L11 = LX11 + LY11

,CALL MODA (Cll,Cl2,C22,Lll,L12,L22,ZZ,ZM,ALFAE,ALFAO,I,Rll,R12.