电力系统课程设计.doc

附图（一）：系统阻尼对暂态稳定性影响曲线(a)delta-t仿真曲线(b)w-t仿真曲线附图（二）：D=0时故障切除时间对系统暂态稳定性影响(a) delta-t曲线(b) w-t曲线附图（三）：D=0.02时故障切除时间对系统暂态稳定性影响(a) delta-t曲线（b）w-t曲线附图（四）：故障严重程度对系统暂态稳定性影响（a）delta-t曲线（b）w-t曲线附图（五）：系统载荷对系统暂态稳定性影响（a）delta-t曲线（b）w-t曲线附图（六）：部分程序清单主程序部分程序清单：%-系统的固定参数-Xd=0.238; Xd2=0.19; Tj=11.28; %发电机的参数XT1=0.13; XT2=0.108; %变压器的参数XL=0.293; XL0=5*XL; %线路的参数fN=50;WN=2*pi*fN;%- 系统的可变更参数-P0=1.2; S0=1.2+0.24j; U0=1.0; I0=1.2-0.24j;T=0.2; %故障持续时间DXL=0.01; %短路点距母线的电抗h=0.01; %迭代步长D=0.02;A=zeros(400,1);W=zeros(400,1);H=zeros(400,1); %-function0;Aclim%-系统阻尼对暂态稳定性影响-H1=figure;D=0.05;function1;plot(H,AA,r-);plot1;主程序里调用的M脚本文件（functiong0.m,function1.m,plot1.m）分别如下：functiong0.m%-设故障类型为单相接地短路-X0=5*DXL+XT1*(XL0+XT2)/(XT1+XL0+XT2); %零序电抗X2=DXL+(Xd2+XT1)*(XL+XT2)/(Xd2+XT1+XL+XT2);%负序电抗DDXL=X0+X2; %等效接地电抗XI=Xd+XT1+XL+XT2; %故障前电抗XII=XI+(Xd+XT1)*(XL+XT2)/(DXL+DDXL); %故障电抗XIII=XI; %故障后电抗E=U0+I0*XI*1.0j; %发电机电势%-临界切除角的计算-absE=abs(E);A0=angle(E);P1=absE*U0/XI;P2=absE*U0/XII;P3=P1;Ab=pi-asin(P0/P3);Aclim=acos(P0*(Ab-A0)+P3*cos(Ab)-P2*cos(A0)/(P3-P2);function1.m%-改进欧拉法的迭代过程-A(1)=A0; W(1)=WN; H(1)=0; A2=0; W2=0;k=1; P4=P2; h1=h;for n=1:399 A2=A(k)+h*(W(k)-WN); W2=W(k)+h*(WN2/(Tj*W(k)*(P0-D*(W(k)-WN)-P4*sin(A(k); A(k+1)=A(k)+0.5*h*(W(k)-WN)+(W2-WN); W(k+1)=W(k)+0.5*h*(WN2/(Tj*W(k)*(P0-D*(W(k)-WN)-P4*sin(A(k)+(WN2/(Tj*W2)*(P0-D*(W2-WN)-P4*sin(A2); H(k+1)=k*h; k=k+1; if (h*kT)&(h*(k-1)=T) P4=P3; endend AA=180*A/pi;plot1.maxis(0 4 -30 180);xlabel(t(s);ylabel(Delta(rad);title(R