RETURN END

FUNCTION CI (A,B)

DOUBLE PRECISION X,CC,CLA,B,CIC

X=DEXP (-I.570796*A/B)

CC=CIC (X)

CI = CC

RETURN

FUNCTION C2 (A,B,H)

DOUBLE PRECISION PIH,AA,BB,X,C2,C1C PIH=3.14I593/H AA=DEXP (P1H*A) BB=DEXP (PIH*B)

X=DSQRT (2.*(AA + BB)/((I+AA)*(1 + BB))) .

C2=CIC (X)

RETURN

FUNCTION C3 (W,S,H)

DOUBLE PRECISION Q,P,Z,T1,T2,T3,T4,X,C3,H,W,S,CCI Q=DEXP (-6.283185*H/(W+S)) P=DSQRT (DSQRT (Q)) Z=1.570796*W/(W + S)

Tl = 2*(P*DSIN (Z)-P**9*DSIN (3*Z)+P**25*DSIN (5*Z)>

T2=2*(P+P* 9+P* 25)

T3= 1 + 2*(Q + Q 4+Q 9)

T4= I-2 (Q*DCOS (2*Z)-Q 4*DC0S (4*Z)+Q**9*DCOS (6*Z))

X = Т1*Т2/ fT3*T4) C3=2 CC1 (X) RETURN END

FUNCTION C4 (W,S,H)

DOUBLE PRECISION Q,P,Z,T1,T2,T3,T4,X,C4,H,W,S CC 1 Q=DEXP (-6.283185*H/(W + S)) P= DSQRT (DSQRT (Q))

Z=1.570796*W/(W+S)

) -P 9* T2=2* (P+P 9+P 25)

T1 = 2*(P*DSIN (Z)-

DSIN (3*Z) + P 25*DSIN (5*Z))

T3= 1+2*(Q+Q*M + Q 9) T4= 1-2*(Q*DCOS (2*Z)-Q 4*DCOS (4*Z)+Q 9*DCOS (6*Z)) X=T1*T3/(T2*T4) C4=2.*CC1 (X) RETURN END

FUNCTION C5 (W,S,H)

DOUBLE PRECISION PIH,A1,B1,X,Y,C5,C1C

PIH=3.141593/H

Al= (PIH*(W+S/2.))

Bl=PIH*S/2.

X=DCOSH (AI)

Y=DCOSH (BI)

X = DSQRT ((l+Y) / (1+X))

C5=C1C (X)

RETURN

FUNCTION C6 (W,S,H)

DOUBLE PRECISION PfH,ALBl,X,Y,C6,ClC PIH=3.I41593/H Al=PIH*(W + S/2.) Bl = PIH*S/2.

X=DCOSH (AI) Y=DCOSH (BI) X= DSQRT ((Y-1 C6=C1C (X) RETURN END

)/(X-L))

FUNCTION C7 (W,S,H)

DOUBLE PRECISION PIH,A1,BLX,Y,A,B,C,C7,C1C

PIH=3.141593/H

AI = PIH*(W + S/2.)

Bl = PIH*S/2.

X = DCOSH (AI)

Y=DCOSH (BI)

A = X-1.

B=X-Y

C=Y+1.

X= DSQRT (C*(A+B + C)/((C+A)*(C+B)))

C7=C1C (X)

RETURN

FUNCTION C8 (W,S,H)

DOUBLE PRECISION PIH,A1,BI,X,Y,A,B,C,C8,C1C

PIH=3.141593/H

A1 = PIH* (W4-S/2.)

Bl = PIH*S/2.

X=DCOSH (AI)

Y=DCOSH (BI)

A=X+1.

B = X-Y C=Y-1.

X = DSQRT(C*(A4B+C) / ((C+A)* (C + B)))

C8=C1C (X)

RETURN

FUNCTION CIC (X)

DOUBLE PRECISION A,B,C,D,X,C1,C1C A=l.

B=DSQRT (1-X*X)

D=(A+B)/2.

B=DSQRT (A*B) A=D

IF ((A-B).GT.1.E-11*A) GOTO 1

C=3.141593/(A+B)

