Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Видаль П. Нелинейные импульсные системы

.pdf
Скачиваний:
24
Добавлен:
24.10.2023
Размер:
10.96 Mб
Скачать

Доказательство. Выбор 4+1—4 = 1 для произвольного k никак не ограничивает общности доказательства. Исходя из уравнения (6-114), получаем:

Л1/=|1Мх)11-]|х||. (6-115)

Функция AV в соответствии с ограничением являет­ ся отрицательно-определенной, что доказывает кри­ терий.

Теперь следует найти условия, которым должны удо­ влетворять коэффициенты разностного уравнения, соот­ ветствующие условию ограниченности.

Условия ограниченности. Рассмотрим систему

 

 

x(4+i) = # [х (4 )]х (4 ),

 

(6-116)

Н — матрица вида

 

 

 

 

 

 

 

И [х (4)] =

Д

 

(6-117)

Предположим, что существуют положительные кон­

станты Ci,

. . ., сп, такие, что

 

 

 

 

 

 

 

1 для любого X,

(6-118)

или

 

 

 

 

 

 

 

шах

Cj

ha (х) | <■ 1

для любого X.

(6-119)

 

/

сi

 

 

 

 

 

Для

 

i = \

из

вышеуказанных

предположений

каждого

Н(х)х

ограничена

для

любого

х, и

система

(6-116)

асимптотически устойчива.

Действительно, для соответствующе выбранного не­

равенства (6-118) определим норму как

 

 

||х||=шах {c,]Xt|}.

 

(6-120)

i

 

 

 

 

Тогда

 

 

 

 

|Н (х)хЦ — шах

Ci

 

<

 

 

/ = 1

 

 

 

<шах

 

 

<

 

i

 

 

 

 

4=1

 

 

 

 

< mf х { s ^ |hij

^ '

m a x {Cj | jc, | } .

(6-121)

 

 

 

■£l

Неравенство, проверяющее согласно (6-118) условия (6-111), является ограничением для ||Н(х)х||. Отметим, что для этой нормы контуры H=const = co-— кубы, реб­

ра которых равны 2с0с~1.

Для уравнения (6-119) выберем норму

II х Ц=

Е Ci |Xi [.

(6-122)

Тогда

 

 

 

 

 

IН (х) х |= У Ci J ] ha ( х )

 

п

 

п

 

/=1

 

 

 

 

 

=S

-Д S

1

■'ji хэ \

 

<x) 1Ci

 

i=l

/=1

 

 

 

m a x

£

-

g - | / , j ( x ) |

(6-123)

 

t= i

 

 

' i= l

 

являются неравенствами, позволяющими получить иско­ мые результаты.

Заметим, что в зависимости от выбора нормы функ­ ция может быть и не ограничивающей.

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

Система m разностных уравнений. Рассмотрим си­ стему, состоящую из m разностных уравнений первогопорядка:

x i (п 4~ 1)=

axxi (к)

#,-х2 (я)

•••

х т(я);

1

Хт (п + 1) =

а'тх т (п) -f- атх2

т(и) -f

... - f aZx ™(п)-

'

 

 

 

 

 

(6-124)

где а1— функции х,(п),

х 2(п),

,

х т(п) и п.

 

Критерий Калмана

и Бертрама

запишется в

виде

двух неравенств

 

 

 

 

 

 

шах

 

 

 

