潮流计算的MATLAB源程序Word下载.docx
《潮流计算的MATLAB源程序Word下载.docx》由会员分享，可在线阅读，更多相关《潮流计算的MATLAB源程序Word下载.docx（11页珍藏版）》请在冰豆网上搜索。
%X为节点号和对地参数矩阵。
其中第一列为节点编号，第二列为节点对地
%参数。
n=input（'
请输入节点数:
n='
）;
n1=input（'
请输入支路数:
n1='
isb=input（'
请输入平衡节点号:
isb='
pr=input（'
请输入误差精度:
pr='
B1=input（'
请输入支路参数:
B1='
B2=input（'
请输入节点参数:
B2='
X=input（'
节点号和对地参数:
X='
Y=zeros（n）;
Times=1;
%置迭代次数为初始值
%创建节点导纳矩阵
fori=1:
n1
ifB1（i,6）==0%不含变压器的支路
p=B1（i,1）;
q=B1（i,2）;
Y（p,q）=Y（p,q）-1/B1（i,3）;
Y（q,p）=Y（p,q）;
Y（p,p）=Y（p,p）+1/B1（i,3）+0.5*B1（i,4）;
Y（q,q）=Y（q,q）+1/B1（i,3）+0.5*B1（i,4）;
else%含有变压器的支路
Y（p,q）=Y（p,q）-1/（B1（i,3）*B1（i,5））;
Y（p,p）=Y（p,p）+1/B1（i,3）;
Y（q,q）二Y（q,q）+1/（B1（i,5F2*B1（i,3））;
end
Y
OrgS=zeros（2*n-2,1）;
DetaS二zeros（2*n-2,1）;
%将OrgS、DetaS初始化
%创建OrgS,用于存储初始功率参数
h=0;
j=0;
ifi~=isb&
B2（i,6）==2
h=h+1;
forj=1:
n
OrgS（2*h-1,1）=OrgS（2*h-1,1）+real（B2（i,3））*（real（Y（i,j））*real（B2（j,3））-imag（Y（i,j））*i
mag（B2（j,3）））+imag（B2（i,3））*（real（Y（i,j））*imag（B2（j,3））+imag（Y（i,j））*real（B2（j,3）））;
OrgS（2*h,1）=OrgS（2*h,1）+imag（B2（i,3））*（real（Y（i,j））*real（B2（j,3））-imag（Y（i,j））*ima
g（B2（j,3）））-real（B2（i,3））*（real（Y（i,j））*imag（B2（j,3））+imag（Y（i,j））*real（B2（j,3）））;
n%对PV节点的处理，注意这时不可再将h初始化为0
B2（i,6）==3
OrgS
%创建PVU用于存储PV节点的初始电压
PVU=zeros（n-h-1,1）;
t=0;
ifB2（i,6）==3
t=t+1;
PVU（t,1）=B2（i,3）;
PVU
%创建DetaS,用于存储有功功率、无功功率和电压幅值的不平衡量
n%对PQ节点的处理
DetaS（2*h-1,1）=real（B2（i,2））-OrgS（2*h-1,1）;
DetaS（2*h,1）=imag（B2（i,2））-OrgS（2*h,1）;
n%对PV节点的处理,注意这时不可再将h初始化为0
DetaS（2*h,1）=real（PVU（t,1））A2+imag（PVU（t,1）F2-real（B2（i,3））A2-imag（B2（i,3）F2;
DetaS
%创建I，用于存储节点电流参数
i=zeros（n-1,1）;
ifi~=isb
I（h,1）=（OrgS（2*h-1,1）-OrgS（2*h,1）*sqrt（-1））/conj（B2（i,3））;
%创建Jacbi雅可比矩阵）
Jacbi=zeros（2*n-2）;
k=0;
ifB2（i,6）==2
ifj~=isb
k=k+1;
ifi==j%对角元素的处理
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（B2（i,3））+real（Y（i,j））*imag（B2（i,3））+imag（I（h,1））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（B2（i,3））+imag（Y（i,j））*imag（B2（i,3））+real（I（h,
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））;
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（B2（i,3））+real（Y（i,j））*imag（B2（i,3））;
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（B2（i,3））+imag（Y（i,j））*imag（B2（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）;
ifk==（n-1）%将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;
n%对PV节点的处理
1））;
Jacbi（2*h,2*k-1）=2*imag（B2（i,3））;
Jacbi（2*h,2*k）=2*real（B2（i,3））;
else%非对角元素的处理
Jacbi（2*h-1,2*k）=real（Y（i,j））*real（B2（i,3））+imag（Y（i,j））*imag（B2（i,3））;
Jacbi（2*h,2*k-1）=0;
Jacbi（2*h,2*k）=0;
ifk==（n-1）%将用于内循环的指针置于初始值，以确保雅可比矩阵换
行
Jacbi
%求解修正方程，获取节点电压的不平衡量
DetaU=zeros（2*n-2,1）;
DetaU=inv（Jacbi）*DetaS;
DetaU
%修正节点电压j=0;
j=j+1;
B2（i,3）=B2（i,3）+DetaU（2*j,1）+DetaU（2*j-1,1）*sqrt（-1）;
B2
%开始循环**********************************************************************whileabs（max（DetaU））>
pr
%!
!
初始功率参数在迭代过程中是不累加的，所以在这里
必须将其初始化为零矩阵
OrgS（2*h-1,1）=OrgS（2*h-1,1）+real（B2（i,3））*（real（Y（i,j））*real（B2（j,3））-imag（Y（i,j））*imag（B2（j,3）））+imag（B2（i,3））*（real（Y（i,j））*imag（B2（j,3））+imag（Y（i,j））*real（B2（j,3）））;
endfori=1:
%创建DetaS
DetaS（2*h,1）=imag（B2（i,2））-OrgS（2*h,1）;
%创建Ii=zeros（n-1,1）;
I
%创建Jacbi
ifi==j
（I（h,1））;
Jacbi（2*h-1,2*k-1）=-imag（Y（i,j））*real（B2（i,3））+real（Y（i,j））*imag（B2（i,3））+imag
else
endifk==（n-1）
ifk==（n-1）
%修正节点电压
Times=Times+1;
%迭代次数加1
Times
一个原始数据的例子
节点数5
支路数5
平衡节点编号5
精度pr0.000001
B1（支路参数矩阵）
[120.04+0.25i0.5i10;
130.1+0.35i010;
230.08+0.30i0.5i10;
420.015i01.051;
53
0.03i01.051]
B2（节点参数矩阵）[0-1.6-0.8i1002;
0-2-1i1002;
0-3.7-1.3i1002;
05+0i1.051.0503;
001.051.0501]X（节点号和对地参数）
[10;
20;
30;
40;
50]
升级会员 