Lipschitz不等约束下不确定分式优化问题鲁棒最优解集的刻画

Some Characterizations of Robust Optimal Solution Set for Uncertain Fractional Optimization Problem with Lipschitz Inequality Constraints

摘要: 本文考虑一类在目标和约束函数中都含有不确定数据的分式优化问题。其中，目标函数是分式函数，约束是满足Lipschitz条件的不等式。利用伪Lagrange型函数给出了鲁棒最优解集的刻画。

Abstract: In this paper, we deal with a class of fractional optimization problem with data uncertainty in both the objective and constraints. The objective function is fractional function and the constraint is inequalities which satisfy Lipschitz conditions. Using pseudo Lagrange-type function, we give those characterizations of robust optimal solution set of uncertain fractional optimization problem.

1. 引言

分式优化问题是最优化理论研究的一类重要问题，它在经济学，金融学，聚类分析及排队选址等问题中有广泛的应用 [1] [2] [3]。然而在实际问题中，由于预测或测量误差，会导致问题输入数据不完整或不确定 [4] [5] [6]，我们需要在知道真实数据和参数前做出决策。因此，含有不确定数据的优化问题引起了研究者的广泛关注。鲁棒优化方法是解决不确定优化问题一个有效方法，见文献 [7] - [16]。

解集特征是不确定规划问题的一个重要研究方向。关于解集特征概念的介绍和研究是由Mangasarian [9] 对可微凸问题提出的，关于凸问题下的鲁棒解集特征可参见文献 [10] [11] [12]。近年来最优问题的解集特征刻画被推广到分式优化问题的最优解集特征的刻画，见文献 [13] [14] [15] [17] [18]。

文献 [13] [14] 在可微的情况下，讨论了凸约束下的鲁棒分式规划问题，得到鲁棒对偶性结论。Sun等 [15] 在不可微的情况下，对目标函数含有不确定数据的凸–凹分式函数，约束函数是含有不确定数据连续的凸–凹函数，利用鲁棒型基本次微分约束规则，给出鲁棒最优解集的刻画；Sisarat，Wangkeeree和Lee [16] 在不可微情况下，考虑了目标函数是不确定凸函数，约束函数是含有不确定数据的Lipschitz函数的优化问题，利用鲁棒型基本约束规则对问题的伪Lagrange型函数和鲁棒最优解集的性质进行了研究。

受文献 [15] [16] 启发，本文主要考虑目标函数是凸–凹分式函数，约束函数是Lipschitz连续函数的问题，在目标函数与约束函数都含有不确定数据且只要求约束集是凸集情况下，利用鲁棒次微分约束规则(RSCQ)对问题的鲁棒最优解集进行刻画。本文的结论推广了文献 [15] [16] 中的相关定理。

文章结构如下：第二部分，介绍基本概念和相关符号；第三部分，利用鲁棒次微分约束规则(RSCQ)对分式优化问题的伪Lagrange型函数和鲁棒最优解集特征进行刻画。

2. 预备知识

在这一节中，我们先回顾几个基本概念和结论，并给出本文要用到的若干引理。

设C是 $R^n$ 的子集，C的指标函数 ${\delta }_c:R^n\to R$，定义为：

