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

小洋洋

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)"}]
7 C/ _3 p/ L0 e

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,/ i- H4 i; H1 y2 v5 c% P
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};2 b\" I* \$ K+ @. {8 a4 V( t
   
3. Current={Ii,Io};
   ' w+ ^  _) s) I0 n& o! l
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   ' I! y3 T. C' }\" }3 `( w- @$ V
5. H=Io/Ii/.CircuitRespones//Together//Simplify;, C\" C7 w+ y) ~
   
6. Print["H(S)=",H]% s, F) [, v$ h' j: p
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},* b9 |7 ~0 ]4 }, _; K
   
8.    den=Denominator[H];5 n! G9 g. C4 E) ^# v; [7 W7 D
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);\" _! z# k8 f' T4 R7 L( V: ?' b
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));0 f) P' Q. h1 W\" O
    
11.    Simplify[{fp,Qp}]];
    - `* W. O7 j3 U1 h; x0 `
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},; z- F, h/ h# ~2 K3 e/ E
    
13.    num=Numerator[H];3 {# X6 l! m, ^7 G( Y
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);8 R1 p$ ^  @) O5 ]
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);7 P, i% Q6 @' Q) p0 J) o
    
16.    Simplify[{fz,Qz0 fz}]];1 k: C\" v\" i6 L- `% w
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];1 Q1 I\" o1 G  [
    
18. Klp=H/.s 0;
    + K- K  R- ^. Z7 ]: Z8 m/ Q
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];$ d) n5 P2 F  {# R  M' `! Y
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    . O$ h7 d: f' q5 l, Z4 j  s
21.    gm1=C1*Wp/Qp;) o# D$ ^$ h( V9 k7 O1 p
    
22.    gm2=C2*Wp*Qp;
    . B, z# ]  }. q5 q! O3 |/ c$ L; S* C
23.    {C1,C2,gm1,gm2}];% U5 E% K4 J% p! s: o
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    8 f0 z2 P7 C' x- w5 ~1 M
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];) W% G) Q- j. f4 ~+ J
    
26. Simplify[H]
    ' e9 q' b0 n' \* `- u. Q
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    9 E4 K2 N\" V7 ]/ H% X
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    ( R- I. f2 t7 c; n. e8 X& R7 a8 K
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],1 }: _7 Y+ U1 V. D5 r# H0 u
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]) O( o\" r- Z0 R; n
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    ) r( ]8 V$ O( E' x1 J+ c
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]. C- U: r6 W5 R, o7 S+ L
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    $ x4 F5 g6 J8 v; m
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    % v6 T2 f! z' V* Z* J7 K
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    1 u) p) _7 @1 c+ W; P/ t/ J
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    5 A) j/ ?3 m0 V; t% e7 |9 @
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];0 W: X. F; t( O( H  X4 ?
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    0 i0 Y. ]. N3 M
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}. i' {8 X+ o. f7 g+ h4 ~
    
40. ,PlotRange All4 o7 _) \4 H1 P0 `
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
1 u, I3 [# b9 D9 Z+ v8 h

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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