function?[delPQ]?=deltaPQ（Ybus）?????????%計(jì)算功率偏差
G=real（Y）;
B=imag（Y）;
[nymy]=size（Y）;
[nbmb]=size（bus）;
nSW=0;????????????????????????%計(jì)數(shù)：平衡節(jié)點(diǎn)個(gè)數(shù)
nPV=0;????????????????????????%計(jì)數(shù)：PV節(jié)點(diǎn)個(gè)數(shù)
nPQ=0;????????????????????????%計(jì)數(shù)：PQ節(jié)點(diǎn)個(gè)數(shù)
for?I=1:nb????????????????????%?nb為總節(jié)點(diǎn)數(shù)
????type=bus（I6）;
????if?type==3
???????nSW=nSW+1;?????????????%?increment?swing?bus?counter
????elseif?type==2
???????nPV=nPV?+1;????????????%?increment?PV?bus?counter
????else
???????nPQ=nPQ+1;?????????????%?increment?PQ?bus?counter
????end
end
ang=bus（:3）;?????????????????%把bus第3列相角存放到矩陣ang中
V=bus（:2）;???????????????????%把bus第2列相角存放到矩陣V中
A=zeros（nymy）;???????????????%矩陣A存放節(jié)點(diǎn)相角差
P=bus（:4）;???????????????????%把bus第4列輸入有功功率存放到矩陣P中
Q=bus（:5）;???????????????????%把bus第5列輸入無功功率存放到矩陣Q中
for?k=1:ny????????????????????%求各節(jié)點(diǎn)相角差存放到矩陣A中
????A（k:）=ang（k）-ang（:）;
end
for?j=1:（nPQ+nPV）?????????????%計(jì)算PQ、PV節(jié)點(diǎn)的有功功率偏差dP存放到矩陣delP中
????delP（j）=P（j）;
????for?k=1:ny
????????delP（j）=delP（j）-V（j）*V（k）*（G（jk）*cos（A（jk））+B（jk）*sin（A（jk）））;
????end
end
for?j=1:nPQ??????????????????%計(jì)算PQ節(jié)點(diǎn)的無功功率偏差dQ存放到矩陣delQ中
????delQ（j）=Q（j）;
????for?k=1:ny
????????delQ（j）=delQ（j）-V（j）*V（k）*（G（jk）*sin（A（jk））-B（jk）*cos（A（jk）））;
????end
end
delPQ=[delPdelQ]‘;??????????%按dP、dQ節(jié)點(diǎn)的順序形成一個(gè)單列矩陣存放到矩陣delPQ中

????
????
????
????