Главная -> Книги (0) (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) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73) (74) (75) (76) (77) (78) (79) (80) (81) (82) (83) (84) ( 85 ) (86) (87) (88) (89) (90) (91) (92) (93) (94) (95) (96) (97) (98) (99) (100) (101) (102) (103) (104) (105) (85) 0007 DATA IFR.IV,IC,IRm,IROUT,IPA,ICH/l,2,3,l,4,6,0/ 0008 xpm=o. 0009 XP = 0. 0010 CALL INPUT/lFR,IV,IC,lRm,lCH,TEXl,m/ ООП 150 XC1=XC11+XC12 0012 XCH1=XC4/XCI 0013 XQHI = XQP*XQP 0014 XPMm = I/(XCHI*36.*XQHI)*(SQRT(l.+12.*XQHI*(XCH[+ + 1.))+1.)**2 0015 IF(XP.NE.O) GOTO 160 0017 ХР=ХРМт 0018 160 XR2=1/(SQRT(XP*XC1*XC4)*(6.2831853*XFR)) 0019 XR3 = XP*XR2 0020 XR5=XR6*(XCHI*(1.+XP)-SQRT(XCHI*XP)/XQR) 0021 XK=XC11/XC1*(1.+XR5/XR6)*XQP*SRRT(1./(XP*XCHI)) 0022 XGSP= (l.+XR5/XR6)**2*XQP*SQRT(l./(XCHI*XP)) 0023 CALL OUTPUT (IPA.ICIROUT.TEXl) 0024 ХР=ХРШ 0025 CALL TEST(IFR,IV,IC,IRIN,TEX1,1N) 0026 XP1N=XP 0027 IF(.NOT.STOP) GOTO 150 0029 RETURN 0030 END 0001 SUBROUTINE HPMQ С ***************************************************************** 0002 CO.MMON/STO/STOP 0003 C0MM0N/DATA1/XFR,XQP,XK,XGSP,XP,XPMIN,XC11,XC12, XR2,XC3,XR4 1,XR5,XR6 0004 DIMENSION TEXl (16),IN(7) e005 LOGICAL*! STOP 0006 DATA TEXl/FREQVQP/K/GSP/P/PMIN/Cl 1/C12,R2, C3 rR4VR5/R6/ V/ 7IN/1,2,5,7,8,10,12/ 0007 DATA IFR,IV,IC,IRIN,IROUT,IPA,ICH/1,2,3,1,4,6,0/ 0008 XPIN=0. 0009 CALL INPUT(IFR,IV,IC,IRIN,ICH,TEX1,IN) 0010 150 XC1=XC11+XC12 ООП XCHI = XC3/XC1 0012 XQHI = XQP*XQP 0013 XPMIN = 1/(XCHI*4.*XQHI)*(SQRT(1.+12*XQHI*(XCHI+ + 1.))-1.)**2 0014 IF(XPIN.NE.O) GOTO 160 0016 XP = XPMIN 0017 160 XR2=1/(SQRT(XP*XC1*XC3)*(6.2831853*XFR)) 0018 XR4=XP*XR2 0019 XR6 = XR5*((1.+ 17XCHI)/XR-(SQRT(L/(XCHI*XP))/XQP)) 0020 XK=XC11*(1.+XR6/XR5)XC1 0021 XGSP = XQP* (1 .+XR6/XR5) **2*SQRT (XP*XCHI) 0022 CALL OUTPUT(IPA,IC,lROUr,TEXl) 0023 XP = XPIN 0024 CALL TEST(IFR,IV,IC,IRIN,TEX1,IN) 0025 XPIN = XP 0026 IF(.NOT.STOP) GOTO 150 0028 RETURN 0029 END О001 SUBROUTINE APMQ С ********«««*««««*«*«***«****»*****«*****«««*««*«*****«•*•*****» §002 COMMON/STO/STOP <to03 COMMON/DATA 1 /XFR,XQP,XK,XGSP,XP,XPMIN,XR 1 ,XG2, XC3,XR4,XR5,XRe,XR7 0004 DIMENSION TEXl (16),Ш(в) 0005 LOGICAL*! STOP 0006 DATA TEX!/FREQ,QP,K/GSP/P/PMIN/R!/C2,C3, R4, rR5/R6,R7, 7IN/1,2,5,8,9,!2/ 0007 DATA IFR,IV,IC,IRIN,IROUT,IPA,ICH/!,2,2,!,5,6,0/ 0008 XP=0. 0009 TYPE 15 0010 CALL INPUT(IFR,IV,IC,IRIN,ICH,TEX1,1N) 0011 150 XPIN=XP 0012 XCH1=XC2/XC3 0013 XQHI = XQP*XQP 0014 XHIL=((6.*(1.+XCHI)**2*XQHI)/(2.*XQHI*(1+XCH1)- XCHI))**2 0015 XHIL = (SQRT (1 .+XHIL/(3.* (1+XCHI) ..2))-1.) • (2.*XQHI *(1+XCHI)-XCHI) 0016 XPMIN=XHIL/(6.*(1+XCH1)**2*XQHI) 0017 IF(XPIN.NE.O) GOTO 160 0019 XP = XPMIN 0020 160 XRI = 1/(SQRT(XP*XC2*XC3)*(6.2831853*XFR)) 0021 XRP=XP*XR1 0022 XR7=XR6* (XP* (1 .+XCHI) +SQRT (XP*XCHI) /XQP) 0023 XHl L = 1 .-2* (SQRT (XP*XCHI)) / (XQP* (1 .+XR7/XR6)) 0024 XR5=XRP/XHIL 0025 XR4=XRP/(1.-XHIL) 0026 XGSP=XQP*( l.+XR7/XR6)**2*XHIL*SQRT(l./(XP*XCHI)J 0027 XK= 1. 0028 TYPE 15 0029 CALL 0UTPUT(IPA,IC,IR0UT,TEX1) 0030 XP = XPIN 0031 CALL TEST(1FR.IV,IC,IRIN,TEX1,1N) 0032 IF(.N0T.ST0P)» GOTO 150 0034 15 FORMAT (/X/ALL-PASS) 0035 RETURN 0036 END 0001 SUBROUTINE BRMQ 0002 COMMON/STO/STOP 0003 COMMON/DAT Al /XFR,XQP,XK,XGSP,XP,XPM1N,XR 1 ,XC2, XC3,XR4,XR5,X R6,XR7 0004 DIMENSION TEXl (16),IN(3) 0005 LOGICAL*! STOP 0006 DATA TEX!/FREQ,QP,K,GSP,P.PMIN.R!,C2,C3, R4, lR5,R6,R7. ., 7IN/!,2.5,8,9,!2/ 0007 DATA IFR,IV,IC.IRIN,IROUT,IPA,ICH/!,2,2,!,5,6.0/ 0008 XP = 0. 0009 TYPE 15 0010 15 FORMAT (/X.BIAND-REJECTION) 0011 CALL INPUT(IFR,IV,IC,IRIN,ICH,TEX1,IN) 0012 150 XPIN=XP 0013 XCHI=XC2/XC3 0014 XPMIN=1./(3.*(1.+XCHI)) 0015 IF(XPIN.NE.O) GOTO 160 0017 XP=XPMIN 0018 160 XR1 = 1/(SQRT(XP*XC2*XC3)*(6.2831853*XFR)) 0019 XRP=XP*XR1 0020 XR7=XP6*XP*(1-J-XCHI) 0021 XHIL= 1.- (SQRT (XP*XCHI)) /(XQP* (1 .+XR7/XR6)) 0022 XR5=XRP/XHIL 0023 XR4=XRP/(1.-XHIL) 0024 XGSP=XQP*(1.+XR7/XR6)**2*XHIL*SQRT(1./(XP*XCH1)) 0025 XK=1. 0026 TYPE 15 0027 CALL OUTPUT (IPA.ICIROUT.TEXl) 0028 XP=XPIN 0029 CALL TEST(IFR,IV,IC,IRIN,TEX1,IN) 0030 IF(.NOT.STOP) GOTO 150 0032 RETURN 0033 END 0001 SUBROUTINE LPHQ С «*«««*♦«•♦«««**«««««»«*««*«•*•♦»*»«**•♦♦««««»»«*»*«♦*»»»»»»»*»**« 0002 COMMON/STO/STOP 0003 COMMON/FIRST/IFIRST 0004 COMMON/DATA 1/XFR,XQP,XK,XR0PT,XR1,XC1,XP2,XR3. XC4,XR6,XR7,XC,XR 0005 DIMENSION TEXl (16),IN(5) 0006 LOGICAL*! STOP 0007 DATA TEX 1 /FREQ/QP.K.ROPTVR 1VC1 .R2/R3.C4,R6, lR7VC,R,,, 7IN/1,2,0,12,13/ 0008 DATA IFR,IV,ICIN,ICOUT,IRIN,IROUT,IPA/1,2,1,2,0,6,3/ 0009 DATA ICH,IRIN1/0,1/ 0010 CALL INPUT(IFR,IV,ICIN,IRIN.ICH,TEX1,IN) 0011 150 ICH=13 0012 IFIRST=0 0013 XROPT=L/(6.2831853*XFR*XC) 0014 TYPE lO.XROPT 0015 10 FORMAT (X,R0PT,T15/=,1PE 10.3) 0016 CALL INPUT(IFR,IV,ICIN,IRIN1,ICH,TEX1,IN) 0017 XK=2 0018 XR2=XR 0019 XR3=XR 0020 XR6=XR 0021 XC1 = XC 0022 XC4=XC 0023 XRI=XQP*XROPT 0024 XR7=XROPT**2/XR 0025 CALL 0UTPUT(IPA,IC0UT,1R0UT,TEX1) 0026 CALL TEST(IFR,IV,IC1N,IRIN.TEX1.IN) 0027 IF(.NOT.STOP) GOTO 150 0029 RETURN 0030 END 0001 SUBROUTINE BPHQ 0002 COMMON/STO/STOP 0003 COMMON/FIRST/IFIRST 0004 COMMON/DATAl/XFR,XQP,XK,XROPT,XRl,XR2,XC3,XR4, (0) (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) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73) (74) (75) (76) (77) (78) (79) (80) (81) (82) (83) (84) ( 85 ) (86) (87) (88) (89) (90) (91) (92) (93) (94) (95) (96) (97) (98) (99) (100) (101) (102) (103) (104) (105) |
|