CN104063852A - 一种基于指数化核范数与混合奇异值截断的张量恢复方法 - Google Patents
一种基于指数化核范数与混合奇异值截断的张量恢复方法 Download PDFInfo
- Publication number
- CN104063852A CN104063852A CN201410321348.0A CN201410321348A CN104063852A CN 104063852 A CN104063852 A CN 104063852A CN 201410321348 A CN201410321348 A CN 201410321348A CN 104063852 A CN104063852 A CN 104063852A
- Authority
- CN
- China
- Prior art keywords
- tensor
- nuclear norm
- matrix
- singular value
- order
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于指数化核范数与混合奇异值截断的张量恢复方法,主要包含以下步骤:首先提出一种新的张量秩定义,张量不同模态下展开矩阵秩的最大值,采用核范数指数和的对数来逼近该张量秩定义,将其转化为凸函数;其次为了消除张量不同模态下展开的矩阵的相关性,引入一系列辅助变量来代替展开矩阵,并将约束条件利用拉格朗日乘子法转化为增广拉格朗日函数;最后采用交替方向法对增广拉格朗日函数中各类变量进行迭代优化,直到收敛。其中,对于核范数的指数和中的优化变量,本发明是一种通用的方法,相对于其它的经典张量恢复方法,该方法能够更好的描述高维数据的内在结构,从而获得更好的恢复结。
Description
技术领域
本发明涉及计算机模式识别技术领域,具体涉及一种基于指数化核范数与混合奇异值截断的张量恢复方法。
背景技术
张量恢复(tensor completion),即高维矩阵的恢复问题,对于一个部分元素缺失的高维矩阵,通过观察其已有位置的元素,从而恢复出缺失部分元素的一般性问题,是计算机模式识别领域中的研究热点之一,被广泛应用于图像去噪、图像恢复、信息推荐系统等众多领域。总的来说,大多数现有的张量恢复方法是基于低秩结构假设,即要求恢复的缺失元素使得整个张量的秩尽可能的小。目前有两种传统的定义张量秩的方法:基于张量的CP(CANDECOMP/PARAFAC)分解方法(CP秩)和Tucker分解方法(Tucker秩)。具体来说,CP秩可以定义为:用秩一张量(rank-one tensor)之和来表示给定张量需要的秩一张量的最小个数。Tucker秩可以定义为:不同模态下展开矩阵的秩的线性加权。无论是CP秩还是Tucker秩,最小化该张量秩的优化问题被证明是一个NP难问题。
为了解决上述问题,Gandy等人采用不同模态下展开矩阵秩的和来表示张量的秩,在计算过程中,将展开矩阵的秩用矩阵的核范数来近似代替。Signoretto等人提出一种Shatten-p范数来代替展开矩阵的核范数,并由此定义了张量的秩,最后讨论了该方法与核范数之间的关系。随后,Liu等人采用不同模态下展开矩阵核范数的线性加权来近似代替Tucker秩,并将该方法应用于图像恢复和医学图像去噪。最后,Tomioka等人对张量恢复方法进行了总结,认为有两种方式可以实现张量恢复:(1)将张量按某一个模态展开成二阶矩阵,可以将张量恢复问题转化为了二阶矩阵的恢复问题;(2)采用不同模态下展开矩阵核范数的线性加权来近似代替Tucker秩。
可以看出,上述方法的目的都是寻找张量Tucker秩的近似。然而Tucker秩定义的几何意义不清晰,并且张量不同模态下的权重难以选择,如果某个模态下展开矩阵的秩很大,而其对应的权重很小,那么上述定义将无法正确地描述张量的秩结构,从而导致张量低秩恢复效果不够理想。张量的CP秩是对矩阵秩的一个推广,它的几何意义比Tucker秩更为明确,然而直接优化CP秩是一个非常困难的问题,甚至连局部最优都难以获得。
发明内容
针对现有技术存在的不足,本发明的目的在于提供一种基于指数化核范数与混合奇异值截断的张量恢复方法,基于CP秩的思想,提出一种新的张量秩定义方式,该张量秩是CP秩的下界,并且能够有效地逼近CP秩,使得张量恢复过程对噪声更为鲁棒。此外,该定义没有任何权重参数需要设置,可以有效的消除参数影响。为了优化问题的可解性,我们用张量不同模态下展开矩阵核范数的指数和的对数来近似张量秩,并提出了一种混合奇异值截断算法来获得优化问题的解析解,从而实现快速有效的张量恢复。
为实现上述目的,本发明提供了如下技术方案:一种基于指数化核范数与混合奇异值截断的张量恢复方法,其特征在于:包括以下三个步骤:
(1)提出一种新的张量秩定义:张量不同模态下展开矩阵秩的最大值;该定义是张量CP秩的下界,能够有效的逼近CP秩,并消除了权重参数的影响,采用核范数指数和的对数来逼近该张量秩定义,将其转化为凸函数;
(2)为了消除张量不同模态下展开的矩阵的相关性,引入一系列辅助变量来代替展开矩阵,并将约束条件利用拉格朗日乘子法转化为增广拉格朗日函数;
(3)采用交替方向法对增广拉格朗日函数中各类变量进行迭代优化,直到收敛;其中,对于核范数的指数和中的优化变量,采用混合奇异值截断算法来获得解析解。
本发明进一步设置,所述的步骤(1)具体包括以下子步骤:
首先,根据张量CP秩和Tucker秩的优缺点,提出一种新的张量秩定义:张量展开矩阵秩的最大值;
其次,将展开矩阵的秩松弛为展开矩阵的核范数,并且利用核范数的指数和的对数来逼近最大值函数,从而将上述张量的秩定义转化为凸函数。
本发明还进一步设置,所述的步骤(2)具体包括以下子步骤:
首先,由于张量在不同模态下的展开矩阵具有相关性,引入一系列辅助矩阵变量来替换不同模态下的展开矩阵,并增加对应的约束条件;
其次,采用拉格朗日乘子法将所有约束条件加入到目标函数中,获得增广拉格朗日函数。
本发明还进一步设置,所述的步骤(3)具体包括以下子步骤:
首先,为了对增广拉格朗日函数中的不同变量进行分别优化,采用交替方向法对增广拉格朗日函数中的各类变量进行迭代优化;
其次,对于核范数指数中的优化变量,采用混合奇异值截断算法来获得解析解。
本发明的优点是:1、本发明所提出一种全新的并且通用的张量秩的定义方法,对任意张量数据都适用。并且,该定义是张量CP秩的下界,能够有效地逼近CP秩,并消除了Tucker秩中权重参数的影响。
2、本发明提出一种有效的方法对张量秩定义进行松弛:将展开矩阵的秩松弛为展开矩阵的核范数,并且利用展开矩阵核范数的指数和的对数来逼近最大值函数,从而将张量恢复中的目标函数转化为凸函数。
3、本发明采用辅助变量来消除张量不同模态下展开的矩阵的相关性,并采用拉格朗日乘子法将所有的约束条件转化为增广拉格朗日函数,从而减小了优化问题的复杂度。
4、本发明采用了一种交替方向法来实现对增广拉格朗日函数的分解迭代优化。对于核范数的指数和中的优化变量,本发明提出了一种混合奇异值截断算法来获得解析解,从而实现快速有效的张量恢复。
下面结合说明书附图和具体实施例对本发明作进一步说明。
附图说明
图1为本发明实施例的整体流程图。
具体实施方式
参见图1,本发明提供一种基于指数化核范数与混合奇异值截断的张量恢复方法,包括以下三个步骤:
(1)提出一种新的张量秩定义:张量不同模态下展开矩阵秩的最大值;该定义是张量CP秩的下界,能够有效的逼近CP秩,并消除了权重参数的影响,采用核范数指数和的对数来逼近该张量秩定义,将其转化为凸函数;
(2)为了消除张量不同模态下展开的矩阵的相关性,引入一系列辅助变量来代替展开矩阵,并将约束条件利用拉格朗日乘子法转化为增广拉格朗日函数;
(3)采用交替方向法对增广拉格朗日函数中各类变量进行迭代优化,直到收敛;其中,对于核范数的指数和中的优化变量,采用混合奇异值截断算法来获得解析解。
作为优选的,本实施例所述的步骤(1)具体包括以下子步骤:
首先,根据张量CP秩和Tucker秩的优缺点,提出一种新的张量秩定义:张量展开矩阵秩的最大值;
其次,将展开矩阵的秩松弛为展开矩阵的核范数,并且利用核范数的指数和的对数来逼近最大值函数,从而将上述张量的秩定义转化为凸函数。
作为优选的,本实施例所述的步骤(2)具体包括以下子步骤:
首先,由于张量在不同模态下的展开矩阵具有相关性,引入一系列辅助矩阵变量来替换不同模态下的展开矩阵,并增加对应的约束条件;
其次,采用拉格朗日乘子法将所有约束条件加入到目标函数中,获得增广拉格朗日函数。
作为优选的,本实施例所述的步骤(3)具体包括以下子步骤:
首先,为了对增广拉格朗日函数中的不同变量进行分别优化,采用交替方向法对增广拉格朗日函数中的各类变量进行迭代优化;
其次,对于核范数指数中的优化变量,采用混合奇异值截断算法来获得解析解。
本发明的方法具体运行的硬件和编程语言并不限制,用任何语言编写都可以完成,为此其它工作模式不再赘述。
本发明的实施例采用一台具有Intel Core-i3中央处理器和4G字节内存的计算机并用Matlab语言编制了基于指数化核范数与混合奇异值截断的张量恢复的工作程序,实现了本发明的方法。
本发明的基于指数化核范数与混合奇异值截断的张量恢复方法主要包括以下三个步骤:张量秩的定义以及该定义的凸松弛设计、辅助变量的引入及转化、交替方向优化策略等模块,具体步骤如下所述:
(1)提出一种有效的张量秩定义,并通过多种松弛策略将其转化为凸函数。其主要包含:
a)对于给定的张量χ,我们给出一种新的张量秩的定义,具体形式如下,
其中X(i)是张量χ在第i模态下的张开矩阵。
b)根据上述定义,张量恢复问题可以描述如下,
s.t.ρΩ(χ)=χΩ
其中ρΩ(·)表示提取张量下标为Ω的元素,这些元素是可以观测到的。
c)由于秩函数rank(·)是非凸的,采用秩函数的凸包络核范数||·||*来对其进行近似表示,则优化函数可转化为
d)根据不等式 ,可以进一步将上述优化问题近似为
s.t.ρΩ(χ)=χΩ
(2)通过引入一系列辅助矩阵变量消除张量不同模态下展开的矩阵的相关性,并采用拉格朗日乘子法将所有的约束条件转化为增广拉格朗日函数。
a)由于X(i)之间存在严重的相关性,使得对各变量的优化无法独立进行。为此,本发明引入一系列辅助矩阵变量
{Mi,i=1,2,…,m}来替换X(i),并增加相应的约束条件。因此,优化问题转化如下,
b)为了将上述约束优化问题转化为非约束优化问题,采用拉格朗日乘子法将所有的约束条件转化为增广拉格朗日函数
其中Qi表示拉格朗日矩阵乘子,<·,·>表示矩阵的内积,μi是损失参数。约束条件ρΩ(χ)=χΩ可以在求解过程中被直接处理,因此不需要加入到增广拉格朗日函数中。
(3)采用了一种交替方向法来实现对增广拉格朗日函数的分解迭代优化,其主要包含:
a)对于上述增广拉格朗日函数,采用交替方向法来对其中的变量{Mi,X(i),Qi}进行迭代优化;
b)对于由于{Mi,i=1,2,…,m}各自独立,因此对各个变量进行单独求解即可,求解过程如下,
对上述问题,本发明提出一种混合奇异值截断算法来获得解析解,即
其中 截断函数 UΣV*是矩阵Y的奇异值分解,tj*通过二分方法求解以下方程获得
其中,是矩阵Y的奇异值。
c)对于求解过程如下,
这是一个简单的最小二乘问题,结合约束ρΩ(χ)=χΩ可以获得如下解,
其中foldi(·)是将矩阵转化成张量的操作,是待恢复的张量元素。
d)对于求解过程如下,
e)步骤b)-d)不断迭代,直到收敛为止。
本发明的有益效果是:
1、本发明所提出一种全新的并且通用的张量秩的定义方法,对任意张量数据都适用。并且,该定义是张量CP秩的下界,能够有效地逼近CP秩,并消除了Tucker秩中权重参数的影响。
2、本发明提出一种有效的方法对张量秩定义进行松弛:将展开矩阵的秩松弛为展开矩阵的核范数,并且利用展开矩阵核范数的指数和的对数来逼近最大值函数,从而将张量恢复中的目标函数转化为凸函数。
3、本发明采用辅助变量来消除张量不同模态下展开的矩阵的相关性,并采用拉格朗日乘子法将所有的约束条件转化为增广拉格朗日函数,从而减小了优化问题的复杂度。
4、本发明采用了一种交替方向法来实现对增广拉格朗日函数的分解迭代优化。对于核范数的指数和中的优化变量,本发明提出了一种混合奇异值截断算法来获得解析解,从而实现快速有效的张量恢复。
上述实施例对本发明的具体描述,只用于对本发明进行进一步说明,不能理解为对本发明保护范围的限定,本领域的技术工程师根据上述发明的内容对本发明作出一些非本质的改进和调整均落入本发明的保护范围之内。
Claims (4)
1.一种基于指数化核范数与混合奇异值截断的张量恢复方法,其特征在于:包括以下三个步骤:
(1)提出一种新的张量秩定义:张量不同模态下展开矩阵秩的最大值;该定义是张量CP秩的下界,能够有效的逼近CP秩,并消除了权重参数的影响,采用核范数指数和的对数来逼近该张量秩定义,将其转化为凸函数;
(2)为了消除张量不同模态下展开的矩阵的相关性,引入一系列辅助变量来代替展开矩阵,并将约束条件利用拉格朗日乘子法转化为增广拉格朗日函数;
(3)采用交替方向法对增广拉格朗日函数中各类变量进行迭代优化,直到收敛;其中,对于核范数的指数和中的优化变量,采用混合奇异值截断算法来获得解析解。
2.根据权利要求1所述的一种基于指数化核范数与混合奇异值截断的张量恢复方法,其特征在于:所述的步骤(1)具体包括以下子步骤:
首先,根据张量CP秩和Tucker秩的优缺点,提出一种新的张量秩定义:张量展开矩阵秩的最大值;
其次,将展开矩阵的秩松弛为展开矩阵的核范数,并且利用核范数的指数和的对数来逼近最大值函数,从而将上述张量的秩定义转化为凸函数。
3.根据权利要求1所述的一种基于指数化核范数与混合奇异值截断的张量恢复方法,其特征在于:所述的步骤(2)具体包括以下子步骤:
首先,由于张量在不同模态下的展开矩阵具有相关性,引入一系列辅助矩阵变量来替换不同模态下的展开矩阵,并增加对应的约束条件;
其次,采用拉格朗日乘子法将所有约束条件加入到目标函数中,获得增广拉格朗日函数。
4.根据权利要求1所述的一种基于指数化核范数与混合奇异值截断的张量恢复方法,其特征在于:所述的步骤(3)具体包括以下子步骤:
首先,为了对增广拉格朗日函数中的不同变量进行分别优化,采用交替方向法对增广拉格朗日函数中的各类变量进行迭代优化;
其次,对于核范数指数中的优化变量,采用混合奇异值截断算法来获得解析解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410321348.0A CN104063852B (zh) | 2014-07-07 | 2014-07-07 | 一种基于指数化核范数与混合奇异值截断的张量恢复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410321348.0A CN104063852B (zh) | 2014-07-07 | 2014-07-07 | 一种基于指数化核范数与混合奇异值截断的张量恢复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104063852A true CN104063852A (zh) | 2014-09-24 |
CN104063852B CN104063852B (zh) | 2016-09-07 |
Family
ID=51551544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410321348.0A Active CN104063852B (zh) | 2014-07-07 | 2014-07-07 | 一种基于指数化核范数与混合奇异值截断的张量恢复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104063852B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104867119A (zh) * | 2015-05-21 | 2015-08-26 | 天津大学 | 基于低秩矩阵重建的结构性缺失图像填充方法 |
CN104966278A (zh) * | 2015-07-17 | 2015-10-07 | 电子科技大学 | 一种去除高光谱图像的噪声的方法及装置 |
CN107515843A (zh) * | 2017-09-04 | 2017-12-26 | 四川易诚智讯科技有限公司 | 基于张量近似的各向异性数据压缩方法 |
CN108107451A (zh) * | 2017-11-30 | 2018-06-01 | 西北工业大学 | 基于张量cp分解的矢量阵列干扰定位方法 |
CN108416750A (zh) * | 2018-03-08 | 2018-08-17 | 闽南师范大学 | 一种图像恢复方法 |
CN109658362A (zh) * | 2018-12-30 | 2019-04-19 | 东北大学秦皇岛分校 | 基于Capped核范数的三维张量补全方法 |
CN110428012A (zh) * | 2019-08-06 | 2019-11-08 | 深圳大学 | 脑网络模型建立方法、脑图像分类方法、装置及电子设备 |
CN110572389A (zh) * | 2019-09-06 | 2019-12-13 | 湖南大学 | 基于增量塔克分解的在线异常检测方法、系统及存储介质 |
CN111147863A (zh) * | 2019-12-30 | 2020-05-12 | 西安电子科技大学 | 基于张量的视频快照压缩成像恢复方法 |
CN116912107A (zh) * | 2023-06-13 | 2023-10-20 | 重庆市荣冠科技有限公司 | 一种基于dct的加权自适应张量数据补全方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005514683A (ja) * | 2001-12-06 | 2005-05-19 | ニューヨーク・ユニバーシティ | 合成、認識および圧縮のための多モード・データ集団の多重線形表現のための論理装置、データ構造、システムおよび方法 |
US8126229B2 (en) * | 2007-08-03 | 2012-02-28 | Siemens Medical Solutions Usa, Inc. | Reduction of lymph tissue false positives in pulmonary embolism detection |
CN101750029B (zh) * | 2008-12-10 | 2011-06-08 | 中国科学院沈阳自动化研究所 | 基于三焦点张量的特征点三维重建方法 |
CN101901472B (zh) * | 2010-07-07 | 2012-12-19 | 清华大学 | 一种基于矩阵秩最小化的非刚性鲁棒批量图像对齐方法 |
-
2014
- 2014-07-07 CN CN201410321348.0A patent/CN104063852B/zh active Active
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104867119B (zh) * | 2015-05-21 | 2018-08-24 | 天津大学 | 基于低秩矩阵重建的结构性缺失图像填充方法 |
CN104867119A (zh) * | 2015-05-21 | 2015-08-26 | 天津大学 | 基于低秩矩阵重建的结构性缺失图像填充方法 |
CN104966278A (zh) * | 2015-07-17 | 2015-10-07 | 电子科技大学 | 一种去除高光谱图像的噪声的方法及装置 |
CN104966278B (zh) * | 2015-07-17 | 2018-01-19 | 电子科技大学 | 一种去除高光谱图像的噪声的方法及装置 |
CN107515843A (zh) * | 2017-09-04 | 2017-12-26 | 四川易诚智讯科技有限公司 | 基于张量近似的各向异性数据压缩方法 |
CN108107451B (zh) * | 2017-11-30 | 2019-05-10 | 西北工业大学 | 基于张量cp分解的矢量阵列干扰定位方法 |
CN108107451A (zh) * | 2017-11-30 | 2018-06-01 | 西北工业大学 | 基于张量cp分解的矢量阵列干扰定位方法 |
CN108416750A (zh) * | 2018-03-08 | 2018-08-17 | 闽南师范大学 | 一种图像恢复方法 |
CN109658362A (zh) * | 2018-12-30 | 2019-04-19 | 东北大学秦皇岛分校 | 基于Capped核范数的三维张量补全方法 |
CN109658362B (zh) * | 2018-12-30 | 2023-09-05 | 东北大学秦皇岛分校 | 基于Capped核范数的数据恢复方法 |
CN110428012A (zh) * | 2019-08-06 | 2019-11-08 | 深圳大学 | 脑网络模型建立方法、脑图像分类方法、装置及电子设备 |
CN110572389A (zh) * | 2019-09-06 | 2019-12-13 | 湖南大学 | 基于增量塔克分解的在线异常检测方法、系统及存储介质 |
CN111147863A (zh) * | 2019-12-30 | 2020-05-12 | 西安电子科技大学 | 基于张量的视频快照压缩成像恢复方法 |
CN111147863B (zh) * | 2019-12-30 | 2021-06-08 | 西安电子科技大学 | 基于张量的视频快照压缩成像恢复方法 |
CN116912107A (zh) * | 2023-06-13 | 2023-10-20 | 重庆市荣冠科技有限公司 | 一种基于dct的加权自适应张量数据补全方法 |
CN116912107B (zh) * | 2023-06-13 | 2024-04-16 | 万基泰科工集团数字城市科技有限公司 | 一种基于dct的加权自适应张量数据补全方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104063852B (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104063852A (zh) | 一种基于指数化核范数与混合奇异值截断的张量恢复方法 | |
Shang et al. | Robust principal component analysis with missing data | |
Feigin et al. | Opers with irregular singularity and spectra of the shift of argument subalgebra | |
Yu et al. | Learning local dependence in ordered data | |
Zhang et al. | Correction of corrupted columns through fast robust Hankel matrix completion | |
CN103886235B (zh) | 一种正面人脸图像生物密钥生成方法 | |
Bibi et al. | High order tensor formulation for convolutional sparse coding | |
CN113746619A (zh) | 预定义时间同步控制的图像加密方法、解密方法及系统 | |
Kaidi et al. | New rank-2 Argyres-Douglas theory | |
Lian et al. | Periodic integrals and tautological systems | |
Hu | Sharp time estimates for exact boundary controllability of quasilinear hyperbolic systems | |
Dinuzzo et al. | Learning low-rank output kernels | |
Biggs et al. | Control affine systems on solvable three-dimensional Lie groups, II | |
Ji et al. | Flag structure for operators in the Cowen–Douglas class | |
CN105160154A (zh) | 一种基于平行因子的多维数据分析方法 | |
Liu et al. | Query-based attention CNN for text similarity map | |
Feng et al. | Bi-iterative algorithm for extracting independent components from array signals | |
Ćaćić et al. | Coupling of gravity to matter, spectral action and cosmic topology | |
Cleju et al. | Analysis-based sparse reconstruction with synthesis-based solvers | |
Kratica | An Electromagnetism-like method for the maximum set splitting problem | |
Shang et al. | Recovering low-rank and sparse matrices via robust bilateral factorization | |
Ali et al. | Deformed complex Hermite polynomials | |
Rodríguez | Alternating optimization low-rank expansion algorithm to estimate a linear combination of separable filters to approximate 2d filter banks | |
Novikov et al. | Reconstruction of vector (signal) by the norms of projections | |
Lu et al. | Fast compressive sensing of high-dimensional signals with tree-structure sparsity pattern |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |