求大神帮忙看一下程序的问题出在哪里，如能帮忙修改万分感谢

小洋洋

CircuitEquations={Ii(1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,  IoIg};Current={Ii,Io};CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];H=Io/Ii/.CircuitRespones//Together//Simplify;Print["H(S)=",H]PoleQpole[H,s]:=Module[{den,fp,Qp},  den=Denominator[H];  fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);  Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));  Simplify[{fp,Qp}]];ZeroQzero[H,s]:=Module[{fz,num,Qz0},  num=Numerator[H];  Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);  fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);  Simplify[{fz,Qz0 fz}]];{fp,Qp}=Simplify[PoleQpole[H,s]];Klp=H/.s0;PrintLabeledList[{fp,Qp},{"fp","Qp"}];DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},  gm1=C1*Wp/Qp;  gm2=C2*Wp*Qp;  {C1,C2,gm1,gm2}];{C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];Simplify[H]h1=Together[H/.{Q6.,WN[2*pi*10^6],      c110.*10^(-12),c210*10^(-12)}//N];Print["gm1=",gm1/.{Q6.,WN[2*pi*10^6],   c110.*10^(-12),c210*10^(-12)}//N,"mS"]Print["gm2=",gm2/.{Q6.,WN[2*pi*10^6],   c110.*10^(-12),c210*10^(-12)}//N,"mS"]h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])PrintLabeledList[N[{Q,W(2*pi)}/.{Q6.,WN[2*pi*10^6],     c110.*10^(-12),c210*10^(-12)}//N],{"Qp","fp(Hz)"}];Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q6.,WN[2*pi*10^6],      c110.*10^(-12),c210*10^(-12)}//N];PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];Plot[{Abs[h]/.sI*2*pi*f*10^6},{f,0.01,3} ,PlotRangeAll ,AxesLabel{"f(MHz)","M(f)"}]


9 t: A5 R- T7 R( M

平凡之不凡

我建议楼主把问题格式调整一下，重发一次，有人会给你解决的

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   ; F% I% C# N: L; a
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};9 Q' R! _1 z! M) o# z
   
3. Current={Ii,Io};9 |9 K3 o1 E6 m2 Z; V
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];4 T$ E6 _, o) N: u$ W3 m8 K& I
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   6 c) U' n( f: c+ Z1 \
6. Print["H(S)=",H]
   \" ~1 @- e; b3 Q4 t% f
7. PoleQpole[H,s]:=Module[{den,fp,Qp},+ f5 y9 O% m& \, x( }+ J( B/ p
   
8.    den=Denominator[H];2 ^, I) K# A- J# }
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   5 e5 M' F) [% G1 Z4 z6 \. g
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    ' h6 J- @/ D4 g) }
11.    Simplify[{fp,Qp}]];! _\" j$ }% @, K5 Q% `( j+ Y
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    2 I2 B# U0 L! V\" u5 L. ~, ?
13.    num=Numerator[H];7 E4 K) t- M5 `7 b
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    6 M) V. b1 x* v9 Y, }
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);# H+ i& l8 j, Q
    
16.    Simplify[{fz,Qz0 fz}]];7 U6 |% T2 K- K\" s$ Y  x
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];) }5 r) q. r! b
    
18. Klp=H/.s 0;
    $ L- l) |( t3 k8 [1 z$ V\" S
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];( H* w0 Q2 f0 ]1 y
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},5 a2 ^& y! M, C
    
21.    gm1=C1*Wp/Qp;; r: m, d\" y, x* G: F
    
22.    gm2=C2*Wp*Qp;
    , b( e$ ^2 v- n- v  r9 b
23.    {C1,C2,gm1,gm2}];4 g: G; w  d9 m  `
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];+ t, W4 d+ h. I4 ^\" x
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];+ B: C0 T& C: A! x: ^: a7 J' q
    
26. Simplify[H]
    ; a0 i  \, n5 Y) j& P. \
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    , o6 z' ~6 t4 q% a2 d) H% N) L3 @' q
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];, \: q9 x6 X; p9 x
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    + B2 C* G0 u, t1 @
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]) O% P/ g7 @6 w! }) A& Q  s
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    ! d( V: ?' R' B4 q+ M
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]0 u* x* Q' q: J. O3 G
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    3 ?5 |& c5 V; A+ C\" m) t; I' E
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    4 c; M' X9 n9 c) k+ W( d; W
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    , @( x8 i8 U; d
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],0 p\" x* u8 j\" f; C; s3 w. o
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];. i  j7 S$ v: w: w9 L9 L9 ?+ g* d
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];5 S2 B* ^0 @8 j\" s9 c- X
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    5 w2 F( K) N$ v2 @
40. ,PlotRange All0 N; ], q( A9 k$ C; N7 s\" X, A
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
: l2 A/ s/ j: p; ^

山林隐逸

同意。。格式改一下。。然后说明你用的什么语言什么软件。。什么环境

guozihu1986

建议你最好能把你的原始.nb文件传上来，这样别人才能帮你看，还有你这个程序的目的是什么，想要达到什么样的结果，你也没有说，另外，关键部分的注释也很重要，这样才能让别人帮你，一点个人意见

guozihu1986

建议你把原始的.nb文件传上来

guozihu1986

把程序的目的说清楚，关键部分的注释也要加上

cnmover

看着眼花啊

guozihu1986

看没有注释的代码，比自己写还要累啊，楼主这个程序是想做什么呢，能不能稍微介绍下，就放下这么一段代码，太难看懂

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教