24.07.2014 Views

Полная версия - Самарский государственный аэрокосмический ...

Полная версия - Самарский государственный аэрокосмический ...

Полная версия - Самарский государственный аэрокосмический ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!