許小勇， 熊臨晨

(東華理工大學(xué) 理學(xué)院，江西 南昌 330013)

非局部條件主要出現(xiàn)在邊界數(shù)據(jù)無法直接測(cè)量的情況下，具有非局部條件的邊值問題具有非常重要的應(yīng)用，最近具有非局部條件的分?jǐn)?shù)階熱傳導(dǎo)方程得到了大量的關(guān)注。Ashyralyev等(2003)利用半群方法求解了一維經(jīng)典的具有非局部條件的熱傳導(dǎo)方程并分析方法的穩(wěn)定性。Zhou等(2009)提出了基于再生核方法求解經(jīng)典的具有非局部條件的拋物型方程。然而利用數(shù)值方法討論含非局部邊界條件分?jǐn)?shù)階熱傳導(dǎo)方程的文獻(xiàn)很少見報(bào)道。Karatay等(2011)利用基于改進(jìn)的高斯消去法求解了非齊次非局部分?jǐn)?shù)階熱傳導(dǎo)方程；Bhrawy等(2014)應(yīng)用Legendre Tau方法求解了此問題。

小波分析理論產(chǎn)生于1980年代初,其思想來源于Fourier分析,但克服了Fourier分析時(shí)頻不能同時(shí)局部化的缺陷, 是Fourier分析理論和應(yīng)用的發(fā)展，在工程中具有非常重要的應(yīng)用(張曉峰等，2016)。由于小波方法可以同時(shí)具備較強(qiáng)的空間和頻譜分析能力，有利于高效的得到高精度解。在數(shù)學(xué)上由正交多項(xiàng)式通過伸縮、平移得到的正交小波近年來被廣泛應(yīng)用于各類方程的數(shù)值計(jì)算, 如：Legendre小波、Chebyshev小波、Bernoulli小波等(Ali et al.，2017；鄭明等，2017；Rahimkhani et al.，2017)。利用第二類Chebyshev小波求解如下具有非局部條件的分?jǐn)?shù)階熱傳導(dǎo)方程：

(1)

非局部條件和邊值條件

u(x,0)=u(x,1)+f(x),0

(2)

u(0,t)=g0(t),u(1,t)=g1(t),0

(3)

其中0<α<1，u(x,t)是關(guān)于x和t的溫度函數(shù)，q(x,t)為已知源項(xiàng)。本研究主要目的是提出一種高精度的數(shù)值方法求解方程(1)，該方法是基于第二類Chebyshev小波分?jǐn)?shù)階微分和積分公式并結(jié)合小波配置法，將所求解問題轉(zhuǎn)化為代數(shù)方程組求解，給出了數(shù)值例子并與相關(guān)結(jié)果進(jìn)行了比較。

1 第二類Chebyshev小波及性質(zhì)

Chebyshev小波ψn,m(x)=ψ(k,n,m,x)包含4個(gè)參數(shù)，在[0,1)區(qū)間上的定義如下：

(4)

(5)

ψn′,m′(y)

對(duì)上式進(jìn)行截?cái)嗫傻茫?/p>

(6)

本次所提算法中需要用到第二類Chebyshev小波函數(shù)任意分?jǐn)?shù)階積分公式，具體如下：

定理1 (Zhou et al.，2017)定義在區(qū)間[(n-1)/2k-1,n/2k-1]上第二類Chebyshev小波函數(shù)ψn,m(x)的任意a(>0)階積分計(jì)算公式如下：

Iαψn,m(x)=

(7)

2 第二類Chebyshev小波函數(shù)任意分?jǐn)?shù)階微分公式

根據(jù)Caputo分?jǐn)?shù)階導(dǎo)數(shù)的定義，給出第二類Chebyshev小波函數(shù)任意分?jǐn)?shù)階微分公式。Caputo分?jǐn)?shù)階導(dǎo)數(shù)定義(Podlubny, 1999)如下：

定理2 定義在區(qū)間[(n-1)/2k-1,n/2k-1]上第二類Chebyshev小波函數(shù)ψn,m(x)的任意α(0<α<1)階微分計(jì)算公式如下：

Dαψn,m(x)=

(8)

為了描述方便，引入下面的記號(hào)

因此，第二類Chebyshev小波函數(shù)可表示成

所以，Dαψn,m(x)

3 算法描述

(9)

