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

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

Рубрика Программирование, компьютеры и кибернетика
Вид курсовая работа
Язык русский
Дата добавления 27.11.2010
Размер файла 213,7 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

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

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

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

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

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

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

студента 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,…,msp, выполняется оценка
  • .
  • Пользуясь этой оценкой и замечая, что в силу (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. Практическая часть по решению линейных интегральных уравнений
  • Для написания функции, находящей решение линейного интегрального уравнения составим алгоритм. Представим алгоритм в виде блок-схемы.
  • Используя данную блок-схему, напишем соответствующую функцию.
  • Функция решения линейных интегральных уравнений будет реализована на С++.
  • 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


Подобные документы

  • Решение уравнения методом половинного деления. Программа в Matlab для уравнения (x-2)cos(x)=1. Решение нелинейных уравнений методом Ньютона. Интерполяция заданной функции. Решение системы линейных алгебраических и обыкновенных дифференциальных уравнений.

    курсовая работа [1,4 M], добавлен 15.08.2012

  • Приведение системы линейных алгебраических уравнений к треугольному виду прямым ходом метода Гаусса. Применение обратного хода метода вращений. Создание алгоритма, блок-схемы и кода программы. Тестовый пример решения уравнения и его проверка в MathCad.

    лабораторная работа [164,3 K], добавлен 02.10.2013

  • Решение систем алгебраических линейных уравнений методом Крамера. Сущность метода прогонки. Программная реализация метода: блок-схема алгоритма, листинг программы. Проверка применимости данного способа решения для конкретной системы линейных уравнений.

    курсовая работа [581,0 K], добавлен 15.06.2013

  • Алгоритм решения систем линейных уравнений методом Гаусса, его этапы. Система уравнений для определения коэффициентов сплайна, представляющая собой частный случай систем линейных алгебраических уравнений. Программная реализация, тестовый пример.

    курсовая работа [431,8 K], добавлен 15.06.2013

  • Характеристика влияния компьютера на здоровье человека. Определение корней уравнения в Microsoft Excel с точностью до шестого знака после запятой. Решение системы линейных уравнений методом вычисления определителей и матричным способом в Microsoft Excel.

    контрольная работа [734,0 K], добавлен 19.03.2012

  • Системы линейных алгебраических уравнений. Матричный метод решения систем линейных уравнений. Решение задачи математическим методом. Блок-схема алгоритма и листинг программы. Расчет трудоемкости разработки программы. Расчет себестоимости и цены программы.

    дипломная работа [144,8 K], добавлен 25.04.2012

  • Методы решения систем линейных алгебраических уравнений. Метод простых итераций и метод Зейделя. разработка программы для решения СЛАУ с произвольным количеством уравнений. Реализация методов Зейделя и простых итераций для получения вектора решений СЛАУ.

    курсовая работа [25,0 K], добавлен 20.11.2008

  • Численные методы решения задач. Решение алгебраических и трансцендентных уравнений. Уточнение корня по методу половинного деления. Решение систем линейных уравнений методом итераций. Методы решения дифференциальных уравнений. Решение транспортной задачи.

    курсовая работа [149,7 K], добавлен 16.11.2008

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

    курсовая работа [939,4 K], добавлен 16.01.2014

  • Отделение корней методом простых интеграций. Дифференцирование и аппроксимация зависимостей методом наименьших квадратов. Решение нелинейного уравнения вида f(x)=0 методом Ньютона. Решение системы линейных уравнений методом Зейделя и методом итераций.

    курсовая работа [990,8 K], добавлен 23.10.2011

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.