Условная устойчивость разностного уравнения третьего порядка в критических случаях | Статья в журнале «Молодой ученый»

Библиографическое описание:

Робочинская А. Я., Шенмаер И. В., Нигматулин Р. М. Условная устойчивость разностного уравнения третьего порядка в критических случаях // Молодой ученый. — 2016. — №25. — С. 113-122. — URL https://moluch.ru/archive/129/34801/ (дата обращения: 16.10.2018).



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

Ключевые слова: разностное уравнение третьего порядка, условная устойчивость, пространство начальных условий, область устойчивости

Введение ипостановка задачи

Исследование устойчивости решений дискретных систем (разностных уравнений) является одной из важнейших задач. Известно, что для линейных уравнений эта задача сводится к выяснению расположения относительно единичной окружности на комплексной плоскости корней характеристического полинома. Наиболее сложными для изучения являются так называемые критические случаи (называемые также граничной устойчивостью), когда некоторые корни характеристического полинома на комплексной плоскости попадают на единичную окружность. В теории устойчивости изучение критических случаев напрямую связано с исследованием границы области асимптотической устойчивости в пространстве параметров [1, 3].

В последнее время наряду с исследованием классической устойчивости (по Ляпунову) большое внимание уделяется условной устойчивости [1, 3]. В отличие от непрерывных систем [4], проблема условной устойчивости для дискретных систем изучена мало [1–3]. Цель нашей работы: провести полное исследование условной устойчивости нулевого решения линейного разностного уравнения третьего порядка

,(1)

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

Область асимптотической устойчивости уравнения (1) в пространстве коэффициентов полностью исследована в [3] и изображена на рис. 1. Она представляет собой тело, ограниченное гиперболическим параболоидом и двумя плоскостями , .

Рисунок1.tif

Рис. 1. Область асимптотической устойчивости уравнения (1)

Границу области асимптотической устойчивости уравнения (1) образуют четыре вершины: , , , , пять ребер: AC, AD и CD, AB и CB, две грани ACD и ABC и часть гиперболического параболоида.

Определение. Нулевое решение уравнения (1) называется условно устойчивым относительно подпространства W пространства начальных условий , если любое решение уравнения (1) с начальными условиями из подпространства W ограничено.

Будем говорить, что уравнение (1) условно устойчиво относительно подпространства W, если его нулевое решение условно устойчиво относительно подпространства W.

В данной статье мы полностью описываем все подпространства W начальных условий, относительно которых нулевое решение уравнения (1) является условно устойчивым в критических случаях.

Основные результаты

Как следует из результатов работы [3], на границе области асимптотической устойчивости потеря устойчивости может произойти только в вершинах , , , и на ребрах AB: , CD: (в остальных случаях решения уравнения (1) ограничены при любых начальных условиях). Рассмотрим все эти критические случаи.

Случай 1. Вершина: .

Общее решение уравнения (1) имеет вид

Зададим начальные условия. Тогда, решая систему

относительно , получаем , , . Для ограниченности решений в этом случае необходимо и достаточно, чтобы . В пространстве начальных значений уравнение задает плоскость , проходящую через . Получаем следующее

Утверждение 1. Уравнение (1) условно устойчиво в критическом случае тогда и только тогда, когда начальные условия уравнения (1) принадлежат подпространству .

В случае условной устойчивости каждое решение уравнения (1) имеет вид . Все решения являются периодическими с периодом 2 (см. пример на рис. 2).

http://www5b.wolframalpha.com/Calculate/MSP/MSP6841b6dfc69d6b3ge65000045cc6a15165b30bd?MSPStoreType=image/gif&s=63http://www5b.wolframalpha.com/Calculate/MSP/MSP3181iic047a1af6g289000060h4a642213e7agi?MSPStoreType=image/gif&s=58

Рис. 2. Примеры решений с начальными условиями, удовлетворяющими (слева) и не удовлетворяющими (справа) уравнению плоскости

На рисунке 2 слева показана условная устойчивость для начальных значений , , удовлетворяющих уравнению плоскости (решение имеет вид ). Справа показана неустойчивость для начальных значений , , , не удовлетворяющих уравнению плоскости (решение имеет вид ).

Случай 2. Вершина : .

Общее решение уравнения (1) имеет вид .

Зададим начальные условия. Тогда, решая систему

