电路Matlab仿真实验报告

1、武汉大学电气工程学院MATLAB电路仿真实验报告姓名：班级：学号：目录实验一 直流电路（1）4一、实验目的4二、预习要求4三、实验例题4四、实验内容6五、实验总结9实验二 直流电路（2）10一、实验目的10二、预习要求10三、实验示例10四、实验内容13五、实验总结18实验三、正弦稳态19二、预习要求19三、实验示例19四、实验内容23实验四、交流分析和网络函数28一、实验目的28二、预习要求28三、实验示例28实验五、动态电路32一、实验目的32二、预习要求32三、实验示例32四、实验内容37五、实验总结45实验六、频率响应46一、实验目的46二、预习要求46三实验示例46实验七 simul

2、ink仿真交流电路54一、实验目的54二、实验内容54实验一 直流电路（1）一、实验目的 1.加深对直流电路的节点电压法和网孔电流法的理解。 2.学习使用MATLAB的矩阵运算的方法。二、预习要求1.复习基尔霍夫KCL和KVL方程及直流电路的相关内容。2.熟悉前面有关矩阵预算的内容。三、实验例题1.节点分析 示例1.1电路如下图所示，求节点电压V1，V2，V3。MATLAB求解： Y = 0.15 -0.1 -0.05;-0.1 0.145 -0.025;-0.05 -0.025 0.075 ;I = 5; 0; 2 ;fprintf('节点V1,V2和V3: n')v = i

3、nv(Y)*I仿真结果： 节点V1,V2和V3: v = 404.2857 350.0000 412.85712、回路分析 示例1.2使用解析分析得到通过电阻RB的电流。另外，求10V电压源提供的功率。解：分析电路得到节点方程，根据节点方程得到矩阵方程，根据矩阵方程，MATLAB求解：Z=40,-10,-30; -10,30,-5; -30,-5,65; V=10,0,0' I=inv(Z)*vV;IRB=I(3)-I(2); fprintf('the current through R 

4、is %8.3f Amps n',IRB) PS=I(1)*10; fprintf('the power supplied by 10V source is %8.4f wattsn',PS) 仿真结果为：the current through R is    0.037 Amps the power supplied by 10V 

5、source is   4.7531 watts四、实验内容编写以下程序，并记录程序和结果，写出简单注释。 1. 电阻电路的计算如图，已知：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,i7；解：（1）解：MATLAB求解：Z = 20 -12 0;-12 32 -12; 0 -12 18;V = 10 0 0'I = inv(Z)*V;i3 = I(1)-I(2);u4 = 8*I(2);u7 = 2*I(3);fp

6、rintf('i3=%f n',i3)fprintf('u4=%f n',u4)fprintf('u7=%f n',u7)仿真结果： i3=0.357143 u4=2.857143 u7=0.476190（2）解：Matlab求解：Z = 0 8 0;-12 32 -12; 0 -12 18;V = 4 0 0'I = inv(Z)*V;Us = 20*I(1)-12*I(2);i3 = I(1)-I(2);i7 = I(3);fprintf('Us=%f n',Us)fprintf('i3=%f n',i

7、3)fprintf('i7=%f n',i7)仿真结果：Us=14.000000 I3=0.500000 i7=0.333333 2.求解电路里的电压，例如V1，V2V5Y = 1  -1   2   -2  0;0   5  - 13  8  0;2  0   4  -11  0 ;176  -5  5 -196  0;0  0  0  0  1;I 

8、;= 0 -200 -120 0 24'V = inv(Y)*I;fprintf('V1=%fVnV2=%fVnV3=%fVnV4=%fVnV5=%fVn',V(1),V(2),V(3),V(4),V(5) 仿真结果：V1=117.479167VV2=299.770833VV3=193.937500VV4=102.791667VV5=24.000000V3.如图所示，已知R1=R2=R3=4 Ohms, R4=4 Ohms,控制系数K1=0.5， K2=4，is=2A,求i1和i2。解：Ma

