广义特征值与极大极小原理
- 格式:doc
- 大小:198.00 KB
- 文档页数:4
广义特征值分解广义特征值分解(Generalized Eigenvalue Decomposition, GED)是一种重要的矩阵分解方法,常常被应用在信号处理、机器学习等领域中。
它能够将两个矩阵同时对角化,得到它们的特征向量和特征值。
在本文中,我们将对广义特征值分解做一个详细的讲解。
步骤一:理解特征值与特征向量在矩阵计算中,特征向量是指在矩阵进行线性变换后仍然保留其方向的向量。
特征值是与特征向量相关的标量,描述了该特征向量在变换中的“伸缩”程度。
一般来说,我们可以通过解决以下方程式来找到一个矩阵的特征向量和特征值:(A−λI)v=0其中,A是一个方阵,λ是一个标量,I是单位矩阵,v是特征向量。
步骤二:理解广义特征值分解在广义特征值分解中,我们要找到两个矩阵A和B的特征向量和特征值。
也就是说,我们需要解决以下方程式:Av=λBv其中,v是特征向量,λ是特征值。
将其转化为标准形式:(A−λB)v=0这样,我们就可以将两个矩阵同时对角化,得到它们的特征向量和特征值。
步骤三:寻找广义特征值分解在实际应用中,我们可以使用数值计算方法来寻找广义特征值分解。
这包括基于迭代算法的方法,如幂法、反幂法和雅可比迭代法等。
其中,幂法是最常用的方法之一,可以用来寻找矩阵的最大特征值和对应的特征向量。
雅可比迭代法则是另一个最常见的方法,可以用来寻找所有特征值和对应的特征向量。
步骤四:应用广义特征值分解广义特征值分解在实际应用中有很多用途。
例如,它可以用来处理分析较大的数据集、图像、信号等。
在信号处理领域中,可以将电磁波分解为多个频率的成分,在图像处理领域中,可以寻找图像的相似性和模式。
在机器学习领域中,广义特征值分解可以被用来进行降维和特征提取。
总之,广义特征值分解是一种非常有用的矩阵分解方法,在很多领域中都被广泛应用。
理解广义特征值分解的原理和寻找方法,将有助于我们更好地应用这种方法来解决实际问题。
有限元结构动力分析的广义特征值的神经计算
李海滨;黄洪钟;赵明扬
【期刊名称】《哈尔滨工业大学学报》
【年(卷),期】2006(38)9
【摘要】广义特征值问题是结构动力分析计算的呼和浩特 010062;3.大连理工大学机械工程学院,辽宁大连 116023)关键之一.应用Reyleigh极小值原理,将神经网络的能量函数的极小点对应于广义特征值问题的最小特征值所对应的特征向量,在神经网络朝着能量函数极小点运动的同时得到了最小特征值所对应的特征向量的精确解答.从特征值的变分特性出发,给出了基于罚函数法的其他特征值的神经网络求解方案,从而在理论上给出了广义特征值问题的所有特征值的神经网络求解方法.仿真计算表明,该方法正确、有效可行.
【总页数】5页(P1523-1527)
【作者】李海滨;黄洪钟;赵明扬
【作者单位】中国科学院,沈阳自动化研究所,沈阳,110015;内蒙古工业大学理学院,呼和浩特,010062;大连理工大学机械工程学院,辽宁,大连,116023;中国科学院,沈阳自动化研究所,沈阳,110015
【正文语种】中文
【中图分类】O3
【相关文献】
1.广义特征值问题与标准特征值问题 [J], 段辉明;李玲;李永红
2.由部分特征值和顺序主子阵构造广义Jacobi矩阵的逆特征值问题 [J], 潘云兰;秦立
3.由特征值和顺序主子阵构造广义Jacobi矩阵的逆特征值问题 [J], 徐秀斌;秦立
4.广义自反矩阵与广义反自反矩阵的广义逆特征值问题 [J], 刘能东;张忠志
5.大型广义特征值问题的部分特征值和特征向量的块迭代求解 [J], 赵小红;陈飞武;吴健;周巧龙
因版权原因,仅展示原文概要,查看原文内容请购买。
极大极小算法原理极大极小算法原理是一种优化算法,广泛应用于数学、计算机科学和工程领域。
它的核心思想是在搜索空间中寻找一个最优解,这个最优解既要满足目标函数的最大值或最小值,又要满足约束条件的可行性。
在这个过程中,极大极小算法通过不断地迭代和更新解的方式,逐步逼近最优解。
一、算法基本原理1.极大极小算法的出发点是:在满足约束条件的前提下,寻找一个使目标函数取得最大或最小值的解。
2.算法的基本步骤如下:(1)初始化:给定问题的参数,如目标函数、约束条件、初始解等。
(2)迭代:在搜索空间中,根据当前解更新解,得到一个新的解。
更新策略可以是基于目标函数的梯度信息,也可以是基于启发式信息等。
(3)评估:评估新的解是否满足停止条件,如目标函数值的改进、约束条件的满意度等。
(4)停止:如果满足停止条件,算法收敛,输出当前解作为最优解;否则,返回步骤2,继续迭代。
二、算法的性质和优点1.极大极小算法具有全局收敛性,即在一定条件下,算法一定能找到最优解。
2.算法具有较强的适应性,可以处理非线性、非凸优化问题,以及具有复杂约束条件的问题。
3.极大极小算法在搜索过程中,可以充分利用目标函数的梯度信息,加快收敛速度。
4.算法具有较好的鲁棒性,对初始解的选择不敏感,适用于不同领域的问题。
三、算法应用与发展1.应用领域:极大极小算法在工程、经济、生物信息学等领域具有广泛的应用,如参数优化、供应链管理、基因表达数据分析等。
2.算法改进:针对原始极大极小算法的不足,研究者们提出了许多改进算法,如自适应步长调整、混合策略、并行计算等。
3.与其他算法比较:极大极小算法与其他优化算法(如梯度下降、牛顿法、遗传算法等)相比,具有更高的收敛速度和更好的全局性能。
总之,极大极小算法作为一种有效的优化方法,在理论上具有严谨的收敛性保证,同时在实际应用中表现出良好的性能。
随着科学技术的不断发展,极大极小算法在未来将继续发挥重要作用,为各个领域的优化问题提供有力支持。
第一章绪论在科学研究和工程应用中,经常需要求解大型稀疏线性方程组血=b(1.1)其中A是n×n的实矩阵,x,6∈R”.目前,求解线性方程组的数值方法可分成两大类,一类是直接法,即通过有限次的运算求出问题的精确解,例如Gauss消去法、列主元及全主元消去法、直接三角分解法等.但是由于直接法计算过程中存储量很大,当需要求解大型稀疏线性方程组时,直接法就不适用了.另一类求解线性方程组的数值方法是迭代法,即通过选取初值,然后用同样的步骤重复计算,求得近似解.在迭代法中,Krylov子空间方法[31是求解大型线性方程组的一类重要方法,国际上有关Krylov子空间方法的研究工作非常活跃.求解对称正定线性方程组的最有效方法是共轭梯度(cO)法【l】及其预处理技术.对称Lanczos方法【13】【161伫9】是解对称不定线性方程组的有效方法之一.在理论上,对称Lanczos方法产生的向量组是正交向量组,但是,在实际计算中,由于舍入误差的影响,Lanczos向量易失去正交性.为了减少存储量和运算量,人们采用重新开始的Lanczos方法,即循环Lanczos迭代法【11I.另外Paigc和Saunders基于对称LRI'ICZOS方法[2]提出了求解对称不定方程组的SYⅣnVJLQ方法【14】和MINI陋S方法【l”,但是对病态线性方程组SYMMLQ方法和MINRES方法常常表现出不稳定性.求解非对称线性方程组的Krylov子空间方法有许多特殊的方法,如Amoldi口]方法、广义极小残量法(GMRES)、双边Lanczos方法、不完全正交化方法等.Y.Saadl31指出,Amoldi【21过程实际上是建立Krylov子空间k,(A,to)=span(to,Aro,...,A”1to)一组标准正交基的过程.将Amoldi回过程用于求解线性方程组可得完全正交化方法(FOM)rss],不完全正交化方法(IOM)H】【3】【loJ是完全正交化方法的一个变形,在理论上它是一种斜投影方法.1991年Freund和Nachfigal基于非对称Lanczos方法提出了求解非对称线性方程组的拟极小残量法(QMR方法)[8】.在用非对称Lanczos方法解非对称线性方程组时,也会发生算法中断或数值不稳定.Krylov子空间方法通常用残量范数作为判断算法终止的条件.若近似解是精确的,残量范数是小的,但是反过来不一定.为克服残量范数作为终止条件的不足,Kasenally[”】在用GMRES方法15删解非对称线性方程组时,考虑求满足扰动方程(A一△。
线性代数中的广义特征值问题线性代数是数学的一个分支,主要研究向量、向量空间和线性变换等概念。
在线性代数的学习中,广义特征值问题是一个重要的概念。
本文将详细介绍广义特征值的概念、求解方法以及其在实际问题中的应用。
一、广义特征值的概念在传统的特征值问题中,我们考虑的是一个方阵A的特征值和特征向量。
特征值是一个标量,而特征向量是方阵A乘以该向量等于特征值乘以该向量。
然而,在某些情况下,方阵A可能是非方阵,这时候就需要考虑广义特征值问题。
广义特征值是非方阵的特征值。
设矩阵A为m×n维,特征向量x为n维列向量,特征值λ为标量。
则广义特征值问题可以表示为Ax = λBx,其中B为m维方阵。
为了求解该问题,需要考虑B的非奇异性。
二、广义特征值问题的求解方法解决广义特征值问题的方法有很多,下面介绍几种常用的方法。
1. 通用特征值问题的转化:将广义特征值问题转化为标准特征值问题。
这种方法适用于一些特殊情况,例如B是正定的或者B很接近一个正定矩阵时。
通过对矩阵进行相似变换,可以将广义特征值问题转化为标准特征值问题,从而利用已有的求解方法求解。
2. 修正特征值问题的求解:对于一些特殊的B矩阵,例如对称正定的B矩阵,可以利用修正特征值问题进行求解。
通过将广义特征值问题转化为修正特征值问题,进一步求解得到广义特征值。
3. 广义特征值问题的迭代法:迭代法是一种常用的数值求解方法,对于广义特征值问题也有相应的迭代算法。
例如广义幂法,可以通过迭代的方式逐渐逼近广义特征值问题的解。
三、广义特征值问题的应用广义特征值问题在实际问题中具有广泛的应用。
以下列举一些常见的应用领域。
1. 物理学中的应用:广义特征值问题在量子力学中有很多的应用。
例如,通过广义特征值问题可以求解量子力学中的能量本征值和波函数等。
2. 工程中的应用:广义特征值问题在结构动力学和振动工程中有着重要的应用。
通过求解广义特征值问题,可以得到结构物的固有频率和振型等信息,从而评估结构物的稳定性和安全性。
极大极小原理
极大极小原理,也被称为最大最小原理或极值原理,是数学分析中经常使用的重要原理。
极大极小原理的核心思想是,如果一个函数在某个区间内取得了最大值或最小值,那么在该区间的边界上也必须取得最大值或最小值。
这个原理对于研究函数的性质和求解最优化问题非常有帮助。
举个例子来说明极大极小原理的应用。
考虑一个函数f(x),在
闭区间[a,b]上定义。
如果在区间内存在一个点x0,使得f(x0)
是f(x)在[a,b]上的最大值,那么可以得出结论:f(x)在区间的
边界上的值必须小于或等于f(x0),即f(x)<=f(x0),其中x属
于[a,b]。
同样地,如果f(x)在[a,b]上取得了最小值f(x0),则
f(x)>=f(x0)。
极大极小原理的应用非常广泛,可以用于证明函数的性质,如连续性、有界性等。
另外,它还可以用于解决一些最优化问题,比如求解函数的最大值或最小值的问题。
基于极大极小原理,可以通过在区间的边界上寻找函数的最大值或最小值,从而找到整个区间上的最大值或最小值。
总之,极大极小原理是数学分析中一项重要的基本原理,通过利用函数在区间内最大最小值的性质,可以帮助我们研究函数的性质和解决最优化问题。
斜拉桥稳定性整体分析【摘要】本文对斜拉桥稳定理论的研究发展概况进行了总结,详细地论述了斜拉桥失稳的两类稳定性问题,并对其稳定问题失稳判别准则进行了分析,探讨了斜拉桥稳定性的两种评价指标。
【关键词】斜拉桥,稳定理论,失稳判别准则,评价指标结构失稳是指在外力作用下结构的平衡状态开始丧失稳定性,稍有扰动,也会引起很大的位移和变形,甚至发生破坏。
此时虽然截面的内力并未超过它的最大抵抗能力,但结构的平衡状态发生了分支,或者是随着变形的发展内外力的平衡己不可能得到,于是结构在外荷载基本不变的情况下可能发生很大的位移最后导致结构的破坏。
一、稳定理论的发展概况与桥梁结构相关的稳定理论已有悠久的历史,同时桥梁失稳事故的发生促进了桥梁稳定理论的发展。
早在1744年欧拉(L.Euler)就进行了弹性压杆屈曲的理论计算。
在国内对于斜拉桥的稳定性问题,李国豪等提出了采用空间杆系屈曲有限元方法进行计算的思路,并给出了计算斜拉桥平面屈曲临界荷载的近似方法。
根据结构经受任意微小外界干扰后,能否恢复初始平衡状态,可把平衡状态分为稳定、不稳定和随遇三种。
研究结构稳定的主要目的就在于防止不稳定平衡状态的发生。
由失稳前后平衡和变形性质,可以把稳定问题分为两大类:第一类稳定,即分支点失稳问题。
见图1;第二类稳定,即极值点失稳问题,见图2。
图1 分支点失稳图2极值点失稳二、斜拉桥的第一类稳定问题在斜拉桥建设的初期,跨径一般较小,再加上计算手段的不成熟,通常只考虑第一类稳定问题,而且常把塔和梁分离开来单独考虑其稳定性。
对斜拉桥稳定性较精确的分析方法是有限元法,这种方法可求得斜拉桥整体的屈曲安全度。
在有限元分析中,斜拉桥被离散为许多单元。
如果知道各个单元的力和位移的关系,则不难推出整体结构的力和位移的关系。
值得注意的是,在压杆刚度矩阵中,需要考虑轴向力对刚度的影响。
对于第一类稳定问题而言,结构失稳时是处于小变形范围,大位移矩阵[KL]较小,通常忽略不计。
第六章 动力问题的有限元法6.1 概述前面几章所研究的问题都属于静力问题,其特点是施加到结构上的外载荷不会使结构产生加速度,且外载荷的大小和方向不随时间变化,因而结构所产生的位移和应力也不随时间变化。
本章将要研究结构分析中另一类重要问题的有限元解法,即动力问题的有限元解法。
动力学问题的特点是,载荷是随时间变化的,因而结构所产生的位移和应力是时间的函数,结构会产生速度和加速度。
由于结构本身的弹性和惯性,结构在动力载荷的作用下,往往呈现出振动的运动形态。
结构振动是工程中一个很普遍很重要的问题。
有些振动对我们有利,例如,振动打桩,振动选料,有些振动对我们有害,例如,机床的振动,仪器与仪表的振动,桥梁、水坝及高层建筑在地震作用下的振动等。
因此,我们必须对振动体本身的振动特性以及它对外部激振力的响应有一个明确的认识,才能更好地利用它有利的一面,而避免它有害的一面,设计出更好的机械和结构。
振动问题主要解决两方面的问题。
1. 寻求结构的固有频率和主振型,从而了解结构的固有振动特性,以便更好地利用或减少振动。
2. 分析结构的动力响应特性,以计算结构振动时动应力和动位移的大小及其变化规律。
6.2 结构的振动方程结构的振动方程可用多种方法建立,这里我们使用达朗伯原理(动静法),仿照前几章建立静力有限元方程的方法,来建立动力问题的有限元方程。
在静力问题中用有限元法建立的平衡方程是}{}]{[F K =δ在振动问题中,对结构的各节点应用达郎伯原理所建立的振动方程仍然具有与上式相同的形式,只不过节点位移是动位移,节点载荷是动载荷,它们都是时间的函数。
上面的方程成为)}({)}(]{[t Q t K =δ (6.1)上式中{})(t δ为节点的动位移,它是时间的函数,)}(]{[t K δ是t 时刻的节点位移产生的弹性恢复力,它与该时刻的节点外力{})(t Q 构成动态平衡。
在动态情况下,结构承受的载荷(集中载荷 ,分布载荷 )可随时间而变化,是时间的函数。
极大极小原理极大极小原理是数学中一个非常重要的概念,它在优化问题中有着广泛的应用。
极大极小原理是指对于一个函数而言,如果在某个点取得极大值或者极小值,那么这个点的导数为零。
极大极小原理是微积分中的基本定理之一,它可以帮助我们求解各种优化问题,包括最大值、最小值以及最优化问题。
首先,我们来看一下极大极小值点的定义。
对于一个函数f(x),如果在点x=a处,f'(a)=0,并且f''(a)>0,那么这个点就是函数f(x)的极小值点;如果在点x=b处,f'(b)=0,并且f''(b)<0,那么这个点就是函数f(x)的极大值点。
通过这个定义,我们可以很容易地找到函数的极大极小值点,从而解决各种优化问题。
在实际应用中,极大极小原理可以帮助我们解决各种实际问题。
比如在工程中,我们经常需要求解一些最优化问题,比如最小成本、最大收益等。
而这些问题往往可以转化为一个函数的极大极小值问题,通过极大极小原理,我们可以很快地找到最优解。
另外,在经济学、物理学、生物学等领域,极大极小原理也有着广泛的应用,可以帮助我们解决各种实际问题。
除了在求解最优化问题中的应用,极大极小原理在数学分析中也有着重要的地位。
通过极大极小原理,我们可以判断一个函数在某个区间内的极值点,从而对函数的性质有更深入的了解。
同时,极大极小原理也为我们提供了一种判断函数极值点的方法,通过这种方法,我们可以更快地找到函数的极值点,从而简化了计算过程。
总之,极大极小原理是数学中一个非常重要的概念,它在优化问题中有着广泛的应用。
通过极大极小原理,我们可以解决各种最优化问题,同时也可以对函数的性质有更深入的了解。
极大极小原理不仅在数学分析中有着重要的地位,同时也在工程、经济学、物理学、生物学等领域有着广泛的应用,可以帮助我们解决各种实际问题。
因此,对极大极小原理的深入理解和掌握,对于我们来说是非常重要的。
结构动态特性分析结构固有特性分析在数学上称之为特征值和特征向量分析,包含固有频率与固有模态分析,是结构动力学中的主要任务之一。
结构固有特性分析是为了研究结构振动的固有规律和内在本质,为结构动力学的进一步分析打下基础,在工程的实际应用以及在求解结构动力响应方面具有很重要的意义。
到目前为止,已经发展了许多求解动态特征问题的数值方法。
在通常的特征值求解方法中,根据解法的特点,可分为四个基本类型:多项式迭代技术、应用特征多项式的Sturm序列的分解法、矢量迭代法和变换方法。
这里不对这些方法做一一介绍,只介绍一些典型常用方法的特点、理论依据以及它们的应用。
、特征值问题的性质结构无阻尼自由振动方程为(6-18)设结构作简谐运动,即(6-19)式中:ω为圆频率,θ为相位角,φ为振幅。
将上式代人式(6-18)得:(6-20)或写成(6-21)其中,;K,M分别为结构的刚度矩阵和质量矩阵。
式(6-21)是结构动力学的广义特征值问题。
显然,由式(6-21)求出的和的值,只取决于结构本身的刚度矩阵K和质量矩阵M,即它们是结构的固有值。
就是结构自振圆频率,称为结构的特征值,与ω相应的空间振动形态(即振型或模态)称为特征向量。
式(6-21)反映的是结构的动态特性,我们的任务就是求解λ和。
在研究特征值问题的具体算法之前,先讨论特征对的一些基本特性。
特征对有如下的特性:(1)如果K和M都对称,且至少有一个矩阵正定,则特征值一定是实数,特征向量也一定是实向量。
如果M正定,并且K为正定或半正定,则所有特征值都是正的实数。
(2)特征向量(或模态向量)关于质量矩阵M和刚度矩阵K正交,即:(6-22)(6-23)在式(6-21)中将特征向量归一化,即:(6-24)式(6-24)称为归一化特征向量。
则式(6-22),(6-23)有:(6-25)(6-26)(3)Ralyeigh商和特征值的极大极小性质定义:(6-27)称为Ralyeigh商。
矩阵论中的广义特征值问题研究矩阵论作为现代代数学的一个重要分支,研究的是关于矩阵的性质和矩阵运算的规律。
其中,广义特征值问题是矩阵论中的一个重要议题。
本文将探讨广义特征值问题的定义、求解方法以及其在实际应用中的作用。
一、广义特征值问题的定义在矩阵论中,一般的特征值问题是指求解方阵A的特征值和特征向量。
而广义特征值问题则是在一般的特征值问题基础上引入了矩阵B,求解的是广义特征值和广义特征向量。
具体而言,给定两个n阶矩阵A和B,广义特征值问题可以表示为:Ax = λBx其中,x为非零向量,λ为广义特征值。
如果存在一个非零向量x和对应的广义特征值λ,满足上述等式,则称x为A和B的广义特征向量和广义特征值。
二、求解广义特征值问题的方法为了求解广义特征值问题,我们需要利用矩阵分解和特征值分析的方法。
以下介绍两种常用的求解方法。
1. 广义特征值问题的化简方法通常情况下,我们可以通过矩阵的相似变换将广义特征值问题转化为一般的特征值问题。
具体步骤如下:步骤一:对矩阵B进行特征值分解,得到矩阵B的特征值和特征向量。
步骤二:将B特征向量作为基,对矩阵A进行相似变换。
步骤三:转化后的矩阵A即为一般的特征值问题,进行求解即可得到广义特征值和广义特征向量。
2. 广义特征值问题的广义奇异值分解方法广义奇异值分解是求解广义特征值问题的一种有效方法,适用于非对称矩阵或奇异矩阵的情况。
具体步骤如下:步骤一:对矩阵A和B同时进行广义奇异值分解。
步骤二:根据广义奇异值分解的结果,得到广义特征值和广义特征向量。
根据以上求解方法,我们可以有效地解决广义特征值问题,从而获得矩阵A和B的广义特征值和广义特征向量。
三、广义特征值问题的应用广义特征值问题在实际应用中有着广泛的应用。
下面介绍一些典型的应用场景。
1. 工程结构振动分析广义特征值问题可用于工程结构的振动分析。
通过求解广义特征值问题,可以确定结构的固有频率和振型,从而评估结构的稳定性和抗震性能。
第二十一讲 广义特征值与极小极大原理
一、 广义特征值问题
1、定义:设A 、B 为n 阶方阵,若存在数λ,使得方程Ax Bx =λ存在非零解,则称λ为A 相对于B 的广义特征值,x 为A 相对于B 的属于广义特征值λ的特征向量。
● 是标准特征值问题的推广,当B =I (单位矩阵)时,广义特征值问题退化为标准特征值问题。
● 特征向量是非零的 ● 广义特征值的求解
()A B x 0-λ= 或者 ()B
A x 0λ-=
→
特征方程 ()det A B 0-λ=
求得λ后代回原方程Ax Bx =λ可求出x
本课程进一步考虑A 、B 厄米且为正定矩阵的情况。
2、等价表述
(1) B 正定,1B -存在
→1
B A x x
-=λ,广义特征值问题化为了标准
特征值问题,但一般来说,1B A -一般不再是厄米矩阵。
(2) B 正定,存在Cholesky 分解,H B G G =,G 满秩 H A x G G x =λ 令H G x y = 则 ()
1
1
H
G A G
y y --=λ 也成为标准特征值问题。
(
)
1
1
H
G A G
--为厄米矩阵,广义特征值是实数,可以按大小顺序
排列12n λ≤λ≤≤λ ,一定存在一组正交归一的特征向量,即存在
12n y ,y ,y 满足
()
1
1
H
i i G A G
y y --=λ
H i
j ij 1i j y y 0
i j
=⎧=δ=⎨≠⎩
还原为()1
H i i x G y -= (i=1,2, ,n),则 ()()
H H H
H i
j i
j i
j ij 1
i j y y x G G x x Bx 0
i j
=⎧===δ=⎨
≠⎩ (带权正交)
二、 瑞利商
A 、
B 为n 阶厄米矩阵,且B 正定,称()()H
H
x A x R x x 0x Bx
=≠为A
相对于B 的瑞利商。
12n x ,x ,x 线性无关,所以,n
x C
∀∈,存在12n a ,a ,a C ∈ ,使
得 n
i
i
i 1
x a x ==
∑
H
n
n n
n
2
H
H
i i i j j j i
j i
i 1j 1i ,j 1
i 1
x Bx a x B a x a a x Bx a ====⎛⎫⎛⎫
==
=
⎪ ⎪⎝⎭⎝⎭
∑∑∑
∑
n
n
n
2
H
H H i i j i
j j i
i j i i
i ,j 1
i ,j 1
i 1
x A x a a x A x a a x Bx a ====
=
λ=
λ∑
∑
∑
∴ ()n
2
i i i 1n
2
i
i 1
a R x a ==λ=
∑
∑
●()1x 0
min R x ≠=λ ()n x 0
max R x ≠=λ
证明:()()()()
()
H
H
H
H
kx A kx x A x R x x Bx
kx B kx =
=
k 为非零常数
可取1k x
=,
kx 1=
∴ ()H
H
x 1
x A x R x x B x
==
(闭区域)
当1x
x =或()i a 0i 2,3,,n == 时,()1R x =λ
i 1λ≥λ ()n
2
i i 11
1n
2
i
i 1
a R x a ==≥λ=λ∑
∑
∴
()1x 0
min R x ≠=λ
另一方面,i n λ≤λ ()n
2
i i 1n
n n
2
i
i 1
a R x a ==≤λ=λ∑
∑
∴ ()n x 0
max R x ≠=λ
[证毕] 当B =I 时,标准特征值问题 A x x =λ (H A A =)
12n
H
i j ij
x x λ≤λ≤≤λ⎧⎨=δ⎩
则 ()
H
1H
x 0x A x m in
x x
≠=λ
()
H
n H
x 0x Ax max
x x
≠=λ
进一步分析可得
()
12x 0
a 0
m in R x ≠==λ ()
n n 1x 0
a 0
m ax R x -≠==λ
()
12k k 1x 0
a a a 0
m in R x +≠=====λ ()
n n 1n k n k 1x 0
a a a 0
m ax R x ----≠=====λ
定理1.设{}r r 1s L span x ,x ,,x += ()r r 1s +λ≤λ≤≤λ ,则 ()r x 0x L
m in R x ≠∈=λ ()s x 0x L
m ax R x ≠∈=λ
这一结果不便于应用,希望对上述结果进行改造,改造成不依赖于i x 的一种表达方式。
1a 0=和n a 0=的情况均对应于
x 在(n-1)维的子空间内变动,
x 在L 中变动是在一个(s-r+1)维子空间中变化。
一般的,x 在n C 的(n-1)维子空间n 1V -中变动时,
()n 1
2x 0
x V m in R x -≠∈≤λ ()n 1
n 1x 0x V m ax R x --≠∈≥λ
即,对于不同的n 1V -,()R x 的最小值及最大值有可能不同,其中各个最小值中最大者为2λ,各个最大值中的最小者为n 1-λ
()n n 1n 12x 0V C x V m ax m in R x --≠∈∈⎡⎤
=λ⎢⎥⎢⎥⎣⎦
()n n 1n 1n 1x 0V C x V m in m ax R x ---≠∈∈⎡⎤
=λ⎢⎥⎢⎥⎣⎦
定理2. 设k V 是n C 的一个k 维子空间,则
()n
k k n k 1x 0V C x V m ax m in R x -+≠∈∈⎡⎤
=λ⎢⎥⎢⎥⎣⎦
()n k k k x 0V C x V m in m ax R x ≠∈∈⎡⎤
=λ⎢⎥⎢⎥⎣⎦
以上两式称为广义特征值的极小极大原理。
● B =I 时,标准特征值问题同样存在上述关系。
● 矩阵奇异值问题:()()2
H
A A A ⎡⎤σ=λ⎣⎦
(非零) ()()
H
H
2
H
2
x
A
A x
A x R x x x
x
=
=
n
k k
2
n k 1
x 0V C x V 2A x m ax m in x -+≠∈∈⎡⎤σ=⎢⎥⎢⎥⎣⎦
n k k
2
k x 0V C x V 2A x m in m ax
x ≠∈∈⎡⎤
σ=⎢⎥⎢⎥⎣⎦。