电力系统分析潮流实验报告(优选.)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
最新文件---------------- 仅供参考--------------------已改成-----------word文本 --------------------- 方便更改
赠人玫瑰,手留余香。
南昌大学实验报告
学生姓名:学号:专业班级:
实验类型:□验证□综合■设计□创新实验日期:实验成绩:
电力系统潮流计算实验
一、实验目的:
本实验通过对电力系统潮流计算的计算机程序的编制与调试,获得对复杂电力系统进行潮流计算的计算机程序,使系统潮流计算能够由计算机自行完成,即根据已知的电力网的数学模型(节点导纳矩阵)及各节点参数,由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解,学会运用电力系统的数学模型,掌握潮流计算的过程及其特点,熟悉各种常用应用软件,熟悉硬件设备的使用方法,加强编制调试计算机程序的能力,提高工程计算的能力,学习如何将理论知识和实际工程问题结合起来。
二、实验内容:
编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型(节点导纳矩阵)及各节点参数,完成该电力系统的潮流计算,要求计算出节点电压、功率等参数。
1、在各种潮流计算的算法中选择一种,按照计算方法编制程序。
2、将事先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。
3、在相应的编程环境下对程序进行组织调试。
4、应用计算例题验证程序的计算效果。
三、实验程序:
function [e,f,p,q]=flow_out(g,b,kind,e,f)
%计算潮流后efpq的终值
s=flow(g,b,kind,e,f);
k=0;
while max(abs(s))>10^-5
J=J_out(g,b,kind,e,f);
J_ni=inv(J);
dv=J_ni*s;
l=length(dv)/2;
for i=1:l
e(i)=e(i)-dv(2*i-1);
f(i)=f(i)-dv(2*i);
end
s=flow(g,b,kind,e,f);
end
l=length(e);
for i=1:l
s1=0;
s2=0;
for j=1:l
s1=s1+g(i,j)*e(j)-b(i,j)*f(j);
s2=s2+g(i,j)*f(j)+b(i,j)*e(j);
end
p(i)=e(i)*s1+f(i)*s2;
q(i)=f(i)*s1-e(i)*s2;
end
function s=flow(g,b,kind,e,f)
%计算当前ef与规定的pqv的差值
l=length(e);
s=zeros(2*l-2,1);
for i=1:(l-1)
s1=0;
s2=0;
for j=1:l
s1=s1+g(i,j)*e(j)-b(i,j)*f(j);
s2=s2+g(i,j)*f(j)+b(i,j)*e(j);
end
s(2*i-1)=kind(2,i)-e(i)*s1-f(i)*s2;
if kind(1,i)==1
s(2*i)=kind(3,i)-f(i)*s1+e(i)*s2;
else
s(2*i)=kind(3,i)^2-f(i)^2-e(i)^2;
end
end
function J=J_out(g,b,kind,e,f)
%计算节点的雅克比矩阵
l=length(e);
J=zeros(2*l-2,2*l-2);
for i=1:(l-1);
if kind(1,i)==1
s=PQ_out(g,b,e,f,i);
for j=1:(2*l-2)
J(2*i-1,j)=s(1,j);
J(2*i,j)=s(2,j);
end
else
s=PV_out(g,b,e,f,i);
for j=1:(2*l-2)
J(2*i-1,j)=s(1,j);
J(2*i,j)=s(2,j);
end
end
end
function pq=PQ_out(g,b,e,f,i)
%计算pq节点的雅克比矩阵
l=length(e);
pq=zeros(2,2*l-2);
for j=1:(l-1)
if j==i
s=0;
for k=1:l
s=s-(g(i,k)*e(k)-b(i,k)*f(k));
end
pq(1,2*i-1)=s-g(i,i)*e(i)-b(i,i)*f(i); s=0;
for k=1:l
s=s-(g(i,k)*f(k)+b(i,k)*e(k));
end
pq(1,2*i)=s+b(i,i)*e(i)-g(i,i)*f(i);
s=0;
for k=1:l
s=s+(g(i,k)*f(k)+b(i,k)*e(k));
end
pq(2,2*i-1)=s+b(i,i)*e(i)-g(i,i)*f(i); s=0;
for k=1:l
s=s-(g(i,k)*e(k)-b(i,k)*f(k));
end
pq(2,2*i)=s+g(i,i)*e(i)+b(i,i)*f(i);
else
pq(1,2*j-1)=-(g(i,j)*e(i)+b(i,j)*f(i)); pq(1,2*j)=b(i,j)*e(i)-g(i,j)*f(i);
pq(2,2*j)=-pq(1,2*j-1);
pq(2,2*j-1)=pq(1,2*j);
end
end
function pv=PV_out(g,b,e,f,i)
%计算pv节点的雅克比矩阵
l=length(e);
pv=zeros(2,2*l-2);
for j=1:(l-1)
if j==i
s=0;
for k=1:l
s=s-(g(i,k)*e(k)-b(i,k)*f(k));
end
pv(1,2*i-1)=s-g(i,i)*e(i)-b(i,i)*f(i); s=0;
for k=1:l
s=s-(g(i,k)*f(k)+b(i,k)*e(k));