function datacheck(putin,n) % putin 为输入信号,n为阶次
[m1,n1]=size(putin);
if m1<n1
putin=putin';
end
l=length(putin);
u=putin; % 输入数据矩阵u为l*1
F=diag(ones(l-1,1),-1);%l*l的对角阵,F 为位移矩阵
Ul=zeros(l,2*n);
for i=1:1:2*n
Ul(:,i)=(F^i)*u;
end
result=Ul'*Ul/l
eig(result)
det(result)