Analysis of errors in some recent numerical quadrature formulas for periodic singular and hypersingular integrals via regularization

被引:8
作者
Sidi, Avram [1 ]
机构
[1] Technion Israel Inst Technol, Dept Comp Sci, IL-32000 Haifa, Israel
关键词
Cauchy Principal Value; Hadamard Finite Part; Circular Hilbert transform; Hypersingular integral; Numerical quadrature; Trapezoidal rule; Euler-Maclaurin expansion; Regularization; EULER-MACLAURIN EXPANSIONS; END-POINT SINGULARITIES;
D O I
10.1016/j.apnum.2014.02.011
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
Recently, we derived some new numerical quadrature formulas of trapezoidal rule type for the singular integrals I(1)[u] = integral(b)(a)(cot pi(x-t)/T)u(x) dx and I-(2)[u] = integral(b)(a)(csc(2) pi(x-t)/T)u(x) dx with b - a = T and u(x) a T-periodic continuous function on R. These integrals are not defined in the regular sense, but are defined in the sense of Cauchy Principal Value and Hadamard Finite Part, respectively. With h = (b - a)/n, n = 1, 2,..., the numerical quadrature formulas Qn((1))[u] for I((1)[)u] and Qn((2))[u] for I-(2)[u] are Qn((1))[u] = h (n)Sigma(j=1) f(t+jh-h/2), f(x) = (cot pi(x-t)/T)u(x), and Qn((2))[u] = h (n)Sigma(j=1) f(t+jh-h/2) - T(2)u(t)h(-1), f(x) = (csc(2) pi(x-t)/T)u(x) We provided a complete analysis of the errors in these formulas under the assumption that u is an element of C-infinity (R) and is T-periodic. We actually showed that, I-(1)[u] - Q(n)((1))[u] = O(n(-mu)) and I-(2)[u] - Q(n)((2))[u] = O(n(-mu)) as n ->infinity, for all mu > 0. In this note, we analyze the errors in these formulas under the weaker assumption that u is an element of C-s (R) for some finite integer s. By first regularizing these integrals, we prove that, if u((s+1)) is piecewise continuous, then I-(1)[u] - Q(n)((1))[u] = O(n(-s-1/2)) as n ->infinity, if s >= 1, and I-(2)[u] - Q(n)((2))[u] = O(n(-s-1/2)) as n ->infinity, if s >= 2 We also extend these results by imposing different smoothness conditions on u((s+1)). Finally, we append suitable numerical examples. (C) 2014 IMACS. Published by Elsevier B.V. All rights reserved.
引用
收藏
页码:30 / 39
页数:10
相关论文
共 14 条
  • [1] Abramowitz M., 1964, NAT BUR STAND APPL M, V55
  • [2] [Anonymous], 1984, Methods of Numerical Integration
  • [3] Davis PJ, 1975, INTERPOLATION APPROX
  • [4] Evans G., 1993, PRACTICAL NUMERICAL
  • [5] Kythe P. K., 2005, HDB COMPUTATIONAL ME
  • [6] Lifanov I.K., 2003, DIFF URAVN, V39, p[1175, 1115]
  • [7] Lifanov I.K., 2006, DIFF URAVN, V42, p[1134, 1263]
  • [8] Sidi A., 1988, Journal of Scientific Computing, V3, P201, DOI 10.1007/BF01061258
  • [9] Euler-Maclaurin expansions for integrals with endpoint singularities: a new perspective
    Sidi, A
    [J]. NUMERISCHE MATHEMATIK, 2004, 98 (02) : 371 - 387
  • [10] Sidi A., 2013, J SCI COMPUT