Решить нелинейное уравнение в excel методом ньютона

Метод Ньютона в Excel

Как видно, процесс нахождения корней нелинейного уравнения методом Ньютона состоит из следующих этапов:

1. Получения шаблона.
2. Уточнение интервалов в ячейках B2 , B3 .
3. Замена в формуле ЕСЛИ запятую ( , ) на точку с запятой ( ; ).
4. Копирование строки итераций до требуемой точности (столбец E ).

Примечание: столбец A — номер итерации, столбец B — корень уравнения X , столбец C — значение функции F(X) , столбец D — значение первой производной dF(X) , столбец E — точность eps .

Трансцендентные уравнения? «Подбор параметра» в Excel!

Нелинейные, трансцендентные уравнения функции одной переменной – это уравнения вида f (x) = 0, в которых нельзя найти алгебраическими методами корни. Функция f (x) – это, как правило, достаточно сложная и громоздкая функция, содержащая в своем составе.

. тригонометрические, логарифмические, степенные и иные нелинейные функции с различной глубиной вложенности. Например: f (x) = sin (3,14^x) + cos (x) = 0. Уравнения такого вида решаются численными методами.

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

Чуть-чуть истории и теории.

Вы задумывались когда-нибудь — откуда и зачем в головах людей, живших в XVI…XVII веках, родились понятия дифференциалов, производных, интегралов? Объяснение, в общем-то, достаточно простое и понятное – эти ученые искали аналитические пути решения прикладных практических задач. И успешно находили.

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

1. Арифметика — сложение, вычитание, умножение, деление.

2. Алгебра – применение элементарных функций (степенной, логарифмической, тригонометрической, …) и алгебраических уравнений функции одной переменной.

3. Гауссовские системы линейных уравнений.

4. Численные методы решения трансцендентных уравнений.

5. Численные методы решения систем трансцендентных уравнений функций нескольких переменных.

6. Дифференцирование и интегрирование функций одной переменной.

7. Дифференцирование и интегрирование функций нескольких переменных.

8. Системы дифференциальных и интегральных уравнений.

9. Масса разнообразных новых и старых специальных методик и подходов мне не известных и известных, но, безусловно, существующих и работающих.

Предлагаю остановиться и разобраться с достаточно высокой четвертой ступенью «лестницы».

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

Каждый из методов имеет свои достоинства и недостатки — они подробно описаны в литературе, и углубляться в них мы не будем. Скажу только, что из вышеперечисленных методов мне на практике довелось использовать все. При решении различных (в основном геометрических и теплотехнических) задач по разным причинам было удобно использовать то один, то другой подход. Метод Ньютона хорош своей быстрой сходимостью и простотой формулы. Комбинированный метод секущих-хорд на основе итерационной формулы Ньютона не требует нахождения производных, быстро «сходится», и главное – не требует анализа функции на сходимость. Метод половинного деления медленно сходится, но не требует никакого предварительного анализа функции.

Трансцендентные уравнения. Два метода решения в Excel.

Если у вас на компьютере нет программы MS Excel, то расчеты можно выполнить в программе OOo Calc из бесплатного пакета Open Office.

Задач, которые требуют для получения ответа составления и решения трансцендентных уравнений, вокруг нас очень много. Это — задачи и физики, и теплотехники, и астрономии, и элементарной геометрии в обычной жизни… Инженерам-конструкторам и программистам в повседневной работе необходимо уметь составлять и быстро решать численными методами нелинейные уравнения. На мой взгляд — это один из критериев профессионализма. Более того, уравнения, которые решаются аналитически, сегодня иногда гораздо проще и быстрее при наличии вычислительной техники решить численными методами, поэтому нужно уметь это делать.

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

Рассмотрим пример из статьи «Расчет геометрии зубчатой передачи». Необходимо найти угол зацепления зубчатой передачи atw . Я обещал в той статье рассказать, как это делается. Выполняю обещание.

Если расстояние между центрами колеса и шестерни не задано, то угол зацепления можно вычислить путем решения трансцендентного уравнения:

inv ( atw )=tg ( atw ) — atw =2* xs *tg ( a )/( z2 + T * z1 )+ tg ( at ) — at

Подставив данные из примера, рассмотренного в вышеупомянутой статье, получим после преобразований следующее уравнение:

inv ( atw )=0,020910

f ( atw )=tg ( atw )— atw -0,020910=0

