《电力系统分析P-Q分解法潮流计算》由会员分享，可在线阅读，更多相关《电力系统分析P-Q分解法潮流计算（11页珍藏版）》请在金锄头文库上搜索。
1、.专业课程设计报告P-Q分解法潮流计算系 别 电气工程系 专业班级 09级电气4班 学生XX 钟剑帆 学号 8指导教师 房大中 提交日期 2011年11月 12日. v.目录P-Q分解法潮流计算1一、原理分析1二、程序流程1三、设计内容23.1 程序设计23.2.1 输入数据73.2.2 输出数据83.2.3 结果分析11四、心得总结12. v.一、原理分析从潮流计算的基本方程出发，采用PQ分解法并通过建立矩阵的修正方程来依次迭代，逐步逼近真值来计算电网的电压和功率分布。二、程序流程三、设计内容3.1 程序设计主函数Sbase_MVA=100.fid=fopen(Nodedata.txt);N
2、=textscan(fid, %s %u %d %f %f %f %f %f %f)fclose(fid);busnumber=size(N1,1)for i=1:busnumber Bus(i).name=N1(i); Bus(i).type=N2(i); Bus(i).no=i; Bus(i).Base_KV=N3(i); Bus(i).PG=N4(i); Bus(i).QG=N5(i); Bus(i).PL=N6(i); Bus(i).QL=N7(i); Bus(i).pb=N8(i); Bus(i).V=1.0; Bus(i).angle=0;endfid=fopen(Aclineda
3、ta.txt);A=textscan(fid, %s %s %f %f %f %f)fclose(fid);aclinenumber=size(A1,1)for i=1:aclinenumber Acline(i).fbname=A1(i); Acline(i).tbname=A2(i); Acline(i).Base_KV=A3(i); Acline(i).R=A4(i); Acline(i).X=A5(i); Acline(i).hB=A6(i); for k=1:busnumber if strcmp(Acline(i).fbname, Bus(k).name) Acline(i).fb
4、no=Bus(k).no; end if strcmp(Acline(i).tbname, Bus(k).name) Acline(i).tbno=Bus(k).no; end end endfid=fopen(Transdata.txt);T=textscan(fid, %s %f %f %s %f %f %f %f)fclose(fid);tansnumber=size(T1,1)for i=1:tansnumber Trans(i).fbname=T1(i); Trans(i).fbBase_KV=T2(i); Trans(i).fbrated_KV=T3(i); Trans(i).tb
5、name=T4(i); Trans(i).tbBase_KV=T5(i); Trans(i).tbrated_KV=T6(i); Trans(i).R=T7(i); Trans(i).X=T8(i); for k=1:busnumber if strcmp(Trans(i).fbname, Bus(k).name) Trans(i).fbno=Bus(k).no; end if strcmp(Trans(i).tbname, Bus(k).name) Trans(i).tbno=Bus(k).no; end end Trans(i).k=Trans(i).tbrated_KV*Trans(i)
6、.fbBase_KV/Trans(i).fbrated_KV/Trans(i).tbBase_KV; tempx=Trans(i).fbrated_KV2/Trans(i).fbBase_KV2; Trans(i).X=tempx*Trans(i).X; Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2)% for Y=G+ matrixG,B,B2=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber); %B:=B;B2:=Bdlmwrite(Gmatrix.tx
7、t, G, delimiter, t,precision, 6);dlmwrite(Bmatrix.txt, B, delimiter, t,precision, 6);GBB2pauseJP,JQ=FormJPQmatrix(Bus,B,B2,busnumber);JPiJP=-inv(JP)JQiJQ=-inv(JQ)pause%maxiteration=0for i=1:busnumber NodeV(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
8、).angle); dQGQL(i)=Bus(i).QG-Bus(i).QL; dPGPL(i)=Bus(i).PG-Bus(i).PL; endNodeV=NodeVNodea=Nodea%VX=VX%VY=VYdQGQL=dQGQLdPGPL=dPGPLpause%for nointer=1:10maxdP=1.;maxdQ=1.;epsilon=0.000001;noiteration=0;while (maxdPepsilon)&(maxdPepsilon)deltaP,deltaQ,maxdP,maxdQ=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPG
9、PL,B,G,busnumber);deltaP;deltaQ;maxdP;maxdQ;da=iJP*deltaP;dV=iJQ*deltaQ;Nodea=Nodea+da;NodeV=NodeV+dV;noiteration=noiteration+1; if noiteration20 break endendfor i=1:busnumber Bus(i).V=NodeV(i); NodeV(i)=NodeV(i)*Bus(i).Base_KV; Bus(i).angle=Nodea(i); Nodea(i)=Nodea(i)*180/pi;endnoiterationNodea=Nod
10、eaNodeV=NodeVClear子函数%生成G、B矩阵function G,B,X=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber)Y=zeros(busnumber);X=zeros(busnumber);for i=1:busnumberY(i,i)=Y(i,i)+Bus(i).pb*j;endfor i=1:aclinenumber f=Acline(i).fbno; t=Acline(i).tbno; Y(f,f)=Y(f,f)+Acline(i).hB*j+1/(Acline(i).R+Acline(i
11、).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; endfor i=1:tansnumber f=Trans(i).fbno; t
12、=Trans(i).tbno; Y(f,f)=Y(f,f)+1/(Trans(i).R+Trans(i).X*j); Y(t,t)=Y(t,t)+1/(Trans(i).R+Trans(i).X*j)/Trans(i).k2; Y(f,t)=Y(f,t)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k; Y(t,f)=Y(t,f)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k; X(f,f)=X(f,f)-1/Trans(i).X; X(t,t)=X(t,t)-1/Trans(i).X; X(f,t)=1/Trans(i).X; X(t,f)=1/Trans(i).X; endG=real(Y);B=imag(Y);end%生成JP、JQ矩阵function JP,JQ=FormJPQmatrix(Bus,B,B2,busnumber)JP=B;JQ=B2;for i=1:busnumber if Bus(i).type=1 for k=1:busnumber JQ(i,k)=0.; JQ(k,i)=0.; JP(i,k)=0.; JP(k,i)=0.; end JQ(i,i)=1.; JP(i,i)=1.; end if Bus(i).type=3 for k=1:busnumber JQ(i,k)=0.;
