Графическое решение нелинейных уравнений в среде excel

Графический способ решения уравнений в среде Microsoft Excel 2007

Тип урока: Обобщение, закрепление пройденного материала и объяснение нового.

Цели и задачи урока:

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

Оборудование: персональные компьютеры, мультимедиапроектор, проекционный экран.

Материалы к уроку: презентация Power Point на компьютере учителя (Приложение 1).

Слайд 1 из Приложения1 ( далее ссылки на слайды идут без указания Приложения1).

Объявление темы урока.

1. Устная работа (актуализация знаний).

Слайд 2 — Соотнесите перечисленные ниже функции с графиками на чертеже (Рис. 1):

у = 6 — х; у = 2х + 3; у = (х + 3) 2 ; у = -(х — 4) 2 ; .

Слайд 3 Графический способ решения уравнений вида f(x)=0.

Корнями уравнения f(x)=0 являются значения х₁, х₂, … точек пересечения графика функции y=f(x) с осью абсцисс (Рис. 2).

Найдите корни уравнения х 2 -2х-3=0, используя графический способ решения уравнений (Рис.3).

Слайд 5 Графический способ решения уравнений вида f (x)=g (x).

Корнями уравнения f(x)=g(x) являются значения х₁, х₂, … точек пересечения графиков функций y=f(x) и у=g(x). (Рис. 4):

Слайд 6 Найдите корни уравнения , используя графический способ решения уравнений (Рис. 5).

2. Объяснение нового материала. Практическая работа.

Решение уравнений графическим способом требует больших временных затрат на построение графиков функций и в большинстве случаев дает грубо приближенные решения. При использовании электронных таблиц, в данном случае – Microsoft Excel 2007, существенно экономится время на построение графиков функций, и появляются дополнительные возможности нахождения корней уравнения с заданной точностью (метод Подбор параметра).

I. Графический способ решения уравнений вида f(x)=0 в Excel.

Дальнейшая работа выполняется учителем в Excel одновременно с учениками с подробными (при необходимости) инструкциями и выводом результатов на проекционный экран. Слайды Приложения 1 используются для формулировки задач и подведения промежуточных итогов.

Пример1: Используя средства построения диаграмм в Excel, решить графическим способом уравнение —х 2 +5х-4=0.

Для этого: построить график функции у=-х 2 +5х-4 на промежутке [ 0; 5 ] с шагом 0,25; \найти значения х точек пересечения графика функции с осью абсцисс.

Выполнение задания можно разбить на этапы:

1 этап: Представление функции в табличной форме (рис. 6):

• в ячейку А1 ввести текст Х, в ячейку A2 — Y;
• в ячейку В1 ввести число 0, в ячейку С1 – число 0,25;
• выделить ячейки В1:С1, подвести указатель мыши к маркеру выделения, и в тот момент, когда указатель мыши примет форму черного крестика, протянуть маркер выделения вправо до ячейки V1 (Рис. 7).

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

После ввода формулы в ячейке окажется результат вычисления по формуле, а в поле ввода строки формул — сама формула (Рис. 8):

• скопировать содержимое ячейки B2 в ячейки C2:V2 за маркер выделения. Весь ряд выделенных ячеек заполнится содержимым первой ячейки. При этом ссылки на ячейки в формулах изменятся относительно смещения самой формулы.

2 этап: Построение диаграммы типа График.

• выделить диапазон ячеек B2:V2;
• на вкладке Вставка|Диаграммы|График выбрать вид График;
• на вкладке Конструктор|Выбрать данные (Рис. 9) в открывшемся окне «Выбор источника данных» щелкнуть по кнопке Изменить в поле Подписи горизонтальной оси — откроется окно «Подписи оси». Выделить в таблице диапазон ячеек B1:V1 (значения переменной х). В обоих окнах щелкнуть по кнопкам ОК;

• на вкладке Макет|Оси|Основная горизонтальная ось|Дополнительные параметры основной горизонтальной оси выбрать:

Интервал между делениями: 4;

Интервал между подписями: Единица измерения интервала: 4;

Положение оси: по делениям;

Выбрать ширину и цвет линии (Вкладки Тип линии и Цвет линии);

• самостоятельно изменить ширину и цвет линии для вертикальной оси;
• на вкладке Макет|Сетка|Вертикальные линии сетки по основной оси выбрать Основные линии сетки.

Примерный результат работы приведен на рис. 10:

3 этап: Определение корней уравнения.

График функции у=-х 2 +5х-4 пересекает ось абсцисс в двух точках и, следовательно, уравнение -х 2 +5х-4=0 имеет два корня: х₁=1; х₂=4.

II. Графический способ решения уравнений вида f(x)=g(x) в Excel.

Пример 2: Решить графическим способом уравнение .

Для этого: в одной системе координат построить графики функций у₁= и у₂=1-х на промежутке [ -1; 4 ] с шагом 0,25; найти значение х точки пересечения графиков функций.

1 этап: Представление функций в табличной форме (рис. 1):

Аналогично Примеру 1, применив приемы копирования, заполнить таблицу. При табулировании функции у₁=воспользоваться встроенной функцией Корень (Рис. 11).

2 этап: Построение диаграммы типа График.

Примерный результат работы приведен на Рис. 12:

3 этап: Определение корней уравнения.

Графики функций у₁= и у₂=1-х пересекаются в одной точке (0;1) и, следовательно, уравнение имеет один корень – абсцисса этой точки: х=0.

III. Метод Подбор параметра.

Графический способ решения уравнений красив, но далеко не всегда точки пересечения могут быть такими «хорошими», как в специально подобранных примерах 1 и 2.

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

Пример 3: Разберем метод Подбор параметра на примере решения уравнения —х 2 +5х-3=0.

1 этап: Построение диаграммы типа График для приближенного определения корней уравнения.

Построить график функции у=—х 2 +5х-3, отредактировав полученные в Примере 1 формулы.

• выполнить двойной щелчок по ячейке B2, внести необходимые изменения;
• с помощью маркера выделения скопировать формулу во все ячейки диапазона C2:V2.

Все изменения сразу отобразятся на графике.

Примерный результат работы приведен на Рис. 13:

2 этап: Определение приближенных значений корней уравнения.

График функции у=-х 2 +5х-3 пересекает ось абсцисс в двух точках и, следовательно, уравнение -х 2 +5х-4=0 имеет два корня.

По графику приближенно можно определить, что х₁≈0,7; х₂≈4,3.

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

1) Начать с поиска более точного значения меньшего корня.

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

• Выделить ячейку Е2;
• перейти на вкладку Данные|Анализ «что-если»|Подбор параметра…;

В открывшемся диалоговом окне Подбор параметра (Рис. 14) в поле Значение ввести требуемое значение функции: 0.

В поле Изменяя значение ячейки: ввести $E$1 (щелкнув по ячейке E1).

Щелкнуть по кнопке ОК.

• В окне Результат подбора (Рис. 15) выводится информация о величине подбираемого и подобранного значения функции:
• В ячейке E1 выводится подобранное значение аргумента 0,6972 с требуемой точностью (0,0001).

Установить точность можно путем установки в ячейках таблицы точности представления чисел – числа знаков после запятой (Формат ячеек|Число|Числовой).

Итак, первый корень уравнения определен с заданной точностью: х₁≈0,6972.

2) Самостоятельно найти значение большего корня с той же точностью. (х₂≈4,3029).

IV. Метод Подбор параметра для решения уравнений вида f(x)=g(x).

При использовании метода Подбор параметров для решения уравнений вида f(x)=g(x) вводят вспомогательную функцию y(x)=f(x)-g(x) и находят с требуемой точностью значения х точек пересечения графика функции y(x) с осью абсцисс.

3. Закрепление изученного материала. Самостоятельная работа.

Задание: Используя метода Подбор параметров, найти корни уравнения с точностью до 0,001.

• ввести функцию у=и построить ее график на промежутке [ -1; 4 ] с шагом 0,25 (Рис. 16):

• найти приближенное значение х точки пересечения графика функции с осью абсцисс (х≈1,4);
• найти приближенное решение уравнения с точностью до 0,001 методом Подбор параметра (х≈1,438).

4. Итог урока.

Слайд 12 Проверка результатов самостоятельной работы.

Слайд 13 Повторение графического способа решения уравнения вида f(x)=0.

Слайд 14 Повторение графического способа решения уравнения вида f(x)=g(x).

5. Домашнее задание.

Используя средства построения диаграмм в Excel и метод Подбор параметра, определите корни уравнения х 2 -5х+2=0 с точностью до 0,01.

Графическое решение уравнений средствами Microsoft Excel

Обращаем Ваше внимание, что в соответствии с Федеральным законом N 273-ФЗ «Об образовании в Российской Федерации» в организациях, осуществляющих образовательную деятельность, организовывается обучение и воспитание обучающихся с ОВЗ как совместно с другими обучающимися, так и в отдельных классах или группах.

Применение табличного процессора Microsoft Excel для графического решения уравнений n-ой степени

Из курса математики известно, что корнями уравнения являются значения точек пересечения графика функции с осью абсцисс. Если же мы решаем систему уравнений, то ее решениями будут координаты точек пересечения графиков функций. Этот метод нахождения корней уравнения называется графическим. Мы уже знаем, что с помощью EXCEL можно строить практически любые графики. Воспользуемся этими знаниями для нахождения корней системы уравнений:

Преобразуем данную систему в приведенную:

Для оценки решений воспользуемся диаграммой, на которой отобразим графики обеих функций. Сначала построим таблицу:

Первая строка – строка заголовков.

При заполнении столбца А: в ячейку А2 заносится начальное значение аргумента Х = – 10, для автоматического заполнения всего столбца в ячейку А3 занести формулу “= А2 + 1” и скопировать ее до ячейки А22.

При заполнении столбца В: в ячейку В2 заносится формула “= А2 * А2”, которая затем копируется до ячейки В22.

При заполнении столбца С: в ячейку С2 заносится формула “ = 2 * А2 + 9”, и также копируется до С22

С помощью Мастера диаграмм построим в одной координатной плоскости графики заданных функций для первоначальной оценки решений/

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

Уточним их, построив два графика в интервалах от – 3 до 0, где находится первое решение, и от 3 до 5, где находится второе решение. Составим новые таблицы. Для первого решения – рисунок 4, для второго – рисунок 5.

Для более точного построения мы уменьшили шаг изменения аргумента. Решением нашей системы будут координаты точек пересечения графиков: Х ₁ = – 2,2; Y ₁ = 4,6; Х ₂ = 4,2; Y ₂ = 17,4. Как вы уже поняли, графическое решение системы дает приблизительные результаты.
Это можно сделать, построив график и определив координаты точек его пересечения с осью OX, либо построив два графика: Y = X3;
Y = 2X2 + 4X – 12 и определив точки их пересечения.

Решение нелинейных уравнений в Excel и Mathcad (стр. 1 )

Из за большого объема этот материал размещен на нескольких страницах: 1 2 3

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ

ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

Кафедра прикладной математики и вычислительной техники

Решение нелинейных уравнений в Excel и Mathcad

к выполнению лабораторных работ

по дисциплине «Вычислительная математика»

Решение нелинейных уравнений в Excel и Mathcad: Метод. указ. / Сост. , — Самара: СГАСУ, 20с.

Методические указания разработаны в соответствии с Государственным образовательным стандартом изучения дисциплины «Вычислительная математика».

Рассмотрена реализация численных методов при решении нелинейных уравнений и систем уравнений в Excel и MathCad. Приведены варианты заданий для индивидуального выполнения и вопросы для самоконтроля и тестирования.

Предназначены для студентов специальности 230201 – «Информационные системы и технологии» всех форм обучения.

Рецензент к. ф-м. н.

Ó , составление, 2012

1 Решение нелинейного уравнения

1.1 Общие сведения о решении нелинейного уравнения

1.2 Отделение корней

1.3 Уточнение корней стандартными средствами Excel и Mathcad

1.4 Метод деления отрезка пополам

1.6 Метод Ньютона (касательных)

