Лаба 2 (Белый шум и случайные блуждания) / White_Gaussian_noise

clc
clear
close all

%% Task №1
matrix = zeros(700,400);
for i = 1:700
for j = 1:400
matrix(i,j) = normrnd(14,8);
end
end

%average of each columns ( time average )
s = zeros(1,400);
for i = 1:400
s(1,i) = mean(matrix(:,i));
end
%average of each rows ( ensemble average )
sk = zeros(1,700);
for i = 1:700
sk(1,i) = mean(matrix(i,:));
end

figure('Name','Ensemble and time average','NumberTitle','off','Position',[500,150,1000,700])
plot(s);
hold on
plot(sk);
hold on
legend('time average',' ensemble average')
xlim([0 400]);

clear sk s
%% Task №2
figure('Name','Scatter plot','NumberTitle','off','Position',[500,150,1000,700])
l = randi(400,3,2);
r(1,:) = xcorr(matrix(:,l(1,1)),matrix(:,l(1,2)),1,'normalized');
r(2,:) = xcorr(matrix(:,l(2,1)),matrix(:,l(2,2)),1,'normalized');
r(3,:) = xcorr(matrix(:,l(3,1)),matrix(:,l(3,2)),1,'normalized');
coeff = table(r(:,2)','VariableNames',{'r-coeff'});

a(1) = scatter(matrix(:,l(1,1)),matrix(:,l(1,2)),20,'r');
hold on
set(a(1),'Visible','off');
a(2) = scatter(matrix(:,l(2,1)),matrix(:,l(2,2)),20,'b');
hold on
set(a(2),'Visible','off');
a(3) = scatter(matrix(:,l(3,1)),matrix(:,l(3,2)),20,'g');
hold on
set(a(3),'Visible','off');
xlim([-20 50])
ylim([-20 50])
lgd = legend(num2str(l(1,:)),num2str(l(2,:)),num2str(l(3,:)));
title('Scatter plot for random realisations')
title(lgd,'Order')


for k = 1:3
cbx(k) = uicontrol('Style','checkbox', 'string', r(k,2) , 'Position',[20 70*k 70 15], 'Callback',{@checkBoxCallback,k,a});
end

clear a matrix vectorx vectory k i j linspace r cbx lgd l