Singular Value Thresholding (SVT) 奇异值阈值

Singular Value Thresholding (SVT) 奇异值阈值这个算法受到压缩感知中迭代算法的启发,在迭代过程中对矩阵进行SVD,然后将较小的奇异值设置为0,生成新的矩阵进行迭代。该算法运算速度快,对于高位低秩矩阵的恢复非常有效。

大家好,又见面了,我是你们的朋友全栈君。

[本文链接:http://blog.csdn.net/shanglianlm/article/details/46009387,转载请注明出处]

为了求解问题

这里写图片描述

因为它是非凸的,我们求解一个它的近似算法

这里写图片描述

对于一个大的 τ 值,它可以用下列等式接近

这里写图片描述

其中第一项为核范式(奇异值的和),第二项为Frobenius范式。

  1. Singular Value Thresholding (SVT) 奇异值阈值

    * 奇异值收缩(singular value shrinkage)*

    首先我们考虑一个秩为 r 的矩阵

    XRn1xn2
    的奇异值分解如下:
    SVD
    其中 U

    V
    分别为 n1×r n2×r 的正交矩阵,奇异值为 ρi 非负的。

    对于每个 τ0 ,我们有软阈值操作 Dτ :
    SVS
    其中 t+ 表示的 t 非负部分,即

    t+=max(0,t)
    。换句话说,这个软阈值操作仅仅应用于矩阵 X <script type=”math/tex” id=”MathJax-Element-56″>X</script> 的奇异值上,使它们趋于零。这也是为什么我们将其成为奇异值收缩(singular value shrinkage)的原因。

    * Singular Value Thresholding (SVT) 奇异值阈值*

    又因为奇异值收缩(singular value shrinkage)是核范式的近似操作(具体证明见[3]),因此上式可以转化为:
    这里写图片描述

    它的迭代方式为:
    这里写图片描述

    这个算法受到压缩感知中迭代算法的启发,在迭代过程中对矩阵进行SVD,然后将较小的奇异值设置为0,生成新的矩阵进行迭代。该算法运算速度快,对于高位低秩矩阵的恢复非常有效。

  2. 用拉格朗日乘子法解释

    原问题为:

    这里写图片描述

    其拉格朗日函数为:

    这里写图片描述

    强对偶成立,且拉格朗日函数的鞍点是原函数与对偶问题的最优解,即

    这里写图片描述

    其迭代解为:

    这里写图片描述

参考或延伸材料
[1] 斯坦福SVT软件
[2] Generalized Singular Value Thresholding
[3] A singular value thresholding algorithm for matrix completion
[4] Exact Matrix Completion via Convex Optimization

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/137392.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • webstorm的永久激活码2021-激活码分享[通俗易懂]

    (webstorm的永久激活码2021)最近有小伙伴私信我,问我这边有没有免费的intellijIdea的激活码,然后我将全栈君台教程分享给他了。激活成功之后他一直表示感谢,哈哈~IntelliJ2021最新激活注册码,破解教程可免费永久激活,亲测有效,下面是详细链接哦~https://javaforall.net/100143.html…

    2022年3月22日
    3.0K
  • particles.js使用及配置

    particles.js使用及配置particles.js使用及配置参考:http://blog.csdn.net/csdn_yudong/article/details/53128570这个项目中有提供demo,可以直接下载这个

    2022年8月5日
    5
  • “国产Linux“的终结

    “国产Linux“的终结

        现今,“国产Linux”不见了,也不再有人提起了,人们的思想获得了大解放。此话根据何在?
     
         6月29日,《2010‘开源中国,开源世界》高峰论坛在北京举行。高峰论坛的议事议程及其特意安排充分显示出,在我国开源发展历史上,这次会议是一个具有特殊意义的发展里程碑。近半年来,我国

    2022年5月13日
    39
  • JAVA异常_什么什么异常

    JAVA异常_什么什么异常异常问题如下:起初xml中返回值类型是这样子,一直在找返回值类型的问题,怎么看都是没有问题的又改为如下,结果还是不对,查询资料反反复复还是出现这个异常突然一下想到会不会是有重复id名字的sql 我用的是idea,直接全局查询,确实查到了一模一样在别的包下的sql,因为我项目是maven的子工程,就给冲突了;最后还是冒着尝试想法去试了一下,把名字一个,哦吼!美滋滋 问题解决 不在报错以后开发还是需要多注意名字重复的问题的!…

    2022年8月19日
    7
  • matlab读取h5文件变慢_matlab载入数据文件

    matlab读取h5文件变慢_matlab载入数据文件filename=’C:\Users\88304\Desktop\Retina-Unet-1\test\test_last_weights.h5′;h5disp(filename)

    2025年8月27日
    6
  • WPFのDecorator 、Adorner和AdornerDecorator

    WPFのDecorator 、Adorner和AdornerDecorator

    2021年6月16日
    99

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号