1.7 Комбинированный метод

1.8 Метод итераций

2 Решение систем нелинейных уравнений

2.1 Общие сведения о решении систем нелинейных уравнений

2.2 Решение систем нелинейных уравнений методом Ньютона

2.3 Решение систем нелинейных уравнений методами итераций

3 Задания к лабораторным работам

Лабораторная № 1. Отделение корней и стандартные инструменты решения нелинейного уравнения

Лабораторная № 2. Сравнение методов уточнения корней нелинейного уравнения

Лабораторная № 3. Решение систем нелинейных уравнений

Лабораторная № 4. Программирование методов решения нелинейных уравнений и систем

4 Вопросы и тесты для самоконтроля

Список рекомендуемой литературы

1 Решение нелинейного уравнения

1.1 Общие сведения о решении нелинейного уравнения

Как правило, нелинейное уравнения общего вида f(х)=0 невозможно решить аналитически. Для практических задач достаточно найти приближенное значение x, в определенном смысле близкое к точному решению уравнения хточн.

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

Достигнутая точность может оцениваться либо «по функции» (в найденной точке x, функция достаточно близка к 0, т. е. выполняется условие |f(x)|≤ ef, где ef требуемая точность по оси ординат), либо «по аргументу» (найден достаточно маленький отрезок [a,b], внутри которого находится корень, т. е. |b–a|≤ ex, где ex требуемая точность по оси абсцисс).

1.2 Отделение корней

Отделение корней может производиться сочетанием графического и аналитического исследования функции. Такое исследование опирается на теорему Вейерштрасса, в соответствии с которой для непрерывной на отрезке [a,b] функции f(х) и любого числа y, отвечающего условию f(a)≤y≤f(b), существует на этом отрезке точка x, в которой функция равна y. Следовательно, для непрерывной функции достаточно найти отрезок, на концах которого функция имеет разные знаки, и можно быть уверенным, что на этом отрезке есть корень уравнения f(х)=0.

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

Пример Найти с точностью до целых все корни нелинейного уравнения y(x)=x3 ‑ 10x + 7=0 а) построив таблицу и б) построив график. Найти корень уравнения на выделенном отрезке, используя опции «Подбор параметра» и «Поиск решения».

Решение Создадим в Excel таблицу, содержащую аргументы и значения функции и по ней построим точечную диаграмму. На рисунке 1 приведен снимок решения.

На графике видно, что уравнение имеет три корня, принадлежащие отрезкам [-4, -3], [0, 1] и [2, 3]. Эти отрезки можно выявить и наблюдая за сменой знаков функции в таблице. По построенному графику можно сделать вывод, что на указанных отрезках функция f(x) монотонна и, следовательно, на каждом из них содержится только по одному корню.

Такой же анализ может быть выполнен и в пакете Mathcad. Для этого достаточно набрать определение функции f(x), используя оператор присваивания (:=) и естественные общепринятые обозначения математических операций и стандартных функций, задать цикл для изменения аргумента, например, а затем вывести на экран таблицу значений функции (располо­жен­ными в одной строке командами x= f(x)=) и график. Цикл можно задать, например, командой x:=-5,-4.5…5. Шаг цикла формируется путем задания начального и следующего за ним значений переменной, а перед конечным значением переменной ставится точка с запятой, которая будет визуально отображена на экране в виде многоточия.

Рисунок 1 – Таблица и график для отделения корней нелинейного уравнения

1.3 Уточнение корней стандартными средствами Excel и Mathcad

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

В Excel для уточнения значений корней можно использовать опции «Подбор параметра» и «Поиск решения». Пример оформления решения приведен на рисунках 2 и 3.

Рисунок 2 – Ввод значений для использования средств решения уравнения в Excel

Рисунок 3 – Результаты использования средств решения уравнения в Excel

В Mathcad для уточнения корней уравнения можно использовать функцию root(….) или блок решения. Пример использования функции root(…) приведен на рисунке 4, а блока решения на рисунке 5. Следует обратить внимание, что в блоке решения (после заголовка блока Given) между левой и правой частями уравнения должен стоять жирный знак равенства (тождества), который можно получить выбором из соответствующей палитры инструментов, либо нажатием одновременно клавиши Ctrl и =.