${\delta }_c:=\{\begin{array}{ll}0\hfill & x\in C,\hfill \\ +\infty \hfill & x\notin C.\hfill \end{array}$

若 $f:R^n\to R$ 对 $\forall x,y\in X\subseteq R^n$， $t\in \left[0,1\right]$ 满足： $f\left(tx+\left(1-t\right)y\right)\le tf\left(x\right)+\left(1-t\right)f\left(y\right)$ 称f为凸函数。当−f是凸函数，称f为凹函数。f在 $\bar{x}$ 处的凸次微分定义为：

$\partial f\left(\bar{x}\right):=\left\{x^{*}\in R^n:\langle x^{*},x-\bar{x}\rangle \le f\left(x\right)-f\left(\bar{x}\right),\forall x\in X\right\}.$

f在 $x\in R^n$ 处沿方向 $d\in R^n$ 的方向导数表示为：

$f^{\prime }\left(x;\cdot \right):=\underset{t\to 0^{+}}{\lim }\frac{f\left(x+td\right)-f\left(x\right)}{t}.$

定义2.1 [16] 称函数 $h:R^n\to R$ 在 $x\in R^n$ 是Lipschitz连续的，若存在 $L>0$ 和x的邻域N满足：

$\left|h\left(y\right)-h\left(z\right)\right|\le L\Vert y-z\Vert ,\forall x,y\in N.$

定义2.2 [19] 设函数 $h:R^n\to R$ 在 $x\in R^n$ 是Lipschitz连续的，h在x处沿方向d的广义Clarke方向导数记为 $h^{\circ }\left(x;d\right)$，定义为：

$h^{\circ }\left(x;d\right):=\underset{\begin{array}{l}y\to x\\ t\to 0^{+}\end{array}}{\lim \sup }\frac{h\left(y+td\right)-h\left(y\right)}{t}.$

定义2.3 [19] 设函数 $h:R^n\to R$ 在 $x\in R^n$ 是Lipschitz连续的，h在x处的广义Clarke次微分记为 ${\partial }^{\circ }h\left(x\right)$，定义为：

${\partial }^{\circ }h\left(x\right):=\left\{\xi \in R^n:h^{\circ }\left(x;d\right)\ge \langle \xi ,d\rangle ,\forall d\in R^n\right\}.$

定义2.4 [16] 设函数 $h:R^n\to R$ 在 $x\in R^n$ 处是Lipschitz连续的，若对每个方向 $d\in R^n$，方向导数 $h^{\prime }\left(x;d\right)$ 存在且等于 $h^{\circ }\left(x;d\right)$，则称h在 $x\in R^n$ 处正则。

本文考虑如下分式优化问题：

$\left(\text{UFP}\right)\underset{x\in R^n}{\min }\left\{\frac{f\left(x,u\right)}{g\left(x,v\right)},x\in C:h_j\left(x,w_j\right)\le 0,\forall w_j\in W_j,j=1,2,\cdots ,m\right\},$

其中 $f:R^n\times R^p\to R$， $g:R^n\times R^p\to R^{+}$， $h_j:R^n\times R^q\to R$ 是给定函数， $C\subseteq R^n$ 是非空凸闭集。 $U\subseteq R^p$， $V\subseteq R^p$， $W_j\subseteq R^q$ 为紧凸的不确定集， $u\in U$， $v\in V$， $w_j\in W_j$ 是不确定参数。(UFP)的鲁棒可行集记为：

$F:=\left\{x\in C,h_j\left(x,w_j\right)\le 0,w_j\in W_j,j=1,2,\cdots ,m\right\}$ .

对于给定紧子集 $W_j\subseteq R^q$ 和给定函数 $h_j:R^n\times R^q\to R$，对 $h_j$ 作如下假设 [16]：

(C₁)对 $\forall x\in R^n$， $w_j\in W_j$，函数 $w_j\mapsto h_j\left(x;w_j\right)$ 是上半连续。

(C₂) $h_j$ 关于 $w_j$ 在x处是一致Lipschitz连续的。

(C₃)对 $\forall \left(x,w_j\right)\in R^n\times W_j$，函数 $h_j\left(\cdot ,w_j\right)$ 是在x处正则。

(C₄)对于 $\left(x,w_j\right)\in R^n\times W_j$，集值映射 $\left(x,w_j\right)\mapsto {\partial }^{\circ }{h}_j\left(\cdot ,w_j\right)\left(x\right)$ 是上半连续。

$J\left(x\right):=\left\{j\in \left\{1,2,\cdots ,m\right\}:\underset{w_j\in W_j}{\max }{h}_j\left(x,w_j\right)=0,\forall x\in F\right\}\ne \varnothing .$

$W_j\left(x\right):=\left\{{\bar{w}}_j\in W_j:h_j\left(x,{\bar{w}}_j\right)=\underset{w_j\in W_j}{\max }{h}_j\left(x,w_j\right),j=1,2,\cdots ,m\right\}\ne \varnothing .$

文章的剩余部分除非特别说明，否则下面假设一直成立。 $f\left(\cdot ,u\right)$ 对 $\forall u\in U$ 是连续凸函数， $f\left(x,\cdot \right)$ 对 $\forall x\in R^n$ 是凹函数； $g\left(x,\cdot \right)$ 对 $\forall x\in R^n$ 是凸函数， $g\left(\cdot ,v\right)$ 对 $\forall v\in V$ 是连续凹函数，且令 $f\ge 0,g>0$，可行集F是非空凸集。

与(UFP)对应的鲁棒优化模型(RUFP)：

$\left(\text{RUFP}\right)\underset{x\in R^n}{\min }\left\{\frac{{\max }_{u\in U}f\left(x,u\right)}{{\min }_{v\in V}g\left(x,v\right)},x\in C:h_j\left(x,w_j\right)\le 0,\forall w_j\in W_j,j=1,2,\cdots ,m\right\}.$

定义2.5若 $\bar{x}$ 是(RUFP)的最优解，则 $\bar{x}\in F$ 称是(UFP)的鲁棒最优解。(UFP)鲁棒最优解集记为：

$S:=\left\{x\in F:\frac{{\max }_{u\in U}f\left(x,u\right)}{{\min }_{v\in V}g\left(x,v\right)}\le \frac{{\max }_{u\in U}f\left(y,u\right)}{{\min }_{v\in V}g\left(y,v\right)},\forall y\in F\right\},$

且假设 $S\ne \varnothing$。

定义2.6 [10] 设 $x\in F$，称鲁棒次微分约束品性(RSCQ)在x处成立，若

$\partial {\delta }_F\left(x\right)\subseteq \partial {\delta }_C\left(x\right)+{\displaystyle \underset{\begin{array}{c}{\lambda }_j\ge 0,w_j\in W_j,\\ {\lambda }_j{h}_j\left(x,w_j\right)=0,j=\left\{1,2,\cdots ,m\right\}\end{array}}{\cup }{\displaystyle \underset{j=1}{\overset{m}{\sum }}{\lambda }_j{\partial }^{\circ }{h}_j\left(\cdot ,w_j\right)\left(x\right)}}.$

引理2.1 [12] 令 $U\subseteq R^p$ 是凸紧集，且令 $f:R^n\times R^p\to R$ 是凸–凹函数，即对于固定的 $u\in U$， $f\left(\cdot ,u\right)$ 是关于 $x\in R^n$ 的凸函数，且对于固定的 $x\in R^n$， $f\left(x,\cdot \right)$ 是关于 $u\in U$ 的凹函数，则

$\partial \left(\underset{u\in U}{\max }f\left(\cdot ,u\right)\right)\left(\bar{x}\right)={\displaystyle \underset{u\in U\left(\bar{x}\right)}{\cup }\partial }f\left(\cdot ,u\right)\left(\bar{x}\right),$

其中 $U\left(\bar{x}\right):=\left\{\bar{u}\in U:f\left(\bar{x},\bar{u}\right)={\max }_{u\in U}f\left(\bar{x},u\right)\right\}$。

引理2.2 [16] 对 $\forall x\in F,\forall j\in J\left(x\right)$， $h_j$ 满足(C₃)，(C₄)，若F是凸集，则：

$F\subseteq \left\{y\in R^n:h_{jx}^{\circ }\left(x,w_j;y-x\right)\le 0,\forall x\in F,\forall j\in J\left(x\right),\forall w_j\in W_j\left(x\right)\right\}.$

3. 鲁棒最优解集的刻画

(UFP)的伪Lagrange型函数 [16] 定义为：

$L^p\left(x,y,\mu ,{\lambda }_j,u,v,w_j\right)=f\left(x,u\right)+\mu \left(-g\right)\left(x,v\right)+{\displaystyle \underset{j\in J\left(y\right)}{\sum }{\lambda }_j}{h}_{jx}^{\circ }\left(y,w_j;x-y\right),x\in R^n.$

其中， $\mu =\frac{{\max }_{u\in U}f\left(x,u\right)}{{\min }_{v\in V}g\left(x,v\right)},y\in R^n,{\lambda }_j\ge 0,u\in U,v\in V,w_j\in W_j,j=1,2,\cdots ,m$。

定理3.1设 $\bar{x}\in S$ 是(UFP)的鲁棒最优解，(RSCQ)在 $\bar{x}$ 处成立，设 $\bar{\mu }=\frac{{\max }_{u\in U}f\left(\bar{x},u\right)}{{\min }_{v\in V}g\left(\bar{x},v\right)}$，则 $\exists {\bar{\lambda }}_j\ge 0,\bar{u}\in U,\bar{v}\in V,{\bar{w}}_j\in W_j$ 和 $j=1,2,\cdots ,m$ 满足：

${\bar{\lambda }}_j{h}_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j;x-\bar{x}\right)=0,\forall j\in J\left(\bar{x}\right);\frac{f\left(x,\bar{u}\right)}{g\left(x,\bar{v}\right)}=\frac{{\max }_{u\in U}f\left(x,u\right)}{{\min }_{v\in V}g\left(x,v\right)},\forall x\in S;$

