Input: , , , , , and . | Output: , , , , and . | Initialization: , , , and . | while not converged and do | Step 1. Given , , , and , computing by | , (a) | where and denote the discrete Fourier transform and the inverse discrete Fourier | transform, respectively, and Fourier transforms of operators and are regarded as the | transforms of their corresponding convolution kernels. | Step 2. Given , , and , update and by the two-dimensional shrinkage | , (b) | , (c) | respectively, where is assumed. | Step 3. Given , , and , update and by | , (d) | , (e) | respectively. | end while |
|