Теория массового обслуживания. Сивохин А.В - 143 стр.

UptoLike

143
0.319168991491612152, 0.219101123495340722, 0.213418946761826100, 0.248310938251221330[ ],
0.319037741491612182, 0.219101123495340722, 0.213468946761826123, 0.248392188251221336[ ]]
P8 := 0.319241786603097210, 0.219101123632489708, 0.213391041580248908, 0.248266048184164534[ ]
Вычисление матрицы переходов после 1, 2, 4 и 8
шагов из заданного начального распределения
P0(1,0,0,0)
> P0:=Vector[row]([1, 0, 0, 0]);
P1:=VectorMatrixMultiply(P0, Pi1);
Pi2:=Pi1^2; # Расчет матрицы вероятностей
переходов Pijk системы за 2 шага
P2:=VectorMatrixMultiply(P1, Pi1);
P2:=VectorMatrixMultiply(P0, Pi2);
Pi4:=Pi1^4; # Расчет матрицы вероятностей
переходов Pijk системы за 4 шага
P4:=VectorMatrixMultiply(P0, Pi4);
Pi8:=Pi1^8; # Расчет матрицы вероятностей
переходов Pijk системы за 8 шагов
P8:=VectorMatrixMultiply(P0, Pi8);
P0 := 1, 0, 0, 0[]
P1 := 0.134000000000000008, 0.234000000000000014, 0.266000000000000014, 0.365999999999999992[ ]
p2
:= [ 0.383599999999999996, 0.218088000000000032, 0.192112000000000005, 0.206200000000000022[ ],
0.357200000000000018, 0.222712000000000021, 0.183888000000000024, 0.236200000000000022[ ],
0.283600000000000019, 0.218088000000000032, 0.242112000000000022, 0.256199999999999984[ ],
0.233600000000000030, 0.218088000000000032, 0.242112000000000022, 0.306200000000000028[ ]]
P2 := 0.383599999999999996, 0.218088000000000032, 0.192112000000000005, 0.206200000000000022[ ]
P2 := 0.383599999999999996, 0.218088000000000032, 0.192112000000000005, 0.206200000000000022[ ]
p4
:= [ 0.327701276800000052, 0.219096438912000036, 0.210234044288000022, 0.242968240000000030[ ],
0.323901603200000054, 0.219117820288000060, 0.211284816512000040, 0.245695760000000042[ ],
0.315201276800000040, 0.219096438912000036, 0.215234044288000026, 0.250468240000000009[ ],
0.307701276800000034, 0.219096438912000064, 0.217734044288000028, 0.255468240000000068[ ]]
P4 := 0.327701276800000052, 0.219096438912000036, 0.210234044288000022, 0.242968240000000030[ ]
p8
:= [ 0.319381491491612158, 0.219101123495340722, 0.213337696761826096, 0.248179688251221332[ ],
0.319313725196562258, 0.219101123952503974, 0.213363429489902112, 0.248221721361032016[ ],
0.319168991491612152, 0.219101123495340722, 0.213418946761826100, 0.248310938251221330[ ],
        [0.319168991491612152, 0.219101123495340722, 0.213418946761826100, 0.248310938251221330],
        [0.319037741491612182, 0.219101123495340722, 0.213468946761826123, 0.248392188251221336]]

    P8 := [0.319241786603097210, 0.219101123632489708, 0.213391041580248908, 0.248266048184164534]


Вычисление матрицы переходов после 1, 2, 4 и 8
шагов из заданного начального распределения
P0(1,0,0,0)
> P0:=Vector[row]([1, 0, 0, 0]);
P1:=VectorMatrixMultiply(P0, Pi1);
Pi2:=Pi1^2;    # Расчет матрицы вероятностей
переходов Pijk системы за 2 шага
P2:=VectorMatrixMultiply(P1, Pi1);
P2:=VectorMatrixMultiply(P0, Pi2);
Pi4:=Pi1^4;    # Расчет матрицы вероятностей
переходов Pijk системы за 4 шага
P4:=VectorMatrixMultiply(P0, Pi4);
Pi8:=Pi1^8;    # Расчет матрицы вероятностей
переходов Pijk системы за 8 шагов
P8:=VectorMatrixMultiply(P0, Pi8);
                                           P0 := [1, 0, 0, 0]

    P1 := [0.134000000000000008, 0.234000000000000014, 0.266000000000000014, 0.365999999999999992]

p2 := [[0.383599999999999996, 0.218088000000000032, 0.192112000000000005, 0.206200000000000022 ],
        [0.357200000000000018, 0.222712000000000021, 0.183888000000000024, 0.236200000000000022],
        [0.283600000000000019, 0.218088000000000032, 0.242112000000000022, 0.256199999999999984],
        [0.233600000000000030, 0.218088000000000032, 0.242112000000000022, 0.306200000000000028]]

    P2 := [0.383599999999999996, 0.218088000000000032, 0.192112000000000005, 0.206200000000000022]

    P2 := [0.383599999999999996, 0.218088000000000032, 0.192112000000000005, 0.206200000000000022]

p4 := [[0.327701276800000052, 0.219096438912000036, 0.210234044288000022, 0.242968240000000030 ],
        [0.323901603200000054, 0.219117820288000060, 0.211284816512000040, 0.245695760000000042],
        [0.315201276800000040, 0.219096438912000036, 0.215234044288000026, 0.250468240000000009],
        [0.307701276800000034, 0.219096438912000064, 0.217734044288000028, 0.255468240000000068]]

    P4 := [0.327701276800000052, 0.219096438912000036, 0.210234044288000022, 0.242968240000000030]

p8 := [[0.319381491491612158, 0.219101123495340722, 0.213337696761826096, 0.248179688251221332 ],
        [0.319313725196562258, 0.219101123952503974, 0.213363429489902112, 0.248221721361032016],
        [0.319168991491612152, 0.219101123495340722, 0.213418946761826100, 0.248310938251221330],

                                                  143