高阶非线性微分方程非局部边值问题的解法 Solving Higher Order Nonlinear Differential Equation with Nonlocal Boundary Value Problem

doi:10.12677/AAM.2017.68124

Advances in Applied Mathematics
Vol.06 No.08(2017), Article ID:22860,5 pages
10.12677/AAM.2017.68124

Solving Higher Order Nonlinear Differential Equation with Nonlocal Boundary Value Problem

Yongfang Zhou, Lijun Ma, Xiangmei Zhang, Dayong Jin, Guozhong Su

●Abstract

●Full-Text PDF

●Full-Text HTML

●Full-Text ePUB

●Linked References

●How to Cite this Article

School of Science, Hebei University of Technology, Tianjin

Received: Nov. 13^th, 2017; accepted: Nov. 23^rd, 2017; published: Nov. 29^th, 2017

ABSTRACT

This paper discusses the numerical method for the higher order nonlinear differential equation with nonlocal boundary value problem. By constructing the reproducing kernel space which satisfies the nonlocal boundary value conditions, the simple reproducing kernel numerical approximate method is established. Convergence of approximate solution and its derivatives is proved, respectively.

Keywords:Nonlocal Boundary Value Problems, Higher Order Nonlinear Differential Equation, Reproducing Kernel Space

高阶非线性微分方程非局部边值问题的解法

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

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

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

摘要

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

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

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

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

1. 引言

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

