matlab在动态电路分析中的应用

1、控制系统仿真作业0引言MATLA是“矩阵实验室” （MATRIXLAB-ORATORY的缩写，由美国 MATHWOR公S 司推出的一种以矩阵运算为基础的交互式程序设计语言和科学计算软件，适用于工程应用和教学研究等领域的分析设计与复杂计算。与其它计算机语言相比，它 具有语句简洁，编程效率高，强大而简易的绘图功能，有效方便的矩阵和数组运 算，尤其是扩充能力强。正因为这些特点，MATLAB已成为教学研究与工程应用的不可缺少的助手，自推出后即流行于欧美。MATLAB是基于矩阵运算的，其基本数据结构是矩阵，也就是说：其变量与常量都是矩阵（标量可看作1X N的矩阵，向量可看作是1 x N或NX 1的矩阵）

2、，其元素可以是复数或任意形成的达 式。MATLAB!将复数看作一个整体处理的，即不区分实部和虚部，而且还具有元素群运算能力。由于这些特点，它有利于分析计算电路的各种问题，主要包括： 直流电阻电路分析；正弦稳态分析；动态电路分析和二端口网络等。 分析电路主 要是求解电路各支路的电压、电流等，具体步骤是先建立适当的数学模型， 然后 通过MATLAB件编程求解，并且编程相比其他高级语言更简便；电路越复杂， 效率越高。1基于MATLAB的电路分析1.1直流电阻电路及正弦稳态分析具体分析按以下几个步骤来实现：（1）建立数学模型。根据所给电路建立适当 的数学模型，对直流电阻电路和正弦稳态交流电路， 可以用

3、同一数学模型，因为 MATLAB的基本元素是复数，其数学模型实际上就是电路教材中的网孔电流方程 和节点电压方程的矩阵形式。例如：三变量的网孔电流方程为：Zi 山+Zi 込+L'si iZiiTj-FZj 站+Z屈二因为MATLAB的变量是复数，所以以上的电流和电压变量上万没有加点。其矩阵形式为：Z X I=Us(2) 编程。由MATLA的语句构成的程序文件叫 M文件，它是以“ .M”作为文件 扩展名的文本文件， 可以直接阅读并可由任何文本编辑器建立。调用 M文件输 入电路元件参数并运行程序后即可得到结果。 实际上该程序的编写相当简单，只 要有电路和程序设计的基础知识即可， 且程序不长。

4、为了简化编程，亦可直接利 用MATLAB勺交互命令，输入电路元件参数后也可得到结果。(3) 例题分析。例题 1:如图 1 所示，已知 R1=R2=R3=4Q， R4=2Q， IS =2A, a =0.5， P =4，求 I 1 和 I 2。解：建模。按图1,建立节点电压万程：(1/R1+1/F2)Un1 + (-1/ R2)Ur2 = IS+a I 2(-1/ R2) Um + (1/ R1+1/R2+1/R3) Ur2 = - a I 2+P I 1/R3I 1 = (Urd Ur2 )/ R2I 2 = Ur2/R4整理以上各式并写成矩阵形式：1风+1 遍-1/fr0-iIsL>Jh

5、01證一闵-10JiQ.0阳01，h *0.编程。利用MATLAB®写的程序如下:clear , format compact %俞入元件参数值R1=4; F2=4; R3=4; R4=2;Is =2; a =0.5 ; B =4;%按 ZX X=Y列写电路的矩阵方程，其中：X= U1; U2;11;12z1 仁 1/R1+1/F2; z12= -1/ R2; z13= 0 ; z14= - a ;z21= -1/ R2; z22= 1/ R2+1/R3+1/R4; z23= - B / R3;z24=a ;z3仁 1/ R2; z32= -1/ R2; z33= -1 ; z34=

6、 0 ;z41= 0 ; z42= 1/ R4 ; z43= 0 ; z44= -1 ;%以上设置系数矩阵X的元素Z= z11, z12, Z13, z14 ; z21, z22, z23, z24 ; z31, z32, z33, z34 ;z41, z42, z43, z44;Y= Is ; 0 ; 0 ; 0;X=Z 丫 ; %求解结果I1=X ，12= X ；I 1 , I2 %输出结果运行程序得：I 1= 1 ,12= 2即正确答案为：I 1= 1A, I2= 2A例题 2:图 2 电路中 L1=3.6H,L2=0.06H, M=0.465H, R1=20Q, R2=0.08 Q,

