Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

数学物理学报, 2019, 39(4): 832-838 doi:

论文

关于差分Riccati方程及时滞微分方程的相关结果

王琼,1, 龙芳,2, 王珺,3

Some Results on Difference Riccati Equations and Delay Differential Equations

Wang Qiong,1, Long Fang,2, Wang Jun,3

通讯作者: 王珺, E-mail: majwang@fudan.edu.cn

收稿日期: 2018-05-14  

基金资助: 国家自然科学基金.  11771090
国家自然科学基金.  11571049
国家自然科学基金.  11271227
国家自然科学基金.  11461070
上海自然科学基金.  17ZR1402900
教育部长江学者和创新团队发展计划项目.  IRT1264
山东大学基本科研业务费专项资金.  2017JC019

Received: 2018-05-14  

Fund supported: the NSFC.  11771090
the NSFC.  11571049
the NSFC.  11271227
the NSFC.  11461070
the Natural Science Foundation of Shanghai.  17ZR1402900
the PCSIRT.  IRT1264
the Fundamental Research Funds of Shandong University.  2017JC019

作者简介 About authors

王琼,E-mail:qiongwangsdu@126.com , E-mail:qiongwangsdu@126.com

龙芳,E-mail:969356135@qq.com , E-mail:969356135@qq.com

摘要

研究了有理系数的差分Riccati方程和常系数的时滞微分方程.当系数满足一定关系时,证明了差分Riccati方程的超越亚纯解具有不小于1的增长级.对于常系数的时滞微分方程,讨论了有理解在z → ∞时的渐近行为.

关键词: 差分Riccati方程 ; 亚纯解 ; 时滞微分方程 ; 有理解

Abstract

We investigate difference Riccati equations with rational coefficients and delay differential equations with constant coefficients. For difference Riccati equations with some relation among coefficients, we prove that every transcendental meromorphic solution is of order no less than one. We also consider the rational solutions for delay differential equations.

Keywords: Difference Riccati equation ; Meromorphic solution ; Delay differential equation ; Rational solution

