Matlab #16: Le produit de convolution 1D

clear all; close all; clc

%%

% Matlab #16: Le produit de Convolution 1D
% www.Electronique-Mixte.fr

%%  Signal bruité  

% Signal 
f0=1; t0=1/f0; N=1e4; n=2; 
t=linspace(0,n*t0,N);
A=1; s_n=A*sin(2*pi*f0*t);

% Bruit
sig=1; 
b_n=sig*randn(1,N);

% Signal + Bruit 
s_nb=s_n+b_n;


%% Le Produit de Convolution 1D

Nf=64;
H=ones(1,Nf)/Nf; 
y_n=zeros(1,N); 

for i=Nf:N
    y0=s_nb(i-Nf+1:i);  
    y_n(i)=sum(y0.*H);
end; 

%% Affichage 

figure; 
plot(t,s_nb,'k','linewidth',2);grid on;hold on;
plot(t,y_n,'r', 'linewidth',2);
plot(t,s_n,'linewidth',2);
legend({'Bruité','Filtré', 'Original'},'fontsize',16);