异质环境下扩散Holling-Tanner模型正稳态解的全局稳定性

Global Stability of the Positive Steady-State Solution of the Diffusion Holling-Tanner Model in a Heterogeneous Environment

摘要: 考虑异质环境对带有扩散项的Holling-Tanner捕食–食饵模型正稳态解全局稳定性的影响。首先利用抛物方程的比较定理等方法证明了该模型正稳态解的存在性，进一步通过迭代过程，得出了该模型正稳态解的下界。最后通过构造一种新的李雅普诺夫函数，研究Holling-Tanner捕食-食饵模型正稳态解的稳定性：该模型在环境异质下的正稳态解是全局稳定的。

Abstract: Consider the influence of heterogeneous environment on the global stability of the positive steady-state solution of the Holling-Tanner predator-prey model with diffusion term. First, the existence of the positive steady-state solution of the model is proved by the comparison theorem of the parabolic equation, and the lower bound of the positive steady-state solution of the model is obtained through the iterative process. Finally, by constructing a new Lyapunov function, the stability of the positive steady-state solution of the Holling-Tanner predator-prey model is studied: the positive steady-state solution of this model is globally stable under heterogeneous environments.

1. 引言

生物数学模型是人们研究种群规模发展变化的重要手段之一 [1]。近年来，捕食–食饵种群模型研究取得了很好的结果。对于Holling-Tanner捕食–食饵模型，很多国内外学者取得了有意义的成果。例如，文献 [2] 用上下解方法讨论了在更简单的参数条件下，该模型的全局稳定性；文献 [3] 通过一种新的方法建立了该模型唯一正平衡解的全局渐进稳定性等等。但是环境异质对物种产生的影响这类研究甚少，在种群竞争中，环境异质可能会改变竞争的结果，理解空间非均质环境如何影响物种间的竞争是一个重要的生物学问题。本文将利用上下解方法，以及一种新构造的Lyapunov函数方法，讨论Holling-Tanner捕食–食饵模型在空间异质情况下正稳态解的全局稳定性。

考虑在异质环境下带有扩散项的Holling-Tanner捕食–食饵模型：

