MATLAB在电路中的应用.ppt
《MATLAB在电路中的应用.ppt》由会员分享,可在线阅读,更多相关《MATLAB在电路中的应用.ppt(59页珍藏版)》请在咨信网上搜索。
1、MATLABMATLAB应用(三)应用(三)Matlab Matlab在电路中的应用在电路中的应用 MAMATLTLAB AB 中的变量与常量都是矩阵中的变量与常量都是矩阵(标量可看做标量可看做1111阶的矩阵阶的矩阵,向量可看做向量可看做n1n1或或1n1n阶的矩阵阶的矩阵),其元素可其元素可以是复数和任意形式的表达式以是复数和任意形式的表达式,它具有元素群运算能力。它具有元素群运算能力。MATLAB MATLAB的的这些优于其他语言的特色这些优于其他语言的特色,有利于分析计有利于分析计算电路的各种问题算电路的各种问题,并且使编程更简便并且使编程更简便,运算效运算效率率更高。更高。2 2 2
2、 2 学习目的:学习目的:通过介绍通过介绍计算电路问题的编程方法和技巧计算电路问题的编程方法和技巧,逐步熟悉逐步熟悉 MATLABMATLAB语言的语言的使使用。用。例题的解法本身例题的解法本身,不一定最佳不一定最佳。求解电路的专用软件:求解电路的专用软件:求解电路的专用软件:求解电路的专用软件:SpiceSpiceSpiceSpice、PSpicePSpicePSpicePSpice等软件等软件等软件等软件3 3 3 3内容内容:电阻电路的求解电阻电路的求解 (例例1-3)1-3)动态电路的求解动态电路的求解 (例例4-7)4-7)例题分析过程:例题分析过程:例题说明例题说明 求解过程:求解
3、过程:建模建模 Matlab Matlab程序说明程序说明 Matlab Matlab程序运行、结果演示程序运行、结果演示4 4 4 4l 电阻电路的求解电阻电路的求解(1)(1)如如u us s=10=10V V,求求i i3 3,u u4 4,u u7 7;(2)2)如已知如已知u u4 4=6=6V V,求求u us s,i,i3 3,u u7 7。图图1 1 例例1 1的电阻的电阻电路电路 例例例例1 1 1 1 如如图图1 1所示的电路所示的电路,己知己知:R R1 1=2=2,R R2 2=4=4,R R3 3=12=12,R R4 4=4=4,R R5 5=12=12,R R6
4、6=4=4,R R7 7=2=2。5 5 5 5 对对图示电路图示电路,用网孔电流法,用网孔电流法列写网孔电流方程如下:列写网孔电流方程如下:A)A)建模建模 解:解:解:解:6 6 6 6写成矩阵形式为:写成矩阵形式为:也可直接列写数字方程为:也可直接列写数字方程为:R R1 1=2=2,R R2 2=4=4,R R3 3=12=12,R R4 4=4=4R R5 5=12=12,R R6 6=4=4,R R7 7=2=27 7 7 7 矩阵方程简写为:矩阵方程简写为:令令u us s=10=10V V,求解矩阵方程得到,求解矩阵方程得到i ia a、i ib b、i ic c。再由再由i
5、i3 3=i=ia a-i-ib b ,u u4 4=R=R4 4i ib b ,u u7 7=R=R7 7i ic c 即可得即可得到问题到问题(1)(1)的解的解8 8 8 8 根据电路的线性性质,可令根据电路的线性性质,可令i i3 3=k=k1 1u us s,u u4 4=k=k2 2u us s,u,u7 7=k=k3 3u us s,由问题由问题(1)(1)的解的解 求得比例系数,进一步使问题求得比例系数,进一步使问题(2)(2)得得 到解答。到解答。具体根据问题具体根据问题(1)(1)的结果可列出以下的表达式:的结果可列出以下的表达式:因此,通过下列表达式即可求得问题因此,通过
6、下列表达式即可求得问题(2)(2)的解:的解:9 9 9 9B)B)MatlabMatlab程序程序(Ex01.m)(Ex01.m)c clear,lear,c close alose allll,format compactformat compactR1=2;R2=4;R3=12;R4=4;R5=12;R6=4;R7=2;R1=2;R2=4;R3=12;R4=4;R5=12;R6=4;R7=2;%为给定元件赋值为给定元件赋值displaydisplay(解问题解问题 (1)(1)%解问题解问题 (1)(1)a11=R1+R2+R3;a11=R1+R2+R3;a12=-R3;a12=-R3;
7、a13=0;a13=0;%将系数矩阵各元素赋值将系数矩阵各元素赋值a21=-R3;a21=-R3;a22=R3+R4+R5;a22=R3+R4+R5;a23=-R5;a23=-R5;a31=0;a31=0;a32=-R5;a32=-R5;a33=R5+R6+R7;a33=R5+R6+R7;b1=1;b2=0;b3=0;b1=1;b2=0;b3=0;us=inputus=input(给定给定 us=),us=),%输入解输入解 (1)(1)的已知条件的已知条件A=a11,a12,a13;a21,a22,a23;a31,a32,a33A=a11,a12,a13;a21,a22,a23;a31,a3
8、2,a33;%列出系数矩阵列出系数矩阵A AB=b1;0;0;B=b1;0;0;I=AB*us;I=AB*us;%I=ia;ib;icI=ia;ib;icia=I(1ia=I(1););ib=I(2);ib=I(2);ic=I(3);ic=I(3);i3=ia-ib,u4=R4*ib,u7=R7*ic i3=ia-ib,u4=R4*ib,u7=R7*ic%解出所需交量解出所需交量display(display(解问题解问题 (2)(2)%利用电路的线性性质及问题利用电路的线性性质及问题 (1)(1)的解求解问题的解求解问题 (2)(2)u u42=input(42=input(给定给定 u42
9、=);u42=);k k1=i3/us;1=i3/us;k2=u4/us;k2=u4/us;k3=u7/us;k3=u7/us;%由问题由问题 (1)(1)得出待求量与得出待求量与 us us 的比例系数的比例系数us2=u42/k2,us2=u42/k2,i32=k1/k2*u42,i32=k1/k2*u42,u72=k3/k2*u42u72=k3/k2*u42%按比例方法求出所需交量按比例方法求出所需交量10101010C)C)程序运行结果程序运行结果解问题解问题 (1)(1)给定给定 us=10us=10i3=0.3704 i3=0.3704 ,u4=2.2222 u4=2.2222 ,
10、u7=0.7407u7=0.7407解问题解问题 (2)(2)给定给定 u42=6u42=6us2=27.0000 us2=27.0000 ,i32=1.0000 i32=1.0000 ,u72=2u72=2 运行结果:运行结果:电路的解:电路的解:(1)(1)i i3 3=0.3704A,=0.3704A,u u4 4=2.2222V,=2.2222V,u u7 7=0.7404V=0.7404V(2)(2)us=27V,us=27V,i i3 3=1A,=1A,u u7 7=2V=2VEx01.m11111111补充说明:补充说明:实际中,实际中,如果熟悉列方程的方法如果熟悉列方程的方法
11、,那么在编写那么在编写MATLABMATLAB程序时可直接写出程序时可直接写出A A和和B B为为:从而可省去给元件和矩阵各元素赋值等语句。从而可省去给元件和矩阵各元素赋值等语句。12121212 例例例例2222 对如对如图图2 2所示的电路所示的电路,已知已知R R1 1=R=R2 2=R=R3 3=4=4,R R4 4=2=2,控制常数控制常数 K K1 1=0.5,k=0.5,k2 2=4,i=4,is s=2A,=2A,求求 i i1 1和和i i2 2。图图2 2 例例2 2的的电路电路13131313 对对图示电路图示电路,用节点电压法列写方程得:,用节点电压法列写方程得:A)A
12、)建模建模 解:解:解:解:uaub14141414 根据根据图示电路图示电路,控制变量,控制变量i i1 1、i i2 2与节点电与节点电压压u ua a、u ub b的关系为:的关系为:整理以上整理以上两式,两式,将将i i1 1、i i2 2也作为未知量也作为未知量,和前面的节点电压共同组成,和前面的节点电压共同组成方程方程 ,并写成矩阵形式并写成矩阵形式有:有:令令 i is s=2A=2A,求,求解上式即解上式即可可得得到到i i1 1和和i i2 2 。uaub15151515B)B)MatlabMatlab程序程序(Ex02.m)(Ex02.m)c clear,lear,form
13、at compactformat compactR R1 1=4=4;R2=4R2=4;R3=4R3=4;R4=2R4=2;%;%设置元件参数设置元件参数is=2is=2;k k1 1=0.5=0.5;k2=4k2=4;%按按A*X=B*isA*X=B*is列写电路的矩阵方程列写电路的矩阵方程,其中其中X=uaX=ua;u;ub b;i i1;1;i2i2。a a1111=1 1/R/R1 1+1 1/R2/R2;a a1 12=-2=-1 1/R2/R2;a a1 13=03=0;a a1 14=-k4=-k1;%1;%设置系数设置系数A Aa2a21 1=-=-1 1/R2/R2;a22=
14、a22=1 1/R2+/R2+1 1/R3+/R3+1 1/R4/R4;a23=-k2/R3a23=-k2/R3;a24=ka24=k1;1;a3a31 1=1 1/R2/R2;a32=-a32=-1 1/R2/R2;a33=-a33=-1;1;a34a34=0 0;a4a41 1=0=0;a42=a42=1 1/R4/R4;a43=0a43=0;a44=-a44=-1;1;A=a1A=a11 1,a,a1 12,a2,a1 13,a3,a1 14 4;a2a21 1,a22,a23,a24,a22,a23,a24;a3a31 1,a32,a33,a34,a32,a33,a34;a4a41 1
15、,a42,a43,a44,a42,a43,a44;B=B=1;1;0 0;0 0;00;%;%设置系数设置系数B BX=AB*is;X=AB*is;i1=X(3),i2=X(4)%i1=X(3),i2=X(4)%显示显示要求的分量要求的分量 16161616Ex02.mC)C)程序运行结果程序运行结果(电路的解电路的解)i i1 1=1 =1,i i2 2=1=117171717 例例例例3333 对如对如图图3 3所示的电路所示的电路,已知已知R R1 1=4=4,R,R2 2=2=2,R,R3 3=4=4,R R4 4=8=8;i is1 s1=2A=2A,i is2s2=0.5A=0.5
16、A。图图3 3 例例3 3的的电路电路(1)(1)负载负载R RL L为何值时能获得最为何值时能获得最 大功率大功率?(2 2)研究研究R RL L在在 001010范围内变范围内变 化时化时 ,其吸收功率的情况。其吸收功率的情况。18181818 解:解:解:解:用戴维南等效电路来求解。对图用戴维南等效电路来求解。对图3(a)3(a)电电路路,断开断开aoao,并在并在aoao端接入外电流源端接入外电流源i ia a,如图如图3 3(b)(b)所示。所示。以以o o为参考点列节点方程为参考点列节点方程得:得:A)A)建模建模图图3 3 例例3 3的的电路电路19191919前面的方程写成矩阵
17、形式为:前面的方程写成矩阵形式为:其中:其中:戴维南等效电路如图戴维南等效电路如图3(c)3(c)所示,其方程为:所示,其方程为:图图3 3 例例3 3的等效的等效电路电路20202020方法方法:令令 i ia a=0,=0,i is1 s1=2A,=2A,i is2s2=0.5A,=0.5A,由由矩阵方程求得矩阵方程求得u u11 11,u u2121,u ua1a1。因。因i ia a=0,=0,由由戴维南等效电路方程得:戴维南等效电路方程得:u uococ=u=ua1a1。再令再令i is1 s1=i=is2s2=0=0,i ia a=1A=1A,仍由仍由矩阵方程可求得另一组矩阵方程可
18、求得另一组u u1212,u u2222,u ua2a2。由于。由于内部电源内部电源i is1 s1=i=is s2 2=0 0,故故u uococ=0=0。从而。从而由由戴维南戴维南等效电路方程有:等效电路方程有:21212121 于是,原电路戴维南等效电路于是,原电路戴维南等效电路如图如图3(d)3(d)所示,负载所示,负载R RL L获得最大功获得最大功率时有率时有:图图3 3 例例3 3的等效的等效电路电路 至于至于问题问题(2),2),由图由图3 3(d d)可得可得R RL L吸收功率为:吸收功率为:再再令令 R RL L=l=l,2 2,3 3,1O1O,即可由上式分别,即可由上
19、式分别求求得得P PL L,并画图。并画图。22222222 可可设设ia ia为一个序列为一个序列(如如ia=0.1ia=0.1,0.20.2,2)2),计算相计算相应的应的u ua a序列序列,再,再用线性拟合用线性拟合,得出得出如下的如下的直线方程直线方程:方法方法:从而求得:从而求得:23232323B)B)MatlabMatlab程序程序 (Ex03-1.m)(Ex03-1.m)c clear,lear,format compactformat compactR R1 1=4;=4;R2=2;R2=2;R3=4;R3=4;R4=8;R4=8;%设置元件参数设置元件参数isis1=2;
20、is2=0.5;1=2;is2=0.5;%按按A*X=B*isA*X=B*is列写此电路的矩阵方程列写此电路的矩阵方程,其中其中X=u1;u2;ua;is=is1;is2;iaX=u1;u2;ua;is=is1;is2;iaa11=1/R1+1/R4;a12=-1/R1;a13=-1/R4;%a11=1/R1+1/R4;a12=-1/R1;a13=-1/R4;%设置系数矩阵设置系数矩阵a21=-1/R1;a22=1/R1+1/R2+1/R3;a23=-1/R3;a21=-1/R1;a22=1/R1+1/R2+1/R3;a23=-1/R3;a31=-1/R4;a32=-1/R3;a33=1/R3
21、+1/R4;a31=-1/R4;a32=-1/R3;a33=1/R3+1/R4;A=a11,a12,a13;a21,a22,a23;a31,a32,a33;A=a11,a12,a13;a21,a22,a23;a31,a32,a33;B=1,1,0;0,0,0;0,-1,1;%B=1,1,0;0,0,0;0,-1,1;%设置系数矩阵设置系数矩阵B B%方法方法:令令 ia=0,ia=0,求求uoc=x1(3);uoc=x1(3);再令再令is1=is2=0,is1=is2=0,设设ia=1,ia=1,求求Req=ua/ia=x2(3).Req=ua/ia=x2(3).Xl=AB*is1;is2;
22、0;uoc=X1(3)Xl=AB*is1;is2;0;uoc=X1(3)X2=AB*0;0;1;Req=X2(3)X2=AB*0;0;1;Req=X2(3)RL=Req;P=uoc2*RL/(Req+RL)2%RL=Req;P=uoc2*RL/(Req+RL)2%求最大负载功率求最大负载功率%也可设也可设RLRL为一数组为一数组,求出的负载功率也为一数组求出的负载功率也为一数组,画出曲线找极大值画出曲线找极大值RL=0:10,p=(RL*uoc./(Req+RL).*uoc./(Req+RL),%RL=0:10,p=(RL*uoc./(Req+RL).*uoc./(Req+RL),%设设RLR
23、L序列序列,求其功率求其功率figure(1),plot(RL,p),gridfigure(1),plot(RL,p),grid%画出功耗随画出功耗随RLRL变化的曲线变化的曲线xlabel(RL),ylabel(p)xlabel(RL),ylabel(p)24242424MatlabMatlab程序程序(Ex03-2.m)(Ex03-2.m)c clear,lear,format compactformat compactR R1 1=4;=4;R2=2;R2=2;R3=4;R3=4;R4=8;R4=8;%设置元件参数设置元件参数isis1=2;is2=0.5;1=2;is2=0.5;%按按
24、A*X=B*isA*X=B*is列写此电路的矩阵方程列写此电路的矩阵方程,其中其中X=u1;u2;ua;is=is1;is2;iaX=u1;u2;ua;is=is1;is2;iaa11=1/R1+1/R4;a12=-1/R1;a13=-1/R4;%a11=1/R1+1/R4;a12=-1/R1;a13=-1/R4;%设置系数矩阵设置系数矩阵a21=-1/R1;a22=1/R1+1/R2+1/R3;a23=-1/R3;a21=-1/R1;a22=1/R1+1/R2+1/R3;a23=-1/R3;a31=-1/R4;a32=-1/R3;a33=1/R3+1/R4;a31=-1/R4;a32=-1/
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 电路 中的 应用
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【a199****6536】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【a199****6536】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。