求大神帮忙看一下程序的问题出在哪里，如能帮忙修改万分感谢 - 数学建模社区-数学中国

CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
* W% t! ^' x# `2 T! V* e
Io (gm1*gm2/s/c1/s/c2)*Ig};
) L6 K. o: ^7 e4 @0 b
Current={Ii,Io};
' V" {+ e& ~, T6 I* S2 a
CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];, S. ], }* c3 r; d- {$ ]
H=Io/Ii/.CircuitRespones//Together//Simplify;
2 s5 Y* o1 K, v
Print["H(S)=",H]- x6 y" P2 B+ ?/ Z
PoleQpole[H,s]:=Module[{den,fp,Qp},
9 B! `* A3 h' k9 M
den=Denominator[H];2 D4 U/ E' h& v
fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
9 b. c& g/ S/ P% c$ ?
Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));9 s" `' L8 V* x
Simplify[{fp,Qp}]];+ K, T8 _* f1 K' |6 q3 V
ZeroQzero[H,s]:=Module[{fz,num,Qz0},+ Q: J, g p1 j% m Y& Q
num=Numerator[H];6 l6 J' m- `: J' K: V9 @
Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]); G+ N" {5 N0 b- v. M
fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);1 B6 ?$ Z0 R5 |$ ?. y# ~
Simplify[{fz,Qz0 fz}]];
7 l9 w) M- m$ h9 G8 N
{fp,Qp}=Simplify[PoleQpole[H,s]];
5 f5 r6 S [7 \/ ]
Klp=H/.s 0;
% \' E" x$ M# w
PrintLabeledList[{fp,Qp},{"fp","Qp"}];3 C z2 L5 k- r" B) j
DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},8 N1 k& S4 n- Z& H
gm1=C1*Wp/Qp;
/ ?% j8 I8 y4 c# {" _$ i- E
gm2=C2*Wp*Qp;
# E2 T* ^# w9 K* X) @3 a1 I
{C1,C2,gm1,gm2}];; L' Q3 ?. g6 [, T \9 x. B
{C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
# j7 y# a6 B5 d: G( T4 h
PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];- U* I; W) q% ~# g+ \ O
Simplify[H]
4 j+ F- D6 m4 h' N7 ^+ l, y. Y
h1=Together[H/.{Q 2.,W N[2*pi*10^6],
! G/ C# V9 y, B6 b, \) f: f* s. j
c1 10.*10^(-12),c2 10*10^(-12)}//N];1 c& H& Z3 m& t) o s7 t
Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
" k! J4 A; ^& R+ W5 H3 i
c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]! i. \0 f d# V! t( ^
Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],- ~1 g! P7 r4 _0 x9 u! R5 {! |
c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]1 I2 y- x6 m: V7 f
h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])2 i( c$ j% | T! P- g. I6 t
PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],0 Z' D6 D6 ~3 H+ ^9 f% W
c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
; n6 ?' z# U- X8 y. Y" t. Z6 z" I. J
Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],1 p; U* j5 a0 c6 ~& H5 U
c1 10.*10^(-12),c2 10*10^(-12)}//N];
$ O* ^2 Q$ O4 V- B
PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
+ a% X6 n" t8 R! v+ G
Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
8 z, N/ E2 P: x1 P' r: x7 V, F
,PlotRange All' u5 [4 y3 b( [$ x/ Z3 ^$ X
,AxesLabel {"f(MHz)","M(f)"}]

复制代码