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

小洋洋

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)"}]

. O; C6 ?: W- X
8 B" U* M! e3 I( b& g- u5 \

平凡之不凡

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

小洋洋

4 N% V1 S- y0 n: o5 x2 \

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,' x& x+ L6 C4 ^1 m
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   8 [0 Q# T3 T( V  g$ K/ i
3. Current={Ii,Io};. O& _) A) a0 W
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   ! m. u5 y  u6 v3 J\" X! l
5. H=Io/Ii/.CircuitRespones//Together//Simplify;2 E; }0 c4 o7 W1 H! f' ^8 W
   
6. Print["H(S)=",H]
   3 p( X. @\" Y8 _
7. PoleQpole[H,s]:=Module[{den,fp,Qp},4 G- w+ K  L1 n+ v9 P
   
8.    den=Denominator[H];0 }5 G- O\" \7 U2 G/ y. `( B
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);* f5 x; S* z1 _
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    7 U' d4 N( S6 z1 ~
11.    Simplify[{fp,Qp}]];7 h8 }  D' {1 [5 ?9 r
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    + S9 U9 A, d' ?) x( s. I
13.    num=Numerator[H];
    ! D$ g* |4 G) n9 t6 @/ c
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    6 u% V0 j6 I/ h\" U\" z0 }
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);) r8 q/ U2 B; b9 e/ s, t
    
16.    Simplify[{fz,Qz0 fz}]];
    $ f' a5 A0 \4 n* w
17. {fp,Qp}=Simplify[PoleQpole[H,s]];, s. w& Y+ e( S& U: J
    
18. Klp=H/.s 0;
    7 @$ b1 ?4 r1 A! `& i
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];\" [1 k9 l0 A4 V8 ^: \
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    * z+ I, S6 [7 e) `! V
21.    gm1=C1*Wp/Qp;! r1 O5 B8 k; s; B6 Q# _1 p$ k
    
22.    gm2=C2*Wp*Qp;( M8 c$ ?6 y  A  ~
    
23.    {C1,C2,gm1,gm2}];
    4 V\" s: \4 B' c3 i
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    $ |9 F\" ]/ _% U, O- ?
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    \" _2 B$ L\" s8 z- d! x4 F# W; T
26. Simplify[H]6 D: f. H7 |' O\" R/ h! J, k5 [
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    % A9 s\" |: d$ K9 Y2 ~  ?
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];2 Z1 _$ b! b; l2 \' f+ [
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    : |# z! P8 U7 `, ^! W: Z
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]. [( K: g, N  a0 I, _
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],/ X# s' Z  f+ o4 `0 z
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    ( Z' h( R: p$ |, A
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    & y+ \: c% [\" F) r& t$ @
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],6 Q3 Z# q) @1 \! l0 R
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];5 a* x1 Q3 h3 U5 Y7 O\" H. L
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    & Z3 T& T: p, x
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];; @5 P9 x' P/ p7 H+ ]
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    ( W/ U5 I+ y7 g' e
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    * h5 A, V7 H% |5 K; s
40. ,PlotRange All
    8 r8 S$ ^0 [9 X& v, k
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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