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

Отправьте статью сегодня! Журнал выйдет 23 ноября, печатный экземпляр отправим 27 ноября.

Опубликовать статью в журнале

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

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



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

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

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

Исследование устойчивости решений дискретных систем (разностных уравнений) является одной из важнейших задач. Известно, что для линейных уравнений эта задача сводится к выяснению расположения относительно единичной окружности на комплексной плоскости корней характеристического полинома. Наиболее сложными для изучения являются так называемые критические случаи (называемые также граничной устойчивостью), когда некоторые корни характеристического полинома на комплексной плоскости попадают на единичную окружность. В теории устойчивости изучение критических случаев напрямую связано с исследованием границы области асимптотической устойчивости в пространстве параметров [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 с.
Основные термины (генерируются автоматически): условная устойчивость, уравнение плоскости, вид, плоскость, решение уравнения, критический случай, общее решение уравнения, ограниченность решений, прямая, уравнение прямой.


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

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

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

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

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

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

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

Асимптотика решения бисингулярной задачи на бесконечной прямой с квадратичной особенностью по времени

В работе построено асимптотическое разложение решения задачи Коши для бисингулярной параболического уравнения, в случае, когда решение соответствующего «вырожденного» уравнения имеет полюс второго порядка по времени в начальной точке. Асимптотика реш...

Решение линейных рекуррентных соотношений второго порядка

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

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

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

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

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

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

Рассматривается вопрос о построении приближенного решения линейных обыкновенных дифференциальных уравнений четвертого порядка. Излагаются два метода: метод конечных разностей и дифференциальной прогонки с модификацией матричного варианта.

Псевдопараболическая регуляризация одной граничной обратной задачи для уравнения теплопроводности

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

Построение локально оптимальных систем с использованием проекционного метода

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

Об одной задаче определения правой части линейного дифференциального уравнения четвертого порядка

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

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

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

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

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

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

Асимптотика решения бисингулярной задачи на бесконечной прямой с квадратичной особенностью по времени

В работе построено асимптотическое разложение решения задачи Коши для бисингулярной параболического уравнения, в случае, когда решение соответствующего «вырожденного» уравнения имеет полюс второго порядка по времени в начальной точке. Асимптотика реш...

Решение линейных рекуррентных соотношений второго порядка

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

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

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

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

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

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

Рассматривается вопрос о построении приближенного решения линейных обыкновенных дифференциальных уравнений четвертого порядка. Излагаются два метода: метод конечных разностей и дифференциальной прогонки с модификацией матричного варианта.

Псевдопараболическая регуляризация одной граничной обратной задачи для уравнения теплопроводности

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

Построение локально оптимальных систем с использованием проекционного метода

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

Об одной задаче определения правой части линейного дифференциального уравнения четвертого порядка

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

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