电力系统分析PQ分解法潮流计算.docVIP

电力系统分析P-Q分解法潮流计算 华南理工大学广州学院 专业课程设计报告 P-Q分解法潮流计算 系 别 电气工程系 专业班级 09级电气4班 学生姓名 钟剑帆 学 号 200930084778 指导教师 房大中 提交日期 2011年11月 12日 华南理工大学广州学院 目录 P-Q分解法潮流计算 .................................... 1 一、原理分析 ......................................... 1 二、程序流程 ......................................... 1 三、设计内容 ......................................... 2 3.1 程序设计 ...................................... 2 3.2.1 输入数据 .................................... 7 3.2.2 输出数据 .................................... 8 3.2.3 结果分析 ................................... 11 四、心得总结 ........................................ 12 华南理工大学广州学院 一、原理分析 从潮流计算的基本方程出发，采用PQ分解法并通过建立矩阵的修正方程来依次迭代，逐步逼近真值来计算电网的电压和功率分布。 二、程序流程 1 华南理工大学广州学院 三、设计 Acline(i).fbno=Bus(k).no; end if strcmp(Acline(i).tbname, Bus(k).name) Acline(i).tbno=Bus(k).no; end end 2 华南理工大学广州学院 end fid=fopen(‘Transdata.txt’); T=textscan(fid, ‘%s %f %f %s %f %f %f %f’) fclose(fid); tansnumber=size(T{1},1) for i=1:tansnumber Trans(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); 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).fbBase_KV/Trans(i).fbrated_KV/Trans(i).tbBase_KV; tempx=Trans(i).fbrated_KV/Trans(i).fbBase_KV; Trans(i).X=tempx*Trans(i).X; Trans(i).R=tempx*Trans(i).R; end %N=0 %Trans(1) %Trans(2) % for Y=G+jB matrix [G,B,B2]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber); %B:=B’;B2:=Bquot; dlmwrite(‘Gmatrix.txt’, G, ‘delimiter’, ‘\t’,’precision’, 6); dlmwrite(‘Bmatrix.txt’, B, ‘delimiter’, ‘\t’,’precision’, 6); G B B2 pause [JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber); JP iJP=-inv(JP) 3 华南理工大学广州学院 JQ iJQ=-inv(JQ) pause %maxiteration=0 for i=1:busn