TLAB及其在理工课程中的应用指南 第三版 十二五陈怀琛 第5－9章 第8章.pptVIP

　　解： 　　◆ 建模 　　先用节点电位法列写其方程，设U1、U2、U3, U4如图 8-16 所示，这四个KCL方程如下： (8-3) (8-4) (8-5) (8-6) 整理成矩阵形式(注意其中s为拉普拉斯算子), 则有 写成 a*U=b 故 U=a＼b 　　给定输入US=1和一个频率ω, 从这个方程组就可解得该频率输出的复数解x, 它的四个分量分别为该频率条件下的U1、U2、U3、U4，对于规定的频率数组作循环，即可求得其频率响应，包括振幅特性和相位特性。为了探讨C2和Ce的影响，可在上述程序的外面，再加两个改变C2和Ce的循环，由此编成程序如下。 　　◆ MATLAB程序 　　w=logspace(0, 3); %规定频率范围及数组值(从 100~103, 按等比取 50 点) 　　C1=1e-5; rs=100; rb=1e4; %给元件赋值 　　hie=1000; hfe=100; us=1; 　　re=200; rc=1000; C2=1e-5; rL=2000; 　　for C2=［C2, 10*C2］ %对C2及10*C2分别循环计算 　　Ce=1e-4; for Ce=［Ce, 10*Ce］ %对Ce及10*Ce分别循环计算 for i=1∶length(w) %对各个频率计算各点输出 s=j*w(i); a11=1/rs+s*C1; a12=-s*C1; %给a矩阵元素赋值 a21=-s*C1; a22=s*C1+1/rb+1/hie; a23=-1/hie; a32=(1+hfe)/hie; a33=-((1+hfe)/hie+1/re+s*Ce); a42=hfe/hie; a43=-hfe/hie; a44=-(1/rc+s*C2./(rL*C2*s+1)); a=[a11, a12, 0, 0; a21, a22, a23, 0; 0, a32, a33, 0; 0, a42, a43, a44]; b=［us/rs; 0; 0; 0］; %给b矩阵元素赋值 x=a＼b; u(∶, i)=x; %求与第i个频率对应的四个输出电压 end s1=j*w; uL=u(4, ∶).*rL./(rL+1./(C2*s1)); %求负载电压 loglog(w, abs(uL)), grid, hold on %绘对数幅频特性图 end end hold off 　　◆ 程序运行结果 　　执行此程序所得的结果如图 8-17 所示, 可以看出，在所给的参数下，把C2加大10 倍可把低频区低端的幅频特性提高将近 10倍，把Ce加大10倍可能在低频区的高端提高其幅频特性。 图 8-17 C2及Ce对放大器低频幅频特性的影响 　　【例 8-2-3】 运算放大器电路如图8－18所示，试分析放大器开环增益和频率响应对整个电路闭环频率响应的影响，并绘出曲线。 图 8-18 运算放大器等效电路 　　解： 　　◆ 建模 　　设运算放大器的开环增益为A，它是频率的函数，则在图示的连接方法下，闭环输出与输入电压之比为 (8-7) 　　在增益A很大时，分母上的第二项可以忽略不计，因而得出理想运放的闭环传递函数 (8-8) 　　式中的s为拉普拉斯算子，将它换成jω，就得出频率响应，所以这两个式子都是复数方程。根据题意，要考虑A=A(ω)对H(ω)的影响，计算将十分冗繁，利用MATLAB可以方便快速地解决这个问题，但必须给出具体数据求数值解。 　　通常，运算放大器的开环传递函数中包括三个实极点，即 其中ω1＜ω2＜ω3，取负号后为其三个极点，A0为直流增益。 　　为了避免自激，通常使ω1和ω2差得很大，例如ω1＜500 1/s，ω2＞1061/s，而且ω2和ω3也要拉开。现设ω1=500，ω2=2×106，ω3