Метод простых итераций для нелинейного уравнения в маткаде

решение уравнений маткад. Нелин_уравнения. Практическая работа 5 Численное решение нелинейных уравнений в Mathcad (метод простых итераций)
НазваниеПрактическая работа 5 Численное решение нелинейных уравнений в Mathcad (метод простых итераций)
Анкоррешение уравнений маткад
Дата24.02.2021
Размер318.5 Kb.
Формат файлаМетод простых итераций для нелинейного уравнения в маткаде
Имя файлаНелин_уравнения.doc
ТипПрактическая работа
#179215
Подборка по базе: Итоговая контрольная работа по обществознанию в 10 классе.docx, Практическая работа №3 «Анализ содержания и методического аппара, Контрольная работа.docx, Лабораторная работа 8.doc, Практическая работа №5.docx, Контрольная работа по русскому языку за 1 полугодие.docx, Практическая работа № 5.docx, Практическая работа №5 (1).docx, Практическая работа № 4..docx, Практическая работа №1 Изучение содержаия обновленных ФГОС ООО.

Видео:Решение нелинейного уравнения методом простых итераций (программа)Скачать

Решение нелинейного уравнения методом простых итераций (программа)

Практическая работа №5

Цель работы:

  • научиться решать нелинейное уравнение в Mathcad, используя вычислительный алгоритм метода простых итераций;

Задание №1.Решить нелинейное уравнение методом итераций.

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

где f(x)  непрерывная функция, требуется определить вещественные корни этого уравнения. Запишем уравнение (1) равносильным уравнением

x= Метод простых итераций для нелинейного уравнения в маткаде(x) (2)

Выберем каким-либо способом, например графически, приближенное значение корня x0 и подставим его в правую часть уравнения (2). Получим некоторое число

x1= Метод простых итераций для нелинейного уравнения в маткаде(x0) (3)

Подставим теперь в правую часть уравнения (3) вместо x0 число x1. Получим новое число x2= Метод простых итераций для нелинейного уравнения в маткаде(x1). Повторяя этот процесс, получим последовательность чисел xn= Метод простых итераций для нелинейного уравнения в маткаде(xn-1) (n=1,2,…). Если эта последовательность сходится, то найдется число Метод простых итераций для нелинейного уравнения в маткаде, которое и будет корнем уравнения (2). Итерационный процесс xn= Метод простых итераций для нелинейного уравнения в маткаде(xn-1) сходится, если Метод простых итераций для нелинейного уравнения в маткаде. Критерием окончания итерационного процесса является выполнение неравенства

