《MATLAB及其在理工课程中的应用指南》课件第8章.pptx

第8章MATLAB在电工和电子

线路中的应用举例

8.1电工原理

8.2晶体管放大电路

8.3电力电子和电机

8.4高频电路

BACK

8.1电工原理

【例8-1-1】梯形直流电路如图8-1所示，问(1)如Us=10V,求Ub,I₇,Ude;

(2)如Ude=4V,求Ubc,I₇,Us。

图8-1例8-1-1的电路图

解：

◆建模

此电路中设节点电压为变量，共有Ua、Ub、Uc、Ua等四个。各支路电流均用这些电压来表示

对b点和c点列出节点电流方程：I₁=I₃+I₇,I₃=I₄+I₆,将上述各

值代入化简，得

(8-1)

(8-2)

由,即可解出问题(1)。

对于问题(2),可以推出类似的矩阵表达式，只是输入输出量不同，请读者自行推导，并与下面的MATLAB程序对照。由于同一系统的系数矩阵有许多相同的项，编程时可以利用这个特点来简化其赋值过程。

读者可从方程(1)、(2)中找到a₁1、a₁2、a₂1、a₂2及a₁3、a₂3的表达

式，并用矩阵左除解出Ub、Uc。

当问题(1)给出Us时，这两个方程中只有两个未知数U₆和U,可写

成矩阵方程

◆MATLAB程序

r1=2;r2=4;r3=4;r4=4;r5=2;r6=12;r7=12;

%为元件赋值

al1=1/(r1+r2)+1/r3+1/r7;a12=-1/r3;a13=1/(r1+r2)

%对各系数赋值

a21=1/r3;a22=-1/r3-1/(r4+r5)-1/r6;a23=0

ubc=u(1)-u(2),ude=u(2)*r5/(r4+r5),i=u(1)/r7

%解问题(1)

%输入问题(1)的已知条件

%列出系数矩阵A1

%u=[ub;uc]

us=input(us=);

A1=[a11,a12;a21,a22]u=A1\[a13*us;0]

ude=input(ude=),%输入问题(2)的已知条件

A2=[A1,-[a13;a23];0,r5/(r4+r5),0]

%列出系数矩阵A2

u=A2\[0;0;ude],%u=[ub;uc;us]

ubc=u(1)-u(2),i=u(1)/r7,

us=(r1+r2)*(u(1)*al1+u(2)*a12)%解问题(2)

◆程序运行结果

(1)输入us=10时，得到ubc=2.2222,ude=0.7407,i7=0.3704;

(2)输入ude=4时，得到ubc=12.0000,i7=2.0000,us=54.0000。

【例8-1-2】如8-2图所示电路，在t0时，开关S位于“1”,电路已处于稳态，t=0时，开关S闭合到“2”,求Uc和IR₂的响应，并画出它们的波形。

IR3

R₃

6QL1F

口

18V

Us

2

Is

3A

1S西安电子科技大学出版社

Ic

十

二Uc

图8-2例8-1-2的电路图

IR2

]R₂

12Ω

十

解：

◆建模

这是一个分析暂态过程的问题，先要找到其初值和终值。

在t=0_时Uc(0_)=-12V,Ic(0_)=0,由于R₂是电流源Is唯一的外部通路，故IR₂(0_)=I,=3;当t=0+时，因为电容器端电压不可能突变，仍有Uc(0+)=Uc(0_)=-12V,电流源向两个电阻和一个电容的并联系统供电，两个电阻的电流应等于电容电压除以电阻，即

电容的充电电流为电流源总电流减去电阻电流，故

Ic(O+)=Is-IR₂(0+)-IR₃(O+)=3-(-1)-(-2)=6A

再分析终值，达到稳态后，电容中将无电流，电流源的全部电流将在两个电阻之间分配，其端电压应相同，它也就是电容上的终电压，结果应为Uc=12V,IR₂f=1A。初值和终值之间的过渡波形u(t)按三要素法计算。所谓三要素是指初值u₀、终值u和该段的充放电时间常数T,表示式为

u(t)=uf+(u₀-uf)e-(t-to)T

MATLAB语句可写成u=uf+(u0-uf)*exp(-(t-t0)/T)。在求多段波

形时，每段的初始时刻t,可能不同，需要注意。

◆MATLAB程序

r1=3;us=18;is=3;r2=12;r3=6;C=1;%给出原始数据

uc0=-12;ir20=uc0/r2;ir30=ucO/r3;%算出初值

icO=is-ir20-ir30;

ir2f=is*r3/(r2+r3);