|
发表于 2016-7-7 18:02:56
|
显示全部楼层
MATLAB脚本:
clear;
clc;
syms Vi Vo V1 V2 V3 R1 R2 C s As;
%[Vout V11 V22 V33]=solve('(Vo-Vi)/R2+(Vo-V3)*s*7*C+(Vo-V1)*s*C=0','(V1-Vo)*s*C+(V1-V2)/3/R1=0','V2/R1+(V2-V1)/3/R1+(V2-V3)/3/R1=0','V3=-V1*As','Vo','V1','V2','V3');
[Vout V11 V22 V33]=solve('(Vo-Vi)/R2+(Vo-V3)*s*7*C+(Vo-V1)*s*C=0','(V1-Vo)*s*C+(V1-V2)/3/R1=0','V2/R1+(V2-V1)/3/R1+(V2-V3)/3/R1=0','V1=0','Vo','V1','V2','V3');
Hf=simplify(Vout/Vi)
clear Vi Vo V1 V2 V3 R1 R2 C s As Vout V11 V22 V33;
默认不考虑放大器传函,解出:1/(105*R1*R2*C^2*s^2 + 8*R2*C*s + 1)
考虑放大器传输函数为As,注释第五行,取消注释第四行,解出:(As + 15*C*R1*s + 4)/(As + 15*C*R1*s + 32*C*R2*s + 105*C^2*R1*R2*s^2 + 8*As*C*R2*s + 105*As*C^2*R1*R2*s^2 + 4)
自己可以假设As分析零极点。 |
|