18525

Анализ многошаговой формулы интегрирования Метод простых итераций. Метод ускоренных итераций Итерации Ньютона-Рафсона

Лекция

Информатика, кибернетика и программирование

Лекция 5 Анализ многошаговой формулы интегрирования Метод простых итераций. Метод ускоренных итераций Итерации НьютонаРафсона. Обратные итерации При неявных методах интегрирования ОДУ возникают нелинейные алгебраические уравнения. Возвратимся к общему виду лине...

Русский

2013-07-08

108.5 KB

6 чел.

Лекция 5

Анализ многошаговой формулы интегрирования Метод простых итераций. Метод ускоренных итераций Итерации Ньютона-Рафсона. Обратные итерации

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

Требуется решить неявное уравнение:

                        (1)

т.к.  член под знаком суммы известен, то заменим его на wn, тогда  выражение (1) принимает вид:

                                           (2)

Необходимо определить . Рассмотрим некоторые варианты решения.

  1.  Метод простых итераций (метод Якоби).

Формула метода простых итераций

                                        (3)

Пусть x* - точное решение (2), тогда

                                        (4)

 Вычитая, получим

                                         (5)

  

Используя теорему о среднем, получим:

                                   (6)

где  

По условию Липшица  тогда

                                      (7)

      По индукции

                                 (6)  

Принимая во внимание теорему о единственности решения, необходимое и достаточное решение о сходимости итерационного процесса Якоби имеет вид:

                                                             (7)

   

Т.к. L £ | lmax|, наибольшее собственное значение матрицы -.                                                                                                                 

Если условие удовлетворено, то итерации Якоби сходятся к единственному решению.

                                                                (8)

Для быстрой сходимости необходимо потребовать :

                                                           (9)

Границы зависят от h, если | lmax | велико, то h должно быть очень мало. Для определения условия окончания итераций рассмотрим случай одного уравнения:

  1.  Метод ускоренных итераций

Метод ускоренных итераций – модификация метода итераций Якоби

                           (10)

где a - параметр ускорения.

Если a=0, то получаем простые итерации.

Условие сходимости введем тем же путем.

Точное решение

                        (11)

Вычитая (11) из (10) и пользуясь теоремой о среднем, получаем:

                       (12)                                             

Условие сходимости

 

  или      ,                    (13)

здесь I – единичная матрица.

3. Итерационный метод Ньютона-Рафсона.

Метод описывается формулой

                   (14)

где An+1(s) – матрица Якоби f `x, оцененная в точке x(tn), однократное применение итерации соответствует решению параметризованной формы. Найдем условие сходимости.

Следуя вышеприведенной последовательности действий, получим:

               (15)

Применение (14) является неэффективной процедурой: необходимо вычислить  на каждой итерации.

4. Обратные итерации.

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

Аналогично можно сформировать уравнения с обратными итерациями в виде:

                                            (16)

которые требуют решения неявных уравнений.

Следуя обычной процедуре, запишем:

Условие сходимости:

    или             –  нижняя граница на h.

Краткие выводы:

Итерации Якоби и ускоренные итерации легко реализуются, но сходимость зависит от максимального собственного значения матрицы Якоби. Если |lmax | велико, то шаг мал.

Условиям сходимости метода Ньютона посвящено много литературы. Итерации Ньютона имеют большую область сходимости, чем простые и ускоренные. Зато обратные итерации имеют громадную область сходимости из-за наличия нелинейной границы на h, но существует проблема решения неявных уравнений.

     Рекомендации:

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

    1. Если корректирующая формула в методе не итерируется, то устойчивость метода зависит как от предсказывающих, так и от корректирующих формул.

    2. Если корректирующая формула итерируется, то нет уверенности, что устойчивость зависит от корректирующей формулы.

 


 

А также другие работы, которые могут Вас заинтересовать

2068. Типы взаимодействия генов 39.77 KB
  Различают следующие основные типы взаимодействия генов: 1) комплементарность, 2) элистаз, 3) полимерия, 4) модифицирующее действие генов.
2069. Роль ДНК как материального носителя наследственности 38.18 KB
  Доказательства роли ДНК как материального носителя наследственности. Структура ДНК, объясняющая ее роль как материального носителя наследственности.
2070. ДНК – основной материальный носитель наследственности 39.4 KB
  Доказательства важнейшего генетического значения ДНК были получены в результате анализа многих факторов и специально поставленных опытов.
2071. ДНК и вирусы 38.46 KB
  Вирусы — это внутриклеточные паразиты животных, растений и бактерий.
2072. Основы работы в AutoCad 1.31 MB
  Рабочее окно графического редактора. Видовые экраны. Двумерные графические примитивы. Редактирование чертежа. Положительные аспекты применения блоков. Редактирование трехмерных поверхностей. Твердотельные объекты и команды их редактирования.
2073. Прокат автомобилей. База данных БД 1.61 MB
  Создание таблиц и заполнение их информацией. Определение условий целостности данных. Обеспечение защиты от несанкционированного доступа. Создание меню проекта и программной оболочки. Инфологическая модель БД.
2074. Проектирование реконструкции участка существующей железной дороги 347.35 KB
  Район проектирования реконструкции участка железной дороги. Характеристика существующего состояния железнодорожной линии. Верхнее строение пути, земляное полотно и искусственные сооружения. Определение существующей наличной провозной способности и её соответствие потребным размерам перевозок.
2075. Арт-терапия с элементами рисунка на тему: Наркотикам – НЕТ! 34 KB
  Узнать отношение к проблеме употребления наркотических веществ. сформировать убеждение в том, что каждый человек несёт персональную ответственность за своё решение. сформировать представления о том, что употребление наркотиков не решает, а усугубляет жизненные проблемы. Место проведения: Учебный класс
2076. Геофизические исследования и работы в скважинах 18.15 MB
  Промышленная геофизика. Исследования геологического разреза скважин. Контроль технического состояния скважин. Вторичное вскрытие пластов и специальные операции. Программно-управляемые и информационно-измерительные системы для ГИРС.