(6-125>

 

i

 

 

 

 

 

182

max

(6-126)

/=i

 

Б частности, этот критерий формулируется так: сум­ ма модулей всех коэффициентов каждой строки матри­ цы А должна быть меньше единицы (и это для всех строк).

Точно так же для столбцов: сумма модулей всех ко­ эффициентов каждого столбца должна быть меньше еди­ ницы. Найдем условия, сформулированные для строк Ж. Лагасса, С. Мира и И. Севели [Л. 6-29], исходя из функции Ляпунова:

 

 

 

 

 

 

 

т

 

 

 

 

 

 

V=E|jCi(n)|.

(6-127)

 

 

 

 

 

 

 

г=1

 

В случае

системы

двух ""уравнений первого порядка

(если принять обозначения аг-3- = а/)

 

 

 

 

 

 

xn^ =

a\xn-\-a\tjrC,

(6-128)

 

 

 

 

 

Уп+i == ^ гМг ""I- ®а1/п*

(6-129)

Достаточные условия устойчивости будут:

 

K l

+

K

l - i <

0;

1

 

 

I

|+

1% 1— К

0,

1

 

 

 

 

 

 

 

(6-130)

 

или

 

 

 

 

 

 

 

 

laJl +

K

l -

К О ;

1

 

 

J а[ I + 1а 11~

1 <

0 .

/

 

 

 

 

 

 

(6-131)

 

 

 

 

 

 

 

 

Рис. 6-23.

 

Для исследователя важно выбрать, исходя из двух областей устойчивости, задаваемых условиями (6-130) и (6-131), наиболее широкую. На практике соотношения (6-130), например, определяют две области Dt и D2. Га­ рантируется устойчивость, если ни в какой момент пере­ менные состояния системы не выходят за поверхность

183

V= |xn|+ |y n |= const, находящуюся внутри области D, полученной в результате пересечения D\ и Dz (рис. 6-23).

Заметим, что система неравенств (6-125) и (6-126) при выборе Ci— Cj= 1 (наиболее частный и наиболее про­ стой случай) не допускает значений коэффициентов

а1, равных единице.

4.Критерий Венгжина—Видаля [Л. 6-30)

Внекоторых задачах случается, что разностные уравнения m-го порядка имеют нелинейные коэффициен­

ты /ь ..., fm, зависимые от xn+m-i, ■■•, хп и п:

 

Хп+тfl-X-n+m—1"Ь ... ~\~fm,Xn = 0

(6-132)

Нельзя судить об устойчивости такого уравнения,

основываясь на замене переменных

 

 

х 1(п) = х(п + 1);

1

 

......................................

 

(6-133)

< < *)= *!■ "-,> + ч

I

 

по системе m уравнений первого порядка, так как по крайней мере — 1) коэффициентов равны единице.

Следует применить другой критерий, пригодный для всех систем, описываемых разностными уравнениями лю­ бого порядка.

Условие асимптотической устойчивости хп = 0 являет­ ся необходимым условием устойчивости для х(^)=0. Достаточно рассмотреть только асиптотическую устой­ чивость для хп = 0 и исследовать разностное уравнение, определяющее

s(t)=0.

Предположим, что с помощью соответствующей за­ мены переменных положение равновесия приведено к на­ чалу координат и определяется как

 

Хп+т—1= •••

= Хп ~ 0 .

(6-134)

Критерий

Функция

п+т—1

 

 

 

 

 

У (•Х п + т >• »

-Хп) == V = =

S

О

V) |JCft |

 

п+т—2

 

N

 

 

\\X*\If ) -

/t—1

lf m ||- * f t |

2

2

 

N

 

N

 

 

184

является функцией Ляпунова, если для n>N возможно найти положительное число v, при котором удовлетво­ ряется неравенство

v < i - 2 \f\.

(6-136)

Ыг

 

Условие (6-136) является достаточным

условием

асимптотической устойчивости.

 

Доказательство

а) Пусть функция V — всегда положительная и стре­ мится к нулю в начале координат. Действительно, пред­ полагаем

п+т

 

^ = 2

(6-137)

N

Если соотношение (6-136) удовлетворено, то функ­ ция W положительна и потому же функция V тоже по­ ложительна, так как она всегда больше W:

 

п+т

 

1

 

 

п+т — 1

 

 

1 /= Ц 7 +

£

, М|

+

 

2

\fm\\Xk\.

(6-138)

 

п+т

2

 

 

 

 

 

 

 

б)

Пусть

конечная

разность

всегда

отрицательная.

Очевидно,

что

п

 

п +

1

п

 

 

 

 

 

 

 

 

 

 

Д 2

*fc=

2

xk— '^lXh =

Xn+1\

 

(6-139)

 

 

N

 

N

 

N

 

 

 

ДУ == |-Хп+т |(1

V)

|fj ||-^n+m-i | '•••

1fm |-Xn j,

 

 

 

 

 

 

 

 

 

(6-140)

и по общему уравнению

(6-132)

 

 

 

|Xn+m |

|/ l I

I Xn+m—i | " t " • • • " Ь

|/ m | |Xn |>

(6-141)

 

 

 

 

 

 

 

 

 

откуда

 

 

AVsS—■\\xn+m\.

 

 

(6-142)

 

 

 

 

 

Этого достаточно, чтобы доказать, что ДУ отрица тельно.

185

З а м е ч а н и е .

В

неравенство

 

(6-136) необходимо ввести про­

извольное положительное число v.

Рассмотрим следующее уравнение:

 

 

Хп+1

п +

1

=0.

(6-143)

 

 

