希尔伯特变换.txt

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;
subplot(311);plot(t,x)
title('带噪声的正弦信号');
subplot(312);plot(t,s);
title('正弦信号');
subplot(313);plot(t,e);
title('滤波结果')





clear all;  
clc;  
close all;  
  
fs = 44100;  
% x = wavread('b.wav');  
t = -5*pi:pi/100:5*pi;