20 PRINT"ANTIMET ELLIP FLTR,CT12/78,1008" 1030 DIM B(16),C(16),D(16),E(30),F(16),R(15),S(15),DB(16),TB(16) 1040 DN=LOG(10)/10:PI=3.1415926 1050 PRINT 1060 PRINT"REJECTION,RIPPLE(DB),1/2-DEG(2-15),TYPE(A,B,OR C):" 1070 INPUT AS,AP,M,T$ 1080 IF AS<=AP THEN STOP 1090 N=2*M 2010 ES=EXP(DN*AS)-1 2020 EP=EXP(DN*AP)-1 2030 V=SQR(ES/EP)+SQR(ES/EP-1) 2040 U=PI*PI/(2*LOG(V+V)) 2050 V=V/(SQR(ES)+SQR(ES+1)) 2060 W=U*LOG(V+SQR(V*V+1))/PI 2070 W=SIN(W)/COS(W):AO=W:W=W*W 2080 Y=EXP(-U):Z=Y:K=M-1 2090 FOR J=1TON 2100 E(J)=1:NEXT J 2110 FORJ=1TO1024 2120 IF K<>M GOTO2150 2130 X=((1-Z)/(1+Z))^2 2140 AO=AO*(W+X)/(1+W*X) 2150 E(K)=E(K)*(1-Z)/(1+Z) 2160 Z=Z*Y:IF Z<.25E-18 GOTO2180 2165 K=K-1 2170 IFK=0 THEN K=N 2175 NEXTJ 2180 E(M)=0:E(N)=E(N)*E(N) 2190 PRINT"U=";U;"AO=";AO;"EP=";E(N) 2200 FOR J=1TOM-1 2210 E(J)=-E(J)*E(N-J) 2220 PRINT"E=";-E(J) 2230 E(N-J)=-E(J):NEXT J 2250 X=SQR(AO*AO+1/(AO*AO)+E(N)*E(N)+1/(E(N)*E(N))) 2260 FOR J=1TOM-1STEP2:K=(J+1)/2 2270 Y=AO*E(J):Y=Y+1/Y 2280 Z=E(N)*E(J) 2290 R(K)=E(M-J)*(1/Z-Z)/Y:S(K)=-X/Y 2300 PRINT"RE";R(K);"SE=";S(K) 2310 R(M-K+1)=R(K) 2320 S(M-K+1)=-S(K):NEXT J 2330 IF K+K=M GOTO3010 2340 R(K+1)=-AO:S(K+1)=0 2350 PRINT"RE=";-AO 3010 IT=2:IF T$="A" THEN IT=1 3020 E8=-E(1):IF T$="A" THEN E8=E(N) 3030 EO=E(N):IF T$="C" THENEO=-E(1) 3040 FP=SQR((E(N)+EO)/(1+E(N)*E8)) 3050 FS=SQR((1+E(N)*EO)/(E(N)+E8)) 3060 D(1)=0 3065 FOR J=ITTOM 3070 D(J)=(E(2*J-1)+E8)/(1+E(2*J-1)*EO) 3080 F(J)=SQR(1/D(J)):NEXT J 3100 SR=0:TQ=0:T0=0:B(1)=0:I=1 3110 FOR J=1TOM 3120 W=(AO^2+E(2*J-1)^2)/(1+(AO*E(2*J-1))^2) 3130 X=(1+EO*E8)*S(J)+EO+E8*W 3140 Y=EO^2+2*EO*S(J)+W 3150 Z=1+2*E8*S(J)+E8^2*W 3160 U=SQR(Y/Z):V=X/Z 3170 R(J)=SQR((U-V)/2):S(J)=SQR((U+V)/2) 3180 PRINT"RF=";-R(J)/FP;"SF=";S(J)/FP 3200 SR=SR+R(J)/U 3210 I=-I:W=I*R(J)/S(J) 3220 TQ=(TQ+W)/(1-TQ*W) 3230 IF T$<>"A" GOTO3270 3240 U=(F(2)-S(J))/R(J):V=(F(2)+S(J))/R(J) 3250 W=I*(V-U)/(1+U*V) 3260 T0=(T0+W)/(1-T0*W) 3270 B(1)=B(1)+R(J):NEXT J 4010 IF T$="A" THEN T0=T0/(1+SQR(1+T0*T0)) 4020 FOR K=ITTOM 4030 DB(K)=0:TB(K)=T0:I=1 4040 FOR J=1TOM 4050 DB(K)=DB(K)+1/(R(J)+(F(K)-S(J))^2/R(J))+1/(R(J)+(F(K)+S(J))^2/R(J)) 4070 I=-I:W=(F(K)-I*S(J))/R(J) 4080 TB(K)=(TB(K)+W)/(1-TB(K)*W):NEXT J:NEXT K 5010 D(M+1)=D(M):F(M+1)=F(M):DB(M+1)=DB(M):TB(M+1)=TB(M):C(1)=0 5020 FOR J=1TOM+1-IT STEP2 5030 TB(M+1-J)=-1/TB(M+1-J):NEXT J 5040 FOR J=ITTOM+1 5050 B(J)=(1+TB(J)^2)*DB(J)/(4*D(J))-TB(J)*F(J)/2 5060 C(J)=TB(J)/F(J):NEXT J 6010 FOR L=1TO2 6020 FOR K=L+2TOM+1 STEP2 6030 FOR J=LTOK-2STEP2 6040 U=C(J)-C(K) 6050 V=1/(U/(B(J)*(D(K)-D(J)))-1) 6060 B(K)=(B(K)-B(J))*V*V-B(J)*(V+V+1) 6070 C(K)=U*V:NEXTJ:NEXTK:NEXTL 7010 W=1:IFT$<>"C" THEN W=((1-TQ*T0)/(TQ+T0))^2 7020 FORJ=1TOM+1 STEP2:B(J)=B(J)*W 7030 C(J)=C(J)*W:NEXT J 7040 PRINT"LD RESIS=";W;" (";1/W;" )" 7050 PRINT" L(C) C(L) PEAK" 7070 IF T$<>"A" THENPRINT" 1 ";FP/B(1) 7085 V=0 7090 FOR J=ITTOM:V=V+C(J) 7100 PRINTJ;FP*C(J);FP/B(J);F(J)/FP:NEXTJ 7110 PRINTM+1;FP*C(M+1);"STPBD EDGE=";FS/FP 7130 PRINT"TESTS";B(M)/B(M+1)-1;(W+1)*SR-V-C(M+1) 7140 STOP 7150 END