относительно , получаем , , . Для ограниченности решений в этом случае необходимо и достаточно, чтобы . Получаем аналогичное случаю 1 следующее

Утверждение 2. Уравнение (1) условно устойчиво в критическом случае тогда и только тогда, когда начальные условия уравнения (1) принадлежат подпространству .

В случае условной устойчивости каждое решение уравнения (1) имеет вид . Как и в случае 1, все решения являются периодическими с периодом 2.

Случай 3. Вершина: .

Общее решение уравнения (1) имеет вид

Зададим начальные условия . Тогда, решая систему

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

Получаем следующее

Утверждение 3. Уравнение (1) условно устойчиво в критическом случае тогда и только тогда, когда начальные условия уравнения (1) принадлежат подпространству .

В случае условной устойчивости каждое решение уравнения (1) имеет вид . Все решения являются периодическими с периодом 2 (см. пример на рис. 3).

http://www5a.wolframalpha.com/Calculate/MSP/MSP57071i90527369d71999000025dah1echh8g6548?MSPStoreType=image/gif&s=60http://www4b.wolframalpha.com/Calculate/MSP/MSP680820dba1dgca86e2ah0000637eabd0f5ac075f?MSPStoreType=image/gif&s=64

Рис. 3. Примеры решений с начальными условиями, удовлетворяющими (слева) и не удовлетворяющими (справа) уравнению прямой

На рисунке 3 слева показана условная устойчивость для начальных значений , , удовлетворяющих уравнению прямой (решение имеет вид ). Справа показана неустойчивость для начальных значений , , , не удовлетворяющих уравнению прямой (решение имеет вид ).

Случай 4. Вершина: .

Общее решение уравнения (1) имеет вид .

Зададим начальные условия . Тогда, решая систему

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

Получаем следующее

Утверждение 4. Уравнение (1) условно устойчиво в критическом случае тогда и только тогда, когда начальные условия уравнения (1) принадлежат подпространству .

В случае условной устойчивости каждое решение уравнения (1) имеет вид : все решения являются постоянными (см. пример на рис. 3).

http://www4b.wolframalpha.com/Calculate/MSP/MSP25211ceh39cf0d7b74fg00001dh0b2df87i24h2e?MSPStoreType=image/gif&s=43http://www5a.wolframalpha.com/Calculate/MSP/MSP1576211f94b0ge27c6g000001h8fe166c1b3gb36?MSPStoreType=image/gif&s=40

Рис. 4. Примеры решений с начальными условиями, удовлетворяющими (слева) и не удовлетворяющими (справа) уравнению прямой

На рисунке 4 слева показана условная устойчивость для начальных значений , , удовлетворяющих уравнению прямой (решение имеет вид ). Справа показана неустойчивость для начальных значений , , , не удовлетворяющих уравнению прямой (решение имеет вид ).

Случай 5. Ребро AB: .

Общее решение уравнения (1) имеет вид .

Зададим начальные условия . Тогда, решая систему

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

Утверждение 5. Уравнение (1) условно устойчиво в критическом случае тогда и только тогда, когда начальные условия уравнения (1) принадлежат подпространству .

В случае условной устойчивости каждое решение уравнения (1) имеет вид и обладает свойством: , Если дополнительно, то .

http://www5b.wolframalpha.com/Calculate/MSP/MSP1231bh9ie68he6fc6dh00003dc7800i5866a8a0?MSPStoreType=image/gif&s=62http://www5b.wolframalpha.com/Calculate/MSP/MSP2393218484fg58a67a3a0000191fad7f8ide076i?MSPStoreType=image/gif&s=48

Рис. 5. Примеры решений с начальными условиями, удовлетворяющими (слева) и не удовлетворяющими (справа) уравнению плоскости β

На рисунке 5 слева показана условная устойчивость при для начальных значений , , , удовлетворяющих уравнению плоскости (решение имеет вид ). Справа показана неустойчивость при для начальных значений , , , не удовлетворяющих уравнению плоскости (решение имеет вид ).

Случай 6. Ребро CD: .

Общее решение уравнения (1) имеет вид .

Зададим начальные условия . Тогда, решая систему

относительно , получаем , , . Для ограниченности решений в этом случае необходимо и достаточно, чтобы . Получаем следующее

Утверждение 6. Уравнение (1) условно устойчиво в критическом случае тогда и только тогда, когданачальные условия уравнения (1) принадлежат подпространству .

