Файл: Элементы математического моделирования в программных средах 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 |
|