基于MATLAB的潮流计算Word格式文档下载.docx
《基于MATLAB的潮流计算Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《基于MATLAB的潮流计算Word格式文档下载.docx（9页珍藏版）》请在冰点文库上搜索。
a=zeros（n,1）;
Ps=zeros（n,1）;
Qs=zeros（n,1）;
P=zeros（n,1）;
Q=zeros（n,1）;
p=zeros（n-1,1）;
q=zeros（m,1）;
aa=zeros（n-1,1）;
u=zeros（m,1）;
k=0;
I=sqrt（-1）;
%由于后面用到ij作为循环变量故重定义虚数单位
Y=[6.311-I*20.3722-2.7523+I*9.1743-3.5587+I*11.387900;
-3.5587+I*11.38798.5587-I*30.9993-5+I*15I*4.98890;
-2.7523+I*9.1743-5+I*157.7523-I*28.75570I*4.9889;
0I*4.98890-I*5.24930;
00I*4.98890-I*5.2493];
%%%%%%%%%%%%%初值的设定%%%%%%%%%%%%%%%%
U=[1;
1;
1];
a=[0;
0;
0];
Ps=[-0.8055;
-0.18;
0.5;
Qs=[-0.532;
-0.12;
3
%%%%%%%%%%%%求失配功率△P（k）和△Q（k）%%%%%%%%%%%%%%
fori=1:
n-1
s=0;
forj=1:
n
s=s+U（j）*（real（Y（i,j））*cos（a（i）-a（j））+imag（Y（i,j））*sin（a（i）-a（j）））;
end
P（i）=U（i）*s;
s=s+U（j）*（real（Y（i,j））*sin（a（i）-a（j））-imag（Y（i,j））*cos（a（i）-a（j）））;
Q（i）=U（i）*s;
detp（i）=Ps（i）-P（i）;
m
detq（i）=Qs（i）-Q（i）;
%%%%%%%求解雅可比矩阵%%%%%%%%%%%
while（max（[detp]）>
=0.0001|max（[detq]）>
=0.0001）%设定精确度
H=zeros（n-1,n-1）;
N=zeros（n-1,m）;
M=zeros（m,n-1）;
L=zeros（m,m）;
ifi==j
4
H（i,j）=U（i）^2*imag（Y（i,j））+Q（i）;
else
H（i,j）=-U（i）*U（j）*（real（Y（i,j））*sin（a（i）-a（j））-imag（Y（i,j））*cos（a（i）-a（j）））;
N（i,j）=-U（i）^2*real（Y（i,j））-P（i）;
N（i,j）=-U（i）*U（j）*（real（Y（i,j））*cos（a（i）-a（j））+imag（Y（i,j））*sin（a（i）-a（j）））;
M（i,j）=U（i）^2*real（Y（i,j））-P（i）;
M（i,j）=U（i）*U（j）*（real（Y（i,j））*cos（a（i）-a（j））+imag（Y（i,j））*sin（a（i）-a（j）））;
L（i,j）=U（i）^2*imag（Y（i,j））-Q（i）;
L（i,j）=-U（i）*U（j）*（real（Y（i,j））*sin（a（i）-a（j））-imag（Y（i,j））*cos（a（i）-a（j）））;
5
%%%%%%求解修正方程得到修正量aa,u%%%%%%%%%%%%%%
J=-inv（[H,N;
M,L]）*[detp'
;
detq'
];
aa=J（1:
n-1,:
u=J（n:
n-1+m,:
n-1%计算a（节点电压相角）,U（节点电压大小）
a（i）=a（i）+aa（i）;
U（i）=U（i）+U（i）*u（i）;
[detp'
]%迭代过程中失配功率变化情况
k=k+1;
n-1%计算功率误差p,q
6
A=a*180/pi;
%对相角进行转换由弧度制装换为角度
display（'
最终的计算结果如下:
'
A
U
各节点的电压为:
fprintf（'
节点%d:
%d∠（%d）°
\n\n'
i,U（i）,A（i））;
失配功率为:
detp
detq
迭代的次数为:
%d\n'
k）%输出迭代次数
2运行结果
>
fengjian
请输入PQ节点个数:
ans=-0.8055
-0.1800
0.5000
-0.3420
0.2575
0.4075
7
ans=
0.0347
0.0022
-0.0094
-0.0178
-0.0138
-0.0407
-0.0447
1.0e-003*
0.5087
-0.1291
-0.2188
-0.3127
-0.4229
-0.4416
A=-7.2049
-5.6096
-5.3697
-0.0056
U=1.0027
1.0263
1.0317
1.0000
8
节点1:
1.002729e+000∠（-7.204950e+000）°
节点2:
1.026321e+000∠（-5.609609e+000）°
节点3:
1.031663e+000∠（-5.369712e+000）°
节点4:
1∠（-5.622073e-003）°
节点5:
1∠（0）°
detp=
1.0e-006*
0.1069-0.0529-0.0548-0.0083
detq=
1.0e-007*
-0.7838-0.5117-0.5255
3
二P—Q解耦迭代
1程序代码
n=5;
m=3;
%%%%%%%%%%%%%%%%各参数初始化%%%%%%%%%%%%%%%%%
%导纳矩阵
%电压矢量
%相角矢量
9
%%%%%%%%%%%节点导纳矩阵的生成%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%初值的设定%%%%%%%%%%%%%%%%%
%%%%%%%%%%%求失配功率△P（k）和△Q（k）%%%%%%%%%%%%%%
10
%%%%%%%%%%%%%%%%生成B'
和B"
矩阵%%%%%%%%%%%%%%
Bp=zeros（n-1,n-1）;
Bpp=zeros（m,m）;
n-1;
Bp（i,j）=imag（Y（i,j））;
end;
m;
Bpp（i,j）=imag（Y（i,j））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
aa=（-1）*inv（Bp）*[detp
（1）/U
（1）;
detp
（2）/U
（2）;
detp（3）/U（3）;
detp（4）/U（4）];
u=（-1）*inv（Bpp）*[detq
（1）/U
（1）;
detq
（2）/U
（2）;
detq（3）/U（3）];
%%%%%%%%%%%%%%求解修正方程得到修正量
aa,u%%%%%%%%%%%%%%%%%%%%%%%%%%
U（i）=U（i）+u（i）;
11
%[detp'
%（记录迭代次数）
生成的节点导纳矩阵为'
Y
12
算得矩阵B’为'
Bp
算得矩阵B”为'
Bpp
2运行结果
fengjian1
A=
-7.2049
-5.6097
-0.0057
U=
1.0027
生成的节点导纳矩阵为
Y=6.3110-20.3722i-2.7523+9.1743i-3.5587+11.3879i00
13
-3.5587+11.3879i8.5587-30.9993i-5.0000+15.0000i0+4.9889i0
-2.7523+9.1743i-5.0000+15.0000i7.7523-28.7557i00+4.9889i
00+4.9889i00-5.2493i0
000+4.9889i00-5.2493i
算得矩阵B’为
Bp=
-20.37229.174311.38790
11.3879-30.999315.00004.9889
9.174315.0000-28.75570
04.98890-5.2493
算得矩阵B”为
Bpp=
-20.37229.174311.3879
11.3879-30.999315.0000
9.174315.0000-28.7557
1.002732e+000∠（-7.204948e+000）°
1.026322e+000∠（-5.609674e+000）°
1.031663e+000∠（-5.369742e+000）°
1∠（-5.676340e-003）°
detp=
1.0e-004*
-0.27470.27610.0562-0.0119
-0.32720.03110.2931
8
14
升级会员 