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

小洋洋

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)"}]
; V" l* X& `! f1 Y

平凡之不凡

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

小洋洋

0 w; t# i) Y' k5 s

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   8 o, C, n) v7 q/ Y3 I' D
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};( B1 B% F! {6 y4 Q! t' Q
   
3. Current={Ii,Io};
   3 ~3 O- p0 g6 f
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   % u4 o# l, @3 i, \+ |% z
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   6 g4 M\" q+ f3 Z3 r\" O8 Q7 Q, p
6. Print["H(S)=",H]/ ~\" o3 B3 [' A) [  z& |% Y: m& }
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},! Q4 a/ M6 f* b+ K: {( v. ?8 i. ^
   
8.    den=Denominator[H];% J- e2 o$ ~. C7 A8 n2 @; Z2 {0 H& \
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);5 v6 A! z3 E5 `. `; h
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));) [2 P7 W! |3 _2 i# f
    
11.    Simplify[{fp,Qp}]];3 d9 G8 f+ s& n7 v  ~% ~
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},7 h7 Y\" V2 {' K0 Z
    
13.    num=Numerator[H];
    ) t6 S0 l: |8 V* T& U6 ~
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    6 z% J7 x9 I\" [$ h8 b2 ~# E
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);6 @) Y( ~$ O8 l4 ]+ K# m
    
16.    Simplify[{fz,Qz0 fz}]];
    ; E/ [8 z+ X: b- G3 X% c& ~  M# |8 q
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    # b1 ^9 x+ i- U! t
18. Klp=H/.s 0;
    * Z$ D/ {2 x5 y9 b% u! g
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];! j2 l) P4 t* n4 |! Q
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},& H7 s: u6 Z* K1 b. K
    
21.    gm1=C1*Wp/Qp;: V$ s/ c  Q9 s
    
22.    gm2=C2*Wp*Qp;
    9 I; d) G7 B  K% T* f
23.    {C1,C2,gm1,gm2}];6 s  d, Z2 l( \/ m
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    7 z6 j$ }5 s% I! }
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    \" b' j& ?; i6 K/ _& l
26. Simplify[H]
    $ k/ q# Z% p' I$ y9 W+ ^\" D. r: y
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    7 l# m, g0 S+ c: V5 O, W
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];/ `' Q( ^- @( a
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],  L+ n9 f9 G  |# Q9 K6 J4 z
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]/ `  Y; f' C0 c& C3 L. u
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    : ]8 }* G  o0 A2 \, Z9 p
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    7 j) N: F. E% o8 [8 c8 N8 X
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    + i# x. c0 m1 c: e& P( s
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    6 A3 B% o) M% i3 h' h
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    2 ^) ?- o& Z' v
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    7 _2 o; r3 Y' y* i, c4 R' E
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];, J) m8 U8 L& k! @
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    ; m, h. e; a/ G, e
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    , @1 e8 z) V3 I) c
40. ,PlotRange All( |1 u& ~; A. \2 `, g) Q
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
  L( q' T' E. E$ _

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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