В пространстве начальных значений для каждого уравнение задает плоскость . В случае условной устойчивости каждое решение уравнения (1) имеет вид и обладает свойством: .

http://www4c.wolframalpha.com/Calculate/MSP/MSP5031ge5f6ee06b1cc8g00004e1iebd831g06fae?MSPStoreType=image/gif&s=63http://www4c.wolframalpha.com/Calculate/MSP/MSP231220cfi4if996a00fi00005eiid3faa0ii5368?MSPStoreType=image/gif&s=55

Рис. 6. Примеры решений с начальными условиями, удовлетворяющими (слева) и не удовлетворяющими (справа) уравнению плоскости γ

На рисунке 6 слева показана условная устойчивость при для начальных значений , , , удовлетворяющих уравнению плоскости (решение имеет вид ). Справа показана неустойчивость при для начальных значений , , , не удовлетворяющих уравнению плоскости (решение имеет вид ).

Сравнение подпространств начальных условий

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

Критический случай (граничные точки)

Геометрическая интерпретация подпространства начальных условий

и

Плоскость:

Прямая :

Прямая :

AB:

Плоскость:

CD:

Плоскость :

Заметим, что все полученные прямые и плоскости проходят через точку .

Выясним взаимное расположение плоскостей , , и прямых , .

  1. Плоскость пересекает по прямой , каноническое уравнение которой , т. е. . Заметим, что при имеем .
  2. Плоскость пересекает по прямой , каноническое уравнение которой , т. е. . Заметим, чтопри имеем.
  3. Рассмотрим прямую и плоскость . Они имеют общую точку и скалярное произведение направляющего вектора прямой на вектор нормали к плоскости равно нулю. Значит .
  4. Рассмотрим прямую и плоскость . Они имеют общую точку и скалярное произведение направляющего вектора прямой на вектор нормали к плоскости равно нулю. Значит
  5. Рассмотрим плоскости и . Координаты векторов нормали к плоскостям не пропорциональны, поэтому плоскость пересекает по прямой Прямая проходит через точку , следовательно, ее каноническое уравнение имеет вид: .
  6. Рассмотрим прямую и плоскость Они имеют общую точку и скалярное произведение направляющего вектора прямой на вектор нормали к плоскости равно нулю. Значит
  7. Рассмотрим прямую и плоскость . При они пересекаются в точке . Если то
  8. Рассмотрим прямую и плоскость . При они пересекаются в точке . Если то
  9. Рассмотрим прямую и плоскость Они имеют общую точку и скалярное произведение направляющего вектора прямой на вектор нормали к плоскости равно нулю. Значит
  10. Рассмотрим прямые и Координаты их направляющих векторов не пропорциональны, значит

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

C:\Users\Равиль\Desktop\untitled.tif

Рис. 7. Взаимное расположение плоскостей (PQRS), (QSU), (PTR) и прямых (QS), (PR).

Выводы

Для всех критических случаев подпространства начальных условий, относительно которых уравнение (1) условно устойчиво, имеет размерность 1 или 2. Причем размерность 1 получается только для граничных точек и . В остальных случаях подпространство начальных условий имеет размерность 2. В общем случае пересечение любой пары подпространств имеет размерность 1, при этом имеют место равенства: , .

Работа поддержана грантом ЮУрГГПУ и КГПУ им. В. П. Астафьева (проект № 16–1022).

