 |
|
|
Главная Волноводные диэлектрические фильтры 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 [ 21 ] 22 23 ПРИЛОЖЕНИЕ 2 56 57 98 59 л0 1 62 63 4 5 б i67 48 .69 7в 71 72 7J 74 75 16 .77 ,78 i79 8в 1 82 в) 4 .85 86 A=Sll+t. B=(Sll-I.}*PE RB=T(A.-8) A = IIB E1*RA/E1 B=l-RB*E .RA/Eil .PO RA=T{AiBl. RBsT{A.-8) A=RB E2.RA/Ea e=(-BB E2tRA/E2)i*PC RA=T(A)a} RBsT(Ai-B) S12=RB EJ+RA/E5 RR(K,l) = CAes(SU) TT(K,I)=CABS(SU) AA(K)=A0 lF{Ae-ABK)4t4,8 8 Ав=РАЯ DD(T)=ALD ALOALD+OALD AL1:AL1-0.5*DaLD ALO = ALO-8.5*DAl.O lF(l,Ee.7) CO TO 9 IFJALD-ALDK)7.7.9 , . 9 PRINT. 5, (ODILl .U=bl) lOOfL) ,L=1,1) DO IB J=tiK 18 PRINT je.AAfJ) ilRRIJ.M) >I1= .1) 1 fTKJ.H) it1=bl) PRINT 5,ALI,ALO ALO=FALD ALO=FALD AL1=FAL1 ALO=ALO.dDB *Ll=ALl*DDB IFAL0-*V.0K)6,6, J END Программа расчета на ЭВМ ЕС-1022 ПРОГРАММА РАСЧЕТА ЭЛЕМЕНТОВ МАТРИЦЫ РАССЕЯНИЯ ДИЭЛЕКТРИЧЕСКОГО ЦИЛИНДРА В ПРЯМОУГОЛЬНОМ ВОЛНОВОДЕ Программа расчета на ЭВМ Мннск-32 элементов матрицы рассеяния диэлектрического цилиндра в прямоугольном волноводе состоит из двух модулей: первый - головная программа, с помощью которой осуществляются ввод исходной информапии и вывод результатов расчета; второй - подпрограмма АВСВЕ, используемая для формирования правых частей дифференциальных уравнений. Исходные данные записываются на трех перфокартах и вводятся операторами 100; ПО; 140. Для контроля правильности ввода исходная информация выводится на АЦПУ с помощью операторов печати 120; 130; 150. Последовательность вывода в этом случае полностью соответствует последовательности ввода. На перфокарту 140 записано условное число L (тип целый), которое в случае вертикальной поляризации принимает значение 1, а в случае горизонтальной поляризации - значение 2. Смысл обозначений, использованных для записи исходных данных, указан в табл. П2.1. и П2.2. Таблица П2.1 Исходные данные (перфокарта 100)
Результаты расчета представляют собой таблицу, которая печатается на АЦПУ с помощью оператора 760 по формату 42. Расположение выходной информации по столбцам дано в табл. П2.3. В последней строке печатается значение DA (оператор 820), соответствующее текущей величине диаметра цилиндра. Численное интегрирование системы дифференциальных уравнений осуществляется по методу Рунге-Кутта четвертого порядка. Обращение к соответствующей подпрограмме RK4 выполняется с помощью оператора 250. В случае сбойной ситуации управление передается на метку 8 и печатается текстовая информация СБОИ в ПР RK4 . Нормированные граничные условия (2.22) за- Таблица П2,2 Исходные данные (перфокарта ПО)
даются операторами 230; 240; 300 и 310, Вычисление коэффициентов 5и и S21 в соответствии с (2.26), (2.27) осуществляется операторами 350-710. После завершения всех циклов головная программа восстанавливает свое начальное состояние и автоматически вводит новые исходные данные. Таблица П2.3 Выходные данные
Подпрограмма ABCDE формирует функции, находящиеся в правых частях дифференциальных уравнений. Конкретные значения этих функций, вычисляемых при обращении к RK4, заносятся в ячейки, отведенные переменным 5(2) и 5(3). Изменение ориентации оси цилиндра в прямоугольном волноводе влияет иа вид функции, задаваемой операторами 120 и 150. В табл. П2.4 приведен образец печати иа АЦПУ контрольного примера расчета для исходных данных, соответствующих табл. П2.1 и П2.2. ОО О О -УО -* <Ьо о о 0 0к о* о О О О: о О око 0 0 с>о оооо qkqk О 0 оо ооо оо о о - кмо tv о л ч> о о> - м ю о о . 0 л 1* 1М INI ♦г оо оОооооо оо ootpniAeo-rae-ooO -.Г4в л м л1<ч *ч 1м тем w ОоООООООООО о . ОвО-* 0 -*fA n, .ol. *-вoe-.or - 00 0000 ооо с ею го к л* А см л г-оо. о- . ..........* .оооооо ооо ол о 1А1Г -ivo i*-o* -moo о aOiAfMiA *OtAOeOO>m о сое П1Л^а> Ю>А|А^А . с о ю * -* п о оо ooooooooqo о г ч|в'*>л^.Г - -0 о iv-orof-rfie-oooo*- ,1Л ь^о v\ <л -а Л'Л >* -л .> -4 <г о ооооооооооо о о -в л в . efw - о 00 00 е о o-in fO о> <> <> о с о ооооооооооо о <ОЧ> Ю < пЛ! -00/М eO w м о .,лйГ1П ЛО вк *Ч)вОсОвООО о о>во л'в fb.fb..r*- Ф> о ю о> > < * о оо о о о о о о о о о о о оо oOOOCjeuo.rtm&.Orar-oO во ,оо1лосм>0вйвовово оо ч. (Ч1>#>А<0<0<0 *0 <0 оо ооооооооооооооо о о >л 4.<ooiov\mia ЛГ4Ю **>л оо г^о-го<м'* *л'л>Л1Л л >л^ о о 1ЛО вo o <cc( <-< * гщ о оо оооо оос оо оо о о о 1ЛО г . -rfM -о оо овО О о <-Kl О'* о w >лл* о О лгм о - мюююглююю о о .оО С0О>О* ОСЬ^О'ОО о -о> оо Ооооооооооо о о с < . 00 i*-. * * в о ЮО # .000 -*- -*- о о оо в ООО ООО оо оС-~ оо оооооооооооооо So SS оооооооооооо оо оооо ООО ооооо
580 И|=(0.,1.> (Г-1.)/ 590 Z=<Z*l.) 600 Zt=o,5 (Z-W) EXP(-0.S zS> 610 Z2=0.5 lZ+W)iEXP(D.5 Z5> 620 Z=Z1 630 W=Z2 640 Zl=YI*Z . 650 Z2 = (*<-Z) B2 660 Z = Zl lll Z2.ft(2) 470 w= ( z 1 ы-гг > I 3 ))/62 680 Zl=0.5 rZ-W) 690 Z2=0>5 lZ W) 700 Z=Zl EXPI-O.5 Z6).Z2 EXPl0,5 2 ) 710 BI 6)=CAbSZ) 720 T=T*ST 730 IT=TT ST 740 iFlTT.GT.rrEND) GO TO 31 750 50 CONTINUE 760 5I WH IТЕ1ПЧ1 ,42)AA,(B( I), 1= 1 ,6) .(b( I ), I=7 . 12 ) iBPoHC,AM 770 T=T1 7 8 0 T T = T T 1 790 AA=AA*SA 800 AN=0. 810 1F(AA-AEND)9.11.H 820 1 1 WRiTECnMl, DDA 830 AA=4A1 840 DA=DA-SD 850 lr(DA-DAEND)34,9,9 860 e WR I T E < ПЧ1 , 2 ) 870 Stop 881 End 010 SUBROUTINE ABCDEcAiS) 020 Dimension А131.si3) 030 Common e . aa ,dA,пи>an,pa>l 040 Г1(ТП=Т1 050 W(V)=SttRT(1.0000 0.1-(2,.V-l.) 2) 060 r2(T2)=R T2 070 X=All) 08 0 SI 2 )=Fir Al 3) ) 090 RR:nM 0A loo GO TOr 1.2) ,L no 1 CONTINUE 120 R=RR 2-l(Z, RR AA) 2).(l,-(l.-E) (0A W(X)+SiN(RR,w(X))/nn)) 130 Go TO 3 140 2 RR=RR>BA 150 R=RR 2-((2. RR AA) 2) (l.t(E-l.l.D4.BA W(X)) 160 3 CONTINUE 170 S(3)=F2(A( 2 ) ) 180 AN=0,001*AN 190 RETURN 200 End Программа расчета на ЭВМ Минск-32 ПРИЛОЖЕНИЕ 3 ПРОГРАММА АНАЛИЗА МНОГОЗВЕННОГО ЗАПРЕДЕЛЬНОГО ВОЛНОВОДНО-ДИЭЛЕКТРИЧЕСКОГО ФИЛЬТРА Настоящая программа предназначена для расчета на ЭВМ ЕС-1022 частотных характеристик коэффициента отражения многозвенных волноводно-диэлектрических фильтров с запредельными связями. Предполагаемся, что в качестве базового звена используется запредельный резонатор с плоским диэлектрическим слоем. В программе реализован алгоритм вычисления S)il путем последовательного преобразования комплексных сопротивлений с учетом дополнительных плоскостей отсчета. |