clc;
clear all;
n = 350; % кількість значень вибірки
x_true = 2; % "правильне" значення вимірюваної величини
error = 0.26;% складова похибки
k = 2.3;
% створюємо вибірку

samples = x_true + randn(1,n) + error*randn(1,n);
istyne_value = mean(samples);
delta_for_x = samples - istyne_value;
sum_delta = sum(delta_for_x);%шум
noise_plus_k = k *(samples + delta_for_x);
ideal = samples * k;
plot(delta_for_x, 'b');hold on;
plot(ideal,'green');hold on;
plot(noise_plus_k,'r');hold off;
legend('Вхідний сигнал', 'Ідеальна мат.модель','З шумами');

experemental_k_ideal = ideal/ samples;
experemental_k_noise = noise_plus_k/samples

% ідеальна модель
istyne_value_ideal = mean(ideal);
delta_for_x_ideal = ideal - istyne_value_ideal;
sum_delta_ideal = sum(delta_for_x_ideal);
dysper_ideal = sum(delta_for_x_ideal.^2)/349;
std_ideal = std(ideal);
cigma_x_ideal = std_ideal/sqrt(350);

% для шуму
istyne_value_noise = mean(noise_plus_k);
delta_for_x_noise = noise_plus_k - istyne_value_noise;
sum_delta_noise = sum(delta_for_x_noise);
dysper_noise = sum(delta_for_x_noise.^2)/349;
std_noise = std(ideal);
cigma_x_noise = std_noise/sqrt(350);