matlab应用举例

《matlab应用举例》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、MATLAB陈嘉2008301130070实验一直流电路（1）一实验目的1加深对直流电路的节点电压法和网孔电流法的理解2学习使用MATLAB的矩阵运算的方法二实验示例1节点分析根据电路图得到矩阵方程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回路分析z=[40,-10,-30;-10,30,-5;-30,-5,65];v=[10,0,0]';I=inv(z)*v;

2、IRB=I(3)-I(2);fprintf('thecurrentthroughRis%8.3fAmps',IRB)ps=I(1)*10;fprintf('thepowersuppliedby10vsourceis%8.4fwatts',ps)输出thecurrentthroughRis0.037Ampsthepowersuppliedby10Vsourceis4.7531watts三实验内容1(1)如果Us=10V,求i3,u4,u7(2)如果U4=4V,求Us,i3,i7使用matlab命令为(1）cle

3、ar%初始化矩阵ZH和向量VZ=[20-120;-1232-12;0-1218];V=[1000]';%解答回路电流I=inv(Z)*V;%I3的计算I3=I(1)-I(2);fprintf('thecurrentI3is%8.2fAmps',I3)%U4的计算U4=8*I(2);fprintf('thevoltageU4is%8.2fVmps',U4)%U7的计算U7=2*I(3);fprintf('thevoltageU7is%8.2fVmps',U7)结果thecurrentI3is0.36Amps

4、thevoltageU4is2.86VmpsthevoltageU7is0.48Vmps（2）clear%初始化矩阵Z和向量UZ=[20-10;-120-12;0018];U=[6-166]';%进行解答I=inv(Z)*U;%计算各要求量Us=I(2)I3=I(1)-0.5I7=I(3)结果Us=14.0000I3=0.5000I7=0.33332clear%初始化节点电压矩阵U和向量IU=[0.725-0.125-0.1-5-1.25;-0.1-0.20.5500;-0.1250.325-0.201.25;10-

5、1-10;00.2-0.201];I=[06500]';%解答节点电压U1，U3，U4与Vb，IaX=inv(U)*I;%所求电压计算V1=X(1)V2=X(1)-10*X(5)V3=X(2)V4=X(3)V5=24结果V1=117.4792V2=299.7708V3=193.9375V4=102.7917V5=243clear%初始化节点电压矩阵Z和向量IZ=[0.5-0.250-0.5;-0.251-10.5;00.50-1;1-1-40];I=[2000]';%解答节点电压V1，V2及电流I1，I2X=inv(

6、Z)*I;%计算未知数V1=X(1)V2=X(2)I1=X(3)I2=X(4)结果如下：V1=6V2=2I1=1I2=1实验二直流电路（2）一实验目的1加深多戴维南定律，等效变换等的了解2进一步了解matlab在直流电路中的作用二实验示例如图所示（图见书本17页2-1），分析并使用matlab命令求解为clear,formatcompactR1=4;R2=2;R3=4;R4=8;%设置元件参数is1=2;is2=0.5;a11=1/R1+1/R4;a12=-1/R1;a13=-1/R4;%设置系数矩阵Aa21=-1/

7、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];%方法1：令ia=0，求uoc=x1（3）；再令is1=is2=0，设ia=1，求Req=ua/ia=x2（3）。X1=AB*[is1;is2;0];uoc=X1(3);X2=AB*[0;0;1];Req=X2(3);RL=Req;P=uoc^2*RL/(Re

8、q+RL)^2;%求最大负载功率RL=0:10,p=(RL*uoc./(Req+RL)).*uoc./(Req+RL),%设RL序列，求其功率figure(1),plot(RL,p),grid%方法2：设一个ia序列，计算一个ia序列，用线性拟合求出其等效开路电压和等效电阻fork=1:21ia(k)=(k-1)*0.1;X=AB*[is1;