Решение линейных интегральных уравнений первого рода

I. ЦЕЛЬ РАБОТЫ

Приобретение навыков решения интегральных уравнений первого рода

II. ТЕОРЕТИЧЕСКИЕ СВЕДЕНИЯ

Задачи для интегральных уравнений первого рода являются некорректно поставленными.

Рассмотрим уравнение Фредгольма первого рода

решение линейных интегральных уравнений первого рода - student2.ru , решение линейных интегральных уравнений первого рода - student2.ru , (1)

где ядро решение линейных интегральных уравнений первого рода - student2.ru является непрерывной функцией по переменным решение линейных интегральных уравнений первого рода - student2.ru ; решение линейных интегральных уравнений первого рода - student2.ru - известная функция; решение линейных интегральных уравнений первого рода - student2.ru - искомая функция.

Будем полагать, что уравнение (1) с точной правой частью решение линейных интегральных уравнений первого рода - student2.ru имеет единственное решение.

Если вместо решение линейных интегральных уравнений первого рода - student2.ru известно лишь ее приближение решение линейных интегральных уравнений первого рода - student2.ru , мало отличающееся (в метрике решение линейных интегральных уравнений первого рода - student2.ru ) от решение линейных интегральных уравнений первого рода - student2.ru , то можем искать лишь приближенное решение уравнения (1). В качестве приближенного решения уравнения (1) нельзя брать точное решение уравнения

решение линейных интегральных уравнений первого рода - student2.ru . (2)

так как такого решения может не существовать. Кроме того, такое решение не обладает свойством устойчивости к малым изменениям правой части уравнения.

Будем искать решение уравнения (2) методом регуляризации. Согласно методу построим сглаживающий функционал решение линейных интегральных уравнений первого рода - student2.ru , выбрав стабилизатор первого порядка решение линейных интегральных уравнений первого рода - student2.ru :

решение линейных интегральных уравнений первого рода - student2.ru решение линейных интегральных уравнений первого рода - student2.ru , (3)

где решение линейных интегральных уравнений первого рода - student2.ru - параметр регуляризации;

решение линейных интегральных уравнений первого рода - student2.ru

решение линейных интегральных уравнений первого рода - student2.ru - заданные неотрицательные непрерывные функции (если нет специальных соображений по выбору функций решение линейных интегральных уравнений первого рода - student2.ru , то обычно полагают решение линейных интегральных уравнений первого рода - student2.ru ). Затем найдем функцию решение линейных интегральных уравнений первого рода - student2.ru , минимизирующую функционал (3), причем параметр решение линейных интегральных уравнений первого рода - student2.ru определим по невязке, т.е. из условия

решение линейных интегральных уравнений первого рода - student2.ru , (4)

где решение линейных интегральных уравнений первого рода - student2.ru - уклонение правой части интегрального уравнения в метрике пространства решение линейных интегральных уравнений первого рода - student2.ru , которое считаем известным:

решение линейных интегральных уравнений первого рода - student2.ru

Решение решение линейных интегральных уравнений первого рода - student2.ru будет устойчиво к малым изменениям в метрике решение линейных интегральных уравнений первого рода - student2.ru правой части уравнения решение линейных интегральных уравнений первого рода - student2.ru .

Функция решение линейных интегральных уравнений первого рода - student2.ru будет являться приближенным решением уравнения (1).

Для нахождения параметра регуляризации будем проводить расчеты с несколькими значениями параметра решение линейных интегральных уравнений первого рода - student2.ru

(например, решение линейных интегральных уравнений первого рода - student2.ru

Для каждого значения решение линейных интегральных уравнений первого рода - student2.ru находим функцию решение линейных интегральных уравнений первого рода - student2.ru , минимизирующую функционал (3).

В качестве искомого значения параметра регуляризации возьмем такое число решение линейных интегральных уравнений первого рода - student2.ru , для которого с требуемой точностью выполняется условие (4), т.е. невязка, полученная при подстановке найденной функции решение линейных интегральных уравнений первого рода - student2.ru в уравнение (2), должна быть сравнима с погрешностью правой части интегрального уравнения.

Обратимся теперь к вариационной задаче

решение линейных интегральных уравнений первого рода - student2.ru

Полагая решение линейных интегральных уравнений первого рода - student2.ru , получим

решение линейных интегральных уравнений первого рода - student2.ru (5)

Проведем дискретизацию сглаживающего функционала, воспользовавшись разностным методом.

Аппроксимируем входящие в функционал решение линейных интегральных уравнений первого рода - student2.ru интегралы квадратурными формулами. Для этого введем на прямоугольнике решение линейных интегральных уравнений первого рода - student2.ru сетку решение линейных интегральных уравнений первого рода - student2.ru решение линейных интегральных уравнений первого рода - student2.ru , так, что решение линейных интегральных уравнений первого рода - student2.ru .

Для простоты рассмотрим равномерную сетку решение линейных интегральных уравнений первого рода - student2.ru

где

решение линейных интегральных уравнений первого рода - student2.ru .

Вычислим решение линейных интегральных уравнений первого рода - student2.ru по формуле средних, заменяя производную разностным отношением

решение линейных интегральных уравнений первого рода - student2.ru ,

где решение линейных интегральных уравнений первого рода - student2.ru .

Таким образом,

решение линейных интегральных уравнений первого рода - student2.ru . (6)

Остальные интегралы вычислим по формуле трапеций

решение линейных интегральных уравнений первого рода - student2.ru (7)

где

решение линейных интегральных уравнений первого рода - student2.ru

решение линейных интегральных уравнений первого рода - student2.ru

где решение линейных интегральных уравнений первого рода - student2.ru

решение линейных интегральных уравнений первого рода - student2.ru (8)

где

решение линейных интегральных уравнений первого рода - student2.ru

Подставляя выражения (6) - (8) в (5), получим

решение линейных интегральных уравнений первого рода - student2.ru (9)

Для решения задачи (9) приравняем к нулю производные от левой части (9) по решение линейных интегральных уравнений первого рода - student2.ru решение линейных интегральных уравнений первого рода - student2.ru . Получим систему уравнений, линейных относительно решение линейных интегральных уравнений первого рода - student2.ru :

решение линейных интегральных уравнений первого рода - student2.ru (10)

где

решение линейных интегральных уравнений первого рода - student2.ru

Систему (10) решим каким-либо методом, например, методом Гаусса.

Параметр решение линейных интегральных уравнений первого рода - student2.ru следует подобрать способом, указанным выше. Заметим, что условие (4) в результате дискретизации запишем в виде:

решение линейных интегральных уравнений первого рода - student2.ru

III. ЗАДАНИЕ

Найти приближенное решение интегрального уравнения

решение линейных интегральных уравнений первого рода - student2.ru

полагая решение линейных интегральных уравнений первого рода - student2.ru

Здесь решение линейных интегральных уравнений первого рода - student2.ru - последняя цифра в номере группы; решение линейных интегральных уравнений первого рода - student2.ru - номер фамилии студента в журнале группы.

IV. ОФОРМЛЕНИЕ ОТЧЕТА

В отчете должны быть представлены:

1. Название работы.

2. Постановка задачи.

3. Описание алгоритма (метода) решения.

4. Текст программы с описанием.

5. Результаты работы программы.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. — М.: Наука, 1979. 286 с.

2. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. — М.: БИНОМ. Лаборатория Знаний, 2007 636с.

3. Калиткин Н.Н. Численные методы. -М.: Наука, 1978. 512 с.

Наши рекомендации