Банк рефератов содержит более 364 тысяч рефератов, курсовых и дипломных работ, шпаргалок и докладов по различным дисциплинам: истории, психологии, экономике, менеджменту, философии, праву, экологии. А также изложения, сочинения по литературе, отчеты по практике, топики по английскому.
Полнотекстовый поиск
Всего работ:
364150
Теги названий
Разделы
Авиация и космонавтика (304)
Административное право (123)
Арбитражный процесс (23)
Архитектура (113)
Астрология (4)
Астрономия (4814)
Банковское дело (5227)
Безопасность жизнедеятельности (2616)
Биографии (3423)
Биология (4214)
Биология и химия (1518)
Биржевое дело (68)
Ботаника и сельское хоз-во (2836)
Бухгалтерский учет и аудит (8269)
Валютные отношения (50)
Ветеринария (50)
Военная кафедра (762)
ГДЗ (2)
География (5275)
Геодезия (30)
Геология (1222)
Геополитика (43)
Государство и право (20403)
Гражданское право и процесс (465)
Делопроизводство (19)
Деньги и кредит (108)
ЕГЭ (173)
Естествознание (96)
Журналистика (899)
ЗНО (54)
Зоология (34)
Издательское дело и полиграфия (476)
Инвестиции (106)
Иностранный язык (62792)
Информатика (3562)
Информатика, программирование (6444)
Исторические личности (2165)
История (21320)
История техники (766)
Кибернетика (64)
Коммуникации и связь (3145)
Компьютерные науки (60)
Косметология (17)
Краеведение и этнография (588)
Краткое содержание произведений (1000)
Криминалистика (106)
Криминология (48)
Криптология (3)
Кулинария (1167)
Культура и искусство (8485)
Культурология (537)
Литература : зарубежная (2044)
Литература и русский язык (11657)
Логика (532)
Логистика (21)
Маркетинг (7985)
Математика (3721)
Медицина, здоровье (10549)
Медицинские науки (88)
Международное публичное право (58)
Международное частное право (36)
Международные отношения (2257)
Менеджмент (12491)
Металлургия (91)
Москвоведение (797)
Музыка (1338)
Муниципальное право (24)
Налоги, налогообложение (214)
Наука и техника (1141)
Начертательная геометрия (3)
Оккультизм и уфология (8)
Остальные рефераты (21697)
Педагогика (7850)
Политология (3801)
Право (682)
Право, юриспруденция (2881)
Предпринимательство (475)
Прикладные науки (1)
Промышленность, производство (7100)
Психология (8694)
психология, педагогика (4121)
Радиоэлектроника (443)
Реклама (952)
Религия и мифология (2967)
Риторика (23)
Сексология (748)
Социология (4876)
Статистика (95)
Страхование (107)
Строительные науки (7)
Строительство (2004)
Схемотехника (15)
Таможенная система (663)
Теория государства и права (240)
Теория организации (39)
Теплотехника (25)
Технология (624)
Товароведение (16)
Транспорт (2652)
Трудовое право (136)
Туризм (90)
Уголовное право и процесс (406)
Управление (95)
Управленческие науки (24)
Физика (3463)
Физкультура и спорт (4482)
Философия (7216)
Финансовые науки (4592)
Финансы (5386)
Фотография (3)
Химия (2244)
Хозяйственное право (23)
Цифровые устройства (29)
Экологическое право (35)
Экология (4517)
Экономика (20645)
Экономико-математическое моделирование (666)
Экономическая география (119)
Экономическая теория (2573)
Этика (889)
Юриспруденция (288)
Языковедение (148)
Языкознание, филология (1140)

Курсовая работа: Решение линейных интегральных уравнений

Название: Решение линейных интегральных уравнений
Раздел: Рефераты по информатике, программированию
Тип: курсовая работа Добавлен 04:06:00 28 ноября 2010 Похожие работы
Просмотров: 148 Комментариев: 2 Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать

