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

小洋洋

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

平凡之不凡

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

小洋洋

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   $ I2 R/ Y' s% d) P7 w
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   6 q6 q- j: ~\" Q  l9 B
3. Current={Ii,Io};
   # m5 c( e\" i) m/ F/ P6 a& s  V' |3 o
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];
   , G; m9 P' Y# C9 w' }( x5 e
5. H=Io/Ii/.CircuitRespones//Together//Simplify;, N# ]# v\" H' u9 f
   
6. Print["H(S)=",H]8 ^$ x: a  {1 n$ e( `\" C
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},
   . ^7 A, p5 p) }. F8 B, p
8.    den=Denominator[H];  C, j: X* A\" a, c
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   . A0 g& n: i- Q\" X2 r: z$ J
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    1 T1 o8 @; K: K9 V5 x% P% d
11.    Simplify[{fp,Qp}]];
    4 \% G1 X3 B/ Q0 g. b! H! A
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},: l. P1 V8 U- n
    
13.    num=Numerator[H];' c% M$ ~  h; R5 S, X$ `
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);- h5 W$ p: N6 M
    
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);$ z$ ?0 O( U$ \
    
16.    Simplify[{fz,Qz0 fz}]];  h/ s2 m* c( N$ u: ]2 n4 T# O
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];
    & s\" W7 L& i' ^( i* U
18. Klp=H/.s 0;
    2 Y) ~3 T; b2 Q
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];) k8 z' \+ H- M& g9 T# v3 T
    
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},
    1 o9 ~' {9 Q& S0 a
21.    gm1=C1*Wp/Qp;3 s8 z/ ^. G6 r! E& R) A
    
22.    gm2=C2*Wp*Qp;
    ) J0 Y- @3 p. X) b7 [
23.    {C1,C2,gm1,gm2}];! e8 `6 \2 ~% n9 H- V. [' x4 [
    
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    2 I% w4 H4 J3 ]: X5 o4 W. N9 l
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];
    . b1 q. E0 E$ K1 ]5 c
26. Simplify[H]
    5 s! I; z: a9 W1 [2 X$ u+ Y: j+ g  L
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    # V- p1 ~0 l, A
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    1 {, B* M# T6 ?% X/ c8 R
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],# X\" I4 N% Q9 x0 m6 p
    
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"], ^% k0 J' B0 ]0 z3 \; t! l& z
    
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],% j/ g% a# Y\" A\" z
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]# B8 [( P) L/ G7 `4 n/ j
    
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])/ ~! b& Z* i5 J, A! P
    
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],# F7 J8 J\" K: o
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];* e2 p! F0 R: \, F6 }  \
    
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],
    ; P4 m+ X# O& X. v
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    0 D2 _; K' B! l
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];
    % |6 {9 h( \5 V
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}- L$ T. c  S$ B2 Y
    
40. ,PlotRange All
    0 `+ p1 ^( S9 u$ ~1 e
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

这是调整完的程序,麻烦大神们帮忙看看吧
4 @3 u& i8 ?6 w3 m% z7 @9 ?9 l

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

我想模拟一个二阶的OTA-C滤波器，也是才学习这个软件，不知道这个程序问题出在哪里了，非常焦急，我把.nb文件穿上去了，望指教
+ m( G2 |2 p3 y: b; `, e