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）;%第一個信號
sig2=exp（-j*0.3*2*pi*（0:N-1）+j*2*pi*rand）;%第二個信號
sig3=exp（j*0.16*2*pi*（0:N-1）+j*2*pi*rand）;%第三個信號
Un=sig1+sig2+sig3+noise;%帶噪聲的信號
un=[zeros（1M-1）Un].‘;%擴展數(shù)據(jù)
A=zeros（MN）;%構建樣本矩陣
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）;%進一步預測信號
w=zeros（MN）;
epsilon=zeros（N-11）;%先驗估計誤差
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（: