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

小洋洋

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)"}]
0 s* q" Z9 T$ O% B" c) |
/ ~2 L6 B" ]5 [1 ?0 L

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
& j( f( [5 B# T- {2 S# @3 S. C

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

guozihu1986

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

guozihu1986

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

山林隐逸

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

小洋洋

0 K. |6 p( f6 w5 c" g. d

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,+ {5 e4 `& \1 P9 I9 T! w: Y
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};7 I9 U$ S9 c% _5 l
   
3. Current={Ii,Io};
   ! q/ K$ a% ]1 \* U& O
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];0 ~1 j4 J& x6 |: N
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;& E) O) a$ M8 ^( m
   
6. Print["H(S)=",H]
   ) W3 R\" D% Z: i9 G
7. PoleQpole[H,s]:=Module[{den,fp,Qp},4 o# j4 {7 e) Y; k. f, X# {4 r
   
8.    den=Denominator[H];- u. `$ `\" T/ C! ~8 e, U  ?7 \7 U. ?
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   ( H2 n2 n# m1 A7 A2 Q5 T* Q7 d6 _
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));) @8 C\" L3 e7 r3 |6 q) K, v; g
    
11.    Simplify[{fp,Qp}]];5 @% u% |+ N' s( L$ e& |
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    $ h' q) L& W9 }; T
13.    num=Numerator[H];
    0 R7 N# d) h% G8 V: q& c5 D+ x
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    , E6 ~5 [% E) d\" s+ O
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    5 L; z$ n4 c& r) E) Q; P7 Q
16.    Simplify[{fz,Qz0 fz}]];
    \" Y8 y, \# I\" e( q0 Q1 L# A
17. {fp,Qp}=Simplify[PoleQpole[H,s]];) o) g0 S, ]+ r! K' _' H
    
18. Klp=H/.s 0;: Z' T4 m; H) w7 Y) c
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];! @: F\" b) X+ Q' h) x
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    ; @0 w9 ^/ j0 @
21.    gm1=C1*Wp/Qp;
    , s  P0 V# m) U0 L7 o8 |0 m  c
22.    gm2=C2*Wp*Qp;
    ) m7 e$ w4 `# `% ~1 a
23.    {C1,C2,gm1,gm2}];
    6 V9 Y( _* d\" p2 ?# `
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    8 t. V$ ^* e0 B1 t6 O: ~
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];2 o7 L& {+ A; m: P
    
26. Simplify[H]9 ?, h3 P* j# h6 {$ i, m) a
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],4 T# b7 l# g6 q
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
      v- _\" a/ M& q8 v
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],0 O9 e& a$ ^1 `7 N% S: h\" g9 C: x
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]! H& j2 E4 _& t
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    9 A: w: P; S' M9 b\" g( A: J
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]' v) j/ h+ X' k; m# y
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    ( C$ z. {0 l  u' \9 {
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    6 q3 R3 J; S* j
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];, l2 M; t, b' F. D4 C; o+ W0 ^
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
      t: s7 C, j# S: E% g9 Q
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];# w# a' ?. b! g& r: J
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];\" X  r) \, B; N8 j* u- B* K
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    0 _4 r. U1 P$ Z3 N- M
40. ,PlotRange All
    . f' V9 r5 X5 s7 L( w! j
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
0 ^& I5 |9 r& Q6 M* Z( I9 F% q

平凡之不凡

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