该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
function [scores bus ] = runpdwga(pop, baseMVA, bus,Main, M, gen, NN, branch, mpopt)
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
baseMVA=1
bus = [
1 3 0 0 0 0 1 1 0 10 1 1.1 0.9;
2 1 0.0831 0.0304 0 0 1 1 0 10 1 1.1 0.9;
3 1 0 0 0 0 1 1 0 10 1 1.1 0.9;
4 1 1.7451 0.6392 0 0 1 1 0 10 1 1.1 0.9;
5 1 0 0 0 0 1 1 0 10 1 1.1 0.9;
6 1 2.9294 1.0729 0 0 1 1 0 10 1 1.1 0.9;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio angle status
Main=[
1 2 0.00084 0.001155 0 250 250 250 0 0 1;
1 6 0.00392 0.00539 0 250 250 250 0 0 1;
2 3 0.0021 0.002888 0 250 250 250 0 0 1;
]
M = [
3 4 0.0007 0.0009625 0 150 150 150 0 0 1;
4 5 0.00042 0.0005775 0 150 150 150 0 0 1;
6 5 0.00028 0.000385 0