我正在对成像进行非线性优化,例如 MRI 和 CT。
我们的问题是的形式。A从未明确形成,因此我们仅限于仅使用Ax和A^Hf 的方法。即便如此,这些评估还是相当昂贵的,而Wx和W^Tx的计算成本却相当低。
我一直在使用非线性共轭梯度 (NCG) 来解决这个问题,但是我的回溯线搜索往往会失败,因为不是 A 的精确共轭转置而是一个近似值。
忽略非线性项并在正规方程上使用共轭梯度,一切都很好,Split-Bregman 算法也是如此。然而,后者比 NCG 慢一个数量级。
我想知道是否有人知道允许梯度(近似的迭代求解器?