结果为

Maximum power occur at 10000.00hms Maximum power dissipation is 0.0025Watts

2、在图示电路里(书本20页2-4)，当R1取0，2，4，6，10，18，24，42，90和186Ω时，求RL的电压UL,电流IL和RL消耗的功率。 解答：

使用matlab命令为： U=48; R=6; Rl=0;

IRl=U/(R+Rl);

fprintf('IRl= %8.1f A \\n',IRl) URl=IRl*Rl;

fprintf('URl= %8.1f V \\n',URl) PRl=IRl*URl;

fprintf('IRl= %8.1f W \\n',PRl)

输出结果：

IRl= 8.0 A URl= 0.0 V IRl= 0.0 W

IRl= 6.0 A URl= 12.0 V IRl= 72.0 W

IRl= 4.8 A URl= 19.2 V IRl= 92.2 W

11

IRl= 4.0 A URl= 24.0 V IRl= 96.0 W

IRl= 3.0 A URl= 30.0 V IRl= 90.0 W

IRl= 2.0 A URl= 36.0 V IRl= 72.0 W

IRl= 1.6 A URl= 38.4 V IRl= 61.4 W

IRl= 1.0 A URl= 42.0 V IRl= 42.0 W

IRl= 0.5 A URl= 45.0 V IRl= 22.5 W

IRl= 0.3 A URl= 46.5 V IRl= 11.6 W

电路分析图见下页

12

即

功率-电阻图为

13

实验三 正弦稳态

一、实验目的

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,I=Ic+Ir,U1=I*Z1,Us=I*Z disp('UcIrIc I U1 Us')

disp('·幅值'),disp(abs([Uc,Ir,Ic,I,U1,Us])) disp('相角'),disp(angle([Uc,Ir,Ic,U1,Us])*180/pi) ha=compass([Uc,Ir,Ic,U1,Us,Uc]); set(ha,'linewidth',3)

Ic =-2.5000 + 4.3301i Ir =1.7321 + 1.0000i I =-0.7679 + 5.3301i U1 =-15.9904 - 2.3038i Us =-7.3301 + 2.6962i UcIrIc I U1 Us 幅值

10.0000 2.0000 5.0000 5.3852 16.1555 7.8102 相角

30.0000 30.0000 120.0000 -171.8014 159.8056

14

90120 2060 15 1015030 51800210330240270300 2、正弦稳态电路，戴维南定理 如图3-3（书本22页），已知C1=0.5F,R2=R3=2Ω,L4=1H,Us(t)=10+10cost,is(t)=5 +5cos2t,求b，d两点之间的电压U(t) 使用matlab命令为：

clear,formatcompact

w=[eps,1,2];Us=[10,10,0];Is=[5,0,5]; Z1=1./(0.5*w*j);Z4=1*w*j; Z2=[2,2,2];Z3=[2,2,2];

Uoc=(Z2./(Z1+Z2)-Z4./(Z3+Z4)).*Us; Zep=Z3.*Z4./(Z3+Z4)+Z1.*Z2./(Z1+Z2); U=Is.*Zep+Uoc; disp('w Um phi')

disp([w',abs(U'),angle(U')*180/pi])

w Um phi

0.0000 10.0000 0 1.0000 3.1623 -18.4349 2.0000 7.0711 -8.1301

由此可以写出U（t）=10=3.1623cos（t-18.4394）+7.0711cos（2t-8.1301） 3、含受控源的电路：戴维南定理 如图3-4-1（书本23页），设Z1=-j250Ω,Z2=250Ω,Is=2∠0°,求负载Zl获得最大功率时的阻抗值及其吸收的功率。 使用matlab命令为

clear,formatcompact

Z1=-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;

15