В частном случае простейших потоков данную СМО можно описать системой обыкновенных дифференциальных уравнений
|
|
|
|
|
|
P = AP, |
|
|
|
|
|
|
|||
где |
|
|
|
|
|
|
0 |
|
0 |
0 |
|
0 |
|
|
|
− λ µ1 |
|
µ2 |
µ3 |
|
|
|
|
|
|
||||||
|
|
− λ |
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
λ |
0 |
0 |
|
µ2 |
|
0 |
µ3 |
|
0 |
|
|
|
||
− µ1 |
|
|
|
|
|
||||||||||
|
|
− λ − |
|
|
|
|
|
|
|
|
|
|
|
||
|
0 |
0 |
|
0 |
|
µ1 |
|
µ3 |
0 |
|
0 |
|
|
|
|
|
|
− µ2 |
|
|
|
|
|
|
|||||||
|
|
|
|
− λ − |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
|
0 |
|
0 |
|
µ2 |
µ1 |
|
0 |
|
|
|
|
|
|
− µ3 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
− |
λ − µ1 |
− |
|
|
|
|
|
|
|
|
A = |
0 |
λ |
|
λ |
0 |
0 |
0 |
|
µ3 |
|
. |
||||
|
|
− |
µ2 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
− λ − µ2 − |
|
|
|
|
|
|
||
0 |
0 |
|
0 |
0 |
|
0 |
|
0 |
|
µ |
|
|
|
||
|
|
|
|
|
|
|
|
|
− µ3 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
− λ − µ1 |
− |
|
|
|
|
|
0 |
0 |
|
0 |
λ |
|
0 |
|
0 |
µ2 |
|
|
|||
|
|
|
− µ3 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
− µ |
− |
µ |
|
|
|
|
0 |
0 |
|
0 |
0 |
|
λ |
|
λ |
λ |
|
− |
|||
|
|
|
|
1 |
|
|
2 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
− µ3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Решим данную систему численно методом Рунге-Кутта (рис 6). |
|||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P000(t) |
0.9 |
|
|
|
|
|
|
P100(t) |
|
|
|
|
|
|
|
P010(t) |
0.8 |
|
|
|
|
|
|
P001(t) |
|
|
|
|
|
|
P110(t) |
|
|
|
|
|
|
|
|
P011(t) |
0.7 |
|
|
|
|
|
|
P101(t) |
|
|
|
|
|
|
P111(t) |
|
|
|
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
|
0.2 |
|
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
|
0 |
0 |
5 |
10 |
15 |
20 |
25 |
30 |
Рис. 6. Аналитическое решение системы уравнений Колмогорова |
|||||||
для СМО с различимыми каналами сигнатуры T(0,15;0,05;0,03;0,0,0) |
|||||||
|
|
|
с диспетчеризацией |
|
|
||
Т.к. данная СМО описывается однородной марковской цепью, к ней применима теорема о транзитивности марковских процессов, т.е. она имеет предельные вероятности состояний. Рассчитаем их, перейдя от системы обыкновенных дифференциальных уравнений (ОДУ) к системе линейных алгебраических уравнений (СЛАУ):
180
|
|
|
− λ µ1 |
|
µ2 |
µ3 |
0 |
|
0 |
|
|
0 |
|
0 |
|
|
|
||||||
|
|
|
|
|
− λ |
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
λ |
0 |
0 |
|
µ2 |
|
0 |
|
µ3 |
|
0 |
|
|
|
|
|||||
|
|
|
− µ1 |
|
|
|
|
|
|
|
|||||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P000 |
|
||||
|
0 |
0 |
|
− λ − |
0 |
|
µ |
|
µ |
|
|
|
0 |
|
0 |
|
|||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
− µ2 |
|
|
1 |
|
|
3 |
|
|
|
|
|
|
|
P100 |
|
||
|
0 |
|
|
0 |
0 |
|
0 |
− λ − |
0 |
|
µ |
|
|
µ |
|
0 |
|
P |
|
||||
|
0 |
|
|
|
|
|
|
|
010 |
|
|||||||||||||
|
|
|
|
|
|
|
− µ |
3 |
|
|
|
2 |
|
|
1 |
|
|
|
|
P |
|
||
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
001 |
|
||
|
0 |
|
|
0 |
λ |
|
λ |
0 |
|
− λ − µ1 |
− |
0 |
|
|
0 |
|
µ |
|
|
P |
|||
|
|
|
|
|
|
− µ2 |
|
|
|
|
3 |
|
|
110 |
|
||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
− λ − µ2 |
− |
|
|
|
|
|
|
011 |
|
|||
0 |
|
0 |
0 |
|
0 |
0 |
|
0 |
|
|
0 |
|
µ1 |
|
P101 |
||||||||
|
|
|
|
|
− µ3 |
|
|
|
|
|
|
||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P111 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− λ − µ1 |
− |
|
|
|
|
|
|||
|
|
|
|
0 |
0 |
|
0 |
λ |
|
0 |
|
0 |
|
µ2 |
|
|
|
|
|||||
|
|
|
|
|
|
|
− µ3 |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
1 |
1 |
|
1 |
|
1 |
|
|
1 |
|
1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Отсюда PT = (0,053; 0,051; 0,062; 0,082; 0,089; 0,178; 0,093; 0,391).
На рис. 7 – 12 представлены результаты имитационного моделирования работы такой СМО в пакете Matlab.
1 |
|
|
|
|
|
|
0.9 |
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
0.2 |
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
0 |
50 |
100 |
150 |
200 |
250 |
300 |
|
||||||
|
|
|
Tj, Tj+1 - Tj=300tau, tau - длина такта |
|
|
|
Рис. 7. Динамика вероятности отказа СМО с различимыми каналами сигнатуры T(0,15;0,05;0,03;0,0,0) с диспетчеризацией
0.85 |
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
0.75 |
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
0.65 |
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
0.55 |
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
0.45 |
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
0.35 |
50 |
100 |
150 |
200 |
250 |
300 |
|
||||||
|
|
|
Tj, Tj+1 - Tj = 300tau, tau - длина такта |
|
||
Рис. 8. Динамика загруженности прибора 1 СМО с различимыми каналами сигнатуры T(0,15;0,05;0,03;0,0,0) с диспетчеризацией
1 |
|
|
|
|
|
|
0.9 |
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
0.2 |
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
0 |
50 |
100 |
150 |
200 |
250 |
300 |
|
||||||
|
|
|
Tj, Tj+1-Tj=300tau, tau - длина такта |
|
|
|
Рис. 9. Динамика загруженности прибора 2 СМО с различимыми каналами сигнатуры T(0,15;0,05;0,03;0,0,0) с диспетчеризацией
181
1 |
|
|
|
|
|
|
0.9 |
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
0.2 |
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
0 |
50 |
100 |
150 |
200 |
250 |
300 |
|
||||||
|
|
|
Tj, Tj+1-Tj=300tau, tau - длина такта |
|
|
|
Рис. 10. Динамика загруженности прибора 3 СМО с различимыми каналами сигнатуры T(0,15;0,05;0,03;0,0,0) с диспетчеризацией
Полагая λ = const, данную систему можно описать системой уравнений Колмогорова, представляющей вероятности состояний СМО Si как функции времени:
dP0 |
(t) |
= P1µ − P0λ; |
|
|
dt |
|
|||
|
|
|||
dPk |
(t) |
= Pk −1λ + Pk +1µ(k +1) − Pk (λ + kµ), i [2,n −1]; |
(2) |
|
dt |
||||
|
|
|||
dpn (t) = Pn−1λ − Pnnµ, dt
при начальных условиях
P0(0)=1, P1(0) = P2(0) = … = Pk(0) = … = Pn(0) = 0.
Так же как и для любой системы, количество состояний которой конечно, и из каждого состояния можно за конечное число шагов перейти в любое другое состояний, для данной СМО существуют предельные вероятности, вычисляемые с помощью СЛАУ вида:
0 = P1µ − P0λ; |
||
0 = P |
− |
λ+Pk +1µ ( k +1) − Pk (λ + kµ ), k [1,n − 1]; |
k |
|
1 |
0 = Pn−1λ − Pnnµ,
Для простейшего эргодического процесса гибели-размножения справедливы следующие формулы предельных вероятностей состояний Pk:орграфом процесса гибели и размножения, изображенным на рис. 11.
Рис. 11. Граф СМО с общей очередью длины m и k приборами равной производительности в классической нотификации
182
P |
= |
λ0λ1 ... λk −1 P , k = 1,2,...,n |
|
|
|
(4) |
|||||||
k |
|
µ1µ2 ... µk |
0 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||
P |
= |
1+ |
λ0 + |
λ0λ1 |
+ ... + |
λ0λ1 ... λn−1 |
−1 |
(5) |
|||||
µ µ |
|
||||||||||||
0 |
|
|
µ |
|
µ µ |
2 |
... |
µ |
n |
|
|||
|
|
|
1 |
1 2 |
|
1 |
|
|
|
|
|||
Стационарное решение для данной многоканальной СМО с однородными приборами и общей очередью можно записать в виде:
|
|
|
ρ |
k |
|
|
|
|
|
ρ k |
|
|
|
||
P |
= |
|
|
k! |
|
|
= |
|
P |
, k = 0,1,2,..., n |
|||||
n |
|
|
|
|
|
|
|
||||||||
|
k |
|
|
|
ρ |
k |
|
|
|
k! |
0 |
|
|||
|
|
|
∑ |
|
|
|
|
|
|
||||||
|
|
|
k! |
|
|
|
|
|
|
|
(6) |
||||
|
|
|
k =0 |
1 |
|
|
|
|
|||||||
|
|
P0 = |
|
|
|
|
|
, |
|
k = 0,1,2,..., n, |
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
n |
|
ρ k |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
∑ |
|
|
|
k! |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
k =0 |
|
|
|
|
|
|
|
|
||
где ρ = λ/µ.
Классическая СМО с однородными приборами без очереди и отказами является частным случаем СМО с различимыми каналами. Рассмотрим одноканальную СМО без очереди и простейшими входящими и выходящими потоками. Результаты имитационного моделирования работы данной СМО сигнатуры T(0,15;0) при входящем потоке интенсивности λ=0,25 (промежутки времени между сигналами имитационной модели распределены по экспоненциальному закону) представлены на рис. 12.
1 |
|
|
|
|
|
|
|
|
|
|
|
|
0.9 |
|
|
|
|
|
|
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
|
|
|
|
|
|
0.7 |
|
|
|
|
|
|
|
|
|
|
|
|
0.6 |
|
|
|
|
|
|
|
|
|
|
|
|
0.5 |
|
|
|
|
|
|
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
|
|
|
|
|
|
0.3 |
|
|
|
|
|
|
|
|
|
|
|
|
0.2 |
|
|
|
|
|
|
|
|
|
|
|
|
0.1 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
50 |
100 |
150 |
200 |
250 |
300 |
350 |
400 |
450 |
500 |
550 |
600 |
|
||||||||||||
|
|
|
|
Tj, Tj+1-Tj =300tau, tau - длина такта |
|
|
|
|
||||
Рис. 12. Динамика вероятности отказа СМО сигнатуры T(0,15;0)
Отклонение средней вероятности отказа СМО согласно имитационной модели после выхода на стационарный режим от предельной вероятности отказа, вычисленной согласно (5) составило 1,59%.
Таким образом, большинство реальных СМО в той или иной степени отклоняется от классической аксиоматике: трудно обеспечить идеальное совпадение пропускных способностей.
183
Литература
1.Котенко, А.П. Моделирование характеристик систем массового обслуживания с неоднородными приборами / Е.А. Каюрин, В.И. Ларина, А.П. Котенко // Сборник научных трудов «Актуальные направления научных исследований XXI века: теория и практика»: [Статья] материалы конференции. Воронеж: – ВГЛТУ, 2017 г. С.179-182
2.Ларина, В.И. Информационно-навигационное обеспечение современных автоматизированных систем диспетчерского управления транспортом / В.И. Ларина, М.С. Щербаков // Сборник научных трудов «Интеллектуальные информационные системы»: [Статья] материалы конференции в 2 ч. – Воронеж: ВГТУ, 2017 г. С. 318-320.
ФГБОУ ВО «Самарский государственный технический университет», Россия
УДК 681.3
В.В. Решетов, Е.В. Кащенко
РАЗРАБОТКА ИНФОРМАЦИОННОЙ СИСТЕМЫ ВНУТРЕННЕГО АУДИТА НАЛОГОВОЙ СЛУЖБЫ
В данной работе описывается необходимость создания информационной системы для внутреннего аудита Федеральной Налоговой Службы, а также основные её требования.
Федеральная налоговая служба (ФНС России) является федеральным органом исполнительной власти, который выполняет функции по контролю и надзору за соблюдением законодательства о налогах и сборах, за правильностью подсчёта, полнотой и своевременностью внесения в надлежащий бюджет налогов и сборов.
Федеральная налоговая служба и ее территориальные органы составляют единую централизованную систему налоговых органов, которым соответственно нужен внутренний аудит.
Внутренний аудит помогает организации достичь своих целей, применяя систематический, дисциплинированный подход для оценки и повышения эффективности процессов управления рисками, контроля и управления.
Необходимость «общения» между ведомствами для проведения определенных подсчётов, сбора статистики и прочих функций и приводит к необходимости создания данной информационной системы.
Информационная система финансового менеджмента это [1]:
– информационная система, которая отслеживает финансовые события и обобщает информацию;
184