PQ分解法潮流计算实验报告.docx

南昌大学实验报告

学生姓名：学号：专业班级：

实验类型：□验证□综合□设计□创新实验日期：实验成绩：

实验二P-Q分解法潮流计算实验

一、实验目的：

本实验针对高压输电网络结构特性，要求根据功率平衡方程的泰勒展开式近似的求取系统的状态变量，使系统潮流计算能够由计算机自行完成，即根据已知的电力网的数学模型（节点导纳矩阵）及各节点参数，由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解，学会运用电力系统的数学模型，掌握P-Q分解法潮流计算的过程及其特点，熟悉应用软件MATLAB，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。

二、实验内容：

编制调试P-Q分解法潮流计算的计算机程序。程序要求根据已知的电网拓扑、线路参数及节点参数，完成电力系统的潮流计算，要求计算出节点电压、线路潮流等参数。

1.根据已知电网拓扑和系统参数生成系统导纳矩阵；

2.采用P-Q分解法编制潮流程序；

3.对程序进行组织调试；

4.应用计算例题验证程序的计算效果；

5.对调试正确的计算程序进行存储、打印；

6.完成本次实验的实验报告。

三、程序清单

1、实验数据

采用课本P48的例2-1作为算例

2、程序

clearall;

clc;

n=5;%节点总数

PQ=3;%PQ节点个数

PV=1;%PV节点个数

edge=10;

error=1.0*10^(-9);%设置误差

ig=sqrt(-1);

loaddata1;

[r1,c1]=size(data1);

A=zeros(n,edge);

forj=1:r1

A(data1(j,2),j)=1;

ifdata1(j,3)0

A(data1(j,3),j)=-1;

else

end

end%生成节点关联矩阵

Y=A*diag(data1(:,4))*A;%生成节点导纳矩阵

G=real(Y);%取实部

B=imag(Y);%取虚部

vol=[1;1;1;1;1.05];%电压初值

ang=[0;0;0;0;0];%相位初值

P_set=[-0.22;-0.18;-0.27;0.35];

Q_set=[-0.14;-0.09;-0.13];%功率初值

mis=1;%误差差值

k=1;%循环次数

while(miserror)

fori=1:n-1

Temp1=0;

forj=1:n

Temp1=Temp1+vol(j,1)*(G(i,j)*cos(ang(i,1)-ang(j,1))+B(i,j)*sin(ang(i,1)-ang(j,1)));

end

Pi(i,1)=vol(i,1)*Temp1;

end

fori=1:n-1

Temp2=0;

forj=1:n

Temp2=Temp2+vol(j,1)*(G(i,j)*sin(ang(i,1)-ang(j,1))-B(i,j)*cos(ang(i,1)-ang(j,1)));

end

Qi(i,1)=vol(i,1)*Temp2;

end

fori=1:1:n-1

forj=1:1:n-1

H(i,j)=vol(i,1)*vol(j,1)*B(i,j);

end

end%生成H矩阵

fori=1:1:PQ

forj=1:1:PQ

L(i,j)=vol(i,1)*vol(j,1)*B(i,j);

e