Главная страница
Навигация по странице:

  • Исходные данные по варианту

  • Основные определения и алгоритмы моделирования случайных величин Нормальное распределени

  • Равномерное распределение

  • Экспоненциальное распределение

  • Основные определения

  • Графики и характеристики в Excel

  • 5485756_Доп_2. Практическая работа 2 Моделирование случайных величин с заданным законом распределения Вариант 1 Проверил 2020


    Скачать 332 Kb.
    НазваниеПрактическая работа 2 Моделирование случайных величин с заданным законом распределения Вариант 1 Проверил 2020
    Дата25.11.2020
    Размер332 Kb.
    Формат файлаdocx
    Имя файла5485756_Доп_2.docx
    ТипПрактическая работа
    #153942

    Министерство науки и высшего образования РФ

    Практическая работа №2

    «Моделирование случайных величин с заданным законом распределения»

    Вариант №1


    Выполнил:

    Проверил:

    2020

    Задание

    1. Необходимо выполнить следующие задания:

    для каждого из следующих трех видов распределений (параметры распределений взять из п.2 данной работы):

    • нормального N(m, σ) (σ > 0),

    • равномерного R(a, b) (a < b),

    • экспоненциальногоE(λ) (λ > 0)

    написать программу вывода в графическом окне графиков функций распределения с использованием функций Excel.

    2. Написать программу, осуществляющую моделирование массивов объемом 300 для законов распределения: нормального, равномерного, экспоненциального.

    Для каждого из этих массивов вывести в графическом окне график гистограммы относительных частот группированной выборки (шесть – восемь интервалов группировки). Продемонстрировать полученные результаты преподавателю.
    Исходные данные по варианту

    Параметры моделируемых массивов:

    для нормального распределения N(m, σ) (σ> 0): m=1, σ= ;

    для равномерного распределения R(a ;b) (a
    для экспоненциального распределения E(λ)( λ> 0 ): .

    Основные определения и алгоритмы моделирования случайных величин
    Нормальное распределение, также называемое распределением Гаусса – распределение вероятностей, которое в одномерном случае задаётся функцией плотности вероятности, совпадающей с функцией Гаусса:



    Математическое ожидание .

    Дисперсия .

    Коэффициент асимметрии: 0;

    Коэффициент эксцесса: 0.

    Аналитическое выражение для моделирования нормальной СВ имеет вид:



    где α1 и α2 - независимые равномерно распределенные случайные числа на интервале [0; 1), m – математическое ожидание СВ, σ – среднеквадратическое отклонение.

    Равномерное распределение – распределение случайной вещественной величины, принимающей значения, принадлежащие интервалу [a, b], характеризующееся тем, что плотность вероятности на этом интервале постоянна.



    Функция распределения:



    Математическое ожидание .

    Дисперсия .

    Коэффициент асимметрии: 0;

    Коэффициент эксцесса: -1,2.

    Аналитическое выражение для моделирования имеет вид:



    где α - независимые равномерно распределенные случайные величины на интервале [0; 1).

    Экспоненциальное распределение - абсолютно непрерывное распределение, моделирующее время между двумя последовательными свершениями одного и того же события.

    Плотность распределения



    Функция распределения:



    Математическое ожидание .

    Дисперсия .

    Коэффициент асимметрии: 2;

    Коэффициент эксцесса: 6.

    Аналитическое выражение для моделирования имеет вид:



    где α - независимые равномерно распределенные случайные числа на интервале [0; 1),

    Основные определения

    Случайная величина – величина которая при испытаниях принимает одно из возможных значений, наперед неизвестно какое. Бывают дискретными и непрерывными.

    Законом распределения дискретной случайной величины называют соответствие между возможными значениями случайной величины и вероятностями их появления. Сумма всех вероятностей Σpi = 1. Закон распределения также может быть задан аналитически (формулой) и графически (многоугольником распределения, соединяющим точки (xi; pi).

    Функция распределения – функция, характеризующая распределение случайной величины или случайного вектора; вероятность того, что случайная величина X примет значение, меньшее или равное х, где х – произвольное действительное число.

    Плотностью распределения вероятностей непрерывной случайной величины Х называется функция f(x) – первая производная от функции распределения F(x). Смысл плотности распределения состоит в том, что она показывает, как часто появляется случайная величина Х в некоторой окрестности точки х при повторении опытов.

    Коэффициент асимметрии – числовая характеризующая степени несимметричности распределения данной случайной величины.

    Коэффициент эксцесса (коэффициент островершинности) – мера остроты пика распределения случайной величины.


    Графики и характеристики в Excel

    Функция СЛЧИС, генерирующая случайные последовательности, равномерно распределенные на отрезке [0,1]. Ее синтаксис СЛЧИС() и она не имеет аргументов – располагаем выборку в столбце А. Далее согласно алгоритму генерируем выборку из заданного распределения объёмом 300 единиц.

    Нормальное распределение

    Выборка из нормального распределения расположена в столбце С:



    Для построения плотности и функции распределения использовали функцию НОРМ.РАСП. Например,

    J2 = =НОРМ.РАСП($H2;F$1;F$2;ЛОЖЬ);

    I2 =НОРМ.РАСП(H2;$F$1;$F$2;ИСТИНА)



    Рис. Функция и плотность распределения

    Для построения гистограммы вычисляем частоты попадания в каждый интервал.



    Здесь Norm – именованный диапазон с нормально распределенной СВ.



    Рис. Гистограмма относительных частот

    Случайные числа для выборок объемом 100, 1000, 10000, 100000 расположены в столбцах A-G.

    Вычисляем их характеристики:



    Графики зависимости оценки мат. ожидания, дисперсии, коэффициента асимметрии и эксцесса от длины выборки:




    Размер выборки

    100

    1000

    10000

    100000

    Среднее

    0,9412

    1,0384

    1,0119

    0,9989

    Дисперсия

    1,0601

    0,9098

    0,9700

    0,9922

    Эксцесс

    -0,0017

    -0,0863

    -0,0956

    0,2451

    Асиметрия

    -0,2838

    -0,0002

    -0,0463

    -0,0226


    Равномерное распределение

    Фрагмент выборки (столбец B):


    Для построения плотности и функции распределения используем аналитическое выражение. Например,

    E6=ЕСЛИ(И(D6>$E$1;D6<$G$1);1/($G$1-$E$1);0)

    F6=ЕСЛИ(И(D6>=$E$1;D6<$G$1);(D6-$E$1)/($G$1-$E$1);ЕСЛИ(D6<$E$1;0;1))



    Рис. Функция и плотность распределения
    Для построения гистограммы вычисляем частоты попадания в каждый интервал.



    Здесь rav – именованный диапазон с равномерно распределенной СВ.


    Рис. Гистограмма относительных частот

    Случайные числа для выборок объемом 100, 1000, 10000, 100000 расположены в столбцах A-G.

    Характеристики:



    Размер выборки

    100

    1000

    10000

    100000

    Среднее

    1,4748

    1,4983

    1,4985

    1,5001

    Дисперсия

    0,0817

    0,0787

    0,0830

    0,0834

    Эксцесс

    -1,2407

    -1,1753

    -1,2001

    -1,2031

    Асиметрия

    0,1916

    0,0235

    0,0118

    -0,0023


    Графики зависимости оценки мат. ожидания, дисперсии, коэффициента асимметрии и эксцесса от длины выборки:


    Экспоненциальное распределение

    Фрагмент выборки (столбец B):



    Для построения плотности и функции распределения использовали функцию ЭКСП.РАСП. Например,

    E6= ЭКСП.РАСП(D6;$E$1;ЛОЖЬ);

    F6= ЭКСП.РАСП(D6;$E$1;ИСТИНА)



    Рис. Функция и плотность распределения

    Для построения гистограммы вычисляем частоты попадания в каждый интервал.

    Интервалы

     

    Левая граница

    Правая граница

    Кол-во, fi

    Середина интервала, xi

    Относительные частоты

    Накопленная частота

    0,006-1,017

    1

    0,006

    1,017

    113

    0,512

    0,3767

    0,3767

    1,017-2,029

    2

    1,017

    2,029

    76

    1,523

    0,2533

    0,6300

    2,029-3,041

    3

    2,029

    3,041

    42

    2,535

    0,1400

    0,7700

    3,041-4,053

    4

    3,041

    4,053

    29

    3,547

    0,0967

    0,8667

    4,053-5,065

    5

    4,053

    5,065

    18

    4,559

    0,0600

    0,9267

    5,065-6,077

    6

    5,065

    6,077

    11

    5,571

    0,0367

    0,9633

    6,077-7,089

    7

    6,077

    7,089

    7

    6,583

    0,0233

    0,9867

    7,089-8,101

    8

    7,089

    8,101

    3

    7,595

    0,0100

    0,9967

    8,101-9,112

    9

    8,101

    9,112

    1

    8,606

    0,0033

    1,0000

     

     

     

     

    300

     

    1,0000

     




    Рис. Гистограмма относительных частот
    Случайные числа в количестве 100000 штук получим с помощью алгоритма Y где x - независимые равномерно распределенные случайные числа на интервале [0; 1). x располагаются в столбце A, а в столбце B-H.

    Характеристики:



    Размер выборки

    100

    1000

    10000

    100000

    Среднее

    2,1353

    1,9577

    1,9947

    1,9863

    Дисперсия

    3,6573

    3,6114

    3,9333

    3,9642

    Эксцесс

    7,0441

    4,4532

    5,6984

    6,0133

    Асиметрия

    2,2426

    1,8382

    1,9926

    2,0250


    Графики зависимости оценки мат. ожидания, дисперсии, коэффициента асимметрии и эксцесса от длины выборки:



    Таблица сравнения теоретических числовых характеристик СВ и характеристик, полученных по моделируемой выборке

     

    M[X]

    D[X]

    Ex[X]

    Sk[X]

     




    Теор.

    Оценка

    Теор.

    Оценка

    Теор.

    Оценка

    Теор.

    Оценка







    Нормальное N(1,1)

    N =100

    1

    0,94124

    1

    1,0601

    0

    -0,00166

    0

    -0,28378







    N =1000

     

    1,03842

     

    0,90985

     

    -0,08634

     

    -0,0002







    N =10000

     

    1,01187

     

    0,97001

     

    -0,09563

     

    -0,04635







    N =100000

     

    0,99886

     

    0,99224

     

    0,245138

     

    -0,02257







    Равномерное R(1, 2)

    N =100

     1,5

    1,47484

     0,08333

    0,08169

    -1,2 

    -1,24071

    0

    0,19158







    N =1000

     

    1,49832

     

    0,07871

     

    -1,17528

     

    0,02345







    N =10000

     

    1,49855

     

    0,08299

     

    -1,20006

     

    0,01185







    N =100000

     

    1,50014

     

    0,08341

     

    -1,20314

     

    -0,00228







    Экспоненциальное E(0,5)

    N =100

    2

    2,13526

     4

    3,65726

     6

    7,044085

    2

    2,24257







    N =1000

     

    1,95769

     

    3,61135

     

    4,45321

     

    1,83816







    N =10000

     

    1,99468

     

    3,93325

     

    5,698422

     

    1,99261







    N =100000

     

    1,98632

     

    3,9642

     

    6,013272

     

    2,025








    Таким образом, с ростом N выборочные значения приближаются к теоретическим.

    Полученные гистограммы распределений приближаются к теоретическим плотностям распределений, т.е. разработанные датчики случайных чисел реализованы верно. Изученные возможности пакета Excel позволяют моделировать и анализировать одномерные случайные величины. Для основных законов распределения имеются встроенные функции, остальные законы распределения можно получить с помощью самостоятельно написанных алгоритмов.


    написать администратору сайта