Решение уравнения с одним неизвестным

Содержание

Слайд 2

пример: f(x) = 3*sin(2*x)-1.5*x-1=0

f=Inline(‘3*sin(2*x)-1.5*x-1‘)
a=input(‘a=‘);
b=input(‘b=‘);
h=input(‘h=‘);
x=a:h:b;
plot(x,f(x)); grid
xlabel(‘x’); ylabel(‘f(x)’)

пример: f(x) = 3*sin(2*x)-1.5*x-1=0 f=Inline(‘3*sin(2*x)-1.5*x-1‘) a=input(‘a=‘); b=input(‘b=‘); h=input(‘h=‘); x=a:h:b; plot(x,f(x)); grid xlabel(‘x’); ylabel(‘f(x)’)

Слайд 3

Уточнение корня на отрезке [a,b], в котором локализован только один корень, осуществляется

Уточнение корня на отрезке [a,b], в котором локализован только один корень, осуществляется
итерационными методами, в которых последовательно, шаг за шагом, производится уточнение начального приближения корня. Итерацией называется совокупность вычислительных операций, приводящих к новому приближенному значению корня. Если каждое последующее значение x(k) (k=1,2,3,…) находится все ближе к точному значению, говорят, что метод сходится. В противном случае метод расходится. Для реализации итерационного процесса должны быть заданы начальное приближение x(0) и точность ε, с которой найти решение уравнения. Условие окончание имеет вид: |x(k)-x(k-1)| ≤·ε.Все методы можно разделить на две группы: с условной и безусловной сходимостью.

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

В этом методе на каждой итерации новое приближение определяется как: x(k)=(a(k-1)+b(k-1))/2, где к – номер итерации.

Алгоритм
Задаем функцию f(x), отрезок [a(0),b(0)], точность ε и k=1.
Вычисляем приближение x(k)=(a(k-1)+b(k-1))/2
Определяем новый отрезок [a(k),b(k)]. Проверяем, если f(a(k-1))*f(x(k))>0, то a(k)=x(k) и b(k)=b(k-1), иначе a(k)=a(k-1) и b(k)=x(k).
Проверяем условие окончания, если |b(k)-a(k)| ≤·2ε, то за ответ принимаем значение равное x=(a(k)+b(k))/2 и переходим на пункт 5, иначе k=k+1 и переходим на пункт 2.
выводим x и f(x).

Методы с безусловной сходимостью

Слайд 4

Блок-схема

начало

f(a)*f(x)>0

x, f(x)

a, b, ε || f(x)

b := x

a=x

x:= (b+a)/2

| b-a | ≤

Блок-схема начало f(a)*f(x)>0 x, f(x) a, b, ε || f(x) b :=

конец

x := (b+a)/2

нет

да

Слайд 5

Решим предыдущий пример при a= -1.6 b= -1.2 и ε= 0.01 т.е.

Решим предыдущий пример при a= -1.6 b= -1.2 и ε= 0.01 т.е.
2ε = 0.02

x= –1,38±0.01 f(x) = –0,038 (невязка)

Слайд 6

Метод простых итераций

В этих методах исходное уравнение f(x)=0 преобразуется к эквивалентному

Метод простых итераций В этих методах исходное уравнение f(x)=0 преобразуется к эквивалентному
виду x=ϕ(x). Тогда на каждой итерации новое приближение будем определять как:
x(1) = ϕ (x(0)), x(2) = ϕ (x(1)), x(3) = ϕ (x(2)),….., т.е. x(k)=ϕ(x(k-1)), k=1,2,3… .
За x(0) принимают любое число на заданном отрезке [a;b]. Вид функции ϕ(x) определим исходя из достаточного условия сходимости, которое записывается как: |ϕ’(x)| < 1, для всех значений x отрезка[a;b], т.е. максимальная производная на заданном отрезке должна быть меньше единицы.

Общий подход для получения итерационной формулы x=ϕ(x)
Помножим обе части уравнения f(x)=0 на множитель, и прибавим к обеим частям по x, тогда итерационная формула будет иметь вид:
x = x + βf(x) = ϕ(x)

Методы с условной сходимостью

Для уравнения x2-5=0 можно положить ϕ(x)=5/x или ϕ(x)=(1/2)(x+5/x) и соответствующие итерационные формулы будут иметь вид x(k)=5/x(k-1) и x(k)=(1/2)(x(k-1)+5/x(k-1)).

А во втором сходится

В первом случаи метод расходится

Слайд 7

Блок-схема

Определить множитель β можно из достаточного условия сходимости.
|ϕ’(x)| < 1 ϕ’(x)

Блок-схема Определить множитель β можно из достаточного условия сходимости. |ϕ’(x)| Мы должны
= 1 + βf’(x) |1 + βf’(x)| < 1 -1 < 1 + βf’(x) < 1 -2 < βf’(x) < 0.
Мы должны выбрать максимальную по модулю производную |f’(x)| на заданном отрезке.
|f’(b)|>|f’(a)| β = -2/f’(b),иначе β = -2/ f’(a)

Слайд 8

Пример: f(x) = 3sin(2x)-1.5x-1 f'(x)=6cos(2x)-1.5 ε=0.01 a = -1,6 b = -1,2

Пример: f(x) = 3sin(2x)-1.5x-1 f'(x)=6cos(2x)-1.5 ε=0.01 a = -1,6 b = -1,2

f'(a) = -7,489 f'(b) = -5,924 β = 0,267 ≈ 0.2
x(k) = x(k-1) + β(3sin(2x(k-1))-1.5x(k-1)-1)

Ответ: x = -1,38±0.01 f(x) = -0,020

Слайд 9

Метод Ньютона или касательных

Пусть известно некоторое приближение x(k-1) к решению x* уравнения

Метод Ньютона или касательных Пусть известно некоторое приближение x(k-1) к решению x*
f(x)=0.

Тогда исходное уравнение можно записать в виде:

f(x(k-1)+∆x(k-1))=0

где ∆x(k-1)= x* -x(k-1) и x* = x(k-1)+ ∆x(k-1)

Разложим функцию в ряд Тейлора и ограничимся линейными членами.

f(x(k-1)+∆x(k-1)) = f(x(k-1))+ f′(x(k-1))∆x(k-1) = 0

откуда

Полученное значение принимаем за новое приближение к решению. Тогда итерационную формулу запишем как:

Слайд 10

Графическая иллюстрация.

За начальное приближение к корню x(0) принимаем одну из границ

Графическая иллюстрация. За начальное приближение к корню x(0) принимаем одну из границ
отрезка [a; b], содержащего один корень.

На каждой итерации, за новое приближение к корню x(k) принимается точка пересечения касательной к графику, построенной в точке f(x(k-1)) с осью абсцисс x:

Слайд 11

алгоритм

Задаем функцию f(x) отрезок [a;b] и точность ε. За начальное приближение x

алгоритм Задаем функцию f(x) отрезок [a;b] и точность ε. За начальное приближение
принимаем одну из границ заданного отрезка [a,b] x=a.
Вычисляем значение шага h= f(x)/f′(x) и новое приближение, как x = x-h.
Проверяем условие окончания если | h | ≤ ε, то выводим последнее значение x и f(x). Иначе перейдем на пункт 2

Блок-схема

Имя файла: Решение-уравнения-с-одним-неизвестным.pptx
Количество просмотров: 41
Количество скачиваний: 0