?? rootmusic.m
字號:
clear all
close all
A=1/10;
N =10;
n = (-(N-1)/2:(N-1)/2)';
Vm = ones(N,1);
u=1:100;
SNR=10.^(20/10);
ui=0.15;
us=0;
Vi = exp(j*n*pi*ui);
Vs=exp(j*n*pi*us);
Sx=zeros(10,10)
for i=1:100
X=sqrt(SNR)*exp(j*2*pi*randn)*(randn+j*randn)*Vi/sqrt(2)+exp(j*2*pi*randn)*(randn+j*randn)'*sqrt(SNR)*Vs/sqrt(2)+(randn+j*randn)/sqrt(2)
Sx=(X*X')+Sx
end
Sx=Sx/100
[U,S,V]=svd(Sx)
G=[];
for k=3:10
G=[G,V(:,k)];
end
GG=G*G'
a=zeros(20,1)
for k=1:10
for l=1:10
a(10-k+l)=a(10-k+l)+GG(k,l);
end
end
z=roots(a)
for k=1:19
if abs(z(k))>1.02|abs(z(k))<0.98
z(k)=0;
end
end
z=sort(z);
anglez=angle(z)/2/pi;
absz=abs(z)
root1=abs(anglez(19)*2)
root2=abs(anglez(18)*2)
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -