新的信号处理方法-TLS-ESPRIT

kevin19821

clear
w1=pi/4;
w2=-pi/3;
w3=pi/3;
w4=-pi/4;
kp=2048;
a1=[exp(-j*pi*sin(w1)*[0:9])].';
a2=[exp(-j*pi*sin(w2)*[0:9])].';
a3=[exp(-j*pi*sin(w3)*[0:9])].';
a4=[exp(-j*pi*sin(w4)*[0:9])].';
A=[a1,a2,a3,a4];
k=1:kp+1;
s=[1.3*cos(k*0.015);sin(k*0.05);cos(k*0.02);sin(k*0.035)];
n=randn(10,kp+1)+j*randn(10,kp+1);
xa=A*s+n;
x=xa(:,1:kp);
y=xa(:,2:kp+1);
r=xa*(xa)'/(kp+1);
rxx=x*(x)'/kp;
rxy=x*(y)'/kp;
[V,D]=eig(r);
[b,index]=sort(diag(D));
lam=b(1);
z=diag(ones(1,9),-1);
cxx=rxx-lam*eye(10);
cxy=rxy-lam*z;
[U,S,Vx]=svd(cxx);
s1=S(1:4,1:4);
u1=U(:,1:4);
v1=(Vx(1:4,:))';
d=eig(s1,(u1)'*cxy*v1);
ww=angle(d);
theta=asin(-ww/pi);
t=theta*180/pi


大家帮看看 程序错在那里了?自己找不出来