0043 CALL TEST(IFR,IV,ICIN.IRIN.TEXLIN)

0044 IF(.NOT.STOP) GOTO 150

0046 RETURN

0047 END

0001 SUBROUTINE GPl

С ««*«««***»•♦♦•«=*♦»»«»«»»«*««««»««««»»«*»«»»»»»»*»»«»*»*»»♦»******

0002 COMMON/STO/STOP

0003 COMMON/FIRST/IFIRST

0004 COMMON/DATA 1 /XFR.XQP.XKH P.XKBP.XKLP.XROPT.XR 1. XR2.XR3.XR4,XR5.XC6

1,XR7.XC8,XC,XR

0005 DIMENSION TEXl (16),IN(5)

0006 LOGICAL*! STOP

0007 DATA TEX1/FREQ/QP,KHP,KBP,KLPVR0PTVR!VR2,

1 ,R4,R5,C6/R7/C8/C.R7IN/1,2,0,15,16/

0008 DATA 1FR.IV.IC1N.ICOUT,IRIN,IROUT,IPAV1,4,1,2,0,7,4/

0009 DATA ICH.IRINl/0.1/

0010 CALL 1NPUT(IFR,IV,1CIN,IRIN,ICH,TEX!,IN)

0011 150 ICH=!6

0012 IFIRST=0

0013 XHI=6.2831835*XFR»XC

0014 XR0PT=!./XHI

0015 TYPE !0,XROpT

0016 10 FORMAT(X.ROPT,T!5,=,!PE10.3)

0017 CALL INPUT(IFR,IV,ICIN.IRIN!,ICH.TEX1,IN)

0018 XR1=XR

0019 XR3=XR

0020 XR5=XR

0021 XR7=XR

0022 XC6=XC

0023 XC8=XC

0024 XR4=XR*(XHI*XR)**2