A=l.

2 D=(A+B)/2 B=DSQRT (A*B) A=D

IF ((A-B).GT.1.E-11*A) GO TO 2

C1 = 3.141593/(A+B)

C1C=C1/C

RETURN

FUNCTION CCl (X)

DOUBLE PRECISION A,B,C,D,X,C1,CCI

A=l.

B=DSQRT (1-X*X)

1 D=(A+B)/2. B=DSQRT(A*B) A=D

IF ((A-B).GT.1.E-llA) GO TO 1

C=3.141593/(A+B)

A=l.

2 D=(A+B)/2 B=DSQRT (A*B) A=D

IF ((A-B).GT.l.E-11*A) GO TO 2

C1 = 3.141593/(A + B)

CC1 = C/C1

RETURN

SUBROUTINE MODA (Cl,C2,C3,Ll,L2,L3,ZZ,ZM,ALFAE,ALFAO,I,Rll,

* R12,R22)

DOUBLE PRECISION Q1,Q2,Q3,Q4,L1,L2,L3,C1,C2 C3,DET E1,E2 A В С

* D,K,Z1I,Z12,Z22,P,Q,ZZ,ZM Q1=(L1*C1-L2*C2)/11.12 Q2=(L2*C3-L1*C2)/11.12 Q3=(L2*C1-L3*C2)/11.12 Q4= (L3*C3- L2*C2) /11.12 DET=DSQRT ((QI-Q4)* 2+4*Q2*Q3) El=(Ql+Q4 + DET)/2. E2=E1-DET

IF (E1.GT.0.AND.E2.GT.0) GO TO 10 E1 = DABS (El) E2=DABS (E2) 10 A=(E1-Q1)/Q2 B=(E2-Q1)/Q2

C=29979E-4*(C1-A*C2)/DSQRT (El)

D= 2.9979E-4*(C1 -B*C2)/DSQRT (E2)

IF (I.EQ.I) PRINT 1,E1,A,C,E2,B,D

FORMAT ( 1Х,МОДАЛЬНЫЕ ПАРАМЕТРЫV2X,E1 =,E15.8,2X,A E15.8,2X,C=,E15.8/2X,E2=,E15.8,2X,B=,E15.8,2X,D=,E15.8) K= A*B/ (C*D*(A-B)) Z11=(C/B-D/A)*K Z12= (C-D)*K Z22=(C*B-D A)*K P=Z11+Z22 Q=Z11*Z22-Z12*Z12 ZZ=P/2.+DSQRT (P*P/4.-Q) ZM=P/2.-DSQRT (P*P/4.-Q) AE=ALFAE/8.686 AO=ALFAO/8.686 K1 = 2./(C*(D/B-D/A)) Rll= (C*AO/B-D*AE/A)*K1 R12=(C*AO-D*AE)*Kl R22= (C*B*AO-D*A*AE)*K1 IF (I.EQ.2) PRINT 2,R11,R12,R22

FORMAT (1Х,МАТРИЦА СОПРОТИВЛЕНИЙ71Х, Rl 1 = ,F10.7/ * 1X,R12=,F10.7/1X,R22=,F10.7) RETURN END

С С С

PROGRAMM KANON АВТОР - ФЕДОРОВ В. Н.

COMPLEX В.и (16) COMMON/EINS/B (4,5).ILL,M,N REAL L,LI1,L12,L22 DIMENSION A(16),FI(16)

COMPLEX ZY1,ZY2,ZY3,ZY4,Z11,Z12,Z22,Y1 1 Y12 Y22

1 KE,KO,Y1E,Y10,Y2E,Y20VE,VO,BD ED

2 Zl,Z2,Z3,Z4,AE,AO,DE,DO Rl = 50.0

R2=50.0 R3=50.0 R4=50.0 C1 = 0.0 C2=0.0 C3=0.0 C4=0.0 1000 TYPE 10

10 FORMAT (IX,ПРОДОЛЖИМ? - Ь, НЕТ - О )

ACCEPT 20,1 20 FORMAT (II )

IF (I.LTl) GO TO 99

TYPE 50

