数学杂志  2022, Vol. 42 Issue (6): 549-561   PDF    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
李宝麟
宋国鑫
一类固定时刻线性脉冲微分方程解的存在唯一性
李宝麟, 宋国鑫    
西北师范大学数学与统计学院, 甘肃 兰州 730070
摘要:本文研究了一类固定时刻线性脉冲微分方程解的存在唯一性问题,利用广义常微分方程理论,证明了此类脉冲微分方程初值问题解的整体存在唯一性定理,推广了具有逐段连续系数的线性脉冲微分方程解的整体存在唯一性结果.
关键词线性脉冲微分方程    Kurzweil积分    广义常微分方程    存在唯一性    
EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR INITIAL VALUE PROBLEMS OF A CLASS OF LINEAR IMPULSIVE DIFFERENTIAL EQUATIONS AT FIXED TIME
LI Bao-lin, SONG Guo-xin    
College of Mathematics and Statistics, Northwest Normal University, Gansu 730070, China
Abstract: In this paper, we study the existence and uniqueness of solutions for a class of linear impulsive differential equations at fixed times. By using the theory of generalized ordinary differential equations, we prove the global existence and uniqueness theorem of solutions for the initial value problem of such impulsive differential equations, and generalize the global existence and uniqueness result of solutions for linear impulsive differential equations with piecewise continuous coefficients.
Keywords: Linear impulsive differential equations     Kurzweil integral     Generalized ordinary differential equation     Existence and uniqueness    
1 引言

脉冲现象是一种瞬时突变现象, 在现代科技各领域的实际问题中是普遍存在的, 其数学模型往往可归结为脉冲微分系统, 其最突出的特点是能够充分考虑到瞬间突发现象对系统状态的影响. 对于脉冲时刻固定的脉冲微分系统, 大量的研究工作已取得一系列重要成果[1][2]. 文献[1], [2]主要应用Riemann积分与经典常微分方程的理论研究脉冲微分方程. 对于齐次线性脉冲微分方程

