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

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ

САМАРСКИЙ ГОСУДАРСТВЕННЫЙ АЭРОКОСМИЧЕСКИЙ УНИВЕРСИТЕТ имени академика С. П. Королева

Кафедра высшей математики

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

К курсовой работе по уравнениям математической физики

САМАРА 2009г.

Реферат

Курсовая работа: пояснительная записка, 30 страниц,8 рисунков, 3 источника, 6 таблиц.

Ключевые слова: МЕТОД КОНЕЧНЫХ РАЗНОСТЕЙ, МЕТОД ЦЕНТРАЛЬНЫХ РАЗНОСТЕЙ, МЕТОД ПРОГОНКИ, МЕТОД ГАЛЕРКИНА, МЕТОД КОЛЛОКАЦИИ, МЕТОД РИТЦА, МЕТОД ЛИБМАНА, ПРИБЛИЖЕННОЕ РЕШЕНИЕ, МЕТОД СЕТОК, ИНТЕГРАЛЬНОЕ УРАВНЕНИЕ, УРАВНЕНИЕ ТЕПЛОПРОВОДНОСТИ, ВОЛНОВОЕ УРАВНЕНИЕ.

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

Необходимо найти приближенное решение линейного дифференциального уравнения второго порядка с помощью методов Галеркина, Ритца и коллокации. Сравнить результаты, построив графики.

Нужно с помощью метода Либмана отыскать приближенное решение задачи Дирихле в квадрате.

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

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

СОДЕРЖАНИЕ

ВВЕДЕНИЕ

I. РЕШЕНИЕ КРАЕВОЙ ЗАДАЧИ

II. МЕТОДЫ ГАЛЕРКИНА, РИТЦА И КОЛЛОКАЦИЙ

III. РЕШЕНИЕ ЗАДАЧИ ДИРИХЛЕ1

IV. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ПЕРВОЙ СМЕШАННОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ НА ОТРЕЗКЕ

V. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ПЕРВОЙ СМЕШАННОЙ ЗАДАЧИ ДЛЯ ВОЛНОВОГО УРАВНЕНИЯ НА ОТРЕЗКЕ

VI. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ИНТЕГРАЛЬНОГО УРАВНЕНИЯ

ЗАКЛЮЧЕНИЕ

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

ВВЕДЕНИЕ

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

Необходимо найти приближенное решение линейного дифференциального уравнения второго порядка с помощью методов Галеркина, Ритца и коллокации. Сравнить результаты, построив графики.

Нужно с помощью метода Либмана отыскать приближенное решение задачи Дирихле в квадрате.

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

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

I. РЕШЕНИЕ КРАЕВОЙ ЗАДАЧИ

Пусть дано дифференциальное уравнение второго порядка

, (1)

Где функция задана таблично

IFi (x)
08,1548
16,8925
25,8327
34,9907
44,3818
54,0188
63,9098
74,0581
84,4615
95,1129
106

Будем искать решение уравнения (1), удовлетворяющее следующим краевым условиям

(2)

Запишем таблицу значений функций

I
0000
10,1-0,20,03
20,2-0,40,12
30,3-0,60,27
40,4-0,80,48
50,5-10,75
60,6-1,21,08
70,7-1,41,47
80,8-1,61,92
90,9-1,82,43
101-23

1. Метод конечных разностей для линейных дифференциальных уравнений второго порядка.

1. Пусть и значения и в каждом узле можно записать конечно-разностными отношениями

(3)

Тогда, используя (3), заменим уравнения (1), (2) системой

(4)

Решая систему (4), получим

2. Пусть тогда, используя (3), заменим уравнения (1), (2) системой:

(5)

Решая систему (5), получим

2. Метод центральных разностей для линейных дифференциальных уравнений второго порядка.

1. Пусть и значения и в каждом узле можно записать центрально-разностными отношениями

(6)

Тогда, используя (6), заменим уравнения (1), (2) системой:

(7)

Решая систему (7), получим:

2. Пусть , тогда, используя (6), заменим уравнения (1), (2) системой:

(8)

Решая систему (8), получим

Рис.1-– решение, полученное с помощью метода конечных разностей (h=0,1), – решение, полученное с помощью метода центральных разностей (h=0,1), – точное решение

Рис.2-– решение, полученное с помощью метода конечных разностей (h=0,2), – решение, полученное с помощью метода центральных разностей (h=0,2) -точное решение

Рис.3- Общий график решений

3. Метод прогонки для линейных дифференциальных уравнений второго порядка.

Конечно-разностные отношения в методе прогонки.

1. Пусть и значения и в каждом узле можно записать конечно-разностными отношениями:

(9)

Тогда, используя (20), заменим уравнения (1), (2), (3) системой:

(10)

Запишем первые n-1 уравнений в виде:

, где (11)

Из системы (21) следует, что (12)

, вычисляются последовательно, но при i=0:

(13)

Остальные , вычисляются по формуле:

(14)

Прямой ход вычислений.

По формулам (11) вычисляем . Далее вычисляем по формулам (13) и по рекуррентным формулам (14) находим .

Обратный ход.

