N=800;
M=16;
noise=0.005*（randn（1N）+j*randn（1N）/sqrt（2））;
sig1=exp（j*0.15*2*pi*（0:N-1）+j*2*pi*rand）;%第一個(gè)信號(hào)
sig2=exp（-j*0.3*2*pi*（0:N-1）+j*2*pi*rand）;%第二個(gè)信號(hào)
sig3=exp（j*0.16*2*pi*（0:N-1）+j*2*pi*rand）;%第三個(gè)信號(hào)
Un=sig1+sig2+sig3+noise;%帶噪聲的信號(hào)
un=[zeros（1M-1）Un].‘;%擴(kuò)展數(shù)據(jù)
A=zeros（MN）;%構(gòu)建樣本矩陣
for?n=1:N
????A（:n）=un（M+n-1:-1:n）;
end
delta=0.004;%調(diào)整參數(shù)
lambda=0.98;%遺忘因子
dn=Un（2:end）;%進(jìn)一步預(yù)測(cè)信號(hào)
w=zeros（MN）;
epsilon=zeros（N-11）;%先驗(yàn)估計(jì)誤差
PI=eye（M）/delta;
for?k=1:N-1%rls算法迭代
????PIn=PI*A（:k）;
????deno=lambda+A（:k）‘*PIn;
????kn=PIn/deno;
????epsilon（k）=dn（k）-w（: