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

小洋洋

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


8 m6 d+ G5 J# `5 |8 b

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   \" d. A# I0 e) r0 ^5 I$ _
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};8 }  {8 t2 H  I/ }6 b; e  F
   
3. Current={Ii,Io};8 M: p: A0 Q5 x) W- j' Z
   
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];$ C/ L$ N/ D4 f: i/ k, ~+ a. o
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   \" p- s% f  g% i) l5 `6 \7 B( j
6. Print["H(S)=",H]/ p4 q# l$ h/ P/ d8 f( [
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},! E5 y/ L4 w% j  T- t$ |5 U( a, y
   
8.    den=Denominator[H];8 Z- H! P7 ]- h7 P1 J
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);2 y/ L5 f' u' D0 m  `\" x
   
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    % k* ^) V  h% }
11.    Simplify[{fp,Qp}]];
    ( S, G4 J\" M! ]1 J$ |
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},
    \" D\" B# j( G; ^  _2 ?
13.    num=Numerator[H];; W$ B6 l% X, R1 |; l1 ^6 V/ J
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    + T* l, @$ [* R5 E% Z) B7 x. I2 U5 Y$ k
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
      Q\" a6 k% Q+ r
16.    Simplify[{fz,Qz0 fz}]];9 Y$ p% |% u. c\" ^
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    . f4 ?  f& `4 y: f% m5 x
18. Klp=H/.s 0;
    , q\" [5 b: |\" C9 X5 [
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];0 M8 u- l) Y% {0 X. }
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},* C7 _8 H+ w4 e* K! ~% Y
    
21.    gm1=C1*Wp/Qp;+ M8 I7 [0 M. Y8 L& c5 O- \) D
    
22.    gm2=C2*Wp*Qp;! y# q! [- F/ `+ v6 W5 V' C
    
23.    {C1,C2,gm1,gm2}];
    : Z2 A: E0 V* B2 w
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];  h2 T# @9 @( b' Z! H( d! W: B3 R
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    3 G( {/ y* ]3 K5 G5 C/ {2 o: d
26. Simplify[H]
    7 y\" r$ i\" _9 c, I- B
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    # J3 e: Z! F, n4 |; r, B, _
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    ; r4 C  d) T2 P( z* T
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],. H+ \# t% [: A( |+ ?
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]\" J, T3 z\" S0 R) [+ r
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],, j1 z0 {+ h) C3 m- @
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]$ e\" N( B7 O9 i% r6 h
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    / u+ y2 ]. l/ x0 S1 L/ d: w( A
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],2 K\" d9 u) Z, i
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];- N. s9 ~: t- K: S' T9 `
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],9 F1 u) C. E# E+ b$ x5 v8 |
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];7 H( e  X: B( ~, s6 L
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    ' c0 k6 s1 X6 d+ m! k: D+ v% d
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    , n9 F8 O- `0 v0 `9 A
40. ,PlotRange All# I% r& }! i/ g% q\" e
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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