base.m

这是一个很适合研究和学习用的电力系统仿真软件

function p = base(p)

global Line Bus Settings

if ~p.n, return, end

p.B = p.Cp./Line.con(p.line,9);
p.y = (1-p.Cp)./Line.con(p.line,9);

fm_errv(p.con(:,6),'Tcsc',p.bus1)

Vb2old = p.con(:,6).*p.con(:,6);
Vb2new = Bus.con(p.bus1,2).*Bus.con(p.bus1,2);

k = Settings.mva*Vb2old./p.con(:,5)./Vb2new;

if p.ty1
  p.con(p.ty1,10) = k(p.ty1).*p.con(p.ty1,10);
  p.con(p.ty1,11) = k(p.ty1).*p.con(p.ty1,11);
end
if p.ty2
  p.con(p.ty2,14) = k(p.ty2).*p.con(p.ty2,14);
  p.con(p.ty2,15) = k(p.ty2).*p.con(p.ty2,15);
end

p.store = p.con;