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

小洋洋

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


6 R  o! o2 X4 j( z# L* E3 d9 S6 y

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   ( j& P4 N& P& c; N
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};, R) Y  h9 _9 w8 }) r
   
3. Current={Ii,Io};( J! v1 v9 f; v% Z  z
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   9 k4 n9 a, h0 q$ l. U5 W
5. H=Io/Ii/.CircuitRespones//Together//Simplify;/ ]: W/ a2 y2 w5 E
   
6. Print["H(S)=",H]
     G8 u2 J; @) \( ?
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   \" }# K( i* _8 o- _( h2 s! |
8.    den=Denominator[H];5 ?5 O' ^$ s. K$ `# n) W8 R% x
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);\" O8 S% G6 M1 v: e* \- |6 r
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));# C( F( P8 `. ~* M- p
    
11.    Simplify[{fp,Qp}]];
    4 L; P6 A, ?: q. B4 Q5 N! q8 |# {: R
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    4 d1 K# P8 P+ u8 C* K
13.    num=Numerator[H];2 {6 F; i- ^5 {5 H2 e% d\" I
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);6 U- n5 e4 f/ A\" ?4 J% b# A) C  k
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);! b( }3 O% y5 t& Q5 {0 s
    
16.    Simplify[{fz,Qz0 fz}]];
    # A: t5 ?' f' ~\" K# X2 p
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    / o/ Y3 E  d. a; J/ h+ w
18. Klp=H/.s 0;
    * Y! U! @. T2 F3 {$ Y/ S\" n) J  Z
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    ( T  l: H+ y; y4 r
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    . V! I8 I( V$ x; e. _/ [
21.    gm1=C1*Wp/Qp;- B* o, V6 p$ t) O, u+ V
    
22.    gm2=C2*Wp*Qp;$ a2 L0 z5 q5 l& R, B2 \) I
    
23.    {C1,C2,gm1,gm2}];1 h9 w7 p' q9 E# C
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];4 r% A; i9 ]/ S8 c' r
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];* y) F9 \* P7 ~4 C  E
    
26. Simplify[H]
    . K# r* M) y9 C/ L) q- G
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],) {* t- H5 w% [8 Y: h
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];7 s( `\" @\" P+ }* {/ s; t& j
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],- h: A# |6 B+ m\" u# Z
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]+ d/ s% Z! F' n: y
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    / y4 ^, K- P7 t$ j3 w6 s. r
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]4 Q' p, ~$ U+ O( Z% ]; r, b8 A
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])9 D9 G/ w) v5 Y3 C- x: Z
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
      |& y# e+ i6 T
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    8 [: o9 T  o, |7 z
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],' K2 A3 V. t1 G/ V
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];$ L, W. J- m; x# {' `\" b: H* I4 o0 }
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];$ b/ t3 `: U# D2 O2 l' Q( K3 v9 H
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    4 ?) t; T5 {6 U6 ?+ `4 D
40. ,PlotRange All. B  a. h( _\" ^& v
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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