Численные методы

Численные методы представляют собой набор алгоритмов, позволяющих получать приближенное (численное) решение математических задач. Два вида погрешностей, возникающих при решении задач. Нахождение нулей функции. Метод половинного деления. Метод хорд.

Рубрика Математика
Вид курс лекций
Язык русский
Дата добавления 06.03.2009
Размер файла 81,2 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

ЛЕКЦИЯ №1

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

Погрешности, возникающие при решении задач, бывают двух видов:

1)абсолютная

p - p? , где p - точное значение, p? - не точное.

2)относительная

Эмпирические данные:

Погрешности Случайные Ошибки

измерительного помехи набора

прибора

1) Нахождение нулей функции;

2) Системы линейных и нелинейных уравнений;

3) Приближение функции. Интерполяция. Экстраполяция.

4) Решение дифференциальных уравнений.

5) Расчет собственных значений и собственных векторов матриц.

НАХОЖДЕНИЕ НУЛЕЙ ФУНКЦИИ

Общая постановка задачи

Дана некоторая функция f(х). Необходимо найти хотя бы одно значение х, при котором f(х)=0.

Этапы:

1) Отделение корней.

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

содержится единственный корень функции.

2) Уточнение корня при помощи одного из численных методов на каждом из выбранных отрезков.

Нуль функции - точка пересечения графика функции с осью Ох.

Непрерывность f(х) в точке х0:

Производная функции: f' =

Физический смысл: f'(х0)- скорость

Геометрический смысл: f'(х0)-тангенс угла наклонной касательной к графику функции, проведенной в данной точке.

Если функция дифференцируема в точке, то она непрерывна. Обратное не верно.

Предел функции в точке:

x: | x-x0| <

е >0 (е)

| f(x) - A| < е

Градиент функции - это вектор.

Геометрический смысл : показывает направление локального возрастания функции в данной точке .

1) Наблюдаем смену знака функции.

2) Исследуем функцию на монотонность.

Теорема №1: если функция f(x) непрерывна на отрезке [a, b] и в концах отрезка принимает значения разных знаков, то на этом отрезке функция имеет хотя бы один корень.

f(x) C[a, b]

f(a) * f(b) < 0 > [a, b] f()=0

Теорема№2: если функция непрерывна и монотонна на отрезке и в концах отрезка принимает значения разных знаков, то на этом отрезке существует только единственный корень функции.

f(x) C[a, b], f ( ) и f(a) * f(b) < 0>[a, b] f() = 0

МЕТОД ПОЛОВИННОГО ДЕЛЕНИЯ

Дано: f(x) непрерывна на [a,b], на [a,b] существует динственный корень f(x)=0, е

1) Делим отрезок пополам. Получаем точку

с= (b + a)/2.

Если f(a) * f(c) < 0,то b:=c.

Если f(b) * f(c) < 0,то а:=с

2) Продолжаем делить [a, b] на 2, пока|b-a| > е, где е- заданная точность.

ЛЕКЦИЯ №2

МЕТОД ХОРД

Дано: 1) f(x) C''[a, b]

2) f(a) * f(b) < 0

3) f'(x) и f''(x) знакопостоянна на отрезке [a, b].

4) е, чтобы получить f(x)=0

1) f(b) 2)

f'(x) >

0 f'(x) > 0

f''(x) >

0 f''(x) < 0

f(a) a x

3) 4)

f'(x)

<0 f'(x) <0

f''(x)

<0 f''(x) > 0

(2.1)

x1(x1,f(x1))

b - неподвижный конец отрезка.

Для случаев 1), 3)

Для случаев 2), 4)

Можем ввести некоторую с:

(2.2)

(2.3)

Алгоритм:

1) Вычисляем неподвижный конец отрезка секущих по формуле(2.3)

2) Находим первое приближение к корню по формуле (2.1)

3) Находим первое приближение к корню по формуле (2.2) до тех пор, пока модуль разности двух последних приближений не станет меньше заданной точности. В этом случае, значением корня является последнее приближение.