1 +

 

 

 

 

 

 

 

 

 

 

Очевидно,

что

если х, < 1,

 

хп -» 0;

1

 

 

 

 

 

 

 

если х, > 1,

 

хп—» 1,

(

 

использование

модифицированного

условия

(6-136)

приведет к

 

 

 

^

л + 1

1

 

 

 

0 <

1

--------—

 

= д т + т г

(6-144>

 

 

 

1

п

 

 

 

 

Более строгое условие асимптотической устойчивости должно

быть:

 

 

 

 

 

 

 

 

 

 

 

( л + 1 ) *

> v-

 

(6' 145>

Проведем сравнение областей устойчивости для линейной и не­ линейной систем [Л. 6-31].

Система первого порядка. Условие, характеризуемое неравенством v < l — \fi\, является достаточным усло­ вием асимптотической устойчивости системы, описывае­ мой уравнением

■^и+1 ~ h f i X n = 0.

(6-146)

Это условие, очевидно, идентично необходимым и до­ статочным условиям устойчивости линейных систем.

Система второго порядка. Достаточное условие асимптотической устойчивости системы, описываемой уравнением

xn+i+fiXn+i~\-fzXn— 0,

(6-147)

может быть охарактеризовано неравенством

 

v < l — 1AI — 1Ы -

(6-148)

Для линейного уравнения с постоянными коэффи­ циентами

Хп+2+ в-Хп+1+ Ьхп= 0

(6-149)

необходимые и достаточные условия устойчивости име­ ют вид:

1 + а +

6 > 0 ;

 

1 —•а +

b ]> 0;

(6-150)

1 - 6 >

0.

 

186

Они определяют область устойчивости линейной си­ стемы, показанную на рис. 6-24, где приведена также область асимптотической устойчивости нелинейной си­ стемы (соответствующей линейным уравнениям с пере­ менными коэффициентами и нелинейному уравнению).

Область асимптотической устойчивости нелинейной системы, конечно, меньшая, чем область устойчивости линейной системы. Невозможно увеличить область

устойчивости нелинейной системы таким образом, чтобы она превосходила область устойчивости линейной систе­ мы (в частности, если fi или /2 являются разрывными).

Укажем выражение условий асимптотической устой­ чивости, часто применяемой на практике. Уравнение (6-147) можно записать:

Azxn+ (fi + 1) Axn+ (fi + /2+ 1) Хп — 0.

(6-151)

Условие асимптотической устойчивости для разност­ ного уравнения

А2хп-Ь/дАхп+ F<iXn= 0

(6-152)

запишется как

 

v < l — \ F i 1| — | F 2 F 11.

(6-153)

Системы третьего и высших порядков. Необходимое и достаточное условие устойчивости для линейного урав­ нения

Хп+з + аХп+ 2 + bxn+l + с х п — 0

(6-154)

187

выражается через коэффициенты а, Ь, с с помощью ряда следующих неравенств:

3 - 3 с + а - 6 > 0 ;

0 ;

3

- j а-

с> ■

1 — а-\-Ь — с у - 0;

(6-155)

1

&-f- b-j- с

0;

 

1 — с2+ ас — Ь> 0.

Эти неравенства определяют область асимптотиче­ ской устойчивости, изображенной на рис. 6-25, где так­ же производится сравнение ее с областью асимптотиче­ ской устойчивости нелинейной системы. При увеличении

порядка разностных уравнений увеличивается и область, заключенная между предложенной областью асимпто­ тической устойчивости и областью устойчивости линей­ ной системы.

З а м е ч а н и е . Отметим, что путем изменений соответствующих координат критерии Калмана — Бертрама, Венгжина — Видаля были обобщены Видалем и Лораном [Л. 6-32, 6-33], Радевым [Л. 6-34] и Бубницким [Л. 6-35]. Однако эти критерии вносят вообще лишь не­ большие усовершенствования, и применение матричной формы запи­ си не меняет смысла метода.

5. Критерий Шеа

Большое затруднение при использовании прямого ме­ тода Ляпунова состоит в том, что рассматриваемая функция Ляпунова не дает всей области асимптотиче­ ской устойчивости исследуемой системы.

188

В. И. Зубов [Л. 6-36] для непрерывных систем свел задачу отыскания функции Ляпунова к нахождению ре­ шения дифференциального уравнения в частных произ­ водных, с помощью которого возможно выбрать функ­ цию Ф, удовлетворяющую определенным условиям. Ре­ шение этого уравнения позволяет получить желаемую функцию Ляпунова. Р. О. Шеа [Л. 6-37] обобщил эти работы на случай разностных уравнений.

Предполагается, что исследуемая система описывает­ ся рядом уравнений с независимыми приращениями

X i ( k + l ) = f i [ X i ( k ) ,

X 2 ( k ) , . . . , X m { k ) \ ,

i= l, 2,

...,

m,

(6-156)

 

или

 

 

(6-157)

x(k + \)=f[x(k)l

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

f i

являются

непрерывными

функциями Xj, причем начало координат является точ­ кой равновесия (нулевое решение).

С другой стороны, Шеа предполагает, что матрица А линейной аппроксимации системы обладает собствен­ ными значениями, меньшими единицы в начале коорди­ нат. В этом случае говорят, что имеет место экспоненци­

альная устойчивость.

 

 

квадра­

Теорема 1. Если положительно-определенная

тичная

форма

 

(6-138), первая

разность

которой

ДУ [*(&)] является решением У(х),

определяющем для

всех х область асимптотической устойчивости, и

 

 

S.V[x(k)} = V[x(k + \))V[x{k)] =

 

при V(0) =0, то

= —Ф'[х(Щ 1— V{x{k)}}

 

(6-158)

 

 

 

 

 

 

 

 

Ц ^ й - = 1 + ф(*)-

 

 

Эту

теорему

легко доказать, разделив оба члена

уравнения (6-158)

на — (1— V[x(k)})

и добавив единицу.

Тогда

достаточно

прологарифмировать

полученное со­

отношение и просуммировать равенства

п-го

порядка

с равенствами

(п+ т)-то порядка.

Суммирование при­

водит к получению некоторого ряда и к получению ре­ шения V[x(n)]. Проверка этого решения может быть произведена подстановкой его в уравнение (6-158).

Теорема

2. Пусть D — область

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

устойчивости

тривиального решения

(6-157). Тогда тео­

18&

рема может быть сформулирована так: если х принад­ лежит D, то 0 ^ Е (х )< 1 .

Доказательство этой теоремы основывается на выра­ жении (6-159) для V[x(n)], найденном при различных вычислениях, произведенных для проверки теоремы 1:

 

V

[ x ( « ) ] =

l - I=-------!------------.

(6-159)

 

 

 

 

П (1+Ф [*(*)!)

 

 

Теорема 3.

Если всюду, за исключением начала ко­

ординат

области

Do,

выполняется неравенство

0 <

< V( x )<\

и если на границах этой области

V(x)=\,

то

все точки D0 лежат в D.

У (х )< 1 —ДН(х)

определенно

Действительно,

если

отрицательная

величина, непрерывность коэффициентов

/ в сочетании

с этим замечанием позволяет доказать,

что V[x{k)] является монотонно убывающей для любой последовательности и что ее предел равен нулю.

Теорема 4. Если D является областью, содержащей начало координат, то условие V(x) = \ соответствует гра­ нице асимптотической устойчивости.

Эта теорема позволяет получить необходимые и до­

статочные условия устойчивости.

решение

уравнения

Теорема

5. Приняв

нулевое

(6-157) экспоненциально

устойчивым, доказывают, что

оно асимптотически устойчиво, если, с одной

стороны,

0<Е (д;)<1

для любого конечного хфО и если, с другой

стороны, предел V(х) для

||х||— >-оо равен единице. Это

условие является необходимым и достаточным.

Теорема

6. Если коэффициенты

могут быть выраже­

ны как коэффициенты ряда Тейлора в окрестности на­ чала координат, то функция V(х) может быть разложе­ на в ряд Тейлора, который сходится в некоторой обла­ сти, расположенной вокруг начала координат.

З а м е ч а н и е .

Эти теоремы касаются существования и свойств

функций Ляпунова

V(x), они определяют целую область асимптоти­

ческой устойчивости

для класса нелинейных импульсных систем.

Граница устойчивости определяется из условия V(x) = l. Если нулевое решение экспоненциально устойчиво, то необходимое и до­

статочное условие устойчивости запишется как

0< Е (х ) < 1 {хфО)

при lim V(x) = l, так как ||х||— >-оо.

 

1, 4, 5. Одна­

Наибольший практический интерес имеют теоремы

ко нужно заметить, что Шеа свел задачу выбора

V(x)

к задаче вы­

бора Ф (х).

190

Соседние файлы в папке книги из ГПНТБ