Его решение назовем приближением x1 решения уравнения (3), т. е.
![]()
Затем то же проделаем с x1 и т. д. Общая формула для построения последовательных приближений по методу Ньютона такова:
(4)
Сравнение с уравнением (Д.28) показывает, что метод Ньютона состоит в применении метода итераций к уравнению (3), переписанному в равносильной форме:

Почему полезна именно такая форма? Для ответа, обозначив правую часть
вычислим производную:

Таким образом, для точного решения
уравнения (17) получаем при
что
А это, в силу п. 4 Добавления, означает, что если только х0 не слишком далеко от
то метод Ньютона сходится со сверхгеометрической скоростью. (Если
что бывает весьма редко, то метод сходится со скоростью геометрической прогрессии.)
В качестве примера рассмотрим то же уравнение (Д. 29), что и в п. 4 Добавления. Здесь формула (4) приобретает вид
![]()
Начав с x0=0, при вычислениях с точностью до 10-6 получаем
(проверьте!). Таким образом, сходимость получилась существенно более быстрая, чем в п. 4 Добавления. Расхождение в последней значащей цифре с п. 4 Добавления объясняется округлениями при вычислениях.
Покажем еще, как применяется метод Ньютона к системе конечных уравнений на примере системы двух уравнений общего вида
(5)
Пусть мы уже имеем некоторое приближение![]()
![]()
Чтобы перейти к следующему приближению, разложим функции f u g около значений
в ряды Тейлора и отбросим в разложениях все нелинейные члены. Тогда взамен (5) мы получим систему уравнений первой степени

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

Таким образом, исходя из нулевого приближения ![]()
мы можем, положив п = 0, найти первое приближение
и т. д. Из последней системы сразу видно также, что если последовательные приближения сходятся при
то в пределе получается одно из решений системы (5). На практике обычно факт сходимости или расходимости распознается легко, так как сходимость в методе Ньютона, если она имеет место, происходит с весьма высокой скоростью.
Метод Ньютона распространяется и на нелинейные уравнения других типов, в частности на краевые задачи для нелинейных дифференциальных уравнений, путем сведения их к последовательному решению линейных задач.
Имеются и другие способы линеаризации уравнений и моделей.
3.4. 5. Детерминированные и вероятностные модели. Другие типы моделей
Математическая модель может включать случайные компоненты — случайные скалярные или векторные величины, случайные последовательности или функции, случайные структуры и т. п., удовлетворяющие статистическим законам. Такие модели называются вероятностными или стохастическими, в отличие от детерминированных моделей, которые таких компонентов не содержат. Так, если какой-либо элемент изучаемого объекта является изделием массового производства и на интересующие нас свойства могут заметно повлиять отклонения параметров от их номинальных значений, то эти параметры часто считают случайными величинами. Случайные функции появляются, например, при рассмотрении воздействия ветра на какие-либо сооружения, сигналов на фоне шума, шероховатых поверхностей, турбулентных движений жидкости и т. д.
Вероятностные модели изучаются с помощью методов теории вероятностей. К сожалению, довольно часто бывает, что вероятностные характеристики случайных компонентов (математические ожидания и дисперсии случайных величин, тем более законы распределения последних, а также аналогичные характеристики случайных функций) оказываются известными с весьма невысокой точностью или даже вовсе неизвестными, т. е. модель не удовлетворяет требованию продуктивности. Методы математической статистики направлены на определение таких характеристик, но и эти методы не всегда удается эффективно применить. Поэтому при построении вероятностных моделей надо уделять существенное внимание источнику таких характеристик. Если они не поддаются определению с необходимой точностью, то можно попытаться поискать другую модель, быть может более грубую, но и более устойчивую относительно пробелов в знании исходных данных. Например, иногда удается провести исследование и вычисления по максимальным отклонениям рассматриваемых параметров.
Приведем пример. Пусть х — решение задачи Коши

где а — случайная функция (переменная ω, как это принято в теории вероятностей, имеет смысл элементарного исхода). Тогда и
— случайная функция, характеристики которой существенно зависят от характеристик функции а. Но пусть не представляется возможным детально узнать характеристики функции а, известно только, что всегда
Тогда, подставляя крайние возможные значения, мы получаем гарантированную оценку решения:![]()
из нее, например, следует, что
с экспоненциальной скоростью при ![]()
Применяется классификация моделей и по другим признакам. Так, различают статические и динамические (эволюционные) модели; для второго типа моделей предметом изучения является изменение рассматриваемого объекта во времени. Промежуточное место занимают квазистатические, стационарные и квазистационарные модели. В квазистатической модели принимается, что изменение объекта происходит столь медленно, что при рассмотрении ситуации в каждый момент можно в первом приближении объект считать статическим (грубо говоря, пренебречь инерционными силами), а время считать добавочным параметром. В стационарной модели считается, что процессы происходят, но изучаемый объект во времени не меняется; простейший пример — электрическая цепь с постоянным током. Естественно определяется и квазистационарная модель.
В связи с перечисленными сейчас типами моделей упомянем еще о применяемых здесь терминах: установившимся процессом обычно называют стационарный или периодический процесс; переходным процессом называют процесс перехода от одного статического состояния или установившегося процесса к другому.
3.5. Модули, альбомы и многообразия математических моделей
3.5.1. Модули и альбомы математических моделей
В этой разделе k — поле, полное относительно нетривиального абсолютного значения.
Пусть Х — топологическое пространство математических моделей. Модулем математической модели ( в дальнейшем – модулем) с пространства X назывем тройку
где
(1) U — открытое подмножество в X,
(2) п — целое неотрицательное число,
(3) φ —отображение U в kп, причем множество φ(U) открыто в kп и отображение
— гомеоморфизм.
Обозначения:
U = О (с) ~ открытое множество модуля с;
φ = М (с) — отображение модуля с;
— размерность модуля с.
Пусть заданы два модуля
и ![]()
|
Из за большого объема этот материал размещен на нескольких страницах:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 |


