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

小洋洋

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

. @' x7 m! k, v8 {' }1 p
/ m" {5 V. w7 S- w

小洋洋

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

guozihu1986

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

guozihu1986

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

山林隐逸

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

小洋洋

- F% m. Q/ O; F

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,+ V% E: |+ v\" K0 j5 O5 }! x4 Z
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};8 X( v4 n: R8 {9 \/ s
   
3. Current={Ii,Io};5 B: k! x+ F, b4 |\" Y  K& W: m
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];: a! q6 \7 W: f/ `. {7 J
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;( w' a* V3 ^  M* F0 v2 I; `
   
6. Print["H(S)=",H]
   + [% |! z0 w: E1 h; p' ^2 M
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   # X2 K3 Y6 B  j5 U! R& ^
8.    den=Denominator[H];
   ! L/ }: k/ T& e0 N$ V
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   # \0 a0 Y6 X$ N( U
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    ( p\" G' x0 [4 ~* K3 D+ V
11.    Simplify[{fp,Qp}]];9 y* t6 F$ ]: X% Z' C
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    / F* {$ Y4 M* q! W\" i  ~3 c
13.    num=Numerator[H];2 E9 s4 M# B. K' Z\" E0 ?
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);! d3 F1 A) v6 H
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    1 X2 L\" z% O4 B' ?4 }1 d0 A) u
16.    Simplify[{fz,Qz0 fz}]];; |  q+ Q9 {/ g) G0 c
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];  Q; E9 E' u, z% L
    
18. Klp=H/.s 0;) S9 Z8 s) Q+ x4 Y
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];7 ~! y% i\" K' k( B( B
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    $ [0 z/ m' N5 p; ]0 g0 s
21.    gm1=C1*Wp/Qp;
    $ Q7 _, K) G0 S5 o. g
22.    gm2=C2*Wp*Qp;0 t& Q  Z; B$ a% Y. f- E$ V
    
23.    {C1,C2,gm1,gm2}];\" Q3 X/ g  D0 u) [
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    + N$ h4 S% e6 Q6 M- c
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];\" g* D- T1 q0 @$ d) v: m4 Y( e\" K- c
    
26. Simplify[H]
    ( e\" y7 p7 z- p$ T( Y: z. C7 o
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    , c& J$ k! T8 c) q9 U4 w1 j
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];3 H9 \& {. b4 v
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],* x8 k! A4 q\" C' C5 l\" I
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    / Q/ D5 Y, y8 a$ C4 `
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    1 t' q- h0 U; _, Q& i  v
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    $ |& F# {4 p0 H\" y
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])# t: J( G7 \' s! x. G3 [
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    : r! u( z* j0 K+ G
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];& k3 ~7 l8 Z( h
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],# B, S9 l( V0 I* L9 c
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];. r( D! u1 F4 L0 d9 o# ^/ `) ~/ f
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];1 ]5 b. P2 }0 N! S. ~
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    - h1 ^, w- p+ X7 f! p( v
40. ,PlotRange All
    7 Q$ k7 f( [- v4 v& M6 H, x
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧

平凡之不凡

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