数学物理学报, 2024, 44(2): 465-475

一阶非线性时标动态方程的 Hyers-Ulam-Rassias 稳定性

邱仰聪,1, 王其如,2,*

1顺德职业技术学院人文学院 广东佛山 528333

2中山大学数学学院 广州 510275

Hyers-Ulam-Rassias Stability of First-Order Nonlinear Dynamic Equations on Time Scales

Qiu Yangcong,1, Wang Qiru,2,*

1School of Humanities, Shunde Polytechnic, Guangdong Foshan 528333

2School of Mathematics, Sun Yat-sen University, Guangzhou 510275

通讯作者: * 王其如, Email:mcswqr@mail.sysu.edu.cn

收稿日期: 2023-03-22   修回日期: 2023-10-7  

基金资助: 国家自然科学基金(12071491)
广东高校科研项目(重点领域专项)(2021ZDZX4114)

Received: 2023-03-22   Revised: 2023-10-7  

Fund supported: NSFC(12071491)
Special Project in Key Fields of Colleges in Guangdong Province(2021ZDZX4114)

作者简介 About authors

邱仰聪,Email:q840410@qq.com

摘要

利用 Picard 算子和动态不等式, 探讨了一类形式更普遍的一阶非线性时标动态方程的 Hyers-Ulam-Rassias 稳定性, 并且提供三个例子说明这些结论的应用.

关键词: 一阶非线性时标动态方程; Hyers-Ulam-Rassias 稳定性; Picard 算子

Abstract

In this paper, by employing the Picard operator and dynamic inequalities, we investigate Hyers-Ulam-Rassias stability of a class of first-order nonlinear dynamic equations on time scales, which is more general than the equations discussed in the references. Three examples are presented to illustrate the applications of the conclusions.

Keywords: First-order nonlinear dynamic equations on time scales; Hyers-Ulam-Rassias stability; Picard operator

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

本文引用格式

邱仰聪, 王其如. 一阶非线性时标动态方程的 Hyers-Ulam-Rassias 稳定性[J]. 数学物理学报, 2024, 44(2): 465-475

Qiu Yangcong, Wang Qiru. Hyers-Ulam-Rassias Stability of First-Order Nonlinear Dynamic Equations on Time Scales[J]. Acta Mathematica Scientia, 2024, 44(2): 465-475

1 引言

考虑形如以下形式的一类一阶非线性时标动态方程

$ (r(t)x(t))^{\Delta}+p_{1}(t)x^{\sigma}(t)+p_{2}(t)x(t)=f(t,x(t)),\quad t\in \mathbb{I}^{\mathcal{K}}$

的 Hyers-Ulam-Rassias 稳定性, 其中 $ \sigma $ 为前跳算子, $ \mathbb{I}=[a,b]_{\mathbb{T}}=[a,b]\cap \mathbb{T} $, $ \mathbb{I}^{\mathcal{K}}=[a,b)_{\mathbb{T}} $, $ \mathbb{T} $ 是一个时标, $ a,b\in \mathbb{T} $, $ a<b $, 并且假定以下两个条件成立

