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

小洋洋

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)"}]
! c) c1 ^6 i% t, k/ n
. ~' ]- ]; z) W% k$ {0 w( @
( n/ ~4 _5 p/ J( {; |, p

平凡之不凡

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

小洋洋

: ^0 z; o5 l& v: |9 c

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,+ O( s6 I6 j) Y
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   8 w# w/ m  t: W% e5 R' }- P  c+ k- t
3. Current={Ii,Io};- I$ U/ k3 ^* @; x. d9 h
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   % ~6 \0 u( ?' T
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   ) ]$ R: f  Z8 i0 }: ^
6. Print["H(S)=",H]
     x! ?, c' o) r4 j. j
7. PoleQpole[H,s]:=Module[{den,fp,Qp},2 B* O7 ^  ?% O$ `
   
8.    den=Denominator[H];
   9 T3 Z% z# Y  G  E* {9 w1 ~( c$ W
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   \" }) A( ]: w4 o5 \
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    * [7 n\" s! S/ `, P  m! u
11.    Simplify[{fp,Qp}]];+ g) Z5 g. d; Y1 e
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    ( r+ s; a. k) U5 Y9 A/ j
13.    num=Numerator[H];% Y, |8 g  Y6 c% B: {
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);  K3 ^5 t7 L, @& i9 a* V# a% x' D
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);, O/ U5 _6 Q( y$ \3 ]& G
    
16.    Simplify[{fz,Qz0 fz}]];7 A+ Q  v1 V7 g, Q, r7 C
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    + \6 f6 C+ N0 y\" J/ w2 s\" Y
18. Klp=H/.s 0;) f. \! Q2 k. C4 u7 Q, j
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    ' K; @2 I* D& Y+ w1 P) n
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    3 y1 b7 A9 |% V0 n# a+ r! g' P
21.    gm1=C1*Wp/Qp;
    # `\" B5 ^; }' g: |
22.    gm2=C2*Wp*Qp;% G% ~: o0 ~0 K/ s/ R: F
    
23.    {C1,C2,gm1,gm2}];  h1 F9 c0 i8 v/ ]
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    5 \; O* ]! X  @% T4 k
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    7 C\" x5 K' V- U1 P3 c+ k
26. Simplify[H]- X- d2 ?$ S& K& l
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    4 y5 V' b4 k2 n: b3 l9 l
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];+ H1 v5 R  K$ @
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    - @/ ^5 k4 p# M. J1 {% H\" G( r
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]# B5 g' |1 ~2 U% a\" f: {/ K4 n9 D
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    # m% Y% e' v% N6 |) f0 X
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    & q/ u/ ~1 j* f
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    ( E( K- C% F9 S+ x  v. c/ Z
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],\" g) n5 G3 }  \\" N4 A
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];\" u$ p, w+ K\" N$ l3 K) z3 S! a
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],+ H# P/ r5 O\" S; H/ h\" d
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    , L3 Y8 z/ j) V* n. S6 |9 z& F  _9 f
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    , s0 k! ]& p) S& v
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    , C( w' f5 E, f3 h
40. ,PlotRange All
    ) R& m  j& y) F0 J\" q9 d
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
9 B" d; [! ~# Z7 L' a7 j

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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