其实在costfunc里有各种的乘以baseMVA,这儿一个 那个一个
比如:
xx = [ Pg; Qg ] * baseMVA;
f = sum( totcost(gencost(ipol, :), xx(ipol)) ); %% cost of poly P or Q
df_dPgQg(ipol) = baseMVA * polycost(gencost(ipol, :), xx(ipol), 1);
其实在最后都乘以了一个系数:
f = f * opt.cost_mult;
df = df * opt.cost_mult;
其中cost_mult = 1e-4 又回到了原来的值