xhat = x + 0.1*randn(1,ntrW);
L = 350;
muW = 0.000001;
h = adaptfilt.filtxlms(L,muS,1,Hhat);
[y,e] = filter(h,xhat,d);

%plot(n,d,'b',n,r,'g',n,e,'r');
plot(n,d,'b',n,y,'g',n,e,'r');
xlabel('Number of iterations');
ylabel('Signal value');
title('Active Noise Control Using the Filtered-X LMS Adaptive Controller');
legend('Original Noise','Anti-Noise','Residual Noise');
sound(e/max(abs(e)),Fs);