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

小洋洋

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

! m$ {; _5 |$ a: R. {
4 E# U, d5 f* U8 _1 r) q; o; `

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   / |8 b/ E$ z4 S6 a
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   3 e1 A( t* q  D& q
3. Current={Ii,Io};
   6 W& Z9 V1 v  Y1 ?* \
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];9 H' X3 u  N$ A( k, ?* W
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;' H& z+ |6 w/ S) F  P
   
6. Print["H(S)=",H]
   0 Z4 B- O. l! l5 r+ D
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   ; l+ {1 H- H7 ]) p( c$ n+ ^; ~& c9 o\" f
8.    den=Denominator[H];
   - W\" E; }. [' R  M$ X
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);$ p% z  o2 l' }7 P. Q
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    3 Z8 q$ s  ]9 ]5 n' S
11.    Simplify[{fp,Qp}]];. ?) T+ y& c1 N3 `  o) K2 `. Y6 I
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},- M  g! Y) G8 ?2 u- l# A
    
13.    num=Numerator[H];% Q& P: E$ U' k5 {* D
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    7 O1 U; S8 ^1 X8 c' ^
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    8 Y% Y, m$ C3 T% U8 G
16.    Simplify[{fz,Qz0 fz}]];
    / ^2 j4 U/ Z5 [% I3 G
17. {fp,Qp}=Simplify[PoleQpole[H,s]];3 I4 t) ?$ d5 }3 R% W: f/ I/ g* a
    
18. Klp=H/.s 0;0 C2 D* }; m( I' H5 Y
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    8 |- E/ {& T' v( e4 o8 R( P. e
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},, U8 s* z4 l& u8 A. A
    
21.    gm1=C1*Wp/Qp;- l) K- Q$ I- G, Y9 O1 \7 g* _# Y( v+ O
    
22.    gm2=C2*Wp*Qp;% V: c$ c6 M6 O\" i5 J$ s8 Z
    
23.    {C1,C2,gm1,gm2}];9 M# [& J; a- F0 {. c' q% r
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];( w5 s, u\" X7 e9 [5 {
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];: @: k! e8 a  p  w! s
    
26. Simplify[H]
    ' a) A% b( t2 T5 o! D
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],, {# V' |( f& a
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    \" Z$ K# s% j( n/ Y\" t
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    ( h+ l% K6 p* k$ w$ U
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    & F0 j8 s7 a6 o5 I
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    6 T( c: j3 }# U2 h5 h6 `
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]$ \# m9 c* S! N9 h* i, ?
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    9 u$ j- _. w% R( C* d# S! y
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    5 h# L  @1 U, s: d
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];( C: |+ M6 @8 s# [7 V' O) M
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    0 z( D3 H0 ?3 ~/ r
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    ! l6 y9 Y5 L\" T- D+ A  A4 s! p3 b
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];& Z: ?+ ^' Q- Y. I3 z; Z4 e% W) C
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}% h3 C% Z2 R) C. ]* S' D+ J
    
40. ,PlotRange All
    & \7 ^+ f: p5 |* x/ s
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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