常矩阵耦合的热方程采样控制系统能控性
Controllability of the Sampled-Data Controlled System of Heat Equations Coupled by Constant Matrices
DOI: 10.12677/aam.2024.134161, PDF, HTML, XML, 下载: 46  浏览: 76  科研立项经费支持
作者: 杨济通, 刘汉兵:中国地质大学数学与物理学院,湖北 武汉
关键词: 采样控制零能控近似能控采样周期Sampled-Data Controlled Null Controllability Approximate Controllability Sampled-Data Period
摘要: 本文研究由常矩阵耦合的热方程采样控制系统的零能控和近似能控问题,所考虑的是给定采样周期情况下系统的两种能控性。我们通过热方程的唯一延拓性,Kalman能控性秩条件和采样周期的选取对系统的能控性进行了刻画。提出在给定采样周期下此类系统近似能控的一个充要条件;并说明当控制施加在整体区域时,这一条件也是在给定采样周期下此类系统零能控的一个充要条件,而当控制只施加在局部区域时,此类系统对任意周期都不是零能控的。
Abstract: This paper studies the null controllability and approximate controllability for sampled-data controlled systems of heat equations coupled by constant matrices. We consider two types of controllability under a given sampled-data period and characterize the controllability of the system through the unique continuation property of the heat equation, the Kalman rank condition and the selection of the sampled-data period. We provide a necessary and sufficient condition for the approximate controllability of such systems under a given sampled-data period; and show that when the control domain is equal to the entire space, this condition is also a necessary and sufficient condition for the null controllability of such systems under a given sampled-data period. However, when the control domain is a proper subset of the entire space, such systems are not null controllable for any period.
文章引用:杨济通, 刘汉兵. 常矩阵耦合的热方程采样控制系统能控性[J]. 应用数学进展, 2024, 13(4): 1696-1708. https://doi.org/10.12677/aam.2024.134161

1. 引言

本文的目标是研究关于一类常矩阵耦合的热方程采样控制系统的能控性,并分别提出此类系统的 T 零能控与 T 近似能控的充要条件,以及研究条件的适用范围。

在现有的研究文献中,连续时间控制系统及其时间离散化系统的能控性问题已取得诸多理论成果。如郭戈等 [1] 针对周期性系统(包括连续和离散两种情形)的能控性问题。对已有文献中有关周期为T的连续周期性系统的能控性结论进行了完善,推导出系统能控等价于单周期内能控的结论,并将这些结论直接推广到离散周期性系统;进一步探讨了离散前后周期性系统能控性保持不变的可能性 [2] ,证明在采用某种不等间距采样方式时,离散前后周期性系统的能控性保持不变。杨磊等 [3] 用循环不变子空间的性质研究了线性切换系统的能控性和能观性。得到线性周期切换系统完全能控和完全能观测的充分必要条件,进一步给出了一般线性切换系统完全能控和完全能观测的充分条件和必要条件。王俭等 [4] 就离散时间系统保持能控性与能观性的充分条件的物理意义,论证了该充分条件实质上就是信号可重构的临界条件,同时说明该充分条件实际上也是必要条件,在两个概念之间建立起类似于等价的联系。随着数字控制器以及计算机的发展,因为大多数控制都是通过数字技术在系统中实现的,这激发了对具有非连续控制的连续时间系统的研究。其中控制只在有限个时刻取到不为零的值的被称为脉冲控制,而控制在作用时间内取分段常值且只进行有限次变化的被称为采样控制 [5] 。有关热方程脉冲控制的能控性问题,Qin S.等 [6] 对Liu X.等 [7] 关于实现常微脉冲控制系统能控性的最小控制个数问题进行了完善,并得到了具有脉冲控制的热方程连续时间系统近似能控的一个充要条件。Wang L.等 [8] 研究了含实矩阵的热方程的约束近似零能控性,分别给出控制作用于系统全局和局部两种情况下系统是全局约束近似零能控的条件。而对于热方程采样控制系统的能控性研究,Wang G.等 [9] 对有界域中热方程的采样周期最优时间和最优控制关于采样间隔做了上下界的误差分析,给出了一种具有一致代价的 L 2 -近似零能控性。在此基础上设计了具有采样数据控制的热方程的时间最优控制,用来近似具有分布式控制的热方程的时间最优控制问题并证明其某种意义上的最优性。Lin P.等 [10] 中给出了一类具有时不变控制热方程的能观不等式和有一致代价的部分近似零能控性,并由此进一步建立了使有界域中采样控制的热方程稳定的输出反馈律。

与以往采样控制系统能控性有关的研究相比,本文所考虑的是关于热方程组的能控性问题,我们借鉴了Qin S.等 [6] 关于系统能控性的一些想法。在本文中我们对由一对常矩阵 ( A , B ) 耦合的热方程采样控制系统进行了研究,定义了两种与采样周期T有关的能控性,结合热方程的唯一延拓性得出了关于系统能控性的两个主要结论。

2. 问题描述与主要结论

2.1. 问题介绍