МЕТОД КАСАТЕЛЬНЫХ

Дано: 1) f(x) C''[a, b];

2) f(a)*f(b) < 0;

3) f'(x) и f''(x) знакопостоянны на [a, b];

4) е, чтобы решить уравнение f(x)=0

т. х0

y=f(x0)+f'(x0)(x-x0) -

уравнение касательной

a x2 x1 b

y=f(b)+f'(b)*(x-b)

(x1,0) : 0= f(b)+ f'(b)(x1-b)

x1=

x2=

xn+1= (2.4)

Второй подход (метод Ньютона):

-приближение

0 = f() = f(xn+hn) ? f(xn)+f'(xn)*hn

x0 = начальное приближение (2.5)

Алгоритм:

1) По формуле (2.5) находим первое приближение к корню х0 (начальное)

2) По формуле (2.4) находим последующее приближение к корню до тех пор, пока модуль разности двух последних приближений не станет заданной точности. В этом случае корень равен последнему приближению.

МЕТОД ИТЕРАЦИЙ

Дано: 1) f(x)C''[a,b]

2)f(a)*f(b)<0

3)f'(x) знакопостоянна

4)е, f(x)=0

Уравнение f(x)=0 заменяется уравнением вида x=ц(x)

ц(x)=x-f(x)*C (2.6)

Пока |xn+1-xn|<е

ц' >0

Cтроим последователь

Выбираем

Находим значение функции

x2= ц(x1), x3= ц(x2)

xn+1= ц(xn) (2.7)

Точка е, для которой выполняется е=f(е), называется неподвижной точкой метода итераций. Очевидно, что эта точка является корнем уравнения f(x)=0.

ц(е) е -f(x)* е

0 f(е)*C

f(е) 0

Достаточное условие: для того, чтобы метод итераций сходился достаточно чтобы:

1) ц(x) (2.8) - Функция является непрерывной и дифференцируемой на [a,b].

2) ц(x) значения - является необходимым условием

3) |ц(x)|<1 для всех

Константа С в формуле(2.6) подбирается таким образом, чтобы функция

ц(x) удовлетворяла условиям сходимости метода итераций.

Скорость сходимости метода Ньютона (касательных) выше сходимости метода секущих (хорд).

ЛЕКЦИЯ №3

МЕТОДЫ РЕШЕНИЯ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ

Общий вид алгебраического уравнения:

а0хn+ а1хn+1+…+ аn-1х+an=0, a00 (3.1)

n=1: а0х+a1=0, x=

n=2: а0х2+a1x+a2=0, x1,2=

Алгебраическое уравнение n- степени имеет ровно n корней.

Теорема Виета (обобщенная):

xn+xn-1+…+x+=0

x1+x2+…+xn=-; (3.2)

x1x2+x1x3+…+xn-1xn=;

x1x2x3…xn=(-1);

Пусть все корни уравнения (3.1) действительны, различны и удовлетворяют соотношениям:

|x1|>>|x2|>>…>>|xn| (3.3)

Преобразуем:

x1(1++…+)= x1=-; (3.4)

Подставим (3.4) : х2=- продолжая получим общую формулу

хk=-, k=1,n (3.5)

Корни уравнения, удовлетворяющие соотношения(3.3), называются отдельными. Задача состоит в том, чтобы по исходному уравнению построить такое уравнение, корни которого будут отделены.

yi=-xim

b0yn + b1yn-1+…+ bn-1y+bn=0 (3.6)

|x1|>|x2|>…>|xn|

Решив уравнение (3.6), корни которого являются отдельными, получим уравнения y1…yn

, i=2,n

Значит |yi-1|>>|xi|

МЕТОД ЛОБАЧЕВСКОГО

Для отделения корней Лобачевский предложил метод квадратирования - способ построения по исходному уравнению нового уравнения, кони которого связаны с корнями исходного следующим образом:

yi=-xi2 (3.7)

Процедура выполнения многократна, пока не достигнем серьёзной разницы модуля разности корней