且 $L^p\left(\cdot ,\bar{x},\bar{\mu },{\bar{\lambda }}_j,\bar{u},\bar{v},{\bar{w}}_j\right)$ 在S上是常数。

证明： $\bar{x}\in S$ 是(UFP)的鲁棒最优解，则有 $\bar{x}$ 是(UFP)鲁棒最优解当且仅当 $\bar{x}$ 是如下问题的鲁棒最优解 [15]：

$\underset{x\in R^n}{\min }\left\{\underset{u\in U,v\in V}{\max }\left\{f\left(x,u\right)+\bar{\mu }\left(-g\right)\left(x,v\right)\right\}:x\in C;h_j\left(x,w_j\right)\le 0,\forall w_j\in W_j,j=1,2,\cdots ,m\right\}.$

对 $\forall x\in R^n$，令

$\varphi \left(x\right)=\underset{u\in U,v\in V}{\max }\left\{f\left(x,u\right)+\bar{\mu }\left(-g\right)\left(x,v\right)\right\}.$

因 $f\left(\cdot ,u\right)$ 对 $\forall u\in U$ 是连续凸函数， $g\left(\cdot ,v\right)$ 对 $\forall v\in V$ 是连续凹函数，故 $\varphi \left(x\right)$ 是连续凸函数。

