%%%%%

%data filtering

Fs=100;
t=[0:Fs]/Fs;
s=10+sin(2*pi*t*5)+sin(2*pi*t*10);
figure;
plot(s);

[b,a]=ellip(8,0.01,40,6/Fs);
[h,w]=freqz(b,a,512*10);
figure;
plot(w*Fs/(2*pi),abs(h));

sf=filter(b,a,s);

figure;

plot(t,sf);