Курсавая. Курсавая 1. Методы решения алгебраических и трансцендентных уравнений
Скачать 0.64 Mb.
|
) 12 f(Zn) f'(Zn) f' ' (Zn ) f(Z ) 1 h2[ f'(Zn)]2 f(Zn) f''Zn)] 0. n f'(Zn) 2 Откуда получаем формулу h 2 f(Zn) f'(Zn) 2[ f'(Zn)]2 f(Zn) f''(Zn) . Для уточнения корня можно воспользоваться методом Эй-ткена – Стеффенсена, который строится как сочетание методов хорд и простой итерации с итерационным процессом вида Можно свести решение системы линейных уравнений n aijxj j1 bi, i 1, n к задаче минимизации (до нуля) функции n n 2 gi aijxj bi , gi 0, i 1, n. i1 j1 В самом деле, если решение (2.32) существует и оно един- ственно, то (2.33) обращается в нуль. В случае противоречивой системы достижение нуля нереально, и существенную роль игра- ет выбор весовых коэффициентов gi, выражающих степень зна- чимости соблюдения того или иного уравнения (2.32). Очевидно, что решение задачи оптимизации не проще использования других рассмотренных здесь методов, требует знакомства с соответ- ствующими методами оптимизации и, может быть, даже с мето- дами случайных испытаний (методами Монте-Карло). Мы не пытаемся рассмотреть все многообразие прямых и итерационных методов решения линейных алгебраических си- стем, большинство из которых потеряло свою практическую зна- чимость в век могущества вычислительной техники. Заметим только, что для систем высокого порядка, где матрица коэффици- ентов может быть разбита на блоки, среди которых есть нулевые, используют клеточные методы, сводящие исходную задачу к по- следовательному перебору клеток (блоков) матрицы и поиску со- ответствующих обратных матриц. Особого внимания заслуживают и «плохообусловленные» системы, для которых определитель матрицы близок к нулю. Здесь малейшее изменение исходных данных ведет к значитель- ному изменению решения (решение неустойчиво по исходным данным). В принципе может обнаружиться и отсутствие решения − несовместность системы. Разрешение такой ситуации, связан- ной обычно с некорректной постановкой задачи, предлагает ме-тод регуляризации А. Н. Тихонова*, сводящий исходную систему AX= Bк всегда совместной системе (ATA+ E) X= ATB, > 0, где параметр выбирается с учетом требуемой точности решения и погрешности исходных данных. К решению системы ATA X = ATB сводится и случай пере-определенных систем, в которых число уравнений m превышает число переменных n (m> n) |