由引理2.1知：

$\begin{array}{l}0\in \partial \left(\varphi +{\delta }_F\right)\left(\bar{x}\right)=\partial \varphi \left(\bar{x}\right)+\partial {\delta }_F\left(\bar{x}\right)\\ ={\displaystyle \underset{\left(\bar{u},\bar{v}\right)\in U^{\prime }\left(\bar{x}\right)\times V^{\prime }\left(\bar{x}\right)}{\cup }\left\{\partial f\left(\cdot ,u\right)\left(\bar{x}\right)+\bar{\mu }\left(-g\right)\left(\cdot ,v\right)\left(\bar{x}\right)\right\}}+\partial {\delta }_F(\; x\; ¯\; )\end{array}$

其中 $\left(\bar{u},\bar{v}\right)\in U^{\prime }\left(\bar{x}\right)\times V^{\prime }\left(\bar{x}\right):=\left\{\left(u,v\right)\in U\times V:f\left(\bar{x},u\right)+\bar{\mu }\left(-g\right)\left(\bar{x},v\right)=\varphi \left(\bar{x}\right)\right\}$。

由(RSCQ)在 $\bar{x}$ 处成立得， $\exists {\bar{\lambda }}_j\ge 0,\bar{u}\in U,\bar{v}\in V,{\bar{w}}_j\in W_j$ 满足：

$\begin{array}{c}0\in \partial f\left(\cdot ,u\right)\left(\bar{x}\right)+\bar{\mu }\left(-g\right)\left(\cdot ,v\right)\left(\bar{x}\right)+\partial {\delta }_C\left(\bar{x}\right)+{\displaystyle \underset{j\in J(\bar{x})}{\sum }{\bar{\lambda }}_j{\partial }^{\circ }{h}_j\left(\cdot ,{\bar{w}}_j\right)\left(\bar{x}\right)}\\ \in \partial f\left(\cdot ,u\right)\left(\bar{x}\right)+\bar{\mu }\left(-g\right)\left(\cdot ,v\right)\left(\bar{x}\right)+\partial {\delta }_C\left(\bar{x}\right)+{\displaystyle \underset{j\in J(\bar{x})}{\sum }{\bar{\lambda }}_j\partial h_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j;\cdot -\bar{x}\right)\left(\bar{x}\right)}\\ \subseteq \partial L^p\left(\cdot ,\bar{x},\bar{\mu },{\bar{\lambda }}_j,\bar{u},\bar{v},{\bar{w}}_j\right)\left(\bar{x}\right).\end{array}$ (1)