Используем метод Ньютона, потому что взять производную представленной выше функции элементарно просто, а итерационная формула очень проста и компактна:

f’( atw )=1/(cos ( atw ))^2—1

atw (i+1) = atw i — f ( atw ) i/ f’( atw ) i

Открываем файл Excel и начинаем работу.

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

1. Инволюту угла зацепления inv( atw ) заносим

в ячейку D3: 0,020910

2. Значение угла зацепления в нулевом приближении atw 0 в радианах записываем

3. Итерационную формулу atw (i+1)= atw i— f( atw )i/ f’( atw )i заносим

в D5: =D4- (TAN (D4) -D4-$D$3)/(1/(COS (D4))^2-1) =0,591706

atw 1= atw 0- (tg ( atw 0) — atw 0- inv ( atw ))/(1/(cos ( atw 0))^2-1)

и копируем в ячейки D6… D14

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

atw =D13- (TAN (D13) -D13-$D$3)/(1/(COS (D13))^2-1) =0,389140

Решение найдено, расчет в Excel завершен!

Решение задачи ландшафтного дизайна с помощью сервиса «Подбор параметра» в Excel

Задача:

Вдоль отмостки стены дома длиной 14 метров необходимо разбить цветник в виде сегмента круга площадью ровно 16 квадратных метров. На сколько метров цветник будет отстоять от края отмостки по центру стены? Каким радиусом необходимо выполнить границу цветника?

1. Длину отмостки стены дома — хорды сегмента круга x в метрах записываем

в ячейку D17: 14,000

2. Площадь цветника – сегмента круга S в квадратных метрах вписываем

в D18: 16,000

3. Предположительное произвольное (не нулевое) значение центрального угла сегмента a в радианах пишем

Трансцендентное уравнение a / sin( a /2 ) -2*cos ( a /2) — (8* S / x ^2) *sin( a /2)=0 вводим

в объединенную ячейку E19F19: =D19/SIN (D19/2) -2*COS (D19/2) — (8*D18/D17^2)*SIN (D19/2)

Включаем сервис «Подбор параметра» в Excel: «Сервис» – «Подбор параметра». Пишем в появившемся окне все как на рисунке слева и нажимаем кнопку OK.

В появившемся новом окне видим, что решение найдено, снова нажимаем на кнопку OK.

Считываем искомое значение центрального угла сегмента a в радианах

в D19: 0,950057

При этом видим, что значение трансцендентного уравнения равно нулю; считываем

в объединенной ячейке E19F19: =D19/SIN (D19/2) -2*COS (D19/2) — (8*D18/D17^2)*SIN (D19/2) =0

4. Радиус наружной границы цветника – радиус сегмента круга r в метрах рассчитывается

в D20: =D17/2/SIN (D19/2) =15,305

r = x /2/sin( a /2)

5. Максимальная ширина цветника – высота сегмента круга h в метрах рассчитывается

в ячейке D21: =D20*(1-COS (D19/2)) =1.695

h = r *(1- cos( a /2))

Ответы получены, вторая задача успешно решена!

Я не приводил вывода использованных формул потому, что это не по теме поста, и, думаю, с геометрией и тригонометрией вы легко разберетесь. Будут вопросы – обращайтесь.

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

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

1. Итерационными численными методами удобно и быстро можно решать трансцендентные уравнения и громоздкие нелинейные алгебраические.

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

3. Использование инструмента «Подбор параметра» в Excel является сегодня, безусловно, наиболее оптимальным и эффективным методом решения нелинейных, трансцендентных уравнений функций одной переменной, а также проведения анализа типа «Что будет? Если…».

Умение применять в работе сервис «Подбор параметра» существенно повышает ваш уровень, как специалиста вообще, так и как пользователя Excel – в частности.

Буду очень рад увидеть ваши комментарии к статье, уважаемые читатели!

1. Решение нелинейных уравнений в MS Excel

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

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

Например , для уравнения выполним преобразование и приведем его к виду f ( x )= 0 т.е. . График этой функции представлен на рисунке 1. Очевидно, что данное уравнение имеет два действительных корня – один на отрезке [-1, 0] , а второй – [1, 2].

Рисунок 1. График функции

1.2 Решение уравнений, используя инструмент “Подбор параметра”

Используя возможности Excel , можно находить корни нелинейного уравнения вида f ( x )=0 в допустимой области определения переменной. Последовательность операций нахождения корней следующая:

1. Производится вычисление значений функции в диапазоне вероятного существования корней от значений аргумента, изменяющегося с определенным шагом;

