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

小洋洋

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 T8 L* s) Q  N( V; Z. w: w# Y

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,6 {! j. M% r) ?' f' d& _3 ]
   
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};$ L; z. r; q& j
   
3. Current={Ii,Io};
   4 s, }% k4 {/ P* v
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];- I) v/ B) P7 P  e/ N' O8 [5 [/ K- B
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   ; N\" Y& w8 W! v
6. Print["H(S)=",H]
   - n0 Y/ w, b, o3 B: c+ X0 p9 J
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   ( j; W9 P/ s1 l# v1 v- ?
8.    den=Denominator[H];
   6 j0 X8 X1 ^$ Q/ [+ k$ ~/ w2 E
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   5 w6 N1 n! M$ v/ G
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));5 Z+ r: K! O% y5 e
    
11.    Simplify[{fp,Qp}]];  ?% \: A9 W+ s6 d: g
    
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},$ c* d\" C8 n6 n& K* J\" ~! j
    
13.    num=Numerator[H];+ Z) k1 I8 ^% Z. E
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);1 k) B* X* S% w  x2 t/ P5 h1 _
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    6 J) S; ?! `1 l1 \* y8 U
16.    Simplify[{fz,Qz0 fz}]];% D5 v9 [* _\" l& b! R
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    ) T\" j; J3 i$ c4 K: Y$ I# R1 w
18. Klp=H/.s 0;
    5 u. y$ _+ |8 H  Q7 L3 l
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];+ K1 o) v$ A$ U/ G* X+ @4 L
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},  [9 v6 f9 B: n. ~* b! @: c
    
21.    gm1=C1*Wp/Qp;
    + I; k) P& k: H
22.    gm2=C2*Wp*Qp;0 ^. _\" E: w7 ~+ k! q$ g5 E$ P5 f9 J  Z
    
23.    {C1,C2,gm1,gm2}];2 W4 x# [6 \+ ^; Q+ B; Y\" y' K4 T9 ^
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    8 q0 Y' y# ^$ |/ K\" w
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];$ v) Y+ p! f9 S* L, C- G
    
26. Simplify[H]
    : t# r/ N, P3 M5 {\" p
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    \" m8 |$ L0 h9 {) C( B5 ?
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];. t9 ~' Q1 E2 T' L
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],) J5 X4 f- H5 A' b! E( z
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]6 b* t! Y. R2 h9 U+ w2 L! \
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],
    - \& q: s0 b0 c& {5 k
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    ) f/ ]4 G6 S& ]3 ~$ d9 I5 a6 w
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])6 X, W. J. v9 o& w6 u: n. H: s5 t
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],5 L( ]9 N/ N; J5 [
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    % l' v\" A: l; H9 ]  ?3 S9 D; A
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],) ^0 f$ w  I6 U\" f3 j
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    ' n0 Y8 m8 r. E. _7 D  I
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    - U( ^# B3 h. U3 a
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}+ [( ]4 r' H/ F
    
40. ,PlotRange All
    6 R# o& i: y- E7 |6 \
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
7 z  x* b9 X  L/ {

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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