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

小洋洋

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

1 K' r# @8 B) d; s; x! V/ D6 d

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,$ b7 r4 J* y6 O9 O; m3 r& T
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   5 R1 W/ `5 _8 N- A# H\" c# U1 h2 B, `
3. Current={Ii,Io};8 ^1 e( c( q4 k' O4 g) [  l/ F3 S
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   $ B  Z8 D2 O1 @9 H/ i6 u) s
5. H=Io/Ii/.CircuitRespones//Together//Simplify;  F. a\" z: \  x/ @7 y
   
6. Print["H(S)=",H]+ y4 B( e6 p. \0 A# y) I8 _
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},! C  O. s2 z, R/ u\" E' H
   
8.    den=Denominator[H];
   - _: U0 Y* {# `. A: W
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);' M, ?- S! J! X/ P6 v5 @
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    8 n5 C( d3 G$ J
11.    Simplify[{fp,Qp}]];
    3 D: Y, u$ A# ^$ K8 E7 P7 a
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},\" ~% p, b! _0 T
    
13.    num=Numerator[H];
      S5 \4 ]1 ~9 b& y+ N! \& _$ u
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    7 G4 Q' \) f\" e) q& H
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);2 Q. Z4 ~: c* F1 ]2 i9 f
    
16.    Simplify[{fz,Qz0 fz}]];
    6 z3 E$ @8 A. i1 n# T1 z
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
      g2 U7 g  @3 S' X
18. Klp=H/.s 0;( Y- g7 O! T5 b+ q0 S( M! j4 U! |
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];9 x9 N5 }  c- P, R2 F8 Z
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    \" [2 R$ K+ o9 V
21.    gm1=C1*Wp/Qp;
    $ w0 ^6 S1 N# z& X1 w
22.    gm2=C2*Wp*Qp;
    & g/ i6 f: u/ x1 V8 }9 o3 c
23.    {C1,C2,gm1,gm2}];2 }1 S; s7 i* ?: ~- T% M
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    ( y5 n7 W7 r* G: a8 s$ P
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];7 g; `( V1 d$ ]0 P! E. a$ w
    
26. Simplify[H]2 P0 t. d8 p; C4 ^# L2 h
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],' a  D2 K  F/ J& t, b) P( o9 ^8 r5 |
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    ) }& H, ^& h/ ?4 q! m
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],  k! `& b9 {. M% Q8 H\" T/ ?  ^/ u
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]$ r* ^8 a/ ^0 k. a/ M: L
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    0 b( |/ Z. Q$ z: W
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    - ^6 f: [# U9 q
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])/ K  \, x* C\" }; P9 a1 B
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
      K, ]0 O$ C& O' L8 Q! U! B
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];- V% o8 y: ^- O, s3 L
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],: R: t+ x( }# C  h3 S& R  q. I0 k
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];$ E9 Y% k\" U8 W  `5 G0 b8 X8 s5 T\" w: `
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    7 f1 N9 s2 F/ D8 }8 C) }# s
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    : W\" n, T) Z' V
40. ,PlotRange All
    7 f3 R1 ~; M5 Z. y! M
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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