angles.m

电力系统分析计算程序

function [s12,c12] = angles(a)

global DAE

t1 = DAE.y(a.bus1);
t2 = DAE.y(a.bus2);

s12 = sin(t1-t2);
c12 = cos(t1-t2);