自强、弘毅、求是、拓新
U2=200;IR=10;R=U2/IR;XC=10;
U=[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
实验总结
学习了正弦交流电路的分析方法，初步了解了MATLAB向量图的绘制，虽然并不能说完全掌握，但是基本有了一定的了解。
实验四交流分析和网络函数
实验目的
1学习交流电路的分析方法
2学习交流电路的MATLAB分析方法
实验内容
1求电流i1(t)和电压uc(t) Z=[10-j*7.5 -6+j*5; -6+j*5 10+j*3]; U2=2*exp(pi*75*j/180); U=[5; U2]; I=Z\\U; I1=I(1);
Uc=-j*10*(I(1)-I(2)); i1_abs=abs(I1);
i1_ang=angle(I1)*180/pi; uc_abs=abs(Uc);
uc_ang=angle(Uc)*180/pi;
fprintf('current I1,magnitude: %f \\n', i1_abs) fprintf('current I1,angle in degree:%f\\n',i1_ang)
11
自强、弘毅、求是、拓新
fprintf('voltage Uc,mangnitude: %f \\n',uc_abs) fprintf('voltage Uc,angle in degree:%f\\n',uc_ang) 结果
current I1,magnitude: 0.597294 current I1,angle in degree:13.266138 voltageUc,mangnitude: 3.202746 voltageUc,angle in degree:-28.980786
2显示一个不平衡wye-wye系统，求相电压VAN，VBN和VCN。 Ua=110;
Ub=110*exp(pi*(-120)*j/180); Uc=110*exp(pi*120*j/180); Za1=1-j;Zb1=1-j*2;Zc1=1-j*0.5; Za2=5-j*12;Zb2=3-j*4;Zc2=5-j*12; Uan=Ua*Za2/(Za1+Za2); Ubn=Ub*Zb2/(Zb1+Zb2); Ucn=Uc*Zc2/(Zc1+Zc2); disp('UanUbnUcn')
disp('幅值'),disp(abs([Uan,Ubn,Ucn]))
disp('相角'),disp(angle([Uan,Ubn,Ucn])*180/pi) ha=compass([Uan,Ubn,Ucn]); set(ha,'linewidth',3) 结果
UanUbnUcn 幅值
99.8755 76.2713 103.1342 相角
-2.1553 -116.8202 116.9789
UanUbnUcn 幅值
99.8755 76.2713 103.1342 相角
-2.1553 -116.8202 116.9789
实验总结
了解运用MATLAB分析交流电路。
12
自强、弘毅、求是、拓新
实验五动态电路
实验目的
1学习动态电路的分析方法
2学习动态电路的MATLAB计算方法
实验内容
1、激励的一阶电路
已知R=2欧姆，C=0.5F, 电容初始电压Uc(0+)=4V，激励的正弦电压Us(t)=Umcoswt，其中w=2rad/s。当t=0时，开关s闭合，求电容电压的全响应，区分其暂态响应与稳态响应，并画出波形
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=6 T=R*C;
ucf=us*Zc/(Zc+R); uc1=uc0*exp(-t/T); figure(1);
subplot(3,1,1); h1=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) 程序运行结果：
13
自强、弘毅、求是、拓新
0.50-0.542050-5012345678910012345678910012345678910
2、二阶欠阻尼电路的零输入响应
如图所示的二阶电路，如L=0.5H，C=0.02F。初始值uc（0）=1V，iL=0，试研究R分别为1Ω，2Ω，3Ω，…，10Ω时，uc（t）和iL(t)的零输入响应，并画出波形。 R=1
clear,format compact L=0.5;R=1;C=0.02; uc0=1;iL0=0;
alpha=R/2/L;wn=sqrt(1/(L*C)); p1=-alpha+sqrt(alpha^2-wn^2); p2=-alpha-sqrt(alpha^2-wn^2); 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),grid subplot(2,1,2)
plot(t(1:end-1),iLn),grid
14
自强、弘毅、求是、拓新
10.50-0.5-100.10.20.30.40.50.60.70.80.910.20.10-0.1-0.200.10.20.30.40.50.60.70.80.91
R=2
10.50-0.5-100.10.20.30.40.50.60.70.80.910.20.10-0.1-0.200.10.20.30.40.50.60.70.80.91
R=3
15