新的信号处理方法-TLS-ESPRIT
clearw1=pi/4;
w2=-pi/3;
w3=pi/3;
w4=-pi/4;
kp=2048;
a1=)].';
a2=)].';
a3=)].';
a4=)].';
A=;
k=1:kp+1;
s=;
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;
=eig(r);
=sort(diag(D));
lam=b(1);
z=diag(ones(1,9),-1);
cxx=rxx-lam*eye(10);
cxy=rxy-lam*z;
=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
大家帮看看 程序错在那里了?自己找不出来
页:
[1]