Рисунок 4 – Решение уравнения с использованием функции root(…) в Mathcad

Рисунок 5 – Решение уравнения с использованием блока решения в Mathcad

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

В то же время, очень просто построить собственную таблицу или написать программу, реализующие один из методов уточнения корней. Здесь можно использовать критерии точности расчета, задаваемые пользователем. При этом достигается и понимание процесса расчетов без опоры на принцип Митрофанушки: «Извозчик есть, довезет».

Далее рассмотрены несколько наиболее распространенных методов. Отметим очевидный момент: при прочих равных условиях тот метод уточнения корней будет более эффективен, в котором результат с той же погрешностью найден с меньшим числом вычислений функции f(x) (при этом достигается и максимальная точность при одинаковом числе вычислений функции).

1.4 Метод деления отрезка пополам

В этом методе на каждом шаге отрезок делится на две равные части. Затем сравнивают знаки функции на концах каждой из двух половинок (например, по знаку произведения значений функций на концах), определяют ту из них, в которой содержится решение (знаки функции на концах должны быть разные), и. сужают отрезок, перенося в найденную точку его границу (а или b). Условием окончания служит малость отрезка, где содержится корень («точность по x»), либо близость к 0 значения функции в средине отрезка («точность по y»). Решением уравнения считают середину отрезка, найденного на последнем шаге.

Пример. Построить таблицу для уточнения корня уравнения x3 –10x+7=0 на отрезке [-4, -3] методом деления отрезка пополам. Определить сколько шагов надо сделать методом деления отрезка пополам и какая при этом достигается точность по х, для достижения точности по y, равной 0,1; 0,01; 0, 001.

Решение Для решения можно использовать табличный процессор Excel, позволяющий автоматически продолжать строки. На первом шаге заносим в таблицу значения левого и правого концов выбранного начального отрезка и вычисляем значение середины отрезка с=(a+b)/2, а затем вводим формулу для вычисления функции в точке a (f(a)) и растягиваем (копируем) её для вычисления f(c) и f(b). В последнем столбца вычисляем выражение (b—a)/2, характеризующего степень точности вычислений. Все набранные формулы можно скопировать во вторую строку таблицы.

На втором шаге нужно автоматизировать процесс поиска той половины отрезка, где содержится корень. Для этого испльзуется логическая функция ЕСЛИ (Меню: ВставкаФункцияЛогические). Для нового левого края отрезка мы проверяем истинность условия f(a)*f(c)>0, если оно верно, то мы в качестве нового значения левого конца отрезка берем число c (т. к. это условие показывает, что корня на отрезке [a, c] нет), иначе оставляем значение a. Аналогично, для нового правого края отрезка мы проверяем истинность условия f(c)*f(b)>0, если оно верно, то мы в качестве нового значения правого конца отрезка берем число c (т. к. это условие показывает, что корня на отрезке [c, b] нет), иначе оставляем значение b.

Вторую строку таблицы можно продолжить (скопировать) на необходимое число последующих строк.

Итерационный процесс завершается, когда очередное значение в последнем столбце становится меньшим, чем заданный показатель точности ex. При этом, значение середины отрезка в последнем приближении, принимается в качестве приближенного значения искомого корня нелинейного уравнения. На рисунке 6 приведен снимок решения. Для построения аналогичного процесса в Mathcad можно использовать бланк, подобный приведенному на рисунке 7. Число шагов N может варьиро­вать­ся до достижения в таблице результатов требуемой точности. При этом таблица будет автоматически удлиняться или укорачиваться.

Итак, одним из трех корней нелинейного уравнения x3 – 10x + 7=0, найденным с точностью e=0,0001, является x= — 3,46686. Как мы видим, он действительно принадлежит отрезку [-4; -3].

Рисунок 6 – Уточнение корня методом деления отрезка пополам в Excel

Рисунок 7 – Уточнение корня методом деления отрезка пополам в Mathcad

1.5 Метод хорд