PDF (307KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

王琼, 龙芳, 王珺. 关于差分Riccati方程及时滞微分方程的相关结果. 数学物理学报[J], 2019, 39(4): 832-838 doi:

Wang Qiong, Long Fang, Wang Jun. Some Results on Difference Riccati Equations and Delay Differential Equations. Acta Mathematica Scientia[J], 2019, 39(4): 832-838 doi:

1 引言与主要结果

本文中,亚纯函数均指复平面上的亚纯函数,并采用Nevanlinna理论中的标准记号和基本结果,见文献[12-13, 15].对于亚纯函数ω, S(r,ω)表示满足S(r,ω)=o(T(r,ω)) (r)的量,允许有一个对数测度有穷的例外集.此外, σ(ω),λ(ω)λ(1/ω)分别表示ω的级,零点收敛指数和极点的收敛指数.

Ablowitz, Halburd和Herbest[1]在2000年首先利用Nevanlinna理论研究了复差分方程的Painlevé性质.随后,利用Nevanlinna理论中对数导数引理的差分对应,人们在复差分方程领域得到了不少研究结果[2, 5-6, 9-11].差分Riccati方程是一类重要的非线性差分方程,它也出现在差分Painlevé方程的研究中,比如文献[7, p197].事实上,文献[7]中的差分Riccati方程具有如下形式

ω(z+1)=A(z)+δω(z)δω(z),
(1.1)

其中A为多项式, δ=±1.A为有理函数且δ=1时, Ishizaki[11]证明了:如果方程(1.1)存在有理解,则其任意超越亚纯解的级一定不小于1/2.显然,方程(1.1)只是下面一般差分Riccati方程的特殊形式

ω(z+1)=a(z)ω(z)+b(z)c(z)ω(z)+d(z),
(1.2)

其中a(z),b(z),c(z),d(z)均为有理函数,满足ac0, adbc0.b0,acd0时,则方程(1.2)退化为Pielou logistic方程.对于Pielou logistic方程,陈宗煊在文献[3]中得到任意有穷级超越亚纯解ω都满足σ(ω)=λ(1/ω)1.此外,他改进了Ishizaki的结论,证明了下面的定理.

定理1.1[4,定理8.2.3]  假设方程(1.2)存在有理解,则有穷级超越亚纯解ω满足σ(ω)=λ(1/ω)1.

陈宗煊还提出了猜想:删除条件"方程(1.2)存在有理解"之后,上面定理的结论依然成立[4, p215].注意到c(z)0时,方程(1.2)退化为一阶多项式系数的线性差分方程,由文献[4,定理5.2.1],可知所有超越亚纯解的级必不小于1.从而,接下来我们仅考虑c(z)0的情况,并不失一般性,假设c(z)1.我们得到:

定理1.2  假设c(z)1, a(z), b(z), d(z)为有理函数满足a(z1)+d(z)0

b(z)+a(z1)a(z)+a(z1)d(z+1)+d(z)d(z+1)0,
(1.3)

则方程(1.2)的任意超越亚纯解ω满足σ(ω)=σ(ω)=λ(1/ω)1.

接下来,我们考虑与微分Painlevé方程有关的时滞微分方程,例如

ω(z)[ω(z+1)ω(z1)]+aω(z)=bω(z),
(1.4)

其中a,b为常数,它连续收敛到第一类Painlevé方程[14]. 2017年, Halburd和Korhonen[8]研究了更一般的方程

ω(z+1)ω(z1)+a(z)ω(z)ω(z)=P(z,ω)Q(z,ω),
(1.5)

其中a(z)为有理函数, P(z,ω)Q(z,ω)为关于ω(z)的非零互素多项式,系数为关于z的有理函数.他们证明了如果方程(1.5)有一个超级小于1的超越亚纯解,则方程右端可退化为简单形式.

很自然地,我们想知道方程(1.5)是否存在有理解,以及有理解在无穷点附近的行为.在P(z,ω)Q(z,ω)均为关于ω的常系数多项式的条件下,我们将讨论方程(1.5)的有理解,并且为了讨论方便,将P(ω),Q(ω)写为

P(ω)=bpωp(z)+bp1ωp1(z)++b0,Q(ω)=cqωq(z)+cq1ωq1(z)++c0,
(1.6)

其中bp0,,b1,b0,cq0,,c1,c0是常数.显然, b0,c0不能同时为零.

定理1.3  假设a(z)为常数, P(z,ω)Q(z,ω)由(1.6)式定义, ω(z)为方程(1.5)的有理解.当z, ω(z)必满足下列性质之一

(1) ω(z)=bp2cqz(1+O(z1)),其中p=q;

(2) ω(z)=d+O(z1),这里d为满足P(d)=0的非零常数;

(3) ω(z)=O(z1)b0=0.

推论1.1  假设a(z)为常数, P(z,ω)Q(z,ω)由(1.6)式定义, p=0q1.那么方程(1.5)不存在有理解.

2 定理1.2的证明

假设ω为方程(1.2)的超越亚纯解.取y=1/ω,则T(r,y)=T(r,ω)+O(1)σ(y)=σ(ω).由方程(1.2)不难看出

y(z+1)=1+d(z)y(z)a(z)+b(z)y(z),
(2.1)

这意味着

P(z,y)=a(z)y(z+1)+b(z)y(z)y(z+1)1d(z)y(z)=0,
(2.2)

P(z,0)10.处理这种方程,我们需要文献[6]中一个引理的特殊情形.

引理2.1[6]  假设ω(z)为差分方程P(z,ω)=0的有穷级超越亚纯解,其中P(z,ω)是关于ω(z)的差分多项式.如果常数a满足P(z,a)0,则

m(r,1ωa)=S(r,ω).

对(2.2)式运用引理2.1,我们可得m(r,1/y)=S(r,y),从而

N(r,ω)=N(r,1/y)=T(r,y)+S(r,y)=T(r,ω)+S(r,ω).
(2.3)

这意味着λ(1/ω)=σ(ω),所以ω必有无穷多个极点.接下来,我们做线性变换ω(z)=A(z)f(z)+B(z),其中A0,B为有理函数.将它代入(1.2)式,我们得到

A(z+1)f(z+1)=˜a(z)f(z)+˜b(z)˜c(z)f(z)+˜d(z),
(2.4)

其中

˜a(z)=a(z)A(z)B(z+1)A(z),˜c(z)=A(z),˜d(z)=B(z)+d(z),
(2.5)

˜b(z)=B(z)a(z)+b(z)B(z+1)B(z)d(z)B(z+1).
(2.6)

˜b(z)0时,由(2.6)式可知B(z)0,并且

B(z+1)=a(z)B(z)+b(z)B(z)+d(z),

从而B是方程(1.2)的一个有理解.与此同时,方程(2.4)将退化为Pielou logistic方程,故有σ(f)=λ(1/f)1.由于ω,f之间存在线性变换,这意味着λ(1/ω)=λ(1/f),σ(ω)=σ(f),我们得到λ(1/ω)=σ(ω)>1,这也给出了定理1.1的另一种证明.

接下来,我们讨论˜b(z)0的情形.为了简便,我们先取A(z)=1,B(z)=a(z1)使得˜a(z)0,则方程(2.4)退化为

f(z+1)=˜b(z)f(z)+˜d(z),
(2.7)

其中

˜b(z)=b(z)a(z)d(z)0,˜d(z)=a(z1)+d(z).
(2.8)

(2.7)式也意味着

P1(z,f)=f(z+1)f(z)+˜d(z)f(z+1)˜b(z)=0

P1(z,0)˜b(z)0.利用引理2.1,不难得到m(r,1/f)=S(r,f),从而N(r,1/f)=T(r,f)+S(r,f),则σ(f)=λ(f),并且f有无穷多个零点.

由于˜b(z)˜d(z)都是有理的,存在R>0使得矩形区域D={z:|Rez|<R,|Imz|<R}包含了˜b(z)˜d(z)的所有零点以及极点.我们将CD分成四个部分

D1={z:RezR},D2={z:RRez},

D3={z:ImzR},D4={z:RImz}.

我们可以在D外部取f(z)一个零点z0,假设其重数为k,然后分两种情形进行讨论.

情形1  a(z1)+d(z)0.由(2.8)式,可知˜d(z)0.那么(2.7)式退化为

f(z+1)=˜b(z)/f(z).
(2.9)

我们就z0所处位置分三种子情形来讨论.

子情形1.1  如果z0D1,则有Rez0R.显然z0+nD1(nN),且z0+n不为˜b的零点和极点.由(2.9)式,可知z0+1f(z)k重极点,且

f(z0+2)=˜b(z0+1)/f(z0+1).

从而z0+2也是f(z)k重零点.通过归纳,容易看出z0+2n(nN)都是f(z)的零点.这意味着,对于任意的r|z0|,我们有

n(r,1f)[r|z0|2]+1r|z0|2,

这里[x]表示为正实数x的整数部分.根据零点收敛指数的定义,我们立刻得到

λ(f)lim sup

因此, \sigma(\omega)=\lambda(1/\omega)=\lambda(f)=\sigma(f)\geq 1.

子情形1.2  如果z_0\in D_3\cup D_4,则|\mbox{Im}\, z_0|\geq R,所以z_0+n(n\in\mathbb{N})不为\tilde{b}(z)的零点和极点.利用类似于情形1.1的证明,我们同样可知z_0+2n(n\in\mathbb{N})f的零点,从而\sigma(\omega)=\lambda(1/\omega)\geq 1也成立.

子情形1.3  如果z_0\in D_2,可知-R\leq \mbox{Re}\, z_0,因此z_0-n(n\in\mathbb{N})不为\tilde{b}(z)的零点和极点.将(2.9)改写为f(z)=\tilde{b}(z)/f(z+1),由此我们知道

f(z-1)=\tilde{b}(z-1)/f(z).
(2.10)

由上式不难看出z_0-1f(z)的极点.继续由(2.10)式可知z_0-2也为f(z)的零点.通过归纳, z_0-2n(n\in\mathbb{N})f的零点.类似于子情形1.1,同理可证\sigma(\omega)=\lambda(1/\omega)\geq 1.

情形2  方程(1.2)的系数满足(1.3)式.假设fz_0附近的局部展开为

f(z)=\delta_k(z-z_0)^k+O(z-z_0)^{k+1},

其中\delta_k为常数.将上式代入(2.7)式,则在z_0点的附近,我们可以得到

f(z+1)=\frac{\tilde{b}(z_0)+O(z-z_0)}{\delta_k(z-z_0)^k+O(z-z_0)^{k+1}+\tilde{d}(z_0)+O(z-z_0)}=\frac{\tilde{b}(z_0)}{\tilde{d}(z_0)}+O(z-z_0),

这意味着f(z_0+1)=\tilde{b}(z_0)/\tilde{d}(z_0).再次利用(2.7)式,可知在z_0的附近,我们有

f(z+2)=\frac{\tilde{b}(z_0+1)+O(z-z_0)}{\tilde{b}(z_0)/\tilde{d}(z_0)+O(z-z_0)+\tilde{d}(z_0+1)+O(z-z_0)}\\~~~~~~~~~~~~~~=\frac{\tilde{b}(z_0+1)\tilde{d}(z_0)+O(z-z_0)}{\tilde{b}(z_0)+\tilde{d}(z_0)\tilde{d}(z_0+1)+O(z-z_0)}.
(2.11)

根据(2.8)式,我们注意到

\tilde{b}(z)+\tilde{d}(z)\tilde{d}(z+1)=b(z)+a(z-1)a(z)+a(z-1)d(z+1)+d(z)d(z+1).

根据(1.3)式,可知\tilde{b}(z)+\tilde{d}(z)\tilde{d}(z+1)\equiv 0.由(2.11)式不难看出z_0+2f(z)的极点.再由(2.7)式,容易看出z_0+3f(z)的零点.利用归纳法,我们知道z_0+3n(n\in\mathbb{N})f(z)的零点.类似情形1,同理可证\sigma(\omega)=\lambda(1/\omega)\geq 1.

3 定理1.3的证明

我们将有理解\omega(z)写为如下形式

\begin{equation}\label{eq-5.1} \omega(z)=\frac{F(z)}{G(z)}+d_0+d_{1}z+d_2z^2+\cdots+d_Nz^N, \end{equation}
(3.1)

其中F, G为关于z的两个互素多项式且\deg F=m<\deg G=n, d_0, d_1, \ldots, d_N是常数.

如果d_N\neq0 (N\geq 1),当z\to\infty时,不难看出

\begin{equation}\label{eq-bu3} \begin{array}{ll} &\omega(z)=d_N z^N(1+o(1)), \, \, \omega'(z)=Nd_Nz^{N-1}(1+o(1)), \\ &\omega(z\pm1)=d_N(z\pm1)^N(1+o(1)), \\ &P(\omega)=b_{p}d^{p}_Nz^{Np}(1+o(1)), \, \, Q(\omega)=c_{q}d^{q}_Nz^{Nq}(1+o(1)).\\ \end{array}\end{equation}
(3.2)

将(3.2)式带入方程(1.5),我们得到

\begin{eqnarray*} c_{q}2Nd_N^{q+1}z^{Nq+N-1}(1+o(1))=b_{p}d_N^{p}z^{Np}(1+o(1)).\end{eqnarray*}

这意味着Nq+N-1=Npc_q 2Nd_N^{q+1}=b_p d_N^p,由此推出N=1, p=q, d_N=b_p/(2c_q),所以

\omega(z)=\frac{F(z)}{G(z)}+d_0+\frac{b_pz}{2c_q}.

接下来,我们考虑d_1=\cdots=d_N=0的情形,并就d_0是否为零区分讨论.首先如果d_0\neq0,则当z\to\infty时, \omega(z)=d_0+o(1),并且

\begin{equation}\begin{array}{ll} \omega'(z)&=o(1), \, \, \omega(z+1)=d_0+o(1), \, \, \omega(z-1)=d_0+o(1), \\ P(\omega)&=P(d_0)+o(1), \, \, Q(\omega)=Q(d_0)+o(1).\\ \end{array}\end{equation}
(3.3)

联立(3.3)和(1.5)式,可得(Q(d_0)+o(1))o(1)=P(d_0)+o(1).这意味着P(d_0)=0.

最后,我们讨论d_0=d_1=\cdots=d_N=0的情形,即

\begin{equation}\label{eq-5.3} \omega(z)=\frac{F(z)}{G(z)}.\end{equation}
(3.4)

将(3.4)式带入方程(1.5)中,我们得到

\begin{equation}\label{eq-5.4} \begin{array}{ll} &Q(\omega)[(\overline{F}\underline{G}-\underline{F}\overline{G})GF+a(F'G-FG')\overline{G}\underline{G}]=P(\omega)\overline{G}\underline{G}GF, \end{array}\end{equation}
(3.5)

其中P(\omega), Q(\omega)由(1.6)式可知

\begin{equation}\label{eq-5.5} P(\omega)=\big(\sum \limits^{p}_{j=0} b_jF^{j}G^{p-j}\big)/G^{p}, \quad Q(\omega)=\big(\sum \limits^{q}_{j=0}c_jF^{j}G^{q-j}\big)/G^{q}.\end{equation}
(3.6)

联立(3.5)和(3.6)式,我们得到

\begin{equation}\label{eq-5.7} C(z)(A(z)+B(z))=D(z), \end{equation}
(3.7)

其中A, B, C, D如下所示

A(z)=GF(\overline{F}\underline{G}-\underline{F}\overline{G}), \, \, B(z)=a(F'G-FG')\overline{G}\underline{G},

C(z)=(\sum \limits^{q}_{j=0}c_jF^{j}G^{q-j})G^p, \, \, D(z)=\overline{G}\underline{G}G^{q+1}F(\sum \limits^{p}_{j=0} b_jF^{j}G^{p-j}).

注意到A+B的最高次项是a(m-n)h_{F, m}h_{G, n}^3z^{3n+m-1}, 这里h_{F, m}, h_{G, n}分别为F, G的最高次项的系数.

我们知道存在b_l\neq 0, c_{k}\neq 0使得b_j=0(j=0, \cdots, l-1), c_t=0(t=0, \cdots, k-1).特别地,如果b_0\not=0c_0\not=0,我们令l=0k=0.由(3.6)式,我们知道

\deg(\sum \limits^{p}_{j=0} b_{j}F^{j}G^{p-j})=ml+n(p-l), \, \, \deg (\sum \limits^{q}_{j=0}c_jF^{j}G^{q-j})=mk+n(q-k).

我们注意到条件a(z)恒为常数,我们设该常数为a.a\not=0,可知\deg(C(A+B))=3n+m-1+np+mk+n(q-k)\deg(D)=n(q+3)+m+ml+n(p-l).由(3.7)式,可知

3n+m-1+np+mk+n(q-k)=n(q+3)+m+ml+n(p-l).

这意味着(n-m)(l-k)=1,从而n=m+1l=k+1.

a=0,我们得到\deg(C(A+B))\leq 3n+m-2+np+mk+n(q-k)\deg(D)=n(q+3)+m+ml+n(p-l).再次由(3.7)式,可知

n(q+3)+m+ml+n(p-l)\leq 3n+m-2+np+mk+n(q-k).

则推出(n-m)(l-k)\geq 2.因此, n\geq m+1, l-k\geq 2同时成立,或者n\geq m+2, l-k\geq 1同时成立.由此l\geq1,也就是b_0=0.所以,当z\to\infty时,我们得到\omega(z)=O(z^{-1}).

参考文献

Ablowitz M J , Halburd R , Herbst B .

On the extension of the Painlevé property to difference equations

Nonlinearity, 2000, 13 (3): 889- 905

DOI:10.1088/0951-7715/13/3/321      [本文引用: 1]

Chen Z X .

Zeros of entire solutions to complex linear difference equations

Acta Mathematica Scientia, 2012, 32 (3): 1141- 1148

DOI:10.1016/S0252-9602(12)60086-1      [本文引用: 1]

Chen Z X .

On growth, zeros and poles of meromorphic solutions of linear and nonlinear difference equations

Science China Math, 2011, 54 (10): 2123- 2133

DOI:10.1007/s11425-011-4265-y      [本文引用: 1]

Chen Z X . Complex Differences and Difference Equations. Beijing: Science Press, 2014

[本文引用: 3]

Chiang Y M , Feng S J .

On the Nevanlinna characteristic of f(z+η) and difference equations in the complex plane

Ramanujan J, 2008, 16 (1): 105- 129

DOI:10.1007/s11139-007-9101-1      [本文引用: 1]

Halburd R G , Korhonen R J .

Difference analogue of the lemma on the logarithmic derivative with applications to difference equations

J Math Anal Appl, 2006, 314 (2): 477- 487

DOI:10.1016/j.jmaa.2005.04.010      [本文引用: 3]

Halburd R G , Korhonen R J .

Existence of finite-order meromorphic solutions as a detector of integrability in difference equations

Physica D, 2006, 218: 191- 203

DOI:10.1016/j.physd.2006.05.005      [本文引用: 2]

Halburd R G , Korhonen R J .

Growth of meromorphic solutions of delay differential equation

Proc Amer Math Soc, 2017, 145 (6): 2513- 2526

DOI:10.1090/proc/13559      [本文引用: 1]

Halburd R G , Korhonen R J .

Nevanlinna theory for the difference operator

Ann Acad Sci Fenn Math, 2006, 31 (2): 463- 478

URL     [本文引用: 1]

Halburd R G , Korhonen R J .

Meromorphic solutions of difference equations, integrability and the discrete Painlevé equations

J Phys, 2007, 40 (6): 1- 38

Ishizaki K .

On difference Riccati equations and second order linear difference equations

Aequationes Mathematicae, 2011, 81 (1): 185- 198

URL     [本文引用: 2]

Laine I. Nevanlinna Theory and Complex Differential Equations. Berlin:Walter de Gruyter, 1993

[本文引用: 1]

Nevanlinna R, Behnke H, Grauert H, et al. Analytic Functions. New York:Springer-Verlag, 1970

[本文引用: 1]

Quispel G R W , Capel H W , Sahadevan R .

Continuous symmetries of differential difference equations:the Kac-van Moerbeke equation and Painlevé reduction

Phys Lett A, 1992, 170 (5): 379- 383

DOI:10.1016/0375-9601(92)90891-O      [本文引用: 1]

Yang C C , Yi H X . Uniqueness Theory of Meromorphic Functions. Dordrecht: Kluwer Academic Publishers, 2004

[本文引用: 1]

/