Гамма-распределение экспериментальной величины

W = [3  4   3   3   3   3   3   3   5   6   6   5   4   2   0   0   1   1   1   1   1   2   1   3   2   2   1   3   3   4   3   2   2   2   1   0   0   0   0   0   2   1   1   2   3   2   2   2   3   4]

alpha = mean(W); %мат.ожидание
lamda = var(W); %дисперсия

clear all, clc, close all
 
W = [3 4 3 3 3 3 3 3 5 6 6 5 4 2 0 0 1 1 1 1 1 2 1 3 2 2 1 3 3 4 3 2 2 2 1 0 0 0 0 0 2 1 1 2 3 2 2 2 3 4];
alpha = mean(W); %мат.ожидание
lamda = var(W); %дисперсия
a = alpha; b = lamda;
 
R = gampdf(W, a, b);
[var, idx] = max(R); %находим максимальную вероятность события и его индекс
Wsort = sort(W); %сортируем массив данных
Wist = W(idx) %находим значение скорости ветра, соответствующее максимальной вероятности
plot(W,R), grid on

x = -9:0.2:9

T = [-37,6 -37,1 -36,2 -33,5 -32,8 -32,9 -34,1 -32,3 -38,7 -30,4 -30,8 -34,3 -38,8 -44,7 -47,7 -49,0 -48,2 ...
    -47,7 -46,2 -45,0 -46,3 -46,5 -43,4 -39,7 -34,6 -27,7 -28,6 -30,1 -36,0 -38,9 -41,9 -44,7 -46,9 -47,4 ...
    -46,9 -45,3 -48,2 -48,5 -47,1 -46,3 -45,1 -45,1 -43,1 -41,3 -41,2 -41,4 -41,1 -40,9 -40,7 -39,7];

clear all, clc, close all
 
T = [-37,6 -37,1 -36,2 -33,5 -32,8 -32,9 -34,1 -32,3 -38,7 -30,4 -30,8 -34,3 -38,8 -44,7 -47,7 -49,0 -48,2 ...
    -47,7 -46,2 -45,0 -46,3 -46,5 -43,4 -39,7 -34,6 -27,7 -28,6 -30,1 -36,0 -38,9 -41,9 -44,7 -46,9 -47,4 ...
    -46,9 -45,3 -48,2 -48,5 -47,1 -46,3 -45,1 -45,1 -43,1 -41,3 -41,2 -41,4 -41,1 -40,9 -40,7 -39,7];
 
% T = -9:0.2:9;
Tsr = mean(T);
sigma = 2;
 
R = normpdf(T, Tsr, sigma);
[var, idx] = max(R);
Tsort = sort(T);
Tist = T(idx)
 
plot(T, R), grid on

M = [1.606, 1.538, 1.905, 2.400, 3.660, 4.536, 2.150, 2.682, 2.773, 2.838, 1.584, 2.181]; %мат.ожидание
Disp = [0.538, 0.595, 0.644, 0.753, 1.057, 1.392, 0.912, 0.941, 0.876, 0.868, 0.594, 0.703]; %дисперсия

gamrnd(alpha, lamda); 
normrnd(mu, sigma)