Из уравнения (12) при i=n-2 и из последнего уравнения системы (10) получаем:

Решив эту систему относительно , получим

(15)

При i=n-2,…,1 используем формулу (12)

вычисляем из второго уравнения системы (10)

(16)

В результате вычислений получим таблицу:

Таблица №1

Прямой ходОбратный ход
IXiPiQiFiMiKiCiDiYi
00008.1548-21-1.1250.0815483.049606
10.1-0.20.036.9025-2.021.0203-1.146580.1626292.744645
20.2-0.40.125.8327-2.041.0412-1.181770.2524762.521233
30.3-0.60.274.9907-2.061.0627-1.243580.3669842.361553
40.4-0.80.484.3818-2.081.0848-1.368060.5388932.250789
50.5-10.754.0188-2.11.1075-1.709770.8566772.176909
60.6-1.21.083.9098-2.121.1308-5.359131.6954012.130132
70.7-1.41.474.0581-2.141.15470.24702410.532052.10254
80.8-1.61.924.4615-2.161.1792-0.40795-3.023272.087729
90.9-1.82.435.1129-2.181.2043-0.59217-1.434182.080518
101-236-2.21.23-0.67952-0.984612.076684

2. Пусть

В результате вычислений по формулам (9)-(16) получим таблицу:

Таблица №2

Прямой ходОбратный ход
IXiPiQiFiMiKiCiDiYi
00008.1548-21-1.1250.0815482.048941
10.2-0.40.125.8327-2.041.0412-1.151210.1560741.844047
20.4-0.80.484.3818-2.081.0848-1.203130.2475191.720701
30.6-1.21.083.9098-2.121.1308-1.316650.4076221.650761
40.8-1.61.924.4615-2.161.1792-1.646360.8359651.619574
51-236-2.21.23-5.714925.9362931.63769

Рис.3-– решение, полученное с помощью метода прогонки с использованием конечно-разностных отношений (h=0,1),– решение, полученное с помощью метода прогонки с использованием конечно разностных отношений (h=0,2) , – точное решение

II. Методы Галеркина, Ритца и коллокаций

Пусть дано дифференциальное уравнение второго порядка и его граничные условия

(17)

1. Метод Галеркина

Введем операторы

На отрезке [a, b] выберем систему базисных функций

Проверим систему на ортогональность

Выбранная система базисных функций является ортогональной и удовлетворяет условию выбора конечной системы базисных функций

Решение краевой задачи (17) ищется в виде

1. Рассмотрим решение задачи (17) с двумя базисными функциями:

Тогда решение

Рассмотрим выражение

(18)

Выражение (18) называется невязкой. Для задачи (1) с двумя базисными функциями

Сi выбирается таким образом, чтобы

Так как ортогональна ко всем базисным функциям, то

Тогда решение задачи (17)

2. Рассмотрим решение задачи (17) с тремя базисными функциями

Тогда решение

Невязка примет вид

Коэффициенты с1 и с2 будем искать из системы

Тогда решение задачи (17)

2. Метод коллокации

Введем операторы

На отрезке [a, b] выберем систему базисных функций

Будем искать решение задачи (17) в виде

1. Рассмотрим решение задачи (17) с двумя базисными функциями

Тогда решение

Составим невязку

На отрезке [-π, π] выберем за точку коллокации 0.

Таким образом, решение задачи (17)

.

2. Рассмотрим решение задачи (17) с тремя базисными функциями

Тогда решение

Составим невязку

На отрезке [-π, π] выберем две точки коллокации: 0 и . Составим систему уравнений

Таким образом, решение задачи (17)

3. Метод Ритца

Составим функционал по формуле

(19)

На отрезке [a, b] выберем систему базисных функций

Будем искать решение задачи (17) в виде

Подставим в (19)

Составим систему уравнений относительно с1 , с2

Таким образом, решение задачи (17)

Рис.4- у1(х)-решение, полученное с помощью метода Галеркина (две базисные функции), у2(х)-решение, полученное с помощью метода коллокации (две базисные функции)

Рис.4-у2(х)- решение, полученное с помощью метода Галеркина (три базисные функции), у4(х)- решение, полученное с помощью метода коллокации (три базисные функции), у5(х)- решение, полученное с помощью метода Ритца (три базисные функции)

Замечание: найти решение методом Ритца для двух базисных функций не удалось, т. к. функция Ф(с1 ) не квадратична относительно переменной с1 и не удовлетворяет условию существования экстремума

III. Решение задачи Дирихле

Применяя метод сеток с шагом , найти решение задачи Дирихле в квадрате с вершинами А(0,0), В(0,1), С(1,1), D(1,0).

(20)

1. Метод Либмана

Найдем значения функции В каждом узле:

На АВ

На ВС

На СD

На АD

Запишем формулу метода последовательных приближений

Пусть , тогда получим

Таблица №3

IU1,1U1,2U2,1U2,2
00000
12,511,49527,56,4952
27,248813,7449,74888,744
38,373215,493411,498210,4934
49,247916,2118512,2166511,21185
59,60712516,6101412,6149411,61014
69,80626916,7995212,8043211,79952
79,90095816,8966512,9014511,89665

