周永芳，马丽君，张相梅，金大永，苏国忠

河北工业大学理学院，天津

收稿日期：2017年11月13日；录用日期：2017年11月23日；发布日期：2017年11月29日

本文讨论高阶非线性微分方程非局部边值问题的数值方法。通过建立满足非局部边值条件的再生核空间，获得简单易行的再生核数值解法。证明近似解及其导数的收敛性。

关键词 :非局部边值问题，高阶非线性微分方程，再生核空间

Copyright © 2017 by authors and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

微分方程非局部边值问题广泛的出现在振动问题、生命科学、化学反应扩散、人口动力学等科学和工程领域之中。在文献 [ 1 ] [ 2 ] [ 3 ] 中，作者讨论了微分方程非局部边值问题解的存在性和唯一性的问题。本文将讨论如下高阶非线性微分方程非局部边值问题

{ y ( n ) ( x ) = F ( x , y ( x ) , y ′ ( x ) , ⋯ , y ( n − 1 ) ( x ) ) y ( i − 1 ) ( x j ) = b i j y ( x k + 1 ) − y ( x k + 2 ) = b n

的数值解法，其中 a < x < b ， 1 ≤ k ≤ n − 1 ， m 1 + m 2 + ⋯ + m k = n − 1 ， b i j 是实数， 1 ≤ i ≤ m j ， 1 ≤ j ≤ k ， a < x 1 < x 2 < ⋯ < x k < x k + 1 < x k + 2 < b 。

将上面方程中的非局部边值条件齐次化，可以将方程转化成如下的形式

{ y ( n ) ( x ) = F ( x , y ( x ) , y ′ ( x ) , ⋯ , y ( n − 1 ) ( x ) ) y ( i − 1 ) ( x j ) = 0 y ( x k + 1 ) − y ( x k + 2 ) = 0 (1)

通过建立包含非局部边值条件的再生核空间，在空间中构造方程(1)的近似解，证明近似解一致收敛于方程精确解，近似解的导数一致收敛于方程精确解的导数。

定义1 W n [ a , b ] = { y ( x ) | y ( n − 1 ) ( x ) 是绝对连续实值函数， y ( n ) ( x ) ∈ L 2 [ a , b ] } 。

W n [ a , b ] 是再生核空间(证明参见文献 [ 4 ] )，对任意 y ( x ) , z ( x ) ∈ W n [ a , b ] ，内积和范数分别为 〈 y , z 〉 W n [ a , b ] = ∑ i = 0 n − 1 y ( i ) ( a ) z ( i ) ( a ) + ∫ a b y ( n ) ( x ) z ( n ) ( x ) d x ， ‖ y ‖ W n [ a , b ] = 〈 y , y 〉 W n [ a , b ] 。

定义2 W n 1 [ a , b ] = { y | y ∈ W n [ a , b ] , y ( i − 1 ) ( x j ) = 0 , y ( x k + 1 ) − y ( x k + 2 ) = 0 } 。

W n 1 [ a , b ] 是 W n [ a , b ] 的闭子空间(证明参见文献 [ 4 ] )， W n 1 [ a , b ] 是再生核空间。设 W n 1 [ a , b ] 的再生核函数为 K ( x , t ) (具体表达式的确定参见文献 [ 4 ] )。

定义3 W [ a , b ] = { y ( x ) | y ( x ) 是绝对连续实值函数， y ′ ( x ) ∈ L 2 [ a , b ] } 。

W [ a , b ] 是再生核空间，对任意的 y ( x ) , z ( x ) ∈ W [ a , b ] ，内积和范数分别为

〈 y , z 〉 W [ a , b ] = y ( a ) z ( b ) + ∫ a b y ′ ( x ) z ′ ( x ) d x , ‖ y ‖ W [ a , b ] = 〈 y , y 〉 W [ a , b ] .

设 W [ a , b ] 的再生核函数为 R ( x , t ) (具体表达式的确定参见文献 [ 4 ] )。

定义线性算子 T ： W n 1 [ a , b ] → W [ a , b ] 。

对任意 y ( x ) ∈ W n 1 [ a , b ] ，令 T y ( x ) = y ( n ) ( x ) ，则方程(1)转化成如下形式：

T y ( x ) = F ( x , y ( x ) , y ′ ( x ) , ⋯ y ( n − 1 ) ( x ) ) , (2)

其中， y ( x ) ∈ W n 1 [ a , b ] ，当 y = y ( x ) ∈ W n 1 [ a , b ] 时， F ( x , y ( x ) , y ′ ( x ) , ⋯ , y ( n − 1 ) ( x ) ) ∈ W [ a , b ] 。

引理1 T ： W n 1 [ a , b ] → W [ a , b ] 是有界线性算子。

设 { x i } i = 1 ∞ 是 [ a , b ] 上的稠密子集。

令 φ i ( x ) = R ( x , x i ) ， Ψ i ( x ) = T ∗ φ i ( x ) ，其中， T ∗ 是 T 的共轭算子。

引理2 假设 { x i } i = 1 ∞ 在 [ a , b ] 上稠密，则 { Ψ i ( x ) } i = 1 ∞ 是 W n 1 [ a , b ] 上的完全系。

证明 由

Ψ i ( x ) = ( T * φ i ) ( x ) = 〈 ( T * φ i ) ( t ) , K ( t , x ) 〉 W n 1 [ a , b ] = 〈 φ i ( t ) , T K ( t , x ) 〉 W n 1 [ a , b ] = T K ( t , x ) ,

有 Ψ i ( x ) ∈ W n 1 [ a , b ] 。

令 〈 y ( x ) , Ψ i ( x ) 〉 W n 1 [ a , b ] = 0 ,   i = 1 , 2 , ⋯ ，其中 y ( x ) ∈ W n 1 [ a , b ] ，即得

〈 y ( x ) , T * φ i ( x ) 〉 W n 1 [ a , b ] = 〈 T y ( x ) , φ i ( x ) 〉 W n 1 [ a , b ] = T y ( x i ) = 0 ,   i = 1 , 2 , ⋯ .

由 { x i } i = 1 ∞ 在 [ a , b ] 上稠密，故 T y ( x ) = 0 。由 T − 1 的存在性可知 y ( x ) ≡ 0 ，定理得证。

我们将 { Ψ i ( x ) } i = 1 ∞ Gram-Schmidt正交化，得到 W n 1 [ a , b ] 上的完全正交系 { Ψ ¯ i ( x ) } i = 1 ∞ ， Ψ ¯ i ( x ) = ∑ k = 1 i β i k Ψ k ( x ) ,   i = 1 , 2 , ⋯ ，其中， β i k 是正交化系数。

定理1 假设 { x i } i = 1 ∞ 在 [ a , b ] 上稠密，如果 y ( x ) ∈ W n 1 [ a , b ] 是方程(2)的解，则

y ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k α k Ψ ¯ k ( x ) , (3)

其中， α k = F ( x k , y ( x k ) , y ′ ( x k ) , ⋯ , y ( n − 1 ) ( x k ) ) ,   k = 1 , 2 , ⋯ 。

证明 y ( x ) ∈ W n 1 [ a , b ] ， { Ψ ¯ i ( x ) } i = 1 ∞ 是 W n 1 [ a , b ] 上的完全正交系，于是

y ( x ) = ∑ i = 1 ∞ ∑ k = 1 i 〈 y ( x ) , Ψ ¯ k ( x ) 〉 W n 1 [ a , b ] Ψ ¯ k ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k 〈 y ( x ) , Ψ k ( x ) 〉 W n 1 [ a , b ] Ψ ¯ k ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k 〈 y ( x ) , T * φ k ( x ) 〉 W n 1 [ a , b ] Ψ ¯ k ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k 〈 T y ( x ) , φ k ( x ) 〉 W n 1 [ a , b ] Ψ ¯ k ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k 〈 T y ( x ) , R ( x , x k ) 〉 W n 1 [ a , b ] Ψ ¯ k ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k T y ( x k ) Ψ ¯ k ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k F ( x k , y ( x k ) , y ′ ( x k ) , ⋯ , y ( n − 1 ) ( x k ) ) Ψ ¯ k ( x ) = ∑ i = 1 ∞ ∑ k = 1 i β i k α k Ψ ¯ k ( x ) .

定理1给出了方程(2)精确解的表达式。

通过截断式(3)中给定的级数，得到方程(2)的近似解

y n ( x ) = ∑ i = 1 n ∑ k = 1 i β i k α k Ψ ¯ k ( x ) ,(4)

显然， ‖ y ( x ) − y n ( x ) ‖ W n 1 [ a , b ] → 0 ,   n → ∞ 。

定理2 假设方程(2)的解存在唯一， y ( x ) 是方程(2)的解， y n ( x ) 是方程的近似解由式(4)给出，则

‖ y ( x ) − y n ( x ) ‖ C → 0 ,   ‖ y ( i ) ( x ) − y n ( i ) ( x ) ‖ C → 0 ,   n → ∞ ,   i = 1 , 2 , ⋯ , n .

证明 注意到 ‖ y ( x ) − y n ( x ) ‖ W n 1 [ a , b ] → 0 ,   n → ∞ ，于是当 n → ∞ 时，有

| y ( x ) − y n ( x ) | = | 〈 y ( x ′ ) − y n ( x ′ ) , K ( x , x ′ ) 〉 W n 1 [ a , b ] | ≤ ‖ y ( x ′ ) − y n ( x ′ ) ‖ W n 1 [ a , b ] ‖ K ( x , x ′ ) ‖ W n 1 [ a , b ] ≤ C ˜ 1 ‖ y ( x ′ ) − y n ( x ′ ) ‖ W n 1 [ a , b ] → 0 ,

所以 ‖ y ( x ) − y n ( x ) ‖ C → 0 ,   n → ∞

| y ( i ) ( x ) − y n ( i ) ( x ) | = | d i d x i ( y ( x ′ ) − y n ( x ′ ) ) | = | 〈 d i d x i ( y ( x ′ ) − y n ( x ′ ) ) , K ( x , x ′ ) 〉 W n 1 [ a , b ] | = | 〈 y ( x ′ ) − y n ( x ′ ) , d i d x i K ( x , x ′ ) 〉 W n 1 [ a , b ] | ≤ ‖ y ( x ′ ) − y n ( x ′ ) ‖ W n 1 [ a , b ] ‖ d i d x i K ( x , x ′ ) ‖ W n 1 [ a , b ] ≤ C ˜ 2 ‖ y ( x ′ ) − y n ( x ′ ) ‖ W n 1 [ a , b ] → 0

所以 ‖ y ( i ) ( x ) − y n ( i ) ( x ) ‖ C → 0 ,   n → ∞ ，这里 C ˜ 1 , C ˜ 2 是常数， i = 1 , 2 , ⋯ , n 。

下面将给出方程(2)的近似解 y n ( x ) 的求解方法。

令

f ( α 1 , α 2 , ⋯ , α n ) = ∑ i = 1 n [ F ( x k , y ( x k ) , y ′ ( x k ) , ⋯ y ( n − 1 ) ( x k ) ) − α k ] 2 ,

若可以获得 α k ( k = 1 , 2 , ⋯ , n ) ，使得 f ( α 1 , α 2 , ⋯ , α n ) 取得最小值，则可以由式(4)获得方程(2)的近似解 y n ( x ) 。

下面给出获得 y n ( x ) 的数值程序。为了获得 y n ( x ) ，只需要确定 α k ( k = 1 , 2 , ⋯ , n ) 即可。

Step1：取初值 α k 0 ( k = 1 , 2 , ⋯ , n ) ；

Step2：计算 y n 0 ( x ) = ∑ i = 1 n ∑ k = 1 i β i k α n 0 Ψ ¯ k ( x ) ；

Step3：计算 f ( α 10 , α 20 , ⋯ , α n 0 ) ；

Step4：如果 f ( α 10 , α 20 , ⋯ , α n 0 ) < ε ( ε 是充分小的正数)，则计算终止；否则，计算

α k 1 = F ( x k , y n 1 ( x k ) , y ′ n 1 ( x k ) , ⋯ , y n 1 ( n − 1 ) ( x k ) ) ,   k = 1 , 2 , ⋯ , n ;

Step5：计算 f ( α 11 , α 21 , ⋯ , α n 1 ) ；

Step6：如果 f ( α 11 , α 21 , ⋯ , α n 1 ) < f ( α 10 , α 20 , ⋯ , α n 0 ) ，则用 α k 1 代替 α k 0 ，计算

y n 1 ( x ) = ∑ i = 1 n ∑ k = 1 i β i k α n 1 Ψ ¯ k ( x ) ,

否则，放弃 α k 1 。取 α k 0 作为初值， f ( α 1 , α 2 , ⋯ , α n ) 的最小值点即是新的 α k 1 ，用 α k 1 代替 α k 0 并且返回到Step2。

类似上面的步骤，寻找 α k 直到获得 α k p ( k = 1 , 2 , ⋯ , n ) ，使得 f ( α 1 p , α 2 p , ⋯ , α n p ) < ε ，则方程(2)的近似解可以由下式得到

y n p ( x ) = ∑ i = 1 n ∑ k = 1 i β i k α n p Ψ ¯ k ( x ) .(5)

文中通过构造包含方程非局部边值条件的再生核空间，获得了一类高阶非线性微分方程非局部边值问题的精确解和近似解，证明了方程近似解及其导数的一致收敛性。该方法可以进一步推广到其他非线性微分方程非局部边值问题的求解。

河北省自然科学基金(A2015202335)。

周永芳,马丽君,张相梅,金大永,苏国忠. 高阶非线性微分方程非局部边值问题的解法Solving Higher Order Nonlinear Differential Equation with Nonlocal Boundary Value Problem[J]. 应用数学进展, 2017, 06(08): 1034-1038. http://dx.doi.org/10.12677/AAM.2017.68124