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

小洋洋

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)"}]
$ X  x6 a' C7 c8 \5 i6 C  B) K
# D! C' {6 t; Y# \7 m

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,, b( {0 W! U. |5 x) I* y
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};& {5 I. D8 ^* z  ]* v. E& @$ `
   
3. Current={Ii,Io};
   , K7 n0 s! }* M6 \! L* }* R
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   9 q; w: ]+ c$ \( Z
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   1 q5 E! j6 D, [8 ^  H
6. Print["H(S)=",H]
   2 W1 K\" d* K' O, e) |; X! u' D
7. PoleQpole[H,s]:=Module[{den,fp,Qp},$ T, z' ?$ S9 e8 o
   
8.    den=Denominator[H];
   ' b' W' I, T& `: ?' e; t9 l
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);. h5 _* a+ `  Q5 z2 z/ B( d# V\" N
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));. O0 [9 g; Z\" a8 P6 A9 s, _
    
11.    Simplify[{fp,Qp}]];. \. i6 c) S  ~( A
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},' B, W5 u: p; l2 c
    
13.    num=Numerator[H];; [2 a& `, P9 G  H7 B+ s
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);# D0 I6 K3 ?/ E% P0 H: R, S
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    - {. _! e, j1 |1 f
16.    Simplify[{fz,Qz0 fz}]];
    \" B) b) H$ X( W& l4 U
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    : ^& ^1 \/ k2 {. ~
18. Klp=H/.s 0;! y( l: P6 e& i; V0 L9 u
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    + v# T\" B9 [  R7 u% E- V' j9 ?+ r
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},; z! L8 ~) r2 D( Y4 R0 o; C' _4 t/ A
    
21.    gm1=C1*Wp/Qp;& z2 g9 U4 {* R1 j( w0 o* _: V
    
22.    gm2=C2*Wp*Qp;# A9 j/ Y& q. ~* M# }
    
23.    {C1,C2,gm1,gm2}];
    , p9 Y\" E8 X6 h8 }4 ~( I
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];% _3 ]# Z4 B( U, E: M4 L0 R
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];# i( k. B8 S7 _4 @& T8 j\" \2 \
    
26. Simplify[H]) r7 O; |1 V; c: t# Q
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],, k' Q- |\" i3 @: F& ]; n( N$ R
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    , r9 i% w# I7 ]% C
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],* ?& D$ `) x( N* `
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    \" Q/ ^9 S  h# c
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    0 |. S4 E1 d% d5 d$ W5 e
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]* _* P\" B\" F\" e5 `9 N& ~$ x
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    . n% m\" {& ^/ L; s( t; S- F/ l
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    - N% X+ R% X# |$ ]. q
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    & S2 S% G8 J6 G3 J) v8 ^8 U
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    ( Q% v' h7 _7 w7 W
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];% ^# b+ o, V: y4 \
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];\" x  {9 ^5 V! @5 [+ u
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}, v- V# b0 @) w4 j) B5 j- w\" k
    
40. ,PlotRange All# y# e, g5 g! \' d9 l. e# S: i
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
* A) l+ e, s- V6 q$ J+ @