高斯赛德尔法潮流计算.docx
- 文档编号:6134941
- 上传时间:2023-01-04
- 格式:DOCX
- 页数:14
- 大小:19.95KB
高斯赛德尔法潮流计算.docx
《高斯赛德尔法潮流计算.docx》由会员分享,可在线阅读,更多相关《高斯赛德尔法潮流计算.docx(14页珍藏版)》请在冰豆网上搜索。
高斯赛德尔法潮流计算
....
高斯——赛德尔法潮流计算
潮流计算高斯——赛德尔迭代法(Gauss一Seidelmethod)是求解电力系统
潮流的方法。
潮流计算高斯——赛德尔迭代法又分导纳矩阵迭代法和阻抗矩阵迭代法两种。
前者是以节点导纳矩阵为基础建立的赛德尔迭代格式;后者是以节
点阻扰矩阵为基础建立的赛德尔迭代格式。
高斯——赛德尔迭代法这是数学上求解线性或非线性方程组的一种常用的迭代方法。
本实验通过对电力网数学模型形成的计算机程序的编制与调试,获得形成电力网数学模型:
高斯---赛德尔法的计算机程序,使数学模型能够由计算机自行形成,即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网
的节点导纳矩阵和各节点电压、功率。
通过实验教学加深学生对高斯---赛德尔法概念的理解,学会运用数学知识建立电力系统的数学模型,掌握数学模型的形成过程及其特点,熟悉各种常用应用软件,熟悉硬件设备的使用方法,加强编制调试计算机程序的能力,提高工程计算的能力,学习如何将理论知识和实际工程问题结合起来。
高斯---赛德尔法潮流计算框图
开始
输入数据,定义数组
给定PQ节点电压初值
给定PV节点电压实部(或虚部)
置迭代计数b=0
计算PQ节点电压实部和虚部
先计算PV节点无功功率
再用其计算PV节点电压实部和
虚部
计算平衡节点的有功和无
功
.....
....
求=+
N
判断所有||是否<0.000001b=b+1
Y
结果输出
结束
[1]系统节点的分类
根据给定的控制变量和状态变量的不同分类如下
①P、Q节点(负荷节点),给定Pi、Qi求Vi、Si,所求数量最多;
②负荷节点,变电站节点(联络节点、浮游节点),给定PGi、QGi的发电机节点,给定QGi的无功电源节点;
③PV节点(调节节点、电压控制节点),给定Pi、Qi求Qn、Sn,所求数量少,可以无有功储备的发电机节点和可调节的无功电源节点;
④平衡节点(松弛节点、参考节点(基准相角)、S节点、VS节点、缓冲节点),给定Vi,δi=0,求Pn、Qn(Vs、δs、Ps、Qs)。
[2]潮流计算的数学模型
1)线性的节点电压方程YV=I
根据S=V可得非线性的节点电压方程(为I的共轭)YV=I==节点功率与节点电流的关系:
2)在国外,对于复数变量不打点,其模要加绝对值符号;在国,对于复数变量,在S、V、I上要打点,Y、Z上不打点,其模不加绝对值符号。
3)式2—5
对于发电机Pi、Qi为正,对负荷来说Pi、Qi为负
4)展开YV=I得
上式代入式2—5得n维的非线性复数电压方程组
式2—6
该式为潮流计算的基本方程
[3]高斯—赛德尔法潮流计算
.....
....
1)高斯法潮流计算
①将式2—6展开成电压方程
式2—7
假设系统节点数是n,PQ节点数为m,m+1及之后的节点是PV节点,第n个节点是平衡节点。
展开式2—7得高斯法潮流计算的基本方程
式2—8
②考虑到i=1时matlab中for语句的使用可写成
③由于平衡节点的电压和相角给定,不用计算,只要计算i=1—n-1节点的电压,但平衡节点的参数和变量要用于其他节点的电压计算.
式2—8的计算过程中有i=1、
2、···n-1
④特点:
在计算i节点的k+1次电压时,所用的i节点前后(包括i节点)
的电压都是k次迭代的结果。
2)高斯—赛德尔法潮流计算
①在高斯法潮流计算中引入赛德尔法迭代方式即为高斯—赛德尔法潮流计
算
②对应式2—8的高斯—赛德尔法潮流计算的方程为
式2—9
在式2—9的计算中有
③特点:
在计算i节点的k+1次电压时,1~i-1节点的电压用的是k+1次时的电压,而i~n-1节点的电压用的是k次时的电压,即在迭代过程中每个被求的电压新值立即被带入到下一个电压新值的计算中。
3)基于导纳矩阵的直角坐标高斯—赛德尔法潮流计算
①设展开式2—6并将实、虚部分列
式2—10
式2—11
②令
式2—12
.....
....
注:
、中不包括j=i的参数和变量;
、中分别有k+1次和k次的变量;
在、中没有单独列出。
()
③
④将2—12代入式2—10和2—11得
式2—13式2—14
⑤将式2—9展开,实、虚部分列,再将式2—12代入,得节点电压的实部、虚
部
式2—15
式2—16
⑥对P、V节点,根据常数
式2—17
4)部分求解方程
对于P、Q节点:
用式2—15求,用式2—16求
对于P、V节点:
用式2—14求用式2—15求,式2—16求
5)为了加速收敛,引入加速因子α,α=1~1.8之间,复数电压:
式2—18
6)实数模型:
式2—19
)式2—20
、是式2—15~式2—17计算出的值,、是考虑到α修正后的值,、是上一次用于迭代的实际值(不一定是式2—15~式2—17计算出的值)
7)三种加速过程
①每次求出的、立即用于求解下一个电压新值;
②每次求出的、同时立即用α进行修正,得到的、同时用于求解下一个电压新值;③每次求出的、分别用α进行修正,得到的、分别用于求解下一个电压新值。
注:
三种加速过程中,速度又快到慢依次为③②①。
8)收敛判据:
复数模型:
实数模型:
,
9)三种收敛判据情况:
①用前后两次经α修正后的电压值;
②用前后两次式2—15~式2—17计算出来的值;
③前一次用α修正的值,后一次用式2—15~式2—17计算出的值。
10)高斯—赛德尔法是用前后两次迭代的最大电压误差作收敛判据,ε取
.....
....
10-5~10-6,牛顿法是用最大功率误差为收敛判据,ε取10-3~10-5,所以后者为好。
[4]编程程序步骤如下
第一步:
设定初值
Vmax
0,i
1定义Z矩阵,s
设定循环次数k100
第二步:
用一判据(Z(i,2)
0)先求PQ节点
用2-15
式求ei(k1)
,再代入2-16替代ei(k)求fi(k1)。
?
(k1)
ei(k1)
jfi(k
1);
Vi
则Vi(k1)
Vi(k1)
Vi(k);
(k1)
arctg(fi
(k1)
(k
1)
);
i
ei
if
Vi(k1)
Vmax;
根据收敛判据Vmax10
5输出代求量,即
Vmax
Vi(k1);
if
Vmax
106
Z(i,6)
Vi(k
1);
Z(i,7)
i(k
1);
第三步:
(Z(i,2)!
0)求PV节点
用2-14求Qi
再用2-16求fi(k1)
,将其代入2-17,求ei(k1)
,
则
?
(k1)
ei(k1)
jfi(k
1);
Vi
Vi(k1)
Vi(k1)
Vi(k);
i(k1)
arctg(fi(k1)
/ei(k1));
if
Vi(k1)
Vmax
根据收敛判据Vmax10
6输出代求量,即
Vi(k1)
Vmax
;
if
Vmax
106
Z(i,3)
Qi(k1);
Z(i,7)
i(k1);
第四步:
求平衡节点n
.....
....
利用式2-13和2-14式求Pi和Qi,然后输出,
即
Z(i,2)
Pi;
Z(i,3)
Qi;
最后输出Z矩阵
试验题目:
用形成Y阵的五节点系统,假定节点
1、2、3为PQ节点,节点4
为PV节点、节点5为平衡节点,试分别用高斯—赛德尔法潮流计算其潮流。
取收敛判据为|△Vmax|<106。
给定:
程序如下:
clear
clc
I=[-2,-3,2,2,3];
J=[4,5,3,1,1];
R=[0,0,0.08,0.04,0.1];
X=[0.015,0.03,0.3,0.25,0.35];
K=[1.05,1.05,0.25,0.25,0];
n=5;
L=5;
Y=zeros(2*n,n);
form=1:
L
i=I(m);j=J(m);r=R(m);x=X(m);k=K(m);
ifi*j==0
Y(2*i-1,i)=Y(2*i-1,i)+r;
Y(2*i,i)=Y(2*i,i)-x;
end
ifi*j>0
Y(2*i-1,j)=Y(2*i-1,j)-r/(r^2+x^2);
.....
....
Y(2*i,j)=Y(2*i,j)+x/(r^2+x^2);
Y(2*j-1,i)=Y(2*i-1,j);
Y(2*j,i)=Y(2*i,j);
Y(2*i-1,i)=Y(2*i-1,i)+r/(r^2+x^2);
Y(2*i,i)=Y(2*i,i)-x/(r^2+x^2)+k;
Y(2*j-1,j)=Y(2*j-1,j)+r/(r^2+x^2);
Y(2*j,j)=Y(2*j,j)-x/(r^2+x^2)+k;
end
ifi*j<0
i=-i;
Y(2*i-1,j)=Y(2*i-1,j)-r/(r^2+x^2)/k;
Y(2*i,j)=Y(2*i,j)+x/(r^2+x^2)/k;
Y(2*j-1,i)=Y(2*i-1,j);
Y(2*j,i)=Y(2*i,j);
Y(2*i-1,i)=Y(2*i-1,i)+r/(r^2+x^2)/k^2;
Y(2*i,i)=Y(2*i,i)-x/(r^2+x^2)/k^2;
Y(2*j-1,j)=Y(2*j-1,j)+r/(r^2+x^2);
Y(2*j,j)=Y(2*j,j)-x/(r^2+x^2);
end
end
Y
P=[-1.6,-2.0,-3.7,5.0,0];
Q=[-0.8,-1.0,-1.3,0,0];
E=[1,1,1,1.05,1.05];
F=[0,0,0,0,0];
k=0;
V=[1,1,1,1.05,1.05];
A=[0,0,0,0,0];
h=3;
m=0.000001;
Vm=1;
whileVm>m
Vm=0;
fori=1:
n-1
j=1;
A1=0;
A2=0;
ifi>j
forj=1:
i-1
g=Y(2*i-1,j);
b=Y(2*i,j);
e=E(j);
f=F(j);
A1=A1+g*e-b*f;
.....
....
A2=A2+g*f+b*e;
end
end
forj=i+1:
n
g=Y(2*i-1,j);
b=Y(2*i,j);
e=E(j);
f=F(j);
A1=A1+g*e-b*f;
A2=A2+g*f+b*e;
end
e=E(i);f=F(i);p=P(i);q=Q(i);g=Y(2*i-1,i);b=Y(2*i,i);
ifi>h
g=Y(2*i-1,i);
b=Y(2*i,i);
Q(i)=-b*(e^2+f^2)-e*A2+f*A1;
q=Q(i);
E(i)=g/(g^2+b^2)*((p*e+q*f)/(e^2+f^2)-A1)+b/(g^2+b^2)*((p*f-q*e)/(e^2+f^2)-A2);
v=V(i);
F(i)=sqrt(v^2-E(i)^2);
A(i)=atan(F(i)/E(i));
A(i)=A(i)*180/pi;
continue
end
E(i)=g/(g^2+b^2)*((p*e+q*f)/(e^2+f^2)-A1)+b/(g^2+b^2)*((p*f-q*e)/(e^2+f^2)-A2);
F(i)=g/(g^2+b^2)*((p*f-q*e)/(e^2+f^2)-A2)+b/(g^2+b^2)*((p*e+q*f)/(e^2+f^2)-A1);
v=sqrt(E(i)^2+F(i)^2);
Vc=v-V(i);
Vc=abs(Vc);
ifVc>Vm
Vm=Vc;
end
V(i)=v;
A(i)=atan(F(i)/E(i));
A(i)=A(i)*180/pi;
end
k=k+1;
end
.....
....
forj=1:
n
e=E(j);f=F(j);g=Y(2*i-1,j);b=Y(2*i,j);
P(n)=P(n)+E(n)*(g*e-b*f);
Q(n)=Q(n)-E(n)*(g*f+b*f);
end
k
P
Q
V
A
运行结果:
Y=
1.3787
-0.6240
-0.7547
0
0
-6.2917
3.9002
2.6415
0
0
-0.6240
1.4539
-0.8299
0
0
3.9002-66.9808
3.1120
63.4921
0
-0.7547
-0.8299
1.5846
0
0
2.6415
3.1120-35.7379
031.7460
0
0
0
0
0
0
63.4921
0-66.6667
0
0
0
0
0
0
0
0
31.7460
0-33.3333
k=
11
P=
-1.6000-2.0000-3.70005.00000.5238
Q=
-0.8000-1.0000-1.30001.38850.5238
V=
0.88851.08171.05791.05001.0500
.....
....
A=
-11.6107-0.41331.17980.00280
>>
.....
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高斯赛德尔法 潮流 计算
![提示](https://static.bdocx.com/images/bang_tan.gif)