短路电流计算

华南理工大学广州学院电气工程系电力系统分析程序报告专业课程设计报告题目：短路电流在matlab中的实现系别电气工程系专业班级***************学生姓名fish_barbar指导教师吴彦祖提交日期***年*月*日华南理工大学广州学院电气工程系电力系统分析程序报告1一、电力系统潮流计算机计算的意义和目的潮流计算是电力系统分析中一种基本的计算，它的任务是对给定的运行条件确定系统的运行状态。采用计算机计算，可以减轻计算人员的工作量，提高计算准确度和计算速度。二、三机九节点的等效电网图279836415(3)f图1-1三机九节点的等效电网图华南理工大学广州学院电气工程系电力系统分析程序报告2三、程序计算的流程图图1-2程序计算的流程图华南理工大学广州学院电气工程系电力系统分析程序报告3四、计算机编程主函数Sbase_MVA=100.fid=fopen('Nodedata.txt');N=textscan(fid,'%s%u%d%f%f%f%f%f%f');fclose(fid);busnumber=size(N{1},1);fori=1:busnumberBus(i).name=N{1}(i);Bus(i).type=N{2}(i);Bus(i).no=i;Bus(i).Base_KV=N{3}(i);Bus(i).PG=N{4}(i);Bus(i).QG=N{5}(i);Bus(i).PL=N{6}(i);Bus(i).QL=N{7}(i);Bus(i).pb=N{8}(i);Bus(i).V=1.0;Bus(i).angle=0;endfid=fopen('Aclinedata.txt');A=textscan(fid,'%s%s%f%f%f%f');fclose(fid);aclinenumber=size(A{1},1);fori=1:aclinenumberAcline(i).fbname=A{1}(i);Acline(i).tbname=A{2}(i);Acline(i).Base_KV=A{3}(i);华南理工大学广州学院电气工程系电力系统分析程序报告4Acline(i).R=A{4}(i);Acline(i).X=A{5}(i);Acline(i).hB=A{6}(i);fork=1:busnumberifstrcmp(Acline(i).fbname,Bus(k).name)Acline(i).fbno=Bus(k).no;endifstrcmp(Acline(i).tbname,Bus(k).name)Acline(i).tbno=Bus(k).no;endendendfid=fopen('Transdata.txt');T=textscan(fid,'%s%f%f%s%f%f%f%f');fclose(fid);tansnumber=size(T{1},1);fori=1:tansnumberTrans(i).fbname=T{1}(i);Trans(i).fbBase_KV=T{2}(i);Trans(i).fbrated_KV=T{3}(i);Trans(i).tbname=T{4}(i);Trans(i).tbBase_KV=T{5}(i);Trans(i).tbrated_KV=T{6}(i);Trans(i).R=T{7}(i);Trans(i).X=T{8}(i);fork=1:busnumberifstrcmp(Trans(i).fbname,Bus(k).name)华南理工大学广州学院电气工程系电力系统分析程序报告5Trans(i).fbno=Bus(k).no;endifstrcmp(Trans(i).tbname,Bus(k).name)Trans(i).tbno=Bus(k).no;endendTrans(i).k=Trans(i).tbrated_KV*Trans(i).fbBase_KV/Trans(i).fbrated_KV/Trans(i).tbBase_KV;tempx=Trans(i).fbrated_KV^2/Trans(i).fbBase_KV^2;Trans(i).X=tempx*Trans(i).X;Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2)%forY=G+jBmatrix[G,B,B2]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber);%B:=B';B2:=Bdlmwrite('Gmatrix.txt',G,'delimiter','\t','precision',6);dlmwrite('Bmatrix.txt',B,'delimiter','\t','precision',6);%G%B%B2YN=G+j*B;%不包括发电机导纳和负荷导纳的节点导纳矩阵%Ychushi=YN;华南理工大学广州学院电气工程系电力系统分析程序报告6%pause[JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber);JP;iJP=-inv(JP);JQ;iJQ=-inv(JQ);%pause%maxiteration=0fori=1:busnumberNodeV(i)=Bus(i).V;Nodea(i)=Bus(i).angle;VX(i)=Bus(i).V*cos(Bus(i).angle);VY(i)=Bus(i).V*sin(Bus(i).angle);dQGQL(i)=Bus(i).QG-Bus(i).QL;dPGPL(i)=Bus(i).PG-Bus(i).PL;endNodeV=NodeV';Nodea=Nodea';%VX=VX'%VY=VY'dQGQL=dQGQL';dPGPL=dPGPL';%pause;%fornointer=1:10maxdP=1.;华南理工大学广州学院电气工程系电力系统分析程序报告7maxdQ=1.;epsilon=0.000001;noiteration=0;while(maxdPepsilon)&(maxdPepsilon)[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,busnumber);deltaP;deltaQ;%maxdP%maxdQda=iJP*deltaP;dV=iJQ*deltaQ;Nodea=Nodea+da;NodeV=NodeV+dV;noiteration=noiteration+1;%pauseifnoiteration20breakendendfori=1:busnumberBus(i).V=NodeV(i);NodeV(i)=NodeV(i)*Bus(i).Base_KV;Bus(i).angle=Nodea(i);Nodea(i)=Nodea(i)*180/pi;end%noiterationNodea=Nodea'华南理工大学广州学院电气工程系电力系统分析程序报告8NodeV=NodeV'%潮流计算节点电压模及相角y=1.0/(j*0.3);%pausefori=1:3YN(i,i)=YN(i,i)+y;endYchushi=YN;fori=4:9y=(Bus(i).PL-j*Bus(i).QL)/Bus(i).V^2;YN(i,i)=YN(i,i)+y;endYN;uM=[0.00.00.01.00.00.00.00.00.0]';fori=1:9fork=1:9ifYchushi(i,k)~=0ifi~=kZ(i,k)=-1/(Ychushi(i,k));endendend华南理工大学广州学院电气工程系电力系统分析程序报告9endZ4=inv(YN)*uM;If=(Bus(4).V*cos(Bus(4).angle)+j*Bus(4).V*sin(Bus(4).angle))/Z4(4);IB230=100.0/sqrt(3)/230.0Ifymzh=abs(If*IB230)fori=1:busnumberV0(i)=Bus(i).V*cos(Bus(i).angle)+j*Bus(i).V*sin(Bus(i).angle);Vf(i)=V0(i)-Z4(i)*If;endV0Vffori=1:9Vf1(i)=abs(Vf(i));endVf1%由V0和Vf计算各支路电流程序补上ZIline=zeros(9,9);fori=1:aclinenumber华南理工大学广州学院电气工程系电力系统分析程序报告10fork=1:busnumberifstrcmp(Acline(i).fbname,Bus(k).name)F1=Bus(k).no;fork2=1:busnumberifstrcmp(Acline(i).tbname,Bus(k2).name)F2=Bus(k2).no;Iline(F2,F1)=(Vf(F2)-Vf(F1))/(Z(F1,F2));endendendendendIline1=Iline;It=zeros(9,9);K=(230/18)/(230/18);fori=1:tansnumberfork2=1:busnumberifstrcmp(Trans(i).fbname,Bus(k2).name)F1=Bus(k2).no;forkk=1:busnumberifstrcmp(Trans(i).tbname,Bus(kk).name)F2=Bus(kk).no;It(F1,F2)=(K*Vf(F1)-Vf(F2))/(K*Z(F2,F1));endendendendend华南理工大学广州学院电气工程系电力系统分析程序报告11It1=It;Iline1Iline=Iline*IB230;fori=1:9forh=1:9shibu=real(Iline(i,h));xubu=imag(Iline(i,h));ifshibu=0&&xubu=0Iline(i,h)=-abs(Iline(i,h));elseIline(i,h)=abs(Iline(i,h));endendendIlineItIt=It*IB230;It=abs(It)I2=It1(1,4)+Iline1(5,4)+Iline1(6,4);dd2=10^(-7);华南理工大学广州学院电气工程系电力系统分析程序报告12dd1=abs(If-I2);ifdd1dd2endclear子函数FormYmatrixfunction[G,B,X]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber)Y=zeros(busnumber);X=zeros(busnumber);fori=1:busnumberY(i,i)=Y(i,i)+Bus(i).pb*j;endfori=1:aclinenumberf=Acline(i).fbno;t=Acline(i).tbno;Y(f,f)=Y(f,f)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(t,t)=Y(t,t)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(f,t)=Y(f,t)-1/(Acline(i).R+Acline(i).X*j);Y(t,f)=Y(t,f)-1/(Acline(i).R+Acline(i).X*j);X(f,f)=X(f,f)-1/Acline(i).X;X(t,t)=X(t,t)-1/Acline(i).X;X(f,t)=1/Acline(i).X;X(t,f)=1/Acline(i).X;endfori=1:tansnumberf=Trans(i).fbno;华南理工大学广州学院电气工程系电