% This is the function program for state space (ESPRIT, Subspace-Rotation)
% the input is the Us and Es
function ssrts=ssprit(Us,Es)
[L,P]=size(Us);
Uu=Us(1:L-1,:);
Ud=Us(2:L,:);
Psi=Uu\Ud; Phi=eig(Psi); ssrts=sort(asin(angle(Phi)/pi))';

