Модель функционирования З П С, оснащённой «n» однородными факторами воздействия на поставки товаров и услуг из другого региона
Вывод основных ДУ
Полагаем, что время обслуживания З П С величина случайная и подчинена показательному закону распределения с параметром ν. Поэтому вероятность того, что время обслуживания поставки не превосходит t определяется выражением
Р(t) = 1 - е 
Вероятность противоположного события равна g(t) = е
может находиться в следующих состояниях:
А
- все факторы ЗПС не проявляют себя.
А
- k – факторов себя проявляют себя, а остальные свободны от облуживания.
А
- все факторы ЗПС проявили себя и обслуживают поставки.
Составим дифференциальные уравнения этих состояний ЗПС.
Обозначим через Δt очень маленький промежуток времени. Составим ДУ для состояния А
. Оно возможно в следующих несовместных случаях:
- в момент времени t все факторы З П С себя не проявляют. За время Δt в области действия З П С не проявилась не одна поставка. Вероятность этого события равна
Р
(t)е
(1)
- в момент времени t один из факторов себя уже проявляет. За время Δt в области действия З П С не проявилась ни одна поставка, а действующий фактор прекратил своё действие. Вероятность этого события равна
Р
(t) (1 - е
) е
(2)
Где Р
(t) –вероятность того, что один фактор себя проявил воздействием на поставку.
Так как рассмотренные события несовместны, то составим уравнение состояния А 
Р
(t+ Δt) = Р
(t) е
+ Р
(t) (1 - е
)е
е
(3)
Где Р
(t+ Δt) - вероятность того, что во время (t + Δt) ни один из факторов
не будет воздействовать на поставки;
Р
(t) - вероятность нахождения З П С в состоянии А
;
е
- вероятность непоявления за Δt в области действия З П С ни
одной поставки;
(1 - е
) - вероятность того, что за время Δt в области действия З П С
один из факторов прекратит свое воздействие на поставки.
Величину е
можно разложить в ряд Тейлора и представить в следующем виде
е
≈ 1 - λ Δt + …,
а
1 - е
≈ νΔt + …,
Учитывая малость величины Δt. Можно уравнение (3) представить в виде, пренебрегая величинами более высокого порядка малости,
′ Р
(t+ Δt) = Р
(t) (1 - λ Δt) + Р
(t) νΔt (1 - λ Δt). (4)
Преобразуем и разделим обе части уравнения (4) на Δt получим следующее соотношение
Р
(t+ Δt) - Р
(t) = - Р
(t) λ Δt + Р
(t) νΔt - Р
(t)∙νΔt∙ λΔt). (5)
Так как в соотношении Р
(t)∙νΔt∙ λΔt Δt∙Δt = Δt
, то можно допустить, что Δt
≈0, то и само соотношение Р
(t)∙νΔt∙ λΔt ≈ 0. Тогда справедливо соотношение
Р
(t+ Δt) - Р
(t) = - Р
(t) λ Δt + Р
(t) νΔt. (6)
Разделим обе части уравнения (6) на Δt получим следующее соотношение
= - Р
(t) λ + Р
(t) ν.
и, перейдя к пределу Δt → 0, получим
(
) =
Р
(t) = - Р
(t) λ + Р
(t) ν. (7)
Составим ДУ для состояния А
. Оно возможно в следующих трёх несовместных случаях:
- в момент времени t «k» факторов З П С занято обслуживанием поставок. За время Δt в область действия З П С не проявилась ни одна поставка и ни один из занятых факторов не освободился. Вероятность этого события равна
Р
(t) (1 - λ Δt)(1 – kνΔt); (8)
- в момент времени t З П С находилась в состоянии А
. За время Δt в области действия З П С проявилась одна поставка, но ни один из действующий факторов не прекратил обслуживание своей поставки. Вероятность этого события равна
Р
(t) λ Δt (1 – kνΔt); (9)
- в момент времени t З П С находилась в состоянии А
. За время Δt в освободился один факторов З П С и в области действия З П С не проявилась ни одна поставка. Вероятность этого события равна
Р
(t) (1 - λ Δt)(1 + k) νΔt. (10)
Тогда
Р
(t+Δt)=Р
(t)(1-λΔt)(1–kνΔt)+ Р
(t)(1– kνΔt)λ Δt+Р
(t)(1-λΔt)(1+k)νΔt. (11)
После аналогичных преобразований получаем
Р
(t) = - (λ+kν)Р
(t) + λ Р
(t)+ Р
(t) (k+1)ν. (12)
Это уравнение справедливо для случая 0≤ k≤ n.
Рассмотрим состояние А
.
- в момент времени t З П С находилась в состоянии А
. За время Δt ни один из занятых факторов не освободился. Вероятность этого события равна
Р
(t) (1 – nνΔt); (13)
После соответствующих преобразований и перехода к пределу при Δt → 0, получим
Р
(t) = - nνР
(t) + λ Р
(t). (14)
В итоге мы получили систему дифференциальных уравнений, которая описывает вероятности переходов состояний З П С, оснащённой n однородными факторами воздействия на поставки товаров и услуг из другого региона.
Р
(t) = - Р
(t) λ + Р
(t) ν.
…………………………………….
Р
(t) = - (λ+kν)Р
(t) + λ Р
(t)+ Р
(t) (k+1)ν. (15)
………………………………………
Р
(t) = - nνР
(t) + λ Р
(t).
Совокупность этих уравнений получила название системы уравнений Эрланга.