b0(m)yn + b1(m)yn-1+…+ bn-1(m)y+ bn(m)=0 (3.8)

Пусть уравнение (3.8) получено в результате m-го шага квадрирования.

m=1 b0(1)=a02, b1(1)= a12=2 a0 a2

bk(1)=ak2-2ak-1ak+1+2ak-2ak+2….,k=0,n

При получении bk коэффициента , который рассчитывается как квадрат соответствующего коэффициента ak минус удвоенное произведение соседних коэффициентов с ak плюс удвоенное произведение следующей пары соседей , чередуя знаки, пока в число соседних коэффициентов не попадут а0 и аn.

m>1b0(m)=( b0(m-1))2, b1(m)=( b1(m-1))2-2b0(m-1)b2(m-1) (3.9)

bk(m)=( b0(m-1))2-2bk-1(m-1)bk-1(m-1)+2bk-2(m-1)bk+2(m-1)

Критерий остановки: bk(m)?( b0(m-1))2, k=0,n (3.10)

Получим корень: yi(m)=-xi2, i=1,n (3.11)

(3.11)-связь корней, полученных на m-шаге процесса квадрирования с корнями исходного уравнения.

yi на m-шаге : , отсюда

, i=1,n (3.12)

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

Признак наличия кратных корней: один или несколько коэффициентов > к половине квадрата коэффициента предыдущего шага.

МЕТОДЫ РЕШЕНИЯ СИСТЕМ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ

СЛАУ

Методы решения СЛАУ делятся на точные и приближенные. К точным методам относятся метод Гаусса, метод Крамера, метод обратной матрицы .

Существуют приближенные методы: метод итераций, Зейделя и т.д.

Общий вид СЛАУ:

(3.13)

Сколько переменных столько и ограничений на них.

пересечение прямых (точка)

пересечение плоскостей (прямая)

точка пересечения трех плоскостей

Т.о. геометрический смысл решения СЛАУ - точка пересечения гиперплоскостей в n-мерном пространстве.

Матрица :=А; B=; X=;

Cn*k*Dk*m=Zn*m , An*n*Bn*1=Xn*1 AX=B (3.14)

ЛЕКЦИЯ №4

МЕТОД ГАУССА

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

1 шаг: Выбираем в матрице А максимальный элемент по всем строкам и столбцам. Путем перестановки строк и столбцов ставим этот элемент на место а11. Теперь а11- главный элемент.

А>А12>…>Аn

Аn должна будет содержать ниже главной диагонали все нули.

, j =1,n ; b1 =b1/a11

Получим систему вида

, i=2,n , j=1,n

Получим А' х=В' и систему

Пусть а221 - максимальный по модулю элемент матрицы А1 по строкам i?2 и столбцам j?2. Если это не так, то добиваемся этого путем перестановки строк и столбцов.

А2:

В2: b12=b11; b22=b21/a221; bi2=bi1-b22-b22ai21

Пусть акк+1 максимальный по модулю элемент матрицы Ак, i?k, j?k.

Пусть на некотором шаге k<n элемент =0, матрица Вк имеет ? множество решений. Причем корни х1,…хк являются зависимыми, а корни хк+1,….xn - независимые.

Если хотя бы один элемент bik при i?k+1 0, то решения у системы нет.

Если была получена матрица Аn, то система имеет единственное решение.

Начинается обратный ход метода Гаусса.

МЕТОД КРАМЕРА

Определитель : det A=

det A==a11a22-a12a21

Минор Hij элемента матрицы aij представляет собой определитель, полученный из матрицы А путем вычеркивания i cтроки и j столбца.

Алгебраическое дополнение Аij элемента аij называется число, равное

Аij=(-1)i+j*Mij

Способы вычисления определителей

1. Привести определитель к треугольному виду (ниже главной диагонали все элементы=0). Достичь этого можно путем вычитания (сложения) строк определителя, умноженных на некоторое число. При перестановки строк/столбцов знак определителя меняется на противоположный. Определитель треугольного вида равен произведению элементов главной диагонали.

