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

