WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     || 2 |
Сибирский математический журнал Май июнь, 2002. Том 43, № 3 УДК 517.962.22 АСИМПТОТИЧЕСКАЯ УСТОЙЧИВОСТЬ РЕШЕНИЙ ВОЗМУЩЕННЫХ ЛИНЕЙНЫХ РАЗНОСТНЫХ УРАВНЕНИЙ С ПЕРИОДИЧЕСКИМИ КОЭФФИЦИЕНТАМИ К. Айдын, А. Я. Булгаков, Г. В. Демиденко Аннотация: Рассматриваются возмущенные линейные системы разностных уравнений y(n + 1) = (A(n) + B(n))y(n), n 0, (1) где {A(n)} T -периодическая матричная последовательность, т. е. A(n + T ) = A(n), n 0, B(n) матрица возмущения. Предполагается, что нулевое решение системы x(n + 1) = A(n)x(n), n 0, асимптотически устойчиво, т. е. все собственные значения матрицы монодромии X(T ) = A(T -1)... A(1)A(0) принадлежат единичному кругу {|| < 1}. Получены условия на возмущение B(n), при которых нулевое решение системы будет асимптотически устойчивым, а также установлена непрерывная зависимость одного класса числовых характеристик асимптотической устойчивости решений системы (1) от коэффициентов системы. Библиогр. 9.

§ 1. Введение В этой работе мы продолжаем рассматривать системы разностных уравнений с периодическими коэффициентами [1]:

x(n + 1) = A(n)x(n), n 0, (1.1) где (N N)-матрица A(n) периодическая с периодом T, т. е.

A(n + T ) = A(n), n 0.

В дальнейшем мы будем предполагать, что нулевое решение системы (1.1) асимптотически устойчиво. Согласно спектральному критерию это означает, что все собственные значения матрицы монодромии системы (1.1):

X(T ) = A(T - 1)... A(1)A(0), принадлежат единичному кругу {|| < 1} (см., например, [2, 3]).

В предыдущей работе [1] мы указали ряд числовых характеристик асимптотической устойчивости решений системы (1.1), не опираясь на спектр матрицы монодромии. Используя эти характеристики, можно получать различные оценки норм решений {x(n)} системы (1.1), из которых вытекает, что норма x(n) Работа выполнена при финансовой поддержке фонда the Scientific and Technical Research Council of Turkey (TUBITAK) и Российского фонда фундаментальных исследований (код проекта 99–01–00533).

© 2002 Айдын К., Булгаков А. Я., Демиденко Г. В.

494 К. Айдын, А. Я. Булгаков, Г. В. Демиденко стремится к нулю при n. В частности, была определена следующая характеристика асимптотической устойчивости решений (1.1):

M(A, T ) = max{ H(0), H(1),..., H(T - 1) }, (1.2) где k-1 k- H(l) = A(j) A(j), l 0, (1.3) k=l j=l j=l и по определению k- A(k - 1)... A(l) при k > l, A(j) = I при k = l.

j=l (Здесь и далее мы используем спектральные нормы матриц.) В работе [1] была доказана следующая Теорема. Для решения {x(n)} начальной задачи x(n + 1) = A(n)x(n), n 0, x(0) = xимеет место n- x(n) 2 1 - H(0) x0 2, n 1.

H(j) j=Из этой теоремы вытекает оценка для решения системы (1.1):

n/ x(n) 1 - H(0) 1/2 x(0), n 0.

M(A, T ) В случае постоянных коэффициентов A(n) A эта оценка совпадает с известной (см., например, [4, 5]). Учитывая, что решение (1.1) имеет вид n- x(n) = A(j) x(0), n 0, j=в силу определения (1.2) из этой оценки получаем n- n/ A(j) 1 - (M(A, T ))1/2.

j=0 M(A, T ) Аналогичным образом доказывается неравенство (n-l-1)/n- A(j) 1 - (M(A, T ))1/2, l n - 1. (1.4) j=l+1 M(A, T ) Отметим, что матричная последовательность {H(l)} является периодической с периодом T и удовлетворяет соотношениям A(l)H(l + 1)A(l) - H(l) = -I, l 0, (1.5) матрицы H(l) эрмитовы положительно определенные, при этом H(l) 1 (см.

[1]).

Асимптотическая устойчивость В частном случае, когда коэффициенты системы (1.1) постоянны, т. е.

A(n) = A, n 0, последовательность {H(l)} стационарна:

H(l) = H = (A)kAk, l 0, k=при этом матрица H решение дискретного уравнения Ляпунова AHA - H = -I.

В этой работе мы будем рассматривать возмущенные линейные системы разностных уравнений y(n + 1) = (A(n) + B(n))y(n), n 0, (1.6) где B(n) (N N)-матрица возмущения. Цель работы получить условия на возмущение в (1.6), при котором нулевое решение системы будет асимптотически устойчивым, а также установить непрерывную зависимость матриц H(l), l 0, от коэффициентов системы разностных уравнений.

