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

小洋洋

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

  N, S2 }1 R6 v& F; f

平凡之不凡

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

小洋洋

7 b2 W; Z$ z' e

1. CircuitEquations={Ii (1+gm1/s/c1+gm1*gm2/s/c1/s/c2)Ig,
   ; Q, |, f7 u1 {4 s\" Z+ r
2.    Io (gm1*gm2/s/c1/s/c2)*Ig};
   & _! Y& y! A, W9 a0 m
3. Current={Ii,Io};
   ) E8 u( q, }* c# V. w1 L
4. CircuitRespones=Together[Flatten[Solve[CircuitEquations,Current]]];/ _\" u; u6 A# S0 H! ^  u* v& {
   
5. H=Io/Ii/.CircuitRespones//Together//Simplify;
   / a5 B$ p9 N0 V, |* m. x
6. Print["H(S)=",H]/ ]8 B/ A0 v9 J# {8 z& X# T+ f
   
7. PoleQpole[H,s]:=Module[{den,fp,Qp},) u7 v4 [, }\" z\" W% e# H6 K0 }/ ]
   
8.    den=Denominator[H];: }% R/ ]& F0 w' c0 X3 F* L6 x$ ~% m
   
9.    fp=Sqrt[Coefficient[den,s,0]/Coefficient[den,s,2]]/(2 pi);
   \" v\" O7 I5 Y% X$ K  ~- Q
10.    Qp=1/(Coefficient[den,s,1]/Coefficient[den,s,2]/(2pi fp));
    . W2 r  D8 i9 X3 e5 [$ H6 D; x6 p# A
11.    Simplify[{fp,Qp}]];
    + V) Z, I0 Y' b& N. l1 @: k/ s- n
12. ZeroQzero[H,s]:=Module[{fz,num,Qz0},2 x/ Y) U4 i) {$ {' ?
    
13.    num=Numerator[H];0 }  Q% ^/ N& j4 {6 u+ f
    
14.    Qz0=(Coefficient[num,s,2]/Coefficient[num,s,1]);
    5 N; q2 Q2 s$ P# |2 v; A& N
15.    fz=Sqrt[Coefficient[num,s,0]/Coefficient[num,s,2]]/(2 pi);
    $ O! D# ]* H- e, J! l& v: C
16.    Simplify[{fz,Qz0 fz}]];$ V  Y$ s1 @% d# A  J6 V5 b' P% w6 r
    
17. {fp,Qp}=Simplify[PoleQpole[H,s]];# t0 M, L) a/ T. M4 d\" m6 I
    
18. Klp=H/.s 0;2 x1 G1 o8 A\" v* H
    
19. PrintLabeledList[{fp,Qp},{"fp","Qp"}];
    1 \$ v, o! h% n* o  W
20. DesignOTA1[Qp_,Wp_,C1_,C2_]:=Module[{gm1,gm2},; U- T! e. M/ i: O( ^
    
21.    gm1=C1*Wp/Qp;
    ; Z/ N\" W6 E, h# b$ M
22.    gm2=C2*Wp*Qp;
    0 o: {! l' M# i( O& W# I
23.    {C1,C2,gm1,gm2}];
    0 V& Z! l8 S7 X; w
24. {C1,C2,gm1,gm2}=Together[DesignOTA1[Q,W,c1,c2]];
    7 G  ~) D4 |9 }* E* D& b
25. PrintLabeledList[{C1,C2,gm1,gm2},{"C1","C2","gm1","gm2"}];& @- x9 D% Q: W* T8 L; L) j: y0 s* {
    
26. Simplify[H]1 \+ K1 C) c, S* U1 G\" a/ ^
    
27. h1=Together[H/.{Q 2.,W N[2*pi*10^6],
    : k, E$ C  p! H$ s
28.       c1 10.*10^(-12),c2 10*10^(-12)}//N];& |- N0 t/ M\" d! N4 c
    
29. Print["gm1=",gm1/.{Q 2.,W N[2*pi*10^6],
    1 j$ x/ f! H3 I* y. t4 r
30.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    , V1 |6 {, O\" t, o% K; p0 i0 T
31. Print["gm2=",gm2/.{Q 2.,W N[2*pi*10^6],! z& X8 ?6 V2 A9 ~7 r
    
32.     c1 10.*10^(-12),c2 10*10^(-12)}//N,"mS"]
    ; d: a( J4 m0 C
33. h=(Numerator[h1]/g)/(Simplify[Denominator[h1]/g])
    / |4 v. y0 q; m! R( e) j3 T
34. PrintLabeledList[N[{Q,W (2*pi)}/.{Q 2.,W N[2*pi*10^6],+ f\" C  K8 m0 [' g1 L1 Q
    
35.       c1 10.*10^(-12),c2 10*10^(-12)}//N],{"Qp","fp(Hz)"}];
    5 J8 U4 Z$ f- B# t0 e8 J
36. Rexample=N[{c1*10^(12),c2*10^(12),gm1,gm2}/.{Q 2.,W N[2*pi*10^6],  e+ V$ z6 k+ Z- f; b
    
37.       c1 10.*10^(-12),c2 10*10^(-12)}//N];
    ! B; `' q  _/ K
38. PrintLabeledList[Rexample,{"C1(pF)","C2(pF)","gm1","gm2"}];8 z1 `3 Y\" \, X
    
39. Plot[{Abs[h]/.s I*2*pi*f*10^6},{f,0.01,8}. z& [8 w, L! c
    
40. ,PlotRange All
    ' Q, o$ a7 \( G5 J
41. ,AxesLabel {"f(MHz)","M(f)"}]

复制代码

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

山林隐逸

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

guozihu1986

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

guozihu1986

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

guozihu1986

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

cnmover

看着眼花啊

guozihu1986

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

小洋洋

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