0025 XR2=XR*(XQP*(1.+XR4/XR) (SQRT(XR4/XR)-I.?

0026 XKHP=(!.H-XR4/XR)/(XR/XR2+!.)

0027 XKBP = XR2/XR

0028 XKLP=XKHP*XR/XR4

0029 CALL 0UTPUT(IPA.IC0UT,IR0UT.TEX1)

0030 CALL TEST(IFR,IV.IGIN,IRIN,TEX!,IN) 003! IF(.NOT.STOP) GOTO 150

0033 RETURN

0034 END

ООО! SUBROUTINE GP2

С ♦•»»»«»»«»»»»««•»••*«»»»»»»»»»•»»»•»»»»»•*»»»*»*»»»»»*•*•»**•»»**

0002 COMMON/STO/STOP

0003 COMMON/FIRST/IFIRST

0004 COMMON/DATAl/XFR,XQP,XKBP,XKLP,XROPT,XRl,XR2, XC3,XR4,XR5,XC6

!,XR7R8,XC.XR

0005 DIMENSION TEXl(!6),IN(7)

0006 LOGICAL*! STOP

0007 DATA TEX!/FREQ,QP,KBP,KLP,ROPT,R!VR2,C3, rR4,R5,Ce.R7/R8/C,R/ 7IN/1 .3,4,0.

0008 DATA IFR.IV.ICIN,1COUT,1RIN,IROUT,IPA/!,4,2,0,7,4/

0009 DATA ICH.IRINl/0.1/

0010 CALL INPUT(IFR.IV,ICIN,IRIN,ICH,TEX1.IN)

0011 150 XKLPIN.= XKLP

0012 1СН=15

0013 IF1RST=0

0014 XROPT=l./(6.2831853*XFR*XC)

0015 TYPE lO.XROPT

0016 10 FORMAT(X,ROPT,T15,=MPE10.3)

0017 CALL 1NPUT(1FR,1V,ICIN,1RIN,1.ICH.TEX1,IN}

0018 XR2=XR

0019 XR7=XR

0020 XR8=XR

0021 XC3=XC

0022 XC6=XC

0023 XR4=XQP*XR0PT

0024 XR5=XROP**2/XR

0025 IF(XKLPIN.EQ.O) GOTO 200

0027 XR1=XR2/XKLP

0028 XKBP=XR4/XR1

0029 GOTO 220

0030 200 XR1 = XR4/XKBP

0031 XKLP=XR2/XR1

0032 220 CALL OUTPUT (1PA,1C0UT,1R0UT,TEX1)

0033 XKLP = XKLP1N

0034 CALL TEST(1FR,IV,1CIN,IR1N.TEX1,1N)

0035 IF(.NOT.STOP) GOTO 150

0037 RETURN

0038 END

0001 SUBROUTINE HLPNMQ

С **»»»»»»»***«*****»***»*»»»****»»♦«»*»*»**»*»♦**»««««««»♦♦♦♦*♦♦♦»

0002 COMMON/STO/STOP

0003 COMMON/FIRST/IFIRST

0004 C0MM0N/DATA1/XFRP,XFRZ.XQP,XK,XGSP,XP,XPMIN,XR1. XR2,XC3,XC4,

1XR5,XR6,XR7,XR8,XR9

0005 DIMENSION TEXl (16),IN (9),XX (3)

0006 LOGICAL* 1 STOP.OPT.RICH

0007 DATA TEX1/FRP,FRZ,QP,K,GSP/P,PMIN.R1,R2,

1 C4 VR5,R6,R7,R8/R97I N/1,2,3,4,6,10,11,16,0/

0008 DATA IFR,IV,IC,IRIN,IROUT,IPA,ICH/2,3,2,1,7,7,0/

0009 XP = 0.

0010 CALL INPUT(IFR,IV,IC,IRIN,1CH,TEX1,IN)

0011 80 XK1 = XK

0012 XPIN=XP

0013 IF((XFRP-XFRZ).NE.O.) GOTO 88

0015 TYPE 5

0016 5 FORMAT(/X,NO BANDREJECT)

0017 XP = XPIN

0018 CALL TEST(IFR,1V,IC,1RIN,TEX1,IN1

0019 GOTO 80

0020 88 XP=0.5

0021 XX (2) =0.

0022 ILAUF=0

0023 OPT==.FALSE.

0024 RICH=.FALSE.

0025 89 ILAUF=ILAUF+1

0026 90 IF(.NOT.OPT) GOTO 95

0028 XP = XPMIN

0029 IF(XPIN.EQ.O) GOTO 100 0031 XP=XPIN

0032 GOTO 100

0033 95 XPMIN=XP

0034 IF (RICH) XP = XP+0.01 0036 IF(.NOT.RICH) XP=XP-0.01

0038 100 XWP=XFRP*6.2831853

0039 XWZ=XFRZ*6.2831853

0040 X=0.

0041 IF((XFRP-XFRZ).GT.O.) X=l

0043 XH1 = 1 .+4.*XQP*XQP*XP* (1 .+XC4/XC3)

0044 XH2=XC3*XWP*(SQRT(XH1)-l.)/(2.*XP*XQP)

0045 XKO= (l.+XP)/(l.+ (l.+XC4/XC3)*XWZ*XWZ*(XC3/XH2) •*2)

0046 IF(XKIN.NE.O.) GOTO 150

0048 IF((1-XKIN/XKO).GE.O) GOTO 150

0050 XK=XKO

0051 XR1 = 1./XH2

0052 XR2=1.E38

0053 GOTO 160

0054 150 XK=XKIN

0055 XR 1 = ХКО/ (XH2*XK)

0056 XR2=1./(XH2*(1.-XKIN/XKO))

0057 160 XR6=XH2*(1.+XP)*(1.-(X/XK))/(XC3*XC4*(XWZ**2-

-XWP**2))

0058 XR5= l./(XC3*XC4*XWP**2/XH2+XP/XR6)

0059 XH4=ABS (1.-ХН2/ (XWP**2*XC3*XC4) * (1 ./XR5+1/XR6) >

0060 XH5=ABS (1.-XQP/XWP* (XH2/XC3f (1 ./XC3-bl./XC4) •(1./XR5+1./XR6)))

0061 XGSP== (l.+XP)/2.*(XH5+XQP*XH4)

0062 XR7=XP*XR9/XK

0063 XR8=XP*XR9/(1.-XK)

0064 XX(ILAUF)=XGSP

0065 IF (XX (1) .LT.XX (2)) RICH = .TRUE. 0067 IF(ILAUF.LE.2) GOTO 89

0069 IF (OPT) GOTO 600

0071 IF(XPMIN.LE.0.02) OPT=.TRUE.

0073 IF (RICH) GOTO 550

0075 IF(XX(2) .LT.XX (3)) OPT=.TRUE.

0077 GOTO 560

0078 550 IF(XX(2).LT.XX(3)) OPT=.TRUE.

0080 560 XX(2)=XX(3)

0081 GOTO 90

0082 600 IF(XFRP-XFPZ)620,80,610

0083 610 TYPE 10

0084 10 FORMAT(X,HIGH-PASS)

0085 GOTO 630

0086 620 TYPE 11

0087 11 FORMAT (X/LOW-PASS)

0085 630 CALL 0UTPUT(IPA,IC,IR0UT,TEX1)

0089 XP=XPIN

0090 CALL TEST(IFR,IV,IC,IRIN,TEX1,IN)

0091 IX(.NOT.STOP) GOTO 80

0093 1000 RETURN

0094 END

0001 SUBROUTINE BRHLPN

С ******•*******•****•**•*••**•*****************•***********•****•*

0002 COMMON/STO/STOP

0003 COMMON/FIRST/IFIRST