函数的水平集
The Level Sets of Functions
DOI: 10.12677/AAM.2021.107237, PDF, HTML, XML, 下载: 515  浏览: 2,204  国家自然科学基金支持
作者: 沈祖沛, 李海泉*:广东金融学院金融数学与统计学院,广东 广州
关键词: 水平集函数性质Level Sets Properties of Function
摘要: 本文利用函数的水平集给出函数的一些基本性质的刻画,包括函数凸性、连续性、可测性、可积性以及一些重要的应用。我们发现,在一定意义下,函数的性质完全由其水平集决定,而且方法上较为统一和直观。
Abstract: In this paper, we give some properties of functions by level sets of functions, including convexity, continuity, measurable, integrability and some important applications. We find that the properties of a function are completely determined by its level set, and the methods are uniform and intuitive.
文章引用:沈祖沛, 李海泉. 函数的水平集[J]. 应用数学进展, 2021, 10(7): 2278-2284. https://doi.org/10.12677/AAM.2021.107237

1. 引言

数学分析是一门很基础且很重要的专业学科,数学分析研究的对象是函数,主要研究函数(包括泛函)的一些基本性质,包括凸性、连续性、可微性、可积性以及一些积分估计等。研究这些性质的方法多种多样,国内也有很多很优秀经典的教材 [1] [2] [3] [4]。水平集方法在现实中有广泛的应用,例如图像分割 [5],计算机视觉 [6] 等领域。然而,介绍水平集方法的几乎都是相应学科的专著。鉴于水平集的重要性,本文主要通过总结和归纳由函数的水平集给出函数性质的刻画。我们发现,在一定意义下,函数的性质完全由其水平集决定,方法上较为统一和直观,而且可读性强。然而,本科数学分析的教材很少涉及这方面的内容。这是写这篇文章的主要动机和意义。

2. 一些定义

我们首先给出一些基本定义,这些定义在正文中经常用到。

定义2.1. 设 f : Ω R Ω 上的实值函数,对每个 t ,我们定义函数上水平集为

S f ( t ) = { x Ω : f ( x ) > t } .

定义函数的下水平集为

L f ( t ) = { x Ω : f ( x ) t } .

上水平集和下水平集我们统称函数的水平集。

定义 2.2. 一个由 Ω 的某些子集构成的集合 Σ ,如果满足下面几条公理,则称它是一个 σ 代数:

i) 若 A Σ ,则 A c Σ

ii) A 1 , A 2 Σ 中的一可数的集合族,则它的可数并 i = 1 A i Σ

iii) Ω Σ

定义2.3. 测度空间是由三部分组成:集合 Ω σ - 代数 Σ 及测度 μ 。我们通常把 Σ 中的元素称为可测集

定义2.4. 设 A R n 是Lebesgue 测度有限的Borel集(由开集生成的最小 σ 代数)。我们定义集合 A 的对称重排 A

A = { x : | x | < r } , ( | S n 1 | / n ) r n = L n ( A ) ,

这里, | S n 1 | 表示 S n 1 的表面积。

定义2.5. 设 A 是一个集合,定义集合 A 的特征函数为

