- 1、本文档共21页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
matlab电力系统快速解耦法潮流计算及短路计算程序
PAGE \* MERGEFORMAT 20
电力系统快速解耦法潮流分析及短路计算
一.程序设计的基本思想:
由于电力系统潮流分析中要利用到矩阵运算,复数运算,故采用matlab编程。采用文件输入,将系统的各个参数以文件的形式输入,便于程序的通用化。
本程序共有两个输入文件,分别为线路参数的文件,和已知的节点状态文件(PQ)
为了使程序不仅仅局限于计算9节点网络,在形成节点导纳的函数Yn()中,利用循环,找出线路首节点中的最大编号数,自动确定节点导纳矩阵的维数。故对于任意n节点网络,均可以计算出节点导纳矩阵
在(3)的前提下,为了使程序支持系统增加节点,增加负荷等造成的PQ参数改变,或者PQ表的加长。对程序做了如下优化。
首先,程序执行的基础是PQ表中平衡节点在第一行,接下来是PV节点,最后是PQ节点,如果系统添加节点,或者删除节点,均在PQ表的末端操作,会造成PQ表的顺序不是平衡节点、PV节点、PQ节点的顺序。故引入了seqencing()函数,其作用就是不论输入的PQ表是什么顺序,在程序读入后均按平衡-PV-PQ的顺序排列。
其次,顺序打乱的PQ表必须与支路参数表对应,故在Yn()函数中加入了两段循环体,使之对应(见相应函数体注释)
(5) 在满足了上述4个条件后,程序便可以通用化了。当然,由于水平有限,且程序未能由大量数据测试,故缺陷在所难免,这里仅是做了通用化的尝试。在本文最末附加了该程序通用化的实例。
二、潮流计算框图
输入系统的支路参数和PQ参数
计算节点导纳矩阵
设置初始电压值U0
设置收敛条件及最大迭代次数Kmax
判断是否收敛
求解修正方程得到ΔU和Δθ
对U0进行修正,U0=U0+ΔU
KKmax ?
N
Y
输出结果
结束
输出潮流不收敛信息
Y
N
开始
三.定义相应的函数
1.形成节点导纳矩阵的函数Yn()
function Y=Yn(x,y)
%定义一名为Yn的函数,其功能是自动识别输入表中节点的个数,形成相应的节点导纳矩阵
[fid,message]=fopen(x,r) ; %从x文件中读入支路参数
if fid==-1; %判断文件是否正确打开
error(message);
end;
[HeadPoint,HeadNumber, EndPoint,EndNumber,R,X,B,k]=textread(x,%s %d %s %d %f %f %f %f);
%将读入的参数处理为以列为向量的数组
fclose(fid);%关闭文件
L=length(HeadNumber); %确定输入表的行数
[Pointstyle,PointNumber,Ps,Qs,Uk,Ok]=seqencing(y);
%调用seqencing函数,引入y文件中的PQ参数
A=PointNumber;
for i=1:L; %通过以下两循环体,实现PQ参数与支路参数的编号对应
for j=1:L;
if HeadNumber(i)==j;
HeadNumber(i)=A(j);
break;
end;
end;
end;
for i=1:L;
for j=1:L;
if EndNumber(i)==j;
EndNumber(i)=A(j);
break;
end;
end;
end;
Y=zeros(L,L); %根据txt文件中数据表的长度建立空的节点导纳矩阵
for i=1:L
m=HeadNumber(i);n=EndNumber(i);
if k(i)==0; %判断是否何种元件,为输电线元件
if n~=0;
Y(m,m)=Y(m,m)+1j*B(i)+1/(R(i)+1j*X(i));
Y(n,n)=Y(n,n)+1j*B(i)+1/(R(i)+1j*X(i));
Y(m,n)=Y(m,n)-1/(R(i)+1j*X(i));
Y(n,m)=Y(n,m)-1/(R(i)+1j*X(i));
else
Y(m,m)=Y(m,m)+R(i)+1j*X(i);
end;
else %为变压器元件
if n~=0;
您可能关注的文档
- Linux内核实验报告实验7.doc
- linux实验报告24和25(AD与DA接口实验).doc
- Linux应用程序分类系统网络和管理工具.doc
- Linux下架构安全邮件服务器之Sendmail.doc
- linux实验报告29(直流电机实验).doc
- linux操作系统下的mysql服务器搭建步骤详解.doc
- linux中的weblogic性能调优.doc
- linux操作系统root管理员用户密码丢失之解决方法.doc
- Linux操作系统-基础操作-教程.doc
- linux的硬件兼容性列表.doc
- 2024.5.20高三英语高考冲刺试题.docx
- 2025高三英语读后续写说课 课件.pptx
- 2020年鲁教版一年级上册数学期末考试卷可打印.docx
- 2018年浙教版一年级上册数学平时训练试卷必考.docx
- 2020年人教版八年级物理上册第一章1.2节运动的描述巩固练习试卷(完整版).docx
- 2022年苏教版一年级上册数学同步训练试卷(word可编辑).docx
- 2020年粤沪版八年级物理上册第一章1.2节运动的描述月考试卷(通用).docx
- 2020学年苏科版一年级上册数学期末试卷【全面】.docx
- 2021年度苏科版八年级物理上册第一章1.2节运动的描述期末模拟试卷全面.docx
- 2021年度人教版一年级上册数学期中试卷【精选】.docx
文档评论(0)