§ 2. Периодические возмущения систем разностных уравнений В этом параграфе мы будем рассматривать возмущенные линейные системы разностных уравнений (1.6) с периодическими коэффициентами и, используя матрицы H(l), укажем условия на матрицу возмущения B(n), при которых нулевое решение системы будет асимптотически устойчивым.

Теорема 2.1. Пусть матрица возмущения B(n) периодическая с периодом T и {H(n)} последовательность матриц, построенная для системы (1.1) по формуле (1.3). Предположим, что для T -периодической последовательности {(n)} такой, что A(n)H(n + 1)B(n) + B(n)H(n + 1)A(n) + B(n)H(n + 1)B(n) I + (n)H(n), 0 n T - 1, (2.1) выполнены условия T - 0 1 + (n), (1 + (l)) < 1. (2.2) l=Тогда нулевое решение системы (1.6) асимптотически устойчиво.

Доказательство. Напомним, что матрицы H(l) эрмитовы положительно определенные. Пусть {y(n)} произвольное решение системы (1.6). Используя для скалярного произведения векторов u, v EN обозначение N u, v = uivi, i=рассмотрим форму H(n + 1)y(n + 1), y(n + 1). Учитывая, что {H(l)} удовлетворяет соотношениям (1.5), будем иметь H(n + 1)y(n + 1), y(n + 1) = H(n + 1)(A(n) + B(n))y(n), (A(n) + B(n))y(n) = (A(n)+B(n))H(n+1)(A(n)+B(n))y(n), y(n) = A(n)H(n+1)A(n)y(n), y(n) 496 К. Айдын, А. Я. Булгаков, Г. В. Демиденко + (A(n)H(n + 1)B(n) + B(n)H(n + 1)A(n) + B(n)H(n + 1)B(n))y(n), y(n) = H(n)y(n), y(n) - y(n), y(n) + (A(n)H(n + 1)B(n) + B(n)H(n + 1)A(n) + B(n)H(n + 1)B(n))y(n), y(n).

Тогда в силу условия (2.1) H(n + 1)y(n + 1), y(n + 1) H(n)y(n), y(n) (1 + (n)).

Oтсюда, применяя (2.2), для любого n 0 получим неравенство n- H(n)y(n), y(n) H(0)y(0), y(0) (1 + (l)). (2.3) l=Следовательно, в силу T -периодичности матрицы H(n) и функции (n) для любого натурального k имеем T - H(0)y(kT ), y(kT ) H(0)y(0), y(0) (1 + (l))k.

l=Пусть n- Y (n) = (A(j) + B(j)) j= матрицант системы (1.6). В силу T -периодичности матрицы (A(n) + B(n)) выполняется равенство k y(kT ) = Y (T )y(0).

Тогда предыдущее неравенство переписывается в виде T - k k H(0)Y (T )y(0), Y (T )y(0) H(0)y(0), y(0) (1 + (l))k.

l=Отсюда имеем T - k Y (T )y(0) 2 cond(H(0)) y(0) 2 (1 + (l))k, l=где cond(H(0)) = H(0) H-1(0) число обусловленности матрицы H(0).

Ввиду условия (2.2) T - (1 + (l))k 0, k.

l=Следовательно, из этой оценки вытекает, что все собственные значения j матрицы монодромии Y (T ) лежат в единичном круге {|| < 1}, т. е. нулевое решение системы (1.6) асимптотически устойчиво.

Теорема доказана.

Сейчас мы приведем аналогичную теорему об асимптотической устойчивости решений возмущенной системы разностных уравнений (1.6), в которой вместо условий (2.1), (2.2) на возмущение B(n) будет фигурировать несколько более жесткое, но гораздо проще проверяемое на практике условие.

Асимптотическая устойчивость Теорема 2.2. Пусть матрица возмущения B(n) периодическая с периодом T и {H(n)} последовательность матриц, построенная для системы (1.1) по формуле (1.3). Предположим, что для норм T -периодической последовательности матриц {C(n)}, C(n) = A(n)H(n + 1)B(n) + B(n)H(n + 1)A(n) + B(n)H(n + 1)B(n), (2.4) выполнены неравенства C(n) < 1, 0 n T - 1. (2.5) Тогда нулевое решение системы (1.6) асимптотически устойчиво. Более того, для решений {y(n)} системы (1.6) имеют место оценки n- H(n)y(n), y(n) H(0)y(0), y(0) (1 + (l)), n > 0, (2.6) l=где 1 - C(l) (l) = -.

H(l) Доказательство. Покажем, что для матрицы возмущения B(n) выполнены условия теоремы 2.1.