(A1) $ r, p_{1}, p_{2}\in {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $, 且

$r^{\sigma}(t)+\mu (t)p_{1}(t)>0,\quad r(t)-\mu (t)p_{2}(t)>0,\quad t\in \mathbb{I}^{\mathcal{K}};$

(A2) $ f\in {\rm C_{rd}}(\mathbb{I}\times\mathbb{R}, \mathbb{R}) $, 且存在某常数 $ L>0 $ 使得

$\left\vert f(t,u(t))-f(t,v(t))\right\vert \leq L\cdot (r(t)-\mu (t)p_{2}(t))\left\vert u(t)-v(t)\right\vert.$

1940 年, Ulam[15]首次提出了泛函方程的一种新型稳定性概念 (后来被称为 Ulam 稳定性). 当微分方程或差分方程的解难以精确求出时, Ulam 稳定性可以有助于推断逼近某个近似解的精确解存在性. 接着, Hyers[10]研究了一类线性泛函方程的稳定性, 并在 1941 年成功求出方程的近似解. 后续的一些研究继续改进了文献 [10] 的结论, 其中包含了 1978 年 Rassias 发表在文献 [12] 的成果.

最近, 一些专家学者推导了富有意义的时标动态方程稳定性结论. 关于时间尺度理论和时标动态方程的研究进展, 详情请参见文献 [1],[2],[4],[5],[8],[9],[11]. 2017 年, Shen[14]采用了积分因子法, 在有限区间 $ [a,b]_{\mathbb{T}} $ 上探讨了两个一阶线性时标动态方程

$y^{\Delta}(t)=p(t)y(t)+f(t)$

$x^{\Delta}(t)=-p(t)x^{\sigma}(t)+f(t)$

的 Ulam 稳定性. 2021 年, 利用动态不等式, Alghamdi 等[3]给出了一阶时标动态方程

$\psi^{\Delta}(t)+p(t)\psi(t)=f(t),\quad t\in \mathbb{I}^{\mathcal{K}}$

的 Hyers-Ulam 和 Hyers-Ulam-Rassias 稳定性结果. 2022 年, Bohner 和 Tikare[6]讨论了非线性时标动态方程

$ x^{\Delta}(t)+p(t)x^{\sigma}(t)=f(t,x(t)),\quad t\in \mathbb{I}^{\mathcal{K}}$

的 Ulam 稳定性. 作者使用了 Picard 算子和动态不等式等方法, 分别得到了方程 (1.2) 新的 Ulam 稳定性结论, 以及另外两个时标动态方程

$x^{\Delta}(t)+q(t)x(t)=g(t,x(t)),\quad t\in \mathbb{I}^{\mathcal{K}}$

$x^{\Delta}(t)=F(t,x(t)),\quad t\in \mathbb{I}^{\mathcal{K}}$

的 Hyers-Ulam-Rassias 稳定性结论.

受到文献 [6] 的启发, 考虑在形式上比上述方程更加普遍的 (1.1) 式的 Hyers-Ulam-Rassias 稳定性, 所得结论将可以覆盖、修正和完善文献 [6] 的相关结果. 除此以外, 由于方程 (1.1) 统合了文献 [6] 所提及的三类方程, 因此相关结论可以更便于使用. 最后, 提供三个例子说明相关结论的应用之处.

2 主要结果

在这一节, 先给出以下定义和引理, 作为相关结论的重要理论基础.

定义 2.1 给 $ {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $ 赋予范数

$\|x\|=\sup_{t\in \mathbb{I}}\left\vert x(t)\right\vert.$

这里, 易知 $ ({\rm C_{rd}}(\mathbb{I}, \mathbb{R}), \|\cdot\|) $ 是一个 Banach 空间.

引理 2.1[16,定理2] 假设 $ x,A\in {\rm C_{rd}}(\mathbb{I}, (0,\infty)) $, $ A^{\Delta}(t)\geq 0 $, $ B,C\in {\rm C_{rd}}(\mathbb{I}, [0,\infty)) $ 且均为正回归函数. 对 $ t\in \mathbb{I} $, 若

$x(t)\leq A(t)+\int_{a}^{t}B(s)\left[x(s)+\int_{a}^{s}C(u)x(u)\Delta u\right]\Delta s$

成立, 则有

$x(t)\leq A(t)e_{B+C}(t,a).$

引理 2.2[7,引理3.1] 函数 $ x $ 是方程 (1.2) 的一个解当且仅当

$x(t)=e_{\ominus p}(t,a)x(a)+\int_{a}^{t}e_{\ominus p}(t,s)f(s,x(s))\Delta s,\quad t\in \mathbb{I}$

成立, 这里

$\ominus p=-\frac{p}{1+\mu p}.$

引理 2.3 (抽象 Gronwall 引理, 见文献 [13,引理 2.1]) 设 $ (X,d,\leq) $ 为有序度量空间. 如果 $ A: X\rightarrow X $ 是一个递增 Picard 算子, 且存在唯一的不动点 $ x^{*} $, 那么对所有 $ x\in X $, 若 $ x\leq A(x) $, 则有 $ x\leq x^{*} $; 若 $ x\geq A(x) $, 则有 $ x\geq x^{*} $.

定义 2.2[6,定义2.8] 称方程 (1.2) 具有 $ \mathcal{N} $ 型 Hyers-Ulam-Rassias 稳定性, 这里的 $ \mathcal{N} $ 为 $ \mathbb{I} $ 上所有非减右连续正函数所组成的族, 如果存在常数 $ K>0 $ 使得: 对任意 $ \varepsilon>0 $ 和 $ \varphi \in \mathcal{N} $, 若函数 $ y $ 满足

$\left\vert y^{\Delta}(t)+p(t)y^{\sigma}(t)-f(t,y(t))\right\vert \leq \varepsilon \varphi(t),\quad t\in \mathbb{I}^{\mathcal{K}},$

则方程 (1.2) 存在一个解 $ x $ 满足

$\left\vert y(t)-x(t)\right\vert \leq K\cdot \varepsilon \varphi(t),\quad t\in \mathbb{I}.$

这里, 正数 $ K $ 被称为 HUR 稳定性常数.

接下来, 给出方程 (1.1) 的 Hyers-Ulam-Rassias 稳定性结论.

定理 2.1 假设条件 (A1) 和 (A2) 均成立. 若 $ EL(b-a)<1 $, 这里

$E=\sup_{t,s\in \mathbb{I},t>s}\left\vert e_{\ominus p}(t,s)\right\vert >0,$

则有

(1) 方程 (1.1) 存在唯一解 $ x^{*} $ 满足初始条件 $ x^{*}(a)=x_{0} $, 这里 $ x_{0} $ 为一常数;

(2) 方程 (1.1) 具有 $ \mathcal{N} $ 型 Hyers-Ulam-Rassias 稳定性, 且 HUR 稳定性常数为

$K=E(b-a)e_{EL}(b,a).$

如果 $ x $ 是方程 (1.1) 的一个解, 且初始条件 $ x(a)=x_{0} $ 成立, 这里 $ x_{0} $ 为一常数, 则由方程 (1.1), 有

$r(t)x^{\Delta}(t)+r^{\Delta}(t)x^{\sigma}(t)+p_{1}(t)x^{\sigma}(t)+p_{2}(t)(x^{\sigma}(t)-\mu (t)x^{\Delta}(t))=f(t,x(t)).$

于是, 可推出

$(r(t)-\mu (t)p_{2}(t))x^{\Delta}(t)+(r^{\Delta}(t)+p_{1}(t)+p_{2}(t))x^{\sigma}(t)=f(t,x(t)),$

$x^{\Delta}(t)+\frac{r^{\Delta}(t)+p_{1}(t)+p_{2}(t)}{r(t)-\mu (t)p_{2}(t)}x^{\sigma}(t)=\frac{f(t,x(t))}{r(t)-\mu (t)p_{2}(t)},\quad t\in \mathbb{I}^{\mathcal{K}}.$

$p(t)=\frac{r^{\Delta}(t)+p_{1}(t)+p_{2}(t)}{r(t)-\mu (t)p_{2}(t)},$

则有

$ (\ominus p)(t)=-\frac{p(t)}{1+\mu (t)p(t)}=-\frac{r^{\Delta}(t)+p_{1}(t)+p_{2}(t)}{r^{\sigma}(t)+\mu (t)p_{1}(t)}.$

由引理 2.2, 可得

$ x(t)=e_{\ominus p}(t,a)x_{0}+\int_{a}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,x(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s.$

定义算子 $ U: {\rm C_{rd}}(\mathbb{I}, \mathbb{R})\rightarrow {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $ 如下

$U(x)(t)=e_{\ominus p}(t,a)x_{0}+\int_{a}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,x(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s.$

对任意 $ x,y\in {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $, 可推出

$\begin{align*} \left\vert U(x)(t)-U(y)(t)\right\vert &\leq \int_{a}^{t}e_{\ominus p}(t,s)\cdot \frac{\left\vert f(s,x(s))-f(s,y(s))\right\vert}{r(s)-\mu (s)p_{2}(s)}\Delta s\\ &\leq E\int_{a}^{t}\frac{\left\vert f(s,x(s))-f(s,y(s))\right\vert}{r(s)-\mu (s)p_{2}(s)}\Delta s\\ &\leq EL\int_{a}^{t}\left\vert x(s)-y(s)\right\vert \Delta s\\ &\leq EL(b-a)\|x-y\|, \end{align*}$

即有

$\|U(x)-U(y)\|\leq EL(b-a)\|x-y\|.$

根据 $ EL(b-a)<1 $ 和压缩映射原理, 可知 $ U $ 在 $ {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $ 是收缩的, 且存在唯一的不动点 $ x^{*} $, 这里的 $ x^{*} $ 不但是方程 (1.1) 的唯一解, 而且还满足初始条件 $ x^{*}(a)=x_{0} $.

接着, 假设 $ y $ 是不等式

$ \left\vert y^{\Delta}(t)+\frac{r^{\Delta}(t)+p_{1}(t)+p_{2}(t)}{r(t)-\mu (t)p_{2}(t)}y^{\sigma}(t)-\frac{f(t,y(t))}{r(t)-\mu (t)p_{2}(t)}\right\vert \leq \varepsilon \varphi(t),\quad t\in \mathbb{I}^{\mathcal{K}}$

的一个解, 这里 $ \varepsilon>0 $ 和 $ \varphi \in \mathcal{N} $. 然后, 假定 $ x $ 为方程 (1.1) 的唯一解且满足 $ x(a)=y(a) $. 类似地, 可得 (2.3) 式, 其中的 $ x_{0} $ 被 $ y(a) $ 替代. 由不等式 (5), 易知 $ y $ 满足方程

$y^{\Delta}(t)+\frac{r^{\Delta}(t)+p_{1}(t)+p_{2}(t)}{r(t)-\mu (t)p_{2}(t)}y^{\sigma}(t)=\frac{f(t,y(t))}{r(t)-\mu (t)p_{2}(t)}+\Phi(t),\quad t\in \mathbb{I}^{\mathcal{K}},$

这里 $ \left\vert \Phi(t)\right\vert \leq \varepsilon \varphi(t) $. 于是, 可得

$y(t)=e_{\ominus p}(t,a)y(a)+\int_{a}^{t}e_{\ominus p}(t,s)\cdot \left[\frac{f(s,y(s))}{r(s)-\mu (s)p_{2}(s)}+\Phi(s)\right]\Delta s,\quad t\in \mathbb{I}.$

由 $ \varphi \in \mathcal{N} $, 可知对 $ t\in \mathbb{I}^{\mathcal{K}} $, 有 $ \varphi(t)>0 $ 和 $ \varphi^{\Delta}(t)\geq 0 $, 所以对一切 $ t\in \mathbb{I} $, 均有

$\begin{align*} &\left\vert y(t)-e_{\ominus p}(t,a)y(a)-\int_{a}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,y(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s\right\vert\\ \leq &\int_{a}^{t}e_{\ominus p}(t,s)\cdot \left\vert \Phi(s)\right\vert \Delta s\leq E\varepsilon \int_{a}^{t}\varphi(s)\Delta s\leq E\varepsilon (b-a)\varphi(t) \end{align*}$

$ \begin{matrix} \left\vert y(t)-x(t)\right\vert &=\left\vert y(t)-e_{\ominus p}(t,a)y(a)-\int_{a}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,x(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s\right\vert \\ &\leq \left\vert y(t)-e_{\ominus p}(t,a)y(a)-\int_{a}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,y(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s\right\vert \\ & +\left\vert \int_{a}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,y(s))-f(s,x(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s\right\vert \\ &\leq E\varepsilon (b-a)\varphi(t)+L\int_{a}^{t}e_{\ominus p}(t,s)\cdot \left\vert y(s)-x(s)\right\vert\Delta s. \end{matrix}$

定义算子 $ V: {\rm C_{rd}}(\mathbb{I}, \mathbb{R})\rightarrow {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $ 如下

$V(x)(t)=E\varepsilon (b-a)\varphi(t)+L\int_{a}^{t}e_{\ominus p}(t,s)x(s)\Delta s,\quad t\in \mathbb{I}.$

对任意 $ x,y\in {\rm C_{rd}}(\mathbb{I},\mathbb{R}) $ 和 $ t\in \mathbb{I} $, 可推出

$\begin{align*} \left\vert V(x)(t)-V(y)(t)\right\vert &=L\left\vert \int_{a}^{t}e_{\ominus p}(t,s)\cdot (x(s)-y(s))\Delta s\right\vert\\ &\leq L\int_{a}^{t}e_{\ominus p}(t,s)\cdot \left\vert x(s)-y(s)\right\vert \Delta s\\ &\leq EL\int_{a}^{t}\left\vert x(s)-y(s)\right\vert\Delta s\\ &\leq EL(b-a)\|x-y\|. \end{align*}$

因为 $ EL(b-a)<1 $, 可知 $ V $ 在 $ {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $ 上是收缩的. 更进一步地, 由 Banach 压缩原理, 可推出 $ V $ 是一个 Picard 算子, 而且存在不动点 $ w^{*} $. 所以, 可得

$w^{*}(t)=E\varepsilon (b-a)\varphi(t)+L\int_{a}^{t}e_{\ominus p}(t,s)w^{*}(s)\Delta s,\quad t\in \mathbb{I}.$

注意到 $ w^{*} $ 是递增的, 且

$w^{*}(t)\leq E\varepsilon (b-a)\varphi(t)+\int_{a}^{t}ELw^{*}(s)\Delta s,\quad t\in \mathbb{I},$

由引理 2.1, 取 $ A(t)=E\varepsilon (b-a)\varphi(t) $, $ B(t)=EL $, $ C(t)=0 $, 有

$w^{*}(t)\leq E\varepsilon (b-a)\varphi(t)e_{EL}(t,a),\quad t\in \mathbb{I}.$

令 $ w(t)=\left\vert y(t)-x(t)\right\vert $, 由 (2.4) 式, 不难看出对一切 $ t\in \mathbb{I} $, 均有 $ w(t)\leq V(w)(t) $, 所以 $ V $ 是 $ {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $ 上的一个递增 Picard 算子. 于是, 利用引理 2.3, 可推导出对所有 $ t\in \mathbb{I} $, $ w(t)\leq w^{*}(t) $ 均成立. 这意味着

$\left\vert y(t)-x(t)\right\vert =w(t)\leq E\varepsilon (b-a)\varphi(t)e_{EL}(t,a)\leq E\varepsilon (b-a)\varphi(t)e_{EL}(b,a),\quad t\in \mathbb{I}.$

因此, 方程 (1.1) 具有 $ \mathcal{N} $ 型 Hyers-Ulam-Rassias 稳定性, 且 HUR 稳定性常数为 $ K=E(b-a)e_{EL}(b,a) $. 证毕.

注2.1 在定义 2.2 中, 如果对 $ t\in \mathbb{I} $, 取 $ \varepsilon=1 $ 或 $ \varphi(t)\equiv 1 $, 那么定理 2.1 中的证明过程可以分别得到简化.

注2.2 对所有 $ t\in \mathbb{I} $, 令 $ r(t)\equiv 1 $, 当分别取 $ p_{2}(t)\equiv 0 $, $ p_{1}(t)\equiv 0 $, $ p_{1}(t)=p_{2}(t)\equiv 0 $ 时, 定理 2.1 的结论将会先后与文献 [6,定理 3.1-3.3] 保持一致.

注2.3 当 $ \mathbb{I} $ 为连续时标 $ [a,b]_{\mathbb{R}} $ 时, 定理 2.1 中的 (2.1) 式可简化为

$ (\ominus p)(t)=-\frac{r^{\Delta}(t)+p_{1}(t)+p_{2}(t)}{r^{\sigma}(t)+\mu (t)p_{1}(t)}=-\frac{r'(t)+p_{1}(t)+p_{2}(t)}{r(t)},$

且 HUR 稳定性常数 $ K=E(b-a)e_{EL}(b,a)=E(b-a){\rm e}^{EL(b-a)} $. 注意到附加条件 $ EL(b-a)<1 $, 可知 $ K<E(b-a)$e.

3 例

在这一节, 给出以下三个例子说明相关结论的意义.

例3.1 令 $ \mathbb{T}=\mathbb{N} $, $ a=2 $, $ b=10 $, $ \mathbb{I}=[2,10]_{\mathbb{T}}=\{2,3,4,5,6,7,8,9,10\} $. 对 $ t\in \mathbb{I}^{\mathcal{K}} $, 考虑方程

$ (tx(t))^{\Delta}-x^{\sigma}(t)+(t-1)x(t)=\frac{1}{50}t\cdot \sqrt[3]{x^{3}(t)+1},$

满足初始条件 $ x(2)=1 $, 以及对应的不等式

$ \left\vert (ty(t))^{\Delta}-y^{\sigma}(t)+(t-1)y(t)-\frac{1}{50}t\cdot \sqrt[3]{y^{3}(t)+1}\right\vert \leq \varepsilon \varphi(t).$

这里, $ r(t)=t $, $ p_{1}(t)=-1 $, $ p_{2}(t)=t-1 $, $ f(t,x)=t\cdot \sqrt[3]{x^{3}+1}/50 $, $ \varepsilon>0 $, $ \varphi \in \mathcal{N} $, $ \mu (t)=1 $.

由于

$r^{\sigma}(t)+\mu (t)p_{1}(t)=t+1-1=t>0,\quad r(t)-\mu (t)p_{2}(t)=t-(t-1)=1>0,\quad t\in \mathbb{I}^{\mathcal{K}},$

以及

$\begin{align*} \left\vert f(t,u(t))-f(t,v(t))\right\vert &=\frac{1}{50}t\cdot \left\vert \sqrt[3]{u^{3}(t)+1}-\sqrt[3]{v^{3}(t)+1}\right\vert\\ &\leq \frac{1}{50}\cdot 10\left\vert u(t)-v(t)\right\vert =\frac{1}{5}\left\vert u(t)-v(t)\right\vert,\quad t\in \mathbb{I}, \end{align*}$

令 $ L=1/5 $, 不难看出方程 (3.1) 的系数函数满足条件 (A1) 和 (A2). 由定理 2.1 的证明过程, 可得 (2.1) 式. 于是, 对一切 $ t\in \mathbb{I} $, $ s\in \mathbb{I}^{\mathcal{K}} $ 且 $ t>s $, 有

$\begin{align*} e_{\ominus p}(t,s)&=\exp\left(\int_{s}^{t}\frac{1}{\mu (u)}\cdot \ln\left(1-\mu (u)\cdot \frac{r^{\Delta}(u)+p_{1}(u)+p_{2}(u)}{r^{\sigma}(u)+\mu (u)p_{1}(u)}\right)\Delta u\right)\\ &=\exp\left(\int_{s}^{t}\ln \frac{r(u)-\mu (u)p_{2}(u)}{r^{\sigma}(u)+\mu (u)p_{1}(u)}\Delta u\right)\\ &=\exp\left(\int_{s}^{t}\ln \frac{1}{u}\Delta u\right)=\exp\left(\sum_{u=s}^{t-1}\ln \frac{1}{u}\right)\\ &=\prod_{u=s}^{t-1}\frac{1}{u}=\frac{(s-1)!}{(t-1)!}. \end{align*}$

当 $ t=s\in \mathbb{I}^{\mathcal{K}} $ 时, 有 $ e_{\ominus p}(t,s)=e_{\ominus p}(t,t)=1 $, 与以上结论保持一致. 所以, 可推出

$E=\sup_{t,s\in \mathbb{I},t>s}\frac{(s-1)!}{(t-1)!}=\frac{1}{2}>0,$

于是

$EL(b-a)=\frac{1}{2}\cdot \frac{1}{5}\cdot (10-2)=\frac{4}{5}<1.$

因此, 可推断出方程 (3.1) 存在唯一解 $ x $ 满足 $ x(2)=1 $. 接着, 由不等式 (3.2), 存在函数 $ \Phi \in {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $, 使得 $ \left\vert \Phi(t)\right\vert \leq \varepsilon \varphi(t) $ 和

$(ty(t))^{\Delta}-y^{\sigma}(t)+(t-1)y(t)=\frac{1}{50}t\cdot \sqrt[3]{y^{3}(t)+1}+\Phi(t),\quad t\in \mathbb{I}^{\mathcal{K}}.$

由定理 2.1 的证明过程, 对 $ t\in \mathbb{I} $, 可得

$\begin{align*} y(t)&=e_{\ominus p}(t,2)y(2)+\int_{2}^{t}e_{\ominus p}(t,s)\cdot \left[\frac{f(s,y(s))}{r(s)-\mu (s)p_{2}(s)}+\Phi(s)\right]\Delta s\\ &=\frac{y(2)}{(t-1)!}+\frac{1}{50}\sum_{s=2}^{t-1}\frac{(s-1)!}{(t-1)!}\cdot s\cdot \sqrt[3]{y^{3}(s)+1}+\sum_{s=2}^{t-1}\frac{(s-1)!}{(t-1)!}\cdot \Phi(s)\\ &=\frac{1}{(t-1)!}\left[y(2)+\frac{1}{50}\sum_{s=2}^{t-1}(s!)\cdot \sqrt[3]{y^{3}(s)+1}+\sum_{s=2}^{t-1}((s-1)!)\cdot \Phi(s)\right]. \end{align*}$

于是, 可知方程 (3.1) 存在唯一解 $ x $ 满足 $ x(2)=y(2) $. 由引理 2.2, 可知

$\begin{align*} x(t)&=e_{\ominus p}(t,2)y(2)+\int_{2}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,x(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s\\ &=\frac{1}{(t-1)!}\left[y(2)+\frac{1}{50}\sum_{s=2}^{t-1}(s!)\cdot \sqrt[3]{x^{3}(s)+1}\right],\quad t\in \mathbb{I}. \end{align*}$

所以, 对 $ t\in \mathbb{I} $, 可推出

$\begin{align*} \left\vert y(t)-x(t)\right\vert &\leq \frac{1}{50}\sum_{s=2}^{t-1}\frac{s!}{(t-1)!}\cdot \left\vert \sqrt[3]{y^{3}(s)+1}-\sqrt[3]{x^{3}(s)+1}\right\vert +\sum_{s=2}^{t-1}\frac{(s-1)!}{(t-1)!}\cdot \left\vert \Phi(s)\right\vert\\ &\leq \frac{1}{50}\sum_{s=2}^{t-1}\frac{s!}{(t-1)!}\cdot \left\vert y(s)-x(s)\right\vert +\varepsilon \sum_{s=2}^{t-1}\frac{(s-1)!}{(t-1)!}\cdot \varphi(s)\\ &\leq \frac{1}{50}\sum_{s=2}^{t-1}\frac{s!}{(t-1)!}\cdot \left\vert y(s)-x(s)\right\vert +\frac{1}{2}\varepsilon \cdot 8 \varphi(t)\\ &=\frac{1}{50}\sum_{s=2}^{t-1}\frac{s!}{(t-1)!}\cdot \left\vert y(s)-x(s)\right\vert +4\varepsilon \varphi(t). \end{align*}$

于是, 采用定理 2.1 的结果, 可知方程 (3.1) 具有 $ \mathcal{N} $ 型 Hyers-Ulam-Rassias 稳定性, 且 HUR 稳定性常数为

$K=\frac{1}{2}\cdot 8\cdot e_{\frac{1}{10}}(10,2)=4\cdot (1.1)^{8}\approx 8.5744.$

例3.2 令 $ \mathbb{T}=3^{\mathbb{N}_{0}} $, $ a=1 $, $ b=81 $, $ \mathbb{I}=[1,81]_{\mathbb{T}}=\{1,3,9,27,81\} $. 对 $ t\in \mathbb{I}^{\mathcal{K}} $, 考虑方程

$ \left(\frac{1}{t}\cdot x(t)\right)^{\Delta}+tx^{\sigma}(t)=\frac{\alpha \left(x^{\beta}(t)+1\right)^{1/\beta}}{t},$

满足初始条件 $ x(1)=2 $, 以及对应的不等式

$ \left\vert \left(\frac{1}{t}\cdot y(t)\right)^{\Delta}+ty^{\sigma}(t)-\frac{\alpha \left(y^{\beta}(t)+1\right)^{1/\beta}}{t}\right\vert\leq \varepsilon \varphi(t).$

这里, $ r(t)=1/t $, $ p_{1}(t)=t $, $ p_{2}(t)=0 $, $ f(t,x)=\alpha \left(x^{\beta}(t)+1\right)^{1/\beta}/t $, $ \alpha>0 $, $ \beta \geq 1 $, $ \varepsilon>0 $, $ \varphi \in \mathcal{N} $, $ \mu (t)=2t $.

由于

$r^{\sigma}(t)+\mu (t)p_{1}(t)=\frac{1}{3t}+2t^{2}>0,\quad r(t)-\mu (t)p_{2}(t)=\frac{1}{t}>0,\quad t\in \mathbb{I}^{\mathcal{K}},$

以及

$\begin{align*} \left\vert f(t,u(t))-f(t,v(t))\right\vert &=\alpha \cdot \frac{1}{t}\cdot \left\vert (u^{\beta}(t)+1)^{1/\beta}-(v^{\beta}(t)+1)^{1/\beta}\right\vert\\ &\leq \alpha \cdot \frac{1}{t}\cdot \left\vert u(t)-v(t)\right\vert,\quad t\in \mathbb{I}, \end{align*}$

可以取 $ L=\alpha $, 则不难看出方程 (3.3) 的系数函数满足条件 (A1) 和 (A2). 类似于例 3.1 的分析过程, 对一切 $ t\in \mathbb{I} $, $ s\in \mathbb{I}^{\mathcal{K}} $ 且 $ t>s $, 均有

$\begin{align*} e_{\ominus p}(t,s)&=\exp\left(\int_{s}^{t}\ln \frac{r(u)-\mu (u)p_{2}(u)}{r^{\sigma}(u)+\mu (u)p_{1}(u)}\Delta u\right)\\ &=\exp\left(\int_{s}^{t}\frac{1}{2u}\ln \frac{\frac{1}{u}}{\frac{1}{3u}+2u^{2}}\Delta u\right)=\exp\left(\int_{s}^{t}\frac{1}{2u}\ln \frac{3}{1+6u^{3}}\Delta u\right)\\ &=\exp\left(\sum_{u=\log_{3}s}^{\log_{3}t-1}\ln \frac{3}{1+6\cdot 27^{u}}\right)=\prod_{u=\log_{3}s}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}. \end{align*}$

当 $ t=s\in \mathbb{I}^{\mathcal{K}} $ 时, 同样有 $ e_{\ominus p}(t,s)=e_{\ominus p}(t,t)=1 $. 于是, 可得

$E=\sup_{t,s\in \mathbb{I},t>s}\left(\prod_{u=\log_{3}s}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)=\frac{3}{7}>0.$

$0<\alpha<\frac{1}{E(b-a)}=\frac{1}{\frac{3}{7}\cdot (81-1)}=\frac{7}{240}$

时, 有 $ EL(b-a)=E\alpha (b-a)<1 $, 这可推导出方程 (3.3) 存在唯一解 $ x $ 满足 $ x(1)=2 $.

接着, 由不等式 (3.4), 存在函数 $ \Phi \in {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $ 使得 $ \left\vert \Phi(t)\right\vert \leq \varepsilon \varphi(t) $ 和

$\left(\frac{1}{t}\cdot y(t)\right)^{\Delta}+ty^{\sigma}(t)=\frac{\alpha \left(y^{\beta}(t)+1\right)^{1/\beta}}{t}+\Phi(t),\quad t\in \mathbb{I}^{\mathcal{K}}.$

类似地, 可得

$\begin{align*} y(t)=\ &e_{\ominus p}(t,1)y(1)+\int_{1}^{t}e_{\ominus p}(t,s)\cdot \left[\frac{f(s,y(s))}{r(s)-\mu (s)p_{2}(s)}+\Phi(s)\right]\Delta s\\ =\ &y(1)\cdot \prod_{u=0}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}+\alpha\int_{1}^{t}\left(\prod_{u=\log_{3}s}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot \frac{\left(y^{\beta}(s)+1\right)^{1/\beta}/s}{1/s}\Delta s\\ &+\int_{1}^{t}\left(\prod_{u=\log_{3}s}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot \Phi(s)\Delta s\\ =\ &y(1)\cdot \prod_{u=0}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}+2\alpha\sum_{v=0}^{\log_{3}t-1}\left(\prod_{u=v}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot 3^{v}\cdot \left(y^{\beta}(3^{v})+1\right)^{1/\beta}\\ &+2\sum_{v=0}^{\log_{3}t-1}\left(\prod_{u=v}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot 3^{v}\cdot \Phi(3^{v}),\quad t\in \mathbb{I}. \end{align*}$

对 $ t\in \mathbb{I} $, 方程 (3.3) 也存在唯一解

$\begin{align*} x(t)&=e_{\ominus p}(t,1)y(1)+\int_{1}^{t}e_{\ominus p}(t,s)\cdot \frac{f(s,x(s))}{r(s)-\mu (s)p_{2}(s)}\Delta s\\ &=y(1)\cdot \prod_{u=0}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}+2\alpha\sum_{v=0}^{\log_{3}t-1}\left(\prod_{u=v}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot 3^{v}\cdot \left(x^{\beta}(3^{v})+1\right)^{1/\beta}, \end{align*}$

满足 $ x(1)=y(1) $, 于是有

$\begin{align*} &\left\vert y(t)-x(t)\right\vert\\ \leq\ &2\alpha\sum_{v=0}^{\log_{3}t-1}\left(\prod_{u=v}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot 3^{v}\cdot \left\vert \left(y^{\beta}(3^{v})+1\right)^{1/\beta}-\left(x^{\beta}(3^{v})+1\right)^{1/\beta}\right\vert\\ &+2\sum_{v=0}^{\log_{3}t-1}\left(\prod_{u=v}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot 3^{v}\cdot \left\vert \Phi(3^{v})\right\vert\\ \leq \ &2\alpha\sum_{v=0}^{\log_{3}t-1}\left(\prod_{u=v}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot 3^{v}\cdot \left\vert \left(y^{\beta}(3^{v})+1\right)^{1/\beta}-\left(x^{\beta}(3^{v})+1\right)^{1/\beta}\right\vert +2\cdot \frac{3}{7}\cdot 40\varepsilon \varphi(t)\\ =\ &2\alpha\sum_{v=0}^{\log_{3}t-1}\left(\prod_{u=v}^{\log_{3}t-1}\frac{3}{1+6\cdot 27^{u}}\right)\cdot 3^{v}\cdot \left\vert \left(y^{\beta}(3^{v})+1\right)^{1/\beta}-\left(x^{\beta}(3^{v})+1\right)^{1/\beta}\right\vert +\frac{240}{7}\varepsilon \varphi(t). \end{align*}$

由定理 2.1, 可知当 $ 0<\alpha<7/240 $ 时, 方程 (3.3) 具有 $ \mathcal{N} $ 型 Hyers-Ulam-Rassias 稳定性, 且 HUR 稳定性常数为

$K=\frac{3}{7}\cdot 80\cdot e_{\frac{3}{7}\alpha}(81,1)=\frac{240}{7}\prod_{u=1}^{4}\left(1+\frac{2}{7}\alpha \cdot 3^{u}\right),$

且易知

$K<\frac{240}{7}\prod_{u=1}^{4}\left(1+\frac{2}{7}\cdot \frac{7}{240}\cdot 3^{u}\right)\approx 77.5169.$

例3.3 令 $ \mathbb{T}=\mathbb{R} $, $ 1<a<b $, $ \mathbb{I}=[a,b]_{\mathbb{R}} $. 对 $ t\in \mathbb{I}^{\mathcal{K}} $, 考虑方程

$ (tx(t))'-x(t)+(t-1)x(t)=\lambda t\cdot \sqrt{x^{2}(t)+1},$

$(tx(t))'+(t-2)x(t)=\lambda t\cdot \sqrt{x^{2}(t)+1},$

满足初始条件 $ x(a)=x_{0}\in \mathbb{R} $, 以及对应的不等式

$ \left\vert (ty(t))'+(t-2)y(t)-\lambda t\cdot \sqrt{y^{2}(t)+1}\right\vert \leq \varepsilon \varphi(t).$

这里, $ r(t)=t $, $ p_{1}(t)=-1 $, $ p_{2}(t)=t-1 $, $ f(t,x)=\lambda t\cdot \sqrt{x^{2}+1} $, $ \lambda>0 $, $ \varepsilon>0 $, $ \varphi \in \mathcal{N} $, $ \mu (t)=0 $.

由于

$r^{\sigma}(t)+\mu (t)p_{1}(t)=r(t)-\mu (t)p_{2}(t)=t>0,\quad t\in \mathbb{I}^{\mathcal{K}},$

以及

$\begin{align*} \left\vert f(t,u(t))-f(t,v(t))\right\vert &=\lambda t\cdot \left\vert \sqrt{u^{2}(t)+1}-\sqrt{v^{2}(t)+1}\right\vert \leq \lambda t\cdot \left\vert u(t)-v(t)\right\vert,\quad t\in \mathbb{I}, \end{align*}$

令 $ L=\lambda $, 同样地, 方程 (3.5) 的系数函数满足条件 (A1) 和 (A2). 根据 (2.5) 式, 有

$(\ominus p)(t)=-\frac{r'(t)+p_{1}(t)+p_{2}(t)}{r(t)}=\frac{1-t}{t}.$

于是, 对一切 $ t\in \mathbb{I} $, $ s\in \mathbb{I}^{\mathcal{K}} $ 且 $ t>s $, 有

$\begin{align*} e_{\ominus p}(t,s)={\rm e}^{\int_{s}^{t}\frac{1-r}{r}dr}=\frac{t}{s}{\rm e}^{s-t}. \end{align*}$

当 $ t=s\in \mathbb{I}^{\mathcal{K}} $ 时, 有 $ e_{\ominus p}(t,s)=e_{\ominus p}(t,t)=1 $, 与以上结论保持一致. 所以, 可推出

$E=\sup_{t,s\in \mathbb{I},t>s}\frac{t}{s}{\rm e}^{s-t}=1>0,$

于是当

$0<\lambda<\frac{1}{b-a}$

时, 可满足 $ EL(b-a)=\lambda (b-a)<1 $. 这时, 方程 (3.5) 存在唯一解 $ x $ 满足 $ x(a)=x_{0} $. 类似地, 由不等式 (12), 存在函数 $ \Phi \in {\rm C_{rd}}(\mathbb{I}, \mathbb{R}) $, 使得 $ \left\vert \Phi(t)\right\vert \leq \varepsilon \varphi(t) $ 和

$(ty(t))'+(t-2)y(t)=\lambda t\cdot \sqrt{y^{2}(t)+1}+\Phi(t),\quad t\in \mathbb{I}^{\mathcal{K}}.$

于是, 有

$\begin{align*} y(t)&=\frac{t}{a}{\rm e}^{a-t}\cdot y(a)+\int_{a}^{t}\frac{t}{s}{\rm e}^{s-t}\left[\lambda \cdot \sqrt{y^{2}(t)+1}+\Phi(s)\right]{\rm d}s\\ &=t{\rm e}^{-t}\left[\frac{{\rm e}^{a}}{a}\cdot y(a)+\lambda \int_{a}^{t}\frac{{\rm e}^{s}}{s}\cdot \sqrt{y^{2}(s)+1}{\rm d}s+\int_{a}^{t}\frac{{\rm e}^{s}}{s}\cdot \Phi(s){\rm d}s\right], \end{align*}$

且方程 (3.5) 存在唯一解 $ x $ 满足 $ x(a)=y(a) $. 由引理 2.2, 可知

$x(t)=t{\rm e}^{-t}\left[\frac{{\rm e}^{a}}{a}\cdot y(a)+\lambda \int_{a}^{t}\frac{{\rm e}^{s}}{s}\cdot \sqrt{x^{2}(s)+1}{\rm d}s\right],\quad t\in \mathbb{I}.$

因此, 对 $ t\in \mathbb{I} $, 均有

$\begin{align*} \left\vert y(t)-x(t)\right\vert &\leq \lambda t{\rm e}^{-t}\int_{a}^{t}\frac{{\rm e}^{s}}{s}\cdot \left\vert \sqrt{y^{2}(s)+1}-\sqrt{x^{2}(s)+1}\right\vert {\rm d}s+t{\rm e}^{-t}\int_{a}^{t}\frac{{\rm e}^{s}}{s}\cdot \left\vert \Phi(s)\right\vert {\rm d}s\\ &\leq \lambda t{\rm e}^{-t}\int_{a}^{t}\frac{{\rm e}^{s}}{s}\cdot \left\vert y(s)-x(s)\right\vert {\rm d}s +\varepsilon t{\rm e}^{-t}\int_{a}^{t}\frac{{\rm e}^{s}}{s}\cdot \varphi(s){\rm d}s\\ &\leq \lambda \int_{a}^{t}\left\vert y(s)-x(s)\right\vert {\rm d}s +\varepsilon (b-a)\varphi(t). \end{align*}$

由定理 2.1, 可知方程 (3.5) 具有 $ \mathcal{N} $ 型 Hyers-Ulam-Rassias 稳定性, 且 HUR 稳定性常数为

$K=(b-a)e_{\lambda}(b,a)=(b-a){\rm e}^{\lambda (b-a)},$

且易知 $ K<(b-a)$e.

参考文献

Agarwal R P, Bohner M.

Basic calculus on time scales and some of its applications

Results Math, 1999, 35: 3-22

[本文引用: 1]

Agarwal R P, Bohner M, O'Regan D, Peterson A.

Dynamic equations on time scales: a survey

J Comput Appl Math, 2022, 141: 1-26

[本文引用: 1]

Alghamdi M A, Aljehani A, Bohner M, Hamza A E.

Hyers-Ulam and Hyers-Ulam-Rassias stability of first-order linear dynamic equations

Publ Inst Math, 2021, 109: 83-93

[本文引用: 1]

Bohner M, Peterson A. Dynamic Equations on Time Scales:An Introduction with Applications. Boston: Birkhäuser, 2001

[本文引用: 1]

Bohner M, Peterson A. Advances in Dynamic Equations on Time Scales. Boston: Birkhäuser, 2003: 824-826

[本文引用: 1]

Bohner M, Tikare S.

Ulam stability for first-order nonlinear dynamic equations

Sarajevo J Math, 2022, 18: 83-96

[本文引用: 6]

Bohner M, Tikare S, dos Santos I L D.

First-order nonlinear dynamic initial value problems

Int J Dyn Syst Differ Equ, 2021, 11: 241-254

[本文引用: 1]

Hilger S. Ein Maßkettenkalkül mit Anwendung auf Zentrumsmannigfaltigkeiten[D]. Würzburg: Universität Würzburg, 1988

[本文引用: 1]

Hilger S.

Analysis on measure chain - a unified approach to continuous and discrete calculus

Results Math, 1990, 18: 18-56

[本文引用: 1]

Hyers D H.

On the stability of the linear functional equation

Proc Nat Acad Sci (USA), 1941, 27: 222-224

[本文引用: 2]

Qiu Y C, Wang Q R.

Existence of nonoscillatory solutions to higher-order nonlinear neutral dynamic equations on time scales

Bull Malays Math Sci Soc, 2018, 41: 1935-1952

DOI:10.1007/s40840-016-0435-y      [本文引用: 1]

Rassias T M.

On the stability of the linear mapping in Banach spaces

Proc Amer Math Soc, 1978, 72: 297-300

[本文引用: 1]

Rus I A.

Gronwall lemmas: ten open problems

Sci Math Jpn, 2009, 70: 221-228

[本文引用: 1]

Shen Y.

The Ulam stability of first order linear dynamic equations on time scales

Results Math, 2017, 72: 1881-1895

DOI:10.1007/s00025-017-0725-1      URL     [本文引用: 1]

Ulam S M.

A Collection of Mathematical Problems

%Interscience Tracts in Pure and Applied Mathematics. New York: Interscience Publishers, 1960

[本文引用: 1]

Wong F H, Yeh C C, Hong C H.

Gronwall inequalities on time scales

Math Inequal Appl, 2006, 9: 75-86

[本文引用: 1]

/