Розподіл і потоки Ерланга
Випадкові величини з експоненціальним розподілом не завжди адекватно описують деякі реальні процеси та події, наприклад час обслуговування і моменти надходження вимог до СМО. Для більш точного моделювання таких процесів доцільніше використовувати гамма-розподілені випадкові величини або ті, що мають розподіл Ерланга. Розподіл Ерланга є результатом підсумовування взаємно незалежних і однаково розподілених експоненціальних випадкових величин і є окремим випадком гамма-розподілу.
Функція щільності розподілу Ерланга k-το порядку з інтенсивністю λ має такий вигляд
Математичне сподівання і дисперсія розподілу Ерланга визначаються як М[х] = 1/kλ, D [x]=1/kλ2.
Для моделювання розподілу Ерланга використовують метод згорток випадкових величин з експоненціальними функціями розподілу. Для цього треба лише обчислити суму k експоненціально розподілених випадкових величин. Зі збільшенням k розподіл Ерланга наближається до нормального.
4.8 Статистична обробка результатів моделювання
Основою для обчислення статистичної оцінки параметра системи є реалізація випадкової величини, яка формується під час прогонів імовірнісної імітаційної моделі. Статистична оцінка також є функцією від випадкових величин, які дістають у результаті прогонів моделі, тому і згадана оцінка є випадковою величиною, закон розподілу якої залежить від закону розподілу досліджуваної випадкової величини та оцінюваного параметра. Чим більше реалізацій випадкової величини, тим точнішу статистичну оцінку параметра системи ми отримуємо. У таких умовах обробка результатів моделювання повинна провадитись лише з використанням методів і алгоритмів, які є оптимальними з погляду затрат часу та використання ресурсів комп'ютера. Під час вибору таких засобів необхідно враховувати, що всі статистичні оцінки мають бути ще й незміщеними, ефективними і спроможними.
Розглянемо деякі методи обчислення основних статистичних оцінок.
4.8.1 Оцінювання ймовірності
Оцінкою ймовірності ρ настання деякої події А є її частість:
де т — кількість випробувань, під час за яких випадкова подія спостерігалась; N — загальна кількість випробувань. Для її використання зазвичай на програмному рівні організовують два лічильники, один з яких призначено для підрахунку загальної кількості випробувань Ν, а другий — кількості успішних випробувань т.
4.8.2 Оцінювання розподілу випадкової величини
Для оцінювання функції розподілу випадкової величини, як звичайно, будується гістограма. Під час її побудови область можливих значень випадкової величини розбивають на η діапазонів і підраховують кількість попадання значень випадкової величини в конкретний інтервал - тk (k = 1,..., n). Оцінка ймовірності попадання випадкової величини в k-й. інтервал має такий вигляд:
Цю величину називають відносною частістю. У процесі моделювання під час підрахунку значень тk кожному інтервалу ставлять у відповідність окремий лічильник.
4.8.3 Оцінювання математичного сподівання
Для оцінювання математичного сподівання випадкової величини використовується формула:
де хk — значення випадкової величини, що належить k-му інтервалу; тk — кількість попадань значень випадкової величини в інтервал; N — загальна кількість випробувань. У більш простому випадку для оцінювання математичного сподівання випадкової величини можна використати звичайне середнє арифметичне:
Щоб запобігти непотрібному завантаженню пам'яті, суму доцільніше підраховувати шляхом поступового накопичення.