Очевидно, {(n)} является T -периодической последовательностью. Покажем, что для нее выполнено условие (2.2) из предыдущей теоремы. Действительно, учитывая, что H(n) 1, имеем 1 - C(n) 1, H(n) а в силу (2.5) 1 - C(n) > 0.

H(n) Из этих неравенств следует, что 1 > 1 + (n) 0, n 0.

Отсюда непосредственно вытекает (2.2).

Покажем теперь, что для матричной последовательности {C(n)} и последовательности {(n)} выполнено условие (2.1). Для этого согласно определению (2.4), очевидно, нужно показать справедливость неравенств C(n)v, v v 2 + (n) H(n)v, v (2.7) для любого вектора v EN.

Выпишем следующую оценку:

C(n)v, v C(n) v 2 = v 2 - (1 - C(n) ) v 2.

Учитывая неравенство H(n)v, v H(n) v 2 и условие (2.5), имеем 1 - C(n) C(n)v, v v 2 - H(n)v, v.

H(n) Отсюда в силу определения (n) приходим к (2.7).

Из предыдущего получаем, что для матрицы B(n) выполнены все условия теоремы 2.1. Поэтому нулевое решение системы (1.6) асимптотически устойчиво. А из оценки (2.3), установленной при доказательстве теоремы 2.1, непосредственно вытекает неравенство (2.6).

Теорема доказана.

498 К. Айдын, А. Я. Булгаков, Г. В. Демиденко Следствие 1. Предположим, что T -периодическая матрица возмущения B(n) удовлетворяет условию B(n) < A(n) 2 + - A(n), 0 n T - 1. (2.8) H(n + 1) Тогда нулевое решение системы (1.6) асимптотически устойчиво и для ее решений имеют место оценки n- H(n)y(n), y(n) H(0)y(0), y(0) (1 + (l)), n > 0, (2.9) l=где 1 - (2 A(l) B(l) + B(l) 2) H(l + 1) (l) = - < 0.

H(l) Доказательство. Вначале заметим, что (2.8) эквивалентно неравенствам 2 A(n) B(n) + B(n) 2 <, 0 n T - 1. (2.10) H(n + 1) Действительно, (2.10) можно переписать в виде ( B(n) + A(n) )2 - A(n) 2 + < H(n + 1) или B(n) + A(n) - A(n) 2 + H(n + 1) B(n) + A(n) + A(n) 2 + < 0.

H(n + 1) Отсюда следует эквивалентность (2.8) и (2.10).

Из неравенств (2.10), очевидно, вытекает условие (2.5), поскольку C(n) (2 A(n) B(n) + B(n) 2) H(n + 1) < 1.

Следовательно, в силу теоремы 2.2 нулевое решение системы (1.6) асимптотически устойчиво.

Из оценки (2.6) ввиду последнего неравенства получаем также оценку (2.9).

Следствие 2. Если выполнены условия следствия 1, то для решений системы (1.6) выполнены неравенства y(n) 2 H(0) y(0) 2(1 - )n, n > 0, где = min {|(l)|} (0, 1).

0lT -Приведем теоремы об асимптотической устойчивости для некоторых частных случаев, вытекающие из доказанных утверждений.

Вначале рассмотрим возмущенную линейную систему разностных уравнений с вещественным параметром µ:

y(n + 1) = (A(n) + µB(n))y(n), n 0. (2.11) Из теоремы 2.2 непосредственно следует Асимптотическая устойчивость Теорема 2.3. Пусть матрица возмущения B(n) периодическая с периодом T и матрицы C(n), n 0, определены в (2.4). Если параметр µ такой, что выполнены неравенства |µ| C(n) < 1, n 0, то нулевое решение системы (2.11) асимптотически устойчиво.

Рассмотрим линейные системы разностных уравнений (1.1) и (1.6) с постоянными коэффициентами. Из следствия 1 вытекает Теорема 2.4. Пусть A(n) = A, B(n) = B, n 0, и матрица H решение дискретного уравнения Ляпунова AHA - H = -I, т. е.

H = (A)kAk, l 0.

k=Если B < A 2 + - A, H то нулевое решение системы (1.6) асимптотически устойчиво и для ее решений имеют место оценки Hy(n), y(n) Hy(0), y(0) (1 - )n, n > 0, где 1 - (2 A B + B 2) H = (0, 1).

H Замечание 1. Аналог теоремы 2.4 содержится в [6]. Подобный результат для систем дифференциальных уравнений с постоянными коэффициентами dy = Ay см. в [6, 7].

dt Замечание 2. Аналогичные результаты для систем дифференциальных dy уравнений с периодическими коэффициентами = A(t)y, A(t + T ) = A(t), dt получены в [8].

§ 3. Непрерывная зависимость характеристик асимптотической устойчивости решений разностных уравнений Рассмотрим систему разностных уравнений (1.1):

