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

小洋洋

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)"}]
1 h0 w: |' U# G9 w* U% V: \6 I6 T
( S; ^; P8 y; }) K

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,. v7 l- m. y: s; U5 h; y6 _
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};5 ]$ u* h4 n0 P& P4 _4 k; @
   
3. Current={Ii,Io};. F3 j8 K6 W. w2 `* c0 [) H. |# j
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];/ J5 B+ p9 N9 H9 {/ Z  }
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;( @1 I% Y. \* F% T
   
6. Print["H(S)=",H]3 F3 S5 |( ^* G
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   8 i$ Y6 Q' J; v
8.    den=Denominator[H];5 ~4 T5 D\" A; a0 [% q
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   . h% m7 U# `# ?4 c2 I/ w0 p  K' p
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    \" C' t. e+ V8 i
11.    Simplify[{fp,Qp}]];- Z  b- B7 J  I9 @- h
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    + F% _\" H- ?4 J% M0 A. {7 q
13.    num=Numerator[H];
    3 ^, m( m$ O) |
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);6 r9 [; m1 b. t$ K1 ]8 x; g
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    2 o& \; T9 K4 ~
16.    Simplify[{fz,Qz0 fz}]];
    & {$ f% @3 q\" S4 j5 P' F
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    6 w6 v/ Z% r\" |/ |6 h
18. Klp=H/.s 0;
    1 J4 S0 R0 l$ v( u& M' g  a2 _; K
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];$ ~* ?9 E& Y8 f2 s\" A2 [
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    ' l& Z# q( m; @7 P) [7 \' e
21.    gm1=C1*Wp/Qp;
    ) ?9 i. i\" a9 g) ]; i
22.    gm2=C2*Wp*Qp;
    5 [1 ~1 w$ z& q  s
23.    {C1,C2,gm1,gm2}];
    $ Y: V4 `! k9 h  {/ ]3 T& p
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];/ f' }/ O: F% D' O4 z  K1 }
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    4 ^2 ?  c8 y3 l
26. Simplify[H]7 \# J+ G4 O2 }
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    7 _4 ]- T6 x4 e  J5 |- |0 K) B
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];/ O  U( B3 s/ Q. d$ I2 u
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    5 g3 V! b! c, P- D1 _
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    4 e1 Q* `; V% V2 O( t1 b
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    / }# T7 n. R' l, U6 z
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    3 W8 u( c8 f- X; r  A8 Z' G2 G% @
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])% M- G$ P% O. C5 V; h
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    : h0 B* H4 I. b% x8 b5 U9 Q
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];6 w8 p3 T6 _' d. \! {
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],- l+ E: n9 z9 n- K! }& w\" c6 o
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];5 Y. k# C7 Q# [( n7 u\" a2 K7 y/ \
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];  U0 C& M1 K\" ^4 [
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}\" ]6 o\" U- e% x, t. }3 z
    
40. ,PlotRange All/ i4 q5 X6 V$ k
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
: d9 w& U: \! C+ h( a

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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