Бюллетень науки и практики — Bulletin of Science and Practice
научный журнал (scientific journal) №12 2017 г.
http://www. bulletennauki. com
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ / PHYSICAL & MATHEMATICAL SCIENCES
УДК 519.64: 519.65
ПРИМЕНЕНИЕ МАТЕМАТИЧЕСКОГО ПАКЕТА К РЕШЕНИЮ ИНТЕГРАЛЬНЫХ
УРАВНЕНИЙ
THE USE OF MATHEMATICAL SOFTWARE PACKAGE FOR SOLVING INTEGRAL
EQUATIONS
©Шувалова Л. Е.
Казанский национальный исследовательский технологический университет г. Нижнекамск, Россия, leshyvalova@yandex.ru
©Shuvalova L.
Kazan National Research Technological University Nizhnekamsk, Russia, leshyvalovayjandex.ru
©Валиуллин А. В. Казанский национальный исследовательский технологический университет г. Нижнекамск, Россия, aidarval111@mail.ru
©Valiullin A.
Kazan National Research Technological University Nizhnekamsk, Russia, aidarval111@mail.ru
Аннотация. Рассмотрен метод механических квадратур решения интегральных уравнений с численной реализацией в математическом пакете MathCad. Построены интерполяционные многочлены, аппроксимирующие искомую функцию. Авторами предлагается программа, которая позволяет автоматически находить приближенные решения интегральных уравнений и при других исходных данных. Расчеты. выполненные в программе дают возможность минимизировать погрешность путем специального выбора узлов сетки.
Abstract. The method of mechanical quadratures for solving integral equations with numerical realization of mathematical package MathCad. Constructed interpolation polynomial approximating the desired function. The authors propose a program that allows you to automatically find approximate solutions of integral equations and other basic data. Calculations, made in the program make it possible to minimize the error by a special selection of mesh nodes.
Ключевые слова: интегральное уравнение, квадратурная формула, математический пакет.
Keywords: integral equation, quadrature formula, mathematical package.
Теория интегральных уравнений Фредгольма II рода хорошо разработана [1]. Как известно, точное аналитическое решение для таких классов задач можно найти лишь в частных случаях. В работе Л. Е. Шуваловой и Л. А. Апайчевой (2013) рассмотрен метод
научный журнал (scientific journal) №12 2017 г.
http://www. bulletennauki. com
механических квадратур для решения нелинейных сингулярных интегральных уравнений, но каждый раз приходится сталкиваться с трудной реализацией вычислительного процесса. Поэтому в данной статье сделана попытка выполнить расчеты в математическом пакете МаШСаё, который позволяет эффективно решать технические задачи.
В общем случае интегральное уравнение Фредгольма II рода представляется в виде
У(Х) -р{а К(х, * уф ^ = /(х) , (1)
где к(х,б), /(х) — известные непрерывные функции в своих областях определения, р — параметр уравнения, а у(х) — искомая функция.
Приближенное решение уравнения (1) ищется методом механических квадратур. Для этого на отрезке [а, Ь] выбираются узлы:а < х1 < х2 < ■■■ < хп = Ь.
Уравнение (1) записывается в узлах сетки:
у(х1)-л{ак(хи8)*у(8) йБ = /ОО , ¿ = 1,П. (2)
Интегралы в равенстве (2) вычисляются приближенно по квадратурным формулам трапеций и Симпсона:
^ <р(х) йх « Т11}=1А] * <р(х¡) . (3)
Получаем систему линейных алгебраических уравнений (СЛАУ):
У1 -Л1$=1А * КИ *У1=П,1 = 1,п. (4)
Решение СЛАУ (4) дает значения у±, / = 1,п, по которым путем интерполяции находится приближенное решение интегрального уравнения (1).
С помощью математического пакета МаШсаё решается следующий тестовый пример
11
У(х) -1f0(x + 2s) * y(s) ds = ex +
■ — X f — У
2->0 ^ ' ' ^ ' 2 Выбираются узлы сетки и вычисляются в них значения правой части
[(х) = ех + х 2 и ядра К(х, б) = х + 2э .
Приведем часть программы:
научный журнал (scientific journal) №12 2017 г.
http://www. bulletennauki. com
Коэффициенты квадратурой формулы (3) А[ вычисляются либо по формулам трапеции, либо по формулам Симпсона. Имеем
Далее формируются матрицы С и D, состоящие из коэффициентов СЛАУ (4) и значений правой части в выбранных узлах /(х^):
научный журнал (scientific journal) №12 2017 г.
http://www. bulletennauki. com
В результате работы данного блока программы получаем приближенные значения искомой функции в узлах xi с применениями формул Симпсона и трапеций соответственно:
Для анализа полученных численных результатов применяем интерполяцию с помощью многочлена Лагранжа и кубического сплайна.
Полученные приближенные решения отразим на графике, сравнив их с точным решением (Рисунок).
научный журнал (scientific journal) №12 2017 г.
http://www. bulletennauki. com
Рисунок. Полученные приближенные и точные решения
Данная программа позволяет автоматически находить приближенные решения интегральных уравнений и при других исходных данных. Кроме того, дает возможность минимизировать погрешность путем специального выбора узлов сетки.
Список литературы:
1. Верлань А. Ф., Сизиков В. С. Интегральные уравнения: методы, алгоритмы, программы. Киев: Наукова Думка, 1986. 543 с.
2. Шувалова Л. Е., Апайчева Л. А. Приближенное решение одного класса нелинейных сингулярных интегральных уравнений // Вестник Казанского технологического университета. 2013. Т. 16. №12. С. 289-292.
Reference:
1. Verlan, A. F., & Sizikov V. S. (1986). Integral equations: methods, algorithms, programs. Kiev, Naukova Dumka, 543. (in Russian)
2. Shuvalova, L. E., & Apaicheva L. A. (2013). The approximate solution of a class of nonlinear singular integral equations. Vestnik Kazanskogo technologichescogo Universiteta, 16, (12), 289-292. (in Russian)
Бюллетень науки и практики — Bulletin of Science and Practice научный журнал (scientific journal)
http://www. bulletennauki. com
№12 2017 г.
Работа поступила в редакцию 22.11.2017 г.
Принята к публикации
26.11.2017 г.
Ссылка для цитирования:
Шувалова Л. Е., Валиуллин А. В. Применение математического пакета к решению интегральных уравнений // Бюллетень науки и практики. Электрон. журн. 2017. №12 (25). С. 13-18. Режим доступа: http://www.bulletennauki.com/shuvalova (дата обращения 15.12.2017).
Cite as (APA):
Shuvalova, L., & Valiullin, A. (2017). The use of mathematical software package for solving integral equations. Bulletin of Science and Practice, (12), 13-18