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

小洋洋

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

平凡之不凡

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

小洋洋

' d6 W) e3 P' {8 b5 O

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   / [+ y  U$ }8 h
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   9 f% g+ M3 i! X# ^* u
3. Current={Ii,Io};5 [& `2 X1 R9 \
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   / H/ _/ n0 b$ ], L) |& [
5. H=Io/Ii/.CircuitRespones//Together//Simplify;' q\" f7 e8 e* u1 e( y, w
   
6. Print["H(S)=",H]
   , d- u, F+ W( d7 e( |
7. PoleQpole[H,s]:=Module[{den,fp,Qp},6 H/ J) V9 r& R/ ?2 s& b0 o8 o* ?
   
8.    den=Denominator[H];
   8 x/ m' H6 V5 m
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);8 z0 \4 k* o! p1 {
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    # B: R1 N! O1 d6 E2 B
11.    Simplify[{fp,Qp}]];
    ; ~, i1 ]# Z8 Y% C2 x4 j
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    * l\" f5 i/ p( `
13.    num=Numerator[H];
    ) ~  D2 P, e3 _
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    / y- I. V. X0 O4 H+ z* I4 R
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    3 v' Z+ G5 S5 ^6 ]1 T
16.    Simplify[{fz,Qz0 fz}]];) C: j+ ~& O, p' w6 {
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];9 Q5 Z2 v( S6 }5 B' ?
    
18. Klp=H/.s 0;  g8 `8 v4 n: F+ m8 T1 w: S
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];: Q4 D6 U+ Z6 F' t: P% I
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    ) h3 x+ M) }. i! n
21.    gm1=C1*Wp/Qp;
    ( k% E8 j8 [& Y3 M2 S
22.    gm2=C2*Wp*Qp;
    9 Z0 n% V! }5 A* e. A4 I' O
23.    {C1,C2,gm1,gm2}];
    0 H! c\" R1 j7 `, t
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    ! c& K; j2 w& I4 s1 w. A$ \
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];1 ~1 A' {\" ^  B
    
26. Simplify[H]0 K. d- R; g: o( U
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],: z( e8 K# v8 |* D1 W, W1 `
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];8 s$ H: u; Q9 v  l( U8 x, K
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],) \2 {- Z8 k, z$ }
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]- [! g# }+ _+ T& K# o; @! N
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    ) i, W# O9 w# M
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]( h) x9 J/ _; p& {! x
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])8 R0 w# I/ p$ M$ k6 r. ~
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],2 {) I. W5 s* w# {& n1 h2 `1 p- T/ O
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];2 n. n' p2 Q. P/ ]7 w
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    5 g0 o, a: u! C
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    * Z% d2 X# y. ]/ Y
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    . U' N/ C\" s, `+ D
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}/ M; {/ a\" g7 C2 `+ I\" Z) P
    
40. ,PlotRange All& `( t( }- t. W3 M# l( k! g- b
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
0 d- S' V( L0 d) q6 w4 e9 B4 A0 N