12章MATLAB在电路中应用课件.ppt

第12章 利用MATLAB计算电路;12.1 MATLAB的概述 ; MATLAB的主要结构和功能;12.2 MATLAB程序设计基础; MATLAB 7.5桌面集成环境;（2）命令窗口;命令窗口环境;（3）工作空间管理窗口;（4）命令历史窗口;（5）当前目录窗口 ;（6）Start菜单 ;（7）编译窗口(MATLAB文本编辑窗口);启动MATLAB文本编辑器的3种方法;利用命令窗口和M文件进行编程。;2.程序控制结构;利用输入函数进行数据的输入。;2.程序控制结构;2.程序控制结构;2.程序控制结构;2.程序控制结构;利用多分支if语句进行编程。;3）switch语句 ;利用switch语句进行数据的选择。;2.程序控制结构;2.程序控制结构;2.程序控制结构;2.程序控制结构;2.程序控制结构;12.3 电路的传递函数及频率特性;频域响应曲线;12.4 MATLAB 在基本电路中的应用例举; 解：这是涉及简单直流电路-戴维南定理的应用电路分析与计算。该电路不含受控源，仅含独立直流电流源。戴维南定理直接应用于简单直流电路的简单电路，可以使用网孔法、节点电压法等方法解决。在这里介绍使用节点电压法解此题。;先求ao端以左的戴维南等效电路。 首先用一个电流为ia的独立直流电流源代替电阻RL，以o点作为参考点，设其电位为零，其次设各个节点电压分别为u1,u2,ua，如图12-13所示;最后根据图12-13可以写出该电路的节点方程组如下所示：;将其写成矩阵形式如式12-4所示 ;编程思想为：令ia=0,is=2A,is2=0.5A，由式（12-4）解得u11，u21,ua1。因为ia=0，由（12-5）可得：uoc=ua1。再令is1=is2=0，ia=1A，仍由式（12-4）解得另一组u21，u22，ua2。由于内部独立电源is1=is2=0，故uoc=0。 ;根据对题目的理论分析，编写MATLAB程序代码如下： clear,format compact R1=4;R2=2;R3=4;R4=8; %设置元件参数 is1=2;is2=0.5; %按A*X=B*is列出此电路的矩阵方程。其中X=[u1;u2;ua],is=[is1;is2;ia] a11=1/R1+1/R4;a12=-1/R1;a13=-1/R4; %设置系数矩阵A a21=-1/R1;a22=1/R1+1/R2+1/R3;a23=-1/R3; a31=-1/R4;a32=-1/R3;a33=1/R3+1/R4; A=[a11,a12,a13;a21,a22,a23;a31,a32,a33];B=[1,1,0;0,0,0;0,-1,1]; %设置系数矩阵B;%方法一：令ia=0,求Uoc=X1(3);再令is1=is2=0,设ia=1,求Req=ua/ia=X2(3) X1=A\B*[is1;is2;0];Uoc=X1(3) X2=A\B*[0;0;1];Req=X2(3) RL=Req;P=Uoc^2*RL/(Req+RL)^2 %求其最大功率 RL=0:0.01:10;P=(RL*Uoc./(Req+RL)).*Uoc./(Req+RL); %设RL序列求其功率 figure(1),plot(RL,P), grid %画出功耗随RL变化的曲线 计算结果为： Uoc = 5.0000 Req = 5.0000 P = 1.2500;功率随负载的变化曲线如下图所示 ;2.正弦稳态电路分析 ; 解：采用网孔电流法求解，设如图12-13所示3个网孔电流分别为，下面列写网孔的回路方程:;% 输入初始参数R1=2;R2=2;R3=2;R4=2;C1=5e-4;C2=8e-4;L1=1e-3;L2=2e-3;w=1000;US=10;ZR1=R1;ZR2=R2;ZR3=R3;ZR4=R4;ZC1=1\(j*w*C1) ;ZC2=1\(j*w*C2);ZL1=j*w*L1;ZL2=j*w*L2; % 计算方程组系数矩阵中各元素的值 Z11=ZR1+ZR4+ZC2+ZL1;Z12=-ZR4;Z21=-ZR4;Z13=-ZC2;Z31=-ZC2;Z22=ZR2+ZR4+ZC1; Z23=-ZR2;Z33=ZR2+ZR3+ZC2+ZL2;Z32=-ZR2;;% 组成方程组、求解结点电流 A=[Z11 Z12 Z13;Z21 Z22 Z23;Z31 Z32 Z33]; B=[0;US;0]; I=A\B。 程序运行结果为: I1=1.4451-1.3276i I2=3.5402-1.7874i I3=1.0821-1.3621i;3.暂态电路分析 ;解：首先设电容电压u为未知量，根据换路定理得出:u(0+)=u(0-)=9v，再根据回路方程和电容