Занятие №5. Изучение моделирования канального уровня и физического уровня модели osi на основе matlab-stateflow

5.3 Исследование модели Петровича для дискретного канала
Необходимо исследовать параметры модели Петровича для дискретного канала:

• 
  длина кодовой комбинации - n = 32 бита;
  
• 
  переходов вероятности:
  
  1. 
     р ₁₂ = 0,21;
     
  2. 
     р24 ₌ 0,76;
     
  3. 
     р ₃₁ = 0,25;
     
  4. 
     р ₄₃ = 0,23;
     
• 
  любовник вероятности:
  
  1. 
     р _{2 о} = 0,49;
     
  2. 
     р _{3 о} = 049.
     

После завершения моделирования в файле error_num . мат будет содержаться выборка значений краткости ошибки. Для определения ее параметров используется программа, приведенная в списке 5.1.

Листинг 5.1 – Программа для определения параметров

Кмакс = 32; % наибольшая возможная кратность ошибки
загрузить('error_num.mat'); % загрузки переменныx из файловой и рабочей области памяти
х = номер_ошибки ( 2, :);
х = сортировать (х);
п = длина (х);
f = п - 1;
хмин = х(1);
хмакс = х(п);
fprintf('Объем выборки = %d\n', n);
fprintf('Число степней свободы = %d\n', f);
fprintf('Минимальное значение = %g\n', xmin);
fprintf('Максимальное значение = %g\n\n', xmax);

k = 0:Kmax; % центральные значения интервалов гистограммы

q = -1:Kmax; % точки переходов эмпирической функции распределения
% определение эмпирической достоверности
N = ист (х, к); % определение количества попадани в интервале
г = Н./н;
% определение значения эмпирической функции распределения
у( 1) = 0;
для I = 2:Kmax+1,
y( I) = y(I-1) + z(I-1);
конец
участок( 2, 1, 1 );
история ( х, к);
title('Гистограмма распределения вероятности ошибки');
xlim( [-1, Кмакс+1]);
подсюжет( 2, 1, 2 );
лестница ( q, y);
title('Эмпирическая функция распределения кратности ошибки');
год([0, 1.1]);
xlim([-1, Кмакс+1]);
disp('Эмпирическая вероятность к-кратной ошибки :' );
для I = 0:Kmax,
fprintf( 'P(%g) = %g\n', I, z(I+1));
конец
Мх = среднее (х);
Sx = станд.( х);
Дх = Sх^2;
Ах = сумма ( (х - Mx).^3)/(n - 1)/Sx^3;
Ex = сумма( (x - Mx).^4)/(n - 1)/Dx^2-3;
дисп('');
disp('Выборочные параметры распределения :' );
fprintf('Математическое ожидание = %g\n', Mx);
fprintf('Среднеквадратичное отклонение = %g\n', Sx);
fprintf('Дисперсия = %g\n', Dx);
fprintf('Асимметрия = %g\n', Ax);
fprintf('Избыток = %g\n\n\n\n', Ex);

В результате выполнения командной программы получаются следующие значения:

Объем выборки = 10000001
Число степней свободы = 10000000
Минимальное значение = 0
Максимальное значение = 17
Математическая формула = 5,16887
Среднеквадратичное отклонение = 2,25987
Дисперсия = 5,10701
Асимметрия = 0,405941
Эксцентриситет = 0,118124