function [w1]=optwbm(bWin,n,w,y,beta,mu,snrThresh,amcInds)
% adapt everyone under no threshold
if isempty(snrThresh)
w = beta*w + mu*bWin*ones(1,size(y,2))*diag(y(n,:)) / …
norm(bWin)^2;
% threshold certain tracks
else
w(:,amcInds) = beta*w(:, amcInds) + mu*bWin*…
ones(1,length(find(amcInds))) * diag(y(n,amcInds)) / …
norm(bWin)^2;
end
w1=w;
end