clear?all
N=256;?
n=0:1:N-1;????????????????
f1=0.2;??????????????????????????????
f2=0.213;
M=10;
snr=3;
?
x1=5*exp（1i*2*pi*f1*n+pi/3）+10*exp（1i*2*pi*f2*n+pi/4）;??
xn=awgn（x1snr）;
?
????M=8;?????????????????????????????????%濾波器階數
????for?k=1:N-M??????????????????????????
????xs（:k）=xn（k:1:k+M-1）.‘;
????end
????R=xs*xs‘/（N-M）;??????????????????????%估計樣本相關矩陣
???
?
????[UE]=svd（R）;??????????????????????