kursovik (675853)
Текст из файла
Министерство образования и науки РФ
ГОУ ВПО “УГТУ-УПИ”
Курсовая работа
по “Вычислительной математике”
на тему: “Некоторые дополнительные вычислительные методы”
Семестр № 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) итерационные методы, позволяющие получать корни системы с заданной точностью путем сходящихся бесконечных процессов. Заметим, что даже результаты точных методов являются приближенными из-за неизбежных округлений. Для итерационных процессов также добавляется погрешность метода.
Пример системы линейных уравнений:
где матрица коэффициентов системы;
- вектор неизвестных;
- вектор свободных членов.
Схема Халецкого
Запишем систему линейных уравнений в матричном виде: ,
где A=[aij] – квадратная матрица порядка n и
Представим матрицу A в виде произведения нижней треугольной матрицы B=[bij] и верхней треугольной матрицы C=[cij] с единичной диагональю , где
Тогда элементы bij и cij определяются по формулам
Отсюда искомый вектор x может быть вычислен из уравнений и
.
Так как матрицы B и C – треугольные, то системы легко решаются:
Из этих двух формул видно, что числа yi выгодно вычислять вместе с коэффициентами cij. Этот метод получил название схемы Халецкого. В схеме применяется обычный контроль с помощью сумм. Если матрица A – симметрическая aij=aji, то
Решение.
В первый раздел таблицы впишем матрицу коэффициентов системы, ее свободные члены и контрольные суммы. Далее так как
, то первый столбец из раздела 1 переносится в первый столбец раздела II. Чтобы получить первую строку раздела II, делим все элементы первой строки раздела I на элемент
, в нашем случае на 3.
Переходим к заполнению второго столбца раздела II, начиная со второй строки. Пользуясь формулами, определяем :
;
;
.
Далее определяя по формулам, заполняем вторую сетку для раздела II:
Затем переходим к третьему столбцу, вычисляя его элементы и
по формулам и т.д., пока не будет заполнена вся таблица раздела II. Таким образом, заполнение раздела II происходит способом “елочки”: столбец - строка, столбец - строка и т.д.
В разделе Ш, пользуясь формулами, определяем и
.
Текущий контроль осуществляется с помощью столбца ∑, над которым производятся те же действия, что и над столбцом свободных членов.
|
|
|
|
|
|
|
|
|
| |||
I |
|
|
|
|
|
| 3 | 1 | -1 | 2 | 6 | 11 |
I |
|
|
|
|
|
| -5 | 1 | 3 | -4 | -12 | -17 |
I |
|
|
|
|
|
| 2 | 0 | 1 | 1 | 1 | 3 |
I |
|
|
|
|
|
| 1 | -5 | 3 | 3 | 3 | -1 |
II |
|
|
|
|
|
| 3 | 0.333333 | -0.333333 | 2 | 2 | 3.666667 |
II |
|
|
|
|
|
| -5 | 2.666667 | -0.25 | 0.25 | -0.75 | 0.5 |
II |
|
|
|
|
|
| 2 | -0.666667 | 2 | -1.25 | -1.75 | -2 |
II |
|
|
|
|
|
| 1 | -5.333333 | 6 | 2.5 | 3 | 4 |
III |
|
| 2 | 1 | ||||||||
III |
|
| -0.75 | -1 | ||||||||
III |
|
| -1.75 | 2 | ||||||||
III |
|
| 3 | 3 |
Метод Зейделя и условия сходимости
Этот метод представляет собой модификацию метода простой итерации. Его смысл заключается в том, что при вычислении (k+1)-го приближения неизвестной xi учитываются уже вычисленные ранее (k+1)-е приближения x1, x2, ..., xi-1. Пусть дана приведенная линейная система (i = 1, 2, …n). Выберем произвольно начальные приближения корней
, стараясь, чтобы они в какой-то мере соответствовали искомым неизвестным x1, x2, x3, ..., xn. Предположим, что k-е приближение
корней известно, тогда в соответствии с идеей метода будем строить (k+1)–е приближение по следующим формулам:
Обычно процесс Зейделя сходится быстрее, чем метод простой итерации. Бывает, что процесс Зейделя сходится, когда простая итерация расходится и т.п. Правда, бывает и наоборот. Во всяком случае, достаточные условия сходимости для метода простой итерации достаточны и для сходимости метода Зейделя. То есть процесс итерации сходится, если выполнено одно из условий
Характеристики
Тип файла документ
Документы такого типа открываются такими программами, как Microsoft Office Word на компьютерах Windows, Apple Pages на компьютерах Mac, Open Office - бесплатная альтернатива на различных платформах, в том числе Linux. Наиболее простым и современным решением будут Google документы, так как открываются онлайн без скачивания прямо в браузере на любой платформе. Существуют российские качественные аналоги, например от Яндекса.
Будьте внимательны на мобильных устройствах, так как там используются упрощённый функционал даже в официальном приложении от Microsoft, поэтому для просмотра скачивайте PDF-версию. А если нужно редактировать файл, то используйте оригинальный файл.
Файлы такого типа обычно разбиты на страницы, а текст может быть форматированным (жирный, курсив, выбор шрифта, таблицы и т.п.), а также в него можно добавлять изображения. Формат идеально подходит для рефератов, докладов и РПЗ курсовых проектов, которые необходимо распечатать. Кстати перед печатью также сохраняйте файл в PDF, так как принтер может начудить со шрифтами.