矩阵补全算法

计算科学 线性代数 凸优化
2021-12-13 06:39:23

我正在尝试实现本文中提出的算法,该算法试图恢复一个矩阵,该矩阵表示一组图像像素强度的噪声较小的数据集。在这种情况下,噪声来自图像中物体材料的镜面反射和阴影。

我在凸优化方面几乎没有(或没有)经验。我的问题是,算法中的停止标准应该是什么?

我使用其中之间的差异,但结果该范数的值保持不变迭代后迭代。||Ai+1Ai||FFA

有人可以给我一些建议以实现该算法吗?

编辑

重新编码算法后,我按照其中一个答案的建议检查了两个条件,这就是绘图的样子。

对于我得到 ||Ei+1Ei||FE的范数之间的差异

对于我得到 ||Ai+1Ai||FA的范数之间的差异

我不知道错误可能出在哪里,我一步一步地检查了代码,这对我来说似乎是正确的。

1个回答

我不确定这个算法是如何工作的,但我猜是一个相当合理的停止标准,而且似乎是需要的。||AiAi+1||FϵA||EiEi+1||FϵE

如果 都没有减少,那么你在算法的两个内部循环中犯了一些大错误。||AiAi+1||F||EiEi+1||F