MATLAB 在电路分析中的应用

1、MATLAB语言课程论文         MATLAB 在电路分析中的应用         姓名：李娜 专业：2010级通信工程 班级：（1）班 指导老师：汤全武 学院：物理电气信息学院    MATLAB 在电路分析中的应用 摘 要 本文将Matlab软件的模拟功能用于电路分析研究，以基本电路理论中典型的直流电阻电路和含有复数运算的正弦稳态电路的计算为例，详述了如何分别运用MATLAB语言编程的方法来对电路进行仿真分析和计算。结论表明，应用这两种方法可以是复杂电路的分析和计算变得非常快捷

2、83;方便，从而为电路分析提供了一个有效的辅助工具。 关键词 MATLAB; 电路分析；模拟；正弦稳态；向量图 一、 问题的提出 MATLAB 语言结构紧凑·语句精炼，指令表达式和数字表达式非常接近，仅需几条简单的语句，就可以完成一大串其他高级语言才能完成的任务，可大大节省编程时间，提高计算效率。基本电路是电类专业非常重要的专业基本课，不仅为后继课程提供了深厚的理论基础，也为电路的分析计算提供了各种方法。其中，在电路分析理论中一般将关于时间的微分方程转化为复数方程求解，在一些电路比较复杂的·方程数量多的情况下，都可以运用MATLAB程序来解决。运用该程序不仅可以节约时间，还

3、可以非常方便的调试电路参数，直观的观察电路中的电流·电压和功率波形。 二、 应用 1 典型直流电阻电路的分析计算 图1所示为典型的直流电阻电路，含有电压控制的受控电流源VCCS，其中，R1=1,R2=2,R3=3,Us=10v,Is=15A,VCCS=0.25U2,现需分析计算电流i1 和电压u2图1 典型直流电阻电路 基本电路分析的基本方法实现建立数学模型，一般是电路方程组。然后通过求解方程组，得到各支路电压和电流。对图1应用回路电流法，可列出如下方程组： R11Im1+R12Im2+R13Im3=Us11 R21Im1+R22Im2+R23Im3=Us22 R31Im1+R32I

4、m2+R33Im3=Us33 其中,R11=R1+R2,R22=R1+R3, R33=R2+R3,R12=R21=-R1,R13=R31=-R2,R23=R32=-R3, US11=Us,Us22=U1,Us33=-U3 而I1=Im1-Im2,Im2=Is=15，Im3=0.25U2,U2=R2(Im1-Im3) 整理以上方程，并写出形如AX=BU的矩阵方程形式，可得R11 R13 0 0 Im1 1 -R12R21 R23 -1 0 Im3 = 0 -R22 US (1)R31 R32 0 1 U1 0 -R32 IS0.25R2 1-0.25R2 0 0 U3 0 0MATLA 语言编程

5、法 应用MATLAB语言编程如下： CLEAR; US=10;IS=15;R1=1;R2=2;R3=3; % 为给定元件赋值 R11=R1+R2;R12=-R1;R21=-R1;R13=-R2;R31=-R2; % 为系数矩阵各元素赋值 R22=R1+R3;R23=-R3;R32=-R3;R33=R2+R3; A=R11 R13 0 0;R21 R23 -1 0;R31 R33 0 1;0.25*R2 1-0.25*R2 0 0; % 列出系数矩阵A B=1 -R12;0 -R22;0 -R23;0 0;USS=US;IS; % 列出系数矩阵B X=AB*USS; % 解出X I1=X(1)-

