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

clc;
clear all;

Y=load('W2_04.txt'); % исходный сигнал ЭКГ
B1=[1/5 1/5 1/5 1/5 1/5];
B2=[1/8 1/4 1/8 0 -1/8 -1/4 -1/8];
M1=length(B1); % Кол-во отсчетов (длина)
M2=length(B2);

for n=1:length(Y) %n=1:250
X1(n)=Y(n,1); %Создаем массивы X1 и X2, пристваивая соот. зн-я из матрицы-столбца Y
X2(n)=Y(n,1);
end

for n=M1:length(Y) %Фильтруем исходный сигнал ЭКГ (1 фильтр) и записываем в массив Y1
Y1(n)= 0.2*X1(n)+0.2*X1(n-1)+0.2*X1(n-2)+0.2*X1(n-3)+0.2*X1(n-4);
end

for n=M2:length(Y) %Фильтруем исходный сигнал ЭКГ (2 фильтр) и записываем в массив Y2
Y2(n)=1/8*X2(n)+1/4*X2(n-1)+1/8*X2(n-2)-1/8*X2(n-4)-1/4*X2(n-5)-1/8*X2(n-6);
end

figure
Fs=250;
T=1/Fs;
tmax=length(Y)*T;
t=0:T:tmax-T;
subplot(3,1,1)
plot(t,Y)
title('ЭКГ сигнал')
ylim([-100,50])
xlim([0,tmax])

subplot(3,1,2)
plot(t,Y1)
title('Обработка сигнала ЭКГ фильтром 1')
ylim([-100 50])
xlim([0,tmax])

subplot(3,1,3)
plot(t,Y2);
title('Обработка сигнала ЭКГ фильтром 2')
ylim([-100 50])
xlim([0,tmax])