clear
clc
>> t=0:1/10000:1-0.0001;
>> s=sin(2*pi*t);
>> n=randn(size(t));
>> x=s+n;
>> w=[0,0.5];
>> u=0.00026;
>> for i=1:9999;
y(i+1)=n(i:i+1)*w';
e(i+1)=x(i+1)-y(i+1);
w=w+2*u*e(i+1)*n(i:i+1);
end
>> figure(1)
>> subplot(3,1,1)
>> plot(t,x);
>> title('帶噪信號')
>> subplot(3,1,2)
>> plot(t,s);
>> title('正弦信號')
>> subplot(3,1,3)
>> plot(t,e)
>> title('濾波結果')