Symmetric trigonometrically-fitted two-step hybrid methods for oscillatory problems

被引:6
作者
Li, Jiyong [1 ,2 ]
Wang, Xianfen [1 ,2 ]
Deng, Shuo [3 ]
Wang, Bin [4 ]
机构
[1] Hebei Normal Univ, Coll Math & Informat Sci, Shijiazhuang 050024, Hebei, Peoples R China
[2] Hebei Key Lab Computat Math & Applicat, Shijiazhuang 050024, Hebei, Peoples R China
[3] Hebei GEO Univ, Sch Math & Sci, Shijiazhuang 050031, Hebei, Peoples R China
[4] Qufu Normal Univ, Sch Math Sci, Qufu 273165, Shandong, Peoples R China
关键词
Trigonometrical-fitting; Two-step hybrid methods; Oscillatory problems; Generalized B2-series; Order conditions; Symmetry; RUNGE-KUTTA METHODS; PERTURBED OSCILLATORS; NUMERICAL-INTEGRATION; COLLOCATION METHODS; SYSTEMS; ORDER; Y'';
D O I
10.1016/j.cam.2018.05.038
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In this paper, we propose a general class of trigonometrically-fitted two-step hybrid (TFTSH) methods for solving numerically oscillatory second-order initial value problems. The TFTSH methods integrate exactly the differential system whose solutions can be expressed as the linear combinations of functions from the set {exp(iwt), exp(-iwt)} or equivalently the set {cos(wt), sin(wt)}, where w represents an approximation of the main frequency of the problem. By introducing the generalized B2-series, the necessary and sufficient conditions for TFTSH methods of up to arbitrarily high order p are derived. We also investigate the symmetry of TFTSH methods and analyze the symmetric conditions of TFTSH methods. Based on the order conditions and symmetric conditions, a diagonally implicit two-stage symmetric TFTSH method with order four is constructed. Some numerical experiments are provided to confirm the theoretical expectations. (C) 2018 Elsevier B.V. All rights reserved.
引用
收藏
页码:115 / 131
页数:17
相关论文
共 25 条
[1]   Sixth-order symmetric and symplectic exponentially fitted modified Runge-Kutta methods of Gauss type [J].
Calvo, M. ;
Franco, J. M. ;
Montijano, J. I. ;
Randez, L. .
COMPUTER PHYSICS COMMUNICATIONS, 2008, 178 (10) :732-744
[2]   Order conditions for a class of two-step methods for y"=f(x,y) [J].
Coleman, JP .
IMA JOURNAL OF NUMERICAL ANALYSIS, 2003, 23 (02) :197-220
[3]   Mixed collocation methods for y" = f(x,y) [J].
Coleman, JP ;
Duxbury, SC .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2000, 126 (1-2) :47-75
[4]   Exponentially fitted two-step hybrid methods for y" = f (x, y) [J].
D'Ambrosio, R. ;
Esposito, E. ;
Paternoster, B. .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2011, 235 (16) :4888-4897
[5]   Two-step hybrid collocation methods for y" = f (x, y) [J].
D'Ambrosio, R. ;
Ferro, M. ;
Paternoster, B. .
APPLIED MATHEMATICS LETTERS, 2009, 22 (07) :1076-1080
[6]   Explicit exponentially fitted two-step hybrid methods of high order for second-order oscillatory IVPs [J].
Franco, J. M. ;
Randez, L. .
APPLIED MATHEMATICS AND COMPUTATION, 2016, 273 :493-505
[7]   Exponentially fitted explicit Runge-Kutta-Nystrom methods [J].
Franco, JM .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2004, 167 (01) :1-19
[8]   ON THE PROBLEM OF ENERGY EQUIPARTITION FOR LARGE SYSTEMS OF THE FERMI-PASTA-ULAM TYPE - ANALYTICAL AND NUMERICAL ESTIMATES [J].
GALGANI, L ;
GIORGILLI, A ;
MARTINOLI, A ;
VANZINI, S .
PHYSICA D, 1992, 59 (04) :334-348
[9]  
Gautschi W., 1961, Numer. Math, V3, P381
[10]  
Hairer E., 2006, GEOMETRIC NUMERICAL