%%-----??run?the?power?flow??-----??????????????????????????????????????????%%%%%%%%%%%%%%%%直流潮流計算固定塊。
%%?initial?state
Va0?=?bus（:?VA）?*?（pi/180）;???????????????????????????????????????????????%%轉(zhuǎn)換為弧度。
%%?build?B?matrices?and?phase?shift?injections
[B?Bf?Pbusinj?Pfinj]?=?makeBdc（baseMVA?bus?branch）;???????????????????%%計算出B矩陣?從而利用makesbus計算出輸入功率的電壓角度和末端功率的輸出角度
%%?compute?complex?bus?power?injections?（generation?-?load）
Pbus?=?real（makeSbus（baseMVA?bus?gen））?-?Pbusinj?-?bus（:?GS）?/?baseMVA;?%%real里面計算母線的復(fù)功率S；取real后相當于有功功率P。
%%?“run“?the?power?flow???
Va?=?dcpf（B?Pbus?Va0?ref?pv?pq）;??????????????????????????????????????%%直