6、IS % 显示要求的分量I1和U2 U2=2*(X(1)-X(2) 程序运行结果 I1=-10.0000 , U2=202 典型的正弦稳态电路的分析与计算 图2所示为典型的正弦稳态电路，其中现需分析该含源一端口在b-o端口间戴维南等效电路  图2 典型的正弦稳态电路图3 在b-o端口间外加电流源后的电路  首先建立数学模型。我们在原含源一端口电路的b-o端子间外加一个正弦电流源，如图3所示。对图3应用结点电压法，并以o点为参考结点，则有如下方程组： Y11úao+Y21úbo=ús11 Y12úao+Y22úbo=ú

7、;s22 其中， 整理以上方程，并转换成形如AX=BU的矩阵方程形式为：  MATLAB 语言编程法实现电路的分析计算 根据式(2)，我们设想，若令íb=0，代入 ús=10-45, 则可求得戴维南等效电源电压ú OC ，它就等于此时的úbo ；然后再令ús=0，将原电路（图2）变成一个无源一端口，并设 íb=10 ，代入式(2)即可求得戴维南等效阻抗,即据此，可设计MATLAB程序。 应用MATLAB 语言编程如下： clear; R1=1;R2=2;L1=4e-4;C1=1e-3;US=5*sqrt(2)-j*5*sqr

8、t(2); % 为给定元件赋值 W=1000; ZR1=1;ZR2=2;ZL1=j*W*L1;ZC1=1/(j*W*C1); Y11=1/(ZR1+ZC1)+1/ZL1+1/ZR2;Y22=1/ZR2; % 为系数矩阵各元素赋值 Y12=-1/ZR2;Y21=-1/ZR2; A=Y11 Y21;(Y12-0.5) Y22;B=1/(ZR1+ZC1) 0;0 1; % 列出各系数矩阵 X0=AB*US;0; % 戴维南等效电源电压UOC等于 b=0, s=2 0 时的Ubo, 是一个复数UOC=X0(2), uoc=abs(UOC),uang=angle(UOC) % 求戴维南等效电源电压的模和

9、辐角 X1=AB*0;1; % 再令 s=0，并设 b=1 0 ，求戴维南等效阻抗Ze Zeq =X1(2) ze=abs(Zeq),zang=angle(Zeq) % 求戴维南等效阻抗Zeq的模和辐角 程序运行结果 UOC =0.0000 + 7.0711i uoc = 7.0711 uang = 1.5708 Zeq = 2.0000 + 1.0000i ze = 2.2361 zang = 0.4636 3 向量与电路 图4 电路图  电路如图4所示，其中的求各支路电流并画向量图。 这是一个交流稳态电路，对二个独立结点列结点电压方程： Y11U1 +Y12U2 =IS1 Y21

10、U1 +Y22U2 =IS2 其中：Y11=G2+G3;Y12=-(G2+G3+G5) Y21=G1+G2+G3+G4;Y22=-(G2+G3) IS1=G5US2;IS2=G1US1 G1=1/R1;G2=1/(R2-jx2);G3=1/-jx3; G4=1/jx4;G5=1/R3. 用Matlab语言编程实现上述计算，程序如下： R1=4 ；R2=3；R3=1；X1=2；X2=0.1；X3=0.8；US1=12;US2=8; 输入初始参数 G1=1 R1；G2=1(R2一j*X2)； G3=1 -j*X3；G4=1j*X1；G5=IR3； Y11=G2+G3 ；Y12=-(G2+G3+G5

11、)； Y21=G1+G2+G3+G4 ；Y22=-(G2+G3)； IS1=G5*US2 ；IS2=G1*US1 % 计算线性方程组系数矩阵中以上各元素的值 A=Y11 ，Y12；Y21，Y22 B=Is1;Is2 % 组成方程组A、BU=A B %解结点电压I1=G1*(U(1) 一US1) %求支路电流I1I2=G2*(U(1) 一U(2) %求支路电流I2I3= G3*(U (1) 一lJ(2) %求支路电流I3 I4=G4*U(1) %求支路电流I4I5=G5*(U(2) 一US2) %求支路电流I5程序运行结果为： I1=-2.9841+0.9862i I2= -1.0625+1.5

12、658i I3= -3.8430-2.4247i I4=7.8895-0.1273i I5=-4.9054-0.8589i 三、MATLAB 应用在电路稳态分析 1 直流稳态分析实例 在图5所示电路中，,求U10 . 图5 直流稳态分析用的实例求解此题的方程组为 对应的M 文件为 A=7 -2 0;-3 2 0;1 0 1; % 定义方程组的系数矩阵A B=16 0 16; % 定义右端矩阵B C=AB % 求解未知变量矩阵C C=4.0000 6.0000 12.0000 % 此为U10 值 2 交流稳态分析 在图6所示的电路中，用2b法求各支路的变量(本例中只比较R6 上的电压)。 图6

13、交流稳态分析的实例与图6对应的2b方程的矩阵形式为 矩阵方程中各子阵的列可写成下面给出其M文件: A=-1 0 0 0 0 1 0 0 0;-1 1 0 0 0 1 0 0;0 1 0 0 0 -1 0 1 1;0 0 0 1 0 0 -1 0 0;0 0 0 0 -1 0 0 0 -1 %输入矩阵AB=1 1 0 0 0 1 0 0 0;0 0 1 -1 0 0 -1 0 0;0 -1 -1 0 0 0 0 1 0; 0 0 0 0 -1 0 0 -1 1 %输入矩阵BC=-1 0 -1 j*1E -4 0 -1 1 1 j*2E -4; %输入矩阵C Ye=diag(C) Ye(5,6)=

14、-1 %产生零矩阵为-1D=j*1 0 0 1 1 -1 1 1 0 -1 -1; Ze=diag(D) Us=0 0 0 0 0 0 1 0+j*0 0 0 ; % 输入电压UsIs=0 -2+j*0 0 0 0 0 0 0 0 ; %输入电流IsE=zeros(5,9) %产生零矩阵EF=zeros(4,9) %产生零矩阵FG=0 0 0 0 0' %输入矩阵GH=0 0 0 0' %输入矩阵 HW=E A;B F;Ye Ze %输入矩阵WN=G;H;Us+Is %输入矩阵NXn=WN %求解支路电压第6条 支路的电压向量为1.0e+002*(-0.0000+0.0004i

15、); 计算其峰值为：0.05656V。 3 MATLAB 应用在电路暂态分析 图7所示的电路中，开关s闭合前已达稳定状态。已知：求开关s在时间t=0瞬时闭合后，电感支路上的电流iL(t).图7 暂态分析所用的电路此题求解的二阶微分方程如下： 对应的M文件为 Desolve('D2y+10*Dy+10*y=1000','Dy(0)=100','Y(0)=0') Ans=100-101.6398*exp(-1.1270*t)+1.6398*exp(-8.8730*t) 其解为 iL(t)=100-101.6398e-1.1270t+1.6398e-8

16、.8370t 由此例看出，用MATLAB自身提供的数值微分函数dsolve求解微分方程简便快捷，大大节省了编程时间,采用同一算法的Fortran语言和C语言程序却多达百条. 四 、结论 本文通过基本电路理论中的典型题目介绍了如何应用MATLAB语言编程的方法来对复杂的电路进行分析和计算。该方法不仅可以节约计算时间，方便的调试电路参数，而且还可以非常直观地观察和测量电路中的电压，电流功率等物理量。结论表明，MATLAB提供了高效简洁的编程方法，其强大而简洁的绘图功能，矩阵和数组的运算能力以及很强的扩充性，能充分满足基本电路分析，计算的需要，从而可以大大的提高计算精度和工作效率，在电路理论学科研究与工程实践中具有很好的应用价值。五、课程体会 经过一学期紧张而有序的课程学习，在忙碌之余也得到了颇多的体会。我深深体会到MATLAB语局简练，功能强大，简单实用，用途广泛，不仅可以大大的提高操作效率，缩短编程时间，是一种简单实用的工具，而且还可以应用于其他学科领域，此次在电路分析中，它有效又简洁地解决了许多的复杂电路问题，给我带来了许多的方便。正是由于我的任课老师汤全武老师的精彩授课和认真的讲解，使我学到了更多的MATLAB语言的知识，并且更好的应用于生活学习中非常感谢汤老师这一学期的教育，愿MATLAB语言有着更广泛的应用前景参考文献： 1 邱关源电路(第三版)M北京：高等教育出版社，1