2. Метод Гаусса

Для нахождения точного решения задачи (20) используем метод Гаусса. Для этого решим систему

Линейный дифференциальный уравнение

(20*)

Введем замену

Тогда (20*) перепишем в виде

Решая систему, получим

Таким образом, получим точное решение задачи (20)

IV. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ПЕРВОЙ СМЕШАННОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ НА ОТРЕЗКЕ

Пусть дано уравнение теплопроводности и его граничные условия

(21)

Решим задачу (21), применяя метод сеток для уравнений параболического типа.

1. Пусть , тогда l=0,02- шаг по оси t, а h=0,2- шаг по оси x. Решение будем искать в виде

(22)

Где (23)

Получим таблицу:

Таблица №4

JTj /xi00,20,40,60,81
0000,040,160,360,641
10,0200,080,20,40,680,72
20,0400,10,240,440,560,74
30,0600,120,270,40,590,61
40,0800,1350,260,430,5050,63
50,100,130,28250,38250,530,5375
JTj /xi1,21,41,61,82
000,80,60,40,20
10,020,80,60,40,20
20,040,660,60,40,20
30,060,670,530,40,20
40,080,570,5350,3650,20
50,10,58250,46750,36750,18250

2. Пусть , тогда l=0,015- шаг по оси t, а h=0,3- шаг по оси x. Решение в виде (22) будем искать по формуле

(24)

В результате получим таблицу

Таблица №5

JTj /xi00,30,60,91,21,51,82
0000,090,360,810,80,50,20
10,01500,120,390,7333330,7516670,50,2166670
20,0300,1450,4022220,6791670,7066670,4947220,2277780
30,04500,1637040,4055090,6375930,6667590,4855560,2343060
40,0600,1767210,4038890,6037730,6316980,4738810,237130
50,07500,1851290,3993420,5751130,6007410,4607250,2370670

Рис.5- Решение, полученное с помощью метода сеток при

Рис.6- Решение, полученное с помощью метода сеток при

Рис.7- График точного решения, полученного аналитически

V. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ПЕРВОЙ СМЕШАННОЙ ЗАДАЧИ ДЛЯ ВОЛНОВОГО УРАВНЕНИЯ НА ОТРЕЗКЕ

Пусть дано волновое уравнение и его граничные условия

(25)

Решим задачу (25), применяя метод сеток для уравнений гиперболического типа.

Заменим производные в (25)

При (26)

Пусть , тогда по формуле (26) получим

Таблица №6

JTj /xi00,10,20,30,40,50,60,7
000-0,14-0,26-0,36-0,44-0,5-0,54-0,56
10,10-0,14-0,26-0,36-0,44-0,5-0,54-0,56
20,20-0,12-0,24-0,34-0,42-0,48-0,52-0,54
30,30-0,1-0,2-0,3-0,38-0,44-0,48-0,5
40,40-0,08-0,16-0,24-0,32-0,38-0,42-0,44
50,50-0,06-0,12-0,18-0,24-0,3-0,34-0,36
JTj /xi0,80,911,11,21,31,41,5
00-0,56-0,54-0,5-0,44-0,36-0,26-0,140
10,1-0,56-0,54-0,5-0,44-0,36-0,26-0,140
20,2-0,54-0,52-0,48-0,42-0,34-0,24-0,120
30,3-0,5-0,48-0,44-0,38-0,3-0,2-0,10
40,4-0,44-0,42-0,38-0,32-0,24-0,16-0,080
50,5-0,36-0,34-0,3-0,24-0,18-0,12-0,060

Рис.7- Решение волнового уравнения методом сеток при

Рис.8- График точного решения, полученного аналитически

VI. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ИНТЕГРАЛЬНОГО УРАВНЕНИЯ

Пусть дано интегральное уравнение

(27)

Будем искать решение уравнения (27) с помощью метода вырожденных ядер.

Представим ядро в виде ряда

Отбросим члены старше пятого порядка

Пусть , тогда

Таким образом, решение задачи (27)

ЗАКЛЮЧЕНИЕ

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

Найдено приближенное решение линейного дифференциального уравнения второго порядка с помощью методов Галеркина, Ритца и коллокации. Сравнение результатов приведено в виде таблиц и графиков.

С помощью метода Либмана получено приближенное решение задачи Дирихле в квадрате. Результаты приведены в виде таблиц.

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

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

Список использованных источников

1. В. Ф. Чудесенко Сборник заданий по специальным курсам высшей математики. Типовые расчеты: Учебное пособие. 4-е изд., стер.-СПб.: Издательство “Лань”, 2007.- 192с.: ил.- (Учебники для вузов. Специальная литература)

2.Вычислительная математика в примерах и задачах. Н. В. Копченова, И. А. Марон. Главная редакция физико-математической литературы изд-ва:Наука, М., 1972.

3. Тихонов, Самарский “Уравнения математической физики”, М.: Наука, 1967.


1 Star2 Stars3 Stars4 Stars5 Stars (No Ratings Yet)
Loading...
Приближенное решение интегрального уравнения