1 лр / Lab1

clear
clc
disp('Вариант 11, Назлуян Давид, МП-42');
variant = 11;
A = [1,0,0,0,0,0,1,0,1,1]
C1 = [0,1,0,1,1,0,1,0,1,1]
C2 = [0,1,1,0,0,0,0,1,0,1]


disp('Сгенерируем М- последовательности:');
M1 = Mfun(A,C1);
M2 = Mfun(A,C2);
Msumma = Msum(M1,M2, variant); %сумма с учетом варианта
figure % график трех последовательностей
subplot(3,1,1) %М1
plot(1:length(M1),M1);
hold on;
title('M1');
subplot(3,1,2)
plot(1:length(M2),M2); %М2
hold on;
title('M2');
subplot(3,1,3)
plot(1:length(Msumma),Msumma); % Сумма М1 и неМ2
hold on;
title('Сумма М1 + неМ2');
figure %все остальные графики с КФ и фильтрацией


%АКФ для М1 с проверкой
subplot(4,2,1);
disp('Построим АКФ для М1');
M1AKF = myCF(M1,M1);
x = 1:length(M1AKF);
plot(x,M1AKF);
hold on;
grid on;
title('АКФ для М1');
subplot(4,2,2);
plot(x,xcorr(M1,M1));
hold on;
grid on;
title('АКФ для М1 проверка');


%ВКФ для М1 и М2 с проверкой
subplot(4,2,3);
disp('Построим ВКФ для М1 и М2');
M1AKF = myCF(M1,M2);
x = 1:length(M1AKF);
plot(x,M1AKF);
hold on;
grid on;
title('ВКФ для М1 и М2');
subplot(4,2,4);
plot(x,xcorr(M2,M1));
hold on;
grid on;
title('ВКФ для М1 и М2 проверка');


%фильтрация М1 из Суммы
subplot(4,2,5);
disp('Построим график фильтрациии М1 из Суммы');
M1AKF = myCF(M1,Msumma);
x = 1:length(M1AKF);
plot(x,M1AKF);
hold on;
grid on;
title('график фильтрациии М1 из Суммы');
subplot(4,2,6);
M1outSum = xcorr(Msumma,M1);
x = 1: length(M1outSum);
plot(x,M1outSum);
hold on;
grid on;
title('проверка');


%фильтрация М2 из Суммы
subplot(4,2,7);
disp('Построим график фильтрациии М2 из Суммы');
M1AKF = myCF(M2,Msumma);
x = 1:length(M1AKF);
plot(x,M1AKF);
hold on;
grid on;
title('график фильтрациии М2 из Суммы');
subplot(4,2,8);
M1outSum = xcorr(Msumma,M2);
x = 1: length(M1outSum);
plot(x,M1outSum);
hold on;
grid on;
title('проверка');



%Задания 10-13
disp('А вот теперь обратимся в вопросам шумоподавления.');
%собственно, генерируем шум
Amplitude = variant; %амплитуда шума
Noise = Amplitude-2*Amplitude*rand(1,length(M1)); %это шум
M3 = M1 + Noise; %смешали
figure
subplot(3,1,1);
hold on;
grid on;
plot(M3); %показали на графике
title('M1+шум');

subplot(3,1,2);
hold on;
grid on;
M3filtered = myCF(M1,M3);
plot(M3filtered); %показали на графике фильтрацию М1 из шума
title('фильтрация М1 из шума');


Db=20*log10(abs(M3filtered/max(M3filtered)));
subplot(3,1,3);
hold on;
grid on;
plot(Db); %показали на графике фильтрацию М1 из шума в дБ
title('фильтрация М1 из шума в Дб');