2. В таблице выделяются ближайшие приближения к значениям корней (пары соседних значений функции с разными знаками);

3. Используя средство Excel Подбор параметра, вычисляются корни уравнения.

2. Работа с матрицами в MS Excel . Решение систем уравнений.

Нахождение определителя матрицы

Перед нахождением определителя необходимо ввести матрицу в диапазон ячеек Excel в виде таблицы.

Для нахождения определителя матрицы в Excel необходимо:

· сделать активной ячейку, в которой в последующем будет записан результат;

· в меню Вставка – Функция в категории Математические выбрать функцию МОПРЕД и нажать OK ;

· на втором шаге задать диапазон ячеек, в котором содержатся элементы матрицы, и нажать OK .

Нахождение обратной матрицы

Для нахождения обратной матрицы необходимо

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

· в меню Вставка – Функция в категории Математические выбрать функцию МОБР и нажать OK ;

· на втором шаге задать диапазон ячеек, в котором содержатся элементы исходной матрицы, и нажать OK .

· после появления значения в левом верхнем углу выделенного диапазона последовательно нажать клавишу F 2 и комбинацию клавиш Ctrl + Shift + Enter .

Для перемножения матриц необходимо

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

· в меню Вставка – Функция в категории Математические выбрать функцию МУМНОЖ и нажать OK ;

· на втором шаге задать два диапазона ячеек с элементами перемножаемых матриц, и нажать OK .

· после появления значения в левом верхнем углу выделенного диапазона последовательно нажать клавишу F 2 и комбинацию клавиш Ctrl + Shift + Enter .

Решение системы уравнений в Excel .

Решение системы уравнений при помощи нахождения обратной матрицы.

Пусть дана линейная система уравнений.

Данную систему уравнений можно представить в матричной форме:

Матрица неизвестных вычисляется по формуле

где A -1 – обратная матрица по отношению к A .

Для вычисления уравнения в Excel необходимо:

· ввести матрицу A;

· ввести матрицу B;

· вычислить обратную матрицу по отношению к А ;

· перемножить полученную обратную матрицу с матрицей B .

Порядок выполнения работы

Задание 1

Найти все корни уравнения 2x 3 -15sin( x )+0,5x-5=0 на отрезке [-3 ; 3].

1. Построить таблицу значений функции f ( x ) для значений x от –3 до 3, шаг 0,2.

Для этого ввести первые два значения переменной x , выделить эти две ячейки, с помощью маркера автозаполнения размножить значения до 3.

Затем ввести формулу для вычисления f ( x ). Скопировать формулу с использованием маркера автозаполнения на весь столбец.

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

2. Выделить цветом пары значений x и f ( x ), где f ( x ) меняет знак (см .р исунок 2).

3. Построить график функции f ( x ).

Рисунок 2. Поиск приближенных значений корней уравнения

4. Скопировать рядом с таблицей произвольную пару выделенных значений x и f ( x ) (см .р исунок 3).

5. Выполнить команду меню Сервис/Подбор параметра. В диалоговом окне (рисунок 3) заполнить следующие поля:

þ Установить в ячейке : в поле указывается адрес ячейки, в которой записана формула правой части функции;

þ Значение : в поле указывается значение, которое должен получить полином в результате вычислений, т.е. правая часть уравнения (в нашем случае 0);

þ Изменяя значение : в поле указывается адрес ячейки (где записано начальное приближение), в которой будет вычисляться корень уравнения и на которую ссылается формула.

Рисунок 3. Диалоговое окно Подбор параметра для поиска первого корня

6. После щелчка на ОК должно получиться значение первого корня -1,65793685 .

7. Выполнить последовательно операции, аналогичные предыдущим, для вычисления значений остальных корней: -0,35913476 и 2,05170101 .

Задание 2

Решить систему уравнений:

1. Ввести значения элементов матриц A и B уравнения в ячейки Excel .

2. Вычислить обратную матрицу с помощью матричной функции МОБР.

3. Перемножить обратную матрицу A -1 на матрицу B с помощью матричной функции МУМНОЖ (Порядок умножения важен ­– первой должна идти матрица A -1 а второй B .)

4. Проверить правильность полученной матрицы корней X .

Контрольные вопросы

1. Порядок действий для решения нелинейного уравнения с помощью инструмента Подбор параметра MS Excel .

2. Порядок действий для решения системы уравнений матричным методом в MS Excel .