第6章 MATLAB在《电路》课程中的应用-4.ppt

本章主要内容： 1．采用M文件编程进行解题。 2．掌握Simpowersystem工具箱常用模块的使用方法。采用simulink进行仿真解题。 例 6.1 电阻电路的计算 如果Us＝10V，求Ubc，I7，Ude 如果Ude＝4V，求Ubc，I7，Us 例 6.1 电阻电路的计算 1．采用M文件编程进行解题。 （1）用基尔霍夫电流定理列方程组 把方程组写成矩阵形式为 Matlab 程序 （2）用基尔霍夫电压定理列方程组 例 6.1 电阻电路的计算 例 6.1 电阻电路的计算 2．掌握Simpowersystem工具箱常用模块的使用方法。采用simulink进行仿真解题。 (1) 打开Simulink工具箱 (2) 新建窗口 (3) 打开simpowersystems工具箱 (4) 把电源模块、电阻模块、Current Measurement 、Voltage Measurement模块 、Sink中的scope模块拖到新建窗口 (5) 修改参数，电阻赋值，电感为inf 、电容为0 。 (6) 连线，注意Current Measurement 、Voltage Measurement模块的连接，其输出接scope 。 (7) 仿真参数设置：直流电源simulation－configuration parameters－ode15s ，交流电源MAX step size取0.001 (8) Scope参数设置：双击scope，在scope窗口中的parameter菜单进行端口设置。右击进行坐标设置。 (9) 左键在scope窗口拖动可放大观看。 例6.2 含受控源的电阻电路 已知is＝1A，R1=100，R2=10，R3=20，R4=20。求i1, i2, Ua, Ub 例6.2 含受控源的电阻电路 图6.2 例6.2的电路图 续例6.2 含受控源的电阻电路 例6.3 戴维南定理 已知is1＝1A，is2＝2A，R1=100，R2=10，R3=20，R4=20。求U1，U2， Ua 例6.3 戴维南定理（续） 写成 例6.4 一阶动态电路 已知us＝12V，is＝0.1A，R1=12，R2=12，R3＝12，C＝0.05uf，T时开关S由1位置打到2位置，求电容C上的电压uc。 例6.4 一阶动态电路（续） 第一段： 电压 电流 初始值 uc(0+)=12V 稳定值 时常数 uc(t )=uc(?)+[uc(0+)-uc(?)]e-t/?1 t≥0 iR(t)=iR(?)+[iR(0+)-iR(?)]e-t/?1 t≥0 例5.4 一阶动态电路（续） 例6.5 变压器电路 例6.6 正弦激励的一阶电路 例5.6 过阻尼零输入响应 方法1 uc的微分方程为 写成 初值为: 例5.6 过阻尼零输入(续) 即有?＞ωn的过阻尼情况。其解为 其中 例5.6 过阻尼零输入(续) 方法2 对方程作L变换，考虑初始条件，可得 整理后得 分解部分分式 求反变换 例5.6 过阻尼零输入(续) p1，p2，r1和r2可用MATLAB中的residue函数求出，其格式为： ? [r，p，k]=residue(num，den) 其中num，den分别为分子、分母多项式系数组成的数组。进而写出： u=r(1)*exp(p(1)*t)+ r(2)*exp(p(2)*t)+… 这样就无需求出其显式，程序特别简明。 例5.7 欠阻尼零输入响应 微分方程同例5.6，不再重复。这里 ，当R=1Ω，2Ω，3Ω，…，10Ω时，? =1，2，3，…，10。显然? =?n=10为临界阻尼，其余为欠阻尼（衰减振荡）情况。 例5.7 欠阻尼零输入(续 方程的解析解为 uc(t)=Ae-?tsin(?t+φ) iL(t)=-t?nCAe-?tsin(?t-?) 其中 例5.7 欠阻尼零输入(续) 方法1：把解析解用MATLAB计算，若不要求解析解，不推荐这种方法，太繁； 方法2：用极点留数方法，其程序与过阻尼的情况相同，只不过出现了复数极点和留数。其核心语句就是两条： % 求极点留数 [r，p，k]=residue(num，den)； % 求时域函数 ucn=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t)； 例5.8 简单正弦稳态电路 方程组为： 设Z1=j?L，Z2=R，Z3=1/j?C，R与C并联 后的阻抗为 ，总阻抗为 Z=Z1+Z23。 利用MATLAB复数运算优势编程 例5.9 正弦稳态：戴维南定理 如图5.9所示电路，已知C1=0.5F，R2=R3=2Ω，L4=1H； Us(t