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

小洋洋

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)"}]
5 t9 y- g/ T% z: s

: E4 o! {1 E- x! x

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,  n) t$ O, k' M\" _2 z' R1 t, u; m/ K
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   2 C\" b3 o5 d% n; i
3. Current={Ii,Io};
   % M6 w! w4 u8 c' H1 A
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   : e/ b! t3 v  x6 Y5 [$ B9 O
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   ! d  b) i& \7 k2 }( B
6. Print["H(S)=",H]- [  b' k* X# S  V$ e2 ?
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   0 W! _; E' T7 ?% H
8.    den=Denominator[H];: H' P\" W1 R! \0 N
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);  K* p+ ]/ e' D6 J
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));) f1 k( u  D$ a# E4 j; }: o
    
11.    Simplify[{fp,Qp}]];: L' J/ I9 D/ d2 ^
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    1 Z0 G+ V' x\" u. I5 W, {* C
13.    num=Numerator[H];2 x5 v+ g1 ^& U& h' J% X* w
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);- r# l& s2 H\" o8 b: }) e
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    $ P) T8 r; i/ C0 C4 s. _: _
16.    Simplify[{fz,Qz0 fz}]];; }- A8 ^4 ?  T, {2 a: b
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];, w5 G% ~) K% P5 ^/ J2 I0 e
    
18. Klp=H/.s 0;5 [* p7 A, U) b1 y
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
      H0 Y) {7 ^( Z8 h  D3 T+ B
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},0 V/ W0 }# z3 |
    
21.    gm1=C1*Wp/Qp;9 W! y5 h; R' |9 R\" {9 \; r( E. V
    
22.    gm2=C2*Wp*Qp;! m$ H2 O/ W9 @: Y% m\" l- }/ m- Q
    
23.    {C1,C2,gm1,gm2}];3 J1 J& o+ B3 {, ?, |# C- |- g8 I
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    8 m- Y\" i4 l  E1 Q6 P9 f4 m
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    6 d* i4 {\" Z0 ]! u: p+ p
26. Simplify[H]
    \" z) L& l9 }. A0 Z& i
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    0 C0 X* S6 r( R7 Z
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    - u: j' d( v\" D% N6 M0 [
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],7 q( o) e9 D% F% r3 {
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]+ h3 l( ^9 v% ?5 ^+ k1 R
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    0 K5 n* c1 P' D% F- J3 ~* w
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    & |% d. Z8 R8 d! i! t# @
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    ) c; B: w* J, p! y2 K# i; I
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    / U5 F4 Y' ?; m4 y
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    ! K- q; a% v8 ~
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    + G8 x; ^- g& ^
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    6 p- e: [5 e  I) j9 K
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];$ P3 ^( @& L( w$ D\" e$ y9 h
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    1 D$ L8 O; o3 J& |9 G4 X6 t' i9 a
40. ,PlotRange All
    ( g- c4 `1 a* l) O/ u
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
6 s/ c. d+ x! c: N7 B8 ^$ w