χ A ( x ) = { 1 x A 0 x A

3. 函数的一些性质的水平集刻画

这一部分我们主要从函数的水平集的角度给出函数的一些基本性质的刻画。

3.1. 函数本身的刻画

我们可以利用函数的水平集的直接构造出函数本身。

性质3.1.1 (层饼表示定理)设 f ( x ) 是任意的非负可测函数,则有

f ( x ) = 0 χ { f > t } ( x ) d t .

证明:这是显然的,实际上

0 χ { f > t } ( x ) d t = 0 f ( x ) 1 d t + f ( x ) 0 d t = f ( x ) .

这个性质在实际应用中极其有用,因为它允许我们在很多情况下把处理一个普通函数的积分简化为处理特征函数的积分(利用Fubini定理)。

3.2. 函数的凸性

我们首先给出凸集和凸函数的定义。

定义3.2.1. 我们称集合 A 为凸集,若对任意的 x , y A , θ [ 0 , 1 ] ,有

θ x + ( 1 θ ) y A .

定义3.2.2. 若 f 的定义域为凸集,我们称 f 是凸函数,若其满足

f ( θ x + ( 1 θ ) y ) θ f ( x ) + ( 1 θ ) f ( y ) x , y d o m f θ [ 0 , 1 ] .

这里 d o m f 表示函数的定义域。

由函数凸性的定义,我们马上得到凸函数的必要条件:

性质3.2.1.凸函数的下水平集是凸集,但是反之不一定成立。

证明:设任意的的 t x , y L f ( t ) ,我们有

f ( θ x + ( 1 θ ) y ) θ f ( x ) + ( 1 θ ) f ( y ) θ t + ( 1 θ ) t t .

因此

θ x + ( 1 θ ) y L f ( t ) .

另一方面,我们很容易验证函数 f ( x ) = e x 的任意下水平集是凸的,但是这个函数不是凸的。然而,这个性质的可以判断一个集合是否是凸集:若某个集合可以描述为一个凸函数的下水平集,那么这个集合为凸集。

这个结论虽然不能从函数的水平集的凸性直接得到函数的凸性,但是我们可以从函数的图像(Epigraph)的角度上看,可以得到函数凸性的一个充分必要条件。

定义3.2.3. 设 f : R n R 。我们定义 f 的图像(Epigraph)为

E p i f = { ( x , t ) R n + 1 | x d o m f , f ( x ) t } .

由这个定义,我们可以得到函数的凸性的一个充要条件:

性质3.2.2. 函数 f 为凸函数当且仅当其图像为凸集。

证明:根据定义,这是显然的。

注:对于凹函数,相关的性质也可以用水平集刻画,只要把下水平集改为上水平集即可。

3.3. 函数的连续性

因为 { ( , b ) : b R } { ( a , + ) : a R } R 上的一组拓扑子基,我们很容易验证:

性质3.3.1. 对应任意的 t ,若函数的水平集 S f ( t ) = { x Ω : f ( x ) > t } 是开的,则函数 f ( x ) 是下半连续的。

类似的,我们可以得到函数的上半连续的刻画:

性质3.3.2. 对应任意的 t ,若函数的水平集 L ˜ f ( t ) = { x Ω : f ( x ) < t } 是开的,则函数 f ( x ) 是上半连续的。

性质3.3.3. 如果一个函数既是上半连续,也是下半连续,那么它一定连续。

3.4. 函数的可测性

我们可以直接根据函数的水平集给出可测函数的定义:

定义3.4.1.对应任意的 t ,若函数的上水平集 S f ( t ) = { x Ω : f ( x ) > t } 是可测集,则称函数 f ( x ) 是可测的。更一般地,如果 f : Ω C 是复函数,那么当它的实部 Re f 和虚部 Im f 都是可测函数时,称 f 可测。

从水平集的角度去定义函数的连续和可测性,我们可以得到可测函数与连续函数的某些关系。例如,如果 Σ R n 上的Borel- σ 代数B,则显然每个连续函数都是Borel可测的,此时函数的水平集

S f ( t ) = { x Ω : f ( x ) > t } 是开集。实际上,上半连续或下半连续函数都是Borel可测的。下半连续函数的可测性是显然的,关于上半连续函数的可测性,因为集合 { x | f ( x ) < t + 1 i } 是可测的,同时

{ x Ω | f ( x ) t } = i = 1 { x | f ( x ) < t + 1 i } ,

所以 { x Ω | f ( x ) t } 是可测的,从而可以得到上半连续的函数的可测性。

通过上面的分析,我们如果从水平集的角度研究函数的连续性与可测性之间的关系,实际上建立一种 σ 代数和拓扑之间的关系的。从这个角度看,这是有意思的。

3.5. 函数的可积性

为了简单起见,本部分都假设 f 是非负的可测函数,且直接讨论一般测度 μ 上的积分。函数的可积性只取决于函数大的值的集合究竟有多大,而以这些大的值在哪里取到没有什么关系。因此,我们很自然可以考虑对应大的 t ,水平集的测度究竟多大。由此直接引出可测函数的分布函数的定义。我们记:

F f ( t ) = μ ( S f ( t ) )

也就是说 F f ( t ) 是使得 f > t 的集合的测度。我们把这个函数称为 f 的分布函数 [7]。显然, F f ( t ) t 是单调不增的。因此我们可以用Riemann积分 0 F f ( t ) d t 来定义 f Ω 上的积分:

.

利用这个定义,我们把层饼表示定理推广为:

性质3.5.1. [8] 设 v 为正实半轴上Borel集上的测度,使得 φ ( t ) = v ( [ 0 , t ) ) 对每个 t > 0 均有限。再设 ( Ω , Σ , μ ) 为一个测度空间且 f ( x ) 是任意的非负可测函数,那么有

证明:为了读者的方便,我们把文献 [8] 中关于此定理的证明摘录于此。根据定义,我们有

0 μ ( S f ( t ) ) v ( d t ) = 0 Ω χ f > t ( x ) μ ( d x ) v ( d t ) .

利用Fubini定理,(3.5)的右端等于

Ω ( 0 χ { f > t } ( x ) v ( d t ) ) μ ( d x ) .

注意到

0 χ { f > t} ( x ) v ( d t ) = 0 f ( x ) v ( d t ) = φ ( f ( x ) ) ,

结论成立。

上述的结论说明函数的可积性完全由其分布函数决定。因此,通过研究函数的水平集,我们完全可以弄清楚可测函数的可积性,而且相对比较自然和直观。

与分布函数密切相关的一个概念:函数的重排。直观上看,分布函数与函数的重排是一对“互逆”函数。函数的可积性完全可以由他们决定。有意思的是,通过对函数重排,得到一些非常有用的积分不等式。例如,根据重排不等式,我们可以得到著名的等周不等式。为了定义函数的重排,我们首先给出特征函数的重排的定义,即:

χ A = χ A

现在假设 f 是一个无穷远处趋于零的Borel可测函数,定义 f 的重排函数为

f ( x ) = 0 χ { | f | > t } ( x ) d t .

函数的重排一个重要的性质就是重排函数 f 与函数 f 就有相同的分布函数,因此,函数的可积性完全由他们决定,这些概念都是从函数的水平集的基础上给出,而且具有很好的性质。从这个方面看,用水平集研究函数的性质具有一定的优越性。

4. 水平集的一些应用

4.1. 几个基本定理

根据层饼表示定理,我们很容易把数学分析中的黎曼积分中极限与积分交换次序的几个重要定理推广到一般的测度空间上。

定理4.1.1 [8] (单调收敛定理)设 f 1 , f 2 , , f n 是测度空间 ( Ω , Σ , μ ) 上递增的可测函数列,且

lim n f n ( x ) = f ( x ) , μ -a . e . x Ω ,

lim n Ω f n ( x ) μ ( d x ) = Ω f ( x ) μ ( d x ) .

证明:不失一般性,我们假设 f n 非负。显然, S f ( t ) = n = 1 S f n ( t ) ,由测度的连续性,我们可以得到

lim n F f n ( t ) = F f ( t ) .

根据一般测度空间上的可测函数积分的定义,我们只需要验证Riemann积分满足

lim n 0 F f n ( t ) d t = 0 F f ( t ) d t .

由于每个分布函数 F f n ( t ) 关于 t n 都是单调的,所以结论显然成立。证明完毕。

作为单调收敛定理的推论,我们很容易得到在一般测度空间的上的Fatou引理和控制收敛定理。

4.2. 水平集在变分法中的应用

无穷维的临界点理论在在微分方程解的存在性的证明中起到很重要的作用,其已经发展成为一门很重要的独立分支。而建立Morse理论、山路定理及更一般的环绕定理等临界点理论,一个主要的思想是主要考虑了泛函水平集的拓扑结构的变化。设 X 是Banach空间, f C 1 ( X , R )

κ = { u X | f ( u ) = 0 } K c = κ f 1 ( c )

定义4.2.1. 设 c R ,称 { u n } X f C 1 ( X , R ) ( P S ) c 序列,若

f ( u n ) c , f ( u n ) 0.

f 满足 ( P S ) c 条件,若它的任意 ( P S ) c 序列都存在收敛子列。若对任意的 c R f 都满足 ( P S ) c 条件,则称 f 满足 ( P S ) 条件。

下面的定理表明,若水平集 f 1 [ a , b ] 中没有 f 的临界点,则 L f ( a ) L f ( b ) 的强形变收缩核。具体证明可以参考文献 [9]。

定理4.2.1.设 f C 1 ( X , R ) 满足(PS)条件,若 f 1 [ a , b ] κ = ,则存在 η C ( [ 0 , 1 ] × L f ( a ) , L f ( b ) ) 使得

i) η ( 0 , · ) = χ L f ( b ) , η ( 1 , L f ( b ) ) L f ( a )

