
- •Е.А. Михайлова, в.А. Михайлова Лабораторный практикум по курсу «вычислительная математика»
- •Предисловие
- •Лабораторная работа № 1 нахождение корня нелинейного уравнения
- •1. Постановка задачи
- •2. Методы решения задачи
- •2.1. Метод деления отpезка пополам
- •2.2. Метод простой итерации.
- •2.3. Метод Ньютона
- •Задание
- •Варианты заданий
- •2. Методы решения системы нелинейных уравнений
- •2.1.Метод простой итерации
- •2.2. Метод Ньютона
- •2. Формулы Ньютона – Котесса
- •2.1. Формулы метода прямоугольников
- •2.2. Формула трапеций
- •2.3. Формула Симпсона (формула парабол)
- •Несобственные интегралы 2-го рода
- •3. Пример
- •2. Методы решения
- •Задание
- •Варианты заданий
- •2. Решить уравнения:
- •Лабораторная работа № 6 краевая задача. Метод стрельбы
- •1. Постановка задачи
- •2. Метод стрельбы
- •Задание
- •Варианты заданий
- •Лабораторная работа № 7 метод наименьших квадратов
- •1. Постановка задачи
- •2. Метод наименьших квадратов
- •3. Пример
- •Задание
- •Варианты заданий
- •Лабораторная работа № 8 системы линейных уравнений
- •1. Постановка задачи
- •Метод прогонки
- •Задание
- •Список рекомендуемой литературы
- •Содержание
- •Лабораторный практикум по курсу «вычислительная математика»
- •4 00062, Г. Волгоград, просп. Университетский, 100.
Несобственные интегралы 2-го рода
Рассмотрим интеграл с конечными пределами
|
|
(4) |
где
подынтегральная функция f(x)
обращается в бесконечность в некоторой
точке
.
По определению можно представить
интеграл (4) в виде следующей суммы
|
|
|
Поэтому, чтобы вычислить (4) с заданной точностью , выбираем из условия
|
|
|
Затем по каким-либо квадратурным формулам приближенно вычисляем определенные интегралы
и
с погрешностью /4 каждый. Исходный интеграл будет равен сумме I1 и I1.
Пример. Вычислим с погрешностью =0,05 интеграл
|
|
(5) |
Подынтегральная функция имеет разрыв при с=2. Представим (5) в виде суммы двух интегралов
|
|
|
|
|
|
и
выбираем
так, чтобы
.
Так, при
интеграл I2
удовлетворяет условию
|
|
|
Следовательно, при =0.1 I2<0.028.
Учитывая полученную оценку интеграла I2, вычислим интеграл
|
|
|
по квадратурной формуле с точностью до 0,022.
Во
многих случаях приближенное вычисление
интеграла (4) облегчается с помощью
метода
выделения особенностей,
предложенного Л.В. Канторовичем. Идея
этого метода состоит в следующем. Из
подынтегральной функции f(x)
выделяют некоторую функцию g(x),
имеющую те же особенности, что и функция
f(x),
но элементарно интегрируемую на данном
промежутке
и такую, чтобы
,
т.е. функция
не должна иметь особенности в точке
х=с.
Запишем (4) в виде
|
|
(6) |
Первый интеграл в (6) берется аналитически, а второй вычисляется численно по квадратным формулам.
Подбор функции g(x) производится различным образом в зависимости от конкретного вида f(x). Рассмотрим правило построения такой функции для f(x), имеющей вид
|
|
|
где
- непрерывная функция вместе со своими
N
производными. Разложим функцию F(x)
в степенной ряд в окрестности точки
х=с:
|
|
|
и удержим столько слагаемых, чтобы выполнялось условие
|
|
(7) |
Тогда интеграл (6) можно записать в виде
|
|
(8) |
причем первый интеграл в (8) берется аналитически
|
|
|
|
|
|
а
второй, в силу условия (7), можно найти
численно, применив одну из известных
квадратных формул. Это возможно, так
как подынтегральная функция
не имеет особенности в точке х=с
и непрерывна вместе со своими N
производными.