- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
电力系统分析潮实验报告
南昌大学实验报告
学生姓名: 学 号: 专业班级:
实验类型:□ 验证 □ 综合 ■ 设计 □ 创新 实验日期: 实验成绩:
电力系统潮流计算实验
一、实验目的:
本实验通过对电力系统潮流计算的计算机程序的编制与调试,获得对复杂电力系统进行潮流计算的计算机程序,使系统潮流计算能够由计算机自行完成,即根据已知的电力网的数学模型(节点导纳矩阵)及各节点参数,由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算计算方法的理解,学会运用电力系统的数学模型,掌握潮流计算的过程及其特点,熟悉各种常用应用软件,熟悉硬件设备的使用方法,加强编制调试计算机程序的能力,提高工程计算的能力,学习如何将理论知识和实际工程问题结合起来。
二、实验内容:
编制调试电力系统潮流计算的计算机程序。程序要求根据已知的电力网的数学模型(节点导纳矩阵)及各节点参数,完成该电力系统的潮流计算,要求计算出节点电压、功率等参数。
在各种潮流计算的算法中选择一种,按照计算方法编制程序。
将事先编制好的电力系统潮流计算的计算程序原代码由自备移动存储设备导入计算机。
在相应的编程环境下对程序进行组织调试。
应用计算例题验证程序的计算效果。
三、实验程序:
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
您可能关注的文档
- 第十章 企业国资产法.doc
- 第七章 自相(计量经济学).ppt
- 第三章 静态电场及其边值问题解(课后题).ppt
- 画法几何(点线)习题解答.ppt
- 数字电路 第2习题解答.ppt
- 正交试验设计(互作用).ppt
- 第六章自测题和题解答.ppt
- 大学物理第13 静电场中导体与电介质.ppt
- 软件工程七章以复习与习题课.ppt
- 《水力学》第一:水静力学.ppt
- 北师版九年级下册数学精品教学课件 期末提分练案 5.1 考点梳理与达标训练.ppt
- 华师版九年级上册数学精品教学课件 第22章 一元二次方程 阶段综合训练【范围22.3】 (2).ppt
- 北师版九年级下册数学精品教学课件 第2章 2.2.1 二次函数y=x2,y=-x2 的图象与性质.ppt
- 华师版九年级上册数学精品教学课件 第23章 阶段综合训练【范围:23.4~23.6】.ppt
- 华师版九年级上册数学精品教学课件 第23章 图形的相似 (4).ppt
- 华师版九年级上册数学精品教学课件 第23章 23.3.2 第3课时 相似三角形的判定定理3.ppt
- 华师版九年级上册数学精品教学课件 第23章 图形的相似 23.3 相似三角形 4.相似三角形的应用 (2).ppt
- 华师版九年级上册数学精品教学课件 第24章 24.1 测量.ppt
- 北师版九年级下册数学精品教学课件 第2章 2.2 二次函数的图象与性质 2.2.1 二次函数y=x2,y=-x2 的图象与性质.ppt
- 苏科版九年级上册数学精品教学课件 第4章 4.2.1概率及概率的计算.ppt
文档评论(0)