$\{\begin{array}{l}{u}_t=d_1\left(x\right)\Delta u+a\left(x\right)u-u^2-\frac{buv}{m+ru},x\in \Omega ,t>0,\\ v_t=d_2\left(x\right)\Delta v+bv-\frac{b{v}^2}{\gamma u},x\in \Omega ,t>0,\\ \frac{\partial u\left(x,t\right)}{\partial \nu }=\frac{\partial v\left(x,t\right)}{\partial \nu }=0,x\in \partial \Omega ,t>0,\\ u\left(x,0\right)=u_0\left(x\right)>0,v\left(x,0\right)=v_0\left(x\right)\ge (\ne )0,x\in \Omega \end{array}$ (1.1)

其中， $u\left(x,t\right)$ 和 $v\left(x,t\right)$ 分别代表捕食者和食饵的密度函数， $\Omega$ 是 $R^n$ 中具有光滑边界 $\partial \Omega$ 的有界区域，施加在 $\partial \Omega$ 上的无通量边界条件表示它们生活在一个封闭的环境里。 $d_1\left(x\right)$ 和 $d_2\left(x\right)$ 分别是u和v的空间依赖扩散系数函数， $a\left(x\right)$ 是空间异质资源函数，其中 $b,r,m,\gamma$ 都是常数。

2. 预备知识

为了证明异质环境下扩散Holling-Tanner模型正稳态解的存在性和全局稳定性，引入如下引理。

引理1 [4] (比较定理)假设 $f\left(s\right)\in C^1\left(\left[0,\infty \right)\right)$ 是正函数，常数 $d>0,\beta >0$， $T\in \left[0,\infty \right)$， $\omega \in C^{2,1}\left(\Omega \times \left(T,\infty \right)\right)\cap C^{1,0}\left(\bar{\Omega }\times \left[T,\infty \right)\right)$ 是正函数。如果 $\omega$ 满足：

$\{\begin{array}{l}\frac{\partial \omega }{\partial t}-d\Delta \omega \le {\omega }^{1+\beta }f\left(\omega \right)\left(\alpha -\omega \right),x\in \Omega ,t\in \left(T,\infty \right)\\ \frac{\partial \omega }{\partial n}=0,x\in \partial \Omega ,t\in \left[T,\infty \right)\end{array}$

且常数 $\alpha >0$，那么

$\underset{t\to \infty }{\lim }\sup \underset{\bar{\Omega }}{\max }\omega \left(\cdot ,t\right)\le \alpha$ (2.1)

引理2 [5] 假设 $\delta >0$，函数 $\psi ,h\in C\left(\left[\delta ,\infty \right)\right)$ 分别满足 $\psi \left(t\right)\ge 0$，且 ${\displaystyle {\int }_{\delta }^{\infty }h\left(t\right)\text{d}t}<\infty$。若 $\psi \in C^1\left(\left[\delta ,\infty \right)\right)$ 是有界的，并且

${\phi }^{\prime }\left(t\right)\le -\psi \left(t\right)+h\left(t\right)$， ${\psi }^{\prime }\left(t\right)\le K$

(其中 $K>0$ )

那么 $\underset{t\to \infty }{\lim }\psi \left(t\right)=0$。

3. 异质环境下扩散Holling-Tanner模型正稳态解的存在性

系统(1.1)对应的稳态问题为：

$\{\begin{array}{l}-d_1\left(x\right)\Delta u=u\left(a\left(x\right)-u-\frac{bv}{m+ru}\right),x\in \Omega \\ -d_2\left(x\right)\Delta v=bv\left(1-\frac{v}{\gamma u}\right),x\in \Omega \\ \frac{\partial u\left(x,t\right)}{\partial \nu }=\frac{\partial v\left(x,t\right)}{\partial \nu }=0,x\in \partial \Omega \end{array}$

本节给出扩散Holling-Tanner模型正稳态解的存在性结果。

定理1 系统(1.1)存在唯一解 $\left(u\left(x,t\right),v\left(x,t\right)\right)$，并满足：对任意小的 ${\epsilon }_1,{\epsilon }_2$ (其中 ${\epsilon }_1>0,{\epsilon }_2\in \left(0,{\epsilon }_1\right]$ )，存在一个常数 $T>0$，当 $t>T$ 时，有

1) ${\epsilon }_2\le u\left(x,t\right),v\left(x,t\right)\le \bar{c}+{\epsilon }_1$ , $\forall x\in \bar{\Omega }.$ (3.1)

2) 存在一个常数 $G>0$，使得

$\underset{t>T_1}{\max }{\Vert u\left(\cdot ,t\right)\Vert }_{C^{2+\alpha }\left(\bar{\Omega }\right)},\underset{t>T_1}{\max }{\Vert v\left(\cdot ,t\right)\Vert }_{C^{2+\alpha }\left(\bar{\Omega }\right)}\le G.$ (3.2)

其中

$\bar{c}:=\underset{\bar{\Omega }}{\max }a\left(x\right)$ , $\underset{\_}{c}:=\underset{\bar{\Omega }}{\min }a\left(x\right).$ (3.3)

证明 由系统(1.1)的第一个方程， $u\left(x,t\right)$ 满足

$u_t-d_1\left(x\right)\Delta u\le u\left(\underset{x\in \bar{\Omega }}{\max }a\left(x\right)-u\right)$

利用引理1，得

$\underset{t\to \infty }{\lim \sup }\underset{\bar{\Omega }}{\max }u\left(x,t\right)\le \underset{x\in \bar{\Omega }}{\max }a\left(x\right)=\bar{c}$

那么对任意的 $\epsilon >0$，存在 $T>0$，使得 $u\left(x,t\right)<\bar{c}+\epsilon$。故模型(1.1)的第二个方程 $v_t=d_2\left(x\right)\Delta v+bv-\frac{b{v}^2}{\gamma u}$ 中， $v\left(x,t\right)$ 满足

$v_t-d_2\left(x\right)\Delta v\le bv\left(1-\frac{v}{\gamma \left(\bar{c}+\epsilon \right)}\right)$

结合系统(1.1)的边界条件 $\frac{\partial u}{\partial \nu }=0$，同理，利用引理1得出

$\underset{t\to \infty }{\lim \sup }\underset{\bar{\Omega }}{\max }v\left(x,t\right)\le \bar{c}+\epsilon$

故存在 $T>0$，使得

$u\left(x,t\right),v\left(x,t\right)\le \bar{c}+{\epsilon }_1$ , $x\in \bar{\Omega }$ , $t>T.$

又因 $u\left(x,T\right),v\left(x,T\right)>0$，我们可以选择一个很小的常数 ${\epsilon }_2$，其中 ${\epsilon }_2\in \left(0,{\epsilon }_1\right]$，使得 $u\left(x,T\right),v\left(x,T\right)>{\epsilon }_2$

故

${\epsilon }_2\le u\left(x,t\right),v\left(x,t\right)\le \bar{c}+{\epsilon }_1.$

因此，问题(1.1)在 $\left[{\epsilon }_2,\bar{c}+{\epsilon }_1\right]\times \left[{\epsilon }_2,\bar{c}+{\epsilon }_1\right]$ 上存在一个正稳态解 $\left(u\left(x,t\right),v\left(x,t\right)\right)$。

4. 异质环境下扩散Holling-Tanner模型正稳态解的下界

对任意的 $T>0$，定义 ${\bar{u}}_1={\bar{v}}_1=\bar{a}+{\epsilon }_1$， ${\underset{\_}{u}}_1={\underset{\_}{v}}_1={\epsilon }_2$。显然有

$\left|a\left(x\right)u_1-u_1^2-\frac{b{u}_1{v}_1}{m+r{u}_1}-a\left(x\right)u_2+u_2^2+\frac{b{u}_2{v}_2}{m+r{u}_2}\right|\le T\left(\left|u_1-u_2\right|+\left|v_1-v_2\right|\right),$ (4.1)

$\left|b{v}_1-\frac{b{v}_1^2}{r{u}_1}-b{v}_2+\frac{b{v}_2^2}{r{u}_2}\right|\le T\left(\left|u_1-u_2\right|+\left|v_1-v_2\right|\right).$

定义如下的迭代序列

$\{\begin{array}{l}{\bar{u}}_{i+1}={\bar{u}}_i+\frac{1}{T}\left(\bar{a}{\bar{u}}_i-{\bar{u}}_i^2-\frac{b{\bar{u}}_i{\underset{\_}{v}}_i}{m+r{\bar{u}}_i}\right),\\ {\underset{\_}{u}}_{i+1}={\underset{\_}{u}}_i+\frac{1}{T}\left(\underset{\_}{a}{\underset{\_}{u}}_i-{\underset{\_}{u}}_i^2-\frac{b{\underset{\_}{u}}_i{\bar{v}}_i}{m+r{\underset{\_}{u}}_i}\right),\\ {\bar{v}}_{i+1}={\bar{v}}_i+\frac{1}{T}\left(b{\bar{v}}_i-\frac{b{\bar{v}}_i^2}{r{\bar{u}}_i}\right),\\ {\underset{\_}{v}}_{i+1}={\underset{\_}{v}}_i+\frac{1}{T}\left(b{\underset{\_}{v}}_i-\frac{b{\underset{\_}{v}}_i^2}{r{\underset{\_}{u}}_i}\right).\end{array}$ (4.2)

定理2 迭代序列(4.2)满足以下的单调性和收敛性：

1) 单调性：常数序列 ${\left\{{\bar{u}}_i\right\}}_{i=1}^{\infty }$， ${\left\{{\underset{\_}{u}}_i\right\}}_{i=1}^{\infty }$， ${\left\{{\bar{v}}_i\right\}}_{i=1}^{\infty }$， ${\left\{{\underset{\_}{v}}_i\right\}}_{i=1}^{\infty }$ 满足

$\{\begin{array}{l}0<{\underset{\_}{u}}_1\le \cdots \le {\underset{\_}{u}}_i\le {\underset{\_}{u}}_{i+1}\le \cdots \le {\bar{u}}_{i+1}\le {\bar{u}}_i\le \cdots \le {\bar{u}}_1,\\ 0<{\underset{\_}{v}}_1\le \cdots \le {\underset{\_}{v}}_i\le {\underset{\_}{v}}_{i+1}\le \cdots \le {\bar{v}}_{i+1}\le {\bar{v}}_i\le \cdots \le {\bar{v}}_1.\end{array}$

2) 收敛性：定义 ${\bar{u}}_{\infty }=\underset{i\to \infty }{\lim }{\bar{u}}_i$， ${\underset{\_}{u}}_{\infty }=\underset{i\to \infty }{\lim }{\underset{\_}{u}}_i$， ${\bar{v}}_{\infty }=\underset{i\to \infty }{\lim }{\bar{v}}_i$， ${\underset{\_}{v}}_{\infty }=\underset{i\to \infty }{\lim }{\underset{\_}{v}}_i$ 满足

$\{\begin{array}{l}\bar{a}{\bar{u}}_{\infty }-{\bar{u}}_{\infty }^2-\frac{b{\bar{u}}_{\infty }{\underset{\_}{v}}_{\infty }}{m+r{\bar{u}}_{\infty }}=0,\\ \underset{\_}{a}{\underset{\_}{u}}_{\infty }-{\underset{\_}{u}}_{\infty }^2-\frac{b{\underset{\_}{u}}_{\infty }{\bar{v}}_{\infty }}{m+r{\underset{\_}{u}}_{\infty }}=0,\\ b{\bar{v}}_{\infty }-\frac{b{\bar{v}}_{\infty }^2}{r{\bar{u}}_{\infty }}=0,\\ b{\underset{\_}{v}}_{\infty }-\frac{b{\underset{\_}{v}}_{\infty }^2}{r{\underset{\_}{u}}_{\infty }}=0.\end{array}$ (4.3)

证明 (1)根据定理1 (3.1)式以及迭代序列(4.2)可得

$\{\begin{array}{l}{\underset{\_}{u}}_1\le {\underset{\_}{u}}_2\le {\bar{u}}_2\le {\bar{u}}_1,\\ {\underset{\_}{v}}_1\le {\underset{\_}{v}}_2\le {\bar{v}}_2\le {\bar{v}}_1.\end{array}$ (4.4)

$\{\begin{array}{l}\bar{a}{\bar{u}}_1-{\bar{u}}_1^2-\frac{b{\bar{u}}_1{\underset{\_}{v}}_1}{m+r{\bar{u}}_1}\le 0,\\ \underset{\_}{a}{\underset{\_}{u}}_1-{\underset{\_}{u}}_1^2-\frac{b{\underset{\_}{u}}_{\infty }{\bar{v}}_1}{m+r{\underset{\_}{u}}_1}\ge 0.\end{array}$

$\{\begin{array}{l}b{\bar{v}}_1-\frac{b{\bar{v}}_1^2}{r{\bar{u}}_1}\le 0\\ b{\underset{\_}{v}}_1-\frac{b{\underset{\_}{v}}_1^2}{r{\underset{\_}{u}}_1}\ge 0\end{array}$

根据(3.3)式有 $\underset{\_}{c}\le \bar{c}$，结合迭代序列(4.1)

利用数学归纳法

当 $i=1$ 时

${\bar{u}}_2-{\underset{\_}{u}}_2\ge {\bar{u}}_1-{\underset{\_}{u}}_1+\frac{1}{T}\left(\bar{a}{\bar{u}}_1-{\bar{u}}_1^2-\frac{b{\bar{u}}_1{\underset{\_}{v}}_1}{m+r{\bar{u}}_1}\right)-\frac{1}{T}\left(\underset{\_}{a}{\underset{\_}{u}}_1-{\underset{\_}{u}}_1^2-\frac{b{\underset{\_}{u}}_1{\bar{v}}_1}{m+r{\underset{\_}{u}}_1}\right)\ge 0$

同理有 ${\bar{v}}_2-{\underset{\_}{v}}_2$，因此(4.4)式是成立的。

假设 $i=N$ 时也成立，当 $i=N+1$ 时有

${\bar{u}}_{N+1}-{\underset{\_}{u}}_{N+1}\ge {\bar{u}}_N-{\underset{\_}{u}}_N+\frac{1}{T}\left(\bar{a}{\bar{u}}_N-{\bar{u}}_N^2-\frac{b{\bar{u}}_N{\underset{\_}{v}}_N}{m+r{\bar{u}}_N}\right)-\frac{1}{T}\left(\underset{\_}{a}{\underset{\_}{u}}_N-{\underset{\_}{u}}_N^2-\frac{b{\underset{\_}{u}}_N{\bar{v}}_N}{m+r{\underset{\_}{u}}_N}\right)\ge 0$

${\bar{u}}_{N+1}-{\bar{u}}_N\le {\bar{u}}_N-{\bar{u}}_{N-1}+\frac{1}{T}\left(\bar{a}{\bar{u}}_N-{\bar{u}}_N^2-\frac{b{\bar{u}}_N{\underset{\_}{v}}_{N-1}}{m+r{\bar{u}}_N}\right)-\frac{1}{T}\left(\bar{a}{\bar{u}}_{N-1}-{\bar{u}}_{N-1}^2-\frac{b{\bar{u}}_{N-1}{\underset{\_}{v}}_{N-1}}{m+r{\bar{u}}_{N-1}}\right)\le 0$

成立。同样可以得到：

${\underset{\_}{u}}_{i+1}\ge {\underset{\_}{u}}_i$ , ${\underset{\_}{v}}_{i+1}\ge {\underset{\_}{v}}_i$ , ${\bar{v}}_{i+1}\ge {\bar{v}}_{i+1}$ , ${\bar{v}}_{i+1}\le {\bar{v}}_i.$

单调性证明完毕。

(2)在证明(1)中的公式已经表明序列 ${\left\{{\bar{u}}_i\right\}}_{i=1}^{\infty }$， ${\left\{{\underset{\_}{u}}_i\right\}}_{i=1}^{\infty }$， ${\left\{{\bar{v}}_i\right\}}_{i=1}^{\infty }$， ${\left\{{\underset{\_}{v}}_i\right\}}_{i=1}^{\infty }$ 是分别收敛于某个常数的，因此收敛性也得到了证明。

本节得出了问题(1.1)正稳态解的下界，通过迭代过程，得到了问题(1.1)正稳态解的更精确估计。

5. 异质环境下扩散Holling-Tanner模型正稳态解的全局稳定性

对于资源分布均匀的捕食–食饵系统，可以采用上下解方法或者构造李雅普诺夫函数来说明正稳态解的全局稳定性 [5] [6] [7] [8] [9]，然而，这并不适用于环境异质的情况。本文采用一种新的李亚普诺夫函数形式： $F\left(t\right)={\displaystyle \underset{\Omega }{\int }{\omega }_{*}}\left(x\right)V\left(u_i\left(x,t\right)\right)\text{d}x$，其中 ${\omega }_i\left(x\right)$ 是权函数。此函数对探索更一般的非同质环境下的扩散捕食–食饵模型 [10] [11] 也是一种有力的工具。这种李雅普诺夫函数也被用于研究扩散SIR传染病模型 [12]。对于模型(1.1)，构造如下收敛的李雅普诺夫函数：

$E\left(t\right):={\displaystyle \underset{\Omega }{\int }{\displaystyle {\int }_{u_{*}\left(x\right)}^{u\left(x,t\right)}\frac{u_{*}\left(x\right)}{d_1\left(x\right)}}}\frac{s-u_{*}\left(x\right)}{s}\text{d}s\text{d}x+\xi {\displaystyle \underset{\Omega }{\int }{\displaystyle {\int }_{v_{*}\left(x\right)}^{v\left(x,t\right)}\frac{v_{*}\left(x\right)}{d_2\left(x\right)}}}\frac{s-v_{*}\left(x\right)}{s}\text{d}s\text{d}x$

定理3 假设在 $\bar{\Omega }$ 上，对任意 $0<\alpha <1$， $a\left(x\right)>0,d_i\left(x\right)>0$，我们定义

$\overline{d_i}={\max }_{x\in \bar{\Omega }}{d}_i\left(x\right)$， $\underset{\_}{d_i}={\min }_{x\in \bar{\Omega }}{d}_i\left(x\right).$

若

$\xi =\sqrt{\frac{\underset{\_}{d_2}\overline{d_2}\left({\underset{\_}{u}}_{\infty }-\epsilon \right)\left({\bar{u}}_{\infty }-\epsilon \right){\underset{\_}{u}}_{\infty }^3{\bar{u}}_{\infty }}{\underset{\_}{d_1}\overline{d_1}{\underset{\_}{v}}_{\infty }{\bar{v}}_{\infty }^3}}\frac{1}{m+r\left({\underset{\_}{u}}_{\infty }-\epsilon \right)}$

那么方程(1.1)有一个唯一的正稳态解 $\left(u_{*},v_{*}\right)$，且正稳态解 $\left(u_{*},v_{*}\right)$ 是全局稳定的。其中 $\underset{t\to \infty }{\lim }u\left(x,t\right)=u_{*}\left(x\right)$， $\underset{t\to \infty }{\lim }v\left(x,t\right)=v_{*}\left(x\right)$。

证明 设 $\left(u\left(x,t\right),v\left(x,t\right)\right)$ 是方程(1.1)的解，定义一个函数

$E\left(t\right):={\displaystyle \underset{\Omega }{\int }{\displaystyle {\int }_{u_{*}\left(x\right)}^{u\left(x,t\right)}\frac{u_{*}\left(x\right)}{d_1\left(x\right)}}}\frac{s-u_{*}\left(x\right)}{s}\text{d}s\text{d}x+\xi {\displaystyle \underset{\Omega }{\int }{\displaystyle {\int }_{v_{*}\left(x\right)}^{v\left(x,t\right)}\frac{v_{*}\left(x\right)}{d_2\left(x\right)}}}\frac{s-v_{*}\left(x\right)}{s}\text{d}s\text{d}x$

其中 $\xi >0$，且当 $t\ge 0$ 时 $E\left(t\right)\ge 0$，根据格林公式

$\begin{array}{c}\frac{\text{d}E\left(t\right)}{\text{d}t}={\displaystyle \underset{\Omega }{\int }\frac{u_{*}\left(u-u_{*}\right)}{d_{1}u}{u}_t\text{d}x}+\xi {\displaystyle \underset{\Omega }{\int }\frac{v_{*}\left(v-v_{*}\right)}{d_{1}v}{v}_t\text{d}x}\\ ={\displaystyle \underset{\Omega }{\int }\frac{u_{*}\left(u-u_{*}\right)}{d_{1}u}\left(d_1\Delta u+f\left(a\left(x\right),u,v\right)-\frac{u}{u_{*}}{d}_1\Delta u_{*}-\frac{u}{u_{*}}f\left(a\left(x\right),u_{*},v_{*}\right)\right)\text{d}x}\\ +\xi {\displaystyle \underset{\Omega }{\int }\frac{v_{*}\left(v-v_{*}\right)}{d_{2}v}\left(d_2\Delta v+g\left(u,v\right)-\frac{v}{v_{*}}{d}_2\Delta v_{*}-\frac{v}{v_{*}}g\left(u_{*},v_{*}\right)\right)\text{d}x}\end{array}$

$\begin{array}{l}={\displaystyle \underset{\Omega }{\int }\left[\frac{u_{*}\left(u-u_{*}\right)}{u}\left(\Delta u-\frac{u}{u_{*}}\Delta u_{*}\right)+\frac{u_{*}\left(u-u_{*}\right)}{d_1}\left(\frac{f\left(a\left(x\right),u,v\right)}{u}-\frac{f\left(a\left(x\right),u_{*},v_{*}\right)}{u_{*}}\right)\right]\text{d}x}\\ +\xi {\displaystyle \underset{\Omega }{\int }\left[\frac{v_{*}\left(v-v_{*}\right)}{v}\left(\Delta v-\frac{v}{v_{*}}\Delta v_{*}\right)+\frac{v_{*}\left(v-v_{*}\right)}{d_2}\left(\frac{g\left(u,v\right)}{v}-\frac{g\left(u_{*},v_{*}\right)}{v_{*}}\right)\right]\text{d}x}\\ \le {\displaystyle \underset{\Omega }{\int }\frac{u_{*}\left(u-u_{*}\right)}{d_1}\left(\frac{f\left(a\left(x\right),u,v\right)}{u}-\frac{f\left(a\left(x\right),u_{*},v_{*}\right)}{u_{*}}\right)\text{d}x}\\ +\xi {\displaystyle \underset{\Omega }{\int }\frac{v_{*}\left(v-v_{*}\right)}{d_2}\left(\frac{g\left(u,v\right)}{v}-\frac{g\left(u_{*},v_{*}\right)}{v_{*}}\right)\text{d}x}\end{array}$

其中 $f\left(a\left(x\right),u,v\right):=a\left(x\right)u-u^2-\frac{buv}{m+ru}$， $g\left(u,v\right):=bv-\frac{b{v}^2}{\gamma u}$。

$\begin{array}{c}\frac{\text{d}E\left(t\right)}{\text{d}t}\le {\displaystyle \underset{\Omega }{\int }\frac{u_{*}\left(u-u_{*}\right)}{d_1}\left(-u-\frac{bv}{m+ru}+u_{*}+\frac{b{v}_{*}}{m+r{u}_{*}}\right)\text{d}x}\\ +\xi {\displaystyle \underset{\Omega }{\int }\frac{v_{*}\left(v-v_{*}\right)}{d_2}b\left(-\frac{v}{\gamma u}+\frac{v}{\gamma u_{*}}\right)\text{d}x}\\ ={\displaystyle \underset{\Omega }{\int }\frac{J_1{\left(u-u_{*}\right)}^2+J_2\left(u-u_{*}\right)\left(v-v_{*}\right)+J_3{\left(v-v_{*}\right)}^2}{d_1{d}_2\gamma u{u}_{*}\left(m+ru\right)\left(m+r{u}_{*}\right)}\text{d}x}\end{array}$

其中

$J_1:=-d_2\gamma u{u}_{*}^2\left[\left(m+ru\right)\left(m+r{u}_{*}\right)-br{v}_{*}\right],$

$J_2:=-b{d}_2\gamma u{u}_{*}^2\left(m+r{u}_{*}\right)+d_1\xi b{v}_{*}^2\left(m+ru\right)\left(m+r{u}_{*}\right),$

$J_3:=-d_1\xi b{u}_{*}{v}_{*}\left(m+ru\right)\left(m+r{u}_{*}\right),$

当 $\xi =\sqrt{\frac{\underset{\_}{d_2}\overline{d_2}\left({\underset{\_}{u}}_{\infty }-\epsilon \right)\left({\bar{u}}_{\infty }-\epsilon \right){\underset{\_}{u}}_{\infty }^3{\bar{u}}_{\infty }}{\underset{\_}{d_1}\overline{d_1}{\underset{\_}{v}}_{\infty }{\bar{v}}_{\infty }^3}}\frac{1}{m+r\left({\underset{\_}{u}}_{\infty }-\epsilon \right)}$ 时，存在 $0<\delta \le 1$，使得 $B^2<4AC$，有

$\frac{\text{d}E\left(t\right)}{\text{d}t}\le -{\displaystyle \underset{\Omega }{\int }\frac{\delta {\left(u-u_{*}\right)}^2+\delta {\left(v-v_{*}\right)}^2}{d_1{d}_2\gamma u{u}_{*}\left(m+ru\right)\left(m+r{u}_{*}\right)}\text{d}x}=:\psi \left(t\right)\le 0$

根据(3.2)式，存在 $G_1>0$，使得 $\left|\psi \left(t\right)\right|<G_1$。利用引理2有：

$\underset{t\to \infty }{\lim }\psi \left(t\right)=-{\displaystyle \underset{\Omega }{\int }\frac{\delta {\left(u-u_{*}\right)}^2+\delta {\left(v-v_{*}\right)}^2}{d_1{d}_2\gamma u{u}_{*}\left(m+ru\right)\left(m+r{u}_{*}\right)}\text{d}x}=0$

由(3.1)式，存在 $t\ge T_1$，有 $0<{\epsilon }_2\le u\left(x,t\right)$，得到

$\underset{t\to \infty }{\lim }u\left(x,t\right)=u_{*}\left(x\right)$ , $\underset{t\to \infty }{\lim }v\left(x,t\right)=v_{*}\left(x\right).$ (5.1)

同时，(3.2)式也表明集合 $\left\{u\left(\cdot ,t\right):t\ge 1\right\}$ 在 $C^2\left(\bar{\Omega }\right)$ 上是相对紧凑的，可以假设

${\Vert u\left(x,t_p\right)-\tilde{u}\left(x\right)\Vert }_{C^2\left(\bar{\Omega }\right)}$， ${\Vert v\left(x,t_p\right)-\tilde{v}\left(x\right)\Vert }_{C^2\left(\bar{\Omega }\right)}\to 0$， $t_p\to \infty$

可以得出，在 $\bar{\Omega }$ 上有

$\tilde{u}\left(x\right)\equiv u_{*}\left(x\right)$ , $\tilde{v}\left(x\right)\equiv v_{*}\left(x\right).$

因此， $\underset{t\to \infty }{\lim }u\left(x,t\right)=u_{*}\left(x\right)$， $\underset{t\to \infty }{\lim }v\left(x,t\right)=v_{*}\left(x\right)$ 在 $C^2\left(\bar{\Omega }\right)$ 上一致成立，那么正稳态解 $\left(u_{*},v_{*}\right)$ 是全局稳定的，证毕。

6. 结束语

针对Holling-Tanner捕食–食饵模型，本文研究了资源空间分散和资源分布不均对种群的影响。本文首先利用抛物方程的比较定理和上下解方法证明了方程(1.1)正稳态解的存在性，进一步通过迭代过程，得出了问题(1.1)正稳态解的下界。最后利用新构造的李雅普诺夫函数，得到Holling-Tanner捕食–食饵模型在异质环境下的正稳态解是全局稳定的。

基金项目

国家自然科学基金青年项目(52004215)；陕西省自然科学基金青年项目(2019JQ-755)；陕西省教育厅基金(19JK0462)。