Министерство образования и науки РФ
ГОУ ВПО “УГТУ-УПИ”
Курсовая работа
по “Вычислительной математике”
на тему: “Некоторые дополнительные вычислительные методы”
Семестр № 3
Преподаватель
Кочнев В.П.
Студент гр. № р-23021д
Логиновских М.А.
Номер зачетной книжки 17309013
Екатеринбург
2004
____________________________________________________________________________
_
Домашнее задание по ________________________________ № ________________
№ записи в книге регистрации __________________ дата регистрации
___________200_г.
Преподаватель _________________________________________
Студент _________________________________________ группа № ________________
Деканат ФДО _______________
СОДЕРЖАНИЕ
1. Решение систем линейных уравнений …………………………………………………… 3 а) Схема Халецкого ……………………………………………………………………....... 3 б) Метод Зейделя и условия сходимости ………………………………………………… 5
2. Методы решения нелинейных уравнений ……………………………………………….. 6 а) Метод хорд ………………………………………………………………………………. 7 б) Метод Ньютона (метод касательных) …………………………………………………. 8 в) Метод итерации ………………………………………………………………………… 9
3. Интерполирование и экстраполирование ……………………………………………….. 11 а) Интерполирование с помощью многочленов ………………………………………… 11 б) Интерполяционный многочлен Лагранжа ……………………………………………. 12 в) Интерполяционные многочлены Стирлинга и Бесселя ……………………………… 13 г) Тригонометрическое интерполирование …………………………………….………... 15
д) Интерполяция сплайнами ……………………………………………………..………... 15
4. Численное дифференцирование и интегрирование ……………………………….…….. 16 а) Постановка задачи численного интегрирования ……………………………………... 16 б) Составные квадратурные формулы ………………………………………………….… 17
5. Приближенное вычисление обыкновенных дифференциальных уравнений …………..
18 а) Метод Рунге-Кутта ……………………………………………………………………… 18 б) Экстраполяционные методы Адамса ………………………………………………….. 20 в) Метод Милна ……………………………………………………………………………. 20 г) Краевые задачи для обыкновенных дифференциальных уравнений ………………...
21
6. Приближенные методы решения дифференциальных уравнений с частными
производными ………………………………………………………………………………… 21 а) Классификация дифференциальных уравнений второго порядка …………………… 22 б) Постановка краевых задач ……………………………………………………………… 23 в) Метод конечных разностей (метод сеток) …………………………………………….. 24 г) Разностные схемы для решения уравнения теплопроводности ……………………… 25
д) Разностные схемы для решения уравнения колебания струны ……………………… 26
7. Список литературы ………………………………………………………………………… 27
1. Решение систем линейных уравнений
Системы линейных уравнений (СЛУ) имеют в вычислениях очень большое
значение, так как к ним может быть приведено приближенное решение широкого
круга задач. Так, основными источниками возникновения СЛУ являются теория
электрических цепей, уравнения балансов и сохранения в механике, гидравлике
и т.д. Существует несколько способов решения таких систем, которые в
основном делятся на два типа: 1) точные методы, представляющие собой
конечные алгоритмы для вычисления корней системы, 2) итерационные методы,
позволяющие получать корни системы с заданной точностью путем сходящихся
бесконечных процессов. Заметим, что даже результаты точных методов являются
приближенными из-за неизбежных округлений. Для итерационных процессов также
добавляется погрешность метода.
Пример системы линейных уравнений: [pic]
Или в матричном виде: [pic],
где [pic]матрица коэффициентов системы;
[pic] - вектор неизвестных; [pic] - вектор свободных членов.
Схема Халецкого
Запишем систему линейных уравнений в матричном виде: [pic],
где A=[aij] – квадратная матрица порядка n и
[pic], [pic] - векторы-столбцы.
Представим матрицу A в виде произведения нижней треугольной матрицы B=[bij]
и верхней треугольной матрицы C=[cij] с единичной диагональю [pic], где
[pic] и [pic].
Тогда элементы bij и cij определяются по формулам
[pic] и [pic]
Отсюда искомый вектор x может быть вычислен из уравнений [pic] и [pic].
Так как матрицы B и C – треугольные, то системы легко решаются:
[pic] и [pic]
Из этих двух формул видно, что числа yi выгодно вычислять вместе с
коэффициентами cij. Этот метод получил название схемы Халецкого. В схеме
применяется обычный контроль с помощью сумм. Если матрица A –
симметрическая aij=aji, то [pic]
Пример. Решить систему [pic]
Решение.
В первый раздел таблицы впишем матрицу коэффициентов системы, ее свободные
члены и контрольные суммы. Далее так как [pic] [pic], то первый столбец из
раздела 1 переносится в первый столбец раздела II. Чтобы получить первую
строку раздела II, делим все элементы первой строки раздела I на
элемент[pic], в нашем случае на 3.
Имеем: [pic]; [pic]; [pic]; [pic]; [pic].
Переходим к заполнению второго столбца раздела II, начиная со второй
строки. Пользуясь формулами, определяем [pic]: [pic]; [pic]; [pic].
Далее определяя по формулам, заполняем вторую сетку для раздела II:
[pic][pic]
[pic]
[pic]
Затем переходим к третьему столбцу, вычисляя его элементы [pic] и [pic] по
формулам и т.д., пока не будет заполнена вся таблица раздела II. Таким
образом, заполнение раздела II происходит способом “елочки”: столбец -
строка, столбец - строка и т.д.
В разделе Ш, пользуясь формулами, определяем [pic] и [pic].
Текущий контроль осуществляется с помощью столбца S, над которым
производятся те же действия, что и над столбцом свободных членов.
|0 |1,200|0,000|0,000|
| |0 |0 |0 |
|1 |1,200|1 |0,948|
| |0 |,0600|0 |
|2 |0,999|1,005|0,999|
| |2 |4 |1 |
|3 |0,999|1.000|1,000|
| |6 |1 |1 |
|4 |1 |1,000|1,000|
| |,0000|0 |0 |
|5 |1 |1,000|1,000|
| |,0000|0 |0 |
Точные значения корней: [pic].
2. Методы решения нелинейных уравнений
Как известно, далеко не всякое уравнение f(x)=0 можно решить точно, т.е. не всегда можно найти число [pic] такое что f([pic])?0. В первую очередь это относится к трансцендентным уравнениям. Кроме того, даже для алгебраических уравнений степени выше четвертой не существуют формулы, выражающей их решения через коэффициенты уравнения при помощи арифметических операций и извлечение корней. Для уравнений третьей и четвертой степени формулы для отыскания корней существуют, но они настолько сложны, что практически не применяются. Поэтому большое значение имеет приближенное вычисление корней уравнения f(x)=0. Для этого существует множество методов некоторые, из которых мы рассмотрим.
Метод хорд
Пусть дано уравнение f(x)=0, где функция f(x) определена и непрерывна на
интервале
[a, b] и f(a)f(b)