clear all, close all, clc % Rimozione di variabili nel workspace, chiusura figure, cancello command Window

N=3e4; % Variabile Scalare
MALength = 1e3; % far vedere 1e3, e poi l'effetto di 1e4, che taglia via anche la parte sinusoidale
n=1:N; % Vettore riga da 1 a N, a intervallo di 1 (1,2,3,4,...,N)
Signal = 0.15*cos(2*pi*n/7000)-n/20000; % Cosinusoide a frequenza 1/7000 meno una retta
x=Signal+randn(1,N); % Aggiungo un segnale casuale di disturbo
plot((1:N-MALength),x(MALength/2 + 1:N-MALength/2)), hold on, grid on % Disegno il segnale, opportunamente ritardato
plot((1:N-MALength),Signal(MALength/2 + 1:N-MALength/2), 'g') % Disegno il segnale rumoroso, opportunamente ritardato



h=ones(1,MALength)/MALength; % Filtro via il rumore
y=conv(x,h); % Convoluzione dei due segnali
pause
plot((1:N-MALength),y(MALength + 1:N),'r') % Disegno il risultato, più ritardato ancora
% per tenere conto dell'effetto del ritardo dovuto al filtro CAUSALE
pause
axis([1, N-MALength, min(Signal)-0.15, max(Signal)+0.15]) % Cambio il punto di osservazione degli assi