Файл: Элементы математического моделирования в программных средах MATLAB 5 и Scilab (Андриевский Фрадков).pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 05.04.2024

Просмотров: 429

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Одношаговый прогноз по AR- модели ( п=5)

 

Uob а,

1

.

!I

'

jI " —!.

 

 

 

 

 

«

 

1

/ /

 

 

 

 

 

уЮ\

 

\

 

 

 

 

 

 

 

/7 \

1

 

 

i

1

 

i j f ^ y ( k \ k - l h

 

 

 

 

 

 

 

 

О

335

340

345

350

 

355

360

сут

Рис. 6.1. Прогноз курса акций на один шаг (р = 1).

189

п > 4 -f- 6 точность практически не возрастает, что позволяет рекомендовать этот диапазон значений п.

6.2.Управление синхронизацией систем на основе адаптивных наблюдателей

Как уже говорилось в п. 3.8, одним из эффективных применений хаотических моделей является современная техника телекоммуникаций, в частности мобильная телефония. Надо сказать, что шумоподобные и широкополосные сигналы применялись в системах связи и раньше, см. [38], но недавние исследования показали возможность использования специфических свойств хаотических систем, таких как их синхронизируемость [33].

В настоящем параграфе рассмотрена задача синхронизации двух нелинейных систем (приемника и передатчика или объекта управления и эталонной модели) в условиях неполноты измерений и при неполной информации о параметрах систем. Строится алгоритм адаптивного управления синхронизацией на основе адаптивного набюдателя. Приводятся пример адаптивной синхронизации двух цепей Чуа и результаты компьютерного моделирования системы связи. Рассматриваемая задача интересна еще и тем, что она иллюстрирует подход к построению нелинейной теории сигналов, намеченный в п. 3.9. Изложение основано на материале работы [143], см. также [6].

6.2.1.Общая постановка задачи и схема решения

Традиционная теория передачи и приема сигналов основана преимущественно на линейных моделях передачи (модуляции) и приема (оценивание параметров) сигналов, а также на стохастических моделях шума (помех). Существующие работы по нелинейной модуляции [16] предполагают медленное изменение параметров линейных моделей. Ниже описывается вариант другого подхода, основанного на синхронизации хаотических сигналов. Такой подход вызывает в последние годы растущий интерес специалистов.

Учитывая сложности перехода к нелинейным моделям, которые необходимы для генерации хаотических сигналов, будем, по возможности, упрощать постановку задачи. В частно-

190



сти, для описания передатчика не будем использовать общую нелинейную модель (3.97), а ограничимся классом систем Лурье, представляемых как соединение линейной динамической и нелинейной статической подсистем (3.83), (3.84). Модулируемые (информационные) параметры, согласно рекомендаций п. 4.1, введем в модель линейно. Таким образом, модель передатчика принимается в виде

 

 

 

т

 

 

 

 

 

 

xd = Axd

Y^

Vd = Cxd)

 

 

(6.26)

 

 

 

i=l

 

 

 

 

где xd G TZn -

вектор состояния передатчика; yd € TZ -

скаляр-

ный выход (передаваемый сигнал); в = col (0Ь ... , 0Ш) -

вектор

параметров

передатчика.

Предполагается, что

нелинейно-

