网站建设网络推广图片,aspx网站配置服务器,石龙镇做网站,微信朋友圈推广方案简介
流形优化是非线性优化的一个分支#xff0c;它主要关注在特定的几何结构下进行优化。在流形优化中#xff0c;优化问题通常是在黎曼流形上进行的#xff0c;而非欧几里得空间。黎曼流形是带有黎曼度量的流形#xff0c;该度量为流形上的每个点都定义了一个内积。这种…
简介
流形优化是非线性优化的一个分支它主要关注在特定的几何结构下进行优化。在流形优化中优化问题通常是在黎曼流形上进行的而非欧几里得空间。黎曼流形是带有黎曼度量的流形该度量为流形上的每个点都定义了一个内积。这种内积结构提供了流形上测量长度和角度的方式这在优化过程中非常重要因为它允许我们定义梯度和Hessian等概念并进行相应的优化操作。在流形优化的背景下流形通常是解的约束集。例如当解必须是正交矩阵、单位范数向量或其它特殊结构时这些集合都可以被看作流形。对于连续优化讨论多面体是考虑线性优化问题的基础而在欧几里得空间中的微积分被用于发展非线性欧几里得优化理论。同样地黎曼流形上的微积分借助微分几何学尤其是黎曼几何学得以促进对于黎曼流形上的优化是至关重要的。
流形优化相对于传统欧几里得空间上的优化方法优势在于
在黎曼流形上迭代减少建模误差提高一些非线性问题的表示能力;充分利用目标函数的内在几何结构通过选取合适的黎曼度量可将一些非凸优化问题转化为凸优化问题;保持最优化问题的尺度不变性。
在各个应用领域的需求牵引下相关学者的研究工作极大地促进了黎曼流形优化理论的发展与工程应用如Hermitian 特征值问题、自适应正则化、迹范数、相机位姿估计、经济负载分配、对称特征值问题、低秩矩阵填充等其技术成果已成功应用于相关学科如航空航天、理论物理、无线通信、医学成像、自动控制、计算机视觉等。
推荐感兴趣的同学阅读参考文献中的书目。
流形优化基础知识
从黎曼流形的概念开始用 M \mathcal{M} M表示一个有限维的黎曼流形并用 T p M T_p\mathcal{M} TpM 表示 M \mathcal{M} M 在点 p p p处的切平面。在 M \mathcal{M} M的每一点 p p p上都有一个内积 ⟨ ⋅ , ⋅ ⟩ p \langle \cdot, \cdot \rangle_p ⟨⋅,⋅⟩p定义在切空间 T p M T_p\mathcal{M} TpM上。黎曼度量 ⟨ ⋅ , ⋅ ⟩ p \langle \cdot , \cdot \rangle_p ⟨⋅,⋅⟩p 相关的范数记作 ∥ ⋅ ∥ p \|\cdot\| _p ∥⋅∥p。用 $l(\gamma) $来表示分段光滑曲线 γ : [ a , b ] → M \gamma : [a, b] \rightarrow \mathcal{M} γ:[a,b]→M的长度。 流形优化的基本问题形式为 min x ∈ M f ( x ) \min_{x \in \mathcal{M}} f(x) x∈Mminf(x)
其中$f: \mathcal{M} \to \mathbb{R} \mathcal{M}$是一个黎曼流形。
常见流形
以下是一些在流形优化中常用的流形例子 球面流形 (Sphere Manifold) 考虑单位球面 S n − 1 S^{n-1} Sn−1 在 R n \mathbb{R}^n Rn 中定义为 S n − 1 { x ∈ R n : ∥ x ∥ 2 1 } S^{n-1} \{ x \in \mathbb{R}^n : \|x\|_2 1 \} Sn−1{x∈Rn:∥x∥21} 此流形上的点是具有单位范数的 R n \mathbb{R}^n Rn 中的向量。这是一个常见的约束特别是在需要正则化或单位范数约束的优化问题中。 正交流形 (Stiefel Manifold) 考虑所有有 n n n行和 p p p 列的列正交矩阵的集合。它可以表示为 S t ( n , p ) { X ∈ R n × p : X T X I p } St(n, p) \{ X \in \mathbb{R}^{n \times p} : X^T X I_p \} St(n,p){X∈Rn×p:XTXIp} 其中 I p I_p Ip 是 p × p p \times p p×p的单位矩阵。这种流形在子空间迭代和降维技术中很有用。 正定锥流形 (Positive Definite Manifold) 考虑所有 n × n n \times n n×n正定矩阵的集合。这是一个在半正定编程和统计中非常重要的约束集。 Grassmann 流形 定义为 n n n 维空间中所有 p p p 维线性子空间的集合。与 Stiefel 流形相比Grassmann 流形不考虑子空间的基的选取只考虑子空间本身。
基础概念
当我们考虑流形优化问题时特别是在黎曼流形上优化的算法需要考虑流形的几何结构。为此我们引入一些特定的概念如黎曼梯度、收缩映射retraction和向量传输vector transport。以下是这些概念的定义 黎曼梯度 (Riemannian Gradient): 对于定义在流形 M \mathcal{M} M上的实值函数 f f f其在点 p p p处的黎曼梯度记作 grad p f \text{grad}_p f gradpf是流形 M \mathcal{M} M上的一个切向量与欧氏梯度 ∇ f ( p ) \nabla f(p) ∇f(p)在该点的内积最大。数学上它满足以下关系 ⟨ grad p f , v ⟩ p D v f ( p ) , ∀ v ∈ T p M \langle \text{grad}_p f, v \rangle_p D_v f(p),\quad \forall v\in{T_p}\mathcal{M} ⟨gradpf,v⟩pDvf(p),∀v∈TpM
其中 D v f ( p ) D_vf(p) Dvf(p)是函数 f f f在点 p p p处沿着方向 v v v的方向导数而 ⟨ ⋅ , ⋅ ⟩ p \langle \cdot, \cdot \rangle_p ⟨⋅,⋅⟩p 是点 p p p处的黎曼度量。 收缩 (Retraction): 让我们从流形上的一点 p p p开始并考虑一个从该点出发的切向量 v v v。收缩映射 R p ( v ) R_p(v) Rp(v) 提供了一个近似的方法来移动到流形上另一点从而在某种意义上对黎曼指数映射进行了近似。 严格地说对于给定的点 p p p和切向量 v v v收缩映射 R p : T p M → M R_p: T_p\mathcal{M} \to \mathcal{M} Rp:TpM→M 满足以下两个条件: - R p ( 0 ) p R_p(0) p Rp(0)p - R p R_p Rp 在 0 0 0处是连续的 传输 (Transport):
传输定义了如何将流形上一点的切空间中的向量“移动”到另一点的切空间中而不改变其“方向”。设 T p M T_p\mathcal{M} TpM和 T q M T_q\mathcal{M} TqM 分别是流形 M \mathcal{M} M 上点 p p p和 q q q的切空间传输是一个映射 T p ( v ) : T p M → T q M T_p(v) : T_p\mathcal{M} \to T_q\mathcal{M} Tp(v):TpM→TqM 其中流形优化中两种常用的传输方式是平行传输parallel transport和向量传输vector transport。
常见算法
有了上面提及的这些概念类似欧式空间中的梯度下降法更新公式 x k 1 x k t k d k x^{k1} x^k t^kd^k xk1xktkdk我们可以简单的得到流形上的梯度下降算法的更新公式 x k 1 R x k ( − t k grad x k f ) x^{k1} R_x^{k}(-t^k \text{grad}_{x^{k}} f) xk1Rxk(−tkgradxkf) 其中 t k t^k tk是步长 R x k R_x^{k} Rxk是在点 x k x^{k} xk的收缩映射。
流形上的梯度下降法的关键在于如何正确地“移动”在流形上。由于流形的几何约束我们不能简单地进行线性更新而是需要利用流形的几何结构进行更新。黎曼梯度和收缩映射是这个过程中的核心组成部分。
近年来多种传统优化算法已成功地被扩展到流形优化领域包括梯度下降共轭梯度法随机梯度下降法牛顿法信赖域法内点法BFGS和ADMM等。需要注意的区别是流形上的梯度及Hessian矩阵信息也与流形本身有关黎曼梯度和黎曼Hessian矩阵同时更新过程也将依赖于收缩映射和传输这些操作确保算法在流形的几何约束内进行迭代并保持算法的稳定性和收敛性。因此对流形及其相关知识的深入了解和熟悉变得尤为重要。
流形优化计算工具
Manopt
Manopt 是一个专门为研究人员和工程师设计的优化工具箱旨在简化和推动定义在流形上的优化问题的解决过程。Manopt 提供了一系列工具和算法帮助用户轻松地在各种流形结构上设置和解决优化问题。Manopt 的优势在于它的灵活性和易用性。它支持多种流形结构如正定对称矩阵、Stiefel流形和球面流形等。此外它的设计目标是为用户提供一个简单、直观的接口使用户无需深入了解流形优化的复杂性就可以轻松应用。此外Manopt 还拥有丰富的文档和教程帮助新手快速上手同时为经验丰富的研究者提供高级功能和自定义选项。
PyManopt是Manopt的Python版本Manopt.jl是为Julia编程语言设计的Manopt版本。
实例
我们考虑Manopt官方给出的实例。计算对称矩阵 A ∈ R n × n A \in \mathbb{R}^{n\times n} A∈Rn×n 的主要特征向量。设 λ 1 ≥ ⋯ ≥ λ n \lambda_1 \geq \cdots \geq \lambda_n λ1≥⋯≥λn 是其特征值。最大的特征值 λ 1 \lambda_1 λ1是以下优化问题的最优值 max x ∈ R n , x ≠ 0 x ⊤ A x x ⊤ x . \max\limits_{x\in\mathbb{R}^n, x \neq 0} \frac{x^\top A x}{x^\top x}. x∈Rn,x0maxx⊤xx⊤Ax. 问题可以改写为 min x ∈ R n , ∥ x ∥ 1 − x ⊤ A x . \min\limits_{x\in\mathbb{R}^n, \|x\| 1} -x^\top A x. x∈Rn,∥x∥1min−x⊤Ax. 函数及其在 R n \mathbb{R}^n Rn 中的梯度为 f ( x ) − x ⊤ A x , ∇ f ( x ) − 2 A x . \begin{align} f(x) -x^\top A x, \nonumber\\ \nabla f(x) -2Ax. \nonumber \end{align} f(x)∇f(x)−x⊤Ax,−2Ax. 向量 x x x 的约束要求 x x x 具有单位2-范数即 x x x 是球面上的一点 S n − 1 { x ∈ R n : x ⊤ x 1 } . \mathbb{S}^{n-1} \{x \in \mathbb{R}^n : x^\top x 1\}. Sn−1{x∈Rn:x⊤x1}. 这是我们应用 Manopt 解决问题所需要的所有信息。
函数 f f f在 S n − 1 \mathbb{S}^{n-1} Sn−1 上是光滑的。它在 S n − 1 \mathbb{S}^{n-1} Sn−1上 x x x 处的黎曼梯度 是 x x x 处球的切线向量。它可以计算为从常规梯度 ∇ f ( x ) \nabla f(x) ∇f(x) 到该切线空间的投影使用正交投影器 P r o j x u ( I − x x ⊤ ) u \mathrm{Proj}_x u (I-xx^\top)u Projxu(I−xx⊤)u g r a d f ( x ) P r o j x ∇ f ( x ) − 2 ( I − x x ⊤ ) A x . \mathrm{grad}\,f(x) \mathrm{Proj}_x \nabla f(x) -2(I-xx^\top)Ax. gradf(x)Projx∇f(x)−2(I−xx⊤)Ax. 这是欧几里得梯度 ∇ f \nabla f ∇f 与黎曼梯度 g r a d f \mathrm{grad}\,f gradf 之间的数学关系的一个示例我们通常已经知道如何从微积分课程中计算 ∇ f \nabla f ∇f并且为优化需要 g r a d f \mathrm{grad}\,f gradf。在 Manopt 中转换是在后台通过一个名为 egrad2rgrad 的函数完成的我们只需要计算 ∇ f \nabla f ∇f。
以下是使用Manopt 中的信赖域法解决这个优化问题的代码 % 生成随机问题数据n 1000;A randn(n);A .5*(AA.);% 创建问题结构manifold spherefactory(n);problem.M manifold;% 定义问题成本函数及其欧几里得梯度。problem.cost (x) -x*(A*x);problem.egrad (x) -2*A*x;% 注意 e 在 egrad 中代表欧几里得% 数值检查梯度一致性可选checkgradient(problem);% 解决[x, xcost, info, options] trustregions(problem);% 显示一些统计数据figure;semilogy([info.iter], [info.gradnorm], .-);xlabel(迭代次数);ylabel(f 的梯度的范数);运行结果 在之后的系列中我们会对上述提到的流形优化中的常见基本概念和算法进行详细的分析。
参考文献
Absil P A, Mahony R, Sepulchre R. Optimization algorithms on matrix manifolds[M]. Princeton University Press, 2008.Boumal N. An introduction to optimization on smooth manifolds[M]. Cambridge University Press, 2023.Sato H. Riemannian optimization and its applications[M]. Berlin: Springer, 2021.潘汉黎曼流形优化及其应用, 北京, 科学出版社, 2020刘浩洋, 户将, 李勇锋文再文最优化建模、算法与理论, 第二版草稿, 7.4 流形约束优化算法 文章转载自： http://www.morning.ygkk.cn.gov.cn.ygkk.cn http://www.morning.drndl.cn.gov.cn.drndl.cn http://www.morning.brsgw.cn.gov.cn.brsgw.cn http://www.morning.ltpzr.cn.gov.cn.ltpzr.cn http://www.morning.3ox8hs.cn.gov.cn.3ox8hs.cn http://www.morning.yfffg.cn.gov.cn.yfffg.cn http://www.morning.lhgqc.cn.gov.cn.lhgqc.cn http://www.morning.bpxmw.cn.gov.cn.bpxmw.cn http://www.morning.hhfwj.cn.gov.cn.hhfwj.cn http://www.morning.yrngx.cn.gov.cn.yrngx.cn http://www.morning.jjzrh.cn.gov.cn.jjzrh.cn http://www.morning.wrdlf.cn.gov.cn.wrdlf.cn http://www.morning.nftzn.cn.gov.cn.nftzn.cn http://www.morning.mrckk.cn.gov.cn.mrckk.cn http://www.morning.xscpq.cn.gov.cn.xscpq.cn http://www.morning.lzrpy.cn.gov.cn.lzrpy.cn http://www.morning.tqdqc.cn.gov.cn.tqdqc.cn http://www.morning.fpjxs.cn.gov.cn.fpjxs.cn http://www.morning.ygkq.cn.gov.cn.ygkq.cn http://www.morning.rjljb.cn.gov.cn.rjljb.cn http://www.morning.ghssm.cn.gov.cn.ghssm.cn http://www.morning.fbzdn.cn.gov.cn.fbzdn.cn http://www.morning.kfbth.cn.gov.cn.kfbth.cn http://www.morning.nnqrb.cn.gov.cn.nnqrb.cn http://www.morning.yqyhr.cn.gov.cn.yqyhr.cn http://www.morning.mtmnk.cn.gov.cn.mtmnk.cn http://www.morning.kdrly.cn.gov.cn.kdrly.cn http://www.morning.dqwykj.com.gov.cn.dqwykj.com http://www.morning.rmdwp.cn.gov.cn.rmdwp.cn http://www.morning.fjzlh.cn.gov.cn.fjzlh.cn http://www.morning.kflbf.cn.gov.cn.kflbf.cn http://www.morning.ydnx.cn.gov.cn.ydnx.cn http://www.morning.aiai201.cn.gov.cn.aiai201.cn http://www.morning.sgfgz.cn.gov.cn.sgfgz.cn http://www.morning.fqklt.cn.gov.cn.fqklt.cn http://www.morning.frqtc.cn.gov.cn.frqtc.cn http://www.morning.ykmkz.cn.gov.cn.ykmkz.cn http://www.morning.nyplp.cn.gov.cn.nyplp.cn http://www.morning.tdmgs.cn.gov.cn.tdmgs.cn http://www.morning.ynjhk.cn.gov.cn.ynjhk.cn http://www.morning.wcgfy.cn.gov.cn.wcgfy.cn http://www.morning.qdzqf.cn.gov.cn.qdzqf.cn http://www.morning.bkqdg.cn.gov.cn.bkqdg.cn http://www.morning.fqzz3.cn.gov.cn.fqzz3.cn http://www.morning.plfy.cn.gov.cn.plfy.cn http://www.morning.ljjph.cn.gov.cn.ljjph.cn http://www.morning.ypwlb.cn.gov.cn.ypwlb.cn http://www.morning.rtzd.cn.gov.cn.rtzd.cn http://www.morning.nqypf.cn.gov.cn.nqypf.cn http://www.morning.tfwg.cn.gov.cn.tfwg.cn http://www.morning.zxqyd.cn.gov.cn.zxqyd.cn http://www.morning.tjkth.cn.gov.cn.tjkth.cn http://www.morning.mfsxd.cn.gov.cn.mfsxd.cn http://www.morning.jxgyg.cn.gov.cn.jxgyg.cn http://www.morning.bfnbn.cn.gov.cn.bfnbn.cn http://www.morning.rxrw.cn.gov.cn.rxrw.cn http://www.morning.hwtb.cn.gov.cn.hwtb.cn http://www.morning.wwkft.cn.gov.cn.wwkft.cn http://www.morning.gxcit.com.gov.cn.gxcit.com http://www.morning.rnzwh.cn.gov.cn.rnzwh.cn http://www.morning.lqrpk.cn.gov.cn.lqrpk.cn http://www.morning.wqpb.cn.gov.cn.wqpb.cn http://www.morning.dlgjdg.cn.gov.cn.dlgjdg.cn http://www.morning.pcgmw.cn.gov.cn.pcgmw.cn http://www.morning.kxqfz.cn.gov.cn.kxqfz.cn http://www.morning.lrnfn.cn.gov.cn.lrnfn.cn http://www.morning.kmwsz.cn.gov.cn.kmwsz.cn http://www.morning.ckhyj.cn.gov.cn.ckhyj.cn http://www.morning.jxzfg.cn.gov.cn.jxzfg.cn http://www.morning.nkyqh.cn.gov.cn.nkyqh.cn http://www.morning.wpmqq.cn.gov.cn.wpmqq.cn http://www.morning.kgnnc.cn.gov.cn.kgnnc.cn http://www.morning.ryglh.cn.gov.cn.ryglh.cn http://www.morning.xmtzk.cn.gov.cn.xmtzk.cn http://www.morning.qlrtd.cn.gov.cn.qlrtd.cn http://www.morning.pzrpz.cn.gov.cn.pzrpz.cn http://www.morning.phjny.cn.gov.cn.phjny.cn http://www.morning.gzgwn.cn.gov.cn.gzgwn.cn http://www.morning.mzhgf.cn.gov.cn.mzhgf.cn http://www.morning.xtqld.cn.gov.cn.xtqld.cn
