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

小洋洋

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)"}]
  o5 }. x* r: P( F
% j, ^8 W7 g$ M  v4 D# K- q
& _6 S+ B3 D: Q! M$ I/ G8 E

小洋洋

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

guozihu1986

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

guozihu1986

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

山林隐逸

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

小洋洋

0 t& {$ ]: e( `6 {0 c1 f, Y) G

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   5 {. m  A4 k2 s% Y/ P
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
     F+ U+ M3 |. `2 W( @# P0 t
3. Current={Ii,Io};
   4 l/ S- ^. g1 H. W4 n, [
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   ; B5 U) W8 }# G: H6 b9 ?
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   , I- A. O  {9 \/ B* R7 }
6. Print["H(S)=",H]
   4 @: L% o4 l. D% _: B# I5 Q, m# S0 T2 S
7. PoleQpole[H,s]:=Module[{den,fp,Qp},3 O  C/ i# E4 H7 |2 H$ @
   
8.    den=Denominator[H];# i* M0 s\" W; k) E( l6 D9 Y! f+ m7 @
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
     @& o6 Q9 W1 G/ n+ b
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    : w4 v- S- L# J0 ]
11.    Simplify[{fp,Qp}]];
    1 r, O7 b& i9 j0 t$ K* V
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    0 d9 s5 J\" F- j* N4 V& G\" Z
13.    num=Numerator[H];4 P4 R4 A6 Q3 e- b; F, [; u
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    & _. c, ^! p# x4 D7 B
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);  p7 U5 X/ s& M0 n7 S6 m
    
16.    Simplify[{fz,Qz0 fz}]];
    % c1 x# t! a, [$ C9 `
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    - q( w3 v5 x/ b5 Y' x5 z3 H
18. Klp=H/.s 0;; V2 p* G\" ~* i1 E
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];( m% {3 d4 P7 m$ @% }) j: ^8 l
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    * f! Q% I8 R/ N5 {: [& ]
21.    gm1=C1*Wp/Qp;
    ! ^\" n5 Q* y+ D' W7 m1 l% a
22.    gm2=C2*Wp*Qp;
    & c1 X$ K7 E1 O; n9 ]1 F$ W) R* X
23.    {C1,C2,gm1,gm2}];6 p5 T/ j# P* O2 P, `
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];, f5 Z2 E1 u* [5 f2 D+ J
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    7 |8 L) ~& V4 j
26. Simplify[H]
    / N) w% b6 T0 R\" c$ g
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    ! ?+ Z6 N, P9 d) l
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];0 D7 P  Z# \( B$ w. g% d3 f& z! V
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    7 R) r9 e9 Y, |& u
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]1 D) D8 H* y6 j+ P/ K
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    # k! J& _& [! p4 q% T
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    : w( E, {% `- `. R/ D& U( k
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
      d' b: M: U. g, V% f
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    ' ~0 _, v* C) `- `
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    : @$ h9 |# T& `- v8 g
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    - }+ _7 M0 O; ?2 Z4 T9 z
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    4 j/ K; F! G\" q: p9 D
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    . _) ^8 L& x. s
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    9 V7 I. Z5 F% b/ ]) j* L\" A3 Z1 _
40. ,PlotRange All
    ! d% f2 h: T! y' j
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
# U8 N& I- h1 T; e5 y! L

平凡之不凡

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