Учебно-методический комплекс по дисциплине «Теоретические основы автоматизированного управления», страница 70

min E{I} =

Если минимум функции I существует, то для него справедливо равенство

(7.26)

Таким образом, операция поиска минимума и вычисления математического ожидания коммутативны. Отсюда можно сделать вывод, что стохастические регуляторы должны описываться уравнениями того же типа, что и детерминированные. Следовательно, управление определяется соотношением

u (N - j) = - KN-j x (N - j).

(7.27)

Устремив N к бесконечности, получим установившееся решение

u (k) = -K×x (k).

(7.28)

Данное выражение описывает стационарный регулятор состояния, принцип действия которого можно пояснить, исходя из следующих соображений. Согласно уравнению (7.21), сигнал управления u(k) должен уменьшать величину x(k+1). Однако как состояние x(k+1), так и управление u(k) зависят только от текущих значений x(k) иv(k) и не зависят от x(k-1),x(k-2), ..., а также от v(k-1), v(k-2), ... . Кроме того, случайное возмущениеv(k) не коррелировано с v(k-1), v(k-2), ... .

Поэтому, при больших значениях N в закон управления должны входить только текущие величины, т. е. u(k) = f(x(k)), и уравнение (7.29) отвечает этому условию. При сравнительно малых значениях N в законе управления необходимо учитывать также случайное начальное состояние x(0) и помеху v(k).

Следовательно, в этом случае необходимо пользоваться нестационарным уравнением оптимального регулятора (7.28). Поскольку оптимизация по отношению к детерминированному начальному состоянию приводит к аналогичному по форме соотношению, оказывается, что уравнение (7.28) описывает регулятор, обеспечивающий оптимальное управление как при детерминированных, так и при случайных помехах (при условии, что в соответствующих критериях качества используются одни и те же весовые матрицы).

Получим теперь выражение для ковариационной матрицы x(k+1) вектора состояния замкнутого контура управления в установившемся режиме. Объединяя соотношения (7.21) и (7.28), запишем уравнение замкнутого контура:

x(k+1) = [A - BK x(k) + F×v(k)

(7.29)

Согласно определению

X(k+1) = [A BK X(k)×[A - B×K]T + F×V×FT.

(7.30)

Это уравнение при k ®¥ дает искомую ковариационную матрицу

X = [A - BK] ×X ×[A - B×K]T + F×V×FT.

(7.31)

Значение критерия качества, которое достигается при включении в контур управления оптимального регулятора состояния, можно определить следующим способом:

E{IN-1,N} = E{xT(N-1)× PN-1,N×x(N-1)+ vT×(N-1)×FT×Q×F×v(N-1)},

E{IN-1} = E{xT(N-1)×PN-1×x(N-1)+ vT×(N-1)×FT×Q×F×v(N-1)}

или для N – 2

E{IN-2} = E{xT(N-2)×PN-2×x(N-2) + vT×(N -2)×FT×Q×F×v(N-2)} +

+ vT×(N-1)×FT×Q×F×v( N-1)}.

Проделаем аналогичным образом все N шагов, учитывая, что v(k) - стационарный случайный процесс:

E{Io} = E{xT (0)×Po×x(0)}+N×E{vT(k)×FT×Q×F×v(k).

(7.32)

Обозначим установившееся значение Po = и возьмем в качестве начального состоянияx(0) случайное возмущение F×v(k). В результате преобразований получим окончательное выражение в виде

(7.33)

7.4. Оптимальные регуляторы с оценкой состояния в случае "белого шума"

В разделе 7.3 предполагалось, что вектор состояния x(k) поддается точному измерению. Однако на практике это далеко не всегда осуществимо, и, как правило, переменные состояния приходится восстанавливать по измеримым координатам системы. Далее будет рассматриваться объект

x(k+1)= A×x(k)+B×u(k)+F×v(k),

(7.34)

у которого вектор измерений выхода определяется уравнением

y(k) = C×x(k) + n(k)

(7.35)

или эквивалентным ему y(k+1)=C×x(k+1)+n(k+1). Считается, что погрешности измерения выходных координат представляют собой векторный "белый шум" со статистиками

E{n(k)} = 0,

                   cov [n(k); t = i -j] = E{n(i)×nT(j)} = N×dij.

(7.36)