50 FORMAT (IX, !.....El......В..........E2......В....

1 !.....ЕЗ.....,В.....!.....Е4......В... П

ACCEPT 40,Е1,Е2,ЕЗ,Е4

40 FORMAT (4F12.5) TYPE 30

30 FORMAT (IX,!.....Rl.....,0М.....!.....R2.....ОМ.

2 !.....R3.....,0М.....!.....R4.....,0М.....i)

ACCEPT 40,R1,R2,R3,R4

С L -ДЛИНА СЕКЦИИ В МЕТРАХ:

РАЗМЕРНОСТЬ - ПФ/М:

L=0.15

С11 = 429.17 С 12=389.01 С22= 468.65

РАЗМЕРНОСТЬ - МКГН/М:

L11 = 0.2584 L12=0.1162 L22=0.1785

RIJ,GIJ - ПОГОННЫЕ ПОТЕРИ, ОМ/М, 1/ОМ*М:

R11 = 0. R12=0.

R22=0. Gl 1 = 0.

1Э.Д.С. El 1Э.Д.С. E2 1Э.Д.С. ЕЗ 1Э.Д.С. E4

(4,5)

G12=0. G22=0.

PRINT 400,C11,C12,C22 PRINT 500.LI 1,L12,L22 400 FORMAT (1X,C1 1 = ,G12.5,1X,C12=,G12.5,1X,

1 C22=,G12.5,3X,-ПФ/М') 500 FORMAT (1X,L11 = ,G12.5,1X,L12=,G12.5,1X,

2 L22=,G12.5,3X.MKrH/M7) Z1 = CMPLX (R1,C1) Z2=CMPLX (R2,C2) Z3=CMPLX (R3,C3) Z4=CMPLX (R4,C4) B(1,5) = CMPLX (El.CI) B(2,5) = CMPLX (E2,C2) B(3,5) = CMPLX (ЕЗ.СЗ) B(4,5) = CMPLX (E4,C4) TYPE 400,C11,C12,C22 TYPE 500,L11,L12,L22 TYPE 900,B(1,5),B (2,5),В (3,5),В TYPE 800,Z1,Z2,Z3,Z4 PRINT 800,Z1,Z2,Z3,Z4 PRINT 900,B (1,5),В (2,5),В (3,5),В (4,5)

900 FORMAT (El = ,2G12.5,E2 =,2012.5,/, 3 E3 =,2012.5,E4 = ,2G 12.5/) ED=(1.,0.) PI = 3 1415926

F = 100 ШАЧАЛЬНАЯ ЧАСТОТА FO МГЦ

N1=1 ! КОЛИЧЕСТВО ЧАСТОТ. ШАГ=РО

WLi=WL*J

D0 1J=1,N1 ! ЦИКЛ В (1,5) = CMPLX (El,CI) В (2,5)=CMPLX (E2,C2) В (3,5)=CMPLX (E3,C3) В (4,5) = CMPLX (E4,C4)

WL = WUJ

F=WL/6.2831852

PRINT 201

PRINT 200,F

W0=WL*1.E-6

Z11=CMPLX (R11,WL*L11)

Z12=CMPLX (R12,WL*L12)

Z22=CMPLX (R22,WL*L22)

Yll = CMPLX (GlLWOCll)

Y12= CMPLX (G12,-W0*C12)

Y22 = CMPLX (G22,W0*C22)

ZY1 = Z11*Y11+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)

VE=CSQRT ((ZYl+ZY4 + BD)/2.)! ПОСТ. РАСПР0СТР. ЧЕТ. ВОЛН AE,DE

VO=CSQRT ((ZYl +ZY4-BD) /2.) ПОСТ. РАСПРОСТР. НЕЧЕТ. ВОЛН AO, DO

ПО ЧАСТОТАМ ! ЭДС. El ! Э.Д.С. Е2 ! Э.Д.С. ЕЗ ! Э.Д.С Е4

10 Заказ 0376

ПРОГРАММА РАСЧЕТА ВОЛНОВЫХ ПРОЦЕССОВ В УСТРОЙСТВАХ

НА СПЛ CANON