2. Рекуррентный способ основан на том, что определитель равен сумме произведений элементов строки/столбца на их алгебраические дополнения. Т.о. задача вычисления определителя n-го порядка сводится к вычислению n определителей n-1 порядка.

Наиболее целесообразно раскладывать определитель по той строке/столбцу, которая содержит максимальное количество нулей. Алгебраическое дополнение 0-го элемента можно не вычислять.

Пусть дана система уравнений вида Ах=В

Если определитель А=0, то система может решений не иметь, либо иметь бесконечное множество решений.

Если определитель А?0, то корни системы могут быть найдены следующим образом.

Пусть Ак-матрица, полученная из матрицы А путем замен к-го столбца на матрицу-столбец В. Тогда решение .

МЕТОД ОБРАТНОЙ МАТРИЦЫ

Пусть дана система Ах=В и detA?0.

Умножим обе части системы на А-1:

А-1*Ах=А-1*В>х=А-1

Способы нахождения обратной матрицы:

1. Способ основан на методе Гаусса.

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

2. Через алгебраические дополнения.

Составить матрицу алгебраических дополнений, в которой на месте aij элементов будут находиться Aij.

Разделить каждый элемент матрицы алгебраических дополнений на detA.

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


Подобные документы

  • Особенности решения алгебраических, нелинейных, трансцендентных уравнений. Метод половинного деления (дихотомия). Метод касательных (Ньютона), метод секущих. Численные методы вычисления определённых интегралов. Решение различными методами прямоугольников.

    курсовая работа [473,4 K], добавлен 15.02.2010

  • Математическая формулировка задачи, существующие численные методы и схемы алгоритмов. Интерполирование функции, заданной в узлах, методом Вандермонда. Среднеквадратичное приближение функции. Вычисление интеграла функций по составной формуле трапеций.

    курсовая работа [3,4 M], добавлен 14.04.2009

  • Изучение способов работы с файлами с помощью автоматического преобразования данных. Решение иррациональных уравнений методами хорд и половинного деления. Вычисление определенного интеграла. Решение систем линейных алгебраических уравнений. Ряды Фурье.

    курсовая работа [759,3 K], добавлен 16.08.2012

  • Методы хорд и итераций, правило Ньютона. Интерполяционные формулы Лагранжа, Ньютона и Эрмита. Точечное квадратичное аппроксимирование функции. Численное дифференцирование и интегрирование. Численное решение обыкновенных дифференциальных уравнений.

    курс лекций [871,5 K], добавлен 11.02.2012

  • Метод Гаусса, метод прогонки, нелинейное уравнение. Метод вращения Якоби. Интерполяционный многочлен Лагранжа и Ньютона. Метод наименьших квадратов, интерполяция сплайнами. Дифференцирование многочленами, метод Монте-Карло и Рунге-Кутты, краевая задача.

    курсовая работа [4,8 M], добавлен 23.05.2013

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

    контрольная работа [58,6 K], добавлен 20.11.2010

  • Решение систем линейных алгебраических уравнений методом простой итерации. Полиномиальная интерполяция функции методом Ньютона с разделенными разностями. Среднеквадратическое приближение функции. Численное интегрирование функций методом Гаусса.

    курсовая работа [2,4 M], добавлен 14.04.2009

  • Общая постановка задачи. Отделение корня. Уточнение корня. Метод половинного деления (бисекции). Метод хорд (секущих). Метод касательных (Ньютона). Комбинированный метод хорд и касательных. Задания для расчётных работ.

    творческая работа [157,4 K], добавлен 18.07.2007

  • Основные понятия теории погрешностей. Приближенное решение некоторых алгебраических трансцендентных уравнений. Приближенное решение систем линейных уравнений. Интерполирование функций и вычисление определенных интегралов, дифференциальных уравнений.

    методичка [899,4 K], добавлен 01.12.2009

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

    методичка [690,6 K], добавлен 26.01.2015

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.