为介绍本文研究的采样控制系统,这里我们先给出一些必要的定义和符号说明。设 Ω N ( N + = { 1 , 2 , } ) 是一具有 C 2 边界 Ω 的有界域, ω Ω 是一非空开子集,其特征函数为 χ ω ω ¯ 表示 ω N 中的闭包,A和B分别是 n × n n × m 实矩阵,视为从 n m n 的线性算子, Δ n = d i a g { Δ , , Δ } Δ 为拉普拉斯算子。

定义 A = Δ n A D ( A ) = H 2 ( Ω ; n ) H 0 1 ( Ω ; n ) ,即对任意 z D ( A ) A z = Δ n z A z ,可以得出 A L 2 ( Ω ; n ) 上生成了一个 C 0 半群 { e A t } t 0

本文考虑如下热方程采样控制系统的零能控与近似能控问题:

{ t y ( t ) A y ( t ) = χ ω B k = 1 p χ [ ( k 1 ) T , k T ) ( t ) u k , Ω × + , y ( 0 ) = y 0 . (2.1)

其中 χ [ ( k 1 ) T , k T ) ( t ) 表示区间 [ ( k 1 ) T , k T ) 上的特征函数, { k T } k = 1 p ( p + ) 被称为采样时刻, T > 0 ( T ) 被称为采样周期, { u k } k = 1 p L 2 ( Ω ; m ) 被称为采样控制,容易得出方程(2.1)是适定的。记 y ( ; y 0 , { u k } k = 1 p ) 为方程(2.1)的解,即

y ( t ; y 0 , { u k } k = 1 p ) = e A t y 0 + k = 1 p ( k 1 ) T k T e A ( t τ ) d τ χ ω B u k , t 0 . (2.2)

在本文中记 , 分别表示 L 2 ( Ω ; n ) 上的范数和内积, σ ( A ) σ ( Δ n ) 分别表示所有A和 Δ n 的特征值的集合, A B A 分别代表A,B和 A 的伴随算子。

本文所考虑关于系统(2.1)的零能控与近似能控定义如下:

定义2.1

给定 T > 0 ,若对任意 y 0 L 2 ( Ω ; n ) ,存在 p + ,和 { u k } k = 1 p L 2 ( Ω ; m ) 使得

y ( p T ; y 0 , { u k } k = 1 p ) = 0. (2.3)

则称系统(2.1)是T-零能控的。

定义2.2

给定 T > 0 ,若对任意 ε > 0 y 0 , y 1 L 2 ( Ω ; n ) ,存在 p + { u k } k = 1 p L 2 ( Ω ; m ) 使得

y ( p T ; y 0 , { u k } k = 1 p ) y 1 ε . (2.4)

则称系统(2.1)是T-近似能控的。

定义2.3

对给定 n × n 实矩阵A,假设有q ( q + ) 个互不相同的特征值 β 1 , β 2 , , β q ,若T对满足

Re [ β i β j ] = 0 ,对任意 i , j = 1 , 2 , , q

的A的特征值成立

T 2 l π Im ( β i β j ) , l = ± 1 , ± 2 , (2.5)

则称采样周期T是非病态的,反之称采样周期T是病态的。

定义2.4

对给定 n × n 实矩阵A,定义如下集合 T 1 T 2

T 1 : = { T | 0 T e ( λ i β j ) τ d τ = 0 , i + , j { 1 , , q } , λ i σ ( Δ n ) , β j σ ( A ) } . (2.6)

T 2 : = { T | 0 T e ( λ i β j ) τ d τ 0 , i + , j { 1 , , q } , λ i σ ( Δ n ) , β j σ ( A ) } . (2.7)

易知若采样周期T是非病态的,则 T T 2

2.2. 主要结论

本文的主要结论阐述为以下两个定理。

定理2.1

(i) 当 Ω \ ω ¯ 时,对于任意采样周期T,系统(2.1)不是T-零能控的。

(ii) ω = Ω 时,系统(2.1)是T-零能控的当且仅当 T T 2 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n

定理2.2

系统(2.1)是T-近似能控的当且仅当 T T 2 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n

推论2.3

(i) ω = Ω 时,若T非病态且 r a n k ( B , A B , , A n 1 B ) = n ,则系统(2.1)是T-零能控的。

(ii) 若T非病态且 r a n k ( B , A B , , A n 1 B ) = n ,则系统(2.1)是T-近似能控的。

3. 定理2.1的证明

下面我们先给出几个证明定理2.1所需要的引理,引理3.1和引理3.2的证明过程分别参见文献 [6] 的定理3.2和引理4.3。

引理3.1设 ω 1 Ω 的非空开子集,V是 n 的子集,则以下两条结论成立:

(i) 若 z L 2 ( Ω ; n ) t > 0 z ( x ) V 对几乎处处的 x Ω e Δ n t z ( x ) V 对于所有 x ω 1

(ii) 若 z L 2 ( Ω ; n ) t > 0 ,则 χ ω 1 e A t z = 0 当且仅当 z = 0

引理3.2 设 ( A , B ) 满足 r a n k ( B , A B , , A n 1 B ) < n 。则存在 z L 2 ( Ω ; n ) \ { 0 } 使得 χ ω B e A t z = 0 对于任意 t > 0

引理3.3 对任意 t 0 有以下等式成立

e A t = e A t e Δ n t .

因为引理3.3的证明要点是 Δ n A 的交换性,方便起见将这部分放在了引理3.5的证明中,这里省略了具体证明。

引理3.4对任意 t > 0 T T 2 则有以下两结论成立:

(i) 对于任意的 y ¯ L 2 ( Ω ; n ) ,存在唯一 y ^ L 2 ( Ω ; n ) 使得

0 T e A τ d τ y ^ = e A t y ¯ .

(ii) 对于任意的 y ¯ L 2 ( Ω ; n ) ,存在唯一 y ^ L 2 ( Ω ; n ) 使得

0 T e A τ d τ y ^ = e A t y ¯ .

引理3.4的证明:(i) 首先我们先证明对于任意的 y ¯ 1 L 2 ( Ω ; n ) T T 2 t > 0 ,存在唯一 y ^ L 2 ( Ω ; n ) 使得

0 T e Δ n A τ d τ y ^ = e Δ n t y ¯ 1 (3.1)

Ψ k j = φ k e j e j n 的单位向量,则 { Ψ i j } i 1 , 1 j n L 2 ( Ω ; n ) 的一组正交基。对于任意的 y ¯ 1 L 2 ( Ω ; n ) 可将其表示为

y ¯ 1 = k = 1 + j = 1 n b k j Ψ k j , b k j .

并且有 y ¯ 1 L 2 ( Ω ; n ) = ( k = 1 + j = 1 n b k j 2 ) 1 / 2 。显然 y ^ L 2 ( Ω ; n ) 是(3.1)的解等价于 y ^ 可唯一表示为 y ^ = k = 1 + j = 1 n a k j Ψ k j , a k j ,其中 k + 1 j n 。满足

i = 1 + j = 1 n 0 T e λ i τ e A τ d τ a i j e j φ i ( x ) = i = 1 + j = 1 n e λ i t b i j e j φ i ( x ) . (3.2)

y ^ L 2 ( Ω ; n ) = ( k = 1 + j = 1 n a k j 2 ) 1 / 2 < ,其中 λ i σ ( Δ n ) 。用 b i 表示 ( b i 1 , , b i n ) ,则(3.1)有唯一解 y ^ 等价于存在唯一 { a i } i + L 2 ( + ; n ) 使得对每一 i +

0 T e λ i τ e A τ d τ a i = e λ i t b i . (3.3)

对任意A存在非奇异矩阵P,使得 A ^ = P 1 ( A ) P 是A的若当形式,用 β j , j = 1 , 2 , , r 表示的A特征值,对每个j, β j 的代数重数记为 m j ,即

A ^ = ( J 1 J 2 J r ) ,

其中 J j 是对角线为 β j 的若当块 j = 1 , 2 , , r 。记 a ^ i = P 1 a i b ^ i = P 1 b i 。则存在唯一 { a i } i + L 2 ( + ; n ) 满足(3.3)等价于存在唯一 { a ^ i } i + L 2 ( + ; n ) 使得每一 i +

0 T e λ i τ e A ^ τ d τ a ^ i = e λ i t b ^ i . (3.4)

而矩阵 0 T e λ i τ e A ^ τ d τ 是上三角形式,又 T T 2 ,故矩阵对角线上元素 0 T e ( λ i β j ) τ d τ 非零, j = 1 , 2 , , r 。因此对每一 i + ,有

det ( 0 T e λ i τ e A ^ τ d τ ) = j = 1 r ( 0 T e ( λ i β j ) τ d τ ) m j 0. (3.5)

这意味着矩阵 0 T e λ i τ e A ^ τ d τ 可逆,由此可得对每一 i + ,存在唯一 a ^ i 满足方程(3.4),即

a ^ i = ( 0 T e λ i τ e A ^ τ d τ ) 1 e λ i t b ^ i . (3.6)

接下来证明 y ^ L 2 ( Ω ; n ) < ,即存在常数 C > 0 使得

i = 1 + P a ^ i n 2 C i = 1 + P b ^ i n 2 . (3.7)

这里 y ^ L 2 ( Ω ; n ) 2 = i = 1 + P a ^ i n 2 y ¯ 1 L 2 ( Ω ; n ) 2 = i = 1 + P b ^ i n 2 。易知矩阵 0 T e λ i τ e A ^ τ d τ 的元素是 0 T e ( λ i β j ) τ d τ 对于某些 j = { 1 , 2 , , r } 或者是 1 k ! 0 T τ k e ( λ i β j ) τ d τ 对于某些 j = { 1 , 2 , , r } k { 1 , 2 , , m ¯ 1 } ,其中 m ¯ = max { m j , j = 1 , 2 , , r } 。存在与 i , j , k 无关的常数 C 1 ( σ ( A ) , T ) 1 C 2 ( σ ( A ) , T ) 1 使得对于所有 i + j = { 1 , 2 , , r } k { 1 , 2 , , m ¯ 1 }

1 k ! 0 T τ k e ( λ i β j ) τ d τ C 1 0 T e ( λ i β j ) τ d τ C 2 0 T e λ i τ d τ . (3.8)

定义集合 O i : = { i + , λ i β j = 0 , j { 1 , 2 , , r } } 。存在与 i , j , k 无关的常数 C 3 ( σ ( A ) , T ) 1 ,使得对于所有 i + \ O i j = { 1 , 2 , , r }

0 T e ( λ i β j ) τ d τ C 3 0 T e λ i τ d τ . (3.9)

A j k 表示矩阵 0 T e λ i τ e A ^ τ d τ 的代数余子式。由(3.8)可知存在与 i , j , k 无关的常数 C 4 1 ,使得 A j k C 4 ( 0 T e λ i τ d τ ) n 1 1 j , k n 。由(3.9)可知存在与i无关的常数 C 5 1 ,使得对于所有 i + \ O i det ( 0 T e λ i τ e A ^ τ d τ ) C 5 ( 0 T e λ i τ d τ ) n 。因此对于所有 i + \ O i

( 0 T e λ i τ e A ^ τ d τ ) 1 ( j = 1 n k = 1 n A j k 2 ( det 0 T e λ i τ e A ^ τ d τ ) 2 ) 1 / 2 C 4 n C 5 0 T e λ i τ d τ . (3.10)

对于 i O i ,可知存在与i无关的常数 C 6 1 ,使得 det ( 0 T e λ i τ e A ^ τ d τ ) C 6 ( 0 T e λ i τ d τ ) n 1 ,故对于 i O i

( 0 T e λ i τ e A ^ τ d τ ) 1 C 4 n C 6 . (3.11)

综上,我们由(3.6),(3.10)和(3.11)推得使不等式(3.7)成立只需取

C = max { max i + \ O i C 4 n e λ i t C 5 0 T e λ i τ d τ , max i O i C 4 n e λ i t C 6 } < .

而不等式(3.7)成立意味着

y ^ L 2 ( Ω ; n ) 2 C P P 1 y ¯ 1 L 2 ( Ω ; n ) 2 < . (3.12)

由此得映射 ( 0 T e Δ n A τ d τ ) 1 e Δ n t : y ¯ 1 y ^ L 2 ( Ω ; n ) 上有界,于是我们证明了存在唯一 y ^ 是(3.1)的解。又对于 t > 0 e A t 可逆,故对任意 y ¯ 1 L 2 ( Ω ; n ) ,存在唯一 y ¯ L 2 ( Ω ; n ) 使得

y ¯ 1 = e A t y ¯ .

结合(3.12)我们得到 y ^ L 2 ( Ω ; n ) 2 C P P 1 e A t y ¯ L 2 ( Ω ; n ) 2 < . 综上我们证明了引理3.4的(i)。

(ii) 同理可证引理3.4的(ii)成立。

引理3.5对任意 n × n 矩阵A, t > 0 T T 2 有下列等式成立

(i) 0 T e A τ d τ e A t = e A t 0 T e A τ d τ e A t ( 0 T e A τ d τ ) 1 = ( 0 T e A τ d τ ) 1 e A t

(ii) 0 T e A τ d τ e A t = e A t 0 T e A τ d τ e A t ( 0 T e A τ d τ ) 1 = ( 0 T e A τ d τ ) 1 e A t

引理3.5的证明:(i)由 A 的定义知 A = Δ n A Δ n = Δ I n × n I n × n 表示n阶单位阵。 A L 2 ( Ω ; n ) 上生成了一个 C 0 半群 { e A t } t 0 e A t = e Δ n t e A t 。于是对任意 n × n 矩阵A, z = ( z 1 , , z n ) T L 2 ( Ω ; n ) ,有

e Δ n t A z = ( e Δ t e Δ t ) ( a 11 a 1 n a n 1 a n n ) ( z 1 z n ) = ( j = 1 n a 1 j e Δ t z 1 j = 1 n a n j e Δ t z n ) .

A e Δ n t z = ( a 11 a 1 n a n 1 a n n ) ( e Δ t e Δ t ) ( z 1 z n ) = ( j = 1 n a 1 j e Δ t z 1 j = 1 n a n j e Δ t z n ) .

其中 a i j ( i = 1 , , n ; j = 1 , , n ) 表示A的元素,由此可得 e A t e Δ n t = e Δ n t e A t 。故有

0 T e A τ d τ e A t = 0 T e Δ n τ e A τ d τ e Δ n t e A t = e Δ n t e A t 0 T e Δ n τ e A τ d τ = e A t 0 T e A τ d τ .

有上述等式成立则有

e A t ( 0 T e A τ d τ ) 1 = ( 0 T e A τ d τ ) 1 0 T e A τ d τ e A t ( 0 T e A τ d τ ) 1 = ( 0 T e A τ d τ ) 1 e A t 0 T e A τ d τ ( 0 T e A τ d τ ) 1 = ( 0 T e A τ d τ ) 1 e A t .

(ii)与(i)同理可证,接下来我们依次证明定理2.1的(i)和(ii)。

证明:(i) 当 Ω \ ω ¯ 时,因为对任意 T > 0 T T 1 T T 2 ,故第一步证明此时对任意 T T 1 系统(2.1)不是T-零能控的,第二步证明此时对任意 T T 2 系统(2.1)也不是T-零能控的。第一步,任取采样周期 T T 1 ,故此时存在 i + j ( 1 , , r ) ,使得

0 T e ( λ i β j ) τ d τ = 0.

结合引理(3.5)可知对此时的 i + j ( 1 , , r )

det ( 0 T e λ i τ e A ^ τ d τ ) = j = 1 r ( 0 T e ( λ i β j ) τ d τ ) m j = 0.

这意味着存在 φ L 2 ( Ω ; n ) \ { 0 } 使得

0 T e A τ d τ φ = 0.

这意味着存在 φ L 2 ( Ω ; n ) \ { 0 } 使得

0 T e A τ d τ φ = 0. (3.13)

采用反证法,假设存在采样周期 T T 1 使系统(2.1)是T-零能控的。由定义2.1知,对于任意

y 0 L 2 ( Ω ; n ) ,存在 p + { u k } k = 1 p L 2 ( Ω ; m ) 使

y ( p T ; y 0 , { u k } k = 1 p ) = e A p T y 0 + k = 1 P e A ( p k ) T 0 T e A τ d τ χ ω B u k = 0 .

于是与(3.13)一起有

φ , e A p T y 0 + k = 1 P e A ( p k ) T 0 T e A τ d τ χ ω B u k = e A p T φ , y 0 + k = 1 p e A ( p k ) T 0 T e A τ d τ φ , χ ω B u k = e A p T φ , y 0 = 0.

由于 y 0 L 2 ( Ω ; n ) 为任取的,这意味着 e A p T φ = 0 ,再由引理3.3与热方程的倒向唯一性我们得 φ = 0 ,与 φ L 2 ( Ω ; n ) \ { 0 } 相矛盾,由此得出,当取采样周期 T T 1 时,系统(2.1)不是T-零能控的。注意到当取采样周期 T T 1 ,此时无论控制域 ω 是否是全空间 Ω ,系统(2.1)都不是T-零能控的。

第二步,采用反证法,我们假设此时系统(2.1)对于某采样周期 T T 2 是T-零能控的。那么根据定义2.1,任取 y 0 L 2 ( Ω ; n ) \ { 0 } ,存在 p + { u k } k = 1 p L 2 ( Ω ; m ) 使得

0 = y ( p T ; y 0 , { u k } k = 1 p ) = e A p T y 0 + k = 1 p ( k 1 ) T k T e A ( p T τ ) d τ χ ω B u k = e A p T y 0 + k = 1 p e A ( p k ) T 0 T e A τ d τ χ ω B u k = e A p T y 0 + k = 1 p 1 e A ( p k ) T 0 T e A τ d τ χ ω B u k + 0 T e A τ d τ χ ω B u P = 0 T e A τ d τ [ ( 0 T e A τ d τ ) 1 e A T ( e A ( p 1 ) T y 0 + k = 1 p 1 e A ( p 1 k ) T 0 T e A τ d τ χ ω B u k ) + χ ω B u p ] .

结合引理3.4,其中

y ¯ = 0 , y ^ = ( 0 T e A τ d τ ) 1 e A T ( e A ( p 1 ) T y 0 + k = 1 p 1 e A ( p 1 k ) T 0 T e A τ d τ χ ω B u k ) + χ ω B u p .

可得

( 0 T e A τ d τ ) 1 e A T ( e A ( p 1 ) T y 0 + k = 1 p 1 e A ( p 1 k ) T 0 T e A τ d τ χ ω B u k ) + χ ω B u p = 0.

Ω \ ω ¯ ,存在 B r ( x 0 ) Ω 使得

B r ( x 0 ) Ω \ ω ¯ . (3.14)

这里 B r ( x 0 ) 表示在 N 上以 x 0 为中心以 r > 0 为半径的邻域,这意味着在 B r ( x 0 ) χ ω B u p = 0 ,结合引理3.5可知在 B r ( x 0 )

0 = ( 0 T e A τ d τ ) 1 e A T ( e A ( p 1 ) T y 0 + k = 1 p 1 e A ( p 1 k ) T 0 T e A τ d τ χ ω B u k ) = e A T 2 ( 0 T e A τ d τ ) 1 e A T 2 ( e A ( p 1 ) T y 0 + k = 1 p 1 e A ( p 1 k ) T 0 T e A τ d τ χ ω B u k ) .

g = ( 0 T e A τ d τ ) 1 e A T 2 ( e A ( p 1 ) T y 0 + k = 1 p 1 e A ( p 1 k ) T 0 T e A τ d τ χ ω B u k )

再由引理3.1的(ii),其中

ω 1 = B r ( x 0 ) , t = T 2 , z = g .

我们得到对几乎处处的 x Ω

g = 0.

由引理3.4和热方程倒向唯一性可知

e A ( p 1 ) T y 0 + k = 1 p 1 e A ( p 1 k ) T 0 T e A τ d τ χ ω B u k = 0. (3.15)

同理重复上述得到(3.15)的过程最终我们可得

e A T y 0 + 0 T e A τ d τ χ ω B u 1 = 0. (3.16)

由(3.16)结合引理3.5得

0 = e A T y 0 + 0 T e A τ d τ χ ω B u 1 = 0 T e A τ d τ [ e A T 2 ( 0 T e A τ d τ ) 1 e A T 2 y 0 + χ ω B u 1 ] .

再由引理3.4得

e A T 2 ( 0 T e A τ d τ ) 1 e A T 2 y 0 + χ ω B u 1 = 0.

结合(3.14)我们发现在 B r ( x 0 ) χ ω B u 1 = 0 ,再根据引理3.1的(ii)

此时 ω 1 = B r ( x 0 ) t = T 2 z = ( 0 T e A τ d τ ) 1 e A T 2 y 0

得到

( 0 T e A τ d τ ) 1 e A T 2 y 0 = 0.

g ¯ = ( 0 T e A τ d τ ) 1 e A T 2 y 0 ,于是有

e A T 2 y 0 = 0 T e A τ d τ g ¯ = 0.

进一步由倒向唯一性我们得出 y 0 = 0 ,这与 y 0 L 2 ( Ω ; n ) \ { 0 } 相矛盾,因此假设不成立,故当 Ω \ ω ¯ 时,系统(2.1)对于任意非病态周期 T T 2 不是T-零能控的。定理2.1的(i)得证。

(ii) 当 ω = Ω 时,我们先证明其充分性,为此,任取 T T 2 ,并设 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n 。我们定义如下 n 的子空间 V k k { 0 , 1 , , n 1 }

V k = { e A k T B a n : a m } . (3.17)

由假设和(3.17)得 n = k = 0 n 1 V k ,因此存在线性映射 P k : n V k 使对任意 a n a = k = 0 n 1 P k a 。从 V k 的定义可以看出对任意 k { 0 , 1 , , n 1 } ,有

e A k T V k R a n g e B . (3.18)

因此根据的 P k 定义与(3.18)可知对于 k { 0 , 1 , , n 1 } ,任意 z L 2 ( Ω ; n )

e A k T P k z R a n g e B . (3.19)

由引理3.1的(i),此时

ω 1 Ω , V = R a n g e B , t = k T , z = e A k T P k z .

我们可推得对于几乎处处的 x Ω ,当 k { 0 , 1 , , n 1 }

e Δ n k T e A k T P k z R a n g e B . (3.20)

因为映射 B : m R a n g e B 是满射,故存在线性映射 C : n m 使得对于所有 a R a n g e B 成立

B C a = a . (3.21)

结合(3.20)和(3.21)可知对于任意 k { 0 , 1 , , n 1 } ,取 u ( k + 1 ) = C e Δ n k T e A k T P k z 可使

e Δ n k T e A k T P k z + B u ( k + 1 ) = 0. (3.22)

ω = Ω ,任取 y 0 L 2 ( Ω ; n ) ,结合(3.21),(3.22)和引理3.5我们得出

y ( n T ; y 0 , { u k } k = 1 n ) = e A n T y 0 + k = 1 n e A ( n k ) T 0 T e A τ d τ B u k = 0 T e A τ d τ [ e A ( n 1 ) T ( 0 T e A τ d τ ) 1 e A T y 0 + k = 1 n e A ( n k ) T B u k ] .

其中令 y ˜ 0 = ( 0 T e A τ d τ ) 1 e A T y 0 ,因 y 0 L 2 ( Ω ; n ) 由引理3.4知 y ˜ 0 L 2 ( Ω ; n ) 于是

y ( n T ; y 0 , { u k } k = 1 n ) = 0 T e A τ d τ [ e A ( n 1 ) T y ˜ 0 + k = 1 n e A ( n k ) T B u k ] = 0 T e A τ d τ [ k = 1 n e A ( n 1 ) T P ( k 1 ) y ˜ 0 + k = 1 n e A ( n k ) T B u k ] = 0 T e A τ d τ k = 1 n [ e A ( n 1 ) T P ( k 1 ) y ˜ 0 + e A ( n k ) T B u k ] = 0 T e A τ d τ k = 1 n [ e A ( n k ) T ( e A ( l 1 ) T P ( k 1 ) y ˜ 0 + B u k ) ] = 0

此时取 u k = C e A ( k 1 ) T P k y ˜ 0 ,由定义2.1我们可知系统(2.1)是T-零能控的,充分性得证。

下面我们证明其必要性,设系统(2.1)是T-零能控的,在定理2.1的(i)中证明了当 T T 1 时,系统(2.1)不是T-零能控的,这意味着若系统(2.1)是T-零能控的,则 T T 2 。于是我们现在只需说明:若系统(2.1)是T-零能控的,那么则有 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n 。为此采用反证法,假设 r a n k ( B , e A T B , , e A ( n 1 ) T B ) < n ,于是根据引理3.2,又 ω = Ω ,存在 z ¯ L 2 ( Ω ; n ) \ { 0 } 使得对于任意 t > 0

B e A t z ¯ = 0 (3.23)

因系统(2.1)是T-零能控的,根据定义2.1,则对任意 z ^ ,存在 p + { u k } k = 1 p L 2 ( Ω ; m ) 使得

0 = y ( p T ; z ^ , { u k } k = 1 p ) = e A p T z ^ + k = 1 P e A ( p k ) T 0 T e A τ d τ B u k . (3.24)

其中取 z ^ = ( 0 T e A τ d τ ) 1 e A ( p + 1 ) T z ¯ ,则有

( 0 T e A τ d τ ) 1 e A ( p + 1 ) T z ¯ 2 = e A p T [ ( 0 T e A τ d τ ) 1 e A ( p + 1 ) T z ¯ ] , ( 0 T e A τ d τ ) 1 e A T z ¯ = k = 1 p e A ( p k ) T 0 T e A τ d τ B u k , ( 0 T e A τ d τ ) 1 e A T z ¯ = k = 1 p u k , B e A ( p k + 1 ) T z ¯ .

结合(3.23)意味着

( 0 T e A τ d τ ) 1 e A ( p + 1 ) T z ¯ 2 = 0. (3.25)

根据引理3.4和(3.25)我们可以得到 e A p T z ¯ = 0 ,再由热方程倒向唯一性得出 z ¯ = 0 而这与 z ¯ L 2 ( Ω ; n ) \ { 0 } 矛盾。由此可知假设错误,故 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n 。必要性得证,综上我们证明了定理2.1。

4. 定理2.2的证明

证明定理2.2之前我们先给出以下两个引理。

引理4.1设 T > 0 { k T } k = 1 p ( p + ) 则有以下两条结论等价

(i) 系统(2.1)是T-近似能控的。

(ii) 若 z 0 L 2 ( Ω ; n ) ,则

χ ω B 0 T e A τ d τ e A ( p k ) T z 0 = 0 对于所有 k { 1 , , p } z 0 = 0 Ω 上。

根据近似能控性与唯一延拓性之间的等价关系结合系统(2.1)的解(2.2)我们容易推得引理4.1,这里省略了具体证明。

引理4.2设 T > 0 满足

r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n .

则有

k = 1 p ker ( B e A ( p k ) T ) = { 0 } .

引理4.2的证明参见文献 [6] 定理5.2的证明。有了以上理论基础,下面我们证明定理2.2。

证明:我们先证明其充分性,设周期 T T 2 且设 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n ,由引理4.2有

k = 1 n ker ( B e A ( n k ) T ) = { 0 } . (4.1)

z L 2 ( Ω ; n ) 满足

χ ω B 0 T e A τ d τ e A ( n k ) T z = 0 对于所有 k { 1 , , n } (4.2)

根据引理3.5与(4.2)这意味着对所有 k { 1 , , n }

χ ω B 0 T e A τ d τ e A ( n k ) T z = B e A ( n k ) T χ ω 0 T e A τ d τ e Δ n ( n k ) T z = 0. (4.3)

由(4.3)推出对所有 k { 1 , , n }

χ ω 0 T e A τ d τ e Δ n ( n k ) T z ker ( B e A ( n k ) T ) . (4.4)

再由引理3.5有

χ ω 0 T e A τ d τ e Δ n ( n k ) T z = χ ω e Δ n ( n k ) T 0 T e A τ d τ z .

再根据(4.4)与引理3.1的(i)一起对所有 k { 1 , , n } ,其中

ω 1 = ω , V = ker ( B e A ( n k ) T ) , t = ( n k ) T .

我们可知几乎处处的 x Ω ,所有 k { 1 , , n }

0 T e A τ d τ z ker ( B e A ( n k ) T ) .

再由(4.1)得出

0 T e A τ d τ z = 0.

由引理3.4我们得到 z = 0 ,进一步由引理4.1我们得出系统(2.1)是T-近似能控的,充分性得证。

下面我们证明其必要性,设系统(2.1)是T-近似能控的,第一步先说明此时 T T 2 ,采用反证法,假设 T T 1 ,故此时存在 z ^ L 2 ( Ω ; n ) \ { 0 } 使得

0 T e A τ d τ z ^ = 0. (4.5)

再由引理3.5得对 p + k { 1 , , p }

χ ω B 0 T e A τ d τ e A ( p k ) T z ^ = χ ω B e A ( p k ) T 0 T e A τ d τ z ^ = 0.

而系统(2.1)是T-近似能控的,结合引理4。1我们推出 z ^ = 0 ,而这与 z ^ L 2 ( Ω ; n ) \ { 0 } 矛盾,故此时 T T 2 。第二步说明此时 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n ,仍采用反证法,假设此时 r a n k ( B , e A T B , , e A ( n 1 ) T B ) < n 。于是由引理3.2可知存在 z L 2 ( Ω ; n ) \ { 0 } 使得 χ ω B e A t z = 0 对于任意 t > 0 ,特殊的有

χ ω B e A k T z = 0 对于任意 k +

任取 T T 2 ,由引理3.4可知存在 z ¯ L 2 ( Ω ; n ) 使得

e A T z = 0 T e A τ d τ z ¯ . (4.6)

于是对任意 k +

χ ω B e A k T z = χ ω B e A ( k 1 ) T e A T z = χ ω B e A ( k 1 ) T 0 T e A t d t z ¯ = 0. (4.7)

由(4.7)结合引理3.5可推得对 p +

χ ω B 0 T e A τ d τ e A ( p k ) T z ¯ = 0 对于所有 k { 1 , , p }

根据假设系统(2.1) T-近似能控,由引理4.1则

χ ω B 0 T e A τ d τ e A ( p k ) T z ¯ = 0 对于所有 k { 1 , , p } z ¯ = 0 Ω 上。

再由(4.6)得 e A T z = 0 ,进一步根据热方程倒向唯一性可得 z = 0 ,与 z L 2 ( Ω ; n ) \ { 0 } 矛盾。这意味着 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n ,必要性得证。至此我们证明了本文的主要定理。

5. 结论

关于主要结论下面给出一些说明:

1) 本文中定义的T-零能控性和T-近似能控性均与采样周期T有关,即考虑的是在给定周期T选择下系统(2.1)的能控性问题。