$f\left(\bar{x},\bar{u}\right)+\bar{\mu }\left(-g\right)\left(\bar{x},\bar{v}\right)=\varphi \left(\bar{x}\right)=\underset{u\in U,v\in V}{\max }\left\{f\left(\bar{x},u\right)+\bar{\mu }\left(-g\right)\left(\bar{x},v\right)\right\}.$ (2)

${\bar{\lambda }}_j{h}_j\left(\bar{x},{\bar{w}}_j\right)=0.$ (3)

由(2)式得 [15]：

$\frac{f\left(\bar{x},\bar{u}\right)}{g\left(\bar{x},\bar{v}\right)}=\frac{{\max }_{u\in U}f\left(\bar{x},u\right)}{{\min }_{v\in V}g\left(\bar{x},v\right)}=\bar{\mu }.$ (4)

由次微分定义和(1)，(2)式得：

$\begin{array}{l}f\left(\bar{x},\bar{u}\right)+\bar{\mu }\left(-g\right)\left(\bar{x},\bar{v}\right)+{\displaystyle \underset{j\in J\left(\bar{x}\right)}{\sum }{\bar{\lambda }}_j{h}_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j;x-\bar{x}\right)}\\ \ge f\left(\bar{x},\bar{u}\right)+\bar{\mu }\left(-g\right)\left(\bar{x},\bar{v}\right)\\ =\underset{u\in U,v\in V}{\max }\left\{f\left(\bar{x},u\right)+\bar{\mu }\left(-g\right)\left(\bar{x},v\right)\right\},\forall x\in R^n.\end{array}$ (5)

对 $\forall x,\bar{x}\in S$ 有：

$\underset{u\in U,v\in V}{\max }\left\{f\left(x,u\right)+\bar{\mu }\left(-g\right)\left(x,v\right)\right\}=\underset{u\in U,v\in V}{\max }\left\{f\left(\bar{x},u\right)+\bar{\mu }\left(-g\right)\left(\bar{x},v\right)\right\}.$ (6)

由(5)，(6)式得：

${\displaystyle \underset{j\in J\left(\bar{x}\right)}{\sum }{\bar{\lambda }}_j{h}_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j,x-\bar{x}\right)}\ge 0,\forall \bar{x}\in S.$ (7)

当 ${\bar{\lambda }}_j>0$ 时，由(3)式得：

$h_j\left(\bar{x},{\bar{w}}_j\right)=0.$

对于 $j\in J\left(\bar{x}\right)$ 结合上式得：

$h_j\left(\bar{x},{\bar{w}}_j\right)=0=\underset{w_j\in W_j}{\max }{h}_j\left(\bar{x},w_j\right).$

故， ${\bar{w}}_j\in {\bar{W}}_j\left(\bar{x}\right)$。

此时由引理2.2得，对于 $\forall j\in J\left(\bar{x}\right),\forall {\bar{w}}_j\in {\bar{W}}_j\left(\bar{x}\right),\forall x\in F$ 有：

${\bar{\lambda }}_j{h}_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j;x-\bar{x}\right)\le 0.$

由上式和(7)式得：

${\bar{\lambda }}_j{h}_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j;x-\bar{x}\right)=0,\forall j\in J\left(\bar{x}\right).$ (8)

下证：

$\frac{f\left(x,\bar{u}\right)}{g\left(x,\bar{v}\right)}=\frac{{\max }_{u\in U}f\left(x,u\right)}{{\min }_{v\in V}g\left(x,v\right)},\forall x\in S.$ (9)

对 $\forall x,\bar{x}\in S$ 得：

$\frac{f\left(x,\bar{u}\right)}{g\left(x,\bar{v}\right)}\le \frac{{\max }_{u\in U}f\left(x,u\right)}{{\min }_{v\in V}g\left(x,v\right)}=\frac{{\max }_{u\in U}f\left(\bar{x},u\right)}{{\min }_{v\in V}g\left(\bar{x},v\right)}=\bar{\mu }.$

由上式得：