ii) η ( t , x ) = x , ( t , x ) [ 0 , 1 ] × L f ( a )

a f 的临界点,上述定理可以推广为

定理4.2.2. 设 f C 1 ( X , R ) 满足(PS)条件, f 只有有限的临界点。 a f [ a , b ) 上唯一的临界值。则存在 η : [ 0 , 1 ] × ( L f ( a ) \ L f ( b ) ) L f ( b ) \ κ b 使得

i) η ( 0 , · ) = χ L f ( b ) \ κ b

ii) η ( t , x ) = x , ( t , x ) [ 0 , 1 ] × L f ( a )

iii) η ( 1 , L f ( b ) \ κ b ) L f ( a )

注:这个定理是建立Morse不等式的关键。

5. 结论

本文研究函数水平集对函数性质刻画,同时也给出了函数水平集一些简单但重要的应用。在一定意义下,函数的性质完全由其水平决定,而且方法上较为统一和直观。另外一方面,函数水平集在计算机视觉理论,图像分割等实际问题当中都有着广泛的应用,鉴于这两点,我们建议在本科数学专业加入关于水平集这一部分的内容。

基金项目

本文受到国家自然科学资金青年资金项目(11701114)资助。

NOTES

*通讯作者。

参考文献

[1] 常庚哲, 史济怀. 数学分析教程[M]. 北京: 高等教育出版社, 2003.
[2] 陈纪修, 於崇华, 金路. 数学分析[M]. 北京: 高等教育出版社, 1999.
[3] 梅加强. 数学分析[M]. 第二版. 北京: 高等教育出版社, 2019.
[4] 华东师范大学数学系. 数学分析[M]. 第四版. 北京: 高等教育出版社, 2010.
[5] 朱付平, 田捷, 林瑶, 葛行飞. 基于Level Set方法的医学图像分割[J]. 软件学报, 2002, 13(9): 1866-1872.
[6] 钱芸, 张英杰. 水平集的图像分割方法综述[J]. 中国图象图形学报, 2008, 13(1): 7-13.
[7] 程民德, 邓东皋, 龙瑞麟. 实分析[M]. 第二版. 北京: 高等教育出版社, 2008.
[8] Ellitt H. Lieb, Michael Loss. 分析学[M]. 第二版. 王斯雷, 译. 北京: 高等教育出版社, 2006.
[9] 张恭庆. 变分法讲义[M]. 北京: 高等教育出版社, 2011.