function?GraphTheory（APGPn_branchn_genn_busFee）
%%?生成節(jié)點?進(jìn)出節(jié)點集
for?i?=?1:n_branch
????if?AP（i3）<0
????????a?=?AP（i3）;
????????AP（i3）?=?AP（i4）;
????????AP（i4）=?a;
????????
????????b?=?AP（i1）;
????????AP（i1）?=?AP（i2）;
????????AP（i2）=?b;
????????
????end
end

AP?=?abs（AP）;
for?i?=?1:n_branch
????AP（i?4）?=?abs（AP（i?4））;
end

for?i?=?1:n_bus
????Node（i）.inNodeSET?=?[];
????Node（i）.outNodeSET?=?[];
end

for??i?=?1:n_branch
????for?j?=?1:n_branch
????????if??AP（j2）?==?i
????????????Node（i）.inNodeSET?=?[?Node（i）.inNodeSET?AP（j1）];
????????????
????????end
????end
end

for??i?=?1:n_branch
????for?j?=?1:n_branch
????????if??AP（j1）?==?i
????????????Node（i）.outNodeSET?=?[?Node（i）.outNodeSET?AP（j2）];
????????end
????end
end

%%?屏顯節(jié)點的進(jìn)出節(jié)點集
for?i?=?1:n_bus
????fprintf（‘%d‘i）;
????fprintf（‘號節(jié)點的進(jìn)出節(jié)點集為：\n‘）;
????disp（Node（i））;
end
%%?main?循環(huán)
result?=?zeros（n_genn_branch）;
Result?=?zeros（n_genn_branch）;
for?j?=?1?:?n_bus-1
????%?找出0進(jìn)線節(jié)點
????for?i?=?1?:?n_bus
????????if（?isempty（Node（i）.inNodeSET））
????????????nowNode?=?i;
????????????break
????????end
????end
????
????temp?=?sum?（GP（:nowNode））;
????temp_result?=?GP（:nowNode）/sum?（GP（:nowNode））;
????
????temp1?=?Node（nowNode）.outNodeSET;
????for?ii?=?1:?length（temp1）
????????for?ii1?=?1:n_branch
????????????if?AP（ii11）==?nowNode?&&?AP（ii12）?==?temp1（ii）
????????????????result（:ii1）?=?temp_result;
????????????end
????????end
????end
????
????%?temp1（ii）?為?nowNode的出節(jié)點
????for?ii?=?1:?length（temp1）
????????for?ii1?=?1:n_branch
????????????if?AP（ii11）==?nowNode?&&?AP（ii1?2）?==?temp1（ii）
????????????????GP（:temp1（ii））?=?GP（:temp1（ii））?+?GP（:nowNode）*（AP（ii14）/temp）;
????????????end
????????end
????????%?從?nowNode的出節(jié)點的進(jìn)線集中刪除當(dāng)前節(jié)點
????????if?length（Node（temp1（ii））.inNodeSET）?==?1
????????????Node（temp1（ii））.inNodeSET?=?[];
????????else
????????????for?ii2?=?1:?length（Node（temp1（ii））.inNodeSET）
????????????????if?Node（temp1（ii））.inNodeSET（ii2）?==?nowNode
????????????????????Node（temp1（ii））.inNodeSET（ii2）=?[];
????????????????????break
????????????????end
????????????end
????????end
????end
????%?將當(dāng)前節(jié)點?排除
????Node（nowNode）.inNodeSET?=?[-1];
end
for?i=1:n_branch
????Result（:i）=result（:i）*Fee（1i）;
end
disp（‘使用圖論方法求得的利用份額結(jié)果為：‘）;
disp（result）;
disp（‘使用圖論方法求得的輸電費用分?jǐn)偨Y(jié)果為：‘）
disp（Result）;

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2486??2015-02-13?17:42??PowerTrace\GraphTheory.m

?????文件?????????98??2015-02-13?17:42??PowerTrace\mainPT.m

?????文件????????391??2015-02-13?17:41??PowerTrace\Pbus5.m

?????文件??????55644??2015-02-13?18:10??PowerTrace\潮流追蹤結(jié)果.png

?????目錄??????????0??2015-03-10?20:21??PowerTrace

-----------?---------??----------?-----??----

????????????????58619????????????????????5