普罗克汝斯忒斯的问题一般、三角和对称托普利兹矩阵

文摘

普罗克汝斯忒斯托普利兹问题是矩阵方程的最小二乘问题在托普利兹矩阵集。本文获得的必要和充分条件的解的存在性和唯一性普罗克汝斯忒斯托普利兹矩阵未知时约束问题的一般,三角形,分别和对称托普利兹矩阵。算法设计和数值例子表明,这些算法是可行的。

1。介绍

考虑约束最小二乘优化问题: 在哪里,,表示弗罗贝尼乌斯标准

这些问题已经被一系列的研究文献。例如,当,然后是一个解决方案(参见[1]),而对称矩阵解案例和正交矩阵的解决方案案例分析(2]和[3分别)。和对称的积极(半)的最小二乘问题讨论了海厄姆(4]。

本文以下托普利兹最小二乘问题分析: 在哪里和托普利兹矩阵的集合,一个矩阵如果其条目满足称为托普利兹矩阵,尽管,,与,,

在这项工作中,我们讨论问题(3)的细节。节2,一般托普利兹的问题讨论,在部分3、三角托普利兹问题和对称托普利兹问题进行了讨论,在部分4,我们给算法和数值例子。

我们首先给出一些定义和引理。

定义1(见[5])。矩阵被称为逆矩阵的当,它可以用或。

引理2(见[5])。让;如果存在非奇异的矩阵和这样 的充分必要条件是 在哪里,,。

备注3。当,,,用。

引理4(见[6])。一个有限维线性赋范线性空间的子空间包含至少一个点固定的最小距离点。

引理5(见[7])。让是一个凸函数连续可微的,然后任何局部最小值点也是一个全局最小值点。和是一个解决问题吗当且仅当满足,在那里意味着函数的梯度。

引理6(见[7])。二次函数凸当且仅当吗是半正定。

2。一般托普利兹问题

我们变换问题(3)下面的等价形式: 在哪里,被定义为 目标函数可以写成 右边的第一项(9)可以表示为 在哪里 右边的第二个任期(9)可以表示为 在哪里,,,。

因此收益率

采取偏导数关于给了

的一阶必要条件, 我们获得 因此下面的定理。

定理7。问题的解决方案(3)存在,它的一般形式可以表示为 在哪里,尽管意味着逆的。

