Processing math: 1%

数学物理学报, 2023, 43(2): 515-530

一类带有弱连续算子的发展方程的最优控制

曾彪,

广西民族大学 数学与物理学院 南宁 530006; 广西应用数学中心 南宁 530006

Optimal Control for a Class of Nonlinear Evolutionary Equations with Weakly Continuous Operators

Zeng Biao,

College of of Mathematics and Physics, Center for Applied Mathematics of Guangxi, Nanning 530006; Guangxi University for Nationalities, Nanning 530006

收稿日期: 2021-08-30   修回日期: 2022-10-18  

基金资助: 国家自然科学基金(12201137)
广西科技基地和人才专项(桂科AD21220103)
广西民族大学引进人才科研启动项目(2019KJQD04)

Received: 2021-08-30   Revised: 2022-10-18  

Fund supported: NSFC(12201137)
Science and Technology Innovation Talent Project of Guangxi(桂科AD21220103)
Start-up Project of Scientific Research on Introducing Talents at School Level in Guangxi University for Nationalities(2019KJQD04)

作者简介 About authors

曾彪,E-mail:biao_zeng@163.com

摘要

该文研究一类带有弱连续算子的发展方程的一个最优控制问题. 通过运用Rothe方法和弱连续算子的一个满射定理, 建立方程的可解性. 然后证明最优控制问题的最优状态-控制对的存在性. 最后把主要结果应用到非平稳的Navier-Stokes-Voigt方程上.

关键词: 发展方程; 最优控制; 弱连续算子; Navier-Stokes-Voigt方程

Abstract

In the paper we study an optimal control problem for a class of nonlinear evolutionary matrixs involving weakly continuous operators. By exploiting the Rothe method and using a surjectivity result for weakly continuous operators, we establish the solvability for the matrix. Then we show the existence of optimal state-control pairs for the optimal control problem. The main results are applied to non-stationary Navier-Stokes-Voigt matrix.

Keywords: Evolutionary matrix; Optimal control; Weakly continuous operator; Navier-Stokes-Voigt matrix

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

本文引用格式

曾彪. 一类带有弱连续算子的发展方程的最优控制[J]. 数学物理学报, 2023, 43(2): 515-530

Zeng Biao. Optimal Control for a Class of Nonlinear Evolutionary Equations with Weakly Continuous Operators[J]. Acta Mathematica Scientia, 2023, 43(2): 515-530

1 引言

众所周知, 如果没有控制系统, 就不可能有制造业、车辆、计算机和环境调节. 近年来, 控制理论在工程、经济、计算机和生态学中得到了广泛的应用, 尤其是在最优控制系统中. 在现代科学和工程计算中, 最优控制问题在各种实际问题中经常遇到. 发展方程的最优控制问题在科学和工程计算中有着广泛的应用. 在过去的几十年中, 对发展方程最优控制问题的研究引起了数学家的极大兴趣. 在许多文献中致力于发展这类问题的数学理论和数值分析. 例如, Ahmed[1-3] 和Papageorgiou[4-6]研究了由发展方程和包含构成的最优控制系统的存在性. 此外, Cesari[7,8]利用不动点技术和抽象算子方程建立了拟线性无限维控制系统解和最优解的存在性. 更多细节, 参见文献[9-12]及其里面的相关文献.

近年来, 非线性发展方程作为在许多物理系统和电磁现象的动力学建模中是一个有用的工具, 已经被数学、科学和工程领域的许多研究人员呈现出来了(参见文献[13-17]). 与线性发展方程相比, 非线性发展方程可以为许多复杂系统提供更真实、更精确的描述. 尽管在线性发展方程的最优控制问题方面已经做了大量的工作, 但在非线性情况还有许多工作需要做. 因此, 建立研究非线性问题的方法是合理和有意义的. 近年来, 非线性发展方程构成的最优控制问题受到越来越多的关注 (参见文献[18-23]). 应该提到的是, 空间的无限维特性使得与非线性发展方程相关的最优控制问题与线性情况相比更难处理. 基于这些事实, 该文主要研究一类非线性发展方程的最优控制问题.

考虑发展三元组框架 (V,H,V),即, V 是一个可分自反Banach空间,H 是一个可分Hilbert空间,嵌入 VH 连续,且 VH 中稠密.此外, 假设嵌入 VH 是紧的. 对于 0<T<,令

V=L2(0,T;V),  V=L2(0,T;V),  W=W1,2(0,T;V)={vV:vV}.

显然 W 是一个Banach空间, 其范数为 , u\in {\cal W}.另外, 令 Y 为一个可分自反Banach空间, {\cal Y}=L^2(0,T;Y), {\cal Y}^*=L^2(0,T;Y^*).

给定两个非线性算子 A,B \colon V \to V^*, 一个线性有界算子 C \colon Y \to V^*, 及 f\in {\cal V}^*,u_0\in V, 考虑如下发展问题.

问题1.1 对于 y\in {\cal Y}, 寻找 u\in {\cal W} 使得

