matlab编程基础与工程应用第八章课件.pptVIP

MATLAb 编程基础与工程应用 机械工业出版社 王敏杰 朱连轩 潘金凤 第八章 MATLAB仿真 本章介绍MATLAB在电路分析、信号与系统、数字信号处理、通信原理、自动控制原理以及数字图像处理课程中的应用。 8.1 MATLAB电路分析仿真 无源元件全部为电阻的电路称为电阻电路。若电阻与电路中的受控源均是线性的，则该电路称为线性电阻电路。在线性电阻电路中，由于电阻的电压与电流之间具有线性性质，所以求解电流、电压问题都可以归类为线性代数方程（组）的求解问题。用户首先需要根据电路原理列出相应节点或回路的电流、电压方程（组），然后再应用MATLAB的数值计算方法，求解待求的电压、电流值。 8.1.1电阻电路 【例8.1-1】电阻电路如图8-1所示，已知R1=R2=R3=1Ω，R4=R5=R6=2Ω，uS1=4V，uS2=-2V，求I3。 电阻电路可用回路电流法、支路电流法、节点电压法等方法求解，本例给出回路电流法。 【解】 解法一：回路电流法。回路电流法以回路电流为变量，根据基尔霍夫电压定律（KVL），列写电路的独立回路组的KVL方程。如图8-1，将3个网孔作为选取的独立回路组，则回路电流法的KVL方程组为 (R1+ R6+ R2) I1- R6 I3- R2 I2=- uS1 (R2+ R4 + R5) I2- R2 I1- R5 I3=- uS2 (8-1) (R3+ R5+ R6) I3- R6 I1- R5 I2= uS2 将各电阻的数值代入方程组（8-1）得 4 I1 - I2- 2 I3=-4 - I1 +5 I2- 2I3=2 (8-2) - 2 I1- 2 I2 + 5 I3 =-2 clear A=[4, -1, -2; -1,5,-2;-2,-2,5]; b=[-4,2,-2]; I=A\b; disp([电流I3为：,num2str(I(3)),A]) 程序运行结果为： 电流I3为：-1.2941A 8.1.2 一阶电路 一阶电路属于动态电路，也称为一阶动态电路。动态电路指含有储能元件（电容或电感）的电路。当动态电路的结构发生变化时，储能元件将存储或释放能量，而且需要一个过程，因此电路也需要一定的时间才能再次达到稳态。一阶电路是指只含有一个电容元件或只含有一个电感元件的动态电路，其方程为一阶线性常微分方程。 当非零初始状态的一阶电路受到激励时，电路的响应称为全响应。一阶电路的全响应f(t) 由初始值f( 0 _+ ) 、稳态解f(∞) 和时间常数τ三个要素决定，其数学表达式为 【例8.1-3】已知电路如图所示，电感无初始储能，U=6V，R1=5Ω，R2=3Ω，t=0时开关S1合上，t=0.1s时开关S2合上，求两次换路后的电感电流。 【解】根据动态电路的三要素解法，求得t=0开关S1合上后电感电流的初始值、稳态解与时间常数，即可求出从t=0开关S1合上后到t=0.1s时开关S2合上之前的电感电流。 当0t0.1s时， 初始值： clear; clc U=6; R1=5; R2=3; L=1; ts=0.1; te=4; t1=0:0.01:ts; %开关S1合上后到开关S2合上之前的时间段 ichu1=0; tao1=L/(R1+R2); iwq1=U/(R1+R2); i1=iwq1+[ichu1-iwq1]*exp(-t1/tao1); t2=ts:0.1:te; %开关S2闭合之后的时间段 ichu2=iwq1+[ichu1-iwq1]*exp(-ts/tao1); tao2=L/R2; iwq2=U/R2; i2=iwq2+[ichu2-iwq2]*exp(-(t2-ts)/tao2); %开关S2闭合之后的电感电流 figure(1), plot([t1,t2],[i1,i2]), axis([0 4 0 3]) xlabel(\it{t}) ylabel(\it{i}_{t}) grid on 运行结果 8.1.3正弦稳态电路 clear clc w=1000; ZR=20; Us=160; ZL=j*w*0.04 ZC=-j*(1/(w*20*1e-6)); Z1=(ZL*ZC)/(ZL+ZC); Z=ZR+Z1; U=Us*(Z1/Z); I=U/ZL; Iy=abs(I) 程序运行结果为 Iy = 3.9801 8.2 MATLAB 信号与系统