$ \left\{ \begin{array}{ll} x'=P(t)x, &t\neq t_{i}, \; i\in\mathbb{Z}, \\ \Delta x|_{t=t_{i}}=C_{i}x, &i\in\mathbb{Z}, \end{array} \right. $

文献[2]考察了其解的存在唯一性与基解矩阵的存在性, 并得出了非齐次线性脉冲微分方程

$ \left\{ \begin{array}{ll} x'=P(t)x+Q(t), &t\neq t_{i}, \; i\in\mathbb{Z}, \\ \Delta x|_{t=t_{i}}=C_{i}x+g_{i}, &i\in\mathbb{Z} \end{array} \right. $

的常数变易公式, 其中

$ \Delta x|_{t=t_{i}}=\lim\limits_{t\to t^{+}_{i}}x(t)-x(t_{i}), \; i\in\mathbb{Z}, $

函数$ P:\mathbb{R}\rightarrow L({\mathbb{C}}^{n}) $$ Q:\mathbb{R}\rightarrow {\mathbb{C}}^{n} $均在$ \mathbb{R}\backslash\{t_{i}\}_{i\in\mathbb{Z}} $上连续, $ L({\mathbb{C}}^{n}) $为复$ n\times n $矩阵全体, $ {\mathbb{C}}^{n} $为复$ n $维欧氏空间, $ \mathbb{Z} $为整数集, $ \{C_{i}\}_{i\in\mathbb{Z}}\subset L({\mathbb{C}}^{n}) $, $ \{g_{i}\}_{i\in\mathbb{Z}}\subset{\mathbb{C}}^{n} $, 数列$ \{t_{i}\}_{i\in\mathbb{Z}} $严格递增且$ \lim\limits_{n\to\infty}t_{n}=+\infty $.

Kurzweil于1957年提出的广义常微分理论为脉冲微分方程的研究提供了新方法, 文献[3]在函数$ f:\overline{B}\times[a, b]\rightarrow {\mathbb{R}}^{n} $关于第二变元具有某些不连续性质的条件下局部地讨论了非线性脉冲微分方程

$ \left\{ \begin{array}{ll} x'=f(x, t)\; a.e.\; 于 \; [a, b]\backslash\{t_{1}, t_{2}, \cdot\cdot\cdot, t_{k}\}, \\ \Delta x|_{t=t_{i}}=I_{i}(x(t_{i})), \; i=1, 2, \cdot\cdot\cdot, k \end{array} \right. $

与一类广义常微分方程之间的等价关系, 其中算子$ I_{i}:\overline{B}\rightarrow {\mathbb{R}}^{n} $连续, $ i=1, 2, \cdot\cdot\cdot, k $,

$ \Delta x|_{t=t_{i}}=\lim\limits_{t\to t^{+}_{i}}x(t)-x(t_{i}), \; i=1, 2, \cdot\cdot\cdot, k, $
$ a\leq t_{1}<t_{2}<\cdot\cdot\cdot<t_{k}\leq b, $

$ \overline{B}=\{x\in{\mathbb{R}}^{n};\|x\|\leq C\} $, $ C>0 $, $ {\mathbb{R}}^{n} $为实$ n $维欧氏空间, $ \|\cdot\| $表示其欧氏范数. 利用上述等价关系及广义常微分方程解关于初值与参数的可微性定理, 文献[4]研究了脉冲微分方程解关于初值与参数的可微性. 文献[5]将文献[3]中广义常微分方程的适定性结果推广至一般Banach空间中的广义常微分方程, 并将所得结论应用于无限滞后型脉冲泛函微分方程, 讨论其解的局部存在唯一性及解相对于初值与参数的连续依赖性.

作为整数阶脉冲微分方程的推广, 分数阶脉冲微分方程已受到广泛关注, 这一领域的研究现已取得了诸多新成果, 文献[6], [7]利用非线性泛函分析理论分别研究了一类在无界区间上具有可列多个脉冲时刻的分数阶脉冲微分方程解的存在性与一类分数阶脉冲微分方程边值问题正解的存在唯一性. 文献[8]研究了一类分数阶脉冲微分方程解的振动性, 并举出相应的算例验证所得结果.

本文将弱化文献[2]中函数$ P $$ Q $连续的条件, 推广文献[3]中的结论, 在整个无界区间$ [t_{0}, +\infty) $上建立线性脉冲微分方程初值问题(IVP)

$ \begin{align} \left\{ \begin{array}{ll} x'=P(t)x+Q(t)\; a.e.\; 于 \; [t_{0}, +\infty)\backslash\{t_{i}\}_{i=1}^{\infty}, \\ \Delta x|_{t=t_{i}}=\lim\limits_{t\to t^{+}_{i}}x(t)-x(t_{i})=C_{i}x(t_{i}), \; i\in\mathbb{N_{+}}, \\ x(t_{0})=x_{0}. \end{array} \right. \end{align} $ (1.1)

与一类广义线性常微分方程初值问题之间的等价关系, 并研究IVP(1.1) 解的整体存在唯一性, 其中函数$ P:[t_{0}, +\infty)\rightarrow L({\mathbb{R}}^n) $$ Q:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $局部Lebesgue可积, $ L({\mathbb{R}}^n) $表示实$ n $阶方阵全体, $ \mathbb{N_{+}} $为正整数集, $ x_{0}\in {\mathbb{R}}^n $, $ \{C_{i}\}_{i=1}^{\infty}\subset L({\mathbb{R}}^n) $, 数列$ \{t_{i}\}_{i=1}^\infty\subset(t_{0}, +\infty) $严格递增且$ \lim\limits_{n\to\infty}t_{n}=+\infty $.

本文共分为三部分, 第二部分介绍文中所用到的基本概念及引理; 第三部分利用Kurzweil积分与广义常微分方程理论证明IVP(1.1) 与一类广义线性常微分方程初值问题之间存在等价关系, 并建立IVP(1.1) 解的整体存在唯一性定理.

2 预备知识

本节将简要介绍Kurzweil积分, 广义常微分方程与正则函数的相关定义和结论.

在本文后续讨论中仍设$ \mathbb{R}=(-\infty, +\infty) $, $ ({\mathbb{R}}^{n}, \|\cdot\|) $为实$ n $维欧氏空间, 且$ L({\mathbb{R}}^n) $为实$ n $阶方阵全体, 其中$ \|\cdot\| $为欧氏范数, 并对任意$ A=(a_{ij})_{n\times n}\in L({\mathbb{R}}^n) $, 规定$ \|A\|_{1}=\sum\limits_{i, j=1}^{n}|a_{ij}|. $

称有限集$ D=\{(\tau_{i}, [\alpha_{i-1}, \alpha_{i}]), \; i=1, 2, \cdot\cdot\cdot, k\} $为区间$ [a, b] $的一个分划, 如果

$ a={\alpha}_{0}<{\alpha}_{1}<\cdot\cdot\cdot<{\alpha}_{k}=b, $
$ {\alpha}_{i-1}\leq{\tau}_{i}\leq{\alpha}_{i}, \; i=1, 2, \cdot\cdot\cdot, k. $

给定正值函数$ \delta:[a, b]\rightarrow (0, +\infty) $, 称$ [a, b] $的分划$ D=\{(\tau_{i}, [\alpha_{i-1}, \alpha_{i}]), \; i=1, 2, \cdot\cdot\cdot, k\} $$ \delta $ –精细的, 如果对每个$ i=1, 2, \cdot\cdot\cdot, k $都有$ [\alpha_{i-1}, \alpha_{i}]\subset(\tau_{i}-\delta(\tau_{i}), \tau_{i}+\delta(\tau_{i})). $

另外, 由文献[3]引理1.4, 若给定区间$ [a, b] $与正值函数$ \delta:[a, b]\rightarrow (0, +\infty) $, 则$ [a, b] $一定存在$ {\delta} $–精细分划$ D=\{(\tau_{i}, [{\alpha}_{i-1}, {\alpha}_{i}]), \; i=1, 2, \cdot\cdot\cdot, k\} $.

定义2.1[3]  称函数$ U:[a, b]\times[a, b] \rightarrow {\mathbb{R}}^{n} $在区间$ [a, b] $上Kurzweil可积, 如果存在向量$ I\in {\mathbb{R}}^{n} $, 使得对任意的$ \varepsilon>0 $, 存在正值函数$ \delta:[a, b]\rightarrow (0, +\infty) $, 使得对$ [a, b] $的任何$ \delta $–精细分划$ D=\{(\tau_{i}, [\alpha_{i-1}, \alpha_{i}]), \; i=1, 2, \cdot\cdot\cdot, k\} $都有

$ \Bigg\|\sum\limits_{i=1}^{k}(U(\tau_{i}, \alpha_{i})-U(\tau_{i}, \alpha_{i-1}))-I\Bigg\|<\varepsilon. $

此时定义$ \int_{a}^{b}DU(\tau, t)=I $$ U $$ [a, b] $上的Kurzweil积分.

特别地, 若$ U(\tau, t)=f(\tau)\cdot t $, 其中$ \tau, \; t\in[a, b] $, 函数$ f:[a, b]\rightarrow {\mathbb{R}}^{n} $, 则记

$ \int_{a}^{b}DU(\tau, t)=\int_{a}^{b}f(s)\text{d}s; $

$ U(\tau, t)=f(\tau)g(t) $, 其中$ \tau, \; t\in[a, b] $, 函数$ f:[a, b]\rightarrow {\mathbb{R}}^{n} $, $ g:[a, b]\rightarrow \mathbb{R} $, 则记

$ \int_{a}^{b}DU(\tau, t)=\int_{a}^{b}f(s)\text{d}g(s), $

并称$ \int_{a}^{b}f(s)\text{d}g(s) $为Kurzwei–Stieltjes积分. 现举以下两例辅助理解Kurzweil积分.

例2.2[9]  考察Kurzweil积分$ \int_{0}^{1}D[t\cdot D(\tau)] $, 其中

$ D(\tau)=\left\{ \begin{array}{ll} 1, \; \tau\in[0, 1]\; \mbox{为有理数}, \\ 0, \; \tau\in[0, 1]\; \mbox{为无理数}, \end{array} \right. $

$ t\cdot D(\tau)=U(\tau, t), \; \tau, \; t\in[0, 1] $, 记$ \{r_{n}\}_{n=1}^{\infty} $$ [0, 1] $中全体有理数所成之集, 对于给定的$ \varepsilon>0 $, 取

$ \delta(\tau)=\left\{ \begin{array}{ll} 1, \;\;\;\;\;\;\;\;\;\;\;\;\; \tau\in[0, 1] 为无理数 , \\ \varepsilon\cdot2^{-(n+1)}, \; \tau=t_{n}, \; n=1, 2, \cdot\cdot\cdot. \end{array} \right. $

又设$ D=\{(\tau_{i}, [{\alpha}_{i-1}, {\alpha}_{i}]), \; i=1, 2, \cdot\cdot\cdot, k\} $$ [0, 1] $的任意$ \delta $–精细分划, 则

$ \begin{align*} &\bigg|\sum\limits_{i=1}^{k}[U(\tau_{i}, \alpha_{i})-U(\tau_{i}, \alpha_{i-1})]-0\bigg|=\bigg|\sum\limits_{{\tau}_{i}\in\{r_{n}\}_{n=1}^{\infty}}D({\tau}_{i})({\alpha}_{i}-{\alpha}_{i-1})\bigg|\\ \leq&\sum\limits_{{\tau}_{i}\in\{r_{n}\}_{n=1}^{\infty}}|D({\tau}_{i})|({\alpha}_{i}-{\alpha}_{i-1})\leq\sum\limits_{{\tau}_{i}\in\{r_{n}\}_{n=1}^{\infty}}2\delta({\tau}_{i})<\sum\limits_{n=1}^{\infty}\frac{\varepsilon}{2^{n}}=\varepsilon, \end{align*} $

因此由定义2.1,

$ \int_{0}^{1}D[t\cdot D(\tau)]=\int_{0}^{1}DU(\tau, t)=0. $

事实上, Kurzweil积分包含Riemann积分与Lebesgue积分.

例2.3[10]  函数

$ f(t)=\left\{ \begin{array}{ll} 0, \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;t=0, \\ 2t\sin\frac{2}{t^{2}}-\frac{2}{t}\cos\frac{2}{t^{2}}, \; t\in(0.1] \end{array}\right. $

$ [0, 1] $上非Lebesgue可积, 但在$ [0, 1] $上Kurzweil可积, 即Kurzweil积分

$ \int_{0}^{1}f(s)\text{d}s=\int_{0}^{1}D[f(\tau)\cdot t] $

存在. 同样, Kurzwei–Stieltjes积分包含Lebesgue–Stieltjes积分.

另外, 由Kurzweil积分的定义不难证明, 若在$ [a, b]\times[a, b] $$ U(\tau, t) $仅是关于变量$ \tau $的函数, 即对任意$ \tau, \; t_{1}, \; t_{2}\in[a, b] $都有$ U(\tau, t_{1})=U(\tau, t_{2}) $, 则$ \int_{a}^{b}DU(\tau, t)=0 $.

定理2.4[3]  设对任意$ c\in(a, b] $, 函数$ U:[a, b]\times[a, b] \rightarrow {\mathbb{R}}^{n} $在区间$ [c, b] $上Kurzweil可积, 且极限

$ \lim\limits_{c\to a^{+}}\bigg[\int_{c}^{b}DU(\tau, t)+U(a, c)-U(a, a)\bigg] $

存在, 则$ U $$ [a, b] $上Kurzweil可积, 且

$ \int_{a}^{b}DU(\tau, t)=\lim\limits_{c\to a^{+}}\bigg[\int_{c}^{b}DU(\tau, t)+U(a, c)-U(a, a)\bigg]. $

定义2.5[3]  设$ O\subset {\mathbb{R}}^n, $区间$ J\subset \mathbb{R} $且函数$ F: O\times J\rightarrow {\mathbb{R}}^n $, 称定义在$ J $的子区间$ P $上的函数$ x: P\rightarrow O $为广义常微分方程

$ \begin{align} \frac{\mathrm{d}x}{\mathrm{d}\tau}=DF(x, t) \end{align} $ (2.1)

的解, 如果对任意$ s_{1}, \; s_{2} \in P $都有$ x(s_{2})-x(s_{1})=\int_{s_{1}}^{s_{2}}DF(x(\tau), t). $特别地, 称

$ \begin{align} \frac{\text{d}x}{\text{d}\tau}=D[A(t)x+g(t)] \end{align} $ (2.2)

为广义线性常微分方程, 其中函数$ A:J\rightarrow L({\mathbb{R}}^n) $, $ g:J\rightarrow {\mathbb{R}}^n $均在$ J $上局部有界变差(即$ A:J\rightarrow L({\mathbb{R}}^n) $$ g:J\rightarrow {\mathbb{R}}^n $均在$ J $的任何闭子区间上有界变差). 事实上, 这是在广义常微分方程(2.1) 中取$ F(x, t)=A(t)x+g(t) $的特殊情形. 在后续讨论中以符号$ \int_{s_{1}}^{s_{2}}\text{d}[A(s)]x(s) $代替$ \int_{s_{1}}^{s_{2}}D[A(t)x(\tau)] $.

引理2.6[3]  若函数$ x:[a, b]\rightarrow {\mathbb{R}}^n $是广义线性常微分方程(2.2) 在$ [a, b] $上的解, 则$ x $$ [a, b] $上的有界变差函数.

定理2.7[3]  设$ J $是一个区间(有界或无界), 函数$ A:J\rightarrow L({\mathbb{R}}^{n}) $$ g:J\rightarrow {\mathbb{R}}^{n} $均在$ J $上局部有界变差, 又设对$ J $的任意内点$ t $都有

$ {\rm det}(I-[A(t)-A(t^{-})])\cdot {\rm det}(I+[A(t^{+})-A(t)])\neq0, $

且当$ \alpha= {\rm inf}J\in J $时,

$ {\rm det}(I+[A(\alpha^{+})-A(\alpha)])\neq0; $

$ \beta= {\rm sup}J\in J $时,

$ {\rm det}(I-[A(\beta)-A(\beta^{-})])\neq0, $

其中$ I $$ n\times n $单位矩阵, 则对任意$ (x_{0}, t_{0})\in{\mathbb{R}}^{n}\times J $, 广义线性常微分方程初值问题

$ \begin{align} \left\{ \begin{aligned} &\frac{\text{d}x}{\text{d}\tau}=D[A(t)x+g(t)], \\ &x(t_{0})=x_{0} \end{aligned} \right. \end{align} $ (2.3)

$ J $上存在唯一解$ x:J\rightarrow {\mathbb{R}}^{n} $, 并且$ x $$ J $上局部有界变差.

引理2.8[3]  设$ [a, b]\subset J $, 函数$ x:[a, b]\rightarrow O $是广义线性常微分方程(2.2) 在$ [a, b] $上的解, 则对任意$ s\in[a, b) $,

$ x(s^{+})-x(s)=[A(s^{+})-A(s)]x(s)+g(s^{+})-g(s), $

对任意$ s\in(a, b] $,

$ x(s)-x(s^{-})=[A(s)-A(s^{-})]x(s)+g(s)-g(s^{-}). $

$ f:[a, b]\rightarrow {\mathbb{R}}^{n} $为正则函数, 如果对于任意$ s\in(a, b] $及任意$ s\in[a, b) $, 极限$ f(s^{-})=\lim\limits_{t\to s^{-}}f(t) $$ f(s^{+})=\lim\limits_{t\to s^{+}}f(t) $分别存在. 特别地, 有界变差函数一定是正则函数, 且定义在某个闭区间上的正则函数一定有界.

引理2.9[3]  设闭球$ \overline{B}=\{x\in{\mathbb{R}}^{n};\|x\|\leq C\} $, 其中$ C>0 $, 函数$ f:\overline{B}\times[a, b]\rightarrow {\mathbb{R}}^{n} $满足:

(ⅰ) 对任意固定的$ x\in \overline{B} $, 函数$ f(x, \cdot) $$ [a, b] $上Lebesgue可测;

(ⅱ) 存在Lebesgue可积函数$ m:[a, b]\rightarrow \mathbb{R} $使得$ \|f(x, s)\|\leq m(s) $对任意$ (x, s)\in \overline{B}\times[a, b] $都成立;

(ⅲ) 存在Lebesgue可积函数$ l:[a, b]\rightarrow \mathbb{R} $使得$ \|f(x, s)-f(y, s)\|\leq l(s)\|x-y\| $对任意$ (x, s), \; (y, s)\; \in \overline{B}\times[a, b] $都成立.

$ \widetilde{F}(x, t)=\int_{a}^{t}f(x, s)\text{d}s, \; (x, t)\in \overline{B}\times[a, b], $

$ x:[a, b]\rightarrow \overline{B} $是正则函数, 则Kurzweil积分$ \int_{a}^{b}D\widetilde{F}(x(\tau), t) $与Lebesgue积分$ \int_{a}^{b}f(x(s), s)\text{d}s $均存在且相等.

引理2.10[3]  设$ A:[a, b]\rightarrow L({\mathbb{R}}^{n}) $$ [a, b] $上的有界变差函数, $ x:[a, b]\rightarrow {\mathbb{R}}^n $$ [a, b] $上的正则函数, 则Kurzweil积分$ \int_{a}^{b}\text{d}[A(s)]x(s) $存在.

3 主要结果

$ x_{0}\in {\mathbb{R}}^n $, $ \{C_{i}\}_{i=1}^{\infty}\subset L({\mathbb{R}}^n) $, 数列$ \{t_{i}\}_{i=1}^\infty\subset(t_{0}, +\infty) $满足

$ t_{1}<t_{2}<\cdot\cdot\cdot<t_{n}<\cdot\cdot\cdot $

$ \lim\limits_{n\to\infty}t_{n}=+\infty $. 又设$ P:[t_{0}, +\infty)\rightarrow L({\mathbb{R}}^n) $$ n\times n $矩阵值函数, $ Q:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $$ n $维列向量函数, 且$ P $, $ Q $均在$ [t_{0}, +\infty) $上局部Lebesgue可积. 本节先建立脉冲微分方程初值问题(IVP)

$ \left\{ \begin{array}{ll} x'=P(t)x+Q(t)\; a.e.\; 于 \; [t_{0}, +\infty)\backslash\{t_{i}\}_{i=1}^{\infty}, \\ \Delta x|_{t=t_{i}}=x(t_{i}^{+})-x(t_{i})=C_{i}x(t_{i}), \; i\in\mathbb{N_{+}}, \\ x(t_{0})=x_{0} \end{array} \right. $

与广义线性常微分方程初值问题

$ \begin{align} \left\{ \begin{aligned} &\frac{\text{d}x}{\text{d}\tau}=D[A(t)x+g(t)], \\ &x(t_{0})=x_{0} \end{aligned} \right. \end{align} $ (3.1)

解之间的等价关系, 其中

$ A(t)=\int_{t_{0}}^{t}P(s)\text{d}s+\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}, \; t\in[t_{0}, +\infty), $
$ g(t)=\int_{t_{0}}^{t}Q(s)\text{d}s, \; t\in[t_{0}, +\infty), $

且对每个$ i\in\mathbb{N_{+}} $,

$ H_{t_{i}}(t)=\left\{ \begin{array}{ll} 0, \; t\leq t_{i}, \\ 1, \; t>t_{i}. \end{array} \right. $

进一步建立IVP(1.1) 在$ [t_{0}, +\infty) $上解的整体存在唯一性定理.

对于IVP(1.1) 的解, 有以下定义.

定义3.1  称函数$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $为IVP(1.1) 在$ [t_{0}, +\infty) $上的解, 如果

(ⅰ) $ x(t_{0})=x_{0} $;

(ⅱ) $ x'(t)=P(t)x(t)+Q(t) $$ [t_{0}, +\infty)\backslash\{t_{i}\}_{i=1}^{\infty} $上几乎处处成立;

(ⅲ) 若$ [a, b]\subset[t_{0}, +\infty) $满足$ [a, b)\cap\{t_{i}\}_{i=1}^{\infty}=\varnothing $, 则$ x(t) $$ [a, b] $上绝对连续;

(ⅳ) $ \Delta x|_{t=t_{i}}=x(t_{i}^{+})-x(t_{i})=C_{i}x(t_{i}), \; i\in\mathbb{N_{+}} $.

引理3.2  设

$ A(t)=\int_{t_{0}}^{t}P(s)\text{d}s+\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}, \; t\in[t_{0}, +\infty), $
$ g(t)=\int_{t_{0}}^{t}Q(s)\text{d}s, \; t\in[t_{0}, +\infty). $

$ A(t) $$ g(t) $均在$ [t_{0}, +\infty) $上局部有界变差.

  由$ \lim\limits_{n\to\infty}t_{n}=+\infty $, 对任意固定的$ t\in[t_{0}, +\infty) $, 存在$ l\in\mathbb{N_{+}} $使得$ t\leq t_{l+1} $, 从而

$ \sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}=\sum\limits_{i=1}^{l}H_{t_{i}}(t)C_{i}, $

$ \sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i} $收敛. 又$ P(t), \; Q(t) $均在$ [t_{0}, +\infty) $上局部Lebesgue可积, 故函数$ A(t) $$ g(t) $均在$ [t_{0}, +\infty) $上有定义. 下证函数$ A(t) $$ g(t) $均在$ [t_{0}, +\infty) $上局部有界变差.

任取$ [a, b]\subset[t_{0}, +\infty) $, 设$ S:a=s_{0}<s_{1}<\cdot\cdot\cdot s_{l}=b $$ [a, b] $的任意分割, 则

$ \sum\limits_{j=1}^{l}\|g(s_{j})-g(s_{j-1})\|=\sum\limits_{j=1}^{l}\|\int_{s_{j-1}}^{s_{j}}Q(s)\text{d}s\|\leq\sum\limits_{j=1}^{l}\int_{s_{j-1}}^{s_{j}}\|Q(s)\|\text{d}s=\int_{a}^{b}\|Q(s)\|\text{d}s, $

$ g(t) $$ [a, b] $上的有界变差函数. 类似地,

$ \sum\limits_{j=1}^{l}\bigg\|\int_{s_{j-1}}^{s_{j}}P(s)\text{d}s\bigg\|_{1}\leq\sum\limits_{j=1}^{l}\int_{s_{j-1}}^{s_{j}}\|P(s)\|_{1}\text{d}s=\int_{a}^{b}\|P(s)\|_{1}\text{d}s, $

$ \begin{align*} &\sum\limits_{j=1}^{l}\bigg\|\sum\limits_{i=1}^{\infty}H_{t_{i}}(s_{j})C_{i}-\sum\limits_{i=1}^{\infty}H_{t_{i}}(s_{j-1})C_{i}\bigg\|_{1} \leq\lim\limits_{n\to\infty}\sum\limits_{j=1}^{l}\sum\limits_{i=1}^{n}[H_{t_{i}}(s_{j})-H_{t_{i}}(s_{j-1})]\|C_{i}\|_{1}\\ =&\lim\limits_{n\to\infty}\sum\limits_{i=1}^{n}\|C_{i}\|_{1}\sum\limits_{j=1}^{l}[H_{t_{i}}(s_{j})-H_{t_{i}}(s_{j-1})] =\lim\limits_{n\to\infty}\sum\limits_{i=1}^{n}\|C_{i}\|_{1}[H_{t_{i}}(b)-H_{t_{i}}(a)]\\ =&\sum\limits_{i=1}^{\infty}\|C_{i}\|_{1}[H_{t_{i}}(b)-H_{t_{i}}(a)], \end{align*} $

因此

$ \begin{align*} \sum\limits_{j=1}^{l}\|A(s_{j})-A(s_{j-1})\|_{1}&=\sum\limits_{j=1}^{l}\bigg\|\int_{s_{j-1}}^{s_{j}}P(s)\text{d}s+\sum\limits_{i=1}^{\infty}H_{t_{i}}(s_{j})C_{i}-\sum\limits_{i=1}^{\infty}H_{t_{i}}(s_{j-1})C_{i}\bigg\|_{1}\\ &\leq\sum\limits_{j=1}^{l}\bigg\|\int_{s_{j-1}}^{s_{j}}P(s)\text{d}s\bigg\|_{1}+\sum\limits_{j=1}^{l}\bigg\|\sum\limits_{i=1}^{\infty}H_{t_{i}}(s_{j})C_{i}-\sum\limits_{i=1}^{\infty}H_{t_{i}}(s_{j-1})C_{i}\bigg\|_{1}\\ &\leq\int_{a}^{b}\|P(s)\|_{1}\text{d}s+\sum\limits_{i=1}^{\infty}\|C_{i}\|_{1}[H_{t_{i}}(b)-H_{t_{i}}(a)], \end{align*} $

即函数$ A(t) $$ [a, b] $上有界变差, 再由$ [a, b]\subset[t_{0}, +\infty) $的任意性, 函数$ A(t) $$ g(t) $$ [t_{0}, +\infty) $上局部有界变差.

$ F_{1}(x, t)=\int_{t_{0}}^{t}[P(s)x+Q(s)]\text{d}s, \; (x, t)\in{\mathbb{R}}^n\times[t_{0}, +\infty), $

则有如下引理.

引理3.3  设$ [a, b]\subset[t_{0}, +\infty) $, $ x:[a, b]\rightarrow {\mathbb{R}}^{n} $为正则函数, 则Kurzweil积分$ \int_{a}^{b}DF_{1}(x(\tau), t) $与lebesgue积分$ \int_{a}^{b}[P(s)x(s)+Q(s)]\text{d}s $均存在且相等.

  因为$ x:[a, b]\rightarrow {\mathbb{R}}^{n} $为正则函数, 所以存在常数$ C>0 $, 使对任意$ t\in[a, b] $都有

$ \|x(t)\|\leq C, $

从而$ x $$ [a, b] $$ \overline{B}=\{x\in{\mathbb{R}}^{n};\|x\|\leq C\} $. 因函数$ P(t) $$ Q(t) $均在$ [t_{0}, +\infty) $上局部Lebesgue可积, 故对一切固定的$ x\in\overline{B} $, 函数

$ f(x, t)=P(t)x+Q(t) $

$ [a, b] $上关于变量$ t $ Lebesgue可测, 且函数$ m(t)=C\|P(t)\|_{1}+\|Q(t)\|, \; t\in[a, b] $$ l(t)=\|P(t)\|_{1}, \; t\in[a, b] $均在$ [a, b] $上Lebesgue可积.

另一方面, 对任意$ (x, s), \; (y, s)\; \in\overline{B}\times[a, b] $

$ \|f(x, s)\|=\|P(s)x+Q(s)\|\leq\|P(s)\|_{1}\cdot\|x\|+\|Q(s)\|\leq m(s), $

$ \|f(x, s)-f(y, s)\|=\|P(s)(x-y)\|\leq l(s)\|x-y\|. $

再令$ \widetilde{F}(x, t)= \int_{a}^{t}f(x, s)\text{d}s, \; (x, t)\in\overline{B}\times[a, b], $由引理2.9, Kurzweil积分$ \int_{a}^{b}D\widetilde{F}(x(\tau), t) $与Lebesgue积分$ \int_{a}^{b}f(x(s), s)\text{d}s $均存在且相等, 于是

$ \int_{a}^{b}DF_{1}(x(\tau), t)=\int_{a}^{b}D\widetilde{F}(x(\tau), t)=\int_{a}^{b}f(x(s), s)\text{d}s= \int_{a}^{b}[P(s)x(s)+Q(s)]\text{d}s. $

$ F_{2}(x, t)=\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}x, \; (x, t)\in{\mathbb{R}}^n\times[t_{0}, +\infty), $

$ A(t)x+g(t)=F_{1}(x, t)+F_{2}(x, t), \; (x, t)\in{\mathbb{R}}^n\times[t_{0}, +\infty), $

且有如下引理.

引理3.4  设区间$ [a, b]\subset[t_{0}, +\infty) $满足$ [a, b)\cap\{t_{i}\}_{i=1}^{\infty}=\varnothing $, 则对任意函数$ x:[a, b]\rightarrow {\mathbb{R}}^n $都有

$ \int_{a}^{b}DF_{2}(x(\tau), t)=0. $

  因$ [a, b)\cap\{t_{i}\}_{i=1}^{\infty}=\varnothing $, 故$ [a, b]\subset[t_{0}, t_{1}] $或存在$ k\in\mathbb{N_{+}} $, 使得$ [a, b]\subset(t_{k}, t_{k+1}] $.

$ [a, b]\subset[t_{0}, t_{1}] $时, 对任意$ t\in[a, b] $都有$ t\leq t_{1} $, 此时对任一函数$ x:[a, b]\rightarrow {\mathbb{R}}^n $

$ F_{2}(x(\tau), t)=\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}x(\tau)=0, $

$ \int_{a}^{b}DF_{2}(x(\tau), t)=0; $

$ [a, b]\subset(t_{k}, t_{k+1}] $, $ k\in\mathbb{N_{+}} $时, 对任意$ t\in[a, b] $$ t_{k}<t\leq t_{k+1} $, 此时对任一函数$ x:[a, b]\rightarrow {\mathbb{R}}^n $

$ F_{2}(x(\tau), t)=\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}x(\tau)=\sum\limits_{i=1}^{k}H_{t_{i}}(t)C_{i}x(\tau)+\sum\limits_{i=k+1}^{\infty}H_{t_{i}}(t)C_{i}x(\tau) =\sum\limits_{i=1}^{k}C_{i}x(\tau), $

于是由Kurzweil积分的定义得$ \int_{a}^{b}DF_{2}(x(\tau), t)=0. $

引理3.5  设$ k\in\mathbb{N_{+}} $, $ s\in(t_{k}, t_{k+1}] $, 则对任意函数$ x:[t_{k}, s]\rightarrow {\mathbb{R}}^n $, 都有

$ \int_{t_{k}}^{s}DF_{2}(x(\tau), t)=C_{k}x(t_{k}). $

  对任意$ k\in\mathbb{N_{+}} $, $ s\in(t_{k}, t_{k+1}] $及任意$ \sigma\in(t_{k}, s) $显然有$ [\sigma, s)\cap\{t_{i}\}_{i=1}^{\infty}=\varnothing $, 因此由引理3.4,

$ \int_{\sigma}^{s}DF_{2}(x(\tau), t)=0 $

对任意$ \sigma\in(t_{k}, s) $都成立. 由定理2.4, 当$ k\in\mathbb{N_{+}} $, $ k\geq2 $时,

$ \begin{align*} \int_{t_{k}}^{s}DF_{2}(x(\tau), t)&=\lim\limits_{\sigma\to t^{+}_{k}}\bigg[\int_{\sigma}^{s}DF_{2}(x(\tau), t)+F_{2}(x(t_{k}), \sigma)-F_{2}(x(t_{k}), t_{k})\bigg]\\ &=\lim\limits_{\sigma\to t^{+}_{k}}\sum\limits_{i=1}^{\infty}H_{t_{i}}(\sigma)C_{i}x(t_{k})-\sum\limits_{i=1}^{\infty}H_{t_{i}}(t_{k})C_{i}x(t_{k})\\ &=\sum\limits_{i=1}^{k}C_{i}x(t_{k})-\sum\limits_{i=1}^{k-1}C_{i}x(t_{k})\\ &=C_{k}x(t_{k}), \end{align*} $

类似地,

$ \int_{t_{1}}^{s}DF_{2}(x(\tau), t)=\lim\limits_{\sigma\to t^{+}_{1}}\sum\limits_{i=1}^{\infty}H_{t_{i}}(\sigma)C_{i}x(t_{1})-\sum\limits_{i=1}^{\infty}H_{t_{i}}(t_{1})C_{i}x(t_{1})=C_{1}x(t_{1}). $

定理3.6  函数$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $是线性脉冲微分方程初值问题(1.1) 在$ [t_{0}, +\infty) $上的解当且仅当$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $是广义线性常微分方程初值问题(3.1) 在$ [t_{0}, +\infty) $上的解.

  必要性设$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $是IVP(1.1) 在$ [t_{0}, +\infty) $上的解, 对任意$ s_{1}, \; s_{2}\in[t_{0}, +\infty) $, $ s_{1}<s_{2} $, 分以下两种情形讨论:

$ [s_{1}, s_{2})\cap\{t_{i}\}_{i=1}^{\infty}=\varnothing $时, 由定义3.1, $ x $$ [s_{1}, s_{2}] $上绝对连续, 且

$ x'(s)=P(s)x(s)+Q(s) $

a.e. 于$ [s_{1}, s_{2}] $, 因此由引理3.3及引理3.4,

$ \begin{align*} x(s_{2})-x(s_{1})&=\int_{s_{1}}^{s_{2}}x'(s)\text{d}s\\ &=\int_{s_{1}}^{s_{2}}[P(s)x(s)+Q(s)]\text{d}s\\ &=\int_{s_{1}}^{s_{2}}DF_{1}(x(\tau), t)\\ &=\int_{s_{1}}^{s_{2}}DF_{1}(x(\tau), t)+\int_{s_{1}}^{s_{2}}DF_{2}(x(\tau), t)\\ &=\int_{s_{1}}^{s_{2}}DF(x(\tau), t). \end{align*} $

$ [s_{1}, s_{2})\cap\{t_{i}\}_{i=1}^{\infty}\neq\varnothing $时, 不妨设存在$ k\in\mathbb{N_{+}} $, 使得$ t_{k}=s_{1}<s_{2}\leq t_{k+1} $, 则由引理3.3, 对任意$ \sigma\in(s_{1}, s_{2}) $

$ x(s_{2})-x(\sigma)=\int_{\sigma}^{s_{2}}x'(s)\text{d}s=\int_{\sigma}^{s_{2}}[P(s)x(s)+Q(s)]\text{d}s=\int_{\sigma}^{s_{2}}DF_{1}(x(\tau), t), $

由引理3.5,

$ x(t_{k}^{+})-x(t_{k})=C_{k}x(t_{k})=\int_{t_{k}}^{s_{2}}DF_{2}(x(\tau), t)=\int_{s_{1}}^{s_{2}}DF_{2}(x(\tau), t). $

再由定理2.4及Lebesgue积分的绝对连续性,

$ \begin{align*} \int_{s_{1}}^{s_{2}}DF_{1}(x(\tau), t)&=\int_{t_{k}}^{s_{2}}DF_{1}(x(\tau), t)\\ &=\lim\limits_{\sigma\to t^{+}_{k}}\bigg[\int_{\sigma}^{s_{2}}DF_{1}(x(\tau), t)+ F_{1}(x(t_{k}), \sigma)-F_{1}(x(t_{k}), t_{k})\bigg]\\ &=\lim\limits_{\sigma\to t^{+}_{k}}\int_{\sigma}^{s_{2}}DF_{1}(x(\tau), t) +\lim\limits_{\sigma\to t^{+}_{k}}\int_{t_{k}}^{\sigma}[P(s)x(t_{k})+Q(s)]\text{d}s\\ &=x(s_{2})-x(t^{+}_{k}), \end{align*} $

因此

$ \begin{align*} x(s_{2})-x(s_{1})&=x(s_{2})-x(t_{k})\\ &=x(s_{2})-x(t^{+}_{k})+x(t^{+}_{k})-x(t_{k})\\ &=\int_{s_{1}}^{s_{2}}DF_{1}(x(\tau), t)+\int_{s_{1}}^{s_{2}}DF_{2}(x(\tau), t)\\ &=\int_{s_{1}}^{s_{2}}DF(x(\tau), t). \end{align*} $

综上所述, 对任意$ s_{1}, s_{2}\in[t_{0}, +\infty) $,

$ x(s_{2})-x(s_{1})=\int_{s_{1}}^{s_{2}}DF(x(\tau), t)=\int_{s_{1}}^{s_{2}}\text{d}[A(s)]x(s)+g(s_{2})-g(s_{1}), $

由定义2.5, 函数$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $是广义线性常微分方程初值问题(3.1) 在$ [t_{0}, +\infty) $上的解.

充分性  设$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $是IVP(3.1) 在$ [t_{0}, +\infty) $上的解, 则由引理2.6, 对任意固定的$ k\in\mathbb{N_{+}} $, $ x $$ [t_{k}, t_{k+1}] $上的有界变差函数, 当然是正则函数, 因此由定义2.5, 引理3.3及引理3.5, 对每个固定的$ k\in\mathbb{N_{+}} $及任意$ s\in(t_{k}, t_{k+1}) $

$ \begin{align*} x(s)-x(t_{k})&=\int_{t_{k}}^{s}\text{d}[A(\tau)]x(\tau)+g(s)-g(t_{k})\\ &=\int_{t_{k}}^{s}DF_{1}(x(\tau), t)+\int_{t_{k}}^{s}DF_{2}(x(\tau), t)\\ &=\int_{t_{k}}^{s}[P(t)x(t)+Q(t)]\text{d}t+C_{k}x(t_{k}), \end{align*} $

$ x(s)=\int_{t_{k}}^{s}[P(t)x(t)+Q(t)]\text{d}t+C_{k}x(t_{k})+x(t_{k}), \; s\in(t_{k}, t_{k+1}), $

因此对每个固定的$ k\in\mathbb{N_{+}} $, $ x'(s)={\varphi}_{k}'(s)=P(s)x(s)+Q(s) $$ (t_{k}, t_{k+1}) $上几乎处处成立, 其中

$ \varphi_{k}(t)=\int_{t_{k}}^{t}[P(s)x(s)+Q(s)]\text{d}s, \; t\in[t_{k}, t_{k+1}], \; k\in\mathbb{N_{+}}. $

类似地, 对任意$ s\in[t_{0}, t_{1}] $

$ \begin{align*} x(s)-x(t_{0})&=\int_{t_{0}}^{s}DF(x(\tau), t)\\ &=\int_{t_{0}}^{s}DF_{1}(x(\tau), t)+\int_{t_{0}}^{s}DF_{2}(x(\tau), t)\\ &=\int_{t_{0}}^{s}DF_{1}(x(\tau), t)\\ &=\int_{t_{0}}^{s}[P(t)x(t)+Q(t)]\text{d}t, \end{align*} $

$ x(s)=\int_{t_{0}}^{s}[P(t)x(t)+Q(t)]\text{d}t+x_{0}, \; s\in[t_{0}, t_{1}], $

因此$ x'(s)=P(s)x(s)+Q(s) $ a.e. 于$ [t_{0}, t_{1}) $.

综合以上讨论可知: $ x'(s)=P(s)x(s)+Q(s) $$ [t_{0}, +\infty)\backslash\{t_{i}\}_{i=1}^{\infty} $上几乎处处成立.

由引理2.8及Lebesgue积分的绝对连续性: 对任意$ k\in\mathbb{N_{+}} $, $ k\geq2 $

$ \begin{align*} x(t_{k}^{+})-x(t_{k})&=[A(t_{k}^{+})-A(t_{k})]x(t_{k})+g(t_{k}^{+})-g(t_{k})\\ &=\lim\limits_{t\to t^{+}_{k}}\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}x(t_{k})-\sum\limits_{i=1}^{\infty}H_{t_{i}}(t_{k})C_{i}x(t_{k})\\ &=\sum\limits_{i=1}^{k}C_{i}x(t_{k})-\sum\limits_{i=1}^{k-1}C_{i}x(t_{k})\\ &=C_{k}x(t_{k}), \end{align*} $

类似可得

$ x(t_{1}^{+})-x(t_{1})=\lim\limits_{t\to t^{+}_{1}}\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}x(t_{1})-\sum\limits_{i=1}^{\infty}H_{t_{i}}(t_{1})C_{i}x(t_{1})=C_{1}x(t_{1}). $

任取闭区间$ [a, b]\subset [t_{0}, +\infty) $满足$ [a, b)\cap\{t_{i}\}_{i=1}^{\infty}=\varnothing $, 下证函数$ x $$ [a, b] $上绝对连续. 设$ (a_{p}, b_{p}), \; p=1, 2, \cdot\cdot\cdot, m $$ [a, b] $中的任意有限多个互不相交的开区间, 由Lebesgue积分的绝对连续性, 对任意$ \varepsilon>0 $, 存在$ \delta>0 $, 使当$ \sum\limits_{p=1}^{m}(b_{p}-a_{p})<\delta $时,

$ \sum\limits_{p=1}^{m}\|x(b_{p})-x(a_{p})\|=\sum\limits_{p=1}^{m}\bigg\|\int_{a_{p}}^{b_{p}}[P(s)x(s)+Q(s)]\text{d}s\bigg\|<\varepsilon, $

$ x(t) $$ [a, b] $上绝对连续.

综上所述, $ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $是IVP(1.1) 在$ [t_{0}, +\infty) $上的解.

定理3.7  若对任意$ i\in\mathbb{N_{+}} $, 矩阵$ I+C_{i} $可逆, 则IVP(1.1) 在区间$ [t_{0}, +\infty) $上存在唯一解$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $.

  易证函数

$ A(t)=\int_{t_{0}}^{t}P(s)\text{d}s+\sum\limits_{i=1}^{\infty}H_{t_{i}}(t)C_{i}, \; t\in[t_{0}, +\infty) $

$ (t_{0}, +\infty) $上左连续, 故对任意$ t\in(t_{0}, +\infty) $

$ I-[A(t)-A(t^{-})]=I. $

另一方面, 对任意$ t\in[t_{0}, +\infty)\backslash\{t_{i}\}_{i=1}^{\infty} $, 有

$ I+[A(t^{+})-A(t)]=I, $

且对任意$ t=t_{i} $, $ i\in\mathbb{N_{+}} $, 有

$ I+[A(t^{+})-A(t)]=I+C_{i}. $

由以上讨论可得

$ {\rm det}(I-[A(t)-A(t^{-})])\cdot {\rm det}(I+[A(t^{+})-A(t)])\neq0, \; t\in(t_{0}, +\infty), $

$ {\rm det}(I+[A(t_{0}^{+})-A(t_{0})])\neq0. $

由引理3.2及定理2.7, IVP

$ \begin{align} \left\{ \begin{aligned} &\frac{\text{d}x}{\text{d}\tau}=D[A(t)x+g(t)], \\ &x(t_{0})=x_{0} \end{aligned} \right. \end{align} $ (3.1)

$ [t_{0}, +\infty) $上存在唯一解$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $, 再由定理3.6, 函数$ x:[t_{0}, +\infty)\rightarrow {\mathbb{R}}^n $是IVP(1.1) 在$ [t_{0}, +\infty) $上的唯一解.

参考文献
[1] 傅希林, 闫宝强, 刘衍胜. 脉冲微分系统引论[M]. 北京: 科学出版社, 2005.
[2] 宋新宇, 郭红建, 师向云. 脉冲微分方程理论及其应用[M]. 北京: 科学出版社, 2011.
[3] Schwabik Š. Generalized ordinary differential equations[M]. Singapore: World Scientific, 1992.
[4] Slavík A. Generalized differential equations: Differentiability of solutions with respect to initial conditions and parameters[J]. Journal of Mathematical Analysis & Applications, 2013, 402(1): 261–274.
[5] Slavík A. Well-posedness results for abstract generalized differential equations and measure functional differential equations[J]. Journal of Differential Equations, 2015, 259(2): 666–707.
[6] 秦宝侠, 刘文斌. 无穷区间上分数阶微分方程脉冲解的存在性[J]. 吉林大学学报(理学版), 2016, 54(1): 8–14. DOI:10.13413/j.cnki.jdxblxb.2016.01.02
[7] Zheng Fengxia. New results for a class of boundary value problems involving impulsive fractional differential equations[J]. Filomat, 2020, 34(3): 707–725. DOI:10.2298/FIL2003707Z
[8] Feng Limei, Han Zhenlai. Oscillation behavior of solution of impulsive fractional differential equations[J]. Journal of Applied Analysis & Computation, 2020, 10(1): 223–233.
[9] 丁传松, 李秉彝. 广义黎曼积分[M]. 北京: 科学出版社, 1989.
[10] Collegari R, Federson M, Frasson M. Linear FDEs in the frame of generalized ODEs: variation-of-constants formula[J]. Czechoslovak Mathematical Journal, 2018, 68(143): 889–920.