武大MATLAB电路仿真实验报告

电路MATLAB仿真报告武汉大学电气工程学院2015级班刘》》2015302540》》MATLAB电路仿真实验报告1实验一直流电路（1）一实验目的1加深对直流电路的节点电压法和网孔电流法的理解2学习使用MATLAB的矩阵运算的方法二实验示例1节点分析电路如图所示（见书本12页），求节点电压V1,V2,V3.根据电路图得到矩阵方程，根据矩阵方程使用matlab命令为Y=0.1500-0.1000-0.0500-0.10000.1450-0.0250-0.0500-0.02500.0750节点v1，v2和v3：v=404.2857350.0000412.85712回路分析电路如图所示（见书本13页），使用解析分析得到同过电阻RB的电流，另外求10V电压源的输出功率。MATLAB电路仿真实验报告2分析电路得到节点方程，根据节点方程得到矩阵方程，根据矩阵方程，使用matlab的命令为z=[40,-10,-30;-10,30,-5;-30,-5,65];v=[10,0,0]';I=inv(z)*v;IRB=I(3)-I(2);fprintf('thecurrentthroughRis%8.3fAmps\n',IRB)ps=I(1)*10;fprintf('thepowersuppliedby10vsourceis%8.4fwatts\n',ps)结果为：thecurrentthroughRis0.037Ampsthepowersuppliedby10Vsourceis4.7531watts三实验内容1根据书本15页电路图，求解电阻电路，已知：R1=2Ω,R2=6Ω,R3=12Ω,R4=8Ω,R5=12Ω,R6=4Ω,R7=2Ω(1)如果Us=10V,求i3,u4,u7(2)如果U4=4V,求Us,i3,i7MATLAB电路仿真实验报告3使用matlab命令为clear%初始化阻抗矩阵Z=[20-120;-1232-12;0-1218];%初始化电压矩阵V=[1000]';%解答回路电流I=inv(Z)*V;%I3的计算I3=I(1)-I(2);fprintf('thecurrentI3is%8.2fAmps\n',I3)%U4的计算U4=8*I(2);fprintf('thevoltageU4is%8.2fVmps\n',U4)%U7的计算U7=2*I(3);fprintf('thevoltageU7is%8.2fVmps\n',U7)MATLAB电路仿真实验报告4结果thecurrentI3is0.36AmpsthevoltageU4is2.86VmpsthevoltageU7is0.48Vmpsclear%初始化矩阵XX=[20-10;-120-12;0018];%初始化矩阵YY=[6-166]';%进行解答A=inv(X)*Y;%计算各要求量Us=A(2)I3=A(1)-0.5I7=A(3)结果MATLAB电路仿真实验报告5Us=14.0000I3=0.5000I7=0.33332求解电路里的电压如图1-4（书本16页），求解V1,V2,V3,V4,V5使用matlab命令为clear%初始化节点电压方程矩阵Z=[0.725-0.125-0.1-5-1.25;-0.1-0.20.5500;-0.1250.325-0.201.25;10-1-10;00.2-0.201];I=[06500]';%解答节点电压U1，U3，U4与Vb，IaA=inv(Z)*I;%最终各电压计算V1=A(1)V2=A(1)-10*A(5)V3=A(2)MATLAB电路仿真实验报告6V4=A(3)V5=24结果V1=117.4792V2=299.7708V3=193.9375V4=102.7917V5=243如图1-5（书本16页），已知R1=R2=R3=4Ω,R4=2Ω,控制常数k1=0.5，k2=4，is=2A，求i1和i2.使用matlab命令为clear%初始化节点电压方程矩阵Z=[0.5-0.250-0.5;-0.251-10.5;00.50-1;1-1-40];I=[2000]';%解答节点电压V1，V2及电流I1，I2MATLAB电路仿真实验报告7A=inv(Z)*I;%计算未知数V1=A(1)V2=A(2)I1=A(3)I2=A(4)结果如下：V1=6V2=2I1=1I2=1实验二直流电路（2）一实验目的1加深多戴维南定律，等效变换等的了解2进一步了解matlab在直流电路中的作用二实验示例如图所示（图见书本17页2-1），分析并使用matlab命令求解为MATLAB电路仿真实验报告8clear,formatcompactR1=4;R2=2;R3=4;R4=8;is1=2;is2=0.5;a11=1/R1+1/R4;a12=-1/R1;a13=-1/R4;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];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:10,p=(RL*uoc./(Req+RL)).*uoc./(Req+RL),figure(1),plot(RL,p),gridfork=1:21ia(k)=(k-1)*0.1;X=A\B*[is1;is2;ia(k)];u(k)=X(3);endfigure(2),plot(ia,u,'x'),gridc=polyfit(ia,u,1);%ua=c(2)*ia=c(1),用拟合函数术,c(1),c(2)uoc=c(1),Req=c(2)MATLAB电路仿真实验报告9RL=012345678910p=Columns1through700.69441.02041.17191.23461.25001.2397Columns8through111.21531.18341.14801.1111A．功率随负载变化曲线01234567891000.20.40.60.811.21.4MATLAB电路仿真实验报告10B.电路对负载的输出特性00.20.40.60.811.21.41.61.8246810121416功率随负载的三实验内容1图见书本19页2-3，当RL从0改变到50kΩ，校验RL为10kΩ的时候的最大功率损耗使用matlab命令为clear%定义电压源和电阻值Us=10;Rs=10000;MATLAB电路仿真实验报告11RL=0:20000;p=(Us^2.*RL)./(RL+Rs).^2;plot(RL,p);输出结果为00.20.40.60.811.21.41.61.82x10400.511.522.5x10-3Maximumpoweroccurat10000.00hmsMaximumpowerdissipationis0.0025Watts2在图示电路里(书本20页2-4)，当R1取0，2，4，6，10，18，24，42，90和186Ω时，求RL的电压UL,电流IL和RL消耗的功率。使用matlab命令为：clear%设置元件参数RL=[02461018244290186];%列出要求的参数同元件间关系式以得出结果UL=48*RL./(RL+6)MATLAB电路仿真实验报告12IL=48./(RL+6)p=2304*RL./(RL+6).^2%画出要求参数随RL变化的曲线plot(RL,UL,'r+')holdonplot(RL,IL,'m*')holdonplot(RL,p,'ks')结果数据UL=Columns1through7012.000019.200024.000030.000036.000038.4000Columns8through1042.000045.000046.5000IL=Columns1through78.00006.00004.80004.00003.00002.00001.6000Columns8through10MATLAB电路仿真实验报告131.00000.50000.2500p=Columns1through7072.000092.160096.000090.000072.000061.4400Columns8through1042.000022.500011.6250UL=Columns1through7012.000019.200024.000030.000036.000038.4000Columns8through1042.000045.000046.5000IL=Columns1through78.00006.00004.80004.00003.00002.00001.6000Columns8through101.00000.50000.2500p=Columns1through7072.000092.160096.000090.000072.000061.4400Columns8through10MATLAB电路仿真实验报告1442.000022.500011.6250UL=Columns1through7012.000019.200024.000030.000036.000038.4000Columns8through1042.000045.000046.5000IL=Columns1through78.00006.00004.80004.00003.00002.00001.6000Columns8through101.00000.50000.2500p=Columns1through7072.000092.160096.000090.000072.000061.4400Columns8through1042.000022.500011.6250UL=Columns1through7012.000019.200024.000030.000036.000038.4000Columns8through10MATLAB电路仿真实验报告1542.000045.000046.5000IL=Columns1through78.00006.00004.80004.00003.00002.00001.6000Columns8through101.00000.50000.2500p=Columns1through7072.000092.160096.000090.000072.000061.4400Columns8through1042.000022.500011.6250UL=Columns1through7012.000019.200024.000030.000036.000038.4000Columns8through1042.000045.000046.5000IL=Columns1through78.00006.00004.80004.00003.00002.00001.6000Columns8through10MATLAB电路仿真实验报告161.00000.50000.2500p=Columns1through7072.000092.160096.000090.000072.000061.4400Columns8through1042.000022.500011.6250实验三正弦稳态一实验目的1学习正弦交流电路的分析方法2学习matlab复数的运算方法二实验示例1如图3-1（书本21页），已知R=5Ω,ωL=3Ω，1/ωc=2Ω，uc=10∠30°V,求Ir，Ic，I和UL,Us，并画出其向量图。使用matlab命令为:Z1=3*j;Z2=5;Z3=-2j;Uc=10*exp(30j*pi/180);Z23=Z2*Z3/(Z2+Z3);Z=Z1+Z23;Ic=Uc/Z3,Ir=Uc/Z2