pmech.m

电力系统分析计算程序

function pm = pmech(p)

global DAE

pm = zeros(p.n,1);

if p.ty1
  pm(p.ty1) =  DAE.x(p.dat1(:,3)) + p.dat1(:,12).* ...
      (DAE.x(p.dat1(:,2)) + p.dat1(:,10).*DAE.x(p.dat1(:,1)));
end

if p.ty2
  pm(p.ty2) = DAE.x(p.dat2(:,1)) - p.dat2(:,7).* ...
      (DAE.x(p.dat2(:,9)) - p.con(p.ty2,3)) + p.dat2(:,5);
  pm(p.ty2) = max(pm(p.ty2),p.dat2(:,4));
  pm(p.ty2) = min(pm(p.ty2),p.dat2(:,3));
end