сти (fi(-)

(г =

1,2, . . . , т ) матрицы Л, С и вектор

В

известны;

0Ь ... , в т

могут изменяться во времени, так как они

содержат

информацию о сообщении, которое подлежит передаче.

На первом этапе синтеза

приемника будем считать век-

тор параметров в постоянным и пренебрежем шумом в канале связи, т.е. будем считать, что принимается сигнал yd. Таким образом, приемник представляет собой другую динамическую систему, которая строит оценки 6{ (г = 1 , . . . , т ) параметров передатчика на основе наблюдения за передава-

емым сигналом yd(t). Задача состоит в получении

уравнений

приемника

 

z = F(z1yd)1

(6.27)

* =

(6.28)

обеспечивающих сходимость

 

Hm (6(t) - в) = 0,

(6.29)

где 6(t) = col (#i(0> • • • >0m(O) ~ вектор оценок параметров. Предлагаемый ниже приемник относится к классу адап-

тивных наблюдателей и имеет вид:

 

х = Ах+В

, у = Сх,

(6.30)

0г\ = ф{{уа,у)у

i = 0 , l , . . . , m ,

(6.31)

191


где х € Кп - вектор состояния адаптивной модели передатчика; в0 € 71 - вспомогательный настраиваемый параметр.

Синтез алгоритма адаптации (6.31) выполним методом скоростного градиента, см. п. 4.3. Состоянием приемника является z = (я,0о>0ь • • • Естественной вспомогательной целью может служить

lim e(t) = 0,

(6.32)

f—оо

 

где e(t) = x(t) — xd(t) - ошибка наблюдения.

 

Хотя условие (6.32) не является необходимым для обеспечения (6.29), оно подсказывает способ выбора подходящей

целевой функции для синтеза алгоритма адаптации (6.31), а

именно выберем

целевую функцию в виде Q{e)

= 1 т Ре, где

т

положительно определенная п

х п-матрица.

Р = Р > 0 -

Очевидно, достижение цели (6.32) эквивалентно соотноше-

нию lim^oo Q(e{t)) =

0, т.е.

асимптотической

(при

t • оо)

минимизации функции

Q(e).

 

 

 

 

 

Лля решения задачи запишем уравнение ошибки:

 

е = Ае + В

 

+

,

У = Се,

 

(6.33)

где в{ = в{-в{ (г = 1 , . . . , т )

-

ошибки

оценки

параметров;

у = yd — У- В соответствии

с методом

скоростного

градиен-

та вычислим скорость изменения Q(e) в силу системы (6.33):

Q(e, в) = еТ Р i^Ae +

В £

0m{yd) + 0оуj

,

(6.34)

а затем вычислим частные производные

(г = 0,1,... га) -

компоненты вектора скоростного градиента

VqQ} где

в =

= col ( 0 o J u . . . , L ) :

 

 

 

 

= ePB<pi{yd),

i = 1,2,... , m,

 

 

 

 

 

(6.35)

= e

PB(yd-y).

 

 

дв0

УУ

У)

 

 

Выражения (6.35) непригодны для непосредственного использования в алгоритме адаптации, поскольку они зависят от

192


вектора ошибки e(t)) недоступного для измерения. Однако эту трудность можно преодолеть, если потребовать выполнения соотношения РВ = С , равносильного тождеству е РВ = у — удТогда алгоритм адаптации, получаемый методом скоростного градиента, имеет вид:

^ = -7.(2/ -

Vd)4>i{Vd)i i = 1,... ,

га,

(6.36)

= ~~7о(у — yd)2-

 

(6.37)

Кроме того, для применимости алгоритмов скоростного градиента требуется так называемое условие достижимости цели: существование значений настраиваемых параметров 0*, таких что Q(e,0*) < 0 при е ф 0. Это равносильно выполнению

при некоторых Р,

0J матричного неравенства РА0 + А0Р

< О,

где Ао = А + ВвцС. Как показано в [63, 106], разрешимость

ма-

тричных соотношений РА0

+ А0 Р < 0, РВ = С эквивалентна

тому, что

линейная часть

системы

(6.33)

гиперминимально-

-фазовая,

т.е. числитель

функции

W(A)

= C(AIn А)"1

В -

гурвицев многочлен степени п — 1. Это и является

условием

достижения цели (6.32) в системе (6.33), (6.36), (6.37).

 

 

Если, кроме того, вектор-функция f(t)

= со1(у?! (^(0)» • • •

••• , ¥\n(l/<f(0)) ограничена и удовлетворяет условию

постоян-

ного возбуждения

(ПВ): существуют а > 0, Т > 0, такие что

 

 

J f{s)f{s)Tds>aI

 

 

(6.38)

 

 

t

 

 

 

 

 

для всех t > 0, то достигается и исходная цель (6.29). Проведенный теоретический анализ позволяет решать ши-

рокий класс задач адаптивной синхронизации. В качестве примера далее рассмотрена задача синхронизации двух цепей Чуа.

6.2.2.Передача сообщений на основе синхронизации с использованием систем Чуа

Пусть в качестве передатчика и приемника используются системы Чуа. Модель передатчика в безразмерной форме имеет вид

7

Б. Р. Андриевский и др.

193