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

小洋洋

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)"}]
5 m* ?# D5 ?* a! G% O$ x
6 Y& @" w/ x' J
* f0 T: c$ C% m0 c( g

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,2 i: p/ z- i- M) D: {5 b
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};! l5 s, X; f6 u7 F1 X; g9 u
   
3. Current={Ii,Io};
   * F) F/ ^: m+ d5 _/ ~9 O5 t
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   % Y9 D+ ?% q( ^- Q
5. H=Io/Ii/.CircuitRespones//Together//Simplify;5 V, k: W1 y5 m$ e; W# R8 o2 \
   
6. Print["H(S)=",H]
   : Q$ ]8 N. O4 e8 I+ J+ }
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   0 d# _1 H* X2 E6 S4 X) j\" E) V
8.    den=Denominator[H];4 I9 f1 \\" ]& R
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);2 E. y9 e  J$ J( f1 a$ |- n( M8 j
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));, R' ?# b1 U4 E- X& s
    
11.    Simplify[{fp,Qp}]];
    * \: \\" L! {1 N/ v; l
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    ! T/ ]% }- I, T3 ?
13.    num=Numerator[H];( V- t# `$ l! `+ G1 S: G1 P4 b1 j3 }
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
      J2 \7 q4 B& n4 q4 g
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);3 Y. E. `2 j* @
    
16.    Simplify[{fz,Qz0 fz}]];* g# s0 u% g9 \- e
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];, u  K. |# G4 S0 C; W! ?& s
    
18. Klp=H/.s 0;
    ( G. N3 ^. {2 o, N7 J
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    % x( `: ~6 o# d: K2 z
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    % C. r  B7 x' r+ h1 K+ v$ I! ^/ ]
21.    gm1=C1*Wp/Qp;
    / r; \+ k* Y4 H# j6 m
22.    gm2=C2*Wp*Qp;( L; H; m5 j$ B\" A- B! N
    
23.    {C1,C2,gm1,gm2}];
    : l: @3 I0 I# r\" D4 v* i0 f0 n
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    # ^: k# s/ ^$ I; r; \
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];1 g, a. Y3 u( x) y! I
    
26. Simplify[H]
    + s& v6 |$ U4 g
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],. p! m* S' I, {. h
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];- p6 _0 k8 B$ r, F7 g+ b
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],; t0 ?0 M( F; t\" ^+ g8 m
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]- m( r\" ?+ c9 e, J$ t5 s
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],  Q9 b4 C( ^% a* p& d
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]( ]8 @5 }+ k/ m& H* ~) b: M
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])7 E3 s& y8 E0 C: m2 ]
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],' {2 @; O8 g+ b/ g8 v3 Y; D
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    ' s6 B) a6 i! L0 x+ a, {# n+ t
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    4 [1 Y' N2 E, u
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    # J4 S8 L0 T; Y- d
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    . [. ^0 }+ F$ `* ]3 [\" M+ e4 ^
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    6 i/ d' Z; V/ B6 B9 P
40. ,PlotRange All
      s8 d0 `4 }3 G. B! g0 Q$ ?
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
7 m$ m3 U1 r8 q0 w; i" J4 ]8 D7 l8 v

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
7 \6 u  _0 I$ U. o  P/ b2 {