使用 Ax=b 中算子 A 的函数的稀疏性迭代算法

计算科学 稀疏矩阵 迭代法
2021-12-11 03:26:35

我要解决一个线性迭代逆问题。我有两个功能,matlab其中一个扮演角色forward,另一个扮演adjoint角色。我对逆问题不太熟悉,我将测试一些算法来检查哪些算法可以帮助我获得稀疏模型。由于我必须使用我的函数来执行正向和反向运算符角色Ax=b,所以我认为我无法使用一些需要A作为矩阵的算法。例如,我可以使用linearized bregman如下:

vk+1=vk+AT(fAuk).

uk+1=δshrink(vk+1,1/μ).

其中,A我使用将模型u转换为数据空间A^T的函数以及将数据转换为模型空间的函数。

我想知道是否有任何其他迭代算法可以以这种方式用于稀疏性?

1个回答

是的,大多数稀疏正则化最小二乘问题的求解器都被设计为使用矩阵向量乘法,而不是直接访问矩阵。参见例如 SPGL1 在

http://www.math.ucdavis.edu/~mpf/spgl1/

可以在 Rice 网站上的压缩传感资源上找到解决此类问题的软件列表(在页面的最底部):

http://dsp.rice.edu/cs