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

clc;
clear all;

% Построение коэффициенетов фильтров 1 и 2
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);

figure
subplot(4,2,1)
stem(B1)
title('Filter 1')
title('Коэффициенты 1-го фильтра')

subplot(4,2,2)
stem(B2)
title('Коэффициенты 2-го фильтра')

%Графики АЧХ
Fs=250;
T=1/Fs;
f=0:Fs/2;
H1=abs(1/5+2/5*cos(2*pi*f*T)+2/5*cos(4*pi*f*T));
subplot(4,2,3)
plot(f,H1)
title('АЧХ 1-го фильтра')

H2=abs(1/4*sin(2*pi*f*T)+1/2*sin(4*pi*f*T)+1/4*sin(6*pi*f*T));
subplot(4,2,4)
plot(f,H2)
title('АЧХ 2-го фильтра')

%Импульсные характеристики
N=20;
for n=1:N
X(n)=0;
end
X(10)=1;
Y1=zeros(1,N);
for n=M1:N
Y1(n)=0.2*X(n)+0.2*X(n-1)+0.2*X(n-2)+0.2*X(n-3)+0.2*X(n-4);
end
subplot(4,2,5)
stem(X,'.')
hold on
stem(Y1,'r')
title('Импульсная характеристика 1-го фильтра')

Y2=zeros(1,N);
for n=M2:N
Y2(n)=1/8*X(n)+1/4*X(n-1)+1/8*X(n-2)-1/8*X(n-4)-1/4*X(n-5)-1/8*X(n-6);
end
subplot(4,2,6)
stem(X,'.')
hold on
stem(Y2,'r')
title('Импульсная характеристика 2-го фильтра')

%Переходные характеристики
N=20;
for n=1:N/2
X(n)=0;
end
for n=N/2:N
X(n)=1;
end
Y1=zeros(1,N);
for n=M1:N
Y1(n)=0.2*X(n)+0.2*X(n-1)+0.2*X(n-2)+0.2*X(n-3)+0.2*X(n-4);
end
subplot(4,2,7)
stem(X,'.')
hold on
stem(Y1,'g')
title('Переходная характеристика 1-го фильтра')

Y2=zeros(1,N);
for n=M2:N
Y2(n)=1/8*X(n)+1/4*X(n-1)+1/8*X(n-2)-1/8*X(n-4)-1/4*X(n-5)-1/8*X(n-6);
end
subplot(4,2,8)
stem(X,'.')
hold on
stem(Y2,'g')
title('Переходная характеристика 2-го фильтра')