рефераты конспекты курсовые дипломные лекции шпоры

Реферат Курсовая Конспект

Пусть уравнение (2.1) можно заменить эквивалентным ему уравнением

Пусть уравнение (2.1) можно заменить эквивалентным ему уравнением - раздел Философия, Математическая модель. Решение нелинейных уравнений. Решение систем линейных алгебраических уравнений X = (X). (2.4) Например, Уравнение - 0.5 = 0 Мо...

x = (x). (2.4)

Например, уравнение - 0.5 = 0 можно заменить эквивалентным ему уравнением x = 0.5sinx.

Выберем каким-либо образом начальное приближение x0. Вычислим значение функции (x) при x = x0 и найдем уточненное значение x1 = (x0). Подставим теперь x1 в уравнение (2.4) и получим новое приближение x2 = (x1) и т. д. Продолжая этот процесс неограниченно, получим последовательность приближений к корню:

xn+1 = (xn). (2.5)

Формула (2.5) является расчетной формулой метода простых итераций.

Если последовательность {xn} сходится при n, т. е. существует

x* = xn , (2.6)

и функция (x) непрерывна, то, переходя к пределу в (2.5) и учитывая (2.6), получим:

x* = xn = (x n -1) = (xn -1) = (x*).

Таким образом, x* = (x*), следовательно, x* - корень уравнения (2.4).

Сходимость метода. Сходимость метода простых итераций устанавливает следующая теорема.

Теорема 2.2. Если в интервале, содержащем корень x* уравнения (2.4), а также его последовательные приближения x0, x1, …, xn, …, вычисляемые по формуле (2.5), выполнено условие:

|(x)| q < 1, (2.7)

то x* = xn.

т. е. итерационный процесс сходится и справедлива следующая оценка погрешности:

|xn - x*| qn|x0 - x*| (2.8)

Оценка (2.8) является априорной. Она показывает, что метод простой итерации сходится со скоростью геометрической прогрессии с знаменателем q. Чем меньше q, тем выше скорость сходимости.

Как следует из теоремы 2.2, условие (2.7) является достаточным для сходимости метода простых итераций. Его выполнение гарантирует сходимость процесса (2.5), но невыполнение условия (2.7), вообще говоря, не означает, что итерационный процесс будет расходиться.

На рис. 2.3 - 2.6 показаны четыре случая взаимного расположения линий y = x и y = (x) и соответствующие итерационные процессы.

Рис. 2.3 и 2.4 соответствуют случаю |(x)| < 1, и итерационный процесс сходится. При этом, если (x) > 0 (рис. 2.3), сходимость носит односторонний характер, а если (x) < 0 (рис. 2.4), сходимость носит двусторонний, колебательный характер. Рис. 2.5 и 2.6 соответствуют случаю |(x)| > 1 - итерационный процесс расходится. При этом может быть односторонняя (рис. 2.5) и двусторонняя (рис 2.6) расходимость.

Рис. 2.3 Рис. 2.4 Рис. 2.5

Рис. 2.6

Погрешность метода. Если известна величина q в условии (2.7), то применима следующая апостериорная оценка погрешности:

|xn - x*| |xn - xn - 1|, n > 1. (2.9)

Критерий окончания. Из оценки (2.9) вытекает следующий критерий окончания итерационного процесса. Вычисления следует продолжать до выполнения неравенства

|xn - xn - 1| < .

Если это условие выполнено, то можно считать, что xn является приближением к x* с точностью .

Если q 0.5, то можно пользоваться более простым критерием окончания:

|xn - xn - 1| < . (2.10)

Пример 2.2.

Используем метод простой итерации для решения уравнения f(x) = sin x - x2 = 0 с точностью = 0.001.

Преобразуем уравнение к виду (2.4):

x = , т. е. (x)= .

Нетрудно убедиться, что корень уравнения находится на отрезке [/6, /3]. Например, вычислив значения f(x) на концах отрезка, получим: f(/6)> 0, а f(/3)< 0, т. е. функция на концах отрезка имеет разные знаки, что в соответствии с теоремой 2.1 указывает на то, что внутри отрезка есть корень. Расположение корня наглядно иллюстрирует рис.2.7.

Рис. 2.7

Подсчитаем, первую и вторую производные функции (x):

(x) = , "(x) = .

Так как "(x) > 0 на отрезке [/6, /3], то производная (x) монотонно возрастает на этом отрезке и принимает максимальное значение на правом конце отрезка, т. е. в точке /3. Поэтому, справедлива оценка:

| (x)| | (/3)| 0.312.

Таким образом, условие (2.7) выполнено, q < 0.5, и можно воспользоваться критерием окончания вычислений в виде (2.10). В табл. 2.2 приведены приближения, полученные по расчетной формуле (2.5). В качестве начального приближения выбрано значение x0 = 1.

Таблица 2.2

 
n xn  
0.8415 0.8861 0.8742 0.8774 0.8765  
     

