MATLAB直流潮流计算

%%-----  run the power flow  -----                                          %%%%%%%%%%%%%%%%直流潮流计算固定块。
%% initial state
Va0 = bus（: VA） * （pi/180）;                                               %%转换为弧度。
%% 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里面计算母线的复功率S；取real后相当于有功功率P。
%% “run“ the power flow   
Va = dcpf（B Pbus Va0 ref pv pq）;                                      %%直