Лаба 5 / filter_ekg

%ЭКГ
global hAxes1 hAxes2 hAxes3 hEda

T=1/Fs;
t=0:T:1-T;

%исходный
axes(hAxes1);
plot(t,x);
grid minor;
ylim([-120 120]);

str=get(hEda,'String');
da=str2num(str);

y=fs50(x,Fs,da);

%отфильтрованный
axes(hAxes2);
plot(t,y);
grid minor;
ylim([-120 120]);

%компенсирующий
for i=1:length(x)
a(i)=x(i)-y(i);
end
axes(hAxes3);
plot(t,a);
grid minor;
ylim([-120 120]);