电力系统潮流计算仿真报告23页word资料.docx
《电力系统潮流计算仿真报告23页word资料.docx》由会员分享，可在线阅读，更多相关《电力系统潮流计算仿真报告23页word资料.docx（35页珍藏版）》请在冰豆网上搜索。
电力系统潮流计算仿真报告23页word资料
BeijingJiaotongUniversity
死记硬背是一种传统的教学方式,在我国有悠久的历史。
但随着素质教育的开展,死记硬背被作为一种僵化的、阻碍学生能力发展的教学方式,渐渐为人们所摒弃;而另一方面,老师们又为提高学生的语文素养煞费苦心。
其实,只要应用得当,“死记硬背”与提高学生素质并不矛盾。
相反,它恰是提高学生语文水平的重要前提和基础。
这个工作可让学生分组负责收集整理,登在小黑板上,每周一换。
要求学生抽空抄录并且阅读成诵。
其目的在于扩大学生的知识面,引导学生关注社会,热爱生活,所以内容要尽量广泛一些,可以分为人生、价值、理想、学习、成长、责任、友谊、爱心、探索、环保等多方面。
如此下去,除假期外,一年便可以积累40多则材料。
如果学生的脑海里有了众多的鲜活生动的材料,写起文章来还用乱翻参考书吗?
电力系统潮流计算
一般说来，“教师”概念之形成经历了十分漫长的历史。
杨士勋（唐初学者，四门博士）《春秋谷梁传疏》曰：
“师者教人以不及，故谓师为师资也”。
这儿的“师资”，其实就是先秦而后历代对教师的别称之一。
《韩非子》也有云：
“今有不才之子……师长教之弗为变”其“师长”当然也指教师。
这儿的“师资”和“师长”可称为“教师”概念的雏形，但仍说不上是名副其实的“教师”，因为“教师”必须要有明确的传授知识的对象和本身明确的职责。
仿真报告
姓名：
TYP
班级：
电气0906
学号：
09291183
指导老师：
吴俊勇
完成日期：
2019.6.24
一、实验内容
电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。
它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。
电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。
对于简单系统，可以将其分为开式网络和闭式网络手工计算。
对于复杂电力系统，根据定解条件，应用牛顿—拉夫逊法进行计算，在手工计算中，由于涉及大量变量、微分方程、矩阵计算，求解很烦琐，而且容易出错，计算不同系统时需要重新计算。
故而我们可以借助计算机来进行潮流计算，方便快捷且准确率高。
二、计算机潮流计算方法
我们常用牛顿—拉夫逊法来进行潮流计算。
牛顿—拉夫逊法（简称牛顿法）在数学上是求解非线性代数方程式的有效方法，其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程，即通常所称的逐次线性化过程。
1、基本原理
从几何意义上，牛顿—拉夫逊法实质上就是切线法，是一种逐步线性化的方法。
2、牛顿—拉夫逊法潮流求解过程
以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。
当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量，由于平衡节点的电压向量是给定的，因此待求量共2（n-1）需要2（n-1）个方程式。
事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。
求解过程大致可以分为以下步骤：
（1）形成节点导纳矩阵；
（2）将各节点电压设初值；
（3）将节点初值代入相关求式，求出修正方程式的常数项向量；
（4）将节点电压初值代入求式，求出雅可比矩阵元素；
（5）求解修正方程，求修正向量；
（6）求取节点电压的新值；
（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步；
（8）计算支路功率分布，PV节点无功功率和平衡节点功率。
其用程序仿真的过程可以由以下流程图来简单表示出来：
三、Matlab程序
1、原代码
%潮流计算
fprintf（'开始潮流计算\n'）;
fprintf（'请输入待求网络的相应参数\n'）;
%参数输入部分
n=input（'网络中的节点数:
n='）;
L=input（'网络中的支路数:
L='）;
ss=input（'平衡节点ss='）;
pr=input（'误差精度:
pr='）;
X1=input（'支路参数:
X1='）;
X2=input（'节点参数:
X2='）;
X=input（'节点号和对地参数:
X='）;
fprintf（'参数输入部分结束\n\n'）;
Y=zeros（n）;
%置迭代次数
mm=1;
%创建节点导纳矩阵
fori=1:
L
ifX1（i,6）==0%不含变压器的支路
p=X1（i,1）;
q=X1（i,2）;
Y（p,q）=Y（p,q）-1/X1（i,3）;
Y（q,p）=Y（p,q）;
Y（p,p）=Y（p,p）+1/X1（i,3）+0.5*X1（i,4）;
Y（q,q）=Y（q,q）+1/X1（i,3）+0.5*X1（i,4）;
else%含有变压器的支路
p=X1（i,1）;
q=X1（i,2）;
Y（p,q）=Y（p,q）-1/（X1（i,3）*X1（i,5））;
Y（q,p）=Y（p,q）;
Y（p,p）=Y（p,p）+1/X1（i,3）;
Y（q,q）=Y（q,q）+1/（X1（i,5）^2*X1（i,3））;
end
end
Y;
OrgS=zeros（2*n-2,1）;
DetaS=zeros（2*n-2,1）;%将OrgS、DetaS初始化
%创建OrgS，用于存储初始功率参数
h=0;
j=0;
fori=1:
n%对PQ节点的处理
ifi~=ss&&X2（i,6）==2
h=h+1;
forj=1:
n
OrgS（2*h-1,1）=OrgS（2*h-1,1）+real（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））+imag（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
OrgS（2*h,1）=OrgS（2*h,1）+imag（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））-real（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
end
end
end
fori=1:
n%对PV节点的处理，注意这时不可再将h初始化为0
ifi~=ss&&X2（i,6）==0
h=h+1;
forj=1:
n
OrgS（2*h-1,1）=OrgS（2*h-1,1）+real（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））+imag（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
OrgS（2*h,1）=OrgS（2*h,1）+imag（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））-real（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
end
end
end
OrgS;
%创建PVU用于存储PV节点的初始电压
PVU=zeros（n-h-1,1）;
t=0;
fori=1:
n
ifX2（i,6）==0
t=t+1;
PVU（t,1）=X2（i,3）;
end
end
PVU;
%创建DetaS，用于存储有功功率、无功功率和电压幅值的不平衡量
h=0;
fori=1:
n%对PQ节点的处理
ifi~=ss&&X2（i,6）==2
h=h+1;
DetaS（2*h-1,1）=real（X2（i,2））-OrgS（2*h-1,1）;
DetaS（2*h,1）=imag（X2（i,2））-OrgS（2*h,1）;
end
end
t=0;
fori=1:
n%对PV节点的处理，注意这时不可再将h初始化为0
ifi~=ss&&X2（i,6）==0
h=h+1;
t=t+1;
DetaS（2*h-1,1）=real（X2（i,2））-OrgS（2*h-1,1）;
DetaS（2*h,1）=real（PVU（t,1））^2+imag（PVU（t,1））^2-real（X2（i,3））^2-imag（X2（i,3））^2;
end
end
DetaS;
%创建I，用于存储节点电流参数
i=zeros（n-1,1）;
h=0;
fori=1:
n
ifi~=ss
h=h+1;
I（h,1）=（OrgS（2*h-1,1）-OrgS（2*h,1）*sqrt（-1））/conj（X2（i,3））;
end
end
I;
%创建Jacbi（雅可比矩阵）
Jacbi=zeros（2*n-2）;
h=0;
k=0;
fori=1:
n%对PQ节点的处理
ifX2（i,6）==2
h=h+1;
forj=1:
n
ifj~=ss
k=k+1;
ifi==j%对角元素的处理
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））+imag（I（h,1））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））+real（I（h,1））;
Jacbi（2*h,2*k-1）=-Jacbi（2*h-1,2*k）+2*real（I（h,1））;
Jacbi（2*h,2*k）=Jacbi（2*h-1,2*k-1）-2*imag（I（h,1））;
else%非对角元素的处理
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h,2*k-1）=-Jacbi（2*h-1,2*k）;
Jacbi（2*h,2*k）=Jacbi（2*h-1,2*k-1）;
end
ifk==（n-1）%将用于内循环的指针置于初始值，以确保雅可比矩阵换行
k=0;
end
end
end
end
end
k=0;
fori=1:
n%对PV节点的处理
ifX2（i,6）==0
h=h+1;
forj=1:
n
ifj~=ss
k=k+1;
ifi==j%对角元素的处理
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））+imag（I（h,1））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））+real（I（h,1））;
Jacbi（2*h,2*k-1）=2*imag（X2（i,3））;
Jacbi（2*h,2*k）=2*real（X2（i,3））;
else%非对角元素的处理
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h,2*k-1）=0;
Jacbi（2*h,2*k）=0;
end
ifk==（n-1）%将用于内循环的指针置于初始值，以确保雅可比矩阵换行
k=0;
end
end
end
end
end
Jacbi;
%求解修正方程，获取节点电压的不平衡量
DetaU=zeros（2*n-2,1）;
DetaU=inv（Jacbi）*DetaS;
DetaU;
%修正节点电压
j=0;
fori=1:
n%对PQ节点处理
ifX2（i,6）==2
j=j+1;
X2（i,3）=X2（i,3）+DetaU（2*j,1）+DetaU（2*j-1,1）*sqrt（-1）;
end
end
fori=1:
n%对PV节点的处理
ifX2（i,6）==0
j=j+1;
X2（i,3）=X2（i,3）+DetaU（2*j,1）+DetaU（2*j-1,1）*sqrt（-1）;
end
end
X2;
%开始循环**********************************************************************
whileabs（max（DetaU））>pr
OrgS=zeros（2*n-2,1）;%!
!
!
初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵
h=0;
j=0;
fori=1:
n
ifi~=ss&&X2（i,6）==2
h=h+1;
forj=1:
n
OrgS（2*h-1,1）=OrgS（2*h-1,1）+real（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））+imag（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
OrgS（2*h,1）=OrgS（2*h,1）+imag（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））-real（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
end
end
end
fori=1:
n
ifi~=ss&&X2（i,6）==0
h=h+1;
forj=1:
n
OrgS（2*h-1,1）=OrgS（2*h-1,1）+real（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））+imag（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
OrgS（2*h,1）=OrgS（2*h,1）+imag（X2（i,3））*（real（Y（i,j））*real（X2（j,3））-imag（Y（i,j））*imag（X2（j,3）））-real（X2（i,3））*（real（Y（i,j））*imag（X2（j,3））+imag（Y（i,j））*real（X2（j,3）））;
end
end
end
OrgS;
%创建DetaS
h=0;
fori=1:
n
ifi~=ss&&X2（i,6）==2
h=h+1;
DetaS（2*h-1,1）=real（X2（i,2））-OrgS（2*h-1,1）;
DetaS（2*h,1）=imag（X2（i,2））-OrgS（2*h,1）;
end
end
t=0;
fori=1:
n
ifi~=ss&&X2（i,6）==0
h=h+1;
t=t+1;
DetaS（2*h-1,1）=real（X2（i,2））-OrgS（2*h-1,1）;
DetaS（2*h,1）=real（PVU（t,1））^2+imag（PVU（t,1））^2-real（X2（i,3））^2-imag（X2（i,3））^2;
end
end
DetaS;
%创建I
i=zeros（n-1,1）;
h=0;
fori=1:
n
ifi~=ss
h=h+1;
I（h,1）=（OrgS（2*h-1,1）-OrgS（2*h,1）*sqrt（-1））/conj（X2（i,3））;
end
end
I;
%创建Jacbi
Jacbi=zeros（2*n-2）;
h=0;
k=0;
fori=1:
n
ifX2（i,6）==2
h=h+1;
forj=1:
n
ifj~=ss
k=k+1;
ifi==j
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））+imag（I（h,1））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））+real（I（h,1））;
Jacbi（2*h,2*k-1）=-Jacbi（2*h-1,2*k）+2*real（I（h,1））;
Jacbi（2*h,2*k）=Jacbi（2*h-1,2*k-1）-2*imag（I（h,1））;
else
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h,2*k-1）=-Jacbi（2*h-1,2*k）;
Jacbi（2*h,2*k）=Jacbi（2*h-1,2*k-1）;
end
ifk==（n-1）
k=0;
end
end
end
end
end
k=0;
fori=1:
n
ifX2（i,6）==0
h=h+1;
forj=1:
n
ifj~=ss
k=k+1;
ifi==j
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））+imag（I（h,1））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））+real（I（h,1））;
Jacbi（2*h,2*k-1）=2*imag（X2（i,3））;
Jacbi（2*h,2*k）=2*real（X2（i,3））;
else
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（X2（i,3））+real（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（X2（i,3））+imag（Y（i,j））*imag（X2（i,3））;
Jacbi（2*h,2*k-1）=0;
Jacbi（2*h,2*k）=0;
end
ifk==（n-1）
k=0;
end
end
end
end
end
Jacbi
DetaU=zeros（2*n-2,1）;
DetaU=inv（Jacbi）*DetaS;
DetaU
%修正节点电压
j=0;
fori=1:
n
ifX2（i,6）==2
j=j+1;
X2（i,3）=X2（i,3）+DetaU（2*j,1）+DetaU（2*j-1,1）*sqrt（-1）;
end
end
fori=1:
n
ifX2（i,6）==0
j=j+1;
X2（i,3）=X2（i,3）+DetaU（2*j,1）+DetaU（2*j-1,1）*sqrt（-1）;
end
end
X2
mm=mm+1;%迭代次数加1
end
mm
2、数据输入
开始潮流计算
请输入待求网络的相应参数
网络中的节点数:
n=22
网络中的支路数:
L=40
平衡节点ss=1
误差精度:
pr=0.0001
支路参数:
X1=[1,7,0.0+0.015i,0.0+0.0i,1,1.05;
2,9.,0.0+0.0217i,0.0+0.0i,1,1.075;
3,2.,0.0+0.0124i,0.0+0.0i,1,1.1;
4,19,0.0+0.064i,0.0+0.0i,1,1.025;
5,18,0.0+0.0375i,0.0+0.0i,1,1.05;
6,17,0.0+0.0337i,0.0+0.0i,1,1;
7,8,0.0106+0.074i,0.0+0.0i,1,1;
7,9,0.0147+0.104i,0.0+0.0i,1,1;
8,9,0.0034+0.0131i,0.0+0.0i,1,1;
8,22,0.0537+0.19i,0.0+0.0i,1,1;
8,22,0.0+0.1653i,0.0+0.0i,1,0;
9,10,0.0+0.002i,0.0+0.0i,1,1;
9,22,0.0559+0.218i,0.0+0.0i,1,1;
9,22,0.0+0.1954i,0.0+0.0i,1,0;
10,11,0.0+0.018i,0.0+0.0i,1,1;
11,11,0.0-1.3675i,0.0+0.0i,1,0;
11,12,0.0033+0.0343i,0.0+0.0i,1,1;
11,12,0.0+1.8797i,0.0+0.0i,1,0;
12,12,0.0-1.3675i,0.0+0.0i,1,0;
12,12,0.0-1.3675i,0.0+0.0i,1,0;
12,13,0.00245+0.0255i,0.0+0.0i,1,1;
12,13,0.0+1.395i,0.0+0.0i,1,0;
12,15,0.0+0.018i,0.0+0.0i,1,1;
13,17,0.0+0.01i,0.0+0.0i,1,1;
14,15,0.0-0.002i,0.0+0.0i,1,1;
14,19,0.0034+0.02i,0.0+0.0i,1,1;
16,16,0.0+0.5018i,0.0+0.0i,1,0;
16,17,0.0+0.001i,0.0+0.0i,1,0.975;
16,18,0.0033+0.0333i,0.0+0.0i,1,1;
16,19,0.0578+0.218i,0.0+0.0i,1,1;
16,19,0.0+0.1887i,0.0+0.0i,1,0;
16,20,0.0165+0.0662i,0.0+0.0i,1,1;
16,20,0.0+0.2353i,0.0+0.0i,1,0;
16,21,0.0374+0.178i,0.0+0.0i,1,1;
16,21,0.0+0.164i,0.0+0.0i,1,0;
19,21,0.0114+0.037i,0.0+0.0i,1,1;
20,22,0.0214+0.0859i,0
升级会员 