matlab 中二阶,MATLAB在二阶动态电路中的应用

0引言在电路分析中,对动态电路的分析无论在理论上还是实际应用中都是十分重要的;关于时间的微分方程,手工求解是相当复杂繁琐的,对于二阶动态电路,特别是一般二阶动态电路的求解就更为复杂.然而利用MATLAB来求解二阶动态电路问题,就可以变得简单,便捷.本文以具体实例讲述了MATLAB在二阶动态电路分析中的应用.利用MATLAB的dsolve函数来求解电路的微分方程(组);用diff函数来求导数,可以根据电压求出电流或根据电流求出电压;利用ex-plot函数绘出电流或电压等物理量的曲线.所有这些只需一个或几个语句即可完成,极大地减少了编程步骤,缩短了编程时间.本文中的实例及一般求解过程参见文献[1,2].本文中的求解更加简便、快捷,由于有了波形,电流、电压的变化规律就更加直观、一目了然.在MATLAB7.0程序中,使用dso1ve函数求解微分方程的具体格式参见文献[3,4],MATLAB程序不支持上下标,但必须区分字母的大小写.1典型二阶动态电路1.1零输入响应如图1所示为零输入响应电路[1].开关断开前,电路已达稳态,t=0时将开关打开.求t>0时电容电压响应uC和电感上的电流iL.图1零输入响应电路微分方程为:d2uCdt2+10ddutC+16uC=0,初始条件:uC(0+)=2.5V,duCdt0+=-i(0C+)=-i(0C-)=-4,(此处文献[1]中有误,duCdt0+=4).应用MATLAB求解的M文件为:uC=dsolve(D2y+10*Dy+16*y=0,Dy(0)=-4,y(0)=2.5),运行结果为:uC=-1/6*exp(-8*t)+8/3*exp(-2*t).即uC=38e-2t-61e-8t,t>0.利用diff函数,对uC求导数,可以求出电流iL:iL=-CddutC=-0.125ddutC.M文件为:symst,%使用syms函数定义符号变量tiL=-0.125*diff(-1/6*exp(-8*t)+8/3*exp(-2*t))%即iL=-0.125*diff(uC).运行结果为:iL=-1/6*exp(-8*t)+2/3*exp(-2*t),即iL=32e-2t-61e-8t,t>0.用MATLAB的ezplot函数画出t>0时(0-5s)uC和iL的响应曲线.M文件为:symst%使用syms函数定义符号变量t;ezplot(-1/6*exp(-8*t)+8/3*exp(-2*t),[05]);ezplot(-1/6*exp(-8*t)+2/3*exp(-2*t),[05]),uC和iL的响应曲线如图2所示.图2uC和iL响应曲线此例中R>2CL,所以是过阻尼状态.由于电阻R较大,电容、电感同时释放能量被电阻消耗,电路中的能量很快被耗尽,所以,电流很快衰减到零.由波形可以观察到电流和电压的变化情况.注意,此例与许多文献中所分析的过阻尼状态的区别.1.2零状态响应如图3所示为零状态响应电路[1].(t)为单位阶跃函数.求电容两端的电压响应uC和电流i.图3零状态响应电路微分方程为:0.25d2uCdt2+ddutC+uC=1.初始条件:duCdt0+,uC(0+)=0.应用MATLAB求解的M文件为:symstuC=dsolve(0.25*D2y+Dy+y=1,Dy(0)=0,y(0)=0).运行结果为:uC=-exp(-2*t)-2*exp(-2*t)*t+1,即uC=-(1+2t)e-2t+1,则电容电压响应为uC=[-(1+2t)e-2t+1](t).利用diff函数,对uC求导数,可以求出电流i:i=Cdd