9、tlab求解：Z = 1 0 0 0; -4 16 -8 -4; 0 0 1 0.5; 0 -8 4 6;V = 2 0 0 0'I = inv(Z)*V;i1 = I(2)-I(3);i2 = I(4);fprintf('i1=%f Vni2=%f Vn',i1,i2)仿真结果：i1=1.000000 Vi2=1.000000 V五、实验总结1、仿真前需进行准确的计算，列出节点或回路方程进而进行矩阵计算。2、熟练矩阵运算公式，即：V=inv（Y）*I实验二 直流电路（2）一、实验目的1.加深对戴维南定理、等效变换的了解。2.进一步了解MATLAB在直流电路中的应用。二

10、、预习要求1.复习戴维南定理等直流电路的相关原理。2.了解MATLAB变量生成的应用。3.了解数组的运算。三、实验示例1.戴维南定理如图所示电路，已知R1=4 Ohms， R2=2 Ohms， R3=4 Ohms；R4=8 Ohms；is1=2A，is2=0.5A.（1）.负载RL为何值时能获得最大功率？（2）.研究RL在010 Ohms范围内变化时，其吸收功率的变化情况。Matlab求解：clear,format compactR1=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/R4;a22=1

11、/R1+1/R2+1/R3;a23=-1/R4;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=AB*is1;is2;0;uoc=X1(3)X2=AB*0;0;1;Req=X2(3)RL=Req;P=uoc2*RL/(Req+RL)2RL=0:10,p=(RL*uoc./(Req+RL).*uoc./(Req+RL),figure(1),plot(RL,p),gridfor k=1:21 ia(k)=(k-1)*0.1; X=AB*is1;is2;ia

12、(k); u(k)=X(3);endfigure(2),plot(ia,u,'x'),gridc=polyfit(ia,u,1);仿真结果：uoc = 2.3333Req = 3.6667P = 0.3712RL = 0 1 2 3 4 5 6 7 8 9 10p = Columns 1 through 7 0 0.2500 0.3391 0.3675 0.3705 0.3624 0.3496 Columns 8 through 110.3350 0.3200 0.3054 0.2915（a）功率随负载的变化曲线（b）电路对负载的输出特性四、实验内容1. 在下图中，当RL从0改

13、变到50千欧，绘制负载功率损耗。校验当RL=10千欧时的最大功率损耗。解：Matlab求解：R=10;U=10;RL=10;P=U2*(RL*1000)/(R+RL)*1000)2RL=0:50;p=(RL*1000*U./(R+RL)*1000).*U./(R+RL)*1000)figure(1),plot(RL,p),grid仿真结果：P = 0.0025p = Columns 1 through 7 0 0.0008 0.0014 0.0018 0.0020 0.0022 0.0023 Columns 8 through 14 0.0024 0.0025 0.0025 0.0025 0.

14、0025 0.0025 0.0025 Columns 15 through 21 0.0024 0.0024 0.0024 0.0023 0.0023 0.0023 0.0022 Columns 22 through 28 0.0022 0.0021 0.0021 0.0021 0.0020 0.0020 0.0020 Columns 29 through 35 0.0019 0.0019 0.0019 0.0018 0.0018 0.0018 0.0018 Columns 36 through 42 0.0017 0.0017 0.0017 0.0016 0.0016 0.0016 0.00

15、16 Columns 43 through 49 0.0016 0.0015 0.0015 0.0015 0.0015 0.0014 0.0014 Columns 50 through 51 0.0014 0.00142. 在图示电路中，当R1取0,2,4,6,10,18,24,42,90和186欧时，求RL的电压UL，电流IL和RL消耗的功率。解：Matlab求解：A=3/4 -1/2 0; 1/2 -33/24 5/6; 0 1 -1;I=15 0 0' U=inv(A)*I; us=U(3); R=6;Z=0 2 4 6 10 18 24 42 90 186;RL=Z(1,:),