证明。让,在那里,然后子空间弗罗贝尼乌斯规范下赋范线性空间。它遵循从引理4至少存在一个矩阵在这样。因此,通过的定义,存在至少一个标量构成矩阵,这意味着我们可以重写作为,在那里。因此问题的解决方案(3)存在。
很明显(13),目标函数是一个二次函数。此外,通过下面的讨论,也是一个凸函数。自是半正定,接下去,对吗,的表达在(10)持有 换句话说,是一种半正定矩阵。然后根据引理6,是一个凸函数。自是连续可微的,然后我们得到任何解决方案(16)也是一个解决方案(7由引理5。显然,线性系统的解决方案(16)是 (见[5])。所需的解决方案的问题(7然后给出了),在那里。定理是等价的问题(3)和问题(7)。

定理8。问题(3当且仅当)有一个独特的解决方案满秩。在这种情况下,唯一的解决方案是 在哪里。

证明(必要性)。如果已满秩,线性系统(16)有一个独特的解决方案;因此,解决方案的问题(7)是唯一决定的,所以是问题(3)。
相反过程的充分性是必要的。

推论9。如果在列满秩,那么问题的解决方案(3)是独一无二的。

3所示。三角托普利兹问题和对称托普利兹的问题

在本节中,我们讨论了三角形托普利兹问题和对称托普利兹的问题。

3.1。三角托普利兹的问题

上三角托普利兹的问题,因为[的定义6),我们的意思是最小化问题 在哪里,是上三角托普利兹矩阵的子空间,一般形式的元素

我们变换问题(21)以下问题 在哪里 证明过程是类似于问题(7)。在这种情况下,相关的函数,未知,线性方程组,矩阵以及标量可以表示为

因此,我们得到一个结论上三角托普利兹问题如下。

定理10。问题的解决方案(21)存在,它的一般形式可以表示为 在哪里,尽管,意味着逆的。

定理11。问题(21当且仅当)有一个独特的解决方案满秩。在这种情况下,唯一的解决方案是 在哪里。

推论12。如果在列满秩,那么问题的解决方案(21)是独一无二的。

同样,我们可以解决下三角托普利兹的问题。

3.2。对称托普利兹的问题

对称托普利兹的问题是以下最小化问题: 在哪里,是对称的托普利兹矩阵的子空间,元素的一般形式: 我们考虑下面的目标函数: 在哪里

在这种情况下,相关的函数,未知,线性方程组,矩阵以及标量可以表示为

因此,我们获得一个结论关于对称托普利兹的问题。

定理13。问题的解决方案(28)存在,它的一般形式可以表示为 在哪里,尽管,意味着逆的。

定理14。问题(28当且仅当)有一个独特的解决方案满秩。在这种情况下,唯一的解决方案是 在哪里。

推论15。如果在列满秩,那么问题的解决方案(28)是独一无二的。

4所示。计算和例子

通用托普利兹的推导问题导致下面的计算算法。

算法16。(1)生成矩阵根据定义(7)。(2)计算矩阵和标量由公式(11)和(12)。(3)如果是满秩的,那么解决线性系统(16)直接获取解决方案;其他计算,形成,尽管。一般来说,让,和在引理2等于,我们可以获得的逆的,这是。(4)计算从,在那里。

其他特殊托普利兹的算法问题可以设计类似。

我们提到,在本文的开始,这两个实例的解决方案实际上是基于算法获得16。我们可以给托普利兹矩阵问题的另一个例子列满秩的如下。根据推论9,12,15,解决方案是独一无二的。

示例17。考虑以下
利用上述算法,通过Matlab计算托普利兹问题,我们得到以下结果。
一般托普利兹矩阵解决方案 上三角托普利兹矩阵解决方案 和对称托普利兹矩阵解决方案

5。结论

我们已经讨论了普罗克汝斯忒斯托普利兹的问题包括,三角形,对称的情况下。改造后的起源问题转化为二次曲面形态,我们与工厂获得一般的或独特的解决方案,通过求解线性方程组。

确认

作者要感谢m . g . Eberle和m . c .斑点首次提出和研究了普罗克汝斯忒斯托普利兹的问题。他们也感谢裁判任何有益的建议或加强定理。

引用

1. g·h·戈卢布和c f . van贷款,矩阵计算,3卷,约翰·霍普金斯大学出版社,巴尔的摩,医学博士,美国,1983年。视图:MathSciNet
2. n·j·海厄姆,”普罗克汝斯忒斯的对称问题,”一些。数值数学,28卷,不。1,第143 - 133页,1988。视图:出版商的网站|谷歌学术搜索|Zentralblatt数学|MathSciNet
3. n·j·海厄姆,“极地decomposition-with应用,计算”工业与应用数学学会的。《科学与统计计算,7卷,不。4、1160 - 1174年,1986页。视图:出版商的网站|谷歌学术搜索|Zentralblatt数学|MathSciNet
4. n·j·海厄姆,”一个最近的对称半正定矩阵,计算”线性代数及其应用卷,103年,第118 - 103页,1988年。视图:出版商的网站|谷歌学术搜索|Zentralblatt数学|MathSciNet
5. x c .他和y . w .太阳,Guangyi倪Juzhen Yinlun,江苏科学技术出版社,南京,中国,1991。视图:MathSciNet
6. e·w·切尼介绍了近似理论,麦格劳-希尔书,纽约,纽约,美国,1966年。视图:MathSciNet
7. l . d .回族和t . x乔数值优化,科学出版社,北京,2006。

版权

版权©2013杨胡安和Yuan-bei邓。这是一个开放的分布式下文章知识共享归属许可,它允许无限制的使用、分配和复制在任何媒介,提供最初的工作是正确引用。