- •Оглавление
- •Предисловие
- •Лекция 1. Понятие погрешности решения систем алгебраических уравнений
- •1.1. Точные и приближенные числа, источники погрешностей, классификация погрешностей
- •Источники погрешностей
- •1.2. Абсолютные и относительные погрешности
- •1.3. Понятие о системе линейных уравнений
- •2.1. Формулы Крамера для решения системы линейных уравнений
- •2.2. Решение системы линейных уравнений методом последовательного исключения неизвестных (метод Гаусса)
- •Пусть дана система
- •1.4. Матричные уравнения
- •Лекция 2. Точные методы решения систем линейных алгебраических уравнений (слау)
- •Обозначим
- •Лекция 3. Приближенные методы решения слау
- •3.1. Понятие предела для векторов и матриц Пусть дана последовательность векторов
- •Аналогично, если имеется последовательность квадратных матриц
- •3.2. Метод простой итерации Пусть дана система линейных уравнений:
- •В матричном виде:
- •Условия сходимости итерационного процесса
- •3.4. Оценка погрешности приближенного процесса методом итерации
- •Лекция 4. Метод зейделя
- •4.1. Условия сходимости процесса Зейделя
- •Пусть дана линейная система
- •4.2. Оценка погрешности процесса Зейделя
- •4.3. Приведение системы линейных уравнений к виду, удобному для итераций
- •Лекция 5. Методы решения нелинейных уравнений
- •5.1. Итерационные методы решения нелинейных уравнений а). Метод простой итерации.
- •Число итераций при использовании этого метода
- •Если корень нас не устраивает, то мы находим
- •Г). Метод Ньютона.
- •5.2. Метод Рыбакова
- •5.3. Метод наискорейшего спуска
- •Лекция 6. Интерполирование и экстраполирование
- •6.1. Интерполяционный многочлен Лагранжа
- •Перепишем этот многочлен в другой форме
- •Окончательно получаем
- •Лекция 7. Конечные разности
- •Рассмотрим некоторые свойства конечных разностей
- •7.1. Первая интерполяционная формула Ньютона для равноотстоящих узлов интерполяции
- •7.2. Вторая интерполяционная формула Ньютона
- •7.3. Оценка погрешностей интерполяционных формул Ньютона
- •Лекция 8. Линейное интерполирование
- •8.1. Интерполирование по Эйткину
- •8.2. Интерполяция и приближение сплайном
- •Приближение линейными сплайнами
- •Лекция 9. Численное интегрирование функций
- •9.1. Постановка задачи
- •9.2. Формула прямоугольников
- •Лекция 10. Численные методы решения задачи коши для обыкновенных дифференциальных уравнений
- •10.1. Задача Коши. Общие замечания. Постановка задачи
- •10.2. Метод Эйлера
- •10.4. Метод Рунге-Кутта
- •Лекция 11. Разностный метод решения уравнений математической физики (метод сеток)
- •11.1. Метод сеток для уравнения параболического типа
- •Лекция 12. Метод сеток для уравнения гиперболического типа
- •12.2. Метод сеток для уравнений Пуассона и Лапласа
- •Образцы выполнения контрольных работ Контрольная работа № 5
- •Контрольная работа № 6
- •Решение уравнения методом Эйлера
- •Решение уравнения модифицированным методом Эйлера
- •Решение уравнения методом Рунге-Кутта
- •Аналитическое решение заданного уравнения
- •Сравнение точного решения и приближенных решений исходного дифференциального уравнения
- •Контрольная работа №7
- •Контрольная работа №8
- •Решение дифференциального уравнения методом Рунге-Кутта
- •Аналитические решения заданного уравнения
- •Примеры заданий для самостоятельной работы Задание 1
- •Задание 2
- •Задание 3
- •Задание 4
- •Задание 5
- •Варианты задания
- •Задание 6
- •Варианты задания
- •Задание 7
- •Варианты задания
- •Задание 8
- •Варианты задания
- •Задача для самостоятельного решения
- •Вопросы для самостоятельной подготовки
- •Библиографический список
Аналитические решения заданного уравнения
Дано
у'' – 3у' + 2у =х, у(0) = у'(0) = 1. (***)
Это линейное неоднородное дифференциальное уравнение 2-го порядка с постоянными коэффициентами. Найдем общее решение уон этого уравнения, которое равно сумме какого-либо частного решения учн неоднородного уравнения и общего решения у∞ соответствующего однородного уравнения
уон = у∞ + учн.
Запишем соответствующее однородное уравнение
у'' – 3у' + 2у = 0,
его характеристическое уравнение k2 – 3k + 2 = 0 имеет корни k1 = 2, k2 = 1. Поэтому
у∞ = с1е2х + с2е2х.
Для определения частного решения неоднородного уравнения (***) воспользуемся методом неопределенных коэффициентов. Так как правая часть уравнения (***) является многочленом первой степени и число ά + βi = 0 не является корнем характеристического уравнения, то частное решение имеет вид:
учн =Ах +В.
Чтобы определить значения коэффициентов А и В, находим производные
учн = А, учн = 0,
подставляем учн, учн, учн в уравнение (19)
-3А +2Ах + 2В = х
и, сравнивая коэффициенты при одинаковых степенях х, получаем систему
из
которой находим
.
Следовательно,
и
Найдем частное решение, удовлетворяющее заданным начальным условиям, для чего продифференцируем общее решение и подставим начальные условия в общее решение и его производную, в результате получим систему
решая
которую, найдем
.Таким
образом, получили аналитическое решение
заданного уравнения, удовлетворяющего
данным начальным условиям
3. Сравним значения точного и приближенного решений заданного уравнения в точках х1, х2, х3. Это сравнение дано в табл.12.
Таблица 12
i |
yi по методу Рунге-Кутта |
у(хi) – точное решение |
|
0 |
1,0000000 |
1,000000 |
|
1 |
1,105349 |
1,105351 |
2∙10-6 |
2 |
1,222951 |
1,222955 |
4∙10-6 |
3 |
1,3555196 |
1,3555301 |
1∙10-5 |
Примеры заданий для самостоятельной работы Задание 1
Решить приближенно методом итерации и методом Зейделя. Сравнить ответы.
Варианты задания:
1.
|
14.
|
2.
|
15.
|
3.
|
16.
|
4.
|
17.
|
5.
|
18.
|
6.
|
19.
|
7.
|
20.
|
8.
|
21.
|
9.
|
22.
|
10.
|
23.
|
11.
|
24.
|
12.
|
25.
|
13.
|
|
Примечание. Для лучшей сходимости итерационного процесса необходимо преобразовать уравнения, выделив наибольшие диагональные элементы при неизвестных. Для этого необходимо или переобозначить неизвестные (циклическая перестановка) или переставить уравнения.
