Jump to content

Жесткое уравнение

(Перенаправлено из «Жесткость (математика)

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

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

Мотивирующий пример

[ редактировать ]
Явные численные методы, демонстрирующие неустойчивость при интегрировании жесткого обыкновенного дифференциального уравнения

Рассмотрим задачу начального значения

( 1 )

Точное решение (показано голубым цветом):

( 2 )

Мы ищем численное решение , которое демонстрирует такое же поведение.

Рисунок (справа) иллюстрирует численные проблемы для различных числовых интеграторов, применяемых в уравнении.

  1. Метод Эйлера с шагом сильно колеблется и быстро выходит за пределы диапазона графика (показано красным).
  2. Метод Эйлера с половиной шага, , дает решение в пределах границ графика, но колеблется около нуля (показано зеленым).
  3. Трапециевидный метод (т. е. двухэтапный метод Адамса – Моултона ) имеет вид
    ( 3 )
    где . Применение этого метода вместо метода Эйлера дает гораздо лучший результат (синий). Численные результаты монотонно убывают до нуля, как и точное решение.

Одним из наиболее ярких примеров жестких обыкновенных дифференциальных уравнений (ОДУ) является система, описывающая химическую реакцию Робертсона: [1]


( 4 )

Если рассматривать эту систему на коротком интервале, например, проблем с численным интегрированием нет. Однако если интервал очень велик (10 11 скажем), то многие стандартные коды не могут правильно его интегрировать.

Коэффициент жесткости

[ редактировать ]

Рассмотрим линейную неоднородную систему с постоянными коэффициентами

( 5 )

где и является постоянной, диагонализируемой, матрица с собственными значениями (предполагается различными) и соответствующие собственные векторы . Общее решение ( 5 ) принимает вид

( 6 )

где являются произвольными константами и является частным интегралом. Теперь предположим, что

( 7 )

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

Устный перевод быть временем (как это часто бывает в физических проблемах), называется переходным решением и стационарное решение . Если велико, то соответствующий член быстро распадется, так как увеличивается и поэтому называется быстрым переходным процессом ; если мал, соответствующий член затухает медленно и называется медленным переходным процессом . Позволять определяться

( 8 )

так что является самым быстрым переходным процессом и самый медленный. Теперь мы определяем коэффициент жесткости как [2]

( 9 )

Характеристика жесткости

[ редактировать ]

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

Дж. Д. Ламберт определяет жесткость следующим образом:

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

Существуют и другие характеристики, которые демонстрируются во многих примерах жестких задач, но для каждого из них есть контрпримеры, поэтому эти характеристики не могут служить хорошим определением жесткости. Тем не менее, определения, основанные на этих характеристиках, широко используются некоторыми авторами и являются хорошим ключом к определению наличия жесткости. Ламберт называет их «утверждениями», а не определениями по вышеупомянутым причинам. Некоторые из них:

  1. Линейная система постоянных коэффициентов является жесткой, если все ее собственные значения имеют отрицательную действительную часть и коэффициент жесткости велик.
  2. Жесткость возникает, когда длина шага ограничивается требованиями стабильности, а не требованиями точности.
  3. Жесткость возникает, когда некоторые компоненты раствора распадаются гораздо быстрее, чем другие. [3]

Этимология

[ редактировать ]

Происхождение термина «жесткость» точно не установлено. По словам Джозефа Окленда Хиршфельдера , термин «жесткий» используется потому, что такие системы соответствуют плотной связи между приводом приводными сервомеханизмами и . [4] По словам Ричарда. Л. Берден и Дж. Дуглас Фейрс,

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

. . .

Проблемы, связанные с быстро затухающими переходными решениями, естественным образом возникают в самых разных приложениях, включая исследование пружин и систем демпфирования, анализ систем управления и проблемы химической кинетики . Все это примеры класса задач, называемых жесткими (математической жесткостью) системами дифференциальных уравнений из-за их применения при анализе движения пружинных и массовых систем, имеющих большие константы пружины (физическая жесткость ). [5]

Например, задача начального значения

( 10 )

с , , , можно записать в виде ( 5 ) с и

( 11 )

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

( 12 )

что довольно велико. Тогда система ( 10 ) заведомо удовлетворяет утверждениям 1 и 3. Здесь жесткость пружины велика, а константа затухания еще больше. [6] (хотя «большой» не является четко определенным термином, но чем больше вышеуказанные величины, тем более выраженным будет эффект жесткости.)Точное решение ( 10 ) есть

( 13 )

Уравнение 13 ведет себя очень похоже на простую экспоненту. , но наличие члена, даже с небольшим коэффициентом, достаточно, чтобы сделать численные вычисления очень чувствительными к размеру шага. Стабильное интегрирование ( 10 ) требует очень маленького шага до достижения гладкой части кривой решения, что приводит к ошибке, намного меньшей, чем требуется для точности. Таким образом, система также удовлетворяет утверждению 2 и определению Ламберта.

А-стабильность

[ редактировать ]

Поведение численных методов решения жестких задач можно проанализировать, применив эти методы к тестовому уравнению при условии начального состояния с . Решение этого уравнения есть . Это решение приближается к нулю, поскольку когда Если численный метод также демонстрирует такое поведение (при фиксированном размере шага), то метод называется A-стабильным. [7] Численный метод, который является L-стабильным (см. ниже), обладает более сильным свойством: решение приближается к нулю за один шаг, когда размер шага стремится к бесконечности. А-стабильные методы не вызывают проблем нестабильности, описанных в мотивирующем примере.

Методы Рунге-Кутты

[ редактировать ]

Методы Рунге – Кутты, примененные к проверочному уравнению принять форму , и по индукции . Функция называется функцией устойчивости . Таким образом, условие, как эквивалентно . Это мотивирует определение области абсолютной устойчивости (иногда называемой просто областью устойчивости ), которая представляет собой множество . Метод A-стабилен, если область абсолютной устойчивости содержит множество , то есть левая полуплоскость.

Пример: методы Эйлера

[ редактировать ]
Розовым диском показана область устойчивости метода Эйлера.

Рассмотрим методы Эйлера, описанные выше. Явный метод Эйлера, примененный к тестовому уравнению является

Следовательно, с . Таким образом, область абсолютной устойчивости для этого метода равна это диск, изображенный справа. Метод Эйлера не является A-стабильным.

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

Пример: Трапециевидный метод

[ редактировать ]
Розовая область — это область устойчивости трапециевидного метода.

Рассмотрим трапециевидный метод

применительно к проверочному уравнению , является

Решение для урожайность

Таким образом, функция устойчивости равна

и область абсолютной устойчивости

Эта область содержит левую полуплоскость, поэтому трапециевидный метод A-стабилен. Фактически область устойчивости идентична левой полуплоскости, и, следовательно, численное решение сходится к нулю тогда и только тогда, когда сходится точное решение. Тем не менее трапециевидный метод не обладает идеальным поведением: он затухает все затухающие компоненты, но быстро затухающие компоненты затухают лишь очень слабо, поскольку как . Это привело к концепции L-стабильности : метод является L-стабильным, если он A-стабилен и как . Трапециевидный метод является A-стабильным, но не L-стабильным. Неявный метод Эйлера является примером L-стабильного метода. [8]

Общая теория

[ редактировать ]

Функция устойчивости метода Рунге–Кутты с коэффициентами и дается

где обозначает вектор со всеми единицами. Это рациональная функция (один многочлен делится на другой).

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

Функция устойчивости неявных методов Рунге – Кутты часто анализируется с использованием звезд порядка . Звезда порядка для метода с функцией устойчивости определяется как набор . Метод A-стабилен тогда и только тогда, когда его функция устойчивости не имеет полюсов в левой плоскости и его порядковая звезда не содержит чисто мнимых чисел. [9]

Многошаговые методы

[ редактировать ]

Линейные многошаговые методы имеют вид

Применительно к проверочному уравнению они становятся

который можно упростить до

где . Это линейное рекуррентное соотношение . Метод A-стабилен, если все решения рекуррентного соотношения сходятся к нулю, когда . Характеристический полином

Все решения сходятся к нулю при заданном значении если все решения из лежат в единичном круге.

Тогда областью абсолютной устойчивости многошагового метода указанного вида является совокупность всех ради чего все такой, что удовлетворить . Опять же, если это множество содержит левую полуплоскость, многошаговый метод называется A-стабильным.

Пример: метод Адамса – Башфорта второго порядка.

[ редактировать ]
Розовая область представляет собой область устойчивости метода Адамса – Башфорта второго порядка.

Определим область абсолютной устойчивости двухшагового метода Адамса–Бэшфорта.

Характеристический полином

у которого есть корни

таким образом, область абсолютной устойчивости

Этот регион показан справа. Он не включает всю левую полуплоскость (фактически он включает только действительную ось между ), поэтому метод Адамса – Башфорта не является A-стабильным.

Общая теория

[ редактировать ]

Явные многошаговые методы никогда не могут быть A-стабильными, как и явные методы Рунге–Кутты. Неявные многошаговые методы могут быть A-стабильными только в том случае, если их порядок не превышает 2. Последний результат известен как второй Далквиста барьер ; это ограничивает полезность линейных многошаговых методов для жестких уравнений. Примером A-стабильного метода второго порядка является упомянутое выше правило трапеций, которое также можно рассматривать как линейный многошаговый метод. [10]

См. также

[ редактировать ]

Примечания

[ редактировать ]
  1. ^ Робертсон, Х.Х. (1966). «Решение системы уравнений скорости реакций». Численный анализ: введение . Академическая пресса. стр. 178–182.
  2. ^ Ламберт (1992 , стр. 216–217)
  3. ^ Ламберт (1992 , стр. 217–220)
  4. ^ Хиршфельдер (1963)
  5. ^ Бремя и ярмарки (1993 , стр. 314)
  6. ^ Крейциг (1972 , стр. 62–68)
  7. ^ Это определение принадлежит Далквисту (1963) .
  8. ^ Определение L-стабильности принадлежит Эле (1969) .
  9. ^ Определение дано Ваннером, Хайрером и Норсеттом (1978) ; см. также Iserles & Nørsett (1991) .
  10. ^ См. Далквист (1963) .
[ редактировать ]
Arc.Ask3.Ru: конец переведенного документа.
Arc.Ask3.Ru
Номер скриншота №: 9fd7af4b887c4643b954091cbb0312cd__1719136740
URL1:https://arc.ask3.ru/arc/aa/9f/cd/9fd7af4b887c4643b954091cbb0312cd.html
Заголовок, (Title) документа по адресу, URL1:
Stiff equation - Wikipedia
Данный printscreen веб страницы (снимок веб страницы, скриншот веб страницы), визуально-программная копия документа расположенного по адресу URL1 и сохраненная в файл, имеет: квалифицированную, усовершенствованную (подтверждены: метки времени, валидность сертификата), открепленную ЭЦП (приложена к данному файлу), что может быть использовано для подтверждения содержания и факта существования документа в этот момент времени. Права на данный скриншот принадлежат администрации Ask3.ru, использование в качестве доказательства только с письменного разрешения правообладателя скриншота. Администрация Ask3.ru не несет ответственности за информацию размещенную на данном скриншоте. Права на прочие зарегистрированные элементы любого права, изображенные на снимках принадлежат их владельцам. Качество перевода предоставляется как есть. Любые претензии, иски не могут быть предъявлены. Если вы не согласны с любым пунктом перечисленным выше, вы не можете использовать данный сайт и информация размещенную на нем (сайте/странице), немедленно покиньте данный сайт. В случае нарушения любого пункта перечисленного выше, штраф 55! (Пятьдесят пять факториал, Денежную единицу (имеющую самостоятельную стоимость) можете выбрать самостоятельно, выплаичвается товарами в течение 7 дней с момента нарушения.)