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

小洋洋

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)"}]
' r7 I$ {5 a/ `/ e. {' ~
5 }+ i4 K) b0 I0 x2 S
* O* c+ L6 w$ i- c. C$ n

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   + v, [, u0 m3 X( [
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   ' K. Q7 u. K1 |  Q. p3 ^& V
3. Current={Ii,Io};3 @; P; _* E* y7 N: B! P
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   0 Q8 [( T0 `- u4 h1 d; J
5. H=Io/Ii/.CircuitRespones//Together//Simplify;* [# S1 Q' Y) X& A- F, z) m$ A
   
6. Print["H(S)=",H]
   : I2 Z3 j6 Y7 p: B! U\" e: C; D
7. PoleQpole[H,s]:=Module[{den,fp,Qp},\" R5 G% y' w+ @- z# o
   
8.    den=Denominator[H];7 G. L5 b; n; V
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   ) U* q/ d7 o- m1 V  B3 D$ N( Q
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));' ?* Y& V/ q2 |8 i6 [4 K
    
11.    Simplify[{fp,Qp}]];* y* r  [- k5 c& U1 G2 k$ w! y* Y; I
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},/ B. j; x! t3 |# s9 x' x\" p! j
    
13.    num=Numerator[H];
    % W7 s: O5 O9 O5 p9 ^7 ], a
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    : F% x; S9 c/ Q* s' e
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);- M/ O% i% k9 w  K! f3 @
    
16.    Simplify[{fz,Qz0 fz}]];
    ; J3 A  @! \( q$ V0 e% c5 U
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    : f  Q. b, ^; S$ `
18. Klp=H/.s 0;
    - W, k. J6 ?2 o8 r# k& u  f! I
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    / f: F; h8 k' x8 r) ]+ A
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    4 k. u6 X( }6 _  F4 }\" W$ u) [
21.    gm1=C1*Wp/Qp;
    - Q, e& H; D1 t; k; x$ b( G  M6 ?
22.    gm2=C2*Wp*Qp;( q( |9 B6 _( s* A7 N6 G
    
23.    {C1,C2,gm1,gm2}];% S( q$ a3 t- h  g& ?
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    9 X: F) `5 p+ S* }
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    $ Y! _2 p/ I# B/ m, ?
26. Simplify[H]3 [, \5 J5 p% u3 q- Y
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],9 k2 c( a\" c, H+ D& N8 I% p- x- N
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    : u! ~0 [  m5 J3 l  q
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    ( @8 _) Z2 Z4 I! R* _5 [5 i0 N$ h
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    \" A* x9 X6 E$ z
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    6 j: [% s, M! t' |& e\" ^
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]' r' F( f1 D1 N3 v: [3 C4 d
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])3 K3 V4 l2 J; |$ I
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    ' R\" o, P& V) c+ |! [. d
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    . i/ a( M3 f  ^4 u4 R9 ]
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],( X' B/ G% q0 Q$ y% e$ ^3 P
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];6 @( F$ r) s) V* |+ {: J
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    9 W0 N5 J  s4 U0 I1 l+ A
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}/ r6 F9 A/ t6 g& Z) O% D
    
40. ,PlotRange All' L7 l: g  w+ n( c
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
1 q( O" i5 s) N, A0 C