\left\{\begin{array}{llll}Au'(t)+Bu(t)= f(t)+Cy(t)\qquad \mbox{ a.e.} \ \ t\in (0,T), \\u(0) = u_0.\end{array}\right.

在这个工作中, 处理问题1.1解的存在性, 并给出其最优控制问题的结果. 主要创新之处如下.

首先, 为了得到非线性发展方程解的存在性, 考虑具有弱连续算子的方程, 而不是用文献中广泛运用的线性、单调性和伪单调性. 尽管弱连续算子不像 M 类算子那样一般, 它也覆盖很多技术上重要的情形(参见文献[24,25]).第二, 考虑问题1.1中参量比其它文献(参见文献[26,27])更少限制的假设. 例如, 不需要算子 B 的强制性假设.第三, 利用具有足够小的时间步长的Rothe方法来建立方程解的存在性, 这允许忽略关于常数的小条件. 第四, 我们将抽象的存在性结果应用于研究非平稳的Navier-Stokes-Voigt方程.

本文结构如下. 第2节回忆一些预备知识并研究发展方程解的存在性.在这里, 定理2.2是主要定理, 它的证明是涉及弱连续算子的静态方程基于Rothe方法的满射结果.第3节考虑一个最优控制问题. 第4节把前面的主要结果应用到给出非平稳的Navier-Stokes-Voigt方程的存在性结果.

2 存在性结果

首先介绍一些有用的预备知识.

(X, \|\cdot\|_{X}) 为一个Banach空间. 记 X^* 为其对偶空间, 并记 \langle\cdot,\cdot\rangle_{X^*,X}X^*X 之间的对偶积. 给定 0 < T < \infty, 令 L^2(0,T; X) 为从 (0,T)X 的所有平方可积函数构成的空间, C([T];X) 为从 [T]X 的所有连续函数构成的空间, 它们的范数分别为

\begin{matrix} \|x\|_{L^2} =\bigg{(}\int_{0}^{T}\|x(t)\|_{X}^2{ d}t \bigg{)}^{\frac{1}{2}}, \ \ \ \ \|x\|_{C} = \sup_{t\in [T]}\|x(t)\|_{X}. \end{matrix}

本文将利用如下弱连续算子的概念及其满射定理.

定义2.1[24,25]X,Y 为Banach空间. 一个算子 F \colon X \to Y 叫做弱连续的, 如果对于任何序列 \{u_{n}\}_{n\geq 1}\subset X 满足 u_{n}\rightharpoonup uX 中, 有 Fu_{n}\rightharpoonup FuY 中.

定义2.2X 为一个Banach空间. 一个算子 F\colon X \to X^* 叫做强制的, 如果

\begin{matrix}\lim_{\|u\|_X\rightarrow \infty}\frac{ \langle Fu,u\rangle_{X^*,X}}{\|u\|_X} = \infty.\end{matrix}

定理2.1[24]X 为一个可分自反的Banach空间, F : X \to X^* 为弱连续和强制的. 那么 F 是满射的, 即, 对于每个 f\in X^*, 存在 u\in X 使得 Fu= f.

以下结果在计算中起着关键作用.

引理2.1 [引理2] 令 \{u_n\}\{w_n\} 为非负序列满足

u_n\le c+\sum\limits_{i=1}^{n-1} w_i u_i \ \ \mbox{ for} \ \ n \ge 1,

其中常数 c>0.

u_n\le c\exp \bigg( \sum\limits_{i=1}^{n-1}w_i \bigg) \ \ \mbox{ for} \ \ n \ge 1.

现在开始研究主要问题.

给定非线性算子 A\colon V \to V^*, B \colon V \to V^*, 和 f\in {\cal V}^*, u_0\in V,考虑如下发展问题.

问题2.1 寻找 u\in {\cal W} 使得

\begin{matrix}\left\{\begin{array}{llll}Au'(t)+Bu(t)= f(t) \qquad \mbox{ a.e.} \ \ t\in (0,T), \\u(0) = u_0.\end{array}\right.\end{matrix}

这一章的目的是证明问题2.1解的存在性.

需要如下关于问题2.1中参量的假设.

H(f):f\in C([T];V^*).

H(A):A\colon V\to V^* 弱连续, 且存在 \alpha_A,\beta_A> 0 使得

\langle Av,v\rangle_{V^*,V}\ge \alpha_A\|v\|_V^2-\beta_A\|v\|_V, \qquad \forall v \in V.

H(B)_1:B\colon V\to V^* 弱连续, 且存在 b_0, b_1\ge 0 使得

\| Bv \|_{V^*} \le b_0 + b_1\| v \|_V, \qquad \forall v \in V.

H(B)_2:B\colon V\to V^* 弱连续, 且存在 \overline{b}_0, \overline{b}_1, \overline{b}_2\ge 0\alpha_B>0 使得

\| Bv \|_{V^*} \le \overline{b}_0 + \overline{b}_1\| v \|_V+ \overline{b}_2\| v \|_V^2, \qquad \forall v \in V

\langle Bv,v\rangle_{V^*,V}\ge \alpha_B\|v\|_V^2, \qquad \forall v \in V,v\neq 0.

本文的主要结果之一如下.

定理2.2 如果 H(f), H(A) 成立, 且 H(B)_1 或者 H(B)_2 其中之一满足, 那么问题2.1有一个解 u\in {\cal W}.

以如下与问题2.1相关联的辅助问题开始定理2.2的证明.

问题2.2 寻找 w\in {\cal V} 使得

\begin{matrix} Aw(t)+B(u_0+\int_0^tw(s){ d}s)= f(t) \qquad \mbox{ a.e.} \ \ t \in (0,T). \end{matrix}

对于a.e. t \in (0, T), 令 w(t)=u'(t).

\begin{matrix} u(t)=u_0+\int_0^tw(s){ d}s, \qquad \forall t \in [T]. \end{matrix}
(2.1)

利用上面的记号, 容易观察到 w\in {\cal V} 是问题2.2的解当且仅当定义为(2.1)式的 u\in {\cal W} 是问题2.1的解. 因此, 在下文中, 为了证明定理2.2, 建立问题2.2解的存在性就够了.

下面通过应用Rothe方法来说明问题2.2的可解性.需要如下记号(见文献[25,29]).固定 N\in {\Bbb N},\tau=\frac{T}{N} 为时间步长, 对于 k=1, 2,\cdots,N, t_k=k\tau.对于任何函数 g\in {\cal V}, 定义其近似

g_\tau^k = \frac{1}{\tau}\int_{t_{k-1}}^{t_k} g(s){ d}s \qquad \ \ k = 1,\cdots,N.

此外, 对于足够小的步长 \tau, 考虑如下给出的 \int_0^tg(s){ d}s 的近似

\begin{matrix}\int_0^{t_n}g(s){ d}s=\sum\limits_{j=1}^{n}\int_{t_{j-1}}^{t_j}g(s){ d}s=\tau\sum\limits_{j=1}^{n} g_\tau^j.\end{matrix}
(2.2)

现在, 考虑如下相应于问题2.2的离散问题, 称为Rothe问题.

问题2.3 寻找一个序列 \{w_\tau^k\}_{k=1}^{N}\subset V 使得

\begin{matrix}\label{Rotheproblem}Aw_\tau^k +Bv_\tau^k= f_\tau^k\end{matrix}
(2.2)

对于 k=1,\cdots,N, 其中 v_\tau^k \in V 定义为v_\tau^k=u_0+\tau\sum\limits_{j=1}^{k}w_\tau^j.

有如下Rothe问题的存在性结果.

定理2.3 如果定理2.2的假设都满足, 那么存在 \tau_0>0 使得对于所有的 \tau\in (0,\tau_0), 问题2.3有一个解 \{w_\tau^k\}_{k=1}^{N}\subset V.

首先, 注意到问题2.3可以等价地表述如下:寻找 w_\tau^{k}\in V 使得

\begin{matrix}F(w_\tau^{k}) =f_{\tau}^k, \qquad \forall k = 1,\cdots,N,\end{matrix}
(2.3)

其中算子 F\colon V \to V^* 定义为

F(v)=Av+B \bigg(u_0+\tau\sum\limits_{j=1}^{k-1}w_\tau^j+\tau v\bigg).

接下来将通过说明 F 是满射的来证明方程(2.3)对于足够小的 \tau 是可解的. 由定理2.1, 只要说明 F 是弱连续和强制的就够了. 显然, 假设 H(A)H(B)_1H(B)_2 推出 F 弱连续.

下面将证明 F 强制. 在 H(B)_1H(B)_2 下分别证明. 首先, 假设 H(B)_1 成立. 对于所有的 v\in V, 由 H(A),H(B)_1

\begin{matrix}\langle Fv,v\rangle_{V^*,V} &\ge&\alpha_A\|v\|_V^{2}-\beta_A\|v\|_V- \big{(}b_0+b_1(K_\tau+\tau\|v\|)_V\big{)}\|v\|_V\\&=&\big{(}\alpha_A-b_1\tau\big{)}\|v\|_V^{2}- \big{(}\beta_A+b_0+b_1K_\tau\big{)}\|v\|_V,\end{matrix}

其中

\begin{matrix}K_\tau=\|u_0\|_V+\tau\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V.\end{matrix}
(2.4)

现在, 取 \tau_0=\frac{\alpha_A}{b_1}. 那么对于所有的 \tau\in (0,\tau_0), 算子 F 是强制的.

其次, 在假设 H(B)_2 下证明 F 的强制性.令z_{\tau}= u_0+\tau\sum\limits_{j=1}^{k-1}w_\tau^jv_\tau=z_{\tau}+\tau v. 那么 v= \frac{1}{\tau}(v_\tau-z_{\tau}),并且由 H(B)_2, 有

\begin{matrix}\langle Bv_\tau,v\rangle_V&=&\langle Bv_\tau,\frac{1}{\tau}(v_\tau-z_{\tau})\rangle_V=\frac{1}{\tau}\langle Bv_\tau,v_\tau\rangle_V-\frac{1}{\tau}\langle Bv_\tau,z_{\tau}\rangle_V\\&\ge&\frac{\alpha_B}{\tau}\| v_\tau \|_V^2-\frac{1}{\tau}\langle Bv_\tau,z_{\tau}\rangle_V\ge \frac{\alpha_B}{\tau}\| v_\tau \|_V^2-\frac{1}{\tau}\|z_{\tau}\|_V\|Bv_\tau\|_{V^*}\\&\ge& \frac{\alpha_B}{\tau}\| v_\tau \|_V^2-\frac{1}{\tau}\|z_{\tau}\|_V(\overline{b}_0 + \overline{b}_1\| v_\tau \|_V+ \overline{b}_2\| v_\tau \|_V^2)\\&=&\frac{\alpha_B}{\tau}\| v_\tau \|_V^2-\frac{1}{\tau}\|z_{\tau}\|_V(\overline{b}_0 + \overline{b}_1\| z_{\tau}+\tau v \|_V+ \overline{b}_2\| z_{\tau}+\tau v \|_V^2)\\&\ge&\frac{\alpha_B}{\tau}\| v_\tau \|_V^2\\&&-\frac{1}{\tau}\|z_{\tau}\|_V(\overline{b}_0 + \overline{b}_1\| z_{\tau}\|_V+\overline{b}_1\tau \|v \|_V+ \overline{b}_2\| z_{\tau}\|_V^2+2\overline{b}_2\tau \| z_{\tau}\|_V\|v \|_V+\overline{b}_2\tau^2\|v \|_V^2)\\&\ge&\frac{\alpha_B}{\tau}\| v_\tau \|_V^2 -\overline{b}_2\tau\|z_{\tau}\|_V\|v \|_V^2-\|z_{\tau}\|_V(\overline{b}_1+2\overline{b}_2\| z_{\tau}\|_V)\|v \|_V\\&&-\frac{1}{\tau}\|z_{\tau}\|_V(\overline{b}_0 + \overline{b}_1\| z_{\tau}\|_V+ \overline{b}_2\| z_{\tau}\|_V^2).\end{matrix}

从而推出

\begin{matrix}\langle Fv,v\rangle_V&\ge& \alpha_A\|v\|_V^{2}-\beta_A\|v\|_V-\overline{b}_2\tau \|z_{\tau}\|_V\|v \|_V^2-\|z_{\tau}\|_V(\overline{b}_1+2\overline{b}_2\| z_{\tau}\|_V)\|v \|_V \\&&+\frac{\alpha_B}{\tau}\| v_\tau \|_V^2-\frac{1}{\tau}\|z_{\tau}\|_V(\overline{b}_0 + \overline{b}_1\| z_{\tau}\|_V+ \overline{b}_2\| z_{\tau}\|_V^2)\\&\ge&\Big{(}\alpha_A-\overline{b}_2\tau K_\tau\Big{)}\|v\|_V^{2}- \Big{(}\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2K_\tau)\Big{)}\|v\|_V\\&&+\frac{\alpha_B}{\tau}\| v_\tau \|_V^2-\frac{1}{\tau}K_\tau(\overline{b}_0 + \overline{b}_1K_\tau+ \overline{b}_2K_\tau^2)\\&\ge&\Big{(}\alpha_A-\overline{b}_2\tau K_\tau\Big{)}\|v\|_V^{2}- \Big{(}\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2K_\tau)\Big{)}\|v\|_V-\frac{1}{\tau}K_\tau(\overline{b}_0 + \overline{b}_1K_\tau+ \overline{b}_2K_\tau^2).\end{matrix}

易证, 可以取 \tau_0 足够小使得对于所有的 \tau\in (0,\tau_0), \alpha_A-\overline{b}_2\tau K_{\tau}> 0, 从而对于所有的\tau\in (0,\tau_0), 算子 F 是强制的.

因此, 通过利用定理2.1推出对于所有的 \tau\in (0,\tau_0) 问题2.3可解. 证毕.

接下来, 建立Rothe问题中出现的序列估计.

引理2.2 如果定理2.2的假设都成立, 那么存在 \tau_0>0 使得对于所有的 \tau\in (0,\tau_0), 如下估计满足

\max_{k=1,\cdots,N}\|w_\tau^k\|_V \le C
(2.5)
max_{k=1,\cdots,N}\|v_\tau^k\|_V \le C
(2.6)

其中 C > 0 不依赖于 \tau.

先开始(2.5)式的证明.

首先, 假设 H(B)_1 成立. 对于所有的 1\le k\le N, 有

\begin{matrix}Aw_\tau^k +Bv_\tau^k= f_\tau^k.\end{matrix}

方程(2.7)两边乘上 w_\tau^k, 与定理2.3的第一部分证明类似, 得到

\langle f_\tau^k,w_\tau^k\rangle=\langle Aw_\tau^k,w_\tau^k\rangle +\langle Bv_\tau^k,w_\tau^k\rangle \ge \Big{(}\alpha_A-b_1\tau\Big{)}\|w_\tau^k\|_V^{2}- \Big{(}\beta_A+b_0+b_1K_\tau\Big{)}\|w_\tau^k\|_V.

从而得到

\|f_\tau^k\|_{V^*}+\beta_A+b_0+b_1K_\tau \ge \Big{(}\alpha_A-b_1\tau\Big{)} \|w_\tau^k\|_V,

其中 K_\tau 定义为(2.4)式. 现在, 取 \tau_0=\frac{\alpha_A}{2b_1}, 这表明对于所有的 \tau\in (0,\tau_0)\alpha_A-b_1\tau>\frac{\alpha_A}{2}. 因此, 推出

\begin{matrix}\|w_\tau^k\|&\le& \frac{2}{\alpha_A}(\|f_\tau^k\|_{V^*}+\beta_A+b_0)+\frac{2}{\alpha_A} b_1K_\tau \\&=& \frac{2}{\alpha_A}\, (\|f_\tau^k\|_{V^*}+\beta_A+b_0)+\frac{2}{\alpha_A}b_1\|u_0\|_V+ \frac{2}{\alpha_A}b_1\tau \sum\limits_{j=1}^{k-1}\|w_\tau^j\|_V.\end{matrix}

由假设 H(f), 存在一个常数 \alpha_f > 0 使得 \|f_\tau^k\|_{V^*}\le \alpha_f 对于所有的 \tau>0k\in {\Bbb N}.

\begin{matrix}m_0=\frac{2}{\alpha_A}\, (\alpha_f+\beta_A+b_0)+\frac{2}{\alpha_A}\, b_1\|u_0\|_V.\end{matrix}

由定理2.1知

\begin{matrix}\|w_\tau^k\|\le m_0 \exp\bigg{(}\frac{2}{\alpha_A}b_1\tau\, \sum\limits_{j=1}^{k-1}\int_{t_{j-1}}^{t_j}{d}s \bigg{)}\le m_0 \exp\Big{(}\frac{2}{\alpha_A}b_1\int_{0}^{t_n}{d}s\Big{)}\le m_0 \exp\Big{(}\frac{2}{\alpha_A}b_1T\Big{)}.\end{matrix}

这证明了当 H(B)_1 满足的情况下估计(2.5)成立.

然后, 假设 H(B)_2 成立. 方程(2.7)两边乘上 w_\tau^k, 与定理2.3证明的第二部分类似, 得到

\begin{matrix}\langle f_\tau^k,w_\tau^k\rangle&=&\langle Aw_\tau^k,w_\tau^k\rangle +\langle Bv_\tau^k,w_\tau^k\rangle\\&\ge& \Big{(}\alpha_A-\overline{b}_2\tau K_\tau\Big{)}\|w_\tau^k\|^{2}_V- \big{(}\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2K_\tau)\big{)}\|w_\tau^k\|_V\\&&+\frac{\alpha_B}{\tau}\|v_\tau\|_V^2-\frac{1}{\tau}K_\tau(\overline{b}_0 + \overline{b}_1K_\tau+ \overline{b}_2K_\tau^2).\end{matrix}

从而得到

\begin{matrix}&&\Big{(}\alpha_A-\overline{b}_2\tau K_\tau\Big{)}\|w_\tau^k\|^{2}_V\\&\le& \big{(}\alpha_f+\beta_A+K_\tau(\overline{b}_1 +2\overline{b}_2K_\tau)\big{)}\|w_\tau^k\|_V+ \frac{1}{\tau}K_\tau(\overline{b}_0 + \overline{b}_1K_\tau+ \overline{b}_2K_\tau^2)-\frac{\alpha_B}{\tau}\|v_\tau\|_V^2.\end{matrix}

可以取 \tau_0>0 使得 \alpha_A-\overline{b}_2\tau K_\tau>\frac{\alpha_A}{2} 对于所有的 \tau\in (0,\tau_0).\tau_0 的选取, 对于 \tau\in (0,\tau_0), 考虑如下两种情形

(i) \overline{b}_0\|u_0\|_V+\overline{b}_1\|u_0\|_{V}^2 +\overline{b}_2\|u_0\|_{V}^3\ge \alpha_B\|v_\tau\|_V^2;

(ii) \overline{b}_0\|u_0\|_V+\overline{b}_1\|u_0\|_{V}^2 +\overline{b}_2\|u_0\|_{V}^3< \alpha_B\|v_\tau\|_V^2.

如果(i)成立, 那么 \{v_\tau\} 有界. 则有

\begin{matrix}\alpha_A\|w_\tau^k\|_V^2&\le&\langle Aw_\tau^k,w_\tau^k\rangle= -\langle f_\tau^k,w_\tau^k\rangle-\langle Bv_\tau^k,w_\tau^k\rangle\\&\le& \|f_\tau^k\|_{V^*}\|w_\tau^k\|_V+\|Bv_\tau^k\|_{V^*}\|w_\tau^k\|_V\\&\le& \alpha_f\|w_\tau^k\|_V+(\overline{b}_0+\overline{b}_1 \|v_\tau^k\|_{V}+\overline{b}_2 \|v_\tau^k\|_{V}^2) \|w_\tau^k\|_V,\end{matrix}

从而

\alpha_A\|w_\tau^k\|_V\le \alpha_f+\overline{b}_0+\overline{b}_1 \|v_\tau^k\|_{V}+\overline{b}_2 \|v_\tau^k\|_{V}^2,

这表明 \{w_\tau^k\} 有界.

如果(ii)成立, 那么