$f\left(x,\bar{u}\right)+\bar{\mu }\left(-g\right)\left(x,\bar{v}\right)\le 0.$ (10)

对 $\forall x\in S$，结合(4)，(5)，(8)式得：

$f\left(x,\bar{u}\right)+\bar{\mu }\left(-g\right)\left(x,\bar{v}\right)\ge 0.$

由上式和(10)式得：

$f\left(x,\bar{u}\right)+\bar{\mu }\left(-g\right)\left(x,\bar{v}\right)=0.$ (11)

故，

$\frac{f\left(x,\bar{u}\right)}{g\left(x,\bar{v}\right)}=\bar{\mu }=\frac{{\max }_{u\in U}f\left(\bar{x},u\right)}{{\min }_{v\in V}g\left(\bar{x},v\right)}=\frac{{\max }_{u\in U}f\left(x,u\right)}{{\min }_{v\in V}g\left(x,v\right)}.$

由上式可知，(9)式成立。

对 $\forall x\in S$，由(8)式和(11)式得：

$\begin{array}{c}{L}^p\left(x,\bar{x},\bar{\mu },{\bar{\lambda }}_j,\bar{u},\bar{v},{\bar{w}}_j\right)=f\left(x,\bar{u}\right)+\bar{\mu }\left(-g\right)\left(x,\bar{v}\right)+{\displaystyle \underset{j\in J\left(\bar{x}\right)}{\sum }{\bar{\lambda }}_j}{h}_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j;x-\bar{x}\right)\\ =f\left(x,\bar{u}\right)+\bar{\mu }\left(-g\right)\left(x,\bar{v}\right)=0.\end{array}$

故 $L^p\left(\cdot ,\bar{x},\bar{\mu },{\bar{\lambda }}_j,\bar{u},\bar{v},{\bar{w}}_j\right)$ 在S上是常数。

注1：若 $h_j\left(\cdot ,w_j\right),j=1,2,\cdots ,m$ 对 $\forall w_j\in W_j$ 是凸函数，则对 $\forall j=1,2,\cdots ,m$，由定理3.1可得：

${\bar{\lambda }}_j{h}_j\left(x,{\bar{w}}_j\right)-{\bar{\lambda }}_j{h}_j\left(\bar{x},{\bar{w}}_j\right)\ge \bar{\lambda }{h^{\prime }}_{jx}\left(\bar{x},{\bar{w}}_j;x-\bar{x}\right)=\bar{\lambda }{h}_{jx}^{\circ }\left(\bar{x},{\bar{w}}_j;x-\bar{x}\right)=0,x\in S.$

对 $x\in F$，由上式和 ${\bar{\lambda }}_j{h}_j\left(\bar{x},{\bar{w}}_j\right)=0$， $j=1,2,\cdots ,m$，得：

且对有：

这就说明伪Lagrange型函数是文献 [15] 中Lagrange型函数。

注2：当关于是连续的凸–凹函数时，定理3.1推广了文献 [15] 中的命题2；当，定理3.1推广了文献 [16] 中命题2。

下面利用给定的(UFP)的一个鲁棒最优解对(UFP)的鲁棒最优解集进行刻画。

定理3.2 设是(UFP)的鲁棒最优解，(RSCQ)在处成立，设，则满足，其中，

证明：首先证明，令，则。由(1)式知及满足：

(12)

由次微分和广义次微分定义得：

(13)

(14)

(14)式两边同乘，联合(8)式得：

(15)

由(12)，(13)，(15)式得：

令，则，

且

(16)

对，由(4)式和(11)式得：

(17)

由(16)，(17)式得：

故，。

下证：。

事实上有：

由(17)式得：

因此，，从而：

故得证。

下证，对有，得：

(18)

由(4)式和(18)式得：

由上式得：

故对有：

此时由，，得：。

故得证。

注3：当关于是连续的凸–凹函数时，定理3.2推广了文献 [15] 中定理3.6；当，定理3.2推广了文献 [16] 中定理4.1。

基金项目

山西省高等学校科技创新项目(NO. 2019L0784)；山西省回国留学人员科研资助项目(NO. 2017-164)；太原师范学院研究生教育创新项目(NO. SYYJSJC-2018)。

NOTES

^*通讯作者。