2) 定理2.1说明系统(2.1)是否T-零能控与控制域 ω 是否为全空间 Ω 有关,若 ω Ω 的真子集时,则对任意采样周期T,系统(2.1)都不是T-零能控的。

3) 当 ω = Ω 时,我们得到了系统(2.1)是T-零能控的一个充要条件,此条件为 T T 2 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n 。在本文定理(2.1)的证明中我们说明了当 T T 1 时,无论控制域 ω 是否为全空间 Ω ,系统(2.1)都不是T-零能控的。

4) 在定理2.2中我们给出了系统(2.1)是T-近似能控的一个充要条件,也是 T T 2 r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n 。与系统(2.1) T-零能控相比系统(2.1) T-近似能控并不依赖于控制域 ω 的大小。

5) 推论2.3分别给出了系统(2.1)是T-零能控的和T-近似能控的充分条件。由定义2.3和定义2.4易知,若采样周期T是非病态的,则 T T 2 ,而采样周期T非病态与 ( A , B ) 满足 r a n k ( B , A B , , A n 1 B ) = n r a n k ( B , e A T B , , e A ( n 1 ) T B ) = n 的一充分条件(参见文献 [11] 的第四章),故推论2.3是定理2.1和定理2.2的特殊情形,因此推论2.3给出的也是一充分条件。