Литература:

  1. Diblík, J.; Halfarová, H.; Šafařík, J. Conditional Stability of Weakly Delayed Planar Linear Discrete Systems.Recent Advances in Mathematical and Computational Methods.Proceedings of the 17th International Conference on Mathematics and Computers in Science and Engineering (MACMESE '15). KualaLumpur, Malaysia: WSEASPress, 2015. — P. 111–117.
  2. ParhiN., Tripathy A. K. On the behavior of solutions of a class third order difference equations // Journal of Difference Equations and Applications. — 2002. — V. 8, No. 5. — P. 415–426.
  3. Нигматулин Р. М., Кипнис М. М. Свойства дискретных систем третьего порядка на границе их областей устойчивости // Фундаментальные исследования. — 2015. — № 9–1.– С. 39–43; URL: http://www.fundamental-research.ru/ru/article/view?id=38962 (дата обращения: 20.10.2016)
  4. Чудинов К. М. Частичная и условная устойчивость линейных функционально-дифференциальных уравнений: диссертация... канд. физ.-мат. наук: 01.01.02.- Пермь, 2005.- 98 с.
Основные термины (генерируются автоматически): условная устойчивость, уравнение плоскости, вид, плоскость, решение уравнения, уравнение прямой, ограниченность решений, общее решение уравнения, критический случай, прямая.


Ключевые слова

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

Похожие статьи

Периодические решения разностного уравнения третьего порядка

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

Решение краевой задачи. Неизвестная функция и постоянные...

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

Интерес к исследованию трещинообразования в композитах напрямую связан с проектированием новых высокопрочных материалов.

Касательная. Задачи на касательную | Статья в журнале...

Рассмотрим решение некоторых задач. Задачи, связанные с определением того, является ли прямая у = kx + b касательной к графику функции у = f(x). Можно указать два способа решения таких задач. Находим общие точки графиков, т. е. решаем уравнение f(x) = kx + b...

Решение методом продолжения задач математической физики...

Общее решение этого уравнения имеет следующий вид

Задача Коши: Найти решение уравнение (2) в области , удовлетворяющее начальные условия: , . (3). Решения задачи Коши задается формулой Даламбера

Исследование статической устойчивости Навоийской ТЭС...

статическая устойчивость, квадратичная форма, характеристическое уравнение, электрическая система, система, уравнение, общий случай, прямой метод, условие нарушения, элемент матрицы.

Методы решения нелинейных уравнений

Статья посвящена изучению методов решения нелинейных уравнений, в том числе, с использованием системы автоматизированного проектирования MathCAD. Рассмотрены шаговый метод, методы половинного деления и Ньютона...

Об одной краевой задаче для нагруженного уравнения...

, , , , , В области для уравнения исследуем следующую задачу: Задачи G. Требуется найти функцию , обладающую следующими свойствами

Если выполнены условия , , , , то в области решение задачи G существует и единственно.

Решение дифференциальных уравнений методом...

тогда это уравнения называется дифференциальным уравнением разрешенное относительно производной. Определение 2([1]). Пусть задано уравнение (1) и функция определена в области Г плоскости R2.

Логические продолжения некоторого типа задач на построение...

Использование общего вида уравнения окружности может помогать решению задач более сложного вида с двумя параметрами, что отражено в последней задаче. Задача1. Даны вершины треугольника ABC.

Обсуждение

Социальные комментарии Cackle

Похожие статьи

Периодические решения разностного уравнения третьего порядка

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

Решение краевой задачи. Неизвестная функция и постоянные...

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

Интерес к исследованию трещинообразования в композитах напрямую связан с проектированием новых высокопрочных материалов.

Касательная. Задачи на касательную | Статья в журнале...

Рассмотрим решение некоторых задач. Задачи, связанные с определением того, является ли прямая у = kx + b касательной к графику функции у = f(x). Можно указать два способа решения таких задач. Находим общие точки графиков, т. е. решаем уравнение f(x) = kx + b...

Решение методом продолжения задач математической физики...

Общее решение этого уравнения имеет следующий вид

Задача Коши: Найти решение уравнение (2) в области , удовлетворяющее начальные условия: , . (3). Решения задачи Коши задается формулой Даламбера

Исследование статической устойчивости Навоийской ТЭС...

статическая устойчивость, квадратичная форма, характеристическое уравнение, электрическая система, система, уравнение, общий случай, прямой метод, условие нарушения, элемент матрицы.

Методы решения нелинейных уравнений

Статья посвящена изучению методов решения нелинейных уравнений, в том числе, с использованием системы автоматизированного проектирования MathCAD. Рассмотрены шаговый метод, методы половинного деления и Ньютона...

Об одной краевой задаче для нагруженного уравнения...

, , , , , В области для уравнения исследуем следующую задачу: Задачи G. Требуется найти функцию , обладающую следующими свойствами

Если выполнены условия , , , , то в области решение задачи G существует и единственно.

Решение дифференциальных уравнений методом...

тогда это уравнения называется дифференциальным уравнением разрешенное относительно производной. Определение 2([1]). Пусть задано уравнение (1) и функция определена в области Г плоскости R2.

Логические продолжения некоторого типа задач на построение...

Использование общего вида уравнения окружности может помогать решению задач более сложного вида с двумя параметрами, что отражено в последней задаче. Задача1. Даны вершины треугольника ABC.

Задать вопрос