实验七：MATLAB在电工与电子线路中的应用.doc

实验七：MATLAB在电工与电子线路中的应用一、实验目的掌握用MATLAB解决电工与电子线路中的难题。二、实验仪器1、计算机2、MATLAB 软件环境三、实验内容1、计算电路中A、B 两点的电位。C点为参考点。AI1B5W5W+15V10W10W15W+-65VI2I4I5C2、电路如图，t=0时合上开关S，合S前电路已处于稳态。试求电容电压uc 和电流ic、i2 并画出它们的波形。S9mA6kW2mF3kWt=0+-C R3、+-4、+-+-5、设将一个二极管与一电阻Rf串接，在此电路的两端加上正向直流电压U0，如图8-15所示，试求出此电路中的电流Idx和电压Udx。四、实验过程1. 建立函数f1clear close allr1=5;r2=10;r3=15;r4=10;r5=5;a11=1/r1+1/r2+1/r5;a12=-1/r2;a13=1/r1;a21=-1/r2;a22=1/r2+1/r3+1/r4;a23=1/r3;us=input(us=)uv=input(uv=)A1=a11,a12;a21,a22A2=a13;a23u=A1a13*us;a23*uvendf1u = 10.0000 20.00002.建立函数f2clear close allr1=6000; r2=3000;c=0.000002;ui=0.009; uc0=54;ir10=uc0/r1;ir20=uc0/r2;ic0=ui-ir10-ir20;ir1f=ui*r2/(r1+r2);ir2f=ui*r1/(r1+r2);ucf=ir1f*r1;icf=0;t=-2:0-eps,0:15;uc(1:3)=54;ir2=0;ic=0;T=r1*r2/(r1+r2)*c;uc(4:19)=ucf+(uc0-ucf)*exp(-t(4:19)/T);ir2(4:19)=ir2f+(ir20-ir2f)*exp(-t(4:19)/T);ic0(4:19)=icf+(ic0-icf)*exp(-t(4:19)/T);subplot(3,1,1);h1=plot(t,uc),subplot(3,1,2);h2=plot(t,ic),subplot(3,1,3);h3=plot(t,ir2);endf2得图：如下 ：3.建立程序f3clear，close allw=input(w=);t=input(t=);z1=50;z2=100+200*j;z3=400*j;U=220*sqrt(2)*sin(w*t);z23=z2*z3/(z2+z3);z=z23+z1;I=U/z;I1=U*z23/(z*z2);I2=U*z23/(z*z3);i=I/sqrt(2),i1=I1/sqrt(2),i2=I2/sqrt(2)disp( i i1 i2 )disp(幅值), disp(abs(i, i1, i2)disp(相角), disp(angle(i, i1, i2)*180/pi)ha=compass(i,i1,i2); set(ha, linewidth, 2);endf3w=2t=1i = 0.6557 - 0.9883ii1 = 0.5322 - 0.5702ii2 = 0.1235 - 0.4181ii i1 i2 幅值 1.1861 0.7800 0.4360相角 -56.4373 -46.9749 -73.5400实验截图如下所示：4建立函数f4clear， close alli3=21.7455-22.6067ix1=0.1+0.5*j;z2=0.1+0.5*j;z3=5+5j;u1=230;u2=227;a11=z1+z3;a12=z3;a13=1;a21=z3;a22=z3+z2;a23=1;A=a11,a12;a21,a22;i=Au1*a13;u2*a23;i3=i(1)+i(2)endf4i3 = 21.7455 -22.6067i5.建立函数f5clear，close allK=1.38e-23; T=300; q=1.6e-19; KT=K*T/q; Is=10e-12; Ud=0:0.01:3.5; Id=Is*(exp(Ud/KT)-1); plot(Ud, Id), grid onaxis(0, max(Ud), 0, 100), hold on line(1.5, 1.8, 76, 76); fill(1.8, 2, 2, 1.8, 76, 72, 80, 76, K); line(1.8, 1.8, 72, 80, linewidth, 2);line(2, 2.5, 76, 76);line (2.5,2.8,2.8,2.5,2.5,74,74,78,78,74,linewidth,2) ; line(2.8, 3.1, 76, 76);plot(1.5, 2.2, 3.1, 76, 76, 76, o);text(1.4, 70, o); text(2.1, 70, Ud); text(2.6, 68, Rf); text(3, 70, U0); U0=input(U0=伏)Rf=input(Rf=欧姆)Id1=1000*(U0-Ud)./Rf plot(Ud, Id; Id1), grid on di, nI=min(abs(Id-Id1); Udx=Ud(nI); Idx=Id1(nI); disp(Udx, Idx=), Udx, Idx, hold offlege