${\begin{cases} y^{(n)} (x) = F (x, y (x), y^{'} (x), \dots, y^{(n - 1)} (x)) \\ y^{(i - 1)} (x_{j}) = b_{i j} \\ y (x_{k + 1}) - y (x_{k + 2}) = b_{n} \end{cases}$

的数值解法，其中 $a < x < b$ ， $1 \leq k \leq n - 1$ ， $m_{1} + m_{2} + \dots + m_{k} = n - 1$ ， $b_{i j}$ 是实数， $1 \leq i \leq m_{j}$ ， $1 \leq j \leq k$ ， $a < x_{1} < x_{2} < \dots < x_{k} < x_{k + 1} < x_{k + 2} < b$ 。

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

${\begin{cases} y^{(n)} (x) = F (x, y (x), y^{'} (x), \dots, y^{(n - 1)} (x)) \\ y^{(i - 1)} (x_{j}) = 0 \\ y (x_{k + 1}) - y (x_{k + 2}) = 0 \end{cases}$ (1)

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

2. 再生核空间

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

$W^{n} [a, b]$ 是再生核空间(证明参见文献 [4] )，对任意 $y (x), z (x) \in W^{n} [a, b]$ ，内积和范数分别为 ${〈 y, z 〉}_{W^{n} [a, b]} = \sum_{i = 0}^{n - 1} y^{(i)} (a) z^{(i)} (a) + \int_{a}^{b} y^{(n)} (x) z^{(n)} (x) d x$ ， ${‖ y ‖}_{W^{n} [a, b]} = \sqrt{{〈 y, y 〉}_{W^{n} [a, b]}}$ 。

定义2 $W^{n 1} [a, b] = {y | y \in 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) \in L^{2} [a, b]}$ 。

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

${〈 y, z 〉}_{W [a, b]} = y (a) z (b) + \int_{a}^{b} y^{'} (x) z^{'} (x) d x$ , ${‖ y ‖}_{W [a, b]} = \sqrt{{〈 y, y 〉}_{W [a, b]}}$ .

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

3. 近似解的构造

定义线性算子 $T$ ： $W^{n 1} [a, b] \to W [a, b]$ 。

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

$T y (x) = F (x, y (x), y^{'} (x), \dots y^{(n - 1)} (x))$ , (2)

其中， $y (x) \in W^{n 1} [a, b]$ ，当 $y = y (x) \in W^{n 1} [a, b]$ 时， $F (x, y (x), y^{'} (x), \dots, y^{(n - 1)} (x)) \in W [a, b]$ 。

引理1 $T$ ： $W^{n 1} [a, b] \to W [a, b]$ 是有界线性算子。

设 ${x_{i}}_{i = 1}^{\infty}$ 是 $[a, b]$ 上的稠密子集。

令 $φ_{i} (x) = R (x, x_{i})$ ， $Ψ_{i} (x) = T^{*} φ_{i} (x)$ ，其中， $T^{*}$ 是 $T$ 的共轭算子。

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

证明由

$\begin{matrix} Ψ_{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) \end{matrix}$ ,

有 $Ψ_{i} (x) \in W^{n 1} [a, b]$ 。

令 ${〈 y (x), Ψ_{i} (x) 〉}_{W^{n 1} [a, b]} = 0, i = 1, 2, \dots$ ，其中 $y (x) \in 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, \dots$ .

由 ${x_{i}}_{i = 1}^{\infty}$ 在 $[a, b]$ 上稠密，故 $T y (x) = 0$ 。由 $T^{- 1}$ 的存在性可知 $y (x) \equiv 0$ ，定理得证。

我们将 ${Ψ_{i} (x)}_{i = 1}^{\infty}$ Gram-Schmidt正交化，得到 $W^{n 1} [a, b]$ 上的完全正交系 ${{\bar{Ψ}}_{i} (x)}_{i = 1}^{\infty}$ ， ${\bar{Ψ}}_{i} (x) = \sum_{k = 1}^{i} β_{i k} Ψ_{k} (x), i = 1, 2, \dots$ ，其中， $β_{i k}$ 是正交化系数。

定理1 假设 ${x_{i}}_{i = 1}^{\infty}$ 在 $[a, b]$ 上稠密，如果 $y (x) \in W^{n 1} [a, b]$ 是方程(2)的解，则

$y (x) = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} α_{k} {\bar{Ψ}}_{k} (x)$ , (3)

其中， $α_{k} = F (x_{k}, y (x_{k}), y^{'} (x_{k}), \dots, y^{(n - 1)} (x_{k})), k = 1, 2, \dots$ 。

证明 $y (x) \in W^{n 1} [a, b]$ ， ${{\bar{Ψ}}_{i} (x)}_{i = 1}^{\infty}$ 是 $W^{n 1} [a, b]$ 上的完全正交系，于是

$\begin{matrix} y (x) = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} {〈 y (x), {\bar{Ψ}}_{k} (x) 〉}_{W^{n 1} [a, b]} {\bar{Ψ}}_{k} (x) = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} {〈 y (x), Ψ_{k} (x) 〉}_{W^{n 1} [a, b]} {\bar{Ψ}}_{k} (x) \\ = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} {〈 y (x), T^{*} φ_{k} (x) 〉}_{W^{n 1} [a, b]} {\bar{Ψ}}_{k} (x) = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} {〈 T y (x), φ_{k} (x) 〉}_{W^{n 1} [a, b]} {\bar{Ψ}}_{k} (x) \\ = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} {〈 T y (x), R (x, x_{k}) 〉}_{W^{n 1} [a, b]} {\bar{Ψ}}_{k} (x) = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} T y (x_{k}) {\bar{Ψ}}_{k} (x) \\ = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} F (x_{k}, y (x_{k}), y^{'} (x_{k}), \dots, y^{(n - 1)} (x_{k})) {\bar{Ψ}}_{k} (x) = \sum_{i = 1}^{\infty} \sum_{k = 1}^{i} β_{i k} α_{k} {\bar{Ψ}}_{k} (x) \end{matrix}$ .

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

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

$y_{n} (x) = \sum_{i = 1}^{n} \sum_{k = 1}^{i} β_{i k} α_{k} {\bar{Ψ}}_{k} (x)$ ,(4)

显然， ${‖ y (x) - y_{n} (x) ‖}_{W^{n 1} [a, b]} \to 0, n \to \infty$ 。

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

${‖ y (x) - y_{n} (x) ‖}_{C} \to 0, {‖ y^{(i)} (x) - y_{n}^{(i)} (x) ‖}_{C} \to 0, n \to \infty, i = 1, 2, \dots, n$ .

证明注意到 ${‖ y (x) - y_{n} (x) ‖}_{W^{n 1} [a, b]} \to 0, n \to \infty$ ，于是当 $n \to \infty$ 时，有

$\begin{matrix} | y (x) - y_{n} (x) | = | {〈 y (x^{'}) - y_{n} (x^{'}), K (x, x^{'}) 〉}_{W^{n 1} [a, b]} | \\ \leq {‖ y (x^{'}) - y_{n} (x^{'}) ‖}_{W^{n 1} [a, b]} {‖ K (x, x^{'}) ‖}_{W^{n 1} [a, b]} \\ \leq {\tilde{C}}_{1} {‖ y (x^{'}) - y_{n} (x^{'}) ‖}_{W^{n 1} [a, b]} \to 0 \end{matrix}$ ,

所以 ${‖ y (x) - y_{n} (x) ‖}_{C} \to 0, n \to \infty$

$\begin{matrix} | y^{(i)} (x) - y_{n}^{(i)} (x) | = | \frac{d^{i}}{d x^{i}} (y (x^{'}) - y_{n} (x^{'})) | \\ = | {〈 \frac{d^{i}}{d x^{i}} (y (x^{'}) - y_{n} (x^{'})), K (x, x^{'}) 〉}_{W^{n 1} [a, b]} | \\ = | {〈 y (x^{'}) - y_{n} (x^{'}), \frac{d^{i}}{d x^{i}} K (x, x^{'}) 〉}_{W^{n 1} [a, b]} | \\ \leq {‖ y (x^{'}) - y_{n} (x^{'}) ‖}_{W^{n 1} [a, b]} {‖ \frac{d^{i}}{d x^{i}} K (x, x^{'}) ‖}_{W^{n 1} [a, b]} \\ \leq {\tilde{C}}_{2} {‖ y (x^{'}) - y_{n} (x^{'}) ‖}_{W^{n 1} [a, b]} \to 0 \end{matrix}$

所以 ${‖ y^{(i)} (x) - y_{n}^{(i)} (x) ‖}_{C} \to 0, n \to \infty$ ，这里 ${\tilde{C}}_{1}, {\tilde{C}}_{2}$ 是常数， $i = 1, 2, \dots, n$ 。

4. 数值算法

下面将给出方程(2)的近似解 $y_{n} (x)$ 的求解方法。

令

$f (α_{1}, α_{2}, \dots, α_{n}) = \sum_{i = 1}^{n} {[F (x_{k}, y (x_{k}), y^{'} (x_{k}), \dots y^{(n - 1)} (x_{k})) - α_{k}]}^{2}$ ,

若可以获得 $α_{k} (k = 1, 2, \dots, n)$ ，使得 $f (α_{1}, α_{2}, \dots, α_{n})$ 取得最小值，则可以由式(4)获得方程(2)的近似解 $y_{n} (x)$ 。

下面给出获得 $y_{n} (x)$ 的数值程序。为了获得 $y_{n} (x)$ ，只需要确定 $α_{k} (k = 1, 2, \dots, n)$ 即可。

Step1：取初值 $α_{k 0} (k = 1, 2, \dots, n)$ ；

Step2：计算 $y_{n 0} (x) = \sum_{i = 1}^{n} \sum_{k = 1}^{i} β_{i k} α_{n 0} {\bar{Ψ}}_{k} (x)$ ；

Step3：计算 $f (α_{10}, α_{20}, \dots, α_{n 0})$ ；

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

$α_{k 1} = F (x_{k}, y_{n 1} (x_{k}), {y^{'}}_{n 1} (x_{k}), \dots, y_{n 1}^{(n - 1)} (x_{k})), k = 1, 2, \dots, n$ ;

Step5：计算 $f (α_{11}, α_{21}, \dots, α_{n 1})$ ；

Step6：如果 $f (α_{11}, α_{21}, \dots, α_{n 1}) < f (α_{10}, α_{20}, \dots, α_{n 0})$ ，则用 $α_{k 1}$ 代替 $α_{k 0}$ ，计算

$y_{n 1} (x) = \sum_{i = 1}^{n} \sum_{k = 1}^{i} β_{i k} α_{n 1} {\bar{Ψ}}_{k} (x)$ ,

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

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

$y_{n p} (x) = \sum_{i = 1}^{n} \sum_{k = 1}^{i} β_{i k} α_{n p} {\bar{Ψ}}_{k} (x)$ .(5)

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

参考文献 (References)

1. Henderson, J. and C.J. Kunkel (2008) Uniqueness of Solutions of Linear Nonlocal Boundary Value Problems. Applied Mathematics Letters, 21, 1053-1056.
https://doi.org/10.1016/j.aml.2006.06.024

2. Henderson, J. (2011) Ex-istence and Uniqueness of Solutions of (k+2)-Point Nonlocal Boundary Value Problems for Ordinary Differential Equations. Nonlinear Analysis: Theory, Methods & Applications, 74, 2576-2584.
https://doi.org/10.1016/j.na.2010.11.048

3. Henderson, J. and Luca, R. (2012) Existence and Multiplicity for Positive Solutions of a Multi-Point Boundary Value Problem. Applied Mathematics and Computation, 218, 10572-10585.
https://doi.org/10.1016/j.amc.2012.04.019

4. 周永芳. 若干微分方程非局部边值问题的一种数值方法[D]: [博士学位论文]. 哈尔滨: 哈尔滨工业大学, 2011: 15-19.

期刊菜单