Федеральное агентство по образованию

Тульский Государственный педагогический университет

имени Л. Н. Толстого

Кафедра информационных технологий

Курсовая работа

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

студента 4 курса группы В

специальности 351500 – МОиАИС

Селиванова Сергея Валериевича

Тула – 2008


Оглавление

Введение

1. Теоретическая часть по решению линейных интегральных уравнений

2. Практическая часть по решению линейных интегральных уравнений

Заключение

Используемые источники


Введение

В данной курсовой работе рассмотрена проблема решения линейных интегральных уравнений. Целью курсовой работы было написание функции, которая по введенным данным (ядру интегрирования, правой части уравнения и отрезку интегрирования) могла бы находить решения линейного интегрального уравнения. Проблема разработки алгоритма решения и написании на его основе функции является практически актуальной, так как решение линейных интегральных уравнений без привлечения ЭВМ является достаточно трудоемким.

Данная курсовая работа состоит двух частей.

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

Во второй главе приводится алгоритм решения линейного интегрального уравнения и, написанной на его основе, функции.


1. Теоретическая часть по решению линейных интегральных уравнений

Существует множество методов решений линейных интегральных уравнений. Рассмотрим один из них – метод итераций.

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

(1)

Будем предполагать, что свободный член и ядро этого уравнения принадлежат соответствующим классам и . Уравнение (1) будем также записывать кратко в виде

, (2)

где интегрирование распространенно на единичный r-мерный куб Gr .

Лемма 1. Если

и (3)

то при решение уравнения (2) удовлетворяет соотношению

,

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


(4)

Принадлежит классу.

Доказательство.

Известно, что при достаточно маломλ решение уравнения (2) можно представить в виде ряда

где Grv -единичный rv-мерный куб. пусть величина Rn определена равенством

.

Тогда пользуясь определением функции F(P,Q1 ,…,Qn ) получим

(5)

Обозначим через С(m1 ,…,mr ) коэффициенты Фурье функции f(P). Так как, по условию, f(P) , то

Аналогичная оценка справедлива, очевидно, и для ядра K(P,Q) уравнения (2) .

Но тогда

и, следовательно,

получим

,

.

Отсюда в силу (5) следует первое из утверждений леммы:

.

Перейдем теперь к доказательству второго утверждения. Так как f(P) и K(P,Q),то, аналогично рассуждениям леммы 12 (1, с.61) легко показать, что

(6)


Где,

В отличие от остальных сомножителей, первый сомножитель в соотношении (6) рассматривается как функция r переменных, соответствующих величине Q1 , а не как функция всех своих переменных.

Далее, рассматривая каждую из функций (v=1,2,…,n)

Как функцию всех rn переменных, соответствующих n величинам Q1 ,…,Qn , согласно первому утверждению леммы 12 (1, с.61) получим, что функция принадлежит классу , где

.

Но в силу (6)

и, следовательно,

.

Чем лемма 1 доказана полностью.

Пусть, как и выше f(P) и K(P,Q),

(7)

и величина γ0 определена равенством (3)

Покажем, что для приближённого решения уравнения (7) можно использовать квадратные формулы с неравномерными сетками.

Теорема 1. Пусть p- простое число, N=p, и величина n определена равенством

Тогда при произвольно малом ε для решения уравнения (7) выполняется асимптотическое равенство

где

Доказательство.

Пусть функция Φ принадлежит классу и σ-сумма модулей её коэффициентов Фурье. Тогда согласно теореме 15 (1, с.94) справедлива квадратурная формула


, (8)

где

(9)

Выберем в лемме 1 . Тогда при для решения уравнения получим

(10)

где согласно (4)функция F(P,Q1 ,…,Qn ) определена равенством F(P,Q1 ,…,Qn )=и принадлежит классу.

Пусть при k=1,2,…,N и v=1,2,…,n точки Mk , v определены равенством

.

Выберем p настолько большим, чтобы выполнялись неравенства n≥1 и N≥rn

Тогда применяя квадратичную формулу (8) получим


(11)

где в силу (9)

(12)

Пользуясь определением n и , получим

.

.

Следовательно,

, , .

В силу (12)

Но тогда из (10) и (11) следует, что


Отсюда, пользуясь оценкой

,

получаем утверждение теоремы.

Результат, полученный в теореме 1, можно усилить, если воспользоваться методом оптимальных коэффициентов.

Лемма 2. Для всякого простого p существуют оптимальные коэффициенты a1 ,…,as такие, что каково бы ни было a>1+ε1 , при любом ε1 (0;1) выполняется оценка

Доказательство.

Пусть z-произвольное целое из интервала Определим функцию Тs (z) равенством

Пусть при z=a достигается минимум этой функции. Тогда, очевидно,


(13)

Так согласно лемме 1(1, с.21)

,

то при произвольном ε > 0 получим из (13),

Отсюда следует, что

(14)

Введём обозначения

Так из (14) в силу определения величины Ts (a) следует оценка

(15)


то пользуясь неравенством, получим

(16)

Чтобы оценить сумму Σ2 , заметим, что для нетривиальных решений сравнения

(17)

Выполняется неравенство

(18)

Действительно, согласно определению величины δp (m) в левой части неравенства (14) отличны от нуля только такие слагаемые, для которых m1 ,…,ms является нетривиальным решением сравнения (5.43). так как любое из этих слагаемых не превосходит всей суммы, то для каждого нетривиального решения сравнения получим

,

Чем неравенство (5.44) доказано.

Пусть функция φ(m1 ,…,ms ) определена равенствами


Тогда пользуясь леммой 18 (1, c.101), получим

. (19)

Обозначим через q минимальное значение произведения , где m1 ,…,ms –произвольное нетривиальное решение сравнения (17).

Тогда, выбирая в лемме 26 (1, с.151)

,

получим, что при любых натуральных m1 ,…,ms , удовлетворяющих условию m1 ,…,ms p, выполняется оценка

.

Пользуясь этой оценкой и замечая, что в силу (18)


при любом ε ≤ a-1 положительном получим из (19)

(20)

Выберем av =av-1 (v=1,2,…,s) (21)

тогда, пользуясь оценками (16) и (20), при получим неравенство, указанное в лемме:

Для завершения доказательства леммы остается убедиться, что величины , определенные равенством (21), являются оптимальными коэффициентами.

Действительно, из (5.39), пользуясь леммой 1(1, c.21) получим

Переписывая эту оценку в виде

убеждаемся, что целые av =av -1 будут оптимальными коэффициентами, чем лемма 2 доказана полностью.

Следствие . Если Ф, то, каково бы ни было для погрешности квадратурной формулы

построенной при N=p с помощью оптимальных коэффициентов, указанных в лемме 2, справедлива оценка

,

Действительно, пользуясь леммами 19 (1, с.106) и 2, получим утверждение следствия

Пусть α>0, , p - простое, N=p, a1 ,…as – оптимальные коэффициенты по модулю p, удовлетворяющие условию леммы 2, и величины γ0 , n определены равенствами

(22)

Теорема 2 Если, то при произвольно малом для решения уравнения


(23)

выполняется равенство

где

(24)

Доказательство.

Выберем в лемме 1 , где γ0 определено первым из равенств (22). Тогда для решения уравнения (23) получим

(25)

где функция F(P,Q1 ,…,Qn ) определена равенством

И принадлежит классу

Пусть при k=1,2,…,N и v=1,2,…,n точки Mk , v определены равенством (24). Тогда согласно квадратурной формуле, указанной в следствии леммы 2, при s=rn и справедливо равенство


(26)

(27)

Пользуясь равенствами (22), получим

Но тогда

и, следовательно

Пользуясь этой оценкой, из (25) и (26) получим


Отсюда, так как в силу выбора n выполняется оценка

Следует утверждение теоремы.

2. Практическая часть по решению линейных интегральных уравнений

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





y[i]=B[i,m];

Используя данную блок-схему, напишем соответствующую функцию.

Функция решения линейных интегральных уравнений будет реализована на С++.

bool solvefredholm2(const double& a,

const double& b,

const int& n,

ap::real_1d_array& y,

const double& epsilon)

{

bool result;

double h;

double t;

double m1;

double x;

ap::real_2d_array smat;

int i;

int j;

int u;

int k1;

int m;

smat.setbounds(1, n, 1, n+1);

y.setbounds(1, n);

h = (b-a)/(n-1);

i = 1;

do

{

x = a+(i-1)*h;

smat(i,n+1) = f(x);

j = 1;

do

{

smat(i,j) = -h*k(x, a+(j-1)*h);

if( j==1||j==n )

{

smat(i,j) = smat(i,j)/2;

}

if( j==i )

{

smat(i,j) = 1+smat(i,j);

}

j = j+1;

}

while(j<=n);

i = i+1;

}

while(i<=n);

y.setbounds(1, n);

result = true;

for(i = 1; i <= n; i++)

{

k1 = i;

m1 = fabs(smat(i,i));

for(j = i+1; j <= n; j++)

{

if( m1<fabs(smat(j,i)) )

{

m1 = fabs(smat(j,i));

k1 = j;

}

}

if( fabs(m1)>=epsilon )

{

for(j = i; j <= n+1; j++)

{

t = smat(i,j);

smat(i,j) = smat(k1,j);

smat(k1,j) = t;

}

for(k1 = i+1; k1 <= n; k1++)

{

t = smat(k1,i)/smat(i,i);

smat(k1,i) = 0;

for(j = i+1; j <= n+1; j++)

{

smat(k1,j) = smat(k1,j)-t*smat(i,j);

}

}

}

else

{

result = false;

break;

}

}

if( result )

{

i = n;

do

{

y(i) = smat(i,n+1);

j = i+1;

while(j<=n)

{

y(i) = y(i)-smat(i,j)*y(j);

j = j+1;

}

y(i) = y(i)/smat(i,i);

i = i-1;

}

while(i>=1);

}

return result;

}


Данная функция решает интегральное уравнение Фредгольма второго рода, заданное ядром интегрирования K(X,S) и правой частью F(X), на отрезке [A, B] методом итераций.

Результат помещается в массив Y с номерами элементов от 1 до N, где 1 соответствует A, N соответствует B.

Epsilon - малое число, передаваемое для сравнения с нолем в ходе решения получаемой системы уравнений.

Для работы этой функции необходима библиотека ap.h


Заключение

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


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

1. Коробов, Н. М. Теоретико-числовые методы в приближенном анализе/ Н. М. Коробов. –М.: 2003. - 316 с.

2. Коробов, Н. М. О приближенном решении интегральных уравнений/

Н. М. Коробов. –ДАН СССР, 1959.

3. http://alglib.sources.ru

Оценить/Добавить комментарий
Имя
Оценка
Комментарии:
Где скачать еще рефератов? Здесь: letsdoit777.blogspot.com
Евгений06:40:53 19 марта 2016
Кто еще хочет зарабатывать от 9000 рублей в день "Чистых Денег"? Узнайте как: business1777.blogspot.com ! Cпециально для студентов!
19:48:19 28 ноября 2015

Работы, похожие на Курсовая работа: Решение линейных интегральных уравнений

Назад
Меню
Главная
Рефераты
Благодарности
Опрос
Станете ли вы заказывать работу за деньги, если не найдете ее в Интернете?

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



Результаты(150916)
Комментарии (1842)
Copyright © 2005-2016 BestReferat.ru bestreferat@mail.ru       реклама на сайте

Рейтинг@Mail.ru