基金项目

湖北省自然科学基金面上项目支持(项目编号:2023AFB622)。

参考文献

[1] 郭戈, 王伟. 保证能控性和能观性的采样方式[J]. 控制与决策, 2003, 18(1): 66-68.
[2] 郭戈. 周期性系统的能控性理论[J]. 控制与决策, 2004, 19(4): 468-470.
[3] 杨磊, 李俊民. 一类线性切换系统的能控性和能观测性的充要条件[J]. 系统工程与电子技术, 2003, 25(5): 588-590.
[4] 王俭, 肖金球, 施颂椒. 离散时间系统能控能观测与信号重构的关系[J]. 电路与系统学报, 2003, 8(6): 136-138.
[5] Chen, T. and Francis, B.A. (2012) Optimal Sampled-Data Control Systems. Springer Science & Business Media, Berlin.
[6] Qin, S. and Wang, G. (2017) Controllability of Impulse Controlled Systems of Heat Equations Coupled by Constant Matrices. Journal of Differential Equations, 263, 6456-6493.
https://doi.org/10.1016/j.jde.2017.07.018
[7] Liu, X. (1995) Impulsive Control and Optimization. Applied Mathematics and Computation, 73, 77-98.
https://doi.org/10.1016/0096-3003(94)00204-H
[8] Wang, L., Yan, Q. and Yu, H. (2021) Constrained Approximate Null Controllability of the Coupled Heat Equation with Impulse Controls. SIAM Journal on Control and Optimization, 59, 3418-3446.
https://doi.org/10.1137/20M1338393
[9] Wang, G., Yang, D. And Zhang, Y. (2017) Time Optimal Sampled-Data Controls for the Heat Equation. Comptes Rendus Mathematique, 355, 1252-1290.
https://doi.org/10.1016/j.crma.2017.11.006
[10] Lin, P., Liu, H. and Wang, G. (2020) Output Feedback Stabilization for Heat Equations with Sampled-Data Controls. Journal of Differential Equations, 268, 5823-5854.
https://doi.org/10.1016/j.jde.2019.11.019
[11] 郑大钟. 线性系统理论[M]. 第2版. 北京: 清华大学出版社, 2002.