Критерий окончания выполняется при n = 5, |x5 - x4| < 0.001. Сходимость двусторонняя, качественный характер такой сходимости представлен на рис. 2.4. Приближенное значение корня с требуемой точностью x* 0.8765.

– Конец работы –

Эта тема принадлежит разделу:

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

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

Если Вам нужно дополнительный материал на эту тему, или Вы не нашли то, что искали, рекомендуем воспользоваться поиском по нашей базе работ: Пусть уравнение (2.1) можно заменить эквивалентным ему уравнением

Что будем делать с полученным материалом:

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

Все темы данного раздела:

Корректность
Определим вначале понятие устойчивости решения. Решение задачи y* называется устойчивым по исходным данным x*, если оно зависит от исходных данны

Вычислительные методы
Под вычислительными методами будем понимать методы, которые используются в вычислительной математике для преобразования задач к виду, удобному для реализации на ЭВМ. Подробнее с различными к

ЛЕКЦИЯ 9
Тема: Элементы теории погрешностей Определение. Пусть u и — точное и приближенное значение некоторой величины соответственно. Тогда абсолютной погрешностью п

Метод деления отрезка пополам является самым простым и надежным способом решения нелинейного уравнения.
Пусть из предварительного анализа известно, что корень уравнения (2.1) находится на отрезке [a0, b0], т. е. x*[a0, b0], так, что f(x

Метод Ньютона (метод касательных)
вать следующий критерий окончания итераций метода Ньютона. При заданной точности > 0 вычисления нужно вести до тех пор, пока не будет выполнено неравенство |xn - xn -

Метод ложного положения
Рассмотрим еще одну модификацию метода Ньютона. Пусть известно, что простой корень x* уравнения f(x) = 0 находится на отрезке [a, b] и на одном из ко

Постановка задачи
Требуется найти решение системы линейных уравнений: a11x1 + a12 x2 + a13x3

Метод исключения Гаусса с выбором главного элемента по столбцу
Хотя метод Гаусса является точным методом, ошибки округления могут привести к существенным погрешностям результата. Кроме того исключение по формулам (3.7) нельзя проводить, если элемент главной ди

Вычисление определителя методом исключения Гаусса
Из курса линейной алгебры известно, что определитель треугольной матрицы равен произведению диагональных элементов. В результате метода исключений Гаусса система линейных уравнений (3.2) с квадратн

Вычисление обратной матрицы методом исключения Гаусса
Обратной матрицей к матрице A называется матрица A-1, для которой выполнено соотношение: A A-1 = E, (3.18) где

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

Метод Зейделя
Модификацией метода простых итераций Якоби можно считать метод Зейделя. В методе Якоби на (k+1)-ой итерации значения x, i = 1, 2, …, n. вычисляются подстановкой

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

Приближение функции многочленами Тейлора
Пусть функция y = f(x) определена в окрестности точки a и имеет в этой окрестности n + 1 производную. Тогда в этой окрестности справедлива формула Тейлора:

Интерполяция функции многочленами Лагранжа
Рассмотрим другой подход к приближению функции многочленами. Пусть функция y = f(x) определена на отрезке [a, b] и известны значения этой функции в некоторой системе узл

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

Постановка задачи численного интегрирования
Далеко не все интегралы можно вычислить по известной из математического анализа формуле Ньютона - Лейбница: I == F(b) - F(a), (5.1) где

Метод прямоугольников
Формулу прямоугольников можно получить из геометрической интерпретации интеграла. Будем интерпретировать интеграл как площадь криволинейной трапеции, ограниченной графиком функции y = f(x

Метод трапеций
Выведем формулу трапеций так же, как и формулу прямоугольников, из геометрических соображений. Заменим график функции y = f(x) (рис.5.1) ломаной линией (рис.5.7), полученной следующим

Метод Симпсона (метод парабол)
Заменим график функции y = f(x) на отрезке [xi, xi+1], i = 0, 2, … , n - 1, параболой, проведенной через точки (xi

Правило Рунге практической оценки погрешности
Оценки погрешности по формулам (5.4), (5.8) и (5.12) являются априорными. Они зависят от длины элементарного отрезка h, и при достаточно малом h справедливо приближенное равенство:

Постановка задачи Коши
Известно, что обыкновенное дифференциальное уравнение первого порядка имеет вид: y (t) = f(t, y(t)). (6.1) Решением уравнения (6.1) являе

Метод Эйлера
Простейшим методом решения задачи Коши является метод Эйлера. Будем решать задачу Коши y (t) = f(t, y(t)). y(t0

Модифицированные методы Эйлера
Первый модифицированный метод Эйлера. Суть этого метода состоит в следующем. Сначала вычисляются вспомогательные значения искомой функции y в точках t = ti + с помощ

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

Хотите получать на электронную почту самые свежие новости?
Education Insider Sample
Подпишитесь на Нашу рассылку
Наша политика приватности обеспечивает 100% безопасность и анонимность Ваших E-Mail
Реклама
Соответствующий теме материал
  • Похожее
  • Популярное
  • Облако тегов
  • Здесь
  • Временно
  • Пусто
Теги