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

小洋洋

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)"}]
" u$ N. u; d4 v$ A
, h3 j* g; Y: L; k1 o: I0 s9 b

平凡之不凡

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

小洋洋

- x. M9 m5 T0 w0 n% c

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,. _: w% c: n8 ?7 o0 |$ B
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   & h! R$ ?2 X9 m- I
3. Current={Ii,Io};: ?5 O+ W$ X4 R# @# Q. A6 ~3 a
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   ' l: v4 N& w; f. @* R* v
5. H=Io/Ii/.CircuitRespones//Together//Simplify;$ q$ h  ?* g+ D0 Q: y
   
6. Print["H(S)=",H]\" f* q. a( p  f
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},$ ~  W  _3 |' K! U0 v- a; J4 W
   
8.    den=Denominator[H];& d4 B# O( a3 c  g1 s7 E
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   2 G9 y* w/ N4 O, D
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    . n$ @& i* x' r3 Y8 ~! Z+ P9 K
11.    Simplify[{fp,Qp}]];
    \" g7 o4 L' Z/ v. b/ ^# W9 z7 K: B; D
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},( ^) P# u- N6 U. H- j% o  H) i
    
13.    num=Numerator[H];
    9 [$ i; |' C5 {* _- J6 [
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    1 ^9 {/ A' w0 X- [8 V7 \2 ^
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);5 r! j\" ^* F' F, U: C
    
16.    Simplify[{fz,Qz0 fz}]];
    \" r) H7 H2 ^4 r+ o3 T8 S, k
17. {fp,Qp}=Simplify[PoleQpole[H,s]];$ b# y; X2 C. j8 _6 b1 g
    
18. Klp=H/.s 0;
    * i; I; I2 X: |6 [) O' B! T3 _
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    ) F3 p; c4 L; o& W; {9 S+ c
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    + B5 g: Z, B* h8 G' E$ m( Z\" L
21.    gm1=C1*Wp/Qp;) W7 i0 V  C, B* C9 E( f+ b
    
22.    gm2=C2*Wp*Qp;
    - I. h6 ^& t0 {# k  m\" {; l
23.    {C1,C2,gm1,gm2}];% D9 @/ ~! D/ z
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    4 Q* g7 Y0 T1 q4 P- b( q7 P! N
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    5 o+ v% Y; W0 M$ R  v, ~
26. Simplify[H]* O7 Z( c8 l. d8 |4 Y
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],8 i' ?! Q1 U5 J+ d4 l5 E
    
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    0 J4 c# ^: ]( Y5 K  o! b- @
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],\" \1 l! ~9 c8 Z# P0 }0 t
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]7 `2 W# x1 a2 Q5 g- R! L* w
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],$ U5 @& ]' p) \8 \! a
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]4 o- W# O- P4 k( I+ b5 u9 f\" I! U' N. A
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])* I* a4 H' @6 W; e% o9 g( Z8 {
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],# V# b. X- Z3 }1 p1 G
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    ! J! N- n$ b+ l( M
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],. l5 V7 u4 f% k3 U% r3 j7 w
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];4 A5 ?* R- ~' e
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    ! t% q$ w( t2 v$ @( a' p* K8 l  }: J
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}\" n3 D4 w) q, p6 P7 z  R
    
40. ,PlotRange All5 e: D8 ^/ f3 E0 ~( c# a+ f4 I
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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