\begin{matrix}\|w_\tau^k \|_V^2&\leq&\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2 K_\tau))\|w_\tau^k\|_V+\frac{2}{\alpha_A\tau}(K_\tau(\overline{b}_0+\overline{b}_1 K_\tau+\overline{b}_2 K_{\tau}^2)-\alpha_B\|v_\tau\|_V^2)\\&=&\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2 K_\tau))\|w_\tau^k\|_V+\frac{2}{\alpha_A\tau}\bigg(\Big(\|u_0\|_V+\tau\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V \Big)\\&&\times\Big(\overline{b}_0+\overline{b}_1 \Big(\|u_0\|_V+\tau\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V \Big)+\overline{b}_2 \Big(\|u_0\|_V+\tau\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V \Big)^2\Big)-\alpha_B\bigg)\\&=&\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2 K_\tau))\|w_\tau^k\|_V\\&&+\frac{2}{\alpha_A\tau}\bigg(\overline{b}_0\|u_0\|_V+\overline{b}_1\|u_0\|_{V}^2+\overline{b}_2\|u_0\|_{V}^3-\alpha_B\|v_\tau\|_V^2\\&&+\tau\Big(m_0+\tau(\overline{b}_1+3\overline{b}_2\|u_0\|_V)\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V+\overline{b}_2\tau^2\Big(\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\Big)^2\Big)\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\bigg)\\&\leq&\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2 K_\tau))\|w_\tau^k\|_V\\&&+\frac{2}{\alpha_A}\Big(m_0+\tau(\overline{b}_1+3\overline{b}_2\|u_0\|_V)\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V+\overline{b}_2\tau^2\Big(\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\Big)^2\Big)\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V,\end{matrix}

其中 m_0=\overline{b}_0+2\overline{b}_1\|u_0\|_V+3\overline{b}_2\|u_0\|_{V}^2. 利用基本事实: “如果对于a,b,x \ge 0, x^2\le ax+b, 那么 x\le a+\sqrt{b},” 得

\begin{matrix}\|w_\tau^k \|_V&\leq&\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2 K_\tau))\\&&+\bigg(\frac{2}{\alpha_A}\Big(m_0+\tau(\overline{b}_1+3\overline{b}_2\|u_0\|_V)\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V+\overline{b}_2\tau^2\Big(\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\Big)^2\Big)\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\bigg)^{\frac{1}{2}}\\&\leq&\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+K_\tau(\overline{b}_1+2\overline{b}_2 K_\tau))+\frac{1}{\alpha_A}\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\\&&+\frac{1}{2}\Big(m_0+\tau(\overline{b}_1+3\overline{b}_2\|u_0\|_V)\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V+\overline{b}_2\tau^2\Big(\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\Big)^2\Big)\\&\leq& m_1+m_2\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V+m_3\Big(\sum\limits_{j=1}^{k-1}\|w_\tau^j \|_V\Big)^2,\end{matrix}

其中

m_{1}=\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+\overline{b}_1\|u_0\|_V+2\overline{b}_2\|u_{0}\|_{V}^2)+\frac{1}{2}m_0,
m_{2}=\frac{2T}{\alpha_A}(\overline{b}_1+4\overline{b}_2\|u_{0}\|_V)+\frac{T}{2}(\overline{b}_1+3\overline{b}_2\|u_{0}\|_V), \quad m_{3}=(\frac{4}{\alpha_A}+\frac{1}{2})\overline{b}_2T^{2}.

此外, 对于 w_\tau^1, 有

\begin{matrix}\alpha_{f}\|w_\tau^1\|_V&\geq&\langle f_\tau^1,w_\tau^1\rangle_V=\langle Aw_\tau^1,w_\tau^1\rangle_V +\langle Bv_\tau^1,w_\tau^1\rangle_V\\&\geq&\alpha_A\|w_\tau^1 \|_V^2+\frac{1}{\tau}\big\langle B(u_0+\tau w_\tau^1),u_0+\tau w_\tau^1\big\rangle_V-\frac{1}{\tau}\big\langle B(u_0+\tau w_\tau^1),u_0\big\rangle_V\\&\geq&\alpha_A\|w_\tau^1 \|_V^2+\frac{1}{\tau}\alpha_B-\frac{1}{\tau}\|u_0\|_V(\overline{b}_0+\overline{b}_1\|u_0+\tau w_\tau^1\|_V+\overline{b}_2\|u_0+\tau w_\tau^1\|_{V}^2)\\&\geq&(\alpha_A-\overline{b}_2\tau \|u_0\|_V) \|w_\tau^1 \|_V^2-\|u_0\|_V(\overline{b}_1+2\overline{b}_2 \|u_0\|_V)\|w_\tau^1\|_V\\&&-\frac{1}{\tau}\|u_0\|_V(\overline{b}_0+\overline{b}_1 \|u_0\|_V+\overline{b}_2 \|u_0\|_{V}^2)+\frac{1}{\tau}\alpha_B,\end{matrix}

\begin{matrix}(\alpha_A-\overline{b}_2\tau \|u_0\|_V) \|w_\tau^1 \|_V^2&\leq&(\alpha_{f}+\beta_A+\|u_0\|_V(\overline{b}_1+2\overline{b}_2 \|u_0\|_V))\|w_\tau^1\|_V\\&&+\frac{1}{\tau}(\|u_0\|_V(\overline{b}_0+\overline{b}_1 \|u_0\|_V+\overline{b}_2 \|u_0\|_{V}^2)-\alpha_B)\\&\leq&(\alpha_{f}+\beta_A+\|u_0\|_V(\overline{b}_1+2\overline{b}_2 \|u_0\|_V))\|w_\tau^1\|_V.\end{matrix}

显然, 对于所有的 \tau\in(0,\tau_0), \alpha_A-\overline{b}_2\tau \|u_0\|_V>\frac{\alpha_A}{2}.

\begin{matrix}\|w_\tau^1 \|_V\leq\frac{2}{\alpha_A}(\alpha_{f}+\beta_A+\|u_0\|_V(\overline{b}_1+2\overline{b}_2 \|u_0\|_V)),\end{matrix}

这表明 w_\tau^1 有界. 由归纳法可以推出在假设 H(B)_2 下估计(2.5)式成立.

接着, 条件(2.6)是(2.5) 式的一个推论. 这就完成了引理的证明. 证毕.

Rothe方法的最后一步, 证明合适定义的逼近函数序列趋于问题2.2的一个解的收敛性.

考虑逐点常数插值 \overline{w}_\tau,\overline{v}_\tau\in {\cal V}\overline{f}_\tau\in {\cal V}^*,它们分别定义为

\begin{matrix}\overline{w}_\tau(t) = w_\tau^k, \qquad\overline{v}_\tau(t) =v_\tau^k, \qquad\overline{f}_\tau(t) = f_\tau^k,\end{matrix}

对于a.e. t\in (t_{k-1}, t_k], k=1,\cdots,N.

定理2.4 如果定理2.2的所有假设都满足, 那么对于一个子序列,有

\begin{matrix}\overline{w}_\tau \rightharpoonup w \ \ \mbox{在${\cal V}$里, 当} \ \tau\rightarrow 0,\end{matrix}

其中 w\in {\cal V} 是问题2.2 的一个解.

考虑分别关于 AB 的Nemytski算子 {\cal A}\colon {\cal V} \to {\cal V}^*{\cal B}\colon {\cal V} \to {\cal V}^*,即

({\cal A} v)(t) = Av(t), \ \ ({\cal B} v)(t) = Bv(t),\qquad \mbox{a.e.} \ t \in (0,T).

那么方程(2.7)变为满足

\begin{matrix}{\cal A}\overline{w}_\tau+ {\cal B}\overline{v}_\tau= \overline{f}_\tau.\end{matrix}
(2.8)

由引理2.2得到的有界性(2.5)式, 序列 \{\overline{w}_\tau\}{\cal V} 中仍有界. 利用 {\cal V} 的自反性, 如果需要, 通过取一个子列, 得到

\begin{matrix}\overline{w}_\tau \rightharpoonup w \ \ \mbox{在$ {\cal V}$里}.\end{matrix}
(2.9)

接下来, 将说明 w\in {\cal V} 是问题2.2的一个解.显然对于a.e. t\in (0,T) 序列 \{\overline{w}_\tau(t)\}V 中有界. 由于 V 自反, 那么对于 \{\overline{w}_\tau(t)\} 一个子列, 也记为 \{\overline{w}_\tau(t)\},有 \overline{w}_\tau(t)\rightharpoonup {\widetilde{w}}(t)V中 a.e. t\in (0,T).这样, 可以推出 w(t) = {\widetilde {w}}(t) 对于 a.e. t\in (0,T). 因此, 有 \overline{w}_\tau(t)\rightharpoonup w(t)V中 a.e. t\in (0,T).由于 A 弱连续, 有 A\overline{w}_\tau(t)\rightharpoonup Aw(t)V^*中对于a.e. t\in (0,T). 那么对于任何 v \in {\mathcal V} 和a.e. t\in (0,T), 有

\begin{matrix}\langle A\overline{w}_\tau(t)-Aw(t),v(t)\rangle_{V} &\ge& -\|A\overline{w}_\tau(t)-Aw(t)\|_{V^*}\|v(t)\|_{V} \\&\ge&-(\|A\overline{w}_\tau(t)\|_{V^*}+\|Aw(t)\|_{V^*})\|v(t)\|_{V}\\&\ge& -(M(\|\overline{w}_\tau(t)\|_V)+M(\|w(t)\|_V))\|v(t)\|_{V}.\end{matrix}

