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

小洋洋

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)"}]
! U2 b8 U. q3 R1 j$ O
" d6 U1 P; D. b$ _
1 P6 m) V- R& h4 u4 ]8 Y% \; r. ]

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
! y: _8 ]3 q; P, z' b, D

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

guozihu1986

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

guozihu1986

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

山林隐逸

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

小洋洋

( C, z8 E7 P6 `  G( H

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
     b' }9 R6 e& H0 x5 {) m6 t
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   + p+ H- h% f8 X
3. Current={Ii,Io};! @+ P& K, O; }& K
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   9 X# S( e$ ^& ^6 R
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   ( s3 B! V8 I) b% w: e
6. Print["H(S)=",H]
   \" S. k) f% d3 u: X
7. PoleQpole[H,s]:=Module[{den,fp,Qp},6 S; ^! p7 e\" c$ t6 T
   
8.    den=Denominator[H];2 p* i% y% m8 L7 c) f3 p
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   , W; _, Y$ Z$ g0 j- |
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));  H; u) z; l' |2 I8 A3 ^
    
11.    Simplify[{fp,Qp}]];
    1 P; r5 Y* p  L5 l9 b# k
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    ( K) s# Z% j' ]3 w, H. }* c0 h
13.    num=Numerator[H];
    ( s\" ~+ E! g- m0 t1 U
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);/ v; a4 S/ r\" O  T
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);1 a. t\" ^+ T( Q  }- h' |
    
16.    Simplify[{fz,Qz0 fz}]];
    3 s/ P, \' q\" Y7 ]3 i\" o6 E
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    # ]8 S; a! x* @$ z# ^; |) E
18. Klp=H/.s 0;
    2 k5 p0 `6 J& D( J6 I  K\" a+ r& \
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    1 ]! q& a  A' [4 u0 N5 q: g% ^
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    , L( ?/ I, }) X8 C# {
21.    gm1=C1*Wp/Qp;/ ]( B0 n2 h$ o  g! H8 k( d7 X/ t( O
    
22.    gm2=C2*Wp*Qp;
    9 I% W. A1 x5 \! M4 n0 C7 R
23.    {C1,C2,gm1,gm2}];6 l; G& f8 ^! h% b
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    ! c, Z, a' N2 V: o; @* g0 q
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
      V+ R' }! T& L: e9 z+ K
26. Simplify[H]
    ! H5 x0 x: Q7 e4 L0 n
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    , x/ ]; w7 t: Y3 ]; D' ]\" W
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];8 H( @0 P9 K, h4 V  P8 f7 W
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    ! R: ~* o- \7 E4 c( \) [# Z2 S
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    \" y( z$ E2 b  Y
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    ' V/ |; _4 M; ^! n
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    $ b2 N, g9 @+ _5 K
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    2 o4 {( f5 L8 f  Z3 K
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    5 B* r5 n3 G9 L  s
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];4 j* b8 e* O0 M
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    $ S8 k: t5 w; t$ d1 |( Q3 V
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];8 r# G/ _4 x0 C
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];, V* |+ ]( e1 y8 p
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}) g2 V\" K. O% [
    
40. ,PlotRange All5 d# \\" w  v* Z* S% o
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
8 _$ Z, P! K# F9 ^& i' h" K. v

平凡之不凡

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