fourth_lab

% b= [0.2 0.3 0.5 0.7];
% a = 1;
% n = 0 : 64;
% x = 0.1 * sin (0.7 * n);
% y = filter (b,a,x);
% figure;
% plot (n, x,'LineWidth', 2);
% hold on;
% plot (n, y);
% grid on;
% xlim([0 64]);
% title("КИХ - фильтр");
% xlabel ('Частота');
% ylabel ('Амплитуда');
% legend('входной сигнал', 'выходной сигнал');
% hold off;

% b = [1 1 1 1];
% a = [1 -0.3 0.25 0.4];
% n = 0 : 64;
% x = 0.1 * sin (0.7 * n);
% y=filter (b,a,x);
% figure;
% plot (n, x,'LineWidth', 2);
% hold on;
% plot (n, y);
% grid on;
% xlim([0 64]);
% title("БИХ - фильтр");
% xlabel ('Частота')
% ylabel ('Амплитуда')
% legend('входной сигнал', 'выходной сигнал')
% hold off;

b= [1 1 1 1];
a = [1 -0.3 0.25 0.4];
n = 0 : 64;
y = dirac(n);
idx = y == Inf;
y(idx) = 1;
y=filter (b,a,y);
impz(y), grid
title("Импульсная характеристика БИХ - фильтра");
xlabel ('Частота')
ylabel ('Амплитуда')

b= [1 1 1 1];
a = [1 -0.3 0.25 0.4];
impz(b, a, 80, 5000), grid on;
title("Импульсная характеристика");
xlabel ('N')
ylabel ('импульсная характеристика БИХ-фильтра (2)')