
- •Содержание
- •Лекция № 1 понятие о численных методах. Решение нелинейных уравнений
- •Понятие о численных методах.
- •1.2 Решение нелинейных уравнений с одним неизвестным.
- •1.2.1. Отделение корня
- •1.2.2. Уточнение корня.
- •1.2.2.1 Метод деления отрезка пополам
- •Пример 1.1
- •1.2.2.2 Метод ньютона (метод касательных)
- •1.2.2.3 Метод итераций (простых)
- •Лекция №2 решение систем линейных алгебраических уравнений
- •2.1 Правило Крамера
- •2.2 Метод обратной матрицы.
- •2.3 Метод Гаусса
- •2.4 Итерационный метод гаусса – зейделя
- •Лекция 3. Интерполирование и аппроксимация функций
- •3.1 Линейная и квадратичная интерполяция
- •3.2 Интерполяция по лагранжу
- •3.3 Апроксимация многочленом по методу наименьших квадратов
- •Лекция № 4. Численное интегрирование
- •4.1 Метод прямоугольников
- •4.2 Метод трапеций
- •4.3 Метод симпсона
- •Если же подынтегральную функциюf(X)проинтерполировать полиномом Лагранжа 2-ойстепени (на отрезке [x0, x2] ):
- •4.4 Метод монте – карло
- •Лекция 5. Численное решение дифференциальных уравнений
- •5.1 Метод эйлера
- •5.2 Метод рунге—кутта
- •5.3 Метод прогноза и коррекции
- •Лекция 6. Методы оптимизации
- •6.1 Метод катящегося шарика
- •6.2 Метод золотого сечения
- •6.3 Метод градиентного спуска (уже для нескольких переменных)
- •6.4 Задача линейного программирования
- •Литература:
- •Вопросы к экзамену по численным методам.
Какую работу нужно написать?
1.2 Решение нелинейных уравнений с одним неизвестным.
Всякое уравнение с одним неизвестным может быть записано в виде:
f(х)=0
Предполагаем, что f(х) — непрерывная функция. И задача решения такого уравнения сводится, таким образом, к поиску х, при котором f(х) обращается в нуль.
Нелинейные уравнения условно делятся на:
— алгебраические уравнения (f(х) – алгебраические функции)
трансцендентные уравнения (тригонометрические, показательные, логарифмические и др.)
Методы решений нелинейных уравнений делят на:
— прямые (когда корни получаются в виде конечной формулы)
— итерационные (методы последовательных приближений)
Об итерационных методах и пойдет речь, (это и есть приближенные – численные – методы решений нелинейных уравнений)
Поиск корня с помощью итерационного метода распадается на три этапа:
1) выбор отрезка [а, b], на котором имеется корень
по теореме Больцано – Коши , если f (а) *f (в) < 0 (т. е. функция на разных концах отрезка имеет разные знаки), то существует с [a,b] и f(c) =0 (это — отделение корня)
2)выяснение того, что корень единственный
(это будет, если f (х))монотонна на [a,b], т. е. при f ' (х)>0 — монотонно возрастает, при f ' (х)<0 — монотонно убывает)
3) построение процесса, позволяющего сузить границы выделенного
отрезка, то есть позволяющего найти приближенные значения корня с любой заданной точностью (это — уточнение корня)
1.2.1. Отделение корня
Рассмотрим на примере: x3 – x – 1 = 0
f (х) = x3 – x – 1 и f ' (х)= 3x2 – 1 — непрерывны.
Определим интегралы монотонности, т. е. найдем точки перегиба (т. е. те точки, где
f ' (х)= 3x2 – 1 = 0 )
Эти
точки
Таким образом, вся область определения функций (— , + ) разобьётся на три интервала:
(— , — 1/3 ] , (— 1/3, +1/3 ) , [1/3, + ).
Вычислим значения функции на краях этих интервалов:
— не
содержит корней (функция на концах
имеет одинаковые знаки)
не
содержит корней
А
этот интервал содержит корни
И поскольку этот интервал является монотонным для f(х), то корень там единственный. Сузим интервал (до конечных точек): [1, 2],
f (1) = 1— 1—1= —1 < 0, f (2) = 8 –2 –1= 5> 0
Итак, на [1, 2] есть корень и единственный.
Отделение корня можно производить и графически.
Для этого уравнение f (х) = 0 сводят к уравнению: 1 (x) = 2 (x)
Затем строят графики y1=1 (x) и у2 =2 (x) и абсцисса точки их пересечения – есть корень. В рассмотренном примере уравнение принимало бы вид:
x3 = x + 1
и строили бы графики кубической параболы (y =x3 ) и прямой ( y=x + 1).
При этом, кстати, может быть обнаружено несколько таких точек пересечения (т.к. уравнение может иметь несколько корней) и, значит, будет отделено несколько отрезков.
Если рассматривать пример уравнения cos(x) =x2, то получим график, приведенныq на рис 1.
y=x2
y=cos
x
Видно, что данное уравнение имеет два корня: на отрезке [- π/2, 0] и [0,π/2]
Но для отделения корня можно активно использовать компьютер. При этом будем использовать такие очевидные положения:
1) если непрерывная на отрезке [a,b] функцияF(x) принимает на его концах значения разных знаков (т.е.F(a)•F(b)<0), то уравнениеF(x)=0имеет на этом отрезке, по крайней мере, один корень
2) если функция F(x)к тому же еще и строго монотонна, то корень на отрезке[a,b]единственный.
Итак, процедуру отделения корней такова. Пусть, все интересующие нас корни уравнения находятся на отрезке [A,B].Будем вычислять значенияF(x),начиная с точкиx=A,двигаясь вправо с некоторым шагомh. Как только обнаружится пара соседних значенийF(x),имеющих разные знаки, так соответствующие значения аргументаx(предыдущее и последующее) можно считать концами отрезка, содержащего корень. Схема такого алгоритма приведена нарис.2.
начало
ввод a,b,h
x1 := a
x2 := x1 + h
y1:= f(x1)
Да
Нет
x2<b
y2:=f(x2)
Да Нет
y1•y2<0
вывод:
x1, x2
x1
:= x2
x2 := x1+h
y1 := y2
конец
Рис.2 . Блок–схема алгоритма отделения корней.