7、RL=42Q, Us=115COS(314t)V,求电流丨1,丨2。图2解：建模。如图2,建立网孔电流方程:(R1+jwL 1) I 1 + jwMI2 = USjwMI1 + ( R2+jwL2+RL)I 2= 0写成矩阵形式： 1 "AUsh0IFfwM(Rr：juf寸利用MATLABS程：clear , format compact%输入元件参数值R1=20; R2=0.08 ； RL=42; L1=3.6 ； L2=0.06 ；M =0. 4650 ； w=10 * pi ； Us=115；%按Z * X=Y列写电路的矩阵方程，其中X= I 1; I 2 z ii=Ri+j

8、* w * L1； z12=j * w * Mz21=j * w * M z 22=R2+j * w * L2+RL；%以上设置系数矩阵X的元素Z=z11 , z12 ; z21 , z22 ； Y=Us ; 0;X=Z 丫 ; %求解结果I 1=X(1) ； I 2=X(2)；I 1，I2，%输出结果r 1=abs ( 11)， angl e仁angl e (I 1 ) * 180/pi ,r2=abs ( I 2)， angl e2=angl e (I 2 ) * 180/pi运行程序得：I 1 =0.0470 - 0.1000 ,I 2 = -0.3501 - 0.0063 I ；r 1

9、 = 0.1105 ，angl e1 = -64.8385 ；r2 = 0.3502，angl e2 = -178.9683上式中r1, r2分别代表电流11, 12的模，angl e1, angle2分别代表电流11, 12的辐角。1.2动态电路分析描述动态电路是用微分方程，一阶动态电路是一阶微分方程来描述的， 可以直接 求解微分方程，但也可以应用三要素法求解。当应用三要素法求解时，电路的全 响应二零输入响应+零状态响应，或者电路的全响应=稳态响应+暂态响应，用公式 表示为：f ( t)= f (s) f (O+)exp (-t / t) + f (s)，该式作为其数学模型。二 阶动态电路用

10、二阶微分方程来描述，要应用拉氏变换来建立S域的数学模型。例3 :如图3所示电路，已知R=2Q, C=0.5F，电容初始电压Uc (0+)= 4 V, 激励为正弦电压US (t)=8cos4t，当t=0时，开关S闭合，求电容电压的全部响应，并绘出波形图。图'解：建模。当t >0时，表征电容电压的微分方程为：dUC /d t+UC / RC=US / RC若用三要素法求解，其解为：uqt)=uqs) + Uc(0+) uqs) exp(- t/ t)式中，Uc(0)为电容的初始电压，Uc(s)为电容的稳 态值，T为时间常数。但由于激励为正弦电压，上述公式适当修改为：Uc(t)=Ucp

11、(t) + Uc(0+)-Ucp(0+) exp(- t / t )式中，Ucp(t)为电容的稳态值，它应是和电源同频率 的正弦电压，只是其幅值和相位不同，记为 UcRt)=Ucpm * cos( t + 9), U cp(0+)为 0+ 时刻的 Ucp(t)值, 即 Ucp(0+)= Ucpm* cos( 9 )。编程。clear allR=2; C= 0.5; T =R * C； Uc0=4; %俞入元件参数Usn=8; w=4; Zc=1/(j * w * C)；t =0:0.1:10 ；Us=Usm * cos( w * t) ; %俞入激励信号UcpFUs * Zc/( R+Zc)

12、; %计算稳态分量Ucp0=Ucp(1) ； %计算稳态分量的初始值Uct= Uc0-Ucp0 * exp(- t /T) ； %t算暂态分量Uc=Uct+Ucp； %计算电路的全响应plot (t ,Uc,'-g',t,Uct,'+r',t,Ucp,'* b' )，grid%绘制稳态分量，暂态分量，全响应的波形图(如图4)legend(' Uc,' Uct',' Ucp)% 波形注释运行程序，得到仿真结果：Uc：电路的全响应Uct :电路的暂态分量Ucp :电路的稳态分量2小结电路分析的基本方法是建立数学模型（一般是方程或者已知电路方程组），并求解方程组，得到各支路电压和电流。当电路规模较大时，求解很复杂，借助 计算机可以大大简化计算量，以前有F0RTRA和BASIC语