function [h,e]=WH(Rss,Rww,M)
e1=10;e0=0;N=0;
while abs(e0-e1)>1e-6
N=N+1;
e0=e1;
Rxs=Rss(M:M+N-1);
Rxx=Rww(M:M+N-1)+Rss(M:M+N-1);
R_xx=zeros(N);

for j=1:N

for n=1:N

R_xx(j,n)=Rxx(abs(j-n)+1);
end
end
h=inv(R_xx)*Rss';
e1=Rss(M)-h'*Rxs';
end,N
e=e1;