第 38 卷第 1 期 辽宁工程技术大学学报(自然科学版) 2019年 2 月
Vo l . 3 8
No.1 Journal of Liaoning Technical University(Natural Science) Feb. 2019
收稿日期:2018-11-01
基金项目:河北省自然科学基金(A2017203100)
作者简介:陈一鸣(1957-),男,黑龙江 齐齐哈尔人,博士,教授,主要从事分数阶微积分方程数值解及应用方面的研究.
辽宁工程技术大学学报(自然科学版) 网址:http://202.199.224.158/ http://xuebao.lntu.edu.cn/
陈一鸣,张兴军.移位 Chebyshev 多项式数值求解分数阶非线性 Sine–Gordon 方程[J].辽宁工程技术大学学报(自然科学版),
2019,38(1):87-91. doi:10.11956/j.issn.1008-0562.2019.01.015
CHEN Yiming, ZHANG Xingjun.Numerical solution of fractional nonlinear Sine-Gordon equation by shift Chebyshev polynomial
[J].Journal of Liaoning Technical University(Natural Science),2019,38(1):87-91. doi:10.11956/j.issn.1008-0562.2019.01.015
移位 Chebyshev 多项式数值求解分数阶非线性
Sine–Gordon 方程
陈一鸣,张兴军
(燕山大学 理学院,河北 秦皇岛 066004)
摘 要:为解决在物理学中有着广泛应用的一种非线性双曲 Sine-Gordon(SG)方程的数值解问题,提出了移位的
Chebyshev 多项式与分数阶微分性质相结合的高效数值算法.首先,我们推导出移位的 Chebyshev 多项式一阶微分算
子矩阵和分数阶微分算子矩阵,然后将 Sine-Gordon(SG)方程转化为线性代数方程组的形式,进而得到分数阶非
线性 SG 方程的数值解.根据所提出的误差校正相关理论,对数值解进行校正以达到更高的精确度.最后用数值算例
及收敛阶数对算法进行验证,表明了本文所提方法的有效性和实用性.
关键字:移位的 Chebyshev 多项式;分数阶非线性 Sine-Gordon(SG)方程;数值解;微分算子矩阵;误差校正
中图分类号:O 241.8 文献标志码:A 文章编号:1008-0562(2019)01-0092-05
Numerical solution of fractional nonlinear Sine-Gordon equation by shift
Chebyshev polynomial
CHEN Yiming, ZHANG Xingjun
(School of Science, Yanshan University, Qinhuangdao 066004, China)
Abstract: In order to solve the numerical solution problem for a class of nonlinear hyperbolic Sine-Gordon(SG)
equation which has been widely used in physics. An efficient numerical algorithm combining the shifted
Chebyshev polynomial with the fractional differential property is proposed. First, we derive the first-order
differential operator matrix and fractional differential operator of the shifted Chebyshev polynomials. Then, we
transform the Sine-Gordon(SG) equation into a form of linear algebraic equations. And the numerical solution of
fractional nonlinear SG equation is obtained. Moreover, correct the numerical solution to achieve higher accuracy
according to the proposed error correction related theory. Finally, numerical examples and convergence order are
given to show the effectiveness and practicability of the proposed method.
Key words:shifted Chebyshev polynomials; fractional nonlinear Sine-Gordon(SG) equation; numerical solution;
differential operator matrix; error correction
0 引言
随着科学的不断发展,分数阶微分方程(FDE)
广泛应用于各个领域中,其建立的分数阶模型可以有
效的模拟各种物理现象,在工程过程、生物系统以及
诸如反常扩散和非指数松弛模式等方面具有更强的
现实性问题
[1]
.一般来说,分数阶微分方程(FDE)的
解析解不易获得,因此获得这些方程的数值解是至关
重要的.近年来,分数阶微分方程数值方法的研究日益
增多,以便通过数值解近似地预测解析解的趋势.许多
研究人员在基于分数阶微分方程(FDE)特性上,通
过 Block-Pulse 函数法
[2]
,Bernstein 多项式法
[3]
,
Bernoulli 小波法
[4]
、Legendre 多项式法
[5]
、谱方法或
谱配置方法
[6-7]
、同伦分析法
[8]
等方法得到分数阶微分
方程(FDE)的数值解.
Sine-Gordon 方程是一个非线性双曲方程,应用
中国煤炭行业知识服务平台www.chinacaj.net