use_doa.m

关于高阶累积量分析的Matlab函数

clf,
clear,
nfft=256;
t=0:1/nfft:1;
signal=[exp(sqrt(-1)*2000*2*pi*t);...
        exp(sqrt(-1)*1500*2*pi*t);...
        exp(sqrt(-1)*3000*2*pi*t)];  
[p,q]=size(signal); 
st=signal;
theta=[10,30,38]*pi/180;
steer=exp(sqrt(-1)*pi*sin(theta));
A=[steer.^(-5);steer.^(-4);steer.^(-3);steer.^(-2);steer.^(-1);steer.^0;...
        steer;steer.^2;steer.^3;steer.^4;steer.^5];
val=1/sqrt(2);
nt=(1/sqrt(2))*(val*randn(size(A*st))+j*val*randn(size(A*st)));
snr=80;
pow=10.^(-snr/10);
xt=A*st;%+pow*nt;
doa(xt',0.5,2.,3,2,1);