CC BY
26
ИII ««I'flATH Ki вычнслнтнл елля mil ныл Н УПРАВЛЕНИЕ
|lt С ||]Ш||Г. 11,4-1
УДК 532.5:621.694:519.6
Горская Т.Ю. - кандидат технических наук, доцент
E-mail: [email protected]
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1
Галимянов А.Ф. - кандидат физико-математических наук, доцент
E-mail: [email protected]
Казанский (Поволжский) федеральный университет
Адрес организации: 420008, Россия, г. Казань, ул. Кремлевская, д. 18
Аппроксимация дробных интегралов частными суммами ряда Фурье
Аннотация
Постановка задачи. Целью настоящего исследования является разработка и применение приближенного метода для вычисления интегралов, являющихся составной частью моделей, использующих интегралы Римана-Лиувилля, а также создание программного продукта, позволяющего проводить подобные вычисления для заданных функций.
Результаты. Основные результаты исследования состоят в построении квадратурной формулы для интеграла, при этом рассматривались случаи, когда плотность интеграла есть функция из пространств непрерывных функций, имеющих обобщенные производные с весом, и гельдеровых классов функций с весом. Для предложенной квадратурной формулы далее исследовалась погрешность ее приближения в пространствах непрерывных функций и квадратично-суммируемых с весом функций. В результате исследования установлены эффективные оценки погрешности аппроксимативного аппарата в предложенных классах функций. Кроме того приближенный метод реализован на компьютере в виде программы на языке СИ.
Выводы. Значимость полученных результатов для строительной отрасли состоит в том, что при решении задач, в том числе задач по нахождению форм конструкций, учитывающих свойства материалов, изменения окружающей среды, модели которых используют интегралы Римана-Лиувиля, можно будет применять приближенный подход, квадратурную формулу, предложенную в статье.
Ключевые слова: интегралы Римана-Лиувилля, ряды Фурье, квадратурные формулы, приближенные вычисления, оценки погрешности.
В последние десятилетия возрос интерес к исследованию дифференциальных уравнений, имеющих дробный порядок, в которых неизвестная функция содержится под знаком производной дробного порядка, а также и к интегральным уравнениям дробного порядка интегрирования. Это обусловлено рядом причин: во-первых, развитием как таковой теории дробного интегрирования и дифференцирования, во-вторых, обширными приложениями применения этого математического аппарата в различных областях науки и производства, особенно, в областях, связанных с нанотехнологиями, задачами диффузии, а также при создании конструкций, учитывающих состояние вещества [1, 2]. Кроме того разными авторами ранее проводились исследования в области приближенных методов дробных интегралов [3-5]. Нами также проводились работы, основанные на известных методологических подходах для решения задач приближения, сначала проекционно-сеточными методами [6, 7], затем по построению квадратурных формул для интегралов [8-10].
Вначале укажем некоторые определения и вспомогательные предложения, которыми будем оперировать при изложении материала.
Напомним, что функция /(х) 6 #$r)(r >0),хЕ (-со; +да), имеющая производные до r-го порядка с весом р (х): /$(l) (х) = р (х) *+/$(l_l) (х),. = 1 , и для нее выполняется условие -да < /$(0) (да) = /$(0) (-да) < да, j = 0 , r - 1. Здесь и далее вес р (х) = ———з-
Функцию /(х) представим рядом Фурье по системе функций {ск(х),s5(х)}, где ск (х) = cos(2karсtgх), s5 (х) = зт(2кагс£дх):
с
/(х) = у + < акСк (х) + Ьк Sk (х),
к = l
2 (f(x)ck Щх
ак =- I ---в-,к = °>1>-
п J 1+ хв
—со ад
2 Г f(x) Sk (х) Ах
bk = - I , 5 в ,к = 1,2,..
2 Cf(x)t п J 1
— сс
Обозначим через:
n
Cnf = Cn(f, х) = у + <акСк (х) + bk Sk (х) - (1)
к = 1
частную сумму ряда Фурье по системе функций {ck (х), Sk (х)}.
Для функции f (х) Е Wp(r\г > 1,х Е (-с; +со) справедливо следующее: Лемма 1.
Пусть функция f (х) EW$r\ г> 1, хЕ (-ж;+со), тогда для нее имеет место представление:
с с
Пх)= п if(t)r+9 + n j^WrMj+B,
где
К% (х, t) = <Ckr (t)Ck (х) + Skr (tfSkte), 1 k=1 1 (t) = — cos H2arctgt +-2-), Skr (t) = —sin (2arctgt +-2-)-
7% у и 2 )7% \ и 2 Из свойства ортогональности функций [с^^х),Зк(х)} с весом р на всей числовой прямой следует следующее: Лемма 2.
. ь(r) & > 1 то E (f) < E Mf(r)) E (К ) если Г(х) Е W(r)
En (f)c<MEn(Kr)Lp.
Следствие 1.
,(r)
то
9пА
En (f)c <
Если fM Е C$r), г > 1, то En (f)c < En{f;r))En(Kr)Lp, если fU) Е W$r\ г > 1, то
Если f(х) Е НП (А), г> 0,0< U <1,
2(п + 1)r+S' Следствие 2.
,(r)
Если Г(х) Е Ш/',& > 1, то Еп(Т)с <
Обозначим через Н = Ь2$р — линейное нормированное пространство квадратично-суммируемых функций с весом р. Напомним, здесь и далее в статье р(х) = Введем
...... Il гс _2(t)dt
норму \\f\\H = Jwf—с—-
Для суммы (1) справедливы следующие оценки, представленные в виде леммы: Лемма 3.
Для любых натуральных п = 1,2,... справедливы соотношения:
\\ип\\с^н < 1, \\ип\\с^с < 2 +1пп.
Справедлива также следующая теорема: Теорема 1.
Если [(х) еС(—ж;+сю), то справедлива оценка погрешности аппроксимации функции частными суммами ряда Фурье (1) в пространстве Н:
М — ипП\н< 2Еп (г)с.
Следствие 1.
Если [(.х) Е (А),г> 0,0< р <1,
то
9?А
\\Г — ипПн<
(п + 1)r+P'
Следствие 2.
''р
Если f(x) е 'P(r), & > 1,
то
пМ
\\/-спП\< (V7D-
Доказательство теоремы 1 следует из известной теоремы Лебега и оценок леммы 3, а следствия 1, 2 теоремы 1 следуют непосредственно из следствий 1, 2 леммы 3, соответственно.
Также справедлива следующая теорема:
Теорема 2.
Если /(х) 6 С(-да;+оо), то справедлива оценка погрешности аппроксимации функции частными суммами ряда Фурье (1) в пространстве непрерывных функций:
\\/-ип/\\с < (3+lnn)J„(/)с.
Следствие 1.
то
Если f(x) е H¡f$ (T), &> 0,0< U < 1,
9п(3 + lnn)T
Wf-Cnf\\H< 2(v + i)r+p •
Следствие 2.
Если f(x) е 'P(r), & > 1,
то
пМ
\\f-Cnf\\
н < 2(п + 1)% (3 +lnv). Утверждения теоремы 2, а равно следствия 1, 2 следуют непосредственно из леммы 3 и ее следствий.
Результаты
Рассмотрим интеграл Римана-Лиувилля на бесконечном промежутке:
(/Cf)(x) =rd) @ ^,—ю <x <ю (2)
—ю
Здесь функция f (x), заданная на числовой оси (—ю; +ю)
Аппроксимируя плотность интеграла (2) частной суммой ряда Фурье (1), придем к следующей приближенной формуле:
1 г(8т + Yn=íakck(x) + bksk (x)gdt ao "
I?(Unf;x) = J-(x — t)1—a-= y + <akCkc(x) + bkskc(x), (3)
—ю k=1
где
Cnc (x) =(2П)? C0S (2na&ctgx +—) , Sn„ (t) = (^Щ)? sin (2na&ctgx +—),
a = B r® f(x)ck(x)dx 7 = 01 b = B [ю f(x)sk(x)dx 7 = 12 k я^—ю 12x2 , , , •• k я-'—ю 12x2 , , , ••
Для формулы (3)-(4) имеет место следующее: Лемма 4.
Если f (x) е С(—ю; +ю), то справедливо соотношение:
\\/?(f — Unf)\\Н < (¿c\\f — Unf \\H•
Скорость сходимости приближенного процесса (3) по метрике пространства H характеризует:
3. Если f (x) е hp р,
\ \ Ic(f — Unf) \ \ н = о{(п + 1)Р+Гпа}
Если f (x) е 'P(r), & > 1, то
у/?(f — Unf)\\H = 0 С(тП+^Т^'пс)•
Утверждения теоремы следуют непосредственно из леммы 4 и следствий 1 и 2 теоремы 1 соответственно.
Теорема 3. Если f (x) е hSP (T), (& > 0,0< U < 1), то справедлива оценка:
Далее рассмотрим оценки погрешности приближенной формулы (3) в равномерной метрике, для доказательства которых воспользуемся утверждением следующей леммы. Лемма 5. Для любых натуральных п = 1,2,... справедливо соотношение:
\\1Сип\\с^с< — (1+Ып).
Теорема 4. Если f(x) Е (А), (г > 0,0< U < 1), то справедлива оценка: , 18?А ( 2 1+ 2~Г~Р )
№(f-cnf)\\c< 2C(^+T)FT%[1+lnn 1 - 2-r-p + i - 2-r-p)2 ln2j.
Если f(x) Е Wp(r\ г > 1, то:
2?М ( 2 1+2~r л
\\It(f -unf)\\c <—-—h + lnn-—— + -——-ln2O.
n+v njjnc~ 2a(n + 1)r( 1 - 2~r (1 - 2~r)2 )
Оценивая остаточный член формулы (3) и последовательно применяя следствие 1,
2 леммы 2 и результаты леммы 5, получаем утверждение теоремы.
Выводы
В работе предложена квадратурная формула (3) с коэффициентами (4) для интеграла Римана-Лиувилля на бесконечном промежутке (2). Установлены скорости сходимости аппроскимационного аппарата (3) в метриках пространств непрерывных функций и квадратично суммируемых с весом функций (соответственно теоремы 4 и 3). Кроме того, данные оценки получены для приближения интеграла (2) в случаях плотностей f(x) Е Нр\(А), (г > 0,0< U< 1) и f (x) Е W$r,г > 1. Считаем, что данные исследования будут интересны в качестве теоретического обоснования применения квадратурных формул для интегралов Римана-Лиувилля на бесконечном промежутке. Кроме того, авторами разработана программа, реализующая приближенный метод на языке программирования СИ.
Список библиографических ссылок
1. Нахушев А. М. Элементы дробного исчисления и их приложения. Нальчик : НИИ ПМА КБНЦ РАН. 2009.
2. Салимов Р. Б. Асимптотическое представление сингулярного интеграла с ядром Гильберта вблизи точки слабой непрерывности плотности // Известия вузов. Математика. № 7. 2015, С. 58-62.
3. Marinov T. M., Ramírez N., Santamaría F. Fractional integration toolbox // Fractional Calculus and Applied Analysis, 2013. V. 16, № 3. P. 670-681.
4. Barton T. A., Purnaras I. K. Zp-solutions of singular integro-differential equations // J. Math. Anal. Appl., 2012, № 386. P. 830-841.
5. Saeed R. K., Ahmed C. Approximate solution for the system of non-linear Volterra integral equations of the second kind by using block-by-block method // Australian Journal of Basic and Applied Sciences, 2008.V. 2, № 1. P. 114-124.
6. Марчук Г. И., Агошков В. И. Введение в проекционно-сеточные методы. М. : Наука, 1981. 416 с.
7. Горская Т. Ю., Ожегова А. В. О сходимости проекционного метода для уравнения задачи движения // Известия КГАСУ, 2013, № 2 (24). С. 112-126.
8. Галимянов А. Ф., Сайфуллина Д. Э. Квадратурный метод решения интегрального уравнения смешанного типа // Известия Вузов. Математика, 2009, № 12. С. 22-27.
9. Галимянов А. Ф., Горская Т. Ю. Обобщенный метод Бубнова-Галеркина для уравнений с дробно-интегральным оператором // Известия КГАСУ, 2013, № 2 (24). С. 112-126.
10. Горская Т. Ю., Галимянов А. Ф., Воронцова В. Л. Квадратурная формула для гиперсингулярного интеграла с логарифмически ослабленным ядром : Тезисы докладов 69 международной научной конференции по проблемам архитектуры и строительства / КГАСУ, Казань. 2017. С. 373.
Gorskaya T.Yu. - candidate of technical sciences, associate professor E-mail: [email protected]
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya st., 1 Galimyanov A.F. - candidate of physical-mathematical sciences, associate professor E-mail: [email protected] Kazan Federal University
The organization address: 420008, Russia, Kazan, Kremlevskaya st., 18
Approximation of fractional integrals by private sums of the Fourier series Abstract
Problem statement. The purpose of this study is the development and application of approximate methods for calculating integrals, which is the part of the models used by the integrals of Riemann-Liouville and the creation of a software product, allowing to carry out similar calculations for the given functions.
Results. Main results of the research consist in the construction of quadrature formulas for the integral, there were cases when the density of the integral is a function from space of continuous functions having generalized derivatives with the weight, and Holder classes of functions with weight. For the proposed quadrature formulas next, we investigated the error of its approximation in spaces of continuous functions and square-summable with weight functions. The study has efficient error estimates approximation of the apparatus in the proposed classes of functions. In addition an approximate method implemented on a computer in the form of a program in the C language.
Conclusions. The significance of the results for the construction industry is that when solving the tasks, including the task of finding shapes of structures, considering material properties, environmental changes, models which use integrals of Riemann-Liouvile, you can apply the approximate approach, the quadrature formula proposed in the article.
Keywords: Riemann-Liouville integrals, Fulier series, quadrature formulas, approximate calculations, error estimates.
References
1. Nakhushev A. M. Elements of fractional calculus and their applications. Nalchik : NII PMA KBSC RAS. 2009.
2. Salimov R. B. Asymptotic representation of the singular integral with the Hilbert kernel near a point of weak continuity of the density of // Izvestiya vuzov. Math. № 7. 2015. P. 58-62.
3. Marinov T. M., Ramirez N., Santamaria F. Fractional integration toolbox // Fractional Calculus and Applied Analysis, 2013. V. 16, № 3. P. 670-681.
4. Barton T. A., Purnaras I. K. Lp-solutions of singular integro-differential equations // J. Math. Anal. Appl., 2012, № 386. P. 830-841.
5. Saeed R. K., Ahmed C. Approximate solution for the system of non-linear Volterra integral equations of the second kind by using block-by-block method. Australian Journal of Basic and Applied Sciences, 2008. V. 2, № 1. P. 114-124.
6. Marchuk G. I., Agoshkov V. I. Introduction to projection-grid methods. M. : Nauka, 1981.416 p.
7. Gorskaya T. Yu., Ozhegova A. V. On the convergence of the projection method for the equation of the motion problem // Izvestiya KGASU, 2013, № 2 (24). P. 112-126.
8. Galimyanov A. F., Sayfullin D. E. Quadrature method for solving an integral equation of mixed type // Izvestiya Vuzov. Matematika, 2009, № 12. P. 22-27.
9. Galimyanov A. F, Gorskaya T. Yu. A generalized Bubnov-Galerkin method for equations with a fractional-integral operator // Izvestiya KGASU, 2013, № 2 (24). P. 112-126.
10. Gorskaya T. Yu., Galimyanov A. F., Vorontsov V. L. Quadrature formula for a hypersingular integral with a logarithmically weakened core : Theses of reports of the 69th International scientific conference on architecture and construction / KGASU, Kazan. 2017. P. 373.
