harj35.m

% Suodatetaan kohinainen signaali

x = load( 'data35.txt' );
hold off
plot( x(:,1), x(:,2), 'b-' )

% Suodatus
dt = x(2,1)-x(1,1);
fc = 10000;
rc = 1/(2*pi*fc);
alpha = dt/(rc+dt);

y(1) = x(1,2);
for i=2:size(x,1)
    y(i) = alpha*x(i,2) + (1-alpha)*y(i-1);
end
hold on;
plot( x(:,1), y, 'r-' )