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

小洋洋

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)"}]

4 K+ |3 K4 K, M( a3 {. C# o3 L

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   ' M$ b( x6 p0 D( [2 Q% p
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};8 v# P! r) ?; U8 A& w0 w1 m6 V
   
3. Current={Ii,Io};
   - l1 O* A+ t) [* ?\" F
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   % I2 ?5 ~; ?5 u* F/ t! m) ^
5. H=Io/Ii/.CircuitRespones//Together//Simplify;5 T- v& [9 M$ M' t' @3 A3 }. Y\" M
   
6. Print["H(S)=",H]
   - n9 t\" I* ?- j9 X8 i1 [5 p. n+ q2 m/ s3 h
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   1 u7 D) g# r# g6 F
8.    den=Denominator[H];
   8 T0 |* G5 L/ M6 B: c2 x3 p
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   7 [' k. o! U6 _0 E6 R; m$ {
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    ; a4 K1 M! s' s! Q4 F, t
11.    Simplify[{fp,Qp}]];5 \3 o9 U' J9 W% b& \7 M
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},9 L* V2 ]6 L% v/ l, h
    
13.    num=Numerator[H];
    * p: D$ X3 @% J\" `; T9 d
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);0 s5 v+ E2 z0 |0 y
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);$ T7 P- v7 F0 B5 I) T' @; e
    
16.    Simplify[{fz,Qz0 fz}]];
    . t* _* D$ G) M6 l
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    4 i3 M1 f/ _8 Z5 k; _( A5 w
18. Klp=H/.s 0;
    % G5 W& J: q' m1 A6 Z( `+ h' m0 a
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    / y7 x7 `# \0 O4 h; X5 U
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},- a5 ?$ U3 I) n2 d. W3 g9 `* O% |
    
21.    gm1=C1*Wp/Qp;
    * i* T9 K' L) v; J
22.    gm2=C2*Wp*Qp;, }$ z# b( q0 H' o
    
23.    {C1,C2,gm1,gm2}];2 j\" z. @1 W* A# B/ i* b9 \
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];; w. G! o6 s8 V! J& f/ G
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    7 M$ r& n# |  F
26. Simplify[H]
    8 {; Y- m7 p, z
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    # }. W+ _5 }7 O
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    $ q; o  }5 y\" H% d0 U* s: `0 z4 P
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
      g1 n0 J\" P9 m% _) |* i; a
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]+ Q& U' ]/ W' \# b* s
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],, |; Q/ p* a6 w# s# r! R
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    + _4 N* a$ n, X( p& I
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    8 I7 G- {9 l  E6 [
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    : `: y/ f7 z2 [9 u
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    6 x. H, X1 m' t: e2 w
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    5 Z) }' m* F, O2 L
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    7 L$ L$ a2 ^6 c8 k1 W! s
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    - _\" Q7 H8 W6 V! [) Y
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}4 B% M& q6 _4 a  k+ o- Z
    
40. ,PlotRange All
    # B$ b7 O$ r  w5 T6 H0 ^
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
5 y1 M: o5 U9 R# B

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
( ]% e, A- F8 [+ i. ^/ B, E5 p