ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.11.2023
Просмотров: 25
Скачиваний: 3
ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
Найдем корни уравнения:
x3+3·x-1=0
ε = 0.001
Используем для этого Метод итераций.
Одним из наиболее эффективных способов численного решения уравнений является метод итерации. Сущность этого метода заключается в следующем. Пусть дано уравнение f(x)=0.
Заменим его равносильным уравнением x=φ(x).
Выберем начальное приближение корня x0 и подставим его в правую часть уравнения. Тогда получим некоторое число x1=φ(x0).
Подставляя теперь в правую часть вместо x0 число x1 получим число x2=φ(x1). Повторяя этот процесс, будем иметь последовательность чисел xn=φ(xn-1)
Если эта последовательность сходящаяся, то есть существует предел ξ = lim(xn), то переходя к пределу в равенстве и предполагая функцию φ(x) непрерывной найдем lim(xn) = φ(lim(xn-1)), n → ∞ или ξ=φ(ξ).
Таким образом, предел ξ является корнем уравнения и может быть вычислен по формуле с любой степенью точности.
Находим первую производную:
Решение.
Представим уравнение в форме:
x = x - λ(x3+3*x-1)
Найдем максимальное значение производной от функции f(x) = x3+3*x-1
max(3•x2+3) ≈ 6
Значение λ = 1/(6) ≈ 0.1667
Таким образом, решаем следующее уравнение:
x-0.1667*(x3+3*x-1) = 0
F(0)=-1; F(1)=3
Поскольку F(0)*F(1)<0 (т.е. значения функции на его концах имеют противоположные знаки), то корень лежит в пределах [0;1].
Остальные расчеты сведем в таблицу.
N | x | F(x) |
1 | 0 | -1 |
2 | 0.1667 | -0.4953 |
3 | 0.2493 | -0.2367 |
4 | 0.2887 | -0.1098 |
5 | 0.307 | -0.05 |
6 | 0.3154 | -0.02257 |
7 | 0.3191 | -0.01015 |
8 | 0.3208 | -0.00455 |
Ответ: x = 0.32080977288594; F(x) = -0.00455
Сходимость:
Найдем корни уравнения:
x3+3·x-1=0
ε = 0.001
Используем для этого Метод половинного деления (метод дихотомии)..
Считаем, что отделение корней произведено и на интервале [a,b] расположен один корень, который необходимо уточнить с погрешностью ε.
Итак, имеем f(a)f(b)<0. Метод дихотомии заключается в следующем.
Определяем половину отрезка c=1/2(a+b) и вычисляем f(c). Проверяем следующие условия:
1. Если |f(c)| < ε, то c – корень. Здесь ε - заданная точность.
2. Если f(c)f(a)<0, то корень лежит в интервале [a,c].
3. Если f(c)f(b)<0, то корень лежит на отрезке[c,b].
Продолжая процесс половинного деления в выбранных подынтервалов, можно дойти до сколь угодно малого отрезка, содержащего корень ξ.
Так как за каждую итерацию интервал, где расположен корень уменьшается в два раза, то через n итераций интервал будет равен:
bn-an=1/2n(b-a)
В качестве корня ξ. возьмем 1/2(an+bn). Тогда погрешность определения корня будет равна (bn – an)/2. Если выполняется условие:
(bn – an)/2 < ε
то процесс поиска заканчивается и ξ = 1/2(an+bn).
Решение.
Число шагов, необходимых для достижения заданной точности определяется неравенством:
F(-10)=-1031; F(10)=1029
Поскольку F(-10)*F(10)<0 (т.е. значения функции на его концах имеют противоположные знаки), то корень лежит в пределах [-10;10].
Итерация 1.
Находим середину отрезка: c = (-10 + 10)/2 = 0
F(x) = -1
F(c) = -1031
Поскольку F(c)*F(b) < 0, то a=0
Итерация 2.
Находим середину отрезка: c = (0 + 10)/2 = 5
F(x) = 139
F(c) = -1
Поскольку F(c)*F(a) < 0, то b=5
Итерация 3.
Находим середину отрезка: c = (0 + 5)/2 = 2.5
F(x) = 22.125
F(c) = 139
Поскольку F(c)*F(a) < 0, то b=2.5
Итерация 4.
Находим середину отрезка: c = (0 + 2.5)/2 = 1.25
F(x) = 4.703
F(c) = 22.125
Поскольку F(c)*F(a) < 0, то b=1.25
Остальные расчеты сведем в таблицу.
N | c | a | b | f(c) | f(x) | ε |
1 | 0 | 0 | 10 | -1031 | -1 | 10 |
2 | 5 | 5 | 10 | -1 | 139 | 5 |
3 | 2.5 | 2.5 | 5 | 139 | 22.125 | 2.5 |
4 | 1.25 | 1.25 | 2.5 | 22.125 | 4.7031 | 1.25 |
5 | 0.625 | 0.625 | 1.25 | 4.7031 | 1.1191 | 0.625 |
6 | 0.3125 | 0.3125 | 0.625 | 1.1191 | -0.03198 | 0.3125 |
7 | 0.4688 | 0.4688 | 0.625 | -0.03198 | 0.5092 | 0.1563 |
8 | 0.3906 | 0.3906 | 0.4688 | 0.5092 | 0.2315 | 0.07813 |
9 | 0.3516 | 0.3516 | 0.3906 | 0.2315 | 0.09814 | 0.03906 |
10 | 0.332 | 0.332 | 0.3516 | 0.09814 | 0.0327 | 0.01953 |
11 | 0.3223 | 0.3223 | 0.332 | 0.0327 | 0.000266 | 0.00977 |
Таким образом, в качестве корня можно принять:
x=(0.3125+0.3223)/2 = 0.3174
Ответ:x = 0.3174; F(x) = 0.000266
Количество итераций, N = 11
Параметр сходимости.
Сходимость метода дихотомии линейная с коэффициентом α = 0.5.
Найдем корни уравнения:
x3+3·x-1=0
ε = 0.001
Используем для этого Метод Ньютона.
Пусть корень ξ уравнения f(x)=0 отделен на отрезке [a,b]. Предположим мы нашли (n-1)-ое приближение корня xn-1. Тогда n-ое приближение xn мы можем получить следующим образом. Положим:
xn = xn-1 + hn-1
Раскладывая в ряд f(x=ξ) в точке xn-1, получим:
f(xn) = f(xn-1+hn-1) = f(xn-1) + f’(xn-1)hn-1=0
Отсюда следует:
Подставим hn-1 в формулу, получим:
Геометрически метод Ньютона эквивалентен замене дуги кривой y=f(x) касательной, проведенной в некоторой точке кривой.
Находим первую производную:
Находим вторую производную:
d2F/dx2 = 6•x
Решение.
F(-10)=-1031; F(10)=1029
Поскольку F(-10)*F(10)<0 (т.е. значения функции на его концах имеют противоположные знаки), то корень лежит в пределах [-10;10].
Вычисляем значения функций в точке a = -10
f(-10) = -1031
f''(-10) = -60
Критерий остановки итераций.
|f(xk)| < εm1
или
где M2 = max|f "(x)|, m1 = min|f'(x)|.
Поскольку f(a)*f''(a) > 0, то x0 = a = -10
Остальные расчеты сведем в таблицу.
N | x | F(x) | dF(x) | h = f(x) / f'(x) |
1 | -10 | -1031 | 303 | -3.4026 |
2 | -6.5974 | -307.9432 | 133.5755 | -2.3054 |
3 | -4.292 | -92.9385 | 58.2631 | -1.5952 |
4 | -2.6968 | -28.704 | 24.8185 | -1.1566 |
5 | -1.5403 | -9.2749 | 10.1172 | -0.9167 |
6 | -0.6235 | -3.113 | 4.1663 | -0.7472 |
7 | 0.1237 | -0.6271 | 3.0459 | -0.2059 |
8 | 0.3296 | 0.02446 | 3.3258 | 0.00735 |
9 | 0.3222 | 5.3E-5 | 3.3114 | 1.6E-5 |
Ответ: x = 0.3222 - 5.3E-5 / 3.3114 = 0.32218535470104; F(x) = 0
Параметр сходимости.
Скорость сходимости метода Ньютона квадратичная с коэффициентом
α = M2/2m1, где M2 = max|f "(x)|, m1 = min|f'(x)|.0>0>0>