ÐÐ¾Ð»Ð½Ð°Ñ Ð²ÐµÑÑÐ¸Ñ - СамаÑÑкий гоÑÑдаÑÑÑвеннÑй аÑÑокоÑмиÑеÑкий ...
ÐÐ¾Ð»Ð½Ð°Ñ Ð²ÐµÑÑÐ¸Ñ - СамаÑÑкий гоÑÑдаÑÑÑвеннÑй аÑÑокоÑмиÑеÑкий ...
ÐÐ¾Ð»Ð½Ð°Ñ Ð²ÐµÑÑÐ¸Ñ - СамаÑÑкий гоÑÑдаÑÑÑвеннÑй аÑÑокоÑмиÑеÑкий ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Кибернетика и информатика<br />
УДК 519.7<br />
ОПТИМИЗАЦИЯ СТОХАСТИЧЕСКОЙ СЛЕДЯЩЕЙ СИСТЕМЫ<br />
С ИСПОЛЬЗОВАНИЕМ ЧИСЛЕННЫХ И ЭВОЛЮЦИОННЫХ МЕТОДОВ<br />
© 2007 М. А. Федорова<br />
Ульяновский государственный университет<br />
В работе исследуется применение эволюционных и численных методов для оптимизации сложных взаимосвязанных<br />
систем фильтрации и управления в условиях априорной неопределенности на примере стохастической<br />
следящей системы (для краткости, – трекера). Для сравнения различных подходов проведены серии<br />
экспериментов на специально разработанном программном продукте. При моделировании трекера обнаружение<br />
нарушений производится на основе метода взвешенных квадратов невязок, адаптация – на основе метода<br />
вспомогательного функционала качества, а в качестве алгоритмов идентификации использованы – для сравнения<br />
их возможностей – метод простой стохастической аппроксимации, метод наименьших квадратов и генетический<br />
алгоритм. В сравнительном аспекте исследуется применение эволюционных и численных методов для<br />
оптимизации сложных взаимосвязанных систем фильтрации и управления на примере стохастической следящей<br />
системы.<br />
Постановка задачи<br />
Рассмотрим заданную в пространстве<br />
состояний линейную инвариантную во времени<br />
стохастическую систему с контуром<br />
управления:<br />
x<br />
n<br />
( ti+ 1)<br />
= Φθ x(<br />
ti<br />
) + Ψθ<br />
u(<br />
ti<br />
) + w(<br />
ti<br />
), x ∈ R ,<br />
y( t ) = H v ),<br />
i<br />
(1)<br />
m<br />
θ<br />
x(<br />
ti<br />
) + ( ti<br />
y ∈ R , (2)<br />
−<br />
+<br />
n<br />
xˆ 0<br />
( ti+ 1)<br />
= Φ xˆ<br />
0 0<br />
( ti<br />
) + Ψ0<br />
u(<br />
ti<br />
), xˆ<br />
0<br />
∈ R , (3)<br />
xˆ<br />
( t<br />
0<br />
ν ( t<br />
i<br />
+<br />
i<br />
) = xˆ<br />
( t<br />
0<br />
) = y(<br />
t<br />
i<br />
−<br />
i<br />
) − H<br />
) + K ν ( t<br />
0<br />
xˆ<br />
( t<br />
0<br />
0<br />
−<br />
i<br />
i<br />
),<br />
),<br />
(4)<br />
+<br />
⎪⎧f<br />
[ ˆ<br />
R<br />
x0<br />
( ti<br />
)],<br />
u ( ti<br />
) = ⎨<br />
* +<br />
q<br />
(5)<br />
⎪⎩−G<br />
0<br />
xˆ0<br />
( ti<br />
),<br />
u ∈R<br />
.<br />
Здесь i ∈ Z ; (1) – объект и (2) – сенсор, параметризованные<br />
параметром неопределенности<br />
θ ; (3)-(4) – фильтр Калмана, спроектированный<br />
для некоторого номинального значения<br />
θ<br />
0<br />
параметра θ ; (5) – управление;<br />
{ w ( ⋅)}<br />
, { v ( ⋅)}<br />
считаются независимыми последовательностями<br />
независимых одинаково<br />
распределенных случайных величин с нулевым<br />
средним значением и ковариациями<br />
Q ≥ 0 и R > 0 соответственно.<br />
θ<br />
θ<br />
Матрицы, присутствующие в системе<br />
(1)-(5), заданы как Φ<br />
0<br />
, Ψ<br />
0<br />
, Q<br />
0 , H<br />
0<br />
и R<br />
0<br />
для<br />
номинального режима работы, т.е. для номинального<br />
значения θ ∈Θ<br />
параметра неопре-<br />
0<br />
деленности θ ∈ Θ , взятого из множества Θ<br />
возможных режимов.<br />
Предполагается, что параметр θ подвержен<br />
внезапным изменениям. Каждое изменение<br />
случается в неизвестный момент<br />
времени t > t c 0<br />
. Это событие можно рассматривать<br />
как переключение θ с θ<br />
0<br />
на некото-о-<br />
рое другое неизвестное значение θ<br />
1<br />
∈ Θ. Чтобы<br />
поддерживать обратную связь (ОС) близкой<br />
к оптимальной, для вновь возникшего<br />
режима (определенного параметром θ ) необходимо<br />
соответствующим образом ее пе-<br />
1<br />
ренастроить. Оптимальной перенастройкой<br />
является альтернативный фильтр Калмана<br />
( KF ), которому соответствует<br />
1<br />
θ с коэффи-<br />
1<br />
циентом K . Таким образом, ОС перенастра-<br />
1<br />
ивается и (отмечена нижним индексом 1<br />
)<br />
подставляется вместо начальной обратной<br />
связи (отмечена нижним индексом 0<br />
).<br />
Проблема заключается в том, что оптимальную<br />
перенастройку нельзя выполнить,<br />
253