本文讨论了非线性差分方程
xn+1 =xn2+ xn-12 的解的渐近性质,给出了零解的收敛域的一个子域以及得到了初始值x
0,x
1 在满足一定的条件下其解发散到无穷大的结论。<br/>
This paper considers a nonlinear difference equation
xn+1 =xn2+ xn-12
with the initial values
邓绍高1,朱立军2,晁越1
1西南交通大学,数学学院,四川 成都
2北方民族大学,数学与信息科学学院,宁夏 银川
收稿日期:2017年10月8日;录用日期:2017年10月24日;发布日期:2017年10月31日
本文讨论了非线性差分方程 x n + 1 = x n 2 + x n − 1 2 的解的渐近性质,给出了零解的收敛域的一个子域以及得到了初始值 x 0 , x 1 在满足一定的条件下其解发散到无穷大的结论。
关键词 :差分方程,平衡点,渐近稳定
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
差分方程在工程技术、经济理论和社会科学中有着十分重要的作用。因此,有许多作者从事差分方程的研究,并取得了很多成果,可参见文献 [
x n + 1 = x n 2 + x n − 1 2 (1.1)
其初始 x 0 , x 1 ∈ ( 0 , + ∞ ) 。
设I为某实数区间,函数 f : I k + 1 → I 连续可微, u = ƒ ( v 0 , v 1 , ⋯ , v k ) ,差分方程:
x n + 1 = ƒ ( x n , x n − 1 , ⋯ , x n − k ) , n ∈ { k , k + 1 , ⋯ } (1.2)
其初始值 x i ∈ I , i ∈ N ( 0 , k ) = { 0 , 1 , ⋯ , k } 。
关于差分方程的平衡点、平衡点的稳定性和渐近稳定性以及渐近平衡点的吸引域的有关概念请参看文献 [
设 β 是方程(1.2)的平衡点,令
p i = f ˙ v i ( β , β , ⋯ , β ) , i ∈ N ( 0 , k )
称方程:
y n + 1 = p 0 y n + p 1 y n − 1 + ⋯ + p k y n − k (1.3)
为方程(1.2)在平衡点 β 处的线性化方程;
称方程:
λ k + 1 − p 0 λ k − p 1 λ k − 1 − ⋯ − p k = 0 (1.4)
为方程(1.3)的特征方程(参看文献 [
引理1.1 [
上述差分方程(1.1)有两个平衡点 β = 0 和 β = 0.5 。在平衡点 β = 0 处,其线性化方程的特征根
λ 1 = λ 2 = 0 ,
所以平衡点 β = 0 是渐近稳定的。在平衡点 β = 0.5 处,其线性化方程有一特征根:
λ = ( 1 + 5 ) / 2 > 1
所以平衡点 β = 0.5 是不稳定的。
本文给出了平衡点 β = 0 的吸引域的子域。即当初始值 x 0 , x 1 在这个子域中时,方程(1.1)的解收敛到零;同时给出了当初始值 x 0 , x 1 在一定的条件下方程(1.1)的解发散到无穷大的结论。
在初始值 x 0 , x 1 满足条件: x 0 2 + x 1 2 ≥ 1 时,方程(1.1)的解发散到无穷。我们在此只给出平衡点 β = 0 的吸引域的一个子域。
定理2.1:设 D = I 2 \ { ( 0.5 , 0.5 ) } ,其中 I = [ 0 , 0.5 ] ,则D是平衡点 β = 0 的吸引域的子域。即 ∀ ( x 0 , x 1 ) ∈ D 方程(1.1)的解 { x n } n = 0 ∞ 满足:
lim n → ∞ x n = 0 .
证明: ∀ ( x 0 , x 1 ) ∈ D ,有 x 2 = x 0 2 + x 1 2 < 0.5 2 + 0.5 2 = 0.5 ,
x 3 = x 2 2 + x 1 2 < 0.5 2 + 0.5 2 = 0.5
令 r = 2 max { x 2 , x 3 } ,则 r < 1 , x 2 ≤ 0.5 r , x 3 ≤ 0.5 r ,
则, x 4 = x 3 2 + x 2 2 ≤ ( 0.5 r ) 2 + ( 0.5 r ) 2 = 0.5 r 2 ≤ 0.5 r ,
x 5 = x 4 2 + x 3 2 ≤ ( 0.5 r ) 2 + ( 0.5 r ) 2 = 0.5 r 2
设 x 2 k ≤ 0.5 r k , x 2 k + 1 ≤ 0.5 r k , k ≥ 1 ,则
x 2 ( k + 1 ) = x 2 k + 1 2 + x 2 k 2 ≤ ( 0.5 r k ) 2 + ( 0.5 r k ) 2 = 0.5 r 2 k ≤ 0.5 r k + 1 ≤ 0.5 r k ,
x 2 ( k + 1 ) + 1 = x 2 ( k + 1 ) 2 + x 2 k + 1 2 ≤ ( 0.5 r k ) 2 + ( 0.5 r k ) 2 = 0.5 r 2 k ≤ 0.5 r k + 1 .
故, ∀ n ∈ N , x 2 n ≤ 0.5 r n , x 2 n + 1 ≤ 0.5 r n 成立。
由此得 lim n → ∞ x 2 n = 0 , lim n → ∞ x ( 2 n + 1 ) = 0 。
从而, lim n → ∞ x n = 0 。证毕。
注记2.1:定理2.1给出了平衡点 β = 0 的吸引域的子域,也暗示了在初始值 x 0 , x 1 满足条件: x 0 2 + x 1 2 < 1 时,方程(1.1)的解有可能收敛到零,也可能发散到无穷大。
定理2.2:当 x 0 ≥ 0.5 , x 1 > 0.5 或 x 0 > 0.5 , x 1 ≥ 0.5 时,方程(1.1)的解 { x n } n = 0 ∞ 发散到无穷大。即:
lim n → ∞ x n = ∞
证明: ,
x 3 = x 2 2 + x 1 2 > 0.5 2 + 0.5 2 = 0.5 ,
令 r = 2 min { x 2 , x 3 } ,
则 r > 1 , x 2 ≥ 0.5 r , x 3 ≥ 0.5 r 。
假设 k ≥ 1 , x 2 k ≥ 0.5 r k , x ( 2 k + 1 ) ≥ 0.5 r k ,
则 x 2 ( k + 1 ) = x 2 k + 1 2 + x 2 k 2 ≥ ( 0.5 r k ) 2 + ( 0.5 r k ) 2 = 0.5 r 2 k ≥ 0.5 r k + 1 ≥ 0.5 r k ,
x 2 ( k + 1 ) + 1 = x 2 ( k + 1 ) 2 + x 2 k + 1 2 ≥ ( 0.5 r k ) 2 + ( 0.5 r k ) 2 ≥ 0.5 r k + 1 .
即 ∀ n ∈ N , n ≥ 1 , x 2 n ≥ 0.5 r n , x ( 2 n + 1 ) ≥ 0.5 r n 。
所以 lim n → ∞ x 2 n = ∞ , lim n → ∞ x ( 2 n + 1 ) = ∞ 。
从而得 lim n → ∞ x n = ∞ 。证毕。
注记2.2:首先,我们可以判断平衡点 β = 0 的吸引域的边界曲线是:
x 1 = x 1 ( x 0 ) , ( 0 ≤ x 0 ≤ c ) ,
且 x 1 是关于 x 0 单调递减的和经过点 ( 0.5 , 0.5 ) ,其中c是某个常数。但要具体地求出c的值是不容易的,若要确定这条边界曲线那就更复杂了。
下面我们来估计一下常数c的取值范围。
设 x 1 = 0 ,则 x 2 = x 1 2 + x 0 2 = x 0 2 ;
x 3 = x 2 2 + x 1 2 = x 0 4 ;
x 4 = x 3 2 + x 2 2 = x 0 8 + x 0 4 < 1 ;
由此得, x 0 < 0.5 5 − 0.5 4 。
所以, c < 0.5 5 − 0.5 4 。
另一方面,由定理2.1知,当 x 1 = 0 , x 0 = 0.5 时, lim n → ∞ x n = 0 。
综上可得, 0.5 ≤ c < 0.5 5 − 0.5 4 。
其次,当 x 0 = 0 时,我们来估计一下 x 1 的取值范围。
x 2 = x 1 2 + x 0 2 = x 1 2 ;
x 3 = x 2 2 + x 1 2 = x 1 4 + x 1 2 < 1 ;
同上可得, 0.5 ≤ x 1 < 0.5 5 − 0.5 。
由于 0.5 5 − 0.5 ≈ 0.7862 , 0.5 5 − 0.5 4 ≈ 0.8867 可知,定理2.1中的吸引域的子域已经是吸引域的大部分了。
国家自然科学基金(61362033);四川省科技厅基础研究计划(2011JYZ002);西南交通大学本科教育教学研究与改革项目(201704010)。
邓绍高,朱立军,晁越. 差分方程xn+1 =xn2+ xn-12解的渐近性质On the Asymptotic Behavior of xn+1 =xn2+ xn-12[J]. 应用数学进展, 2017, 06(07): 892-895. http://dx.doi.org/10.12677/AAM.2017.67107