Метод простых итераций для нелинейного уравнения в маткаде, где ε 3-x-1=0

  • Метод простых итераций для нелинейного уравнения в маткаде
  • x5-x-0.2=0
  • Видео:Метод простых итераций пример решения нелинейных уравненийСкачать

    Метод простых итераций пример решения нелинейных уравнений

    Практическая работа №6

    Цель работы:

    • научиться решать нелинейное уравнение в Mathcad, используя вычислительный алгоритм метода бисекции (половинного деления).

    Задание №1.Решить нелинейное уравнение методом бисекции.

    Метод бисекции (половинного деления) состоит в следующем.

    Пусть дано уравнение

    Деление продолжать до получения отрезка длиной ε ≤ 10-3.

    • Определите левую часть нелинейного уравнения как функцию f(x).
    • Построим график функции f(x), точка пересечения которого с осью и даст приближенное значение корня исходного уравнения. Определим по графику значения концов отрезка, содержащего корень уравнения: a=0, b=1.
    • Задайте значения концов отрезка [a,b] в векторной форме. Для этого, используя палитру векторов и матриц, получите, оператор вектора на два элемента и заполните его идентификаторами концов отрезка a0 и b0. При помощи клавиши пробела выведите курсор за знак оператора, нажмите клавишу :, введите еще один оператор вектора на два элемента и заполните его значениями концов отрезка: 0 и 1.
    • Задайте дискретную переменную i (от 0 до 10), определяющую количество этапов деления отрезка [a,b], за которое мы рассчитываем получить решение уравнения.
    • Определите формулы вычисления концов отрезка на каждом этапе деления его пополам. На отрезке [ai,bi] серединой будет точка Метод простых итераций для нелинейного уравнения в маткаде.

    Если Метод простых итераций для нелинейного уравнения в маткаде(знаки функции на концах разделенного отрезка разные), то начало нового отрезка для вычислений ai+1 совпадет с ai, т.е. ai+1=ai. В противном случае (знаки функции на концах разделенного отрезка одинаковые) ai+1 совпадет с точкой половинного деления отрезка, т.е. Метод простых итераций для нелинейного уравнения в маткаде.

    Подобная операция необходима для другого конца отрезка [ai, bi], для bi. Если Метод простых итераций для нелинейного уравнения в маткаде, то конец нового отрезка для вычислений bi+1 совпадет с точкой половинного деления отрезка, т.е. Метод простых итераций для нелинейного уравнения в маткаде; в противном случае bi+1 совпадет с концом отрезка [ai, bi], т.е. bi+1=bi. Если же Метод простых итераций для нелинейного уравнения в маткаде, то точка Метод простых итераций для нелинейного уравнения в маткаде будет являться корнем нелинейного уравнения. Для реализации описанного процесса в Mathcad запишем формулы для определения начала и конца разделенного отрезка, воспользовавшись условным оператором if в векторной форме. Такая форма записи позволяет в вычислениях Mathcad использовать результаты предыдущих операций. Введем оператор вектора на три элемента и заполним его идентификаторами искомых величин ai+1, bi+1, gi (gi  для вычисления Метод простых итераций для нелинейного уравнения в маткаде). Введем знак присвоения и еще один оператор вектора на три элемента. Последний оператор заполните следующими формулами:

    Метод простых итераций для нелинейного уравнения в маткаде

    Условный оператор if работает следующим образом. Если логическое выражение, записанное в скобках, истинно, то переменная слева от знака := принимает значение выражения или величины, стоящей первой после логического выражения, в противном случае переменной присваивается значение выражения, или величины, стоящей на втором месте после логического выражения.

    • Чтобы увидеть результат вычислений, введите последовательно a=,b=,g=. Вектор a покажет начала всех разделенных отрезков, b  концы этих отрезков. Вектор g содержит значения функции f(x) в серединах всех отрезков [a,b]. Среди значений вектора a или вектора b будет корень исходного уравнения.
    • Определим корень уравнения следующим образом. Вычислим длину отрезка [ai,bi]i=0,…,10, определив ее как Метод простых итераций для нелинейного уравнения в маткаде. Проанализируем значения вектора ε. Найдем элемент вектора, значение которого меньше заданной в условии задачи минимальной длины отрезка [ai,bi]. В нашем примере это значение элемента ε10. Корнем уравнения будет значение Метод простых итераций для нелинейного уравнения в маткаде.

    Результат на экране:

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

    Локализация корня уравнения.

    Метод простых итераций для нелинейного уравнения в маткаде

    Д Метод простых итераций для нелинейного уравнения в маткаде
    Метод простых итераций для нелинейного уравнения в маткадеискретная переменная i.

    Формулы для определения начала и конца отрезка и определения значения f(x) в середине отрезка.
    Метод простых итераций для нелинейного уравнения в маткаде
    Метод простых итераций для нелинейного уравнения в маткадеЗначения длины вычисляемых отрезков.
    Р Метод простых итераций для нелинейного уравнения в маткаде
    езультаты вычислений.
    Метод простых итераций для нелинейного уравнения в маткадеДлина отрезка [a10, b10] -3 .
    З Метод простых итераций для нелинейного уравнения в маткаденачение f(x) наиболее близкое к нулю.

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

    Методом бисекции (половинного деления) уточнить корень уравнения.

    1. x3+x-1000=0 с точностью до 10-4 .
    2. tg(x)=x на отрезке Метод простых итераций для нелинейного уравнения в маткадес точностью до 0.0001.

    Видео:1 3 Решение нелинейных уравнений методом простых итерацийСкачать

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

    Практическая работа №7

    Цель работы:

    • научиться решать нелинейное уравнение в Mathcad, используя вычислительный алгоритм метода Ньютона (метод касательных).

    Задание №1.Решить нелинейное уравнение методом Ньютона.

    Знаменитый метод Ньютона является одним из наиболее эффективных способов численного решения нелинейных уравнений. Этот метод состоит в следующем. Пусть дано уравнение

    где f(x)  непрерывная функция, требуется определить вещественные корни этого уравнения с точностью ε. Расчетная формула метода Ньютона имеет следующий вид

    Метод простых итераций для нелинейного уравнения в маткаде(2)

    Выберем каким-либо способом, например, графически, приближенное значение корня x0 и, подставляя его в правую часть уравнения (2), можем начать итерационный процесс вычисления корня уравнения.

    Условием завершения итерационного процесса является выполнение условия

    Метод простых итераций для нелинейного уравнения в маткаде(3)

    В случае выполнения неравенства (3), корнем уравнения (1) будем считать значение xk+1.
    Рассмотрим реализацию алгоритма метода Ньютона в Mathcad на следующем примере.

    Видео:4 Метод простой итерации Mathcad Решение системы линейных уравнений СЛАУСкачать

    4 Метод простой итерации Mathcad Решение системы линейных уравнений СЛАУ

    Найти действительные корни уравнения

    Метод простых итераций для нелинейного уравнения в маткаде

    с точностью до трех значащих цифр.

    • Установим формат чисел так, чтобы результаты вычислений отображались с нужным количеством знаков после десятичной точки (запятой).
    • Преобразуем заданное уравнение к виду (1) и определим его правую часть как функцию f(x).
    • Определим приближенное значение корня уравнения графически. Для этого построим график функции Метод простых итераций для нелинейного уравнения в маткаде. Этим значением будет значение x=a в точке пересечения графика функции с осью абсцисс.
    • Определим начальное значение для вычислений следующим образом:

    x0:a.

    • Определим дискретный аргумент k:=0..9, задающий количество итераций (вычисляемых значений xk).
    • Определим первую производную функции f(x) как функцию d(x):

    Метод простых итераций для нелинейного уравнения в маткаде

    • Определим вычислительную формулу Ньютона в векторной форме:

    Метод простых итераций для нелинейного уравнения в маткаде

    • Для получения вектора результатов в свободном месте рабочего пространства Mathcad напишите: x=. Введите условие оценки погрешности вычислений Метод простых итераций для нелинейного уравнения в маткаде. Результат вычислений погрешности выведите на экран. Корнем уравнения будет x2=1.171.

    Результат на экране:

    Метод простых итераций для нелинейного уравнения в маткаде

    Нелинейное уравнение определено в виде функции.

    Локализация корней уравнения графическим способом.

    Метод простых итераций для нелинейного уравнения в маткаде

    Приближенное значения корня уравнения.

    Метод простых итераций для нелинейного уравнения в маткаде

    К Метод простых итераций для нелинейного уравнения в маткадеоличество итераций.

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

    Метод простых итераций для нелинейного уравнения в маткаде

    Вычислительная формула Ньютона.

    Результат вычислений. Расчет погрешности вычислений.

    Метод простых итераций для нелинейного уравнения в маткаде Метод простых итераций для нелинейного уравнения в маткаде

    Метод простых итераций для нелинейного уравнения в маткаде

    Погрешность вычислений становится допустимой на второй итерации Метод простых итераций для нелинейного уравнения в маткаде, поэтому допустимое решение получено на второй итерации Метод простых итераций для нелинейного уравнения в маткаде.
    Задачи для самостоятельного решения.

    Используя метод Ньютона, найти с точностью ε=10-6 положительный корень уравнений

    1. Метод простых итераций для нелинейного уравнения в маткаде
    2. Метод простых итераций для нелинейного уравнения в маткаде

    Видео:3 Шаговый метод Mathcad Численные методы решения нелинейного уравненияСкачать

    3 Шаговый метод Mathcad Численные методы решения нелинейного уравнения

    Решение нелинейных уравнений в 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) (при этом достигается и максимальная точность при одинаковом числе вычислений функции).

    Видео:7 Метод половинного деления Mathcad Численные методы решения нелинейного уравненияСкачать

    7 Метод половинного деления Mathcad Численные методы решения нелинейного уравнения

    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). В последнем столбца вычисляем выражение (ba)/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

    Видео:11 Метод Ньютона (Метод касательных) Mathcad Численные методы решения нелинейного уравненияСкачать

    11 Метод Ньютона (Метод касательных) Mathcad Численные методы решения нелинейного уравнения

    1.5 Метод хорд

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

    Метод простых итераций для нелинейного уравнения в маткаде.

    Видео:Mathcad Prime. Урок 5 - Способы решения уравненийСкачать

    Mathcad Prime. Урок 5 - Способы решения уравнений

    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.

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

    Метод простых итераций для нелинейного уравнения в маткаде

    Видео:10 Численные методы решения нелинейных уравненийСкачать

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

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

    Видео:Метод Ньютона (метод касательных) Пример РешенияСкачать

    Метод Ньютона (метод касательных) Пример Решения

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

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

    Метод простых итераций для нелинейного уравнения в маткаде

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

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

    Метод простых итераций для нелинейного уравнения в маткаде

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

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

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

    Видео:2.2 Итерационные методы решения СЛАУ (Якоби, Зейделя, релаксации)Скачать

    2.2 Итерационные методы решения СЛАУ (Якоби, Зейделя, релаксации)

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

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

    Метод простых итераций для нелинейного уравнения в маткаде

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

    Метод простых итераций для нелинейного уравнения в маткаде,

    Метод простых итераций для нелинейного уравнения в маткаде

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

    Метод простых итераций для нелинейного уравнения в маткаде.

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

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

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

    Метод простых итераций для нелинейного уравнения в маткаде.

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

    Видео:Mathcad-09. Пример: уравненияСкачать

    Mathcad-09. Пример: уравнения

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

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

    Метод простых итераций для нелинейного уравнения в маткаде.

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

    Метод простых итераций для нелинейного уравнения в маткаде.

    Если одно из решений системы и начальные значения x0 и y0 лежат в области D, задаваемой неравенствами: axb, cyd, то расчет по методу простых итераций сходится при выполнении в области D соотношений:

    Видео:3.Системы нелинейных уравнений MathcadСкачать

    3.Системы нелинейных уравнений Mathcad

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

    Дата добавления: 2015-07-23 ; просмотров: 4914 ; Нарушение авторских прав

    Ознакомиться с основными методами решения нелинейных уравнений и их реализацией в пакете MathCAD.

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

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

    Широко распространен графический способ определения приближенных значений действительных корней – строят график функции f(x) и отмечают точки пересечения его с осью ОХ. Построение графиков часто удается упростить, заменив уравнение f(x)=0 равносильным ему уравнением Метод простых итераций для нелинейного уравнения в маткаде, где функции f1(x) и f2(x) — более простые, чем функция f(x). В этом случае следует искать точку пересечения этих графиков.

    Пример 1.Графически отделить корни уравнения x lg x = 1. Перепишем его в виде равенства lg x=1/xи найдем абсциссы точек пересечения логарифмической кривой y = lg x и гиперболы y = 1/x (рис. 5). Видно, что единственный корень уравнения Метод простых итераций для нелинейного уравнения в маткаде.

    Метод простых итераций для нелинейного уравнения в маткадеРис. 5. Метод графического отделения корней

    Реализация классических приближенных методов решения в пакете MathCAD.

    Метод половинного деления

    Отрезок, на концах которого функция принимает значения разного знака, делится пополам и, если корень лежит правее центральной точки, то к центру подтягивается левый край, а если – левее, то правый край. Новый суженный отрезок снова делится пополам и процедура повторяется. Этот метод прост и надежен, всегда сходится (хотя часто медленно – расплата за простоту!). Программная реализация его в пакете MathCAD рассмотрена в лабораторной работе №7 данного пособия.

    Метод хорд

    В качестве последовательных приближений к корню уравнения принимаются значения х1, х2, . хn точек пересечения хорды АВ с осью абсцисс (рис. 6).

    Уравнение хорды AB имеет вид: Метод простых итераций для нелинейного уравнения в маткаде. Для точки пересечения ее с осью абсцисс (х=х1, y=0) имеем:

    Метод простых итераций для нелинейного уравнения в маткаде.

    Пусть для определенности кривая у = f(x) будет выпукла вниз и, следовательно, расположена ниже своей хорды АВ, т.е. на отрезке [a,b] f ²(x)>0. Возможны два случая: f(а)>0 (рис. 6, а) и f(а) 3 –0,2x 2 –0,2х–1,2 = 0 с точностью e= 0,01. (Точный корень уравнения x = 1,2).

    Для организации итерационных вычислений в MathCAD документе используется функция until(a, z),котораявозвращает значение величины z, пока выражение a не становится отрицательным.

    Метод простых итераций для нелинейного уравнения в маткаде

    Метод простых итераций для нелинейного уравнения в маткаде

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

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

    Метод простых итераций для нелинейного уравнения в маткадеРис. 7. Метод касательных Ньютона

    При этом не обязательно задавать отрезок [а, b], содержащий корень уравнения), а достаточно лишь задать начальное приближение корня x=х0, которое должно находиться на том же конце интервала [а, b], где знаки функции и ее второй производной совпадают.

    Уравнение касательной, проведенной к кривой y = f(x) через точку В0 с координатами х0 и f(х0), имеет вид: Метод простых итераций для нелинейного уравнения в маткаде

    Отсюда найдем следующее приближение корня х1как абсциссу точки пересечения касательной с осью Ох (y = 0): Метод простых итераций для нелинейного уравнения в маткаде

    Аналогично могут быть найдены и последующие приближения как точки пересечения с осью абсцисс касательных, проведенных в точках В1, В2 и так далее. Формула для (i +1) приближения имеет вид:

    Метод простых итераций для нелинейного уравнения в маткаде

    Условием окончания итерационного процесса является неравенство ïf(xi)ï 1), то процесс итерации может быть расходящимся (рис. 8, в).

    Пример 4. Решитьуравнение x 3 – x – 1 = 0 методом простой итерации с точностью e = 10 -3 . Реализация этой задачи представлена следующим MathCAD документом.

    Метод простых итераций для нелинейного уравнения в маткаде

    Реализация приближенных методов решения встроенными функциями MathCAD

    Использование функции root

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

    📸 Видео

    1.1 Решение нелинейных уравнений метод деления отрезка пополам (бисекций) Мathcad15Скачать

    1.1 Решение нелинейных уравнений метод деления отрезка пополам (бисекций) Мathcad15

    8 Метод простой итерации Ручной счет Решение системы линейных уравнений СЛАУСкачать

    8 Метод простой итерации Ручной счет Решение системы линейных уравнений СЛАУ

    Средство для решения систем уравнений в MathCAD 14 (29/34)Скачать

    Средство для решения систем уравнений в MathCAD 14 (29/34)
    Поделиться или сохранить к себе: