main1.m

可用传输容计算程序

function [P,PL1]=main1(D,P,Xl)%计算正常情况下的ATC；
PL1=zeros(length(D),1);
PLL=D(:,5);
area1=[1 2 5 7 28 3 4 6 8 9 11];
area2=[13 12 16 14 18 17 15 19 20 23];
area3=[10 27 22 21 30 24 25 29 26];
%例如计算区域1到区域3的ATC则增加区域1发电，增加区域3负荷；
while sum(abs(PL1)>PLL)==0&sum(P(:,2)>P(:,4))==0&sum(P(:,3)>P(:,5))==0
    P(area1,2)=P(area1,2)+(10/sum(P(area1,2)))* P(area1,2);%增加区域1发电;
    P(area3,3)=P(area3,3)+(10/sum(P(area3,3)))* P(area3,3);%增加区域3负荷；
    [O,PL]=chaoliu(D,P,Xl);
    PL1=PL;
end