Ключевые слова:
Цель работы: изучить методы решения нелинейных уравнений с одним неизвестным и апробировать их в опытно-экспериментальной работе.
Задачи работы:
- Проанализировать специальную литературу и выбрать наиболее рациональные способы решения нелинейных уравнений, позволяющие глубоко изучить и усвоить данную тему всем выпускникам средней школы.
- Разработать некоторые аспекты методики решения нелинейных уравнений с применением ИКТ.
- Изучить методы решения нелинейных уравнений:
‒ Шаговый метод
‒ Метод деления пополам
‒ Метод Ньютона
Введение.
Без математической грамотности невозможно успешное освоение методов решения задач по физике, химии, биологии и другим предметам. Весь комплекс естественных наук построен и развивается на базе математических знаний. Например, исследование ряда актуальных задач математической физики приводит к необходимости решения нелинейных уравнений. Решение нелинейных уравнений необходимо в нелинейной оптике, физике плазмы, теории сверхпроводимости и физике низких температур. По этой теме есть достаточное количество литературы, но во многих учебниках и статьях трудно разобраться ученику средней школы. В данной работе рассмотрены методы решения нелинейных уравнений, которые можно использовать при решении прикладных задач физики, химии. Интересным представляется аспект применения информационных технологий к решению уравнений и задач по математике.
Шаговый метод.
Пусть требуется решить нелинейное уравнение вида уравнение F(x)=0. Предположим также, что нам задан некоторый интервал поиска . Требуется найти интервал [а,b] длиной h, содержащий первый корень уравнения, начиная с левой границы интервала поиска.
Рис. 1. Шаговый метод
Решить подобную задачу можно несколькими способами. Шаговый метод является наиболее простым из численных методов решения неравенств, но для достижения большой точности необходимо существенно уменьшить шаг, а это сильно увеличивает время расчётов. Алгоритм решения уравнений с помощью данного метода состоит из двух этапов.
I этап. Отделение корней.
На этом этапе определяются участки, на каждом из которых находится только один корень уравнения. Есть несколько вариантов реализации этого этапа:
- Подставляем значения X (желательно с каким-то достаточно мелким шагом) и смотрим где функция сменит знак. Если функция сменила знак, это значит, что на участке между предыдущим и текущим значением X лежит корень (если функция не меняет характер возрастания/убывания, то можно утверждать, что корень на этом интервале один).
- Графический метод. Строим график и оцениваем на каких интервалах лежит один корень.
- Исследуем свойства конкретной функции.
II этап. Уточнение корней.
На данном этапе значение корней уравнения, определенных ранее, уточняется. Как правило на этом этапе используются итерационные методы. Например, метод половинного деления (дихотомии) или метод Ньютона.
Метод половинного деления
Быстрый и достаточно простой численный метод решения уравнений, основанный на последовательном сужении интервала, содержащего единственный корень уравнения F(x)=0 до того времени, пока не будет достигнута заданная точность Е. Данный метод обычно используется при решении квадратных уравнений и уравнений высших степеней. Однако у данного метода есть существенный недостаток - если на отрезке [а,b] содержится более одного корня, то с его помощью не удастся добиться хороших результатов.
Рис. 2. Метод дихотомии
Алгоритм данного метода следующий:
‒ Определить новое приближение корня х в середине отрезка [а;b]: х=(а+b)/2.
‒ Найти значения функции в точках а и х: F(a) и F(x).
‒ Проверить условие F(a)*F(x)
‒ Перейти к пункту 1 и вновь поделить отрезок пополам. Алгоритм продолжить до того времени, пока не будет выполнено условие |F(x)|
Метод Ньютона
Самый точный из численных методов решения; подходит для решения очень сложных уравнений, но усложняется необходимостью вычисления производных на каждом шаге. заключается в том, что если x n - некоторое приближение к корню уравнения 
, то следующее приближение определяется как корень касательной к функции f(x), проведенной в точке x n .
Уравнение касательной к функции f(x) в точке x n имеет вид:
В уравнении касательной положим y = 0 и x = x n +1 .
Тогда алгоритм последовательных вычислений в методе Ньютона состоит в следующем:
Сходимость метода касательных квадратичная, порядок сходимости равен 2.
Таким образом, сходимость метода касательных Ньютона очень быстрая.
Без всяких изменений метод обобщается на комплексный случай. Если корень x i является корнем второй кратности и выше, то порядок сходимости падает и становится линейным.
К недостаткам метода Ньютона следует отнести его локальность, поскольку он гарантированно сходится при произвольном стартовом приближении только, если везде выполнено условие 
, в противной ситуации сходимость есть лишь в некоторой окрестности корня.
Метод Ньютона (метод касательных) обычно применяется в том случае, если уравнение f(x) = 0 имеет корень , и выполняются условия:
1) функция y= f(x) определена и непрерывна при ;
2) f(a)·f(b) (функция принимает значения разных знаков на концах отрезка [a;b ]);
3) производные f"(x) и f""(x) сохраняют знак на отрезке [a;b ] (т. е. функция f(x) либо возрастает, либо убывает на отрезке [a;b ], сохраняя при этом направление выпуклости);
Смысл метода заключается в следующем: на отрезке [a;b ] выбирается такое число x 0 , при котором f(x 0) имеет тот же знак, что и f""(x 0), т. е. выполняется условие f(x 0)·f""(x) > 0 . Таким образом, выбирается точка с абсциссой x 0 , в которой касательная к кривой y=f(x) на отрезке [a;b ] пересекает ось Ox . За точку x 0 сначала удобно выбирать один из концов отрезка.
Рассмотрим данный алгоритм на конкретном примере.
Пусть нам дана возрастающая функция y = f(x) =x 2– 2, непрерывная на отрезке (0;2), и имеющая f "(x) =2x>0 и f ""(x) = 2> 0 .
В нашем случае уравнение касательной имеет вид: y-y 0 =2x 0 ·(x-x 0). В качестве точки x 0 выбираем точку B 1 (b; f(b)) = (2,2). Проводим касательную к функции y = f(x) в точке B 1 , и обозначаем точку пересечения касательной и оси Ox точкой x 1 . Получаем уравнение первой касательной:y-2=2·2(x-2), y=4x-6. Ox: x 1 =
Рис. 3. Построение первой касательной к графику функции f(x)
y=f(x) Ox через точку x 1 , получаем точку В 2 =(1.5; 0.25) . Снова проводим касательную к функции y = f(x) в точке В 2 , и обозначаем точку пересечения касательной и Ox точкой x 2 .
Уравнение второй касательной: y-2.25=2*1.5(x-1.5), y = 3x - 4.25. Точка пересечения касательной и оси Ox: x 2 = .
Затем находим точку пересечения функции y=f(x) и перпендикуляра, проведенного к оси Ox через точку x 2 , получаем точку В 3 и так далее.
Рис. 4. Построение второй касательной к графику функции f(x)
Первое приближение корня определяется по формуле:
= 1.5.
Второе приближение корня определяется по формуле:
=
Третье приближение корня определяется по формуле:
Таким образом, i -ое приближение корня определяется по формуле:
Вычисления ведутся до тех пор, пока не будет достигнуто совпадение десятичных знаков, которые необходимы в ответе, или заданной точности e - до выполнения неравенства |xi-xi-1|
В нашем случае, сравним приближение, полученное на третьем шаге с реальным ответом. Как видно, уже на третьем шаге мы получили погрешность меньше 0.000002.
Решение уравнения при помощи САПР MathCAD
Для простейших уравнений вида f (x ) = 0 решение в MathСAD находится с помощью функции root .
root(f (х 1 , x 2 , … ) , х 1 , a, b ) - возвращает значение х 1 , принадлежащее отрезку [ a, b ] , при котором выражение или функция f (х ) обращается в 0. Оба аргумента этой функции должны быть скалярами. Функция возвращает скаляр.
Рис. 5. Решение нелинейного уравнения в MathCAD (функция root)
Если в результате применения данной функции возникает ошибка, то это может означать, что уравнение не имеет корней, или корни уравнения расположены далеко от начального приближения, выражение имеет локальные max и min между начальным приближением и корнями.
Чтобы установить причину ошибки, необходимо исследовать график функции f (x ). Он поможет выяснить наличие корней уравнения f (x ) = 0 и, если они есть, то определить приблизительно их значения. Чем точнее выбрано начальное приближение корня, тем быстрее будет найдено его точное значение.
Если начальное приближение неизвестно, то целесообразно использовать функцию solve . При этом если уравнение содержит несколько переменных, нужно указать после ключевого слова solve список переменных, относительно которых решается уравнение.
Рис. 6. Решение нелинейного уравнения в MathCAD (функция solve)
Заключение
В ходе исследования были рассмотрены как математические методы, так и решение уравнений с использованием программирования в САПР MathCAD. Различные методы имеют свои достоинства и недостатки. Следует отметить, что применение того или иного метода зависит от начальных условий заданного уравнения. Те уравнения, которые хорошо решаются известными в школе методами разложения на множители и т. п., не имеет смысла решать более сложными способами. Прикладные задачи математики, важные для физики, химии и требующие сложных вычислительных операций при решении уравнений успешно решаются, например, с помощью программирования. Их же хорошо решать методом Ньютона.
Для уточнения корней можно применять несколько методов решения одного и того же уравнения. Именно это исследование и легло в основу данной работы. При этом легко проследить, какой метод наиболее удачен при решении каждого этапа уравнения, а какой метод на данном этапе лучше не применять.
Изученный материал, с одной стороны, способствует расширению и углублению математических знаний, привитию интереса к математике. С другой стороны, задачи реальной математики важно уметь решать тем, кто собирается приобрести профессии технического и инженерного направления. Поэтому данная работа имеет значение для дальнейшего образования (например, в высшем учебном заведении).
Литература:
- Митяков С. Н. Информатика. Комплекс учебно-методических материалов. - Н. Новгород: Нижегород. гос. техн. ун-т.,2006
- Вайнберг М. М., Треногин В. А. Теория ветвления решений нелинейных уравнений. М.: Наука, 1969. - 527 с.
- Бронштейн И. Н., Семендяев К. А. Справочник по математике для инженеров и учащихся ВТУЗов - М.: Наука, 1986.
- Омельченко В. П., Курбатова Э. В. Математика: учебное пособие. - Ростов н/Д.: Феникс, 2005.
- Савин А. П. Энциклопедический словарь юного математика. - М.: Педагогика, 1989.
- Корн Г., Корн Т. Справочник по математики для научных работников и инженеров. - М.: Наука, 1973.
- Кирьянов Д. Mathcad 15/MathcadPrime 1.0. - С-Пб.: БХВ-Петербург, 2012.
- Черняк А., Черняк Ж., Доманова Ю. Высшая математика на базе Mathcad. Общий курс. - С-Пб.: БХВ-Петербург, 2004.
- Поршнев С., Беленкова И. Численные методы на базе Mathcad. - С-Пб.: БХВ-Петербург, 2012.
Ключевые слова: нелинейные уравнения, прикладная математика, САПР MathCAD, метод Ньютона, шаговый метод, метод дихотомии. .
Аннотация: Статья посвящена изучению методов решения нелинейных уравнений, в том числе, с использованием системы автоматизированного проектирования MathCAD. Рассмотрены шаговый метод, методы половинного деления и Ньютона, приведены подробные алгоритмы применения данных методов, а также проведен сравнительный анализ указанных методов.
Метод Ньютона (также известный как метод касательных) - это итерационный численный метод нахождения корня (нуля) заданной функции. Метод был впервые предложен английским физиком, математиком и астрономом Исааком Ньютоном (1643-1727), под именем которого и обрёл свою известность.
Метод был описан Исааком Ньютоном в рукописи De analysi per aequationes numero terminorum infinitas (лат.О б анализе уравнениями бесконечных рядов), адресованной в 1669 году Барроу , и в работе De metodis fluxionum et serierum infinitarum (лат.Метод флюксий и бесконечные ряды) или Geometria analytica (лат.Аналитическая геометрия) в собраниях трудов Ньютона, которая была написана в 1671 году. Однако описание метода существенно отличалось от его нынешнего изложения: Ньютон применял свой метод исключительно к полиномам. Он вычислял не последовательные приближения x n , а последовательность полиномов и в результате получал приближённое решение x.
Впервые метод был опубликован в трактате Алгебра Джона Валлиса в 1685 году, по просьбе которого он был кратко описан самим Ньютоном. В 1690 году Джозеф Рафсон опубликовал упрощённое описание в работеAnalysis aequationum universalis (лат. Общий анализ уравнений). Рафсон рассматривал метод Ньютона как чисто алгебраический и ограничил его применение полиномами, однако при этом он описал метод на основе последовательных приближений x n вместо более трудной для понимания последовательности полиномов, использованной Ньютоном.
Наконец, в 1740 году метод Ньютона был описан Томасом Симпсоном как итеративный метод первого порядка решения нелинейных уравнений с использованием производной в том виде, в котором он излагается здесь. В той же публикации Симпсон обобщил метод на случай системы из двух уравнений и отметил, что метод Ньютона также может быть применён для решения задач оптимизации путём нахождения нуля производной или градиента.
В соответствии с данным методом задача поиска корня функции сводится к задаче поиска точки пересечения с осью абсцисс касательной, построенной к графику функции .
Рис.1 . График изменение функции
Проведенная в любой точке касательная линия к графику функции определяется производной данной функции в рассматриваемой точке, которая в свою очередь определяется тангенсом угла α (). Точка пересечения касательной с осью абсцисс определяется исходя из следующего соотношения в прямоугольном треугольнике: тангенс угла в прямоугольном треугольнике определяется отношением противолежащего катета к прилежащему катету треугольнику. Таким образом, на каждом шаге строится касательная к графику функции в точке очередного приближения . Точка пересечения касательной с осью Ox будет являться следующей точкой приближения . В соответствии с рассматриваемым методом расчет приближенного значения корня на i -итерации производится по формуле:
Наклон прямой подстраивается на каждом шаге наилучшим образом, однако следует обратить внимание на то, что алгоритм не учитывает кривизну графика и следовательно в процессе расчета остается неизвестно в какую сторону может отклониться график.
Условием окончания итерационного процесса является выполнение следующего условия:
где ˗ допустимая погрешность определения корня.
Метод обладает квадратичной сходимостью. Квадратичная скорость сходимость означает, что число верных знаков в приближённом значении удваивается с каждой итерацией.
Математическое обоснование
Пусть дана вещественная функция , которая определена и непрерывна на рассматриваемом участке. Необходимо найти вещественный корень рассматриваемой функции.
Вывод уравнения основано на методе простых итераций, в соответствии с которым уравнение приводят к эквивалентному уравнению при любой функции . Введем понятие сжимающего отображения, которое определяется соотношением .
Для наилучшей сходимости метода в точке очередного приближения должно выполняться условие . Данное требование означает, что корень функции должен соответствовать экстремуму функции .
Производная сжимающего отображения определяется в следующем виде:
Выразим из данного выражение переменную при условии принятого ранее утверждения о том, что при необходимо обеспечить условие . В результате получим выражение для определения переменной :
С учетом этого сжимающая функция прием следующий вид:
Таким образом, алгоритм нахождения численного решения уравнения сводится к итерационной процедуре вычисления:
Алгоритм нахождения корня нелинейного уравнения по методу
1. Задать начальную точку приближенного значения корня функции , а также погрешность расчета (малое положительное число ) и начальный шаг итерации ( ).
2. Выполнить расчет приближенного значения корня функции в соответствии с формулой:
3. Проверяем приближенное значение корня на предмет заданной точности, в случае:
Если разность двух последовательных приближений станет меньше заданной точности , то итерационный процесс заканчивается.
Если разность двух последовательных приближений не достигает необходимой точности , то необходимо продолжить итерационный процесс и перейти к п.2 рассматриваемого алгоритма.
Пример решения уравнений
по методу Ньютона для уравнения с одной переменной
В качестве примера, рассмотрим решение нелинейного уравнения методом Ньютона для уравнения с одной переменной . Корень необходимо найти с точностью в качестве первого приближения .
Вариант решения нелинейного уравнения в программном комплексе MathCAD представлен на рисунке 3.
Результаты расчетов, а именно динамика изменения приближенного значения корня, а также погрешности расчета от шага итерации представлены в графической форме (см. рис.2).
Рис.2 . Результаты расчета по методу Ньютона для уравнения с одной переменной
Для обеспечения заданной точности при поиске приближенного значения корня уравнения в диапазоне необходимо выполнить 4 итерации. На последнем шаге итерации приближенное значение корня нелинейного уравнения будет определяться значением: .
Рис.3 . Листинг программы в MathCad
Модификации метода Ньютона для уравнения с одной переменной
Существует несколько модификаций метода Ньютона, которые направлены на упрощение вычислительного процесса.
Упрощенный метод Ньютона
В соответствии с методом Ньютона требуется вычислять производную функции f(x) на каждом шаге итерации, что ведет к увеличению вычислительных затрат. Для уменьшения затрат, связанных с вычислением производной на каждом шаге расчета, можно произвести замену производной f’(x n ) в точке x n в формуле на производную f’(x 0) в точке x 0 . В соответствии с данным методом расчета приближенное значение корня определяется по следующей формуле: Модифицированный метод Ньютона
Разностный метод Ньютона
В результате приближенное значение корня функции f(x) будет определяться выражением разностного метода Ньютона:
Двух шаговый метод Ньютона
В соответствии с методом Ньютона требуется вычислять производную функции f(x) на каждом шаге итерации, что не всегда удобно, а иногда практически невозможно. Данный способ позволяет производную функции заменить разностным отношением (приближенным значением):
В результате приближенное значение корня функции f(x) будет определяться следующим выражением:
где
Рис.5 . Двух шаговый метод Ньютона
Метод секущих является двух шаговым, то есть новое приближение определяется двумя предыдущими итерациями и . В методе необходимо задавать два начальных приближения и . Скорость сходимости метода будет линейной.
- Назад
- Вперёд
Для того, чтобы добавить Ваш комментарий к статье, пожалуйста, зарегистрируйтесь на сайте.
2. Метод Ньютона решения систем нелинейных уравнений.
Этот метод обладает гораздо более быстрой сходимостью, чем метод простой итерации. В основе метода Ньютона для системы уравнений (1.1) лежит использование разложения функций
, где 
(2.1)
в ряд Тейлора, причём члены, содержащие вторые и более высокие порядки производных, отбрасываются. Такой подход позволяет решение одной нелинейной системы (1.1) заменить решением ряда линейных систем.
Итак, систему (1.1) будем решать методом Ньютона. В области D выберем любую точку 
и назовём её нулевым приближением к точному решению исходной системы. Теперь функции (2.1) разложим в ряд Тейлора в окрестности точки . Будем иметь
Т.к. левые части (2.2) должны обращаться в ноль согласно (1.1), то и правые части (2.2) тоже должны обращаться в ноль. Поэтому из (2.2) имеем
Все частные производные в (2.3) должны быть вычислены в точке .
(2.3) есть система линейных алгебраических уравнений относительно неизвестных Эту систему можно решить методом Крамера, если её основной определитель будет отличен от нуля и найти величины
Теперь можно уточнить нулевое приближение , построив первое приближение с координатами
т.е.
. (2.6)
Выясним, получено ли приближение (2.6) с достаточной степенью точности. Для этого проверим условие
, 
(2.7)
где 
наперёд заданное малое положительное число (точность, с которой должна быть решена система (1.1)). Если условие (2.7) будет выполнено, то за приближённое решение системы (1.1) выберем (2.6) и закончим вычисления. Если же условие (2.7) выполняться не будет, то выполним следующее действие. В системе (2.3) вместо 
возьмём уточнённые значения
, (2.8)
т.е. выполним следующие действия
. (2.9)
После этого система (2.3) будет системой линейных алгебраических уравнений относительно величин Определив эти величины, следующее второе приближение 
к решению системы (1.1) найдём по формулам
Теперь проверим условие (2.7) 
Если это условие выполняется, то заканчиваем вычисления, приняв за приближённое решение системы (1.1) второе приближение 
. Если же это условие не выполняется, то продолжаем строить следующее приближение, приняв в (2.3) 
Строить приближения нужно до тех пор, пока условие на не будет выполнено.
Рабочие формулы метода Ньютона для решения системы (1.1) можно записать в виде.
Вычислить последовательность
Здесь 
являются решением системы
Сформулируем алгоритм вычислений по формулам (2.11)-(2.13).
1. Выберем нулевое приближение , принадлежащее области D.
2. В системе линейных алгебраических уравнений (2.13) положим 
,а .
3. Решим систему (2.13) и найдём величины 
.
4. В формулах (2.12) положим 
и вычислим компоненты следующего приближения .
5. Проверим условие (2.7) на : (См. алгоритм вычисления максимума нескольких величин.)
6. Если это условие выполняется, то заканчиваем вычисления, выбрав за приближённое решение системы (1.1) приближение . Если же это условие не выполняется, то перейдём к п.7.
7. Положим 
для всех .
8. Выполним п.3, положив 
.
Геометрически этот алгоритм можно записать в виде.
Алгоритм. Вычисления максимума нескольких величин .
Пример . Рассмотрим использование метода Ньютона для решения системы двух уравнений.
Методом Ньютона с точностью до решить следующую систему нелинейных уравнений
, (2.14)
здесь 
. Выберем нулевое приближение 
, принадлежащее области D. Построим систему линейных алгебраических уравнений (2.3). Она будет иметь вид
(2.15)
Обозначим
Решим систему (2.15) относительно неизвестных 
, например методом Крамера. Формулы Крамера запишем в виде
(2.17)
где основной определитель системы (2.15)
(2.18)
а вспомогательные определители системы (2.15) имеют вид
.
Найденные значения подставим в (2.16) и найдём компоненты первого приближения 
к решению системы (2.15).
Проверим условие
, (2.19)
если это условие выполняется, то заканчиваем вычисления, приняв за приближённое решение системы (2.15) первое приближение, т. е. 
. Если условие (2.19) не выполняется, то положим 
, 
и построим новую систему линейных алгебраических уравнений (2.15). Решив её, найдём второе приближение 
. Проверим его на . Если это условие будет выполняться, то за приближённое решение системы (2.15) выберем 
. Если условие на не будет выполняться, положим 
, 
и построим следующую систему (2.15) для нахождения 
и т. д.
Задания
Во всех заданиях требуется:
Составить программу численной реализации метода, согласно предложенному алгоритму.
Получить результаты вычислений.
Проверить полученные результаты.
Задана система двух нелинейных уравнений.
1. 
2. 
3. 
4. 
5. 
6. 
7. 
8. 
9. 
10.
11. 
12. 
13. 
14. 
15. 
.
Глава 3. Численные методы решения систем линейных алгебраических уравнений (СЛАУ).
Цель работы . Знакомство с некоторыми приближёнными методами решения СЛАУ и их численной реализацией на ПК.
Предварительные замечания. Все методы решения СЛАУ обычно разделяют на две большие группы. К первой группе относятся методы, которые принято называть точными. Эти методы позволяют для любых систем найти точные значения неизвестных после конечного числа арифметических операций, каждая из которых выполняется точно.
Ко второй группе относятся все методы, не являющиеся точными. Их называют итерационными, или численными, или приближёнными. Точное решение, при использовании таких методов, получается в результате бесконечного процесса приближений. Привлекательной чертой таких методов является их самоисправляемость и простота реализации на ПК.
Рассмотрим некоторые приближённые методы решения СЛАУ и построим алгоритмы их численной реализации. Приближённое решение СЛАУ будем получать с точностью до , где некоторое очень маленькое положительное число.
1. Метод итерации.
Пусть СЛАУ задана в виде
(1.1)
Эту систему можно записать в матричном виде
, (1.2)
где 
- матрица коэффициентов при неизвестных в системе (1.1), 
- столбец свободных членов, 
- столбец неизвестных системы (1.1).
. (1.3)
Решим систему (1.1) методом итерации. Для этого выполним следующие действия.
Во-первых. Выберем нулевое приближение
(1.4)
к точному решению (1.3) системы (1.1). Компонентами нулевого приближения могут быть любые числа. Но удобнее за компоненты нулевого приближения взять либо нули 
, либо свободные члены системы (1.1) 
Во-вторых. Компоненты нулевого приближения подставим в правую часть системы (1.1) и вычислим
(1.5)
Величины, стоящие слева в (1.5) являются компонентами первого приближения 
Действия, в результате которых получилось первое приближение, называются итерацией.
В-третьих. Проверим нулевое и первое приближения на
(1.6)
Если все условия (1.6) выполняются, то за приближённое решение системы (1.1) выберем, либо , либо всё равно, т.к. они отличаются друг от друга не больше чем на и закончим вычисления. Если хотя бы одно из условий (1.6) не будет выполнено, то перейдём к следующему действию.
В-четвёртых. Выполним следующую итерацию, т.е. в правую часть системы (1.1) подставим компоненты первого приближения и вычислим компоненты второго приближения 
, где
В-пятых. Проверим 
и на , т.е. проверим условие (1.6) для этих приближений. Если все условия (1.6) будут выполнены, то за приближённое решение системы (1.1) выберем, либо , либо всё равно, т.к. они отличаются друг от друга не больше чем на . В противном случае будем строить следующую итерацию, подставив компоненты второго приближения в правую часть системы (1.1).
Итерации нужно строить до тех пор, пока два соседних приближения 
и будут отличаться друг от друга не больше, чем на .
Рабочую формулу метода итерации решения системы (1.1) можно записать в виде
Алгоритм численной реализации формулы (1.7) может быть таким.
Достаточные условия сходимости метода итерации для системы (1.1) имеют вид
1. 
, .
2. 
, 
.
3. 
2. Метод простой итерации.
Пусть система линейных алгебраических уравнений (СЛАУ) задана в виде
(2.1)
Чтобы систему (2.1) решить методом простой итерации, её сначала надо привести к виду
(2.2)
В системе (2.2) 
-ое уравнение представляет собой -ое уравнение системы (2.1), разрешённое относительно –ой неизвестной (
).
Метод решения системы (2.1), состоящий в сведении её к системе (2.2) с последующим решением системы (2.2) методом итерации, называется методом простой итерации для системы (2.1).
Таким образом, рабочие формулы метода простой итерации решения системы (2.1) будут иметь вид
(2.3)
Формулы (2.3) можно записать в виде
Алгоритм численной реализации метода простой итерации для системы (2.1) по формулам (2.4) может быть таким.
Этот алгоритм можно записать геометрически.
Достаточные условия сходимости метода простой итерации для системы (2.1) имеют вид
1. 
, .
2. 
, 
.
3. 
3. Стационарный метод Зейделя.
Метод Зейделя решения СЛАУ отличается от метода итерации тем, что найдя какое-то приближение для -той компоненты, мы сразу же используем его для отыскания следующих 
, 
, …, 
-ой компонент. Такой подход позволяет обеспечить более высокую скорость сходимости метода Зейделя по сравнению с методом итерации.
Пусть СЛАУ задана в виде
(3.1)
Пусть 
- нулевое приближение к точному решению 
системы (3.1). И пусть найдено 
-ое приближение 
. Определим компоненты 
-ого приближения по формулам
(3.2)
Формулы (3.2) можно записать в компактном виде
, 
, 
(3.3)
Алгоритм численной реализации метода Зейделя решения системы (3.1) по формулам (3.3) может быть таким.
1. Выберем , например, 
, 
2. Положим .
3. Для всех вычислим .
4. Для всех проверим условия 
.
5. Если все условия в п.4 будут выполнены, то за приближенное решение системы (3.1) выберем либо , либо и закончим вычисления. Если хотя бы одно условие в п.4 не будет выполнено, перейдем к п.6.
6. Положим и перейдем к п.3.
Этот алгоритм можно записать геометрически.
Достаточное условие сходимости метода Зейделя для системы (3.1) имеет вид 
, .
4. Нестационарный метод Зейделя.
Этот метод решения СЛАУ (3.1) обеспечивает еще более высокую скорость сходимости метода Зейделя.
Пусть каким-либо образом для системы (3.1) найдены компоненты -ого приближения и -ого приближения .
Вычислим вектор поправки
Подсчитаем величины
, (4.2)
Расположим величины 
, в порядке их убывания.
В таком же порядке перепишем уравнения в системе (3.1) и неизвестные в этой системе., : Линейная алгебра и нелинейные ... Руководство для лабораторных работ по ... методические указания для практических работ по для студентов ...
Учебная литература (естественные науки и технические) 2000-2011 цикл опд – 10лет цикл сд – 5 лет
Литература... Естественные науки в целом 1. Астрономия [Текст] : пособие для ... Численные методы : Линейная алгебра и нелинейные ... Руководство для лабораторных работ по ... методические указания для практических работ по дисциплине "Экономика транспорта" для студентов ...
- естественные науки (1)
Учебное пособие... руководство для студентов и преподавателей, предназначенное для использования не только при изучении методов работы ... выработке практических навыков с использованием реальных данных. Методические рекомендации по выполнению зачетной работы по данному...
- естественные науки - физико-математические науки - химические науки - науки о земле (геодезические геофизические геологические и географические науки)
Документ... для студентов естественно - ... работ по дисциплине "Генетика и селекция", посвященных актуальным проблемам этой науки . Систематизирована самостоятельная работа студентов по теоретическому и практическому ... линейного , нелинейного , динамического. Все методы ...
- естественные науки - физико-математические науки - химические науки - науки о земле (геодезические геофизические геологические и географические науки) (7)
Список учебниковОпределитель Еремина в линейной и нелинейной алгебре : линейное и нелинейное программирование: новый метод / Еремин, Михаил... Для студентов и преподавателей геологических специальностей вузов. кх-1 1794549 99. Д3 П 693 Практическое руководство по ...
Решение нелинейных уравнений методом Ньютона
Для решения электроэнергетических задач существует несколько моди-фикаций метода. Они позволяют увеличить скорость сходимости итераци-онного процесса и уменьшить время расчета.
Основное достоинство метода – он обладает быстрой сходимостью.
Идея метода состоит в последовательной замене на каждой итерации расчета исходной нелинейной системы уравнений некоторой вспомогатель-ной линейной системой уравнений, решение которой позволяет получить очередное приближение неизвестных, более близкое к искомому решению (линеаризация ).
Рассмотрим нелинейное уравнение в общем виде:
Искомое решение уравнения – точка, в которой кривая пересекает ось абсцисс.
Задаем начальное приближение неиз-вестной х (0) . Определяем значение функции в этой точке w(х (0)) и проводим касательную к кривой в точке В. Точка пересечения этой касательной с осью абсцисс определяет сле-дующее приближение неизвестной х (1) и т.д.
Разложим уравнение (1) в ряд Тейлора в окрестностях точки х (0) . Рас-смотрим члены разложения, содержащие только 1-ю производную:
(2)
х – х (0) = Δх - поправка к неизвестной. Если определим её, то сможем определить и следующее приближение.
Из (2) определяем поправку 
(3)
Тогда следующее приближение: 
(5)
Аналогично получаем к -е приближения:
Это рекуррентная формула метода Ньютона для решения нелинейных уравнений. Она позволяет определять очередные приближения неизвестных.
Формулу (6) можно получить другим способом из рисунка:
Итерационный процесс сходится, если уменьшается и приближается к 0 . Результат достигнут, если .
Комментарий к геометрической интерпретации
Итерационный шаг метода сводится к замене кривой на прямую, ко-торая описывается левой частью уравнения (2). Она является касательной к кривой в точке . Этот процесс называется линеаризацией . Точка пере-сечения касательной к кривой с осью х дает очередное приближение неиз-вестной . Поэтому этот метод называется методом касательных .
Пример:
Пример:
Для того, чтобы определить этим методом все корни нелинейного урав-нения, нужно любым способом определить приблизительное расположение этих корней и задать начальные приближения в близи них.
Простой способ определения области расположения корней - табуляция .
Итерационный процесс Ньютона не сходится , если начальные приближения выбраны так, что:
Процесс или не сходится или сходится очень плохо.
Метод Ньютона-Рафсона для решения СНАУ
Рафсон показал, что итерационный метод Ньютона, предложенный для решения одного нелинейного уравнения , можно использовать для решения систем нелинейных уравнений.
При этом, для решения систем нелинейных уравнений нужно вместо од-ной неизвестной рассматривать совокупность(вектор) неизвестных :
вместо одной невязки уравнения, рассматриваем вектор невязок уравнений системы:
Одна производная в (6) замещается матрицей производных . Операция деления в (6) замещается умножением на обратную матрицу производных. В этом случае метод Ньютона-Рафсона отличается от метода Ньютона пере-ходом от одномерной задачи к многомерной .
Рассмотрим систему действительных нелинейных алгебраических уравне-ний:
(7)
В матричном виде ее можно записать:
где Х = х 2 – вектор – столбец неизвестных;
w 1 (х 1 , х 2 , … х n)
W = w 2 (х 1 , х 2 , … х n) – вектор-функция.
w n (х 1 , х 2 , … х n)
Пусть 
- начальные приближения неизвестных. Разложим каждое уравнение системы (7) в ряд Тейлора в окрестности точки Х (0)
, то есть выполним приближенную замену исходных нелинейных уравнений линей-ными, в которых сохраняется только 1-я производная (линеаризация). В ре-зультате система уравнений (7) принимает вид:
(9)
В результате получили систему линейных уравнений (линеаризованная система), в которой неизвестными являются поправки . Коэф-фициенты при неизвестных в этой системе – первые производные от урав-нений w j исходной нелинейной системы по всем неизвестным Х i . . Они обра-зуют матрицу коэффициентов – матрицу Якоби :
= 
Каждая строка матрицы состоит из первых производных от очередного урав-нения нелинейной системы по всем неизвестным.
Запишем линеаризованную систему (9) в матричной форме:
(10)
Здесь - вектор невязок уравнений исходной системы. Его эле-менты получаем при подстановке в уравнения нелинейной системы очеред-ных приближений неизвестных;
- матрица Якоби . Ее элементами являются первые частные про-изводные от всех уравнений исходной системы по всем неизвестным;
- вектор поправок к искомым неизвестным. На каждой итерации он может быть записан:
Систему (10) с учетом принятых обозначений можно записать:
(12)
Эта система линейна относительно поправок ΔХ (к) .
Система (13) - линеаризованная система уравнений, которой заменяется исходная СНАУ на каждом шаге итерационного процесса.
Система (13) решается любым известным способом, в результате находим вектор поправок . Затем из (11) можем найти очередные приближения неизвестных:
Т.о. каждый шаг итерационного процесса состоит в решении линейной сис-темы (13) и определении очередного приближения из (14).
Из (11) и (12) можно получить общую рекуррентную формулу (в матричном виде), соответствующую методу Ньютона–Рафсона:
(15)
Она имеет структуру, соответствующую формуле (6).
Формула (15) в практических расчетах используется редко , так как здесь нужно обращать матрицу Якоби (большой размерности) на каждой итерации расчетов. В реальных расчетах поправки определяются в результате решения линейной системы (13).
Контроль завершения итерационного процесса выполняем по вектору невязок:
Это условие должно выполняться для невязок всех уравнений системы.
Алгоритм решения СНАУ методом Ньютона-Рафсона
1. Задание вектора начальных приближений неизвестных .
Задание точности расчета є , других параметров расчета
2. Определение невязок нелинейных уравнений в точке приближения ;
2.3. Определение элементов матрицы Якоби в точке очередного прибли-жения неизвестных ;
2.4. Решение линеаризованной системы (13) любым известным методом. Определение поправок к неизвестным .
2.5. Определение очередного приближения неизвестных в соответ-ствии с (14).
2.6. Контроль завершения итерационного процесса в соответствии с (16). Если условие не выполняется, то возврат к пункту 2.
Примерчик:
Решить СЛАУ методом Ньютона-Рафсона:
(решение Х 1 =Х 2 =2)
Запишем уравнения в виде невязок:
Определяем элементы матрицы Якоби:
Матрица Якоби:
Реализуем алгоритм метода Ньютона-Рафсона:
1) Первая итерация:
Начальные приближения 
Невязки 
Матрица Якоби: 
Линеаризованная система уравнений:
1-е приближение неизвестных:
2) Вторая итерация
3) Третья итерация:
… ……… …… …… …… ……..
Решение систем уравнений установившегося режима методом Ньютона-Рафсона
Нелинейное уравнение установившегося режима в форме баланса мощ-ности для -го узла имеет вид:
(17)
Это уравнение с комплексными неизвестными и коэффициентами. Для того, чтобы такие уравнения вида (17) можно было решать методом Ньюто-на-Рафсона, их преобразуют: разделяют действительные и мнимые части. В результате этого каждое комплексное уравнение вида (17) распадается на два действительных уравнения, которые соответствуют балансу активной и ре-активной мощности в узле:
Здесь -заданные мощности в узле;
Неизвестные составляющие напряжения в узлах. Их нужно
определить в результате расчета.
В правой части уравнений (18) - расчетная суммарная мощность пере-токов в ветвях, подходящих к -му узлу.
Запишем эти уравнения (18) в виде невязок :
Невязки уравнений (19) соответствует расчетному небалансу активной и реактивной мощности в -ом узле.
Невязки описывают режим узла і и являются нелинейными функциями от неизвестных напряжений в узлах . Нужно, чтобы -> 0.
Будем решать методом Ньютона-Рафсона систему 2n уравнений вида (19), то есть для решения задачи расчета установившегося режима электри-ческой сети методом Ньютона - Рафсона нужно:
1) сформировать систему 2n уравнений вида (19) для всех узлов электрической сети, кроме балансирующих;
2) организовать итерационный процесс метода Ньютона-Рафсона
для решения этой системы уравнений. В результате решения
получаем искомые составляющие напряжений в узлах .
Запишем эту систему уравнений в общем виде:
(20)
Получили систему 2 нелинейных уравнений невязок с 2 неизвест-ными, которыми. Неизвестными в ней являются составляющие напряжения - модули и углы .
Для решения системы (20) методом Ньютона-Рафсона нужно составить вспомогательную линеаризованную систему уравнений вида (13), решая ко-торую на каждой итерации, определяем поправки к неизвестным:
(21)
С учетом принятых обозначений система (21) может быть записана:
(22)
где -матрица Якоби, её элементами являются частные производные от уравнений системы (20) по всем неизвестным - составляющим напряже-ний 
Вектор невязок уравнений системы (20). Их значения получаем при подстановке в уравнения очередных приближений неизвестных;
Вектор поправок к неизвестным:
; ΔӨ i
= Ө i (к+1) - Ө i (к)
, ΔU i = U i (к+1) - U i (к) .
Для определения элементов матрицы Якоби применяем аналитическое дифференцирование , т.е. дифференцируем каждое уравнение системы (20) по искомым величинам – углам и модулям напряжений. Чтобы сформировать матрицу Якоби, нужно получить аналитические выражения для производных следующих видов :
1) Производная от уравнения невязки активной мощности го узла по углу напряжения этого же узла: ;
2) Производная от уравнения невязки активной мощности го узла по углу напряжения смежного j- го узла: ;
3) Производная от невязки активной мощности го узла по модулю напряжения этого же узла: ;
4) Производная от невязки активной мощности го узла по модулю напряжения смежного узла: ;
Аналогично определяются ещё четыре вида производных – производные от уравнений невязки реактивной мощности го узла по всем неизвестным:
5) ; 6) ; 7) ; 8) .
С учетом этих производных матрицу Якоби можно записать в общем виде:
(23)
Определим аналитические выражения для производных, дифференци-руя уравнения системы (20) по неизвестным величинам. Они имеют вид:
(24)
Матрица Якоби в общем случае - квадратная матрица, симметричная, размерностью , её элементами являются частные производные от невязок уравнений (небаланса мощностей) по всем неизвестным.
Если узлы не связаны между собой, то соответствующие произ-водные в матрицы матрице Якоби, расположенные вне диагонали, будут равны нулю (аналогично матрице проводимостей) – т.к. в соответствующих форму-лах (24) взаимная проводимость y ij является сомножителем и. y ij =0.
Каждая строка матрицы – это производные от очередного уравнения системы (20).
Наличие в схеме моделируемой сети особых узлов (опорные и балансирую-щие узлы, узлы ФМ) сказывается на структуре системы уравнений устано-вившегося режима и на структуре матрицы Якоби:
1. Для узлов с фиксацией модуля напряжения (ФМ), в которых заданы и неизвестными являются и , из матрицы Якоби исключается стро-ка производных (т.к. Q i не задана, то и уравнение баланса реак-тив-ной мощности (18), (19) составить нельзя) и столбец производных (т.к. модуль напряжения U i известен и он исключается из состава неизвест-ных).
2. Для узлов опорных и балансирующих – соответствующие строки и столбцы матрицы исключаются;
3. Если узлы не связаны непосредственно – соответствующие произ-водные в матрице равны нулю.
Матрицу Якоби можно разбить на четыре блока :
1) - производные от уравнений небаланса активной мощности (20) по углам напряжений;
2) - производные от уравнений небаланса активной мощности по модулям напряжений;
3) - производные от уравнений небаланса реактивной мощности (20) по углам напряжений;
4) - производные от уравнений небаланса реактивной мощности по модулям напряжений.
Это матрицы-клетки частных производных небалансов активной и реактив-ной мощностей по неизвестным углам и модулям напряжений. В общем случае, это квадратные матрицы размерностью n×n.
С учетом этого, матрица Якоби может быть представлена в виде блочной мат-рицы:
Где 
субвектора неизвестных величин.
С учетом этого,Тогда линеаризованную систему уравнений (22) можно запи-сать в ви-де:
. (25)
Решая эту линейную систему уравнений (любым известным методом) на
кКаждой итерации метода, находим поправки к неизвестным , а затем и
очередные приближения неизвестных:
(26)
Очередное приближение неизвестных можно, также, получить с использо-ванием итерационной формулы метода Ньютона-Рафсона, аналогичной (15):
- 
· (27)
Тут требуется обращение матрицы Якоби на каждой итерации – громоздкая вычислительная операция.
Алгоритм решения систем уравнений установившегося режима методом Ньютона - Рафсона
1. Задание начальных значений неизвестных напряжений . В ка-честве начальных приближений принимаем: , т.е. номинальные напряжения узлов;
2. Задание условий расчета: точность ε , предельное количество итера-ций , ускоряющие коэффициенты и др.
3. Определение невязок уравнений в соответствии с уравнениями (20) при очередных приближениях неизвестных;
4. Определение элементов матрицы Якоби в соответствии с (24) при очередных приближениях неизвестных;
5. Решение линеаризованной системы уравнений (25) и определение поправок к неизвестным ;
6. Определение очередных приближений неизвестных в соответствии с (26);
7. Проверка завершения итерационного процесса:
Значения невязок уравнений для всех узлов должны быть меньше задан-ной точности.
Если условие не выполняется, то возврат к пункту 3 и повторение рас-чета при новых приближениях неизвестных.
Существует ряд модификаций метода Ньютона-Рафсона. В том числе:
1. Модифицированный метод Ньютона-Рафсона.
Матрицу Якоби рассчитывают один раз при начальных значениях неизвест-ных. На последующих итерациях она принимается постоянной . Это значи-тельно сокращает объем вычислений на каждой итерации, но увеличивает ко-личество итераций.
2. Разделенный метод Ньютона-Рафсона.
Производные вида очень малы и их значениями можно прине-бречь. В результате, в матрице Якоби остаются два блока - 1-й и 4-й, и сис-тема (25), состоящая из уравнений, распадается на две независимые сис-темы размерностью . Каждая из этих систем решается отдельно от другой. Это приводит к сокращению объема вычислений и необходимой памяти ЭВМ.
