second_lab

X = zeros (1,10000);
Y = zeros (1,10000);
X(1:100) = 5;
Y(301:400) = 1;


% conv_R = conv(X,Y);
% conv_R = conv_R(1:1000);
% plot(conv_R);
% title('Автокорреляция одного импульса');
% xlabel('Время t');
% ylabel('Значение');
%
% X(3500:3600)=5;
% Y(3801:4000)=1;

% xcorr_R = xcorr(X);
% figure;
% plot(xcorr_R);
% xlim([6000 14000]);
% ylim([0 5200]);
% xticks (6000:1000:14000);
% title('Корреляция двух импульсов');
% xlabel('Время t');
% ylabel('Значение');

t = 0:0.1:100;
s=3;
p=0;
for i = 1:10
p = p + rectpuls(t-s, 3);
s = s + 5 + 10 * rand;
end
% figure;
% plot(t,p),grid, set(gca,'FontSize',12);
% title('Сигнал без периодичных повторов импульсов');
% xlabel('Время t');
% ylabel('U');

delayed_pulse=zeros(1,300);
delayed_pulse=[delayed_pulse,p];

xcorr_p = xcorr(p,delayed_pulse);
figure;
plot(xcorr_p);
xlim([0 2000]);
ylim([-20 320]);
title('Корреляция');
xlabel('Время t');
ylabel('Значение');

xcov_p = xcov(p,delayed_pulse);
figure;
plot(xcov_p);
xlim([0 2400]);
ylim([-100 220]);
xticks (0:400:2400);
title('Ковариация');
xlabel('Время t');
ylabel('Значение');

% figure;
% out = awgn(delayed_pulse,1.5);
% plot(delayed_pulse, 'LineWidth',3);
% hold on;
% plot(delayed_pulse+out);
% legend('Меандр','Меандр с шумом');
% title('Сигнал без периодичных повторов импульсов c шумом');
% xlabel('Время t');
% ylabel('U');
% hold off;
%
% f = 2;
% t = 1:0.01:5;
% sin_x = sin(2*pi*f*t);

% figure
% hold on
% plot(t,sin_x, 'LineWidth',3);
% sin_x_noise = awgn(sin_x, 10);
% plot(t,(sin_x+sin_x_noise));
% legend('Синусоида','Синусойда с шумом');
% title('Синусоида с шумом');
% xlabel('Время t');
% ylabel('U');
% hold off;

% t=1:0.01:9;
% xcorr_R3=xcorr(sin_x,(sin_x+sin_x_noise));
% figure
% plot(t,xcorr_R3),grid;
% title('Корреляция');
% xlabel('Время t');
% ylabel('Значение');

% barker = comm.BarkerCode('Length',13,'SamplesPerFrame',13);
% Bar = barker();
% Bar = Bar';
% t = 0:(length(Bar)*13) - 1;
% s0 = square(t);
% s(1:169)=0;
% for j=1:13
% for i=1:13
% s(:,i*j) = s0(:,i*j)*Bar(j);
% end
% end
% figure
% hold on
% plot(t, s,'red', 'LineWidth',1.5);
% w_noise = awgn(s, 12.5);
% plot(t,w_noise,'blue');
% legend('Код Баркера','Код Баркера с шумом');
% title('Код Баркера');
% xlabel('Время t');
% ylabel('U');
% hold off

% xcorr_R4=xcorr(Bar);
% plot(xcorr_R4), grid;
% title('Автокорреляция');
% xlabel('Время t');
% ylabel('Значение');

%
% K=0:1:13;
% X=randi([-1,1],1,14);
% w_noise2 = awgn(X,11);
% figure
% hold on
% plot(K,X,'g'), grid;
% plot(K,w_noise2,'r'), grid;
% legend('Шумоподобный сигнал','Шумоподобный сигнал с шумом');
% xlabel('Время t');
% ylabel('U');
% hold off;
%
% R5=xcorr(X);
% plot(R5), grid;
% title('Автокорреляция');
% xlabel('Время t');
% ylabel('Значение');

% t =0:1:10000;
% fn = 100; % несущая частота
% fc = 0.5; % частота сигнала
% m = 0.5; % индекс модуляции
%
%
% mono = sin(2*pi*fc*t);
% mat1 = mean(mono);
% dispers1 = var(mono);
% dispers12 = dispers1^(1/2);
% % histogram(mono);xlabel('Моногармонический');grid on;
%
% lin = sin(2*pi*fc*t.^2);
% mat2 = mean(lin);
% dispers2 = var(lin);
% dispers22 = dispers2^(1/2);
% % histogram(lin);xlabel('С линейной частотой модуляцией');grid on;
%
% ampl = sin(2*pi*fn*t).*(1+m*sin(2*pi*fc*t));
% mat3 = mean(ampl);
% dispers3 = var(ampl);
% dispers32 = dispers3^(1/2);
% % histogram(ampl);xlabel('С амплитудной модуляцией');grid on;
%
% chastot = sin(2*pi*fc*(cos(0.5*t)).*t);
% mat4 = mean(chastot);
% dispers4 = var(chastot);
% dispers44 = dispers4^(1/2);
% % histogram(chastot);xlabel('С частотной модуляцией');grid on;
%
% phase = sin(2*pi*fn*t+5*sin(2*pi*t));
% mat5 = mean(phase);
% dispers5 = var(phase);
% dispers55 = dispers5^(1/2);
% % histogram(phase);xlabel('С фазовой модуляцией');grid on;