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

CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig, Y' e- t9 a1 l) t# O& ]0 x6 r
Io (gm1*gm2/s/c1/s/c2)*Ig};
0 p$ j' h' u- K# c5 C6 Y7 r
Current={Ii,Io};0 G/ b2 s7 ?4 J$ h5 j# Z# R
CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
$ H( K8 J+ K8 z8 H, [
H=Io/Ii/.CircuitRespones//Together//Simplify;
2 Y' z7 T3 n6 g+ a5 O
Print["H(S)=",H]4 E- b+ ^# a: m* ]
PoleQpole[H,s]:=Module[{den,fp,Qp},
4 N- W) ~0 v+ S% j$ r
den=Denominator[H];. y! R$ X7 p* [7 Z% j* a" e$ x) ~9 d
fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
- k0 }. e! q2 L" s+ F3 R
Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));! J& e; C4 q7 V S
Simplify[{fp,Qp}]];5 ?. r$ ~; W2 g- K
ZeroQzero[H,s]:=Module[{fz,num,Qz0},
7 E( o9 `. G _% v+ z6 R/ {7 z
num=Numerator[H];; O# ]" e% x" ]& `
Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);/ F0 ?4 ~: q; z6 n! t9 S7 c
fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);- B, h6 l1 ]: w) u/ i- X
Simplify[{fz,Qz0 fz}]];
! U" k1 M# v2 X2 [' [# S9 U' S7 B
{fp,Qp}=Simplify[PoleQpole[H,s]];
. |$ I8 T9 t# O `$ i# B
Klp=H/.s 0;
/ n1 [4 a9 {+ m v1 P% \
PrintLabeledList[{fp,Qp},{"fp","Qp"}];
" M( _1 ?' i3 X+ w1 _6 x4 r# R
DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
1 _- J8 C# ?( ^: @& u4 Z; v2 G: h
gm1=C1*Wp/Qp;
( S, `9 c7 N \+ i4 Q) S& N
gm2=C2*Wp*Qp;
7 Q4 p: t; h5 x; `. ~ Z e
{C1,C2,gm1,gm2}]; e6 ^) D) U) ?# U3 _* O
{C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
! b0 q( n2 J: o; Z/ C+ \( U) W9 n3 F
PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];+ f& J- i7 c, g6 a; h" n d4 H
Simplify[H]
; O+ G& L. j( N) o3 F, |( ]
h1=Together[H/.{Q 2.,W N[2*pi*10^6],% t( V8 |8 W( e9 J0 F
c1 10.*10^(-12),c2 10*10^(-12)}//N];
% C7 w/ e/ e( A, R
Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
% M. x. t3 Z7 u. @" e3 z& z
c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
" r- x, E& O# X% x9 R4 k
Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],) w# I4 d/ q) k; F0 e
c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]9 s, j5 M9 |& K$ n
h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
7 ?; f; O$ y s5 n( s
PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
2 t$ g+ D& @# B5 E+ d4 n+ n
c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
# \ @0 g" i% Q$ y9 Z# C+ g% O- R
Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
% P! `- g( y1 t& W
c1 10.*10^(-12),c2 10*10^(-12)}//N];, Z, p0 z2 o* [$ S9 _
PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
: ]+ R# h4 G9 Q1 K9 A7 Z
Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
0 P! l2 M$ V4 q- G$ G4 ?5 b( ?
,PlotRange All6 }: `3 Y5 Q& D! I, v' c
,AxesLabel {"f(MHz)","M(f)"}]

复制代码