Лабы_Вариант 4_5 сем / ЛР_5 / filt_program

global hAxes1 hAxes2 hAxes3 hEda

Y=load('W5_04.txt');
LY=length(Y);
Fs=250;
T=1/Fs;
tmax=LY*T;
t=0:T:tmax-T;
N=length(Y);
axes(hAxes3)
cla('reset');
hold on
set (hAxes3,'YLim',[-100 160])
plot (t,Y)
title ('Тестовая синусоида')

x = Y.';

axes(hAxes2)
str=get(hEda,'String');
da=str2double(str);
y=fs50 (x,Fs,da);
cla('reset');
hold on
set (hAxes2,'YLim',[-100 160])
plot (t,y)
title ('Сигнал после фильтрации')

y2=zeros(1,N);
axes(hAxes1)
for j=1:N
y2(1,j)=x(1,j)-y(1,j);
end
cla('reset');
hold on
set (hAxes1,'YLim',[-100 160])
plot (t,y2)
title ('График компенсирующего сигнала')