Input：$\mathit{D},{\mathit{W}}_{\mathrm{d}},{\mu }^0,\lambda =0.7/\sqrt[]{\mathrm{m}\mathrm{a}\mathrm{x}(n_1,n_2)\mathrm{*}{n}_3},\epsilon ={10}^{-7}$

Initialization：${\mathit{B}}^0={\mathit{T}}^0=y^0=0,{\mathit{W}}_{\mathrm{s}\mathrm{w}}=1,{\mathit{W}}^0={\mathit{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}\odot {\mathit{W}}_{\mathrm{s}\mathrm{w}},{\mu }^0=3\times {10}^{-6},\rho =1.1,c=1,k=0$

while not converge do

1.Fix the others and update ${\mathit{T}}^{k+1}$：${\mathit{T}}^{k+1}={\mathit{S}}_{\lambda {\mathit{W}}^k{\mu }^{-1}}\left(\mathit{D}-{\mathit{B}}^k-\frac{y^k}{\mu }\right)$

2.Fix the others and update ${\mathit{B}}^{k+1}$：${\mathit{B}}^{k+1}=D_{\mathit{W}{\mu }^{-1}}\left(\mathit{D}-{\mathit{T}}^k-\frac{y^k}{\mu }\right)$

3.Fix the others and update $y^{k+1}$：$y^{k+1}=y^k+{\mu }^k(\mathit{D}-{\mathit{B}}^{k+1}-{\mathit{T}}^{k+1})$

4.Fix the others and update ${\mathit{W}}^{k+1}$：${\mathit{W}}_{\mathrm{s}\mathrm{w}}^{k+1}=\frac{c}{\left|{\mathit{T}}^k\right|+\epsilon };{\mathit{W}}^{k+1}={\mathit{W}}_{\mathrm{r}\mathrm{e}\mathrm{c}}\odot {\mathit{W}}_{\mathrm{s}\mathrm{w}}^{k+1}$

5.Update$\mu$：${\mu }^{k+1}=\rho {\mu }^k$

6.Check the convergence conditions：$\frac{{‖{\mathit{B}}^{k+1}+{\mathit{T}}^{k+1}-\mathit{D}‖}_F^2}{{‖\mathit{D}‖}_F^2}<\epsilon \mathrm{o}\mathrm{r}{‖{\mathit{T}}^{k+1}‖}_0={‖{\mathit{T}}^k‖}_0$

7.Update$k$：$k=k+1$

end while

Output：${\mathit{B}}^k,{\mathit{T}}^k$

19.946

2.835

88.916

inf

inf

inf

inf