انجام محاسبات آماری با متلب

%% Run_Me

%% Clearing
clc;
clear;
close all;

%% Make points
N = 100;
mu = input('Please insert MU  : ');
sigma = input('Please insert sigma : ');
data = normrnd(mu,sigma,[1,N]);

min_data = min(data);
max_data = max(data);
range_data = (max_data - min_data)/10;
start_point = min_data;
end_point = min_data + range_data;
for i = 1:  10
    idx1 = data >= start_point;
    idx2 = data <= end_point;
    idx = and(idx1,idx2);
    mycell{i} = data(idx);
    start_point = min_data + (i * range_data);
    if i == 9
        end_point = max_data;
    else
        end_point = min_data + ((i+1) * range_data);
    end
end
Last_Freq = 0;
for i = 1 : 10
   arr =  mycell{i};
   Table(i,1) = length(arr);
   Freq(i) = length(arr);
   Table(i,2) =  mean(arr);
   Center(i) = mean(arr);
   Last_Freq = Last_Freq + Freq(i);
   Table(i,3) = Last_Freq;
   Freq_Accum(i) = Last_Freq;
   Table(i,4) = Freq(i)/N;
   Freq_Ratio(i) = Freq(i)/N;
end

%% Plot Results in table
Class = {'c1';'c2';'c3';'c4';'c5';'c6';'c7';'c8';'c9';'c10'};
Table_Contact = array2table(Table,'RowNames',Class,'VariableNames',{'Freq';'Center';'Freq_Accum';'Freq_Ratio'});
disp(Table_Contact);


Table2(1) = mean(data);
Table2(2) = var(data);
Table2(3) = std(data);
Table2(4) = mode(data);
Table2(5) = data(ceil(end/2));
Table2(6) = Table2(3) / Table2(1);
Table2(7) = prctile(data,1);
Table2(8) = prctile(data,2);
Table2(9) = prctile(data,3);
Table2(10) = prctile(data,4);
Table2(11) = prctile(data,10);
Table2(12) = prctile(data,20);
Table2(13) = prctile(data,30);
Table2(14) = prctile(data,40);
Table2(15) = quantile(data,0.25);
Table2(16) = quantile(data,0.5);
Table2(17) = quantile(data,0.75);
Table2(18) = quantile(data,1);


%% Plot Results in table
Raw_Name = {'Mean';'Var';'STD';'Mode';'Middle';'CV';'Perc_1';'Perc_2';'Perc_3';'Perc_4' ...
    ;'Dec_1';'Dec_2';'Dec_3';'Dec_4';'Quan_1';'Quan_2';'Quan_3';'Quan_4'};
Table_Contact2 = array2table(Table2','RowNames',Raw_Name,'VariableNames',{'Value'});
disp(Table_Contact2);

figure;
bar(data,0.4);
title('Bar chart of data');

figure;
histogram(data,30);
title('Histogram of data');

fprintf('Finish\n');