目录

  • 算法
  • k的选择
  • \(x\)的初始化
  • 代码

抱歉,真的没怎么看懂,当然,估计和我现在没法静下心来好好看也有关系。

算法

想法非常非常简单吧,就是在原来幂法的基础上,每次迭代的时候再加个截断。当然,论文里给出了,为什么这么做的理由,把我弄得晕晕的,但是思想就是这么朴素。现在的问题是:
1.k怎么选?
2.初始\(x\)的选择

k的选择

这个我没在论文里找到,但是看数值实验,感觉在k上面是有操作空间的。

\(x\)的初始化

\(x\)的初始化,是这篇论文的大头,讲了怎么样怎么样就能怎么样怎么样。
总结就是有如下3种方案:

  • \(x=e_j,j=argmax\{A_{ii}\}\)实在是简单粗暴啊。
  • 分俩步,第一步先把\(k\)放大一些,然后进行迭代(初始化估计就用第一种吧),迭代几步之后,把\(k\)变回来,再继续迭代。
  • 当\(k\approx p\)的时候,采用Moghaddam et al. 2006后向选取的方法。

注:文章还提到,当\(A\)具有唯一的稀疏的主特征向量,那么,TPower方法能够从\(A\)的带噪声的观测中,讲该特征向量近似地恢复出来。

代码

def You_eig_value(C, x, k):  #幂法p = C.shape[1]x1 = x     #初始化while True:x2 = C @ x1truncate(x, k)x2 = x2 / np.sqrt(x2 @ x2)if np.sum(np.abs(x2-x1)) < 0.0001:breakelse:x1 = x2return x1def truncate(x, k): #截断p = len(x)label = np.argsort(np.abs(x))[:p-k]x[label] = 0

转载于:https://www.cnblogs.com/MTandHJ/p/10527942.html

Truncated Power Method for Sparse Eigenvalue Problems相关推荐

  1. Power Method for dominate eigenvalue

    Power method to solve eigenvalues #matlab code A = []; x = [];%%input A and intial x k = 1; [xp p] = ...

  2. A Physics-Data-Driven Bayesian Method for Heat Conduction Problems

    A Physics-Data-Driven Bayesian Method for Heat Conduction Problems Xinchao Jiang*, Hu Wang†, Yu li S ...

  3. Robust PCA 学习笔记

    很久没有写学习笔记了,年初先后忙考试,忙课程,改作业,回家刚安定下来,读了大神上学期给的paper,这几天折腾数学的感觉很好,就在这里和大家一起分享一下,希望能够有所收获.响应了Jeffrey的建议, ...

  4. scipy.sparse求稀疏矩阵前k个特征值

    背景: 要在python中处理7000*7000的稀疏矩阵,计算前k小的特征值和相应的特征向量.不想在matlab中做这件事了,所有的数据预处理和展现工作都想在python中完成.然而一般的linal ...

  5. 【Sparse R-CNN】《Sparse R-CNN:End-to-End Object Detection with Learnable Proposals》

    arXiv-2020 作者的博客:Sparse R-CNN 文章目录 1 Background and Motivation 2 Related Work 3 Advantages / Contrib ...

  6. Sparse R-CNN细节剖析

    文章目录 一:2D目标检测--问题探讨 二:Sparse细节探讨 三:Sparse R-CNN讲解 一:2D目标检测--问题探讨 当前目标检测已经有许多经典框架,大致可以分为三大类:Dense met ...

  7. CVPR 2011 全部论文标题和摘要

    CVPR 2011 Tian, Yuandong; Narasimhan, Srinivasa G.; , ■Rectification and 3D reconstruction of curved ...

  8. ICLR2020国际会议焦点论文(Spotlight Paper)列表(内含论文源码)

    来源:AINLPer微信公众号(点击了解一下吧) 编辑: ShuYini 校稿: ShuYini 时间: 2020-02-21     2020年的ICLR会议将于今年的4月26日-4月30日在Mil ...

  9. UBIFS - UBI File-System

    参考:http://www.linux-mtd.infradead.org/doc/ubifs.html#L_raw_vs_ftl UBIFS - UBI File-System Table of c ...

最新文章

  1. linux内核oom,linux OOM killer分析
  2. VS2010中经常使用的快捷键
  3. mongo 改字段名_一日一技:修改MongoDB集合中的字段名
  4. spring_security权限应用
  5. 基础理论:集合的Hausdorff距离
  6. hdu4821 字符串hash(有多少(M*L长的,M个不相同)子串))
  7. K8s普通用户配置权限解决User “system:anonymous“ cannot list nodes at the cluster scope
  8. 2021牛客暑期多校训练营7 B-xay loves monotonicity(线段树+不降子序列)
  9. linux算法平台,Linux实时调度算法与测试平台的研究与实现
  10. linux python3_在Linux上安装Python 3
  11. ubuntu 关闭IPv6
  12. 我对jeecg-boot项目的理解、使用心得和改进建议
  13. 13-[LVI-SAM]visual_odometry_初探
  14. TMS320F28035 中断中使用DINT,无法关闭中断的原因
  15. html怎么清除背景颜色,怎么去除Word复制网页内容的背景色?
  16. oracle让电脑变卡以及解决方法
  17. 头条号想过新手,但指数却一直达不到650,该怎么过新手?
  18. php slik转wav_node.js安装wx-voice转换silk和mp3
  19. OTFS-SCMA学习笔记
  20. 配置windows系统中 PHP的环境变量

热门文章

  1. Linux目录结构与Windows对比
  2. GCN-LSTM预测道路交通车辆速度 英文 Vehicle Speed Forecasting Based On GCN-LSTM Combined Model
  3. TOOD: Task-aligned One-stage Object Detection 原理与代码解析
  4. Oracle 错误总结及解决方法
  5. LCD1602 液晶
  6. Visual Studio 2017下载地址和安装教程(图解版)
  7. [附源码]java毕业设计在线学习网站的设计与实现
  8. ubuntu下mongodb及其可视化工具robomongo安装
  9. Oracle数据库:约束条件:主键约束、唯一约束、检查约束、非空约束、外键约束、默认值填写
  10. 不错的每日站会的一个实践