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

UptoLike

67
10. Оценка, интерполяция и визуализация
нормированной корреляционной
функции Rx(Tk, Tl) значений в сечениях
случайной функции X(t)
> Rx:=Matrix(1..m, 1..m);
# Определение матрицы для нормированной
корреляционной функции Rx(Tk, Tl)
whattype(Kx);
for k from 1 to m do
for l from 1 to m do
Rx[k, l]:=Kx[k, l]/(sqrt(Kx[k,
k])*sqrt(Kx[l, l]));# Расчет нормированной
корреляционной функции Rx(Tk, Tl) с помощью
циклов
end do;
end do;
Rx:=[seq([seq(Kx[k, l]/(sqrt(Kx[k,
k])*sqrt(Kx[l, l])), k=1..m)], l=1..m)];#
Аналитический расчет нормированной
корреляционной функции
whattype(Rx);
matrixplot(Rx, heights=histogram, axes=normal,
gap=0.25, style=patchcontour);# Визуализация
корреляционных моментов Kx(Tk, Tl)
Rx :=
000000
000000
000000
000000
000000
000000
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
?
10. Оценка, интерполяция и визуализация
     нормированной корреляционной
 функции Rx(Tk, Tl) значений в сечениях
         случайной функции X(t)
> Rx:=Matrix(1..m, 1..m);
# Определение матрицы для нормированной
корреляционной функции Rx(Tk, Tl)
whattype(Kx);

for k from 1 to m do
    for l from 1 to m do
        Rx[k, l]:=Kx[k, l]/(sqrt(Kx[k,
k])*sqrt(Kx[l, l]));# Расчет нормированной
корреляционной функции Rx(Tk, Tl) с помощью
циклов
    end do;
end do;

Rx:=[seq([seq(Kx[k, l]/(sqrt(Kx[k,
k])*sqrt(Kx[l, l])), k=1..m)], l=1..m)];#
Аналитический расчет нормированной
корреляционной функции
whattype(Rx);

matrixplot(Rx, heights=histogram, axes=normal,
gap=0.25, style=patchcontour);# Визуализация
корреляционных моментов Kx(Tk, Tl)
                   ?0   0   0    0   0   0?
                   ?                      ?
                   ?0   0   0    0   0   0?
                   ?                      ?
                   ?                      ?
                   ?0   0   0    0   0   0?
             Rx := ?                      ?
                   ?0   0   0    0   0
                                          ?
                                         0?
                   ?
                   ?                      ?
                   ?0   0   0    0   0   0?
                   ?                      ?
                   ?                      ?
                   ?0   0   0    0   0   0?



                            67