x(n + 1) = A(n)x(n), n 0, и возмущенную систему (1.6):

y(n + 1) = (A(n) + B(n))y(n), n 0.

Будем предполагать, что матрица возмущения B(n) удовлетворяет условиям теоремы 2.2 и, значит, нулевое решение (1.6) асимптотически устойчиво. Тогда наряду с матричной последовательностью {H(l)}, определенной в (1.3), можно ввести аналогичную T -периодическую последовательность {H(l)} по формулам k-1 k- H(l) = (A(j) + B(j)) (A(j) + B(j)), l 0.

k=l j=l j=l 500 К. Айдын, А. Я. Булгаков, Г. В. Демиденко Матрицы {H(l)} эрмитовы положительно определенные, H(l) 1 и для них выполняются соотношения (A(l) + B(l))H(l + 1)(A(l) + B(l)) - H(l) = -I, l 0. (3.1) По аналогии с (1.2) определим величину M(A + B, T ) = max{ H(0), H(1),..., H(T - 1) }, (3.2) которая будет являться числовой характеристикой асимптотической устойчивости решений возмущенной системы разностных уравнений (1.6). Для решений этой системы, очевидно, справедливо неравенство n/ y(n) 1 - H(0) 1/2 y(0), n 0.

M(A + B, T ) Из общих соображений ясно, что при достаточно малых возмущениях B(n), n 0, величина модуля разности |M(A, T ) - M(A + B, T )| должна быть малой. В этом параграфе мы приведем оценки для норм H(l) - H(l), l 0, из которых будет вытекать, что при малых возмущениях |M(A, T ) - M(A + B, T )| c max { B(j) }.

0jT -Теорема 3.1. Предположим, что для матриц C(n), определенных в (2.4), выполнено условие (2.5), т. е. C(n) < 1, n 0. Тогда имеют место оценки cmax H(l) - H(l) H(l), l = 0, 1,..., T - 1, (3.3) 1 - cmax где cmax = max { C(j) } (0, 1).

0jT - Доказательство. Поскольку для последовательностей {H(l)}, {H(l)} выполняются соотношения (1.5), (3.1) соответственно, то для T -периодической последовательности {V (l)}, V (l) = H(l) - H(l), l 0, имеем равенства (A(l) + B(l))V (l + 1)(A(l) + B(l)) - V (l) = -C(l), l 0, (3.4) где матрицы C(l) определены в (2.4). В следующей лемме мы покажем, что матрицы V (l), l = 0, 1,..., T - 1, являются решениями дискретных уравнений Ляпунова со специальными правыми частями. Этот факт существенно будет использован при получении оценок (3.3).

Введем следующие обозначения:

n- A(n) = A(n) + B(n), Y (n, l) = A(j), n l 0, j=l D(l) = C(l) + Y (l + 1, l)C(l + 1)Y (l + 1, l) + Y (l + 2, l)C(l + 2)Y (l + 2, l) + · · · + Y (l + T - 1, l)C(l + T - 1)Y (l + T - 1, l), l = 0, 1,..., T - 1.

Асимптотическая устойчивость Лемма 3.1. Для каждой матрицы V (l), l = 0, 1,..., T - 1, выполнено соотношение Y (l + T, l)V (l)Y (l + T, l) - V (l) = -D(l). (3.5) Доказательство. В силу (3.4), используя наши обозначения, имеем V (l) = C(l) + A(l)V (l + 1)A(l) = C(l) + A(l)(C(l + 1) + A(l + 1)V (l + 2)A(l + 1))A(l) = C(l) + A(l)C(l + 1)A(l) + (A(l + 1)A(l))V (l + 2)A(l + 1)A(l) = C(l) + Y (l + 1, l)C(l + 1)Y (l + 1, l) + Y (l + 2, l)V (l + 2)Y (l + 2, l).

Если T > 2, то, вновь применяя (3.4) для матрицы V (l + 2), получим V (l) = C(l) + Y (l + 1, l)C(l + 1)Y (l + 1, l) + Y (l + 2, l)(C(l + 2) + A(l + 2)V (l + 3)A(l + 2))Y (l + 2, l) = C(l) + Y (l + 1, l)C(l + 1)Y (l + 1, l) + Y (l + 2, l)C(l + 2)Y (l + 2, l) + Y (l + 3, l)V (l + 3)Y (l + 3, l), и т. д. Наконец, V (l) = C(l) + Y (l + 1, l)C(l + 1)Y (l + 1, l) + Y (l + 2, l)C(l + 2)Y (l + 2, l) + · · · + Y (l + T - 1, l)C(l + T - 1)Y (l + T - 1, l) + Y (l + T, l)V (l + T )Y (l + T, l) = D(l) + Y (l + T, l)V (l + T )Y (l + T, l).

Pages:     || 2 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.