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

小洋洋

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 A5 W( v' N% v) q

平凡之不凡

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

小洋洋

. o8 U4 X) Q# c) F

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   % \: |8 o& u8 g  O  V; s
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   6 N( [( u; Z, P0 J
3. Current={Ii,Io};/ w0 \' p' ^' q/ q( w5 T. e
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   ' u4 Z# C6 a/ `
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   3 J# i, F1 U# i' x2 G\" ~) G
6. Print["H(S)=",H]
   $ {) I' [, N; g; ~, b; ]6 I
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   ( t: l) r1 y9 T/ c+ ^8 ^
8.    den=Denominator[H];
   2 z9 e0 i: e3 S; v) H: T; K
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);7 `' e3 W, p0 d6 u: o\" H2 `! U: f
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));+ G: V. M; z3 c6 Z  `6 ]
    
11.    Simplify[{fp,Qp}]];
    1 ~, m% g- l# j# u/ U
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},8 y0 ^% F4 `/ S/ M: B9 ~
    
13.    num=Numerator[H];
    $ o* N: S9 A7 Y, J' m$ E9 h- m
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);2 {) k! C& @; r/ G
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);: o4 ]\" S  c' F: n3 f& ~
    
16.    Simplify[{fz,Qz0 fz}]];
    . f  N# s$ W/ U# g* X
17. {fp,Qp}=Simplify[PoleQpole[H,s]];! M# R7 X% U, Z! b' H' }
    
18. Klp=H/.s 0;
    3 V9 N- X/ v7 \) b1 F; e
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    ( f7 }) t6 l6 }% p
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},, K9 b0 @0 H, s
    
21.    gm1=C1*Wp/Qp;
    # [! }( {8 |, \& @8 h' Z
22.    gm2=C2*Wp*Qp;1 m: Y6 E0 [2 B  {\" b
    
23.    {C1,C2,gm1,gm2}];
    / V, J+ G; H  O$ [' Q4 O
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];5 u6 s  Q5 J3 H; s4 u
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    % Y) x$ x' N* f# Z
26. Simplify[H]
    ) ]& r& N/ X6 t; ]' j
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    ) e. w4 b$ \: X8 x; A7 D+ n7 d
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    # h' `  S# o3 b$ [6 [\" x& F7 F
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    + O- c) Z, E, D3 d) a
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]. F  I3 J  B- x
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    3 l5 Z+ R\" e( t- e
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    1 c6 a1 R3 {2 x/ ~) f
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])) t& e2 D9 x! B# H! z\" A
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    , G+ z; @0 a6 ]4 C. C0 a/ w7 L7 e
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    & ~7 M\" B5 Y! E3 s* a* q
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],1 _  T* y0 u* k4 v; @
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];/ v) A6 e$ k) L- U3 T% q% x
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];8 k\" J# k( D& |, x( L
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    + B1 N! d5 _0 d2 m& r
40. ,PlotRange All% U8 W4 I- h; U\" f# K( _8 k
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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