- •Гоу впо «Кемеровский государственный университет»
- •Методы вычислений
- •Предисловие
- •Тема 1. Элементы теории погрешностей
- •1.1. Определения
- •1.2. Абсолютная и относительная погрешности
- •1.3. Значащие цифры и число верных знаков
- •1.4. Погрешности арифметических действий
- •1.5. Погрешность вычисления функции
- •Задания
- •Варианты
- •Тема 2. Интерполирование
- •2.1. Постановка задачи интерполирования
- •2.2. Интерполяционный многочлен Лагранжа
- •2.3. Интерполяционная формула Гаусса
- •2.4. Сплайн- интерполяция
- •2.5. Линейный сплайн
- •Таким образом, линейный сплайн имеет вид
- •2.6. Параболический сплайн
- •2.7. Кубический сплайн
- •Задания
- •Варианты функций
- •Тема 3. Численное интегрирование
- •3.1. Постановка задачи интегрирования
- •3.2. Квадратурные формулы
- •3.3. Выбор шага интегрирования
- •3.4. Квадратурная формула Гаусса
- •Тогда формула Гаусса будет иметь вид
- •Задания
- •Варианты
- •Тема 4. Решение трансцендентных (нелинейных) уравнений
- •4.1. Отделение корней
- •4.2. Метод последовательных приближений (метод простой итерации)
- •4.3. Метод половинного деления (метод проб, метод дихотомии)
- •4.4. Метод пропорционального деления (метод хорд)
- •4.5. Метод Ньютона (метод касательных)
- •4.6. Метод Ньютона модифицированный
- •4.7. Метод Чебышева
- •Задания
- •Варианты уравнений
- •Тема 5. Решение спектральной задачи
- •5.1. Метод скалярных произведений
- •5.2. Метод вращения
- •Задания
- •Варианты матриц
- •Тема 6. Решение систем линейных алгебраических уравнений
- •6.1. Обусловленность матрицы
- •6.2. Метод Гаусса
- •6.3. Метод Гаусса с выбором главного элемента
- •6.4. Нахождение определителя и обращение матрицы с помощью метода Гаусса
- •6.5. Итерационные методы (метод Якоби, метод Зейделя, метод релаксации)
- •6.6. Оптимизация скорости сходимости итерационного процесса
- •6.7. Итерационные методы вариационного типа
- •6.8. Методы сопряженных направлений
- •Задания
- •Тема 7. Решение систем нелинейных алгебраических уравнений
- •7.1. Определения
- •Интеграл в правой части (8.4) вычисляется с помощью численной квадратуры
- •8.2. Метод Эйлера
- •8.3. Методы Рунге-Кутта
- •8.4. Метод Адамса
- •8.5. Метод Милна
- •Задания
- •Варианты
- •8.6. Краевая задача для обыкновенного дифференциального уравнения второго порядка
- •8.6.1. Интегро-интерполяционный метод
- •8.6.2. Метод прогонки
- •Задания
- •Варианты заданий
- •Тема 9. Методы решения дифференциальных уравнений в частных производных
- •9.1. Простейшие приемы построения разностных схем
- •9.2. Сходимость, аппроксимация, устойчивость разностных схем
- •9.3. Решение уравнения параболического типа
- •9.3.1. Явная разностная схема
- •9.3.2. Неявная разностная схема
- •9.3.3. Реализация метода разностной прогонки для уравнения параболического типа
- •Задание
- •Варианты
- •9.4. Решение уравнения эллиптического типа
- •9.4.1. Метод матричной прогонки
- •Задание:
- •Варианты
- •Содержание
- •Тема 7. Решение систем нелинейных алгебраических уравнений 40
- •Тема 8. Решение обыкновенных дифференциальных уравнений 43
- •Тема 9. Методы решения дифференциальных уравнений в частных производных 49
3.4. Квадратурная формула Гаусса
В формулах Гаусса
коэффициенты Ai абсциссы ti подбираются так, чтобы формула была точной для всех многочленов наивысшей возможной степени 2n-1 (табл. 3.1).
При вычислении интеграла следует сделать замену
xi= .
Тогда формула Гаусса будет иметь вид
= ,
где .
Таблица 3. 1
N |
ti |
Ai |
Rn |
4 |
-t1=t4=0,86113631 -t2=t3=0,33998104 |
A1=A4=0,34785484 A2=A3=0,65214516 |
R42,88*10-7f(8)(ξ), -1< ξ<1 |
5 |
-t1=t5=0,90617985 -t2=t4=0,53846931 t3=0 |
A1=A5=0,23692688 A2=A4=0,47862868 A3=0,56888889 |
R50,08*10-8f(10)(ξ), -1< ξ<1 |
7 |
-t1=t7=0,94910791 -t2=t6=0,74153119 -t3=t5=0,40584515 t4=0 |
A1=A7=0,12948496 A2=A6=0,27970540 A3=A5=0,38183006 A4=0,41795918 |
R72,17*10-15f(14)(ξ), -1< ξ<1 |
Задания
Вычислить приближенное значение интеграла , используя формулы трапеции, Симпсона, трех-восьмых, прямоугольников и Гаусса (n=4, 5 или 7). Оценить остаточный член формул.
Вычислить значение интеграла с заданной точностью , используя формулу трапеции или Симпсона, двумя способами:
- выбрать шаг интегрирования из оценки остаточного члена,
- использовать метод последовательного удвоения числа шагов.
Вычислить значение интеграла с заданной точностью , если функция f(x) имеет разрыв второго рода.
Вычислить значение интеграла , заменив функцию f(x) кубическим сплайном.
Вычислить двойной интеграл .
Варианты
Для заданий 1,2,4:
f(x)=x3 e2x ; a=0 ; b=1.
f(x)= ; a=0 ; b=4.
f(x)= ; a= -2 ; b= -1.
f(x)= ; a=0 ; b=1.
f(x)= ; a=0 ; b=1.
f(x)= ; a=1 ; b=3.
f(x)= ; a=1 ; b=2.
f(x)= ; a=0,5 ; b=2,5.
f(x)= ; a=5 ; b=7.
f(x)= ; a=0 ; b=5.
f(x)=cos(x)/(x+2) ; a=0,4 ; b=1,2.
f(x)= ; a=0,4 ; b=1.2.
f(x)=(x+1)sin(x) ; a=1,6 ; b=2.4.
f(x)=(x+1)cos(x2) ; a=0,2 ; b=1.
f(x)=sin(x2-0,4)/(x+2) ; a=0,8 ; b=1,2.
f(x)=ln(1+x2)/(1+x2) ; a=0; b=1.
f(x)=ln(5+4cos(x)) ; a=0; b=3,1416.
f(x)=x*ln(1+x); a=0; b=1.
Для задания 3:
f(x)= ; a= 0; b= 0,5.
f(x)= ; a= 0; b= 1.
f(x)= ; a= 0; b= 1.
f(x)= ; a= 0; b= 2.
f(x)= ; a= 0; b= 1.
f(x)= ; a= 1; b= 2.
f(x)= ; a= -1; b= 1.
f(x)= ; a= -1; b= 1.
f(x)= ; a= -1; b= 1.
f(x)= ; a= 0; b= 1.
Для задания 5:
.
.
.
.
.
.
правильный шестиугольник, вписанный в единичный круг.
.
.
ромб с центром в начале координат и с вершинами в точках (0, -4); (0, 4); (-2, 0); (2, 0).
Тема 4. Решение трансцендентных (нелинейных) уравнений
Пусть задана функция f(x) действительного переменного. Требуется найти корни уравнения
f(x)=0. (4.1)
Задача нахождения корней уравнения (4.1) обычно решается в два этапа. На первом этапе изучается расположение корней и проводится их разделение, то есть выделяются области, содержащие только один корень. На втором этапе, используя начальное приближение, строится итерационный процесс для уточнений корня.