untitled.asv

可用传输容计算程序

clc;clear;
D=[1	1	2	0.0575	120	-17.3913; 
2	1	3	0.1852	120	-5.3996; 
3	2	4	0.1737	80	-5.7571; 
4	3	4	0.0379	100	-26.3852; 
5	2	5	0.1983	100	-5.0429; 
6	2	6	0.1763	100	-5.6721; 
7	4	6	0.0414	100	-24.1546; 
8	5	7	0.116	100	-8.6207; 
9	6	7	0.082	100	-12.1951; 
10	6	8	0.042	120	-23.8095;
11	6	9	0.208	100	-4.8077; 
12	6	10	0.556	60	-1.7986; 
13	9	11	0.208	60	-4.8077; 
14	9	10	0.11	100	-9.0909; 
15	4	12	0.256	100	-3.9063; 
16	12	13	0.14	120	-7.1429; 
17	12	14	0.2559	100	-3.9078; 
18	12	15	0.1304	80	-7.6687; 
19	12	16	0.1987	80	-5.0327; 
20	14	15	0.1997	100	-5.0075; 
21	16	17	0.1932	60	-5.1760; 
22	15	18	0.2185	60	-4.5767; 
23	18	19	0.1292	60	-7.7399; 
24	19	20	0.068	100	-14.7059; 
25	10	20	0.209	100	-4.7847; 
26	10	17	0.0845	60	-11.8343; 
27	10	21	0.0749	80	-13.3511; 
28	10	22	0.1499	100	-6.6711; 
29	21	22	0.0236	120	-42.3729; 
30	15	23	0.202	100	-4.9505; 
31	22	24	0.179	60	-5.5866; 
32	23	24	0.27	100	-3.7037; 
33	24	25	0.3292	60	-3.0377; 
34	25	26	0.38	80	-2.6316; 
35	25	27	0.2087	100	-4.7916; 
36	28	27	0.396	80	-2.5253; 
37	27	29	0.4153	60	-2.4079; 
38	27	30	0.6027	100	-1.6592; 
39	29	30	0.4533	80	-2.2060; 
40	8	28	0.2	80	-5.0000; 
41	6	28	0.0599	120	-16.6945];
P=[1	99.2	0;
2	80	61.7;
3	0	72.4;
4	0	27.6;
5	0	0;
6	0	0;
7	0	52.8;
8	78	56;
9	0	0;
10	0	35.8;
11	0	14;
12	0	41.2;
13	50	0;
14	70	0;
15	0	88.2;
16	0	23.5;
17	0	9;
18	0	3.2;
19	0	0;
20	20	0;
21	0	57.5;
22	56	0;
23	42	0;
24	0	8.7;
25	0	23;
26	0	8.5;
27	28	0;
28	68	0;
29	0	25.4;
30	32	0];
Xl=[0.0575
0.1852
0.1737
0.0379
0.1983
0.1763
0.0414
0.116
0.082
0.042
0.208
0.556
0.208
0.11
0.256
0.14
0.2559
0.1304
0.1987
0.1997
0.1932
0.2185
0.1292
0.068
0.209
0.0845
0.0749
0.1499
0.0236
0.202
0.179
0.27
0.3292
0.38
0.2087
0.396
0.4153
0.6027
0.4533
0.2
0.0599];


B(30,30)=0;%定义节点导纳矩阵为一个30X30的零矩阵；
for i1=1:length(D)
    B(D(i1,2),D(i1,3))=D(i1,6); B(D(i1,3),D(i1,2))=D(i1,6);%定义导纳矩阵非对角元素即互导纳；
end
for i2=1:30%定义导纳矩阵对角元素；
    for i3=1:length(D)
        if i2==D(i3,2)|i2==D(i3,3)
            B(i2,i2)=B(i2,i2)-D(i3,6);
        end
    end
end
B(13,:)=[];B(:,13)=[];%消去导纳矩阵平衡节点所在的行和列；
P(13,:)=[];%消去平衡节点注入功率；
P1=P(:,2)-P(:,3);%节点注入功率减去负荷功率；
O=inv(B)*P1;%直流潮流模型算出各节点相角；
O=[O(1:12);0;O(13:29)];
Xm=diag([Xl]);%Xm为线路电抗对角矩阵；
for i4=1:length(D)%M为支路节点关系矩阵起始节点为1，终止节点为-1；
    M(i4,D(i4,2))=1;M(i4,D(i4,3))=-1;
end
Pl=inv(Xm)*M*O;