16、 i=us./(R+RL) u=us.*RL./(R+RL) p=(RL.*us./(R+RL).*us./(R+RL)figure(1),plot(RL,i),gridfigure(2),plot(RL,u),gridfigure(3),plot(RL,p),grid仿真结果：RL = 0 2 4 6 10 18 24 42 90 186i = Columns 1 through 7 8.0000 6.0000 4.8000 4.0000 3.0000 2.0000 1.6000 Columns 8 through 10 1.0000 0.5000 0.2500u = Columns 1 t

17、hrough 7 0 12.0000 19.2000 24.0000 30.0000 36.0000 38.4000 Columns 8 through 10 42.0000 45.0000 46.5000p = Columns 1 through 7 0 72.0000 92.1600 96.0000 90.0000 72.0000 61.4400 Columns 8 through 10 42.0000 22.5000 11.6250功率P、电压U、电流I随电阻R的变化曲线分别如下图所示五、实验总结1、经过这次实验基本了解了MATLAB变量生成的应用。2、经过这次实验更加深刻了戴维南等效电

18、路的原理。3、了解了MATLAB中图像的生成。实验三、正弦稳态一、实验目的1. 学习正弦交流电路的分析方法2. 学习MATLAB复数的运算方法二、预习要求3. 复习有关正弦交流电路的有关概念4. 了解有关相量图的绘制三、实验示例 1.如图所示电路，已知欧姆，欧姆，1/wC=2欧姆，uc=1030 V, 求Ir，Ic，I，Ul，Us，并画出其相量图。解：Matlab求解：Z1=3j;Z2=5;Z3=-2j;Uc=10*exp(30j*pi/180);Z23=Z2*Z3/(Z2+Z3);Z=Z1+Z23;Ic=Uc/Z3,Ir=Uc/Z2,I=Ic+Ir,U1=I*Z1,Us=I*Z;disp(&

19、#39;Uc Ir Ic I u1 Us')disp('·ùÖµ'),disp(abs(Uc,Ir,Ic,I,U1,Us)disp('Ïà½Ç'),disp(angle(Uc,Ir,Ic,I,U1,Us)*180/pi)ha=compass(Uc,Ir,Ic,I,Us,Uc);set(ha,'linewidth',3)仿真结果：Ic = -2.5000 + 4.3301iIr = 1.7321 + 1.0000iI = -0.7679 + 5.3301iU1

20、= -15.9904 - 2.3038iUc Ir Ic I u1 Us幅值 10.0000 2.0000 5.0000 5.3852 16.1555 7.8102相角 30.0000 30.0000 120.0000 98.1986 -171.8014 159.80562.正弦稳态电路，戴维南定理 如图所示电路，已知C1=0.5F，R2=R3=2欧姆，L4=1H；US（t）=10+10cost,Is(t)=5+5cos2t,求b，d两点之间的电压U（t）。解：Matlab求解：clear,format compactw=eps,1,2;Us=10,10,0;Is=5,0,5;Z1=1./(0

21、.5*w*j);Z4=1*w*j;Z2=2,2,2;Z3=2,2,2;Uoc=(Z2./(Z1+Z2)-Z4./(Z3+Z4).*Us;Zeq=Z3.*Z4./(Z3+Z4)+Z1.*Z2./(Z1+Z2);U=Is.*Zeq+Uoc;disp(' w Um phi ')disp(w',abs(U'),angle(U')*180/pi)仿真结果：w Um phi 0.0000 10.0000 01.0000 3.1623 -18.43492.0000 7.0711 -8.13011. 含受控源的电路，戴维南定理如图所示电路，设Z1=-J250 Ohms，

22、Z2=250欧姆，Is=20 A，求负载ZL获得最大功率时的阻抗值及其吸收功率。解：Matlab求解：clear,format compactZ1=-j*250;Z2=250;ki=0.5;Is=2;a11=1/Z1+1/Z2;a12=-1/Z2;a13=0;a21=-1/Z2;a22=1/Z2;a23=-ki;a31=1/Z1;a32=0;a33=-1;A=a11,a12,a13;a21,a22,a23;a31,a32,a33;B=1,0;0,1;0,0;X0=AB*Is;0;Uoc=X0(2),X1=AB*0;1;Zeq=X1(2),PLmax=(abs(Uoc)2/4/real(Zeq)

23、仿真结果：Uoc = 5.0000e+002 -1.0000e+003iZeq = 5.0000e+002 -5.0000e+002iPLmax = 625四、实验内容 1.如图所示电路，设R1=2欧姆，R2=3欧姆，R3=4欧姆，jXL=j2，-jXC1=-j3，-jXC2=-j5，Us1=80 V,Us2=60 V, Us3=80, Us4=150, 求各支路的电流相量和电压相量。解：Matlab求解：clear,format compactR1=2;R2=3;R3=4;ZL=2*j;ZC1=-3*j;ZC2=-5*j;US1=8;US2=6;US3=8;US4=15;Y1=1/R1+1/

24、ZL;Y2=1/ZC1+1/R2;Y3=1/R3+1/ZC2;a11=1/Y1;a12=1/Y2;a13=1/Y3;a21=0;a22=-1;a23=1;a31=-1;a32=1;a33=0;b1=0;b2=US2/R2-US3/R3-US4/ZC2;b3=-US1/ZL-US2/R2;A=a11,a12,a13;a21,a22,a23;a31,a32,a33;B=b1;b2;b3;I=inv(A)*B;I1=I(1),I2=I(2),I3=I(3),ua=I1/Y1,ub=I3./(-Y3),I1R=ua/R1,I1L=(US1-ua)./ZL,I2R=(US2-ua+ub)/R2,I2C=

25、(ua-ub)./ZC1,I3R=(US3-ub)/R3,I3C=(US4-ub)./ZC2仿真结果：I1 = 1.2250 - 2.4982iI2 = -0.7750 + 1.5018iI3 = -0.7750 - 1.4982iua = 3.7232 - 1.2732iub = 4.8135 + 2.1420iI1R = 1.8616 - 0.6366iI1L = 0.6366 - 2.1384iI2R = 2.3634 + 1.1384iI2C = 1.1384 - 0.3634iI3R = 0.7966 - 0.5355iI3C = 0.4284 + 2.0373i2.含互感的电路：复

26、功率如图电路，已知R1=4 欧姆，R2=R3=2欧姆，XL1=10欧姆，XL2=8欧姆，XM=4欧姆，Xc=8欧姆，Us=100 V, Is=100 A。求电压源、电流源发出的复功率。解：Matlab求解：clear,format compactR1=4;R2=2;R3=2;XL1=10;XL2=8;XM=4;XC=8;US=10;IS=10;Y1=1/R1+1/(-j*XC);Y2=1/(j*(XL1-XM);Y3=1/(j*XM);Y4=1/(j*(XL2-XM)+R2);Y5=1/R3;a11=1;a12=-1;a13=0;a14=0;a15=0;a21=0;a22=0;a23=0;a2

27、4=1;a25=-1;a31=0;a32=1;a33=-1;a34=-1;a35=0;a41=1/Y1;a42=1/Y2;a43=1/Y3;a44=0;a45=0;a51=0;a52=0;a53=-1/Y3;a54=1/Y4;a55=1/Y5;A=a11,a12,a13,a14,a15;a21,a22,a23,a24,a25;a31,a32,a33,a34,a35;a41,a42,a43,a44,a45;a51,a52,a53,a54,a55;B=-US/R1;-IS;0;0;0;I=inv(A)*B;I1=I(1);I2=I(2);I3=I(3);I4=I(4);I5=I(5);ua=-I1

28、/Y1;ub=I3/Y3;uc=I5/Y5;Ii=US/R1+ua/R1;Pus=US*IiPis=uc*IS仿真结果： Pus = 54.0488 - 9.3830iPis = 1.7506e+002 +3.2391e+001i2.正弦稳态电路，求未知参数如图所示电路，已知Us=100V，I1=100mA，电路吸收的功率P=6W，XL1=1250欧姆，Xc=750欧姆，电路呈感性，求R3及XL3.3.、正弦稳态电路，利用模值求解如图所示电路，已知IR=10A，Xc=10欧姆，并且U1=U2=100V，求XL。Matlab求解：clearU2=200;IR=10;R=U2/IR;XC=10;U

29、=200*exp(-150j*pi/180);200*exp(-30j*pi/180);I=(U-200)./(-j*XC);X=200./(I-10);XL=imag(X)仿真结果：XL =5.3590 74.6410实验四、交流分析和网络函数一、实验目的1. 学习交流电路的分析方法。2. 学习交流电路的MATLAB分析方法。二、预习要求1.复习有关网络函数的相关概念三、实验示例1.在图里，如果R1=20欧姆，R2=100欧姆，R3=50欧姆，并且L1=4H，L2=8H以及C1=250微法，求v3(t)，其中w=10 rad/s. 解：Matlab求解：Y=0.05-0.0225*j 0.0

30、25*j -0.0025*j;0.025*j 0.01-0.0375*j 0.0125*j;-0.0025*j 0.0125*j 0.02-0.01*j; c1=0.4*exp(pi*15*j/180); i=c1;0;0;V=inv(Y)*i; v3_abs=abs(V(3); v3_ang=angle(V(3)*180/pi; fprintf('voltage v3, magnitude: %f n voltage v3, angle in degree:%f',v3_abs,v3_ang)仿真结果：voltage v3, magnitude: 1.850409 volta

31、ge v3, angle in degree:-72.453299四、实验内容 1.电路显示如图，求电流i1(t)和电压VC（t）.将本题中电路图改成频域模型，产生的电路如下图所示，电阻为欧姆。Matlab求解：Y=1 1 -1;6-5*j 0 4-2.5*j;6-5*j -10-8*j 0;c2=5;c3=2*exp(pi*75*j/180);v=0;c2;c3;i=inv(Y)*v;it_abs=abs(i(3);it_ang=angle(i(3)*180/pi;Vc_abs=abs(i(1)*-10*j);Vc_ang=angle(i(1)*-10*j)*180/pi;fprintf(&

32、#39;voltage it,magnitude: %f n voltage it,angle in degree: %f ',it_abs,it_ang)voltage it,magnitude: 0.387710 voltage it,angle in degree: 15.019255 >> fprintf('voltage Vc,magnitude: %f n voltage Vc,angle in degree: %f ',Vc_abs,Vc_ang)voltage Vc,magnitude: 4.218263 voltage Vc,angle i

33、n degree: -40.861691 >>3. 在图里，显示一个不平衡wye-wye系统，求相电压Van，Vbn，Vcn.Matlab求解Y=6+13*j 0 0;0 4+6*j 0;0 0 6-12.5*j;c1=110;c2=110*exp(pi*(-120)*j/180);c3=110*exp(pi*120*j/180);v=c1;c2;c3;i=inv(Y)*v;Van_abs=abs(i(1)*(5+12*j);Van_ang=angle(i(1)*(5+12*j)*180/pi;Vbn_abs=abs(i(2)*(3+4*j);Vbn_ang=angle(i(2)*

34、(3+4*j)*180/pi;Vcn_abs=abs(i(3)*(5-12*j);Vcn_ang=angle(i(3)*(5-12*j)*180/pi;Y=6+13*j 0 0;0 4+2*j 0;0 0 6-12.5*j; c1=110;c2=110*exp(pi*(-120)*j/180);c3=110*exp(pi*120*j/180);i=inv(Y)*v; Van_abs=abs(i(1)*(5+12*j);Van_ang=angle(i(1)*(5+12*j)*180/pi;Vbn_abs=abs(i(2)*(3+4*j); Vbn_ang=angle(i(2)*(3+4*j)*1

35、80/pi;Vcn_abs=abs(i(3)*(5-12*j);Vcn_ang=angle(i(3)*(5-12*j)*180/pi; fprintf('voltage Van,magnitude: %f n voltage Van,angle in degree: %f ',Van_abs,Van_ang);voltage Van,magnitude: 99.875532 voltage Van,angle in degree: 2.155276 >> fprintf('voltage Vbn,magnitude: %f n voltage Vbn,ang

36、le in degree: %f ',Vbn_abs,Vbn_ang);仿真结果如下：voltage Vbn,magnitude: 122.983739 voltage Vbn,angle in degree: -93.434949 fprintf('voltage Vcn,magnitude: %f n voltage Vcn,angle in degree: %f ',Vcn_abs,Vcn_ang);voltage Vcn,magnitude: 103.134238 voltage Vcn,angle in degree: 116.978859 实验五、动态电路一

37、、实验目的1. 学习动态电路的分析方法.2. 学习动态电路的MATLAB的计算方法。二、预习要求3. 复习一阶、二阶动态电路的基本概念4. 了解MATLAB动态电路的计算方法。三、实验示例1.一阶动态电路，三要素公式电路如图所示，已知R1=3 Ohms, R2=12 Ohms, R3=6 Ohms, C=1F; Us=18V, is=3A, 在t<0时， 开关s位于“1”，电路已处于稳态。（1） t=0时。开关s闭合到“2”，求Uc（t），iR2（t），并画出波形。（2） 若经10s，开关s又复位到“1”，求Uc（t），iR2（t），并画出波形。解：Matlab求解：clear allR

38、1=3;us=18;is=3;R2=12;R3=6;C=1;uc0=-12;ir20=uc0/R2;ir30=uc0/R3;ic0=is-ir20-ir30;ir2f=is*R3/(R2+R3);ir3f=is*R2/(R2+R3);ucf=ir2f*R2;icf=0;t=-2-eps,0-eps,0+eps,0:9,10-eps,10+eps,11:20figure(1),plot(t),griduc(1:3)=-12;ir2(1:3)=3;T=R2*R3/(R2+R3)*C;uc(4:14)=ucf+(uc0-ucf)*exp(-t(4:14)/T);ir2(4:14)=ir2f+(ir2

39、0-ir2f)*exp(-t(4:14)/T);uc(15)=uc(14);ir2(15)=is;ucf2=-12;ir2f=is;T2=R1*R3/(R1+R3)*C;uc(15:25)=ucf2+(uc(15)-ucf2)*exp(-(t(15:25)-t(15)/T2);ir2(15:25)=is;figure(2)subplot(2,1,1);h1=plot(t,uc);grid,set(h1,'linewidth',2)subplot(2,1,2);h2=plot(t,ir2);grid,set(h2,'linewidth',2)仿真结果：1. 二阶过

40、阻尼电路的零输入响应下图是典型的二阶动态电路，其零输入响应有过阻尼、临界阻尼。欠阻尼三种情况。本例讨论过阻尼的情况，如已知L=0.5H，C=0.02F，R=12.5欧姆，初始值UC（0）=1v，iL（0）=0，求t=0时的Uc（t）和IL（t）的零输入响应，并画出波形。解：Matlab求解：方法一：clear,format compactL=0.5;R=12.5;C=0.02;uc0=1;iL0=0;alpha=R/2/L;wn=sqrt(1/(L*C);p1=-alpha+sqrt(alpha2-wn2);p2=-alpha-sqrt(alpha2-wn2);dt=0.01;t=0:dt:1

41、;uc1=(p2*uc0-iL0/C)/(p2-p1)*exp(p1*t);uc2=-(p1*uc0-iL0/C)/(p2-p1)*exp(p2*t);iL1=p1*C*(p2*uc0-iL0/C)/(p2-p1)*exp(p1*t);iL2=-p2*C*(p1*uc0-iL0/C)/(p2-p1)*exp(p2*t);uc=uc1+uc2;iL=iL1+iL2;subplot(2,1,1),plot(t,uc),gridsubplot(2,1,2),plot(t,iL),grid仿真结果：方法二：Matlab求解：clear,format compactL=0.5;R=12.5;C=0.02

42、;uc0=1;iL0=0;alpha=R/2/L;wn=sqrt(1/(L*C);p1=-alpha+sqrt(alpha2-wn2);p2=-alpha-sqrt(alpha2-wn2);dt=0.01;t=0:dt:1;num=uc0,R/L*uc0+iL0/C;den=1,R/L,1/L/C;r,p,k=residue(num,den);ucn=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);iLn=C*diff(ucn)/dt;figure(1),subplot(2,1,1),plot(t,ucn),gridsubplot(2,1,2)plot(t(1:end-1)

43、,iLn),grid仿真结果：四、实验内容1.正弦激励的一阶电路如图所示，已知R2=2欧姆，C=0.5F，电容初始电压Uc（0）=4V，激励的正弦电压Us（t）=umcoswt，其中w=2 rad/s。当t=0时，开关s闭合，求电容电压的全响应，区分其暂态响应的稳态响应，并画出波形。解：Matlab求解：uc0=4;w=2;R=2;C=1;Zc=1/(j*w*C);dt=0.1;t=0:dt:10;us=6*cos(w*t);%È¡Um=6T=R*C;ucf=us*Zc/(Zc+R);uc1=uc0*exp(-t/T);figure(1);subplot(3,1,1);h1

44、=plot(t,ucf);grid,set(h1,'linewidth',2)subplot(3,1,2);h2=plot(t,uc1);grid,set(h2,'linewidth',2);uc=ucf+uc1;subplot(3,1,3);h3=plot(t,uc);grid,set(h3,'linewidth',2)仿真结果：2.二阶欠阻尼电路的零输入响应如图所示的二阶电路，如L=0.5H，C=0.02F，初始值uc（0）=1v，IL=0，试研究及分别为1欧姆、2欧姆、，10欧姆时，uc（t）和IL（t）的零输入响应，并画出波形。解：Mat

45、lab求解：1、R=1欧姆clear,format compactL=0.5;R=1;C=0.02;uc0=1;iL0=0;alpha=R/2/L;wn=sqrt(1/(L*C);p1=-alpha+sqrt(alpha2-wn2);p2=-alpha-sqrt(alpha2-wn2);dt=0.01;t=0:dt:1;num=uc0,R/L*uc0+iL0/C;den=1,R/L,1/L/C;r,p,k=residue(num,den);ucn=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);iLn=C*diff(ucn)/dt;figure(1),subplot(2,1

46、,1),plot(t,ucn),gridsubplot(2,1,2)plot(t(1:end-1),iLn),grid仿真结果如下：2、R=2欧姆代码同上，仿真结果如下图3、 R=34、 R=45、 R=56、 R=67、 R=78、 R=89、 R=910、R=10五、实验总结1、熟悉并了解了MATLAB暂态电路的计算方法2、通过MATLAB的仿真更加深刻理解了一阶二阶动态电路的基本工作状态。实验六、频率响应一、实验目的1.学习有关频率响应的相关概念2.学习的频率计算二、预习要求1.频率响应函数（即正弦稳态网络函数）定义为输出相量与激励相量之比，即（）（）（j（w）=Y/F.复习相关内容2.

47、用MATLAB中的abs（H）和angle（H）语句可直接计算幅频响应和相频响应，而且其图像的频率坐标（横坐标）可以是线性的，也可以是半对数的，这给计算和绘制特性和相频特性带来很大方便，了解相关函数的用法。三实验示例1.一阶低通电路的频率响应图示电路是一阶低通电路，若以为响应，求频率响应函数，画出其幅频特性（）和相频特性（）。解：MATLAB求解：clear,format compactww=0:0.2:4;H=1./(1+j*ww);figure(1)subplot(2,1,1),plot(ww,abs(H),grid,xlabel('ww'),ylabel('ang

48、le(H)');subplot(2,1,2),plot(ww,angle(H),grid,xlabel('ww'),ylabel('angle(H)');figure(2)subplot(2,1,1),semilogx(ww,20*log10(abs(H),grid,xlabel('ww'),ylabel('abs（H）');subplot(2,1,2),semilogx(ww,angle(H),grid,xlabel('ww'); ylabel('angle(H)')结果为A：线性频率特性

49、B：对数频率特性频率响应：二阶低通电路二阶低通函数的典型形式为（）（）令H0=1，画出Q=1/3,1/2,1/2,1,2,5的幅频相频响应，当Q=1/2时，成为最平幅度特性，即在通带内其幅频特性最为平坦。Matlab求解：clear,format compactfor Q=1/3,1/2,1/sqrt(2),1,2,5ww=logspace(-1,1,50);H=1./(1+j*ww/Q+(j*ww).2);figure(1)subplot(2,1,1),plot(ww,abs(H),hold onsubplot(2,1,2),plot(ww,angle(H),hold onfigure(2)subplot(2,1,1),semilogx(ww,20*log10(abs(H),hold onsubplot(2,1