對(duì)方程(9)中變量從0到x進(jìn)行兩次積分,可得

(10)

=(I2Ψ(x))TUΨ(t)

(11)

在方程(11)中令x=1，得

(12)

u(x,t)=(1-x)g0(t)+xg1(t)

+((I2Ψ(x))T-x(I2Ψ(1))T)UΨ(t)

(13)

對(duì)u(x,t)中變量t求α次導(dǎo)數(shù)，得

+((I2Ψ(x))T-x(I2Ψ(1))T)UDαΨ(t)

(14)

在實(shí)際計(jì)算中I2Ψ(1)≈I2Ψ(1-10-15)，根據(jù)配置法，為了計(jì)算方程(9)中的未知系數(shù)矩陣U，選取如下配置點(diǎn)

xi=(2i-1)/(2kM),tj=(2j-1)/(2kM),i=1,2,…,2k-1M,j=1,2,…,2k-1M-1。將表達(dá)式(9)、(14)代入方程(1)中，并考慮上述配置點(diǎn)，得到2k-1M×(2k-1M-1)個(gè)線性方程

(1-xi)Dαg0(tj)+xiDαg1(rj)+((I2Ψ(xi))T-xi(I2Ψ(1))TUDαΨ(tj)-ΨT(xi)UΨ(tj)

=q(xi,tj)

(15)

另外考慮非局部邊值條件，將配置點(diǎn)xii=1,2,…,2k-1M代入方程(2)，可得2k-1M個(gè)方程

u(xi,0)=u(xi,1)+f(xi)

(16)

由此可以得到含2k-1M×2k-1M個(gè)變量的線性方程組，求此方程組可解出系數(shù)矩陣U，將U代入u(x,t)的表達(dá)式(13)中，便可求得方程(1)的近似解。

4 數(shù)值算例

0

非局部條件u(x,0)=u(x,1)-ln(1+x(1-x)),0

方程精確解為u(x,t)=t2ln(1+x(1-x))。

Bhrawy等(2014)應(yīng)用Legendre tau-譜方法求解了此問題，并給出了α=0.45和0.95時(shí)不同參數(shù)N,M情況下的最大絕對(duì)誤差(MAEs)。表1給出了α=0.45和0.95的最大絕對(duì)誤差，并與Bhrawy的結(jié)果進(jìn)行了比較。圖1分別畫出了k=1時(shí),M=10,α=0.45(圖1a)和M=12，α=0.95(圖1b)的整體誤差函數(shù)圖。為了說明所提出方法的收斂性，筆者畫出了當(dāng)k=1,α=0.45,0.95時(shí)最大絕對(duì)誤差的半對(duì)數(shù)曲線圖(圖2)，可以看出誤差呈指數(shù)遞減，這也是期待的譜精度。本研究的方法比Bhrawy介紹的方法結(jié)果更精確。

例2 考慮如下具有非局部條件的時(shí)間分?jǐn)?shù)階熱傳導(dǎo)方程

+4π2t2sin(2πx)，0

非局部條件u(x,0)=u(x,1)-sin(2πx)，0

表1 k=1時(shí)不同M值下最大絕對(duì)誤差值

Bhrawy等(2014)應(yīng)用Legendre tau方法求解了此問題，為了進(jìn)一步說明本次算法比文獻(xiàn)結(jié)果更精確，在表2中比較了兩種算法的最大絕對(duì)誤差值，表明在參數(shù)相同的情況下，本次結(jié)果具有更高的精度。

表2 k=1時(shí)不同M值下最大絕對(duì)誤差值

5 結(jié)論

在Caputo分?jǐn)?shù)階導(dǎo)數(shù)定義下，利用移位的第二類Chebyshev 多項(xiàng)式解析形式，推導(dǎo)了第二類Chebyshev小波函數(shù)的任意α(0<α<1)階微分公式。建立了一種求解具有非局部條件的分?jǐn)?shù)階熱傳導(dǎo)方程有效的小波配置法。通過使用有限數(shù)量的小波基函數(shù)就可以得到較高精度的解。所給全局誤差與多項(xiàng)式次數(shù)的半對(duì)數(shù)曲線圖說明本研究方法具有指數(shù)收斂特點(diǎn)。所得數(shù)值結(jié)果與精確解和Legendre tau方法近似解的比較，說明了本次所提方法的有效性。