Nb = 7;   % НОМЕР БРИГАДЫ
b = [0.64, 0.71, 0.77];     % ВЕКТОР КОЭФФИЦИЕНТОВ ЧИСЛИТЕЛЯ ПЕРЕДАТОЧНОЙ ФУНКЦИИ
a = [1, -0.95, 0.68];     % ВЕКТОР КОЭФФИЦИЕНТОВ ЗНАМЕНАТЕЛЯ ПЕРЕДАТОЧНОЙ ФУНКЦИИ
eps = 10^(-3);
N2 = 1000;   % ДЛИНА ВОЗДЕЙСТВИЯ
Fs = 7000;   % ЧАСТОТА ДИСКРЕТИЗАЦИИ

disp('% п.1. ВЫЧИСЛЕНИЕ ИМПУЛЬСНОЙ ХАРАКТЕРИСТИКИ - функция impz')
disp('%')
disp('%')

disp('% п.2. ВЫЧИСЛЕНИЕ ИМПУЛЬСНОЙ ХАРАКТЕРИСТИКИ - функция filter')
disp('%')
disp('%')
disp('% Для вывода ГРАФИКА ИМПУЛЬСНОЙ ХАРАКТЕРИСТИКИ нажмите <ENTER>')
pause
n = 0:(N2-1);
u0 = [1 zeros(1,(N2-1))];    % ЦИФРОВОЙ ЕДИНИЧНЫЙ ИМПУЛЬС
h2 = filter(b,a,u0);         % ИМПУЛЬСНАЯ ХАРАКТЕРИСТИКА
subplot(2,1,2), stem(n,h2,'fill','MarkerSize',3), grid
xlabel('n'), ylabel('h(n)'), title('Impulse Response h(n) - filter')
disp('%')
disp('%')
disp('% Для продолжения нажмите <ENTER>')
pause
disp('%')
disp('%')
disp('% п.3. ВЫЧИСЛЕНИЕ РЕАКЦИИ ПО ФОРМУЛЕ СВЕРТКИ')
disp('%')
disp('%')
disp('% Для вывода ГРАФИКОВ ВОЗДЕЙСТВИЯ И РЕАКЦИИ, вычисленной по ФОРМУЛЕ СВЕРТКИ,  нажмите <ENTER>')
pause
h1 = eps + h2;
x = input_1(N2);    % ВОЗДЕЙСТВИЕ (ДИСКРЕТНЫЙ ПРЯМОУГОЛЬНЫЙ ИМПУЛЬС)
y1 = conv(x,h1);    % РЕАКЦИЯ ДЛИНЫ, РАВНОЙ ДЛИНЕ СВЕРТКИ

ls = length(y1); %длина свёртки
N1 = -N2+1+ls