应用Fatou引理得到

\begin{matrix}\liminf_{\tau\rightarrow 0}\langle {\cal A}\overline{w}_\tau-{\cal A}w,v\rangle_{{\cal V}}&=& \liminf_{\tau\rightarrow 0}\int_{0}^{T}\langle A\overline{w}_\tau(t)-Aw(t),v(t)\rangle_{V}{ d}t\\&\ge&\int_{0}^{T}\liminf_{\tau\rightarrow 0}\langle A\overline{w}_\tau(t)-Aw(t),v(t)\rangle_{V} { d}t = 0.\end{matrix}

由于 v \in {\mathcal V} 任意, 推断出

\begin{matrix}\label{Aw}{\cal A}\overline{w}_\tau\rightharpoonup {\cal A}w \ \\mbox{在 {\cal V}^*中,\ 当} \ \tau \to 0.\end{matrix}
(2.10)

此外, 对于a.e. t \in (0,T), 令 (I_0v)(t)=\int_0^tv(s){d}s.由于算子 I_0\colon {\cal V}\to C([T];V) 是线性连续的, 它也是弱连续的.因此, 由(2.9)式推出

\begin{matrix}I_0\overline{w}_\tau\rightharpoonup I_0w \ \ \mbox{在 C([T];V)中, 当} \ \tau\rightarrow 0.\end{matrix}

从而也可以得到

\begin{matrix}I_0\overline{w}_\tau(t) \rightharpoonup I_0w(t) \ \ \mbox{在V中, 当} \ \tau\rightarrow 0\end{matrix}
(2.11)

对于所有的 t\in [T].实际上, 令 v^*\in V^*t\in [T].选取 \varphi\in (C([T];V))^*=M_{ca}(\sum\limits_I,V^*) 定义为 \varphi(y)=\langle v^*,y(t)\rangle_{V^*,V} 对于 y\in C([T];V).

\begin{matrix}\langle v^*,I_0\overline{w}_\tau(t)\rangle_{V^*,V}=\varphi(I_0\overline{w}_\tau)\to \varphi(I_0w)=\langle v^*,I_0w(t)\rangle_{V^*,V}, \qquad \mbox{as} \ \tau\rightarrow 0.\end{matrix}

由于 v^* \in V^* 任意, 推出(2.11)式满足.

另一方面,由有界性(2.5)式知, 存在一个常数 C>0 使得 \|\overline{w}_\tau(s)\|\le C对于a.e. s\in (0,T).那么, 对于 t\in (t_{k-1},t_k], 有

\begin{matrix}\|\overline{v}_\tau(t)-u_0-I_0\overline{w}_\tau(t)\|&=&\bigg\|\int_{0}^{t_k}\overline{w}_\tau(s){ d}s-\int_{0}^{t}\overline{w}_\tau(s){ d}s\bigg\|\nonumber\\&\le&\int_{t}^{t_k}\|\overline{w}_\tau(s)\|{ d}s\le \|\overline{w}_\tau\|_{{\cal V}}\sqrt{t_k-t}.\label{convergence2}\end{matrix}
(2.12)

由于当 \tau\rightarrow 0t\rightarrow t_k, 推出对于所有的 t\in [T], 当 \tau\rightarrow 0\overline{v}_\tau(t)-u_0-I_0\overline{w}_\tau(t) \to 0.

联立(2.11)和(2.12)式, 有

\begin{matrix}\overline{v}_\tau(t)\rightharpoonup u_0+I_0w(t)\ \ \mbox{在V中, 当} \ \tau\rightarrow 0\end{matrix}
(2.13)

对于所有的 t\in [T].类似地可以推出

\begin{matrix}{\cal B}\overline{v}_\tau\rightharpoonup {\cal B}(u_0+I_0w) \ \\mbox{在 {\cal V}^*中,\ \ 当} \ \tau \to 0.\end{matrix}
(2.14)

另外,文献[引理3.3]推出, 当 \tau \to 0\overline{f}_{\tau}\rightarrow f{\cal V}^*中. 从而由收敛性(2.9)和(2.14)式, 通过在(2.8)式中取极限得到

{\cal A}w + {\cal B}(u_0+I_0w) = f.

这表明 w\in {\cal V} 是问题2.2的一个解. 证毕.

这一章最后, 给出问题2.1的一个特殊情形.

问题2.4 寻找 u\in {\cal W} 使得

\left\{\begin{array}{llll}u'(t)+Bu(t)= f(t) \ \ \ \mbox{ a.e.} \ \ t\in (0,T), \\u(0) = u_0.\end{array}\right.

定理2.5 如果 H(f) 成立, 以及 H(B)_1 或者 H(B)_2 其中之一满足, 那么问题2.4有一个解 u\in {\cal W}.

由于对偶算子 {\cal J}:V\rightarrow V^*, 定义为 {\cal J}v=v 对于 v\in V, 是线性连续的, 它也是弱连续的. 这个结果是定理2.2的一个推论. 证毕.

3 最优状态-控制对的存在性

这一章考虑限制Lagrange最优控制问题的最优状态-控制对的存在性, 它的成本泛函受问题1.1 约束.

{\cal Y}_{ ad} {\cal Y} 中的一个非空闭凸子集. 记 u^{y}\in {\cal W} 为问题1.1关于控制 y\in {\cal Y}_{ ad} 的解. 对于每个 y \in {\cal Y}_{ ad}, 令

{\cal S}(y)=\{u\in {\cal W}| (u, y)\ \mbox {是问题1.1中系统方程的解}\}.

考虑如下Lagrange问题.

问题3.1 寻找 y^0\in {\cal Y}_{{ ad}},u^0\in {\cal S}(y^0) 使得

J(u^0,y^0)\le J(u^y,y),\qquad \forall y\in {\cal Y}_{ ad},\forall u^y\in {\cal S}(y),

其中 J(u,y):= \int _0^T{\mathfrak L}(t,u(t),y(t)){ d}t.

一对 (u,y)\in {\cal W}\times {\cal Y}_{ad} 叫做可行的如果它满足问题1.1. 若 (u,y) 可行, 则 u\in {\cal S}(y).

需要如下关于 {\mathfrak L} 的条件.

H({\mathfrak L}): 函数 {\mathfrak L}:[T]\times V\times Y\rightarrow \mathbb{R} \bigcup \{\infty \} 满足

(i) 函数 {\mathfrak L}:[T]\times V\times Y\rightarrow \mathbb{R} \bigcup \{\infty \} Borel可测;

(ii) 对于a.e. t\in [T], {\mathfrak L}(t,\cdot,\cdot )V\times Y 上是序列下半连续的;

(iii) 对于每个 x\in X 和a.e. t\in [T], {\mathfrak L}(t,x,\cdot )Y 上是凸的;

(iv) 存在一个常数 d>0 和一个函数 \psi \in L^1([T], \mathbb{R} _+) 使得

{\mathfrak L}(t,u,y)\ge \psi (t)+d\Vert y\Vert _{Y}^2.

需要如下关于 A 的条件.

H(A)_1:A\colon V\to V^*有界弱连续, 且存在 \overline{\alpha}_A,\overline{\beta}_A,\lambda_1,\lambda_2>0使得

\langle Av,v\rangle_{V^*,V}\ge \overline{\alpha}_A\|v\|_V^2-\overline{\beta}_A\|v\|_V,\qquad \forall v \in V

\int_0^t\langle Aw'(s),w(s)\rangle_{V^*,V}{d}s\ge \lambda_1\|w(t)\|_V^2-\lambda_2\|w(0)\|_V^2,\qquad \forall w \in {\cal W}.

定理3.1 如果定理2.2的条件和 H({\mathfrak L}),H(A)_1 都成立, 那么问题3.1至少有一个最优可行对.

定义 J_0(y)=\displaystyle \inf\limits_{u^y\in {\cal S}(y)}J(u^y,y), \forall y\in {\cal Y}_{ ad}.{\cal S}(y) 只有有限元素, 则存在某个 \tilde{u}^y\in {\cal S}(y) 使得 J(\tilde{u}^y,y)=\displaystyle \inf\limits_{u^y\in {\cal S}(y)}J(u^y,y)=J_0(y). 若集合 {\cal S}(y) 有无限元素且 \displaystyle \inf\limits_{u^y\in {\cal S}(y)}J(u^y,y)=+\infty, 则无需证明. 现在, 假设J_0(y)=\displaystyle \inf\limits_{u^y\in {\cal S}(y)}J(u^y,y)<+\infty .H({\mathfrak L}), 有 J_0(y)>-\infty . 现在把证明分成几步.

第1步 由下确界的定义, 存在一个序列 \{u_n^y\}\!\subseteq\! {\cal S}(y) 满足当 n\rightarrow \infty 时, J(u_n^y,y)\!\rightarrow \!J_0(y). 考虑一个可行对序列 \{u_n^y,u\}, 有

\begin{matrix}{\cal A}{u^{y}}'+{\cal B}u^{y}= f+{\cal C}y.\end{matrix}
(3.1)

第2步 现在说明存在某个 \tilde{u}^y\in {\cal S}(y) 使得 J(\tilde{u}^y,y)=\displaystyle \inf\limits_{u^y\in {\cal S}(y)}J(u^y,y)=J_0(y).

实际上, 对于 t\in[\frac{jT}{k},\frac{(j+1)T}{k}),~0\leq j\leq k-1, 有

\begin{matrix}&&\int_0^t\langle A{u_n^y}'(s),u_n^y(s)\rangle_{V^*,V}{d}s+\int_0^t\langle Bu_n^y(s),u_n^y(s)\rangle_{V^*,V}{ d}s\\&=&\int_0^t\langle f(s),u_n^y(s)\rangle_{V^*,V}{d}s+\int_0^t\langle Cy(s),u_n^y(s)\rangle_{V^*,V},\end{matrix}

那么

\begin{matrix}&&\lambda_1\|u_n^y(t)\|_V^2-\lambda_2\|u_0\|_V^2+\alpha_B\|u_n^y(t)\|_V^2\\&\le& \|f\|_{{\cal V}^*}\|u_n^y\|_{L^2(0,t;V)}+\int_0^t\langle Cy(s),u_n^y(s)\rangle_{V^*,V}\\&\le& \|f\|_{{\cal V}^*}\Big{(}\int_{0}^{t}\|u_n^y(s)\|_V^2{d}s\Big{)}^{\frac{1}{2}} +\|C\|\|y\|_{L^2(0,t;Y)}\Big{(}\int_{0}^{t}\|u_n^y(s)\|_V^2{d}s\Big{)}^{\frac{1}{2}}.\end{matrix}

易证存在 r_{0}>0 使得

\|u_n^y\|_{{\cal V}}\leq r_{0}.

类似地, 可以得到 \{(u_n^y)'\}{\cal V} 中有界. 因此, \{u_n^y\}{\cal W} 中有界.

不失一般性, 可以假设当 n\rightarrow \infty u_n^y\rightharpoonup \tilde{u}^y{\cal W}中对于 y\in {\cal Y}_{ad}. 在(3.1)式两边令 n\rightarrow \infty, 由Lebesgue控制收敛定理得到

{\cal A}(\tilde{u}^y)'+{\cal B}\tilde{u}^y= f+{\cal C}y,

这表明 \tilde{u}^y\in {\cal S}(y). 现在验证对于任何 y\in {\cal Y}_{ ad}, J(\tilde{u}^y,y)=\displaystyle \inf\limits_{u^y\in {\cal S}(y)}J(u^y,y)=J_0(y). 实际上, 由于 {\cal W} 连续嵌入到 L^1(0,T;V) 中, 通过可行对的定义, 假设 H({\mathfrak L}) 和Balder定理[30], 有

J_0(y)=\displaystyle \lim _{n\rightarrow \infty }\int _0^T{\mathfrak L}(t,u_n^y(t),y(t)){d}t\ge \displaystyle \int _0^T{\mathfrak L}(t,\tilde{u}^y(t),y(t)){d}t =J(\tilde{u}^y,y)\ge J_0(y),

J(\tilde{u}^y,y)=J_0(y). 这表明对每个 y\in {\cal Y}_{ad}, J_0(y)\tilde{u}^y\in {\cal W} 达到最小.

第3步 证明存在 y^0\in {\cal Y}_{ad} 使得对于所有的 u\in {\cal Y}_{ ad}, J(y^0)\!\le\! J_0(y).\inf\limits_{u\in Y_{{ ad}}}J_0(y)\!=\!+\infty , 则无需证明. 假设 \inf\limits_{u\in {\cal Y}_{ ad}}J_0(y)<+\infty . 与第1步类似, 可以证明 \inf\limits_{y\in {\cal Y}_{ ad}}J_0(y)>-\infty , 且存在一个序列 \{y_n\}\subseteq {\cal Y}_{ ad} 使得当 n\rightarrow \infty 时, J(u_n)\rightarrow \inf\limits_{u\in {\cal Y}_{ ad}}J_0(y).H({\mathfrak L})(iv) 得到集合 \{y_n\}{\cal Y} 中有界, 由于 {\cal Y} 是一个自反Banach空间, 存在一个子列, 还记为 \{y_n\}, 当 n\rightarrow \infty 时, 弱收敛于某个 y^0\in {\cal Y}. 注意到 {\cal Y}_{ ad} 是闭凸的, 得到 y^0\in {\cal Y}_{ ad}.

假设 \tilde{u}^{y_n} 是问题1.1关于 y_n 的解, 其中 J(y_n) 达到最小值. 那么, (\tilde{u}^{y_n},y_n) 是一个可行对且能验证对于 t\in (0,T),有

\begin{matrix}{\cal A}(\tilde{u}^{y_n})'+{\cal B}\tilde{u}^{y_n}= f+{\cal C}y_n.\end{matrix}
(3.2)

类似可以得到集合 \{\tilde{u}^{y_n}\}\subset {\cal W} 有界, 且存在子列, 还记为 \{\tilde{u}^{y_n}\}, 和 \tilde{u}^{y^0}\in {\cal W}, 使得当 n\rightarrow \infty 时, \tilde{u}^{y_n}\rightharpoonup \tilde{u}^{y^0}{\cal W}里. 若(3.2) 式两边 n\rightarrow \infty , 则

{\cal A}(\tilde{u}^{y^0})'+{\cal B}\tilde{u}^{y^0}= f+{\cal C}y^0,

这表明 (\tilde{u}^{y^0},y^0) 是一个可行对. 由于 {\cal W} 连续嵌入到 L^1(0,T;V) 中, 由假设H({\mathfrak L}) 和Balder定理[30], 有

\begin{matrix}\inf\limits_{y\in {\cal Y}_{ ad}}J_0(y)&=& \lim _{n\rightarrow \infty }\int _0^T{\mathfrak L}(t,\tilde{u}^{y_n}(t),y_n(t)){ d}t \ge \int _0^T{\mathfrak L}(t,\tilde{u}^{y^0}(t),y^0(t)){ d}t\\&=& J(\tilde{u}^{y^0},y^0)\ge \inf\limits_{y\in {\cal Y}_{ ad}}J_0(y).\end{matrix}

因此, J(\tilde{u}^{y^0},y^0)=J(y^0)=\displaystyle \inf\limits_{x^{y^0}\in {\cal S}{(y^0)}}J(x^{y^0},y^0). 此外, J(y^0)=\displaystyle \inf\limits_{y\in {\cal Y}_{ ad}}J_0(y),Jy^0\in {\cal Y}_{ ad} 达到最小值. 证毕.

定理3.2 假设定理2.2的假设和 H(A) 都成立. 若 H({\mathfrak L})(iv) 被替换为

(iv)' 存在常数 c,d>0 和一个函数 \psi \in L^1([T], \mathbb{R} _+) 使得

{\mathfrak L}(t,u,y)\ge \psi (t)+c\Vert u\Vert _{V}^2+d\Vert y\Vert _{Y}^2,

那么问题3.1至少有一个最优可行对.

与定理3.1的证明类似, 由 (iv)' 得到 \{u_k\}{\cal V} 中有界. 其余证明略.

4 不稳定的Navier-Stokes-Voigt方程

这一章把定理2.2中主要结果应用到证明Navier-Stokes-Voigt方程解的存在性.

过去一些年, Navier-Stokes-Voigt方程解的存在性吸引了许多数学家的注意(参见文献[32-37]). 注意到在很多实际情况中, Navier-Stokes-Voigt方程中的一阶时间导数可以经常被忽略不计, 之后一些时间流体流动变为几乎静止. 因此, 一些应用只研究稳态解, 另一方面, 这可能是不够的和不准确的.用这个办法, 得到一个发展方程, 它是一个Kelvin-Voigt粘弹性不可压缩流体动力学的弱形式.该方程整体解的存在性定理在许多方面得到了推广, 在文献[32]中得到了Navier-Stokes-Voigt 方程的一些结果.

考虑如下Navier-Stokes-Voigt方程.

\begin{matrix}\label{NSV}\left\{\begin{array}{llll}u'_t-\nu\Delta u-\alpha^2\Delta u'+(u\cdot \nabla)u+\nabla p=f+y, \qquad& \mbox{在$ Q$里}, \\\nabla\cdot u=0, \qquad &\mbox{在$ Q$里}, \\u(x,t)=0,\qquad &\mbox{在$ \partial\Omega\times (0,T)$上}, \\u(x,0)=u_0(x),\qquad& \mbox{在$ \Omega$里}.\end{array}\right.\end{matrix}
(4.1)

下面, 首先给出方程(4.1)的变分形式, 它是一个发展方程, 然后说明怎样利用定理2.2证明发展方程解的存在性.

定义如下两个内积

(u,v)_1:=\int_{\Omega}\sum\limits_{j=1}^3u_jv_j{ d}s,\qquad u=(u_1,u_2,u_3),v=(v_1,v_2,v_3)\in (L^2(\Omega))^3

(u,v)_2:=\int_{\Omega}\sum\limits_{j=1}^3\nabla u_j\cdot\nabla v_j{ d}s,\qquad u=(u_1,u_2,u_3),v=(v_1,v_2,v_3)\in (H_0^1(\Omega))^3,

及其对应的范数\|v\|_1:=\sqrt{(v,v)_1},\|v\|_2:=\sqrt{(v,v)_2}.

\begin{matrix}V_0=\{v\in (C_0^{\infty}(\Omega))^3:\nabla \cdot v=0\}.\end{matrix}

HV_0(L^2(\Omega))^3 中的闭包, 记 VV_0(H_0^1(\Omega))^3 中的闭包. 显然,HV 是Hilbert空间, 内积分别为 (\cdot,\cdot)_1(\cdot,\cdot)_2. 因此 V\subset H \subset V^*, 其中的嵌入是稠密、连续、紧的. 那么存在一个常数 \overline{c}>0 使得 \|v\|_H\le \overline{c}\|v\|_V. 显然空间 V 是一个自反Banach空间, 范数为\| v \| = \| v \|_{H^{1}(\Omega;\mathbb{R} ^d)},v \in V.

{\cal V}=L^2(0,T;V), \ \ {\cal V}^*=L^2(0,T;V^*), \ \ {\cal H}=L^2(0,T;H), \ \ {\cal W}=W^{1,2}(0,T;V).

众所周知, 嵌入 {\cal W}\hookrightarrow C([T];V) 是连续, 且嵌入 {\cal W}\hookrightarrow (L^2(Q))^3 是紧的(参见文献[38]). {\cal V}^*{\cal V} 之间的对偶积为

\begin{matrix}\langle v^*,v\rangle_{{\cal V}^*,{\cal V}}:=\int_0^T\langle v^*(t),v(t)\rangle_{V^*,V}{ d}t,\qquad v\in {\cal V},v^*\in {\cal V}^*.\end{matrix}

现在定义三线性形式 b:V\times V\times V\rightarrow \mathbb{R}

\begin{matrix}b(u,v,w):=\sum\limits_{i,j}^3\int_{\Omega} u_i\frac{\partial v_j}{x_i}w_j{ d}s,\end{matrix}

其中的积分有意义. b 的性质可以参见文献[39]第II章引理1.1和引理1.3或者文献[40,41].

现在给出方程(4.1)弱解的定义.

定义4.1[31,32] 对于 f\in C([T];V^*),u_0\in V, u\in {\cal W} 叫做方程(4.1)在 (0,T) 的一个弱解, 如果

\left\{\begin{array}{llll}( u'(t),v)_1+\nu(u(t),v)_2+\alpha^2(u'(t),v)_2+b(u(t),u(t),v)\\=\langle f(t)+y(t),v\rangle_{V^*,V},\ \forall v\in V, \ \mbox{a.e.} \ t\in (0,T),\\u(0)=u_0.\end{array}\right.

下面想要给出一个方程的等价形式.为此, 介绍一个线性连续算子 E:V\rightarrow V^*

\langle Eu,v\rangle_{V}:=(u,v)_2,\qquad \forall u,v\in V,

和一个非线性算子 F:V\rightarrow V^*

\langle Fu,v\rangle_{V}:=b(u,u,v),\qquad \forall u,v\in V.

现在, 有定义4.1的等价形式.

定义4.2 一个函数 u\in {\cal W} 叫做方程(4.1)在区间 (0,T) 上的一个弱解, 如果

\begin{matrix}\left\{\begin{array}{llll}u'+\nu Eu+\alpha^2Eu'+Fu=f+y,\qquad &\mbox{在$ {\cal V}^*$里},\\u(0)=u_0,\qquad &\mbox{在$ V$里}.\end{array}\right.\end{matrix}

如下方程(4.1)的存在性结果是定理2.2的一个推论.

定理4.1 方程(4.1)有一个弱解.

用如下泛函结构来应用定理2.2:

\begin{matrix} A\colon V\to V^*,~~&&A={\cal J}+\alpha^2E,\\B \colon V\to V^*,~~&&B=\nu E+F,\\C\colon (L^{2}(\Omega))^3\to V^*,~~&&Cy=y, \qquad \forall y\in (L^{2}(\Omega))^3.\end{matrix}

接下来验证 H(f),H(A),H(B)_2.

H(f)显然. 易知 J 是有界弱连续的. E,F 的有界性和弱连续性的验证可以参见文献[31,41], 从而 A,B 是有界弱连续的. 此外, 对于 v\in V,w \in {\cal W}, 有

\langle Av,v\rangle_{V^*,V}=\langle {\cal J}v,v\rangle_{V^*,V}+\langle \alpha^2Ev,v\rangle_{V^*,V}=(v,v)_1+\alpha^2(v,v)_2\ge \alpha^2\|v\|_V^2,
\langle Bv,v\rangle_{V^*,V}=\langle \nu Ev,v\rangle_{V^*,V}+\langle Fv,v\rangle_{V^*,V}=\nu(v,v)_2= \nu\|v\|_V^2,

这表明 H(A),H(B)_2 成立且\alpha_A=\alpha^2,\beta_A=0,\alpha_B=\nu. 证毕.

考虑如下二次目标泛函的最小化

J(u,y)=\frac{\alpha_Q}{2}\int_0^T\int_Q|u(x,t)-u_Q(x,t)|^2{ d}x{ d}t +\frac{\gamma}{2}\int_0^T\int_Q|y(x,t)|^2{ d}x{ d}t.

其中 u_Q\in (L^2(Q))^3, 系数 \alpha_Q 是一个正常数, 正则参数 \gamma 用来调节控制的损耗, 也是一个正常数.

容许控制集记为 {\cal Y}_{ ad}, 在 (L^2(Q))^3 中是非空凸闭的. 考虑如下最优控制问题:

问题4.1 最小化上面的泛函 J(u,y), 满足方程(4.1)和限制 y\in {\cal Y}_{ ad}.

最后有如下结果, 它是定理3.1的一个推论.

定理4.2 问题4.1最少有一个最优可行对.

{\mathfrak L}(t,u(t),y(t))=\frac{\alpha_Q}{2}\int_Q|u(x,t)-u_Q(x,t)|^2{ d}x+\frac{\gamma}{2}\int_Q|y(x,t)|^2{ d}x.

易知 H({\mathfrak L}) 成立. 现在验证 H(A)_1.

对于 w \in {\cal W}, 有

\begin{matrix}\int_0^t\langle Aw'(s),w(s)\rangle_{V^*,V}{ d}s&=&\int_0^t\langle {\cal J}w'(s),w(s)\rangle_{V^*,V}{ d}s+\int_0^t\langle \alpha^2Ew'(s),w(s)\rangle_{V^*,V}{ d}s\\&=&\int_0^t(w'(s),w(s))_1{ d}s+\alpha^2\int_0^t(w'(s),w(s))_2{ d}s\\&=& \frac{1}{2}(\|w(t)\|_H^2-\|w(0)\|_H^2)+\frac{\alpha^2}{2}(\|w(t)\|_V^2-\|w(0)\|_V^2)\\&\ge& -\frac{\overline{c}^2}{2}\|w(0)\|_V^2+\frac{\alpha^2}{2}(\|w(t)\|_V^2-\|w(0)\|_V^2),\end{matrix}

这表明 H(A)_1 成立, 且\overline{\alpha}_A=\alpha^2,\overline{\beta}_A=0,\lambda_1=\frac{\alpha^2}{2},\lambda_2=\frac{\overline{c}^2}{2}+\frac{\alpha^2}{2}.然后, 通过应用定理3.1, 问题4.1最少有一个最优可行对.证毕.

参考文献

Ahmed N U.

Existence of optimal controls for a class of systems governed by differential inclusions in Banach spaces

J Optim Theory Appl, 1986, 50: 213-237

DOI:10.1007/BF00939270      URL     [本文引用: 1]

Ahmed N U.

Existence of optimal relaxed controls for differential inclusions on Banach space//Lakshmi- kantham

Nonlinear Analysis and Applications. New York: CRC Press, 1987: 39-49

[本文引用: 1]

Ahmed N U.

Properties of the relaxed trajectories for a class of nonlinear evolution matrixs on a Banach space

SIAM Journal on Control and Optimization, 1983, t2: 953-968

[本文引用: 1]

Papageorgiou N S.

Existence of optimal controls for nonlinear systems in Banach spaces

J Optim Theory Appl, 1987, 53(3): 1581-1600

[本文引用: 1]

Papageorgiou N S.

Properties of the relaxed trajectories of evolution matrixs and optimal control

SIAM J Contr Optim, 1989, 27: 267-288

DOI:10.1137/0327014      URL     [本文引用: 1]

Papageorgiou N S.

Relaxation and existenee of optimal controls for systems governed by evolution inclusions in separable Banach spaces

J Optim Theory Appl, 1990, 64: 573-594

DOI:10.1007/BF00939425      URL     [本文引用: 1]

Cesari L.

Existence of Solutions and Existence of Optimal Solutions

Mathematical Theories of Optimization, 1983, 979: 88-107

[本文引用: 1]

Cesari L. Optimization:Theory and Applications. New York: Springer-Verlag, 1983

[本文引用: 1]

Ahmed N U, Teo K L. Optimal Control of Distributed-Parameter Systems. New York: North-Holland, 1981

[本文引用: 1]

Lions J L. Optimal Control of Systems Governed by Partial Differential Equations. New York: Springer-Verlag, 1971

[本文引用: 1]

Zeng B.

Feedback control systems governed by evolution matrixs

Optimization, 2019, 68: 1223-1243

DOI:10.1080/02331934.2019.1578358      URL     [本文引用: 1]

Zeng B, Liu Z H.

Existence results for impulsive feedback control systems

Nonlinear Analysis: Hybrid Systems, 2019, 33: 1-16

DOI:10.1016/j.nahs.2019.01.008      URL     [本文引用: 1]

Aizicovici S, Pavel N H.

Anti-periodic solutions to a class of nonlinear differential matrixs in Hilbert space

J Functional Analysis, 1991, 99: 387-408

DOI:10.1016/0022-1236(91)90046-8      URL     [本文引用: 1]

Barbu V.

Nonlinear Semigroup and Differential Equations in Banach Spaces

Leyden: Noordhoff, 1976

[本文引用: 1]

Carstensen C, Gwinner J.

A theory of discretization for nonlinear evolution inequalities applied to parabolic Signorini problems

Ann Mat Pura Appl, 1999, 177: 363-394

DOI:10.1007/BF02505918      URL     [本文引用: 1]

Denkowski Z, Migórski S, Papageorgiou N S. An Introduction to Nonlinear Analysis:Theory. New York: Kluwer Academic/Plenum Publishers, 2003

[本文引用: 1]

Li X J, Yong J M.

Optimal Control Theory for Infinite Dimensional Systems

Boster: Birkhäuser, 1995

[本文引用: 1]

Guo B Z, Xu Y S, Yang D H.

Optimal actuator location of minimum norm controls for heat matrix with general controlled domain

J Differential Equations, 2016, 261: 3588-3614

DOI:10.1016/j.jde.2016.05.037      URL     [本文引用: 1]

Reyes J C D L, Dhamo V.

Error estimates for optimal control problems of a class of quasilinear matrixs arising in variable viscosity fluid flow

Numer Math, 2016, 132: 691-720

DOI:10.1007/s00211-015-0737-2      URL     [本文引用: 1]

Sattayatham P.

Strongly nonlinear impulsive evolution matrixs and optimal control

Nonlinear Analysis, 2004, 57: 1005-1020

DOI:10.1016/j.na.2004.03.025      URL     [本文引用: 1]

Tanabe H.

Equations of Evolution

London: Pitman, 1979

[本文引用: 1]

Xiang X L.

Optimal controls for a class of strongly nonlinear evolution matrixs with constraint

Nonlinear Analysis, 2001, 47: 57-66

DOI:10.1016/S0362-546X(01)00156-0      URL     [本文引用: 1]

Yamazaki N.

Convergence and optimal control problems of nonlinear evolution matrixs governed by time-dependent operator

Nonlinear Analysis, 2009, 70: 4316-4331

DOI:10.1016/j.na.2008.09.015      URL     [本文引用: 1]

Francŭ J.

Weakly continuous operators, applications to differential matrixs

Application of Mathematics, 1984, 39(1): 45-56

DOI:10.21136/AM      URL     [本文引用: 3]

Roubiček T.

Nonlinear Partial Differential Equations with Applications

Berlin: Birkhäuser, 2005

[本文引用: 3]

Zeng B, Migórski S.

Evolutionary subgradient inclusions with nonlinear weakly continuous operators and applications

Comput. Math. Appl., 2018, 75: 89-104

[本文引用: 1]

Zeng B.

Feedback control for nonlinear evolutionary matrixs with applications

Nonlinear Analysis: Real World Applications, 2022, 66: 103535.

[本文引用: 1]

Shen S, Liu F, Chen J, et al.

Numerical techniques for the variable order time fractional diffusion matrix

Appl Math Comput, 2012, 218: 10861-10870

DOI:10.1016/j.amc.2012.04.047      URL    

Kačur J.

Method of Rothe in Evolution Equations

Teubner-Texte zur Mathematik 80. Leipzig: B G Teubner, 1985

[本文引用: 1]

Balder E J.

Necessary and sufficient conditions for L^1-strong weak lower semicontinuity of integral functionals

Nonlinear Anal, 1987, 11: 1399-1404

DOI:10.1016/0362-546X(87)90092-7      URL     [本文引用: 2]

Anh C T, Nguyet T M.

Optimal control of the instationary three dimensional Navier-Stokes-Voigt matrixs

Numerical Functional Analysis and Optimization, 2016, 37(4): 415-439

DOI:10.1080/01630563.2015.1136891      URL     [本文引用: 2]

Anh C T, Trang P T.

Pull-back attractors for three-dimensional Navier-Stokes-Voigt matrixs in some unbounded domains

Proc Royal Soc Edinburgh Sect A, 2013, 143: 223-251

DOI:10.1017/S0308210511001491      URL     [本文引用: 3]

We study the first initial–boundary-value problem for the three-dimensional non-autonomous Navier–Stokes–Voigt equations in an arbitrary (bounded or unbounded) domain satisfying the Poincaré inequality. The existence of a weak solution to the problem is proved by using the Faedo–Galerkin method. We then show the existence of a unique minimal finite-dimensional pull-back \\smash{\\mathcal D_\\sigma}-attractor for the process associated with the problem, with respect to a large class of non-autonomous forcing terms. We also discuss relationships between the pull-back attractor, the uniform attractor and the global attractor.

Celebi A O, Kalantarov V K, Polat M.

Global attractors for 2D Navier-Stokes-Voight matrixs in an unbounded domain

Appl Anal, 2009, 88: 381-392

DOI:10.1080/00036810902766682      URL     [本文引用: 1]

García-Luengo J, Marín-Rubio P, Real J.

Pullback attractors for three-dimensional non-autonomous Navier-Stokes-Voigt matrixs

Nonlinearity, 2012, 25: 905-930

DOI:10.1088/0951-7715/25/4/905      URL     [本文引用: 1]

Kalantarov V K, Titi E S.

Global attractor and determining modes for the 3D Navier-Stokes-Voight matrixs

Chin Ann Math Ser B, 2009, 30: 697-714

DOI:10.1007/s11401-009-0205-3      URL     [本文引用: 1]

Kalantarov V K, Titi E S.

Gevrey regularity for the attractor of the 3D Navier-Stokes-Voight matrixs

J Nonlinear Sci, 2009, 19: 133-152

DOI:10.1007/s00332-008-9029-7      URL     [本文引用: 1]

Yue G, Zhong C K.

Attractors for autonomous and nonautonomous 3D Navier-Stokes-Voight matrixs

Discrete Cont Dyna Syst Ser B, 2011, 16: 985-1002

[本文引用: 1]

Robinson J C. Infinite-Dimensional Dynamical Systems. Cambridge: Cambridge University Press, 2001

[本文引用: 1]

Temam R.

Navier-Stokes Equations: Theory and Numerical Analysis, 2nd edition

North-Holland: Amsterdam, 1979

[本文引用: 1]

Contantin P, Foias C. Navier-Stokes Equations, Chicago Lectures in Mathematics. Chicago: University of Chicago Press, 1988

[本文引用: 1]

Dudek S, Kalita P, Migórski S.

Stationary flow of non-Newtonian fluid with nonmonotone frictional boundary conditions

Z Angew Math Phys, 2015, 66: 2625-2646

DOI:10.1007/s00033-015-0545-7      URL     [本文引用: 2]

/