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

小洋洋

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)"}]
! e2 x/ P! T. Z" n
: u& }; e5 A2 T

平凡之不凡

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

小洋洋

# u; e3 z0 e2 o4 _3 n

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,2 n3 Y+ J% e6 [+ v
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   1 L* e3 K! \& |& Q\" R
3. Current={Ii,Io};
   & e& i9 }* a. e9 x
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];8 O5 U$ U5 J  v5 x( ?\" Z
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   ! x3 f; B- f/ J+ Z3 J. W$ y! W# f
6. Print["H(S)=",H]6 N1 d% t, q\" w# I* s% [9 S
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   ! Y, D9 _- l* Y- K
8.    den=Denominator[H];
   : u% w0 [$ @7 Z! r: f, |! V\" j. f
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   % m3 S% }; K\" T. B& }3 v
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));8 v! E* x$ \2 c4 m' B# Y
    
11.    Simplify[{fp,Qp}]];
    4 u# c* l$ k  b3 q
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},: \  g' s9 ^8 ~  ]
    
13.    num=Numerator[H];+ C% s* h4 W, r& }! d
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);6 |3 i; L/ n1 ^$ f8 e2 S
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    6 H1 J! S; a( k5 X1 |
16.    Simplify[{fz,Qz0 fz}]];& s* |- m; |0 j' x. ?
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    % M4 d$ X\" C1 R  v  g
18. Klp=H/.s 0;
    ( ~; I( k2 K, w& M
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    ; o5 [+ h' Y\" b2 N\" d& k
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
      F0 X6 H) f% I* G
21.    gm1=C1*Wp/Qp;1 b9 f0 [- R; ?
    
22.    gm2=C2*Wp*Qp;5 p2 T3 L1 T1 Z, E3 O, b) }) x0 f
    
23.    {C1,C2,gm1,gm2}];! ]# t\" g' a, n! v
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];' T' r! L! }5 h3 ^3 T  k, z6 @- v
    
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    : H4 g% d5 ?' F3 i, Q1 B5 l
26. Simplify[H]
    ; J; p4 i5 A0 ]- D( p6 X
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    ) w* [' f8 E: V* p3 O* u3 R5 G
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    % g* P& e+ c/ v  h0 e7 C2 G# v
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    1 o: Y3 k5 W% I, N7 h8 [) h
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    7 I5 \* V\" T- o1 D& a4 b
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    \" ~7 d1 T\" @4 E& K' G2 \' A
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]2 J! I6 H3 k6 S$ b# b
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    ( F- L0 D, |+ c+ E
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],
    7 f! K+ z* r+ v- `\" _
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    , n! K7 x. k8 F* i5 H+ P: {# l
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    5 C3 G' q- f4 q
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];% U8 x\" m( w* K, H
    
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    - ~/ d: C) J\" |# K
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}
    1 @\" P* l5 O' ]0 b' J
40. ,PlotRange All% o; {, A) j% z\" b+ b. e
    
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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