В этом методе нелинейная функция f(x) на отделенном интервале [а, b] заменяется линейной – уравнением хорды, т. е. прямой соединяющей граничные точки графика на отрезке. Условие применимости метода – монотонность функции на начальном отрезке, обеспечивающая единственность корня на этом отрезке. Расчет по методу хорд аналогичен расчету методом деления отрезка пополам, но теперь на каждом шаге новая точка x внутри отрезка [a,b] рассчитывается по любой из следующих формул:

.

1.6 Метод Ньютона (касательных)

Идея, на которой основан метод, аналогична той, которая реализована в методе хорд, только на каждом шаге кривая f(x) заменяется касательной к ней, проведенной в предыдущей найденной точке. В качестве начальной точки в зависимости от свойств функции берется или левая граница отрезка, содержащего корень – x0 = а (если f(а) f»(х) > 0), или правая его граница: x0 = b (если f(b) f»(х)>0). Расчет нового приближения на следующем шаге i+1 производится по формуле:

.

Алгоритм применим для монотонных функций, сохраняющих выпуклость или вогнутость в промежутке между начальным приближением и корнем уравнения (т. е. должен сохраняться знак первой и второй производных функции f(x)). работоспособен при выпуклых и монотонных функциях f(x). В расчетах нет необходимости отслеживать две границы отрезка, поэтому достаточно на каждом шаге вычислять значения x, f(x) и f′(x). При этом легко оценить «точность по y», по значению левой части уравнения на очередном шаге. Для оценки «точности по x» нужно отслеживать разницу приближений на предыдущем и последующих шагах, которая связана с разницей между найденным приближением и точным значением корня.

Следует обратить внимание на следующую особенность метода: последовательность x1, x2, x3,… приближается к корню с другой стороны, в отличие от использования метода хорд при прочих равных условиях.

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

Уточнить корень уравнения tg (0,55x+0,1) – x2=0 на отрезке [0.6, 0.8] методом касательных до точности 0,001.

Точность вычислений можно оценить из соотношения

2 Решение систем нелинейных уравнений

2.1 Общие сведения о решении систем нелинейных уравнений

Систему n нелинейных уравнений с n неизвестными x1, x2, . xn записывают в виде:

где F1, F2,…, Fn – функции независимых переменных, среди которых есть нелинейные.

Как и в случае систем линейных уравнений, решением системы является такой вектор X*, который при подстановке обращает одновременно все уравнения системы в тождества.

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

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

Рассмотрим две группы таких методов: метод Ньютона с различными его модификациями и методы итераций (простых итераций и Зейделя).

2.2 Решение систем нелинейных уравнений методом Ньютона

Будем рассматривать этот метод на примере системы двух нелинейных уравнений с двумя неизвестными:

Начальные значения x0 и y0 определяются графически. Для нахождения каждого последующего приближения (xi+1, yi+1) используют вектор значений функций и матрицу значений их первых производных, рассчитанные в предыдущей точке (xi, yi).

,

Для расчета новых приближений на шаге i+1 используется матричная формула

.

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

Расчет останавливают при выполнении одного (а иногда и обоих) из двух условий. Первое из них заключается в том, что на очередном шаге максимальное по модулю из изменений аргументов x и y становится меньше заданная погрешность по аргументам. В соответствии со вторым из условий, на очередном шаге максимальное по модулю значение левых частей уравнений должно отличаться от нуля меньше, чем заданная погрешность по функциям.

В упрощенном методе Ньютона матрица производных и матрица, обратная к ней вычисляются только один раз (в начальной точке) и для расчетов используется матричная формула

.

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

2.3 Решение систем нелинейных уравнений методами итераций

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

.

Для решения такой системы задаются начальным приближением x0, y0. Уточненные решения получают по шагам, подставляя в правые части уравнений значения, найденные на предыдущем шаге. В методе простых итераций для уточнения решения используют формулы:

.

Если одно из решений системы и начальные значения x0 и y0 лежат в области D, задаваемой неравенствами: a ≤ x ≤ b, c ≤ y ≤ d, то расчет по методу простых итераций сходится при выполнении в области D соотношений: