Matlab-statistics da eksperiment natijalariga ishlov berish va baholash

Amaliy mashg’ulot № 6
MATLAB-STATISTICS DA EKSPERIMENT NATIJALARIGA ISHLOV BERISH VA BAHOLASH

Телекоммуникация соҳасининг теҳник воситаларини синаб кўришда уларнинг қувватини билдирувчи (carrier power) катталиклар ўлчаниб кўрилади, бундай тажрибаларда йигирмадан кам бўлмаган ўлчаш натижалари лабаратория журналларига қайд қилинади

% The following section perform measurement of carrier power
% The example test is a conducted measurement. The measurement was repeated
% 45 times. The following results were obtained (before correcting for
% cabling and attenuator network insertion loss):
% 0,57 mW; 1,76 mW; 1,88 mW; 2,54 mW; 2,98 mW; 3,02 mW; 3,13 mW; 3,52 mW;
% 3,72 mW; 3,88 mW; 4,05 mW; 4,10 mW; 4,29 mW; 4,44 mW; 4,56 mW; 4,72 mW;
% 4,85 mW; 4,93 mW; 5,05 mW; 5,16 mW; 5,23 mW; 5,37 mW; 5,45 mW; 5,50 mW;
% 5,54 mW; 5,64 mW; 5,76 mW; 5,86 mW; 5,95 mW; 6,06 mW; 6,15 mW; 6,36 mW;
% 6,54 mW; 6,72 mW; 6,88 mW; 7,13 mW; 7,21 mW; 7,38 mW; 7,62 mW; 7,91 mW;
% 8,23 mW; 8,46 mW; 8,62 mW; 9,35 mW; 9,98 mW
>> x = [0.57, 1.76, 1.88, 2.54, 2.98, 3.02, 3.13, 3.52, 3.72, 3.88, 4.05, 4.10, 4.29, 4.44, 4.56, 4.72, 4.85, 4.93, 5.05, 5.16, 5.23, 5.37, 5.45, 5.50, 5.54, 5.64, 5.76, 5.86, 5.95, 6.06, 6.15, 6.36, 6.54, 6.72, 6.88, 7.13, 7.21, 7.38, 7.62, 7.91, 8.23, 8.46, 8.62, 9.35, 9.98];


Тақсимот қонунини визуал усулда аниқлаш
Авваламбор гистограммани чизиб оламиз
% Vertical Bar Graph Plot
[m, xout] = hist(x, 10);
plot_bargraph = bar(xout, m/length(x));

Сўнгра гистограмма устида қуйидаги тақсимот қонунларини чизамиз.
Нормал тақсимот қонуни (Гаусс қонуни)
Экспоненциал тақсимот қонуни
Тўғрибурчакли тақсиомт қонуни
Релея тақсимот қонуни

Шу билан бирга ҳар тақсимот қонунинг параметрлари (MU, SIGMA, LAMBDA, A, B) лар эканлиги ҳисобга оламиз.

holdа on
% Perform Normal Distribution
[MU, SIGMA] = normfit(x);
y = normpdf(x, MU, SIGMA);
plot_normdst = plot(x, y, '--ko', 'MarkerEdgeColor', 'k', 'MarkerFaceColor', 'g', 'MarkerSize', 6);
% Perform Exponential Distribution
MU = expfit(x);
y = exppdf(x, MU);
plot_expdst = plot(x, y, '-.r*');
% Perform Uniform (Continuous) Distribution
[A, B] = unifit(x);
y = unifpdf(x, A, B);
plot_unifdst = plot(x, y, 'b', 'LineWidth', 2);
% Perform Rayleigh Distribution
B = raylfit(x);
y = raylpdf(x, B);
plot_rayldst = plot(x, y, '.m');
grid on
hold off

Қабул қилинган натижаларни таҳлил қиламиз.