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

小洋洋

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)"}]
& l( F9 e( p3 c/ l

平凡之不凡

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

小洋洋

1 G+ q- q" x& Z$ H1 N. o

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   + M0 U* [0 M; J! t
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   9 |5 m  A* _- U9 B( T! V5 i5 C1 V
3. Current={Ii,Io};, j3 \0 @' _+ _9 X( O4 \) O\" K
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   $ \/ C6 c, E3 s. r9 i/ |
5. H=Io/Ii/.CircuitRespones//Together//Simplify;+ c+ r# H% }' X# e1 q
   
6. Print["H(S)=",H]$ @# T) e# F7 }# w
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},: V* M; B  {+ C9 X0 W% |# w% l' s
   
8.    den=Denominator[H];$ l4 u/ V  \  ?- F# O1 H
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);2 b\" l! }0 X9 O+ b: j2 P
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    : ]2 N, k/ J+ ]' F2 Y; J( M\" n
11.    Simplify[{fp,Qp}]];6 p; t3 l1 e  p, j- `
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    % R1 N1 E0 }: o+ i
13.    num=Numerator[H];, n: ]$ v) V% v: I8 \+ a9 Y
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);0 z1 L' a) O4 X3 f7 B/ ~
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    5 i3 g: E0 W$ Y9 s+ A& }' j) K
16.    Simplify[{fz,Qz0 fz}]];% n+ P& I) ~2 n$ O$ T6 p: F
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];7 z+ g& ?/ z8 L4 ]( Q6 o, `4 |! k. W) C
    
18. Klp=H/.s 0;
    0 }+ C' i8 H3 X1 z* b) V, U, N
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    / @( ~5 ~. n9 x
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},8 Y2 ]# ~- h7 A& F6 y
    
21.    gm1=C1*Wp/Qp;& A0 L; K3 m+ a0 k
    
22.    gm2=C2*Wp*Qp;
    ( B6 z4 {3 Q3 X& E
23.    {C1,C2,gm1,gm2}];' b. ^' }, a* q3 L
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    - x6 s* @0 y; B, s
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    7 h\" _2 ?6 P9 P% M- U9 J
26. Simplify[H]
    ' B. M+ v& U: W/ N% h
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],/ c6 z: r- e' b3 o5 V- [
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    ) d\" Z( K* i\" ]' D
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    / b; T6 W/ D; H* o- B# N
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    9 V, W$ F3 u1 l7 w7 C: o
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],\" V4 E3 i5 Z6 C0 K. m
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    : F\" g* e  A( ~  q- s) E  P\" s+ f
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    5 Y' a7 e& q& A* n* e, r/ A1 y8 X. X' G
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    7 G4 u2 ?9 x% r0 c; M) ~5 j
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    4 _# [5 s, z  h% G/ G' N
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    # ^, h+ a. L9 T
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];; g/ D3 ]\" N* u+ u2 ^3 n
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
      i( f! N7 F5 w( d. u- R
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    8 e6 `4 P3 u- F* X; ?5 a
40. ,PlotRange All( w! |0 ^/ W\" D( v
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
) B, s( z8 T1 u# u

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
% v8 ]/ @9 p) o6 r( M+ ?