收稿日期：2020-11-25基金项目：福建省自然科学基金(2016J01032)作者简介：程金发(1966-)，男，江西省乐平市人，博士，教授，博士生导师.*通信作者.E-mail ：***************.cn非一致格子上离散分数阶差分与分数阶和分程金发*（厦门大学数学科学学院福建厦门，361005）摘要：众所周知，一致格子上分数阶和分与分数阶差分的思想概念也是最近几年才兴起的，并且在该邻域得到了很大的发展.但是在非一致格子x ()z =c 1z 2+c 2z +c 3或者x ()z =c 1q z +c 2q -z +c 3上，分数阶和分与分数阶差分的定义是什么，这是一个十分复杂和有趣的问题.本文首次提出非一致格子上分数阶和分与Riemann-Liouville 分数阶差分、Caputo 分数阶差分的定义以及非一致格子上广义Abel 积分方程的求解等基础性结果.关键词：超几何差分方程；非一致格子；分数阶和分；分数阶差分；特殊函数中图分类号：33C45；33D45；26A33；34K37文献标志码：A文章编号：2095-7122（2021）01-0001-013On the fractional sum and fractional difference on nonuniform latticesCHENG Jinfa *(School of Mathematical Sciences,Xiamen University,Xiamen,Fujian 361005,China )Abstract:As is well known,the idea of a fractional sum and difference on uniform lattice is more current,and gets a lot of development in this field.But the definitions of fractional sum and fractional difference of f ()z on nonuniform lattices x ()z =c 1z 2+c 2z +c 3or x ()z =c 1q z +c 2q -z +c 3seem much more complicated andinteresting.In this article,for the first time we propose the definitions of the fractional sum and fractional difference on nonuniform lattices.The solution of the generalized Abel equation is obtained etc.Key words:special function;orthogonal polynomials;adjoint difference equation;difference equation of hy-pergeometric type;nonuniform lattice第34卷第1期2021年3月闽南师范大学学报（自然科学版）Journal of Minnan Normal University （Natural Science ）Vol.34Ｎｏ．1Mar.20211背景回顾及问题提出正如我们在本文序言指出的，分数阶微积分的概念几乎与经典微积分同时起步，可以回溯到Euler 和Leibniz 时期.经过几代数学家的努力，特别是近几十年来，分数阶微积分已经取得了惊人的发展和广阔的应用，有关分数阶微积分的著作层出不穷，例如文献[1-4],但是在一致格子x ()z =z 和x ()z =q z 或者q -z ,z ∈C 上关于离散分数阶微积分的思想，仍然是最近才兴起的.虽然关于一致格子x ()z =z 和x ()z =q z 的离散分数微积分出现和建立相对较晚,但是该领域目前已经做出了大量的工作，且取得了很大的发展[5-8].在最近十年的学术著作中，程金发[9],Goodrich 和Peterson [10]相继出版了两本有关离散分数阶方程理论、离散分数微积分的著作，其中全面系统地介绍了离散分数微积分的基本定义和基本定理，以及最新的参考资料.有关q -分数阶微积分方面的著作可参见Annaby 和Mansour [11].非一致格子的定义回溯到超几何型微分方程[12-13]:σ()z y ′′()z +τ()z y ′()z +λy ()z =0，(1)的逼近，这里σ()z 和τ()z 分别是至多二阶和一阶多项式，λ是常数.Nikiforov 等[14-15]将式(1)推广到如下最一般的复超几何差分方程σˉ[]x ()s ΔΔx ()s -12éëêùûú∇y ()s ∇x ()s +12τˉ[]x ()s éëêùûúΔy ()s Δx ()s +∇y ()s ∇x ()s +λy ()s =0,(2)这里σˉ()x 和τˉ()x 分别是关于x ()s 的至多二阶和一阶多项式，λ是常数，Δy ()s =y ()s +1-y ()s ,∇y ()s =y ()s -y ()s -1，并且x ()s 必须是以下非一致格子.定义1[16-17]两类格子函数x ()s 称之为非一致格子，如果它们满足x ()s =-c 1s 2+-c 2s +-c 3,(3)x ()s =c 1q s +c 2q -s +c 3,(4)这里c i ,-c i 是任意常数，且c 1c 2≠0,-c 1-c 2≠0.当c 1=1,c 2=c 3=0，或c 2=1,c 1=c 3=0或者-c 2=1,-c 1=-c 3=0时，这两种格子函数x ()s ：x ()s =s ，(5)x ()s =q s 或x ()s =q -s(6)称之为一致格子.给定函数F ()s ，定义关于x γ()s 的差分或差商算子为∇γF ()s =∇F ()s ∇x γ()s ,且∇k γF ()z =∇∇x γ()z ()∇∇x γ+1()z ⋯()∇()F ()z ∇x γ+k -1()z .()k =1,2,⋯关于差商算子，命题1是常用的.命题1给定两个复函数f ()s ,g ()s ，成立恒等式Δυ()f ()s g ()s =f ()s +1Δυg ()s +g ()s Δυf ()s =g ()s +1Δυf ()s +f ()s Δυg ()s ,Δυ()f ()s g ()s =g ()s +1Δυf ()s -f ()s +1Δυg ()s g ()s g ()s +1=g ()s Δυf ()s -f ()s Δυg ()s g ()s g ()s +1,Δυ()f ()s g ()s =f ()s -1Δυg ()s +g ()s Δυf ()s =g ()s -1Δυf ()s +f ()s Δυg ()s ,(7)Δυ()f ()s g ()s =g ()s -1Δυf ()s -f ()s -1Δυg ()s g ()s g ()s -1=闽南师范大学学报（自然科学版）2021年2g ()s Δυf ()s -f ()s Δυg ()s g ()s g ()s -1.我们必须指出,在非一致格子式(3)或者式(4),即使当n ∈N ，如何建立非一致格子的n -差商公式，也是一件很不平凡的工作，因为它是十分复杂的，也是难度很大的.事实上,在文献[14-15]中,Nikiforov 等利用插值方法得到了如下n -阶差商∇()n 1[]f ()s 公式：定义2[12-13]对于非一致格子式(3)或式(4)，让n ∈N +，那么∇()n 1[]f ()s =∑k =0n ()-1n -k[]Γ()n +1q[]Γ()k +1q[]Γ()n -k +1q×∏l =0n∇x []s +k -()n -12∇x []s +()k -l +12f ()s -n +k =∑k =0n()-1n -k[]Γ()n +1q[]Γ()k +1q[]Γ()n -k +1q×∏l =0n ∇x n +1()s -k ∇x []s +()n -k -l +12f ()s -k ,(8)这里[]Γ()s q 是修正的q -Gamma 函数，它的定义是[]Γ()s q=q -()s -1()s -24Γq ()s ，并且函数Γq ()s 被称为q -Gamma 函数;它是经典Euler Gamma 函数Γ()s 的推广.其定义是Γq ()s =ìíîïïïï∏k =0∞(1-q k +1)()1-q s -1∏k =0∞(1-q s +k),当||q <1;q -()s -1()s -22Γ1q ()s ,当||q >1.(9)经过进一步化简后，Nikiforov 等在文献[14]中将n 阶差分∇()n 1[]f ()s 的公式重写成下列形式:定义3[14]对于非一致格子式(3)或式(4)，让n ∈N +，那么∇()n 1[]f ()s =∑k =0n ()[]-n qk[]k q ![]Γ()2s -k +c q[]Γ()2s -k +n +1+c qf ()s -k ∇x n +1()s -k ,这里[]μq=γ()μ=ìíîïïïïq u2-q -u 2q 12-q -12如果x ()s =c 1q s +c 2q -s +c 3;μ,如果x ()s =-c 1s 2+-c 2s +c 3,(10)且c =ìíîïïïïïïïïlog c 2c 1log q ,当x ()s =c 1q s +c 2q -s +c 3,-c 2-c 1,当x ()s =-c 1s 2+-c 2s +c 3.程金发：非一致格子上离散分数阶差分与分数阶和分第1期3现在存在两个十分重要且具有挑战性的问题需要进一步深入探讨:1）对于非一致格子上超几何差分方程式（2），在特定条件下存在关于x ()s 多项式形式的解，如果用Rodrigues 公式表示的话，它含有整数阶高阶差商.一个新的问题是：若该特定条件不满足，那么非一致格子上超几何差分方程式（2）的解就不存在关于x ()s 的多项式形式，这样高阶整数阶差商就不再起作用了.此时非一致格子超几何方程的解的表达形式是什么呢？这就需要我们引入一种非一致格子上分数阶差商的新概念和新理论.因此,关于非一致格子上α-阶分数阶差分及α-阶分数阶和分的定义是一个十分有趣和重要的问题.显而易见，它们肯定是比整数高阶差商更为难以处理的困难问题，自专著[14-15]出版以来，Nikiforov 等并没有给出有关α-阶分数阶差分及α-阶分数阶和分的定义，我们能够合理给出非一致格子上分数阶差分与分数阶和分的定义吗?2)另外，我们认为作为非一致格子上最一般性的离散分数微积分，它们也会有独立的意义，并可以导致许多有意义的结果和新理论.本文的目的是探讨非一致格子上离散分数阶和差分.受文章篇幅所限，本文我们仅合理给出非一致格上分数阶和分与分数阶差分的基本定义，其它更多结果例如：非一致格子离散分数阶微积分的一些基本定理，如:Euler Beta 公式，Cauchy Beta 积分公式，Taylor 公式、Leibniz 公式在非一致格子上的模拟形式，非一致格子上广义Abel 方程的解，以及非一致格子上中心分数差分方程的求解，离散分数阶差和分与非一致格子超几何方程之间联系等内容，请参见笔者新专著[16].2非一致格子上的整数和分与整数差分设x ()s 是非一致格子，这里s ∈ℂ.对任意实数γ，x γ()s =x ()s +γ2也是一个非一致格子.让∇γF ()s =f ()s .那么F ()s -F ()s -1=f ()s []x γ()s -x γ()s -1.选取z ,a ∈ℂ，和z -a ∈N .从s =a +1到z ，则有F ()z -F ()a =∑s =a +1zf ()s ∇x r()s .因此，我们定义∫a +1z f ()s d ∇x γ()s =∑s =a +1zf ()s ∇xγ()s .容易直接验证下列式子成立.命题2给定两个复变函数F ()z ,f ()z ，这里复变量z ,a ∈C 以及z -a ∈N ，那么成立1）∇γéëêùûú∫a +1zf ()s d ∇x γ()s =f ()z ;2）∫a +1z∇γF ()s d ∇x γ()s =F ()z -F ()a .现在让我们定义非一致格子上的广义n -阶幂函数[]x ()s -x ()z ()n 为[]x ()s -x ()z ()n =∏k =0n -1[]x ()s -x ()z -k ,()n ∈N +,当n 不是正整数时，需要将广义幂函数加以进一步推广，它的性质和作用是非常重要的，非一致格子上广义幂函数[]x γ()s -x γ()z ()α的定义如下：闽南师范大学学报（自然科学版）2021年4定义4[17-18]设α∈C ，广义幂函数[]x γ()s -x γ()z ()α定义为[]x γ()s -x γ()z ()α=ìíîïïïïïïïïïïïïïïïïïïïïΓ()s -z +a Γ()s -z ,如果x ()s =s ,Γ()s -z +a Γ()s +z +γ+1Γ()s -z Γ()s +z +γ-α+1,如果x ()s =s 2,()q -1αq α()γ-α+12Γq ()s -z +αΓq ()s -z ,如果x ()s =q s ,12α()q -12αq -α()s +γ2Γq ()s -z +αΓq ()s +z +γ+1Γq ()s -z Γq ()s +z +γ-α+1,如果x ()s =q s +q -s 2.(11)对于形如式(4)的二次格子，记c =-c 2-c 1，定义[]x γ()s -x γ()z ()α=-c 1αΓ()s -z +a Γ()s +z +γ+c +1Γ()s -z Γ()s +z +γ-α+c +1；(12)对于形如式(3)的二次格子，记c =logc 2c 1log q，定义[]xγ()s -x γ()z ()α=éëùûc 1()1-q 2αq -α()s +γ2Γq()s -z +a Γq()s +z +γ+c +1Γq()s -z Γq()s +z +γ-α+c +1,(13)这里Γ()s 是Euler Gamma 函数，且Γq ()s 是Euler q -Gamma 函数，其定义如式(9).命题3[17-18]对于x ()s =c 1q s +c 2q -s +c 3或者x ()s =-c 1s 2+-c 2s +-c 3，广义幂[]x γ()s -x γ()z ()α满足下列性质：[]x γ()s -x γ()z []x γ()s -x γ()z -1()μ=[]x γ()s -x γ()z ()μ[]xγ()s -x γ()z -μ=(14)[]xγ()s -x γ()z ()μ+1；(15)[]xγ-1()s +1-x γ-1()z ()μ[]xγ-μ()s -x γ-μ()z =[]x γ-μ()s +μ-x γ-μ()z []x γ-1()s -x γ-1()z ()μ=[]x γ()s -x γ()z ()μ+1；(16)ΔzΔx γ-μ+1()z []xγ()s -x γ()z ()μ=-∇s∇x γ+1()s []x γ+1()s -x γ+1()z ()μ=(17)-[]μq []x γ()s -x γ()z ()μ-1；(18)∇z∇x γ-μ+1()z {}1[]xγ()s -x γ()z ()μ=-ΔsΔx γ-1()s ìíîïïüýþïï1[]x γ-1()s -x γ-1()z ()μ=(19)[]μq[]xγ()s -x γ()z ()μ+1(20)这里[]μq 定义如式（10）.程金发：非一致格子上离散分数阶差分与分数阶和分第1期5现在让我们详细给出非一致格子x γ()s 上整数阶和分的定义，这对于我们进一步给出非一致格子x γ()s 上分数阶和分的定义是十分有帮助的.设γ∈R ，对于非一致格子x γ()s ，数集{}a +1,a +2,⋯,z 中f ()z 的1-阶和分定义为y 1()z =∇-1γf ()z =∫a +1z f ()s d ∇x γ()s ，（21）这里y 1()z =∇-1γf ()z 定义在数集{}a +1,mod ()1中.那么由命题2，我们有∇1γ∇-1γf ()z =∇y 1()z ∇x γ()z =f ()z ，（22）并且对于非一致格子x γ()s ，数集{}a +1,a +2,⋯,z 中f ()z 的2-阶和分定义为y 2()z =∇-2γf ()z =∇-1γ+1[]∇-1γf ()z =∫a +1z y 1()s d ∇x γ+1()s =∫a +1z d ∇x γ+1()s ∫a +1s f ()t d ∇x γ()t =∫a +1z f ()t d ∇x γ()t ∫tz d ∇x γ+1()s =∫a +1z []x γ+1()z -x γ+1()t -1f ()s d ∇x γ()s .(23)这里y 2()z =∇-2γf ()z 定义在数集{}a +1,mod ()1中.同时，可得∇1γ+1∇1γ-1y 1()z =∇y 2()z ∇x γ+1()z =y 1()z ，∇2γ∇-2γf ()z =∇∇x γ()z ()∇y 2()z ∇x γ+1()z =∇y 1()z ∇x γ()z =f ()z .（24）更一般地，由数学归纳法，对于非一致格子x γ()s ，数集{}a +1,a +2,⋯,z 中函数f ()z ，我们可以给出函数f ()z 的n -阶和分定义为y k ()z =∇-kγf ()z =∇-1γ+k -1[]∇-()k -1γf ()z =∫a +1z y k -1()s d ∇x γ+k -1()s =1[]Γ()k q∫a +1z []xγ+k -1()z -x γ+k -1()t -1()k -1f ()t d ∇x γ()t ,()k =1,2,⋯(25)这里[]Γ()k q=ìíîïïq -()k -1()k -2Γq ()k ,如果x ()s =c 1q s +c 2q -s +c 3;Γ()α,如果x ()s =-c 1s 2+-c 2s +c 3,这满足下式[]Γ()k +1q=[]k q []Γ()k q ,[]Γ()2q =[]1q []Γ()1q =1.那么成立∇kγ∇-k γf ()z =∇∇x γ()z ()∇∇x γ+1()z ⋯()∇y k ()z ∇x γ+k -1()z =f ()z .()k =1,2,⋯（26）需要指出的是，当k ∈C 时，式（25）右边仍然是有意义的，因此自然地，我们就可以对非一致格子x γ()s 闽南师范大学学报（自然科学版）2021年6给出函数f ()z 的分数阶和分定义如下：定义5（非一致格子分数阶和分）对任意Re α∈R +，对于非一致格子式（3）和式（4），数集{}a +1,a +2,⋯,z 中的函数f ()z ，我们定义它的α-阶分数阶和分为∇-αγf ()z =1[]Γ()αq∫a +1z []xγ+α-1()z -x γ+α-1()t -1()α-1f ()s d ∇x γ()s ，（27）这里[]Γ()αq=ìíîïïq -()s -1()s -2Γq ()α,如果x ()s =c 1q s +c 2q -s +c 3;Γ()α,如果x ()s =-c 1s 2+-c 2s +c 3,这满足下式[]Γ()α+1q=[]αq []Γ()αq .3非一致格子上的Abel 方程及分数阶差分非一致格子x γ()s 上f ()z 的分数阶差分定义相对似乎更困难和复杂一些.我们的思想是起源于非一致格子上广义Abel 方程的求解.具体来说，一个重要的问题是：让m -1<Re α≤m ，定义在数集{}a +1,a +2,⋯,z 的f ()z 是一给定函数，定义在数集{}a +1,a +2,⋯,z 的g ()z 是一未知函数，它们满足以下广义Abel 方程∇-αγg ()z =∫a +1z []x γ+α-1()z -x γ+α-1()t -1()α-1[]Γ()αqg ()t d ∇x γ()t =f ()t ，（28）怎样求解该广义Abel 方程式（28）？为了求解方程式（28），我们需要利用重要的Euler Beta 公式在非一致格子下的基本模拟.定理1[16]（非一致格子上Euler Beta 公式）对于任何α,β∈C ，那么对非一致格子x ()s ，我们有∫a +1z []x β()z -x β()t -1()β-1[]Γ()βq[]x ()t -x ()αα[]Γ()α+1qd ∇x 1()t =[]x β()z -x β()α()α+β[]Γ()α+β+1q.（29）定理2（Abel 方程的解）设定义在数集{}a +1,mod ()1中的函数f ()z 和函数g ()z 满足∇-αγg ()z =f ()z ,0<m -1<Re α≤m ，那么g ()z =∇m γ∇-m +αγ+αf ()z （30）成立.证明我们仅需证明∇-m γg ()z =∇-m +αγ+αf ()z ，即∇-()m -αγ+αf ()z =∇-()m -αγ+α∇-αγg ()z =∇-m γg ()z .事实上，由定义5可得程金发：非一致格子上离散分数阶差分与分数阶和分第1期7∇-()m -αγ+af ()z =∫a +1z []xγ+m -1()z -x γ+m -1()t -1()m -α-1[]Γ()m -αqf ()t d ∇x γ+α()t =∫a +1z []x γ+m -1()z -x γ+m -1()t -1()m -α-1[]Γ()m -αqd ∇x γ+α()t ⋅∫a +1z []xγ+α-1()t -x γ+α-1()s -1()α-1[]Γ()αqg ()s d ∇x γ()s =∫a +1zg ()s ∇x γ()s ∫sz []xγ+m -1()z -x γ+m -1()t -1()m -α-1[]Γ()m -αq⋅[]xγ+α-1()t -x γ+α-1()s -1()α-1[]Γ()αqd ∇x γ+α()t .在定理1中，将α+1替换成s ；α替换成α-1；β替换成m -α，且将x ()t 替换成x γ+α-1()t ，那么x β()t 替换成x γ+m -1()t ，则我们能够得出下面的等式∫sz []xγ+m -1()z -x γ+m -1()t -1()m -α-1[]Γ()m -αq[]xγ+α-1()t -x γ+α-1()s -1()α-1[]Γ()αqd ∇x γ+α()t =[]xγ+m -1()z -x γ+m -1()s -1()-m -1[]Γ()m q,因此，我们有∇-()m -αγ+af ()z =∫a +1z []x γ+m -1()z -x γ+m -1()s -1()-m -1[]Γ()m qg ()s d ∇x γ()s =∇-mγg ()z ,这样就有∇m γ∇-()m -αγ+a f ()z =∇m γ∇-m γg ()z =g ()z .由定理2得到启示，很自然地我们给出关于f ()z 的Riemann-Liouville 型α-阶()0<m -1<Re α≤m 分数阶差分的定义如下：定义6（Riemann-Liouville 分数阶差分）让m 是超过Re α的最小正整数，对于非一致格子x γ()s ，数集{}α,mod ()1中f ()z 的Riemann-Liouville 型α-阶分数阶差分定义为∇αγf ()z =∇m γ()∇α-mγ+αf ()z .（31）形式上来说，在定义5中，如果α替换成-α，那么式（27）的右边将变为∫a +1z []xγ-α-1()z -x γ-α-1()t -1()-α-1[]Γ()-αqf ()t d ∇x γ()t =∇∇x γ-α()t ()∇∇x γ-α+1()t ⋯()∇∇x γ-α+n -1()t ⋅∫a +1z[]xγ+n -α-1()z -x γ+n -α-1()t -1()n -α-1[]Γ()n -αqf ()t d ∇x γ()t =∇n γ-α∇-n +αγf ()z =∇αγ-αf ()z .（33）闽南师范大学学报（自然科学版）2021年8从式（33），我们也可以得到f ()z 的Riemann-Liouville 型α-阶分数阶差分如下:定义7（Riemann-Liouville 型分数阶差分2）对任意Re α>0，对于非一致格子x γ()s ，数集{}a +1,a +2,⋯,z 中f ()z 的Riemann-Liouville 型α-阶分数阶差分定义为∇αγ-αf ()z =∫a +1z x γ-α-1()z -x γ-α-1()t -1()-α-1[]Γ()-αqf ()t d ∇x γ()t ，（34）将∇γ-α()t 替换成∇γ()t ，那么∇αγf ()z =∫a +1z []x γ-1()z -x γ-1()t -1()-α-1[]Γ()-αqf ()t d ∇x γ+α()t ,（35）这里假定[]Γ()-αq ≠0.4非一致格子上Caputo 型分数阶差分在本节，我们将给出非一致格子上Caputo 型分数阶差分的合理定义.定理3（分部求和公式）给定两个复变函数f (s )，g (s )，那么∫a +1z g (s )∇γf (s )d ∇x γ(s )=f (z )g (z )-f (a )g (a )-∫a +1z f (s -1)∇γg (s )d ∇x γ(s )，这里z ,a ∈C ，且假定z -a ∈N .证明应用命题1，可得g (s )∇γf (s )=∇γ[f (z )g (z )]-f (s -1)∇γg (s )，这样就有g (s )∇r f (s )=∇r [f (z )g (z )]-f (s -1)∇r g (s ).关于变量s ，从a +1到z 求和，那么可得∫a +1z g (s )∇γf (s )d ∇x γ(s )=∫a +1z ∇γ[f (z )g (z )]∇x γ(s )-∫a +1z f (s -1)∇γg (s )d ∇x γ(s )=f (z )g (z )-f (a )g (a )-∫a +1z f (s -1)∇γg (s )d ∇x γ(s ).与非一致格子上Riemann-Liouville 型分数阶差分定义的思想来源一样，对于非一致格子上Caputo 型分数阶差分定义思想，也是受启发于非一致格子上广义Abel 方程式（28）的解.在本文第3节，借助于非一致格子上的Euler Beta 公式，我们已经求出广义Abel 方程∇-αγg (z )=f (z )，0<m -1<α≤m ，是g (z )=∇αγf (z )=∇m γ∇-m +αγ+αf (z ).（36）现在我们将用分部求和公式，给出式（36）的另一种新的表达式.事实上，我们有∇a γf (z )=∇m γ∇-m +aγ+a f (z )=∇mγ∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α-1)[Γ(m -α)]qf (s )d ∇x γ+α(s ).（37）应用恒等式∇(s )[x γ+m -1(z )-x γ+m -1(s )](m -α)∇x γ+α(s )=∇(s )[x γ+m -1(z )-x γ+m -1(s -1)](m -α)∇x γ+α(s -1)=-[m -α]q [x γ+m -1(z )-x γ+m -1(s -1)](m -α-1)，那么以下表达式程金发：非一致格子上离散分数阶差分与分数阶和分第1期9∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α-1)[Γ(m -α)]qf (s )d ∇x γ+α(s )，可被改写成∫a +1zf (s )∇(s ){-[x γ+m -1(z )-x γ+m -1(s )](m -α)[Γ(m -α+1)]q}d ∇s =∫a +1z f (s )∇γ+α-1{-[x γ+m -1(z )-x γ+m -1(s )](m -α)[Γ(m -α+1)]q}d ∇x γ+α-1(s ).应用分部求和公式，可得∫a +1zf (s )∇γ+α-1{-[x γ+m -1(z )-x γ+m -1(s )](m -α)[Γ(m -α+1)]q}d ∇x γ+α-1(s )=f (a )[x γ+m -1(z )-x γ+m -1(a )](m -α)[Γ(m -α+1)]q+∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α)[Γ(m -α+1)]q∇γ+α-1[f (s )]d ∇x γ+α-1(s ).因此，这可导出∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α-1)[Γ(m -α)]q}f (s )d ∇x γ+α(s )=f (a )[x γ+m -1(z )-x γ+m -1(a )](m -α)[Γ(m -α+1)]q+∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α)[Γ(m -α+1)]q∇γ+α-1[f (s )]d ∇x γ+α-1(s ).(38)进一步，考虑∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α)[Γ(m -α+1)]q∇γ+α-1[f (s )]d ∇x γ+α-1(s )，(39)利用恒等式∇(s )[x γ+m -1(z )-x γ+m -1(s )](m -α+1)∇x γ+α-1(s )=∇(s )[x γ+m -1(z )-x γ+m -1(s -1)](m -α+1)∇x γ+α-1(s -1)=-[m -α+1]q [x γ+m -1(z )-x γ+m -1(s -1)](m -α)，表达式(39)能被改写成∫a +1z∇γ+α-1[f (s )]∇(s ){-[x γ+m -1(z )-x γ+m -1(s -1)](m -α+1)[Γ(m -α+2)]q}d ∇s =∫a +1z∇γ+α-1[f (s )]∇γ+α-2{-[x γ+m -1(z )-x γ+m -1(s -1)](m -α+1)[Γ(m -α+2)]q}d ∇x γ+α-2(s ).由分部求和公式，我们有∫a +1z ∇γ+α-1[f (s )]∇γ+α-2{-[x γ+m -1(z )-x γ+m -1(s -1)](m -α+1)[Γ(m -α+2)]q}d ∇x γ+α-2(s )=∇γ+α-1f (a )[x γ+m -1(z )-x γ+m -1(a )](m -α+1)[Γ(m -α+2)]q +∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α+1)[Γ(m -α+2)]q[∇γ+α-2∇γ+α-1]f (s )d ∇x γ+α-2(s )=闽南师范大学学报（自然科学版）2021年10∇γ+α-1f (a )[x γ+m -1(z )-x γ+m -1(a )](m -α+1)[Γ(m -α+2)]q+∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α+1)[Γ(m -α+2)]q∇2γ+α-2f (s )d ∇x γ+α-2(s ).因此，我们得到∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α)[Γ(m -α+1)]q∇γ+α-1[f (s )]d ∇x γ+α-1(s )=∇γ+α-1f (a )[x γ+m -1(z )-x γ+m -1(a )](m -α+1)[Γ(m -α+2)]q+∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α+1)[Γ(m -α+2)]q∇2γ+α-2f (s )d ∇x γ+α-2(s ).(40)同理，用数学归纳法，我们可得∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α+k -1)[Γ(m -α+k )]q∇kγ+α-k [f (s )]d ∇x γ+α-k (s )=∇kγ+α-kf (a )[x γ+m -1(z )-x γ+m -1(a )](m -α+k )[Γ(m -α+k +1)]q+∫a +1z[x γ+m -1(z )-x γ+m -1(s -1)](m -α+k )[Γ(m -α+k +1)]q∇k +1γ+α-(k +1)f (s )d ∇x γ+α-(k +1)(s ).(k =0,1,⋯,m -1)(41)将式（38）,（40）和（41）代入式（37），则有∇αγf ()z =∇m γìíîïïf ()a []x γ+m -1()z -x γ+m -1()a ()m -α[]Γ()m -α+1q +∇γ+α-1f ()a []xγ+m -1()z -x γ+m -1()a ()m -α+1[]Γ()m -α+2q+∇kγ+α-kf ()a []x γ+m -1()z -x γ+m -1()a ()m -α+k []Γ()m -α+k +1q+⋯+∇m -1γ+α-()m -1f ()a []x γ+m -1()z -x γ+m -1()a ()2m -α-1[]Γ()2m -αq+üýþïï∫a +1z []xγ+m -1()z -x γ+m -1()s -1()2m -α-1[]Γ()2m -αq∇m γ+α-mf ()s d ∇x γ+α-m ()s =∇m γ{}∑k =0m -1∇kγ+α-kf ()a []x γ+m -1()z -x γ+m -1()a ()m -α+k []Γ()m -α+k +1q+∇α-2m γ+α-m ∇mγ+α-m f ()z =∑k =0m -1∇kγ+α-kf ()a []x γ-1()z -x γ-1()a ()-α+k []Γ()-α+k +1q+∇α-m γ+α-m ∇mγ+α-m f ()z .总之，我们有下面的程金发：非一致格子上离散分数阶差分与分数阶和分第1期11定理4（广义Abel 方程解2）假设定义在数集{}a +1,a +2,⋯,z 上的函数f ()z 和g ()z 满足∇-αγg ()z =f ()z ,0<m -1<Re α≤m ,那么g ()z =∑k =0m -1∇k γ+α-kf ()a []xγ-1()z -x γ-1()a ()-α+k []Γ()-α+k +1q+∇α-m γ+α-m ∇mγ+α-m f ()z .受到定理4的启示，我们很自然地给出函数f ()z 的α-阶()0<m -1<Re α≤m Caputo 分数阶差分如下：定义8（Caputo 分数阶差分）让m 是Re α超过的最小整数，非一致格子上定义在数集{}a +1,a +2,⋯,z 函数f ()z 的α-阶Caputo 分数阶差分定义为C∇αγf ()z =∇α-m γ+α-m ∇mγ+α-m f ()z .最后，本文再强调指出：对于非一致格子上超几何差分方程式（2），在特定条件下存在关于x ()s 多项式形式的解，如果用Rodrigues 公式表示的话，它含有整数阶高阶差分.一个重要的问题是：若该特定条件不满足，那么非一致格子超几何差分方程的解就不存在关于x ()s 的多项式形式，这样高阶整数阶差分将不再起作用了，这就迫切需要我们引入一种非一致格子上分数阶差分的新概念和新理论.因此，关于非一致格子上阶分数阶差分及阶分数阶和分的定义是一个十分有趣和重要的问题.有关非一致格子超几何差分方程与离散分数阶差和分的联系，更深入的内容参见笔者著作[16]及文献[19-21].(42)(43)参考文献：[1]Kilbas A A,Srivastava H M,Trujillo J J.Theory and applications of fractional differential equations[M].Holland:North-Hol-land Mathatics Studies,Elsevier,2006.[2]Miller S,Ross B.An introduction to the fractional calculus and fractional differential equations[M].NewYork:JohnWiley andSons,1993.[3]Podlubny I.Fractional Differential Equations[M].San Diego,CA:Academic Press,1999.[4]Samko S G,Kilbas A A,Marichev O I.Fractional integrals and derivatives:theory and applications[M].London:Gordon andBreach,1993.[5]Anastassiou G A.Nabla discrete fractional calculus and nalba inequalities[J].Mathematical and Computer Modelling,2010,51:562-571.[6]Atici F M,Eloe P W.Discrete fractional calculus with the nable operator[J].Electronic Journal of Qualitative Theory of Differ-ential Equations,Spec.Ed.I,2009(3):1-12.[7]Atici F M,Eloe P W.Initialvalue problems in discrete fractional calulus[J].Pro.Amer.Math.Soc,2009,137:981-989.[8]Ferreira A C,Torres F M.Fractional h-differences arising from the calculus of variations[J].Appl Anal Discrete Math,2011(5):110-121.[9]程金发.分数阶差分方程理论[M].厦门:厦门大学出版社,2011.[10]Goodrich C,Peterson A C.Discrete fractional discrete fractional discrete fractional calculus[M].Switzerland:Springer Inter-national Publishing,2015.[11]Annaby M H,Mansour Z S.q-Fractional Calculus and Equations[M].NewYork:Springer-Verlag,2012.[12]Andrews G E,Askey R,Roy R.Special functions.Encyclopedia of Mathematics and its Applications[M].Cambridge:Cam-bridge University Press,1999.[13]Wang Z X,Guo D R.Special Functions[M].Singapore:World Scientific Publishing,1989.闽南师范大学学报（自然科学版）2021年12[14]Nikiforov A F,Suslov S K,Uvarov V B.Classical orthogonal polynomials of a discrete variable[M].Berlin:Springer-Verlag,1991.[15]Nikiforov A F,Uvarov V B.Special functions of mathematical physics:a unified introduction with applications[M].Basel:Birkhauser Verlag,1988.[16]程金发.非一致格子超几何方程与分数阶差和分[M].北京:科学出版社,2021.[17]Atakishiyev N M,Suslov S K.Difference hypergeometric functions,in:progress in approximation theory[M].New York:Springer-Verlag,1992:1-35.[18]Suslov S K.On the theory of difference analogues of special functions of hypergeo-metric type[J].Russian Math Surveys,1989,44:227-278.[19]Cheng J F,Jia L K.Generalizations of rodrigues type formulas for hypergeometric difference equations on nonuniform[J].Journal of Difference Equations and Applications,2020,26(4):435-457.[20]Cheng J F,Dai W Z.Adjoint difference equation for a Nikiforov-Uvarov-Suslov difference equation of hypergeometric typeon non-uniform Lattices[J].Ramanujan Journal,2020,53:285-318.[21]Cheng J F.On the complex difference equation of hypergeometric type on non-uniform lattices[J].Acta Mathematical Sinica,English Series,2020,36(5):487–511.[责任编辑：钟国翔]程金发：非一致格子上离散分数阶差分与分数阶和分第1期13。
