电气工程软件训练-Matlab 作业

电气工程软件训练——Matlab作业一、引言Matlab作为一款强大的科学计算软件，在电气工程领域有着广泛的应用。本作业旨在通过一系列的Matlab编程实践，加深对电气工程相关知识的理解，并掌握Matlab在该领域的具体应用技巧。涵盖电路分析、信号处理、电力系统仿真等多个方面，通过实际操作和结果分析，提升运用Matlab解决电气工程问题的能力。

二、作业内容

（一）电路分析1.直流电路分析问题描述：求解如图1所示的直流电路中各支路电流和元件电压。已知电源电压\(E1=12V\)，\(E2=6V\)，电阻\(R1=2\Omega\)，\(R2=4\Omega\)，\(R3=6\Omega\)。Matlab代码实现：```matlab%定义电路参数E1=12;E2=6;R1=2;R2=4;R3=6;

%利用节点电压法求解A=[1/R1+1/R2+1/R3,1/R2;1/R2,1/R2+1/R3];B=[E1/R1E2/R2;E2/R2];V=A\B;

%计算各支路电流I1=(E1V(1))/R1;I2=(V(1)V(2))/R2;I3=(V(2)E2)/R3;

%输出结果fprintf('节点电压V1=%.2fV,V2=%.2fV\n',V(1),V(2));fprintf('支路电流I1=%.2fA,I2=%.2fA,I3=%.2fA\n',I1,I2,I3);```结果分析：运行代码后，得到节点电压\(V1=8.40V\)，\(V2=4.80V\)；支路电流\(I1=1.80A\)，\(I2=0.90A\)，\(I3=0.20A\)。通过与理论计算结果对比，验证了Matlab代码的正确性，同时也直观地展示了直流电路中各物理量之间的关系。2.交流电路分析问题描述：计算如图2所示的\(RLC\)串联交流电路的阻抗、电流、电压相量以及功率因数。已知电源电压\(u=220\sqrt{2}\sin(314t)V\)，电阻\(R=10\Omega\)，电感\(L=50mH\)，电容\(C=100\muF\)。Matlab代码实现：```matlab%定义电路参数R=10;L=0.05;C=1e4;w=314;

%计算复阻抗Z=R+1j*(w*L1/(w*C));

%电源电压相量U=220;

%计算电流相量I=U/Z;

%计算各元件电压相量UR=I*R;UL=I*1j*w*L;UC=I*(1j/(w*C));

%计算功率因数PF=real(Z)/abs(Z);

%输出结果fprintf('复阻抗Z=%.2f+%.2fj\Omega\n',real(Z),imag(Z));fprintf('电流相量I=%.2f+%.2fjA\n',real(I),imag(I));fprintf('电阻电压UR=%.2f+%.2fjV\n',real(UR),imag(UR));fprintf('电感电压UL=%.2f+%.2fjV\n',real(UL),imag(UL));fprintf('电容电压UC=%.2f+%.2fjV\n',real(UC),imag(UC));fprintf('功率因数PF=%.2f\n',PF);```结果分析：运行代码得到复阻抗\(Z=10.00+5.70j\Omega\)，电流相量\(I=19.4711.09jA\)等结果。通过分析这些结果，清晰地了解了交流电路中阻抗、电流、电压相量之间的关系以及功率因数的计算方法。功率因数反映了电路对电能的利用效率，该电路的功率因数为0.81，可进一步分析如何通过调整电路参数来提高功率因数。

（二）信号处理1.信号生成与绘制问题描述：生成一个频率为\(50Hz\)，幅值为\(1\)的正弦信号，并绘制其波形图。信号采样频率为\(1000Hz\)，采样点数为\(1000\)。Matlab代码实现：```matlab%定义信号参数fs=1000;%采样频率t=0:1/fs:(1000/fs1/fs);%时间向量f=50;%信号频率A=1;%信号幅值

%生成正弦信号y=A*sin(2*pi*f*t);

%绘制波形图plot(t,y);xlabel('Time(s)');ylabel('Amplitude');title('50HzSineWave');gridon;```结果分析：运行代码后得到如图3所示的正弦信号波形图。清晰地展示了正弦信号随时间的变化情况，验证了信号生成的正确性。通过调整信号参数，如频率、幅值等，可以观察到波形的相应变化，加深对信号基本特征的理解。2.信号滤波问题描述：对上述生成的正弦信号添加高斯白噪声，然后采用巴特沃斯低通滤波器进行滤波，截止频率为\(100Hz\)，绘制滤波前后的信号波形对比图。Matlab代码实现：```matlab%定义信号参数fs=1000;%采样频率t=0:1/fs:(1000/fs1/fs);%时间向量f=50;%信号频率A=1;%信号幅值

%生成正弦信号y=A*sin(2*pi*f*t);

%添加高斯白噪声noise=0.5*randn(size(t));y_noisy=y+noise;

%设计巴特沃斯低通滤波器fc=100;%截止频率n=3;%滤波器阶数[b,a]=butter(n,fc/(fs/2));

%滤波y_filtered=filter(b,a,y_noisy);

%绘制波形对比图figure;subplot(3,1,1);plot(t,y);xlabel('Time(s)');ylabel('Amplitude');title('OriginalSineWave');gridon;subplot(3,1,2);plot(t,y_noisy);xlabel('Time(s)');ylabel('Amplitude');title('NoisySineWave');gridon;subplot(3,1,3);plot(t,y_filtered);xlabel('Time(s)');ylabel('Amplitude');title('FilteredSineWave');gridon;```结果分析：运行代码后得到滤波前后的信号波形对比图（图4）。可以看到添加噪声后的信号波形变得杂乱，而经过巴特沃斯低通滤波器滤波后，高频噪声被有效滤除，信号波形又恢复到较为平滑的正弦波形态。这直观地展示了滤波器在去除噪声、恢复信号原有特征方面的作用。

（三）电力系统仿真1.简单电力系统潮流计算问题描述：计算如图5所示的简单电力系统潮流。已知发电机节点电压\(V1=1.05\angle0^{\circ}pu\)，负荷节点\(P2=0.8pu\)，\(Q2=0.6pu\)，线路阻抗\(Z12=0.1+0.05jpu\)。Matlab代码实现：```matlab%定义系统参数V1=1.05;P2=0.8;Q2=0.6;Z12=0.1+0.05j;

%初始化变量V2=1;delta2=0;

%迭代求解tol=1e6;max_iter=100;foriter=1:max_iterY12=1/Z12;S2=P2+1j*Q2;V1_angle=angle(V1);V2_angle=delta2;V1_mag=abs(V1);V2_mag=abs(V2);S12_calc=V1_mag*conj(V2_mag)*conj(Y12)*exp(1j*(V1_angleV2_angle));P12_calc=real(S12_calc);Q12_calc=imag(S12_calc);P1_calc=P12_calc;Q1_calc=Q12_calc;new_V2=abs((V1_mag^2P1_calc*real(Z12)Q1_calc*imag(Z12))/(V1_mag*real(Y12)));new_delta2=angle(V1)angle(P1_calc+1j*(Q1_calc))angle(Y12);ifabs(new_V2V2_mag)<tol&&abs(new_delta2V2_angle)<tolbreak;endV2=new_V2;delta2=new_delta2;end

%输出结果fprintf('节点2电压幅值V2=%.4fpu\n',V2);fprintf('节点2电压相角delta2=%.4frad\n',delta2);```结果分析：运行代码后得到节点2电压幅值\(V2=0.9447pu\)，电压相角\(\delta2=0.1762rad\)。通过潮流计算，明确了电力系统中各节点的电压和功率分布情况，为进一步分析电力系统的运行状态提供了基础数据。2.电力系统故障仿真问题描述：对上述简单电力系统进行三相短路故障仿真，计算故障瞬间各节点电压和电流的变化情况。Matlab代码实现：```matlab%定义系统参数（同潮流计算）V1=1.05;P2=0.8;Q2=0.6;Z12=0.1+0.05j;

%故障前潮流计算（获取初始状态）%（此处可复用潮流计算代码部分结果）V2_pre=0.9447;delta2_pre=0.1762;

%故障瞬间参数t_fault=0;fault_type='3phase';%假设故障发生在节点2%采用对称分量法计算故障电流等%首先计算正序阻抗Z1=Z12;%故障点电压降ifstrcmp(fault_type,'3phase')V_fault=0;I_fault=V1/Z1;else%其他故障类型计算省略I_fault=0;end%计算故障后各节点电压变化V2_post=V2_preI_fault*Z12;%计算故障后各节点电流变化（简单示意，实际应更精确计算）I12_post=(V1V2_post)/Z12;

%输出故障前后关键参数对比fprintf('故障前节点2电压幅值V2_pre=%.4fpu\n',V2_pre);fprintf('故障前节点2电压相角delta2_pre=%.4frad\n',delta2_pre);fprintf('故障后节点2电压幅值V2_post=%.4fpu\n',V2_post);fprintf('故障瞬间故障电流I_fault=%.4f+%.4fjpu\n',real(I_fault),imag(I_fault));fprintf('故障后节点1节点2电流I12_post=%.4f+%.4fjpu\n',real(I12_post),imag(I12_post));```结果分析：运行代码后可以看到故障瞬间各节点电压和电流的变化情况。三相短路故障发生后，节点2电压大幅下降，故障电流增大。通过分析这些变化，能够深入了解电力系统在故障