潮流计算的计算机算法Word格式.docx
- 文档编号:20534992
- 上传时间:2023-01-23
- 格式:DOCX
- 页数:24
- 大小:101.04KB
潮流计算的计算机算法Word格式.docx
《潮流计算的计算机算法Word格式.docx》由会员分享,可在线阅读,更多相关《潮流计算的计算机算法Word格式.docx(24页珍藏版)》请在冰豆网上搜索。
q=B1(i,2);
Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j);
%
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;
Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;
end
非对角元
对角元
disp('
导纳矩阵Y='
disp(Y)
%
%下面是求P,Q,V,O矩阵
V=zeros(1,n);
O=zeros(1,n);
P=zeros(1,n);
Q=zeros(1,n);
G=real(Y);
B=imag(Y);
n
P(i)=B2(i,3);
Q(i)=B2(i,4);
V(i)=B2(i,5);
O(i)=B2(i,6);
end
B3=B(1:
n-1,1:
n-1);
B4=B(1:
n2,1:
n2);
%%
不含平衡节点,由节点导纳虚部构成
所有PQ节点
%下面是求ΔP,ΔQ矩阵
DX=0;
ICT=1;
Mp=1;
Mq=1;
whileICT~=0
m1=1;
m2=1;
所有节点数
节点数
求DP,DQ
对V矩阵求逆
ΔP/V
ΔP/V/B3
Δ角=-ΔP/V/V/B3
ΔV=-ΔQ/V/B
ifi~=isb
C(i)=0;
D(i)=0;
forj1=1:
C(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1))+B(i,j1)*sin(O(i)-O(j1)));
D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*sin(O(i)-O(j1))-B(i,j1)*cos(O(i)-O(j1)));
DP(m1)=P(i)-C(i);
m1=m1+1;
ifB2(i,2)==1DQ(m2)=Q(i)-D(i);
m2=m2+1;
endendm1=m1-1;
%m2=m2-1;
%PQDPQ=[DP'
;
DQ'
];
%V1=V(:
1:
m1);
V2=diag(V1);
V3=inv(V2);
%H=V3*DP'
%K=-inv(B3)*H;
%-deltO=V3*K;
%max1=max(abs(DP));
fori=1:
m1
ifmax1<
pr
Mp=0;
else
O(i)=O(i)+deltO(i)'
endendV4=V(:
m2);
V5=diag(V4);
V6=inv(V5);
L=V6*DQ'
N=-inv(B4)*L;
deltV=N;
%max2=max(abs(DQ));
m2
ifmax2<
Mq=0;
ifB2(i,2)==1;
V(i)=V(i)+deltV(i)'
ifMp==0&
&
Mq==0
ICT=0;
DX=DX+1;
%迭代结束,开始输出结果
'
迭代次数为:
disp(DX);
E(i)=V(i)*cos(O(i))+1j*V(i)*sin(O(i));
o(i)=180*angle(E(i))/pi;
修正后各节点电压标么值为(节点号从小到大排列):
'
disp(V);
修正后各节点电压相角为(节点号从小到大排列):
disp(o);
%计算各个节点的功率
各节点的功率为:
forp=1:
C(p)=0;
forq=1:
C(p)=C(p)+conj(Y(p,q)*conj(E(q)));
S(p)=E(p)*C(p);
disp(S);
%计算各支路的功率
n1p=B1(i,1);
Si(p,q)=E(p)*(conj(E(p))*conj(Y(p,p)-Y(p,q))+(conj(E(p))-conj(E(q)))*conj(Y(p,q)));
各条支路的首端功率为:
disp(Si(p,q));
Si(q,p)=E(q)*(conj(E(q))*conj(Y(q,q)-Y(p,q))+(conj(E(q))-conj(E(p)))*conj(Y(p,q)));
各条支路的末端功率为:
disp(Si(q,p));
DS(i)=Si(p,q)+Si(q,p);
各条支路的功率损耗为:
disp(DS(i));
%计算平衡节点功率
Sp=0;
Sp=Sp+V(n)*conj(Y(n,i))*conj(V(i));
平衡节点功率为:
disp(Sp);
、对支路参数(B1)、节点参数(B2)的说明
1.节点数:
14
2.支路数:
20
3.支路矩阵B1的各支路参数:
起点编号,终点编号,电阻,电抗,电导,电纳
[1
2
0.013350.042110
0;
1
3
00.2091200
4
00.5561800
10
0.058110.176320
0.034
11
0.067010.171030
0.0128
0.056950.173880
0.0346
12
00.2520200
14
0.054030.223040
0.0492
00.1100100
13
00.1761500
5
0.031810.08450
9
0.127110.270380
6
0.082050.192070
0.094980.19890
7
8
0.220920.199880
0.122910.255810
0.170930.348020
0.066150.130270
0.046990.197970
0.0438
0.019380.059170
0.0528
4.
节点参数矩阵B2的各节点参数:
(对应的每一列为)节点编号,类型,注入有功,注入无功,电压幅值,电压相位其中节点类型:
1=PQ节点,2=PV节点,0=平衡节点
>
clear
请输入节点数:
n1=input('
请输入支路数:
isb=input('
请输入平衡节点号:
pr=input('
请输入误差精度:
请输入支路参数:
B1='
请输入节点参数:
n2=input('
请输入PQ节点个数:
%非对角元
%对角元
%下面是求P,Q,V,O矩阵
%不含平衡节点,由节点导纳虚部构成
%所有PQ节点
%下面是求ΔP,ΔQ矩阵
C(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1))+B(i,j1)*sin(O(i)-O(j1)));
D(i)=D(i)+V(i)*V(j1)*(G(i,j1)*sin(O(i)-O(j1))-B(i,j1)*cos(O(i)-O(j1)));
max1=max(abs(DP));
V4=V(:
V5=diag(V4);
V6=inv(V5);
L=V6*DQ'
N=-inv(B4)*L;
deltV=N;
%ΔV=-ΔQ/V/B
max2=max(abs(DQ));
elseifB2(i,2)==1;
V(i)=V(i)+deltV(i)'
Mp=1;
o(i)=180*angle(E(i))/pi;
%计算各个节点的功率
enddisp(S);
Si(p,q)=E(p)*(conj(E(p))*conj(Y(p,p)-Y(p,q))+(conj(E(p))-conj(E(q)))*conj(Y(p,q)));
disp('
disp(Si(p,q));
Si(q,p)=E(q)*(conj(E(q))*conj(Y(q,q)-Y(p,q))+(conj(E(q))-conj(E(p)))*conj(Y(p,q)));
disp(Si(q,p));
n=14
n1=20
请输入平衡节点号:
isb=14
pr=0.00001
B1=[120.013350.0421100;
00.209120
00.556180
0.05811
0.17632
0.06701
0.17103
0.05695
0.17388
00.252020
0.05403
0.22304
00.110010
00.176150
0.03181
0.0845
0.12711
0.27038
0.08205
0.19207
0.09498
0.1989
0.22092
0.19988
0.12291
0.25581
0.17093
0.34802
0.06615
0.13027
0.04699
0.19797
0.0438;
0.01938
0.05917
0.0528;
B2=[11
-0.4780.039
-0.076
-0.016
0;
00
10;
-0.295
-0.166
-0.09
-0.058
-0.035
-0.018
-0.061
-0.135
-0.149
-0.05
0.183
01.045
-0.942
01.01
-0.112
0.047
1.7
00.1741.9
1.060;
]
请输入PQ节点个数:
n2=9
导纳矩阵
Y=
Columns1through5
10.5130-38.2963i
-6.8410+21.5786i
0.0000+4.7819i
0.0000+1.7980i
0.0000+
0.0000i
9.5680-34.8916i
0.0000+0.0000i
0.0000-19.5490i
0.0000+9.0901i
5.3261-24.2825i
-3.9020
+10.3654i
-3.9020+10.3654i
5.7829
-14.7683i
-1.8809+
4.4029i
Columns6through10
-1.6860+
-1.7011+
9.5213
-1.1350+
-4.9991
-1.6860+5.1158i
-1.9860+5.0688i
5.1158i
5.1939i
-1.8809+4.4029i
3.8359-8.4970i
0.0000+0.0000i-30.1895i
4.7819i
-1.9550+4.0941i
-1.7011+5.1939i
0.0000+3.9679i
-1.0259+4.2350i
4.0150-5.4279i
-2.4890+2.2520i
-1.5260+3.1760i
0.0000+5.6770i
6.7249-10.6697i
-1.1370+2.3150i
-3.0989+6.1028i
-1.4240+3.0291i
2.5610-5.3440i
+15.2631i
Columns11through14
迭代次数为:
42
修正后各节点电压标么值为
Columns1through10
(节点号从小到大排列)
1.21281.2148
1.58091.0450
1.5627
1.5364
1.5602
1.6264
1.6792
1.6654
1.01001.7000
1.9000
1.0600
修正后各节点电压相角为(节点号从小到大排列)
-1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 潮流 计算 计算机 算法