1.输入. 红外图像序列，$\mathcal{O}$，${\mathcal{B}}^0={\mathcal{S}}^0=0$，L，H，${\lambda }_1$，${\lambda }_2$，${\rho }_0=0.01$，${\rho }_{\mathrm{m}\mathrm{a}\mathrm{x}}=1e7$，$\mu =1.5$，$\zeta =1e-6$

  2.重复以下步骤，直到满足停止条件：

  步骤1. 根据式（11）更新${\mathrm{U}}_{\mathrm{n}}$

  步骤2. 根据式（14）更新${\mathrm{W}}_{\mathrm{n}}$

  步骤3. 根据式（16）更新$\mathcal{B}$

  步骤4. 根据式（19）更新$\mathcal{G}$

  步骤5. 根据式（21）更新$\mathcal{S}$

  步骤6. 根据式（16）更新$\mathcal{B}$

  步骤7. 根据式（14）更新${\mathrm{W}}_{\mathrm{n}}$

  步骤8. 根据式（11）更新${\mathrm{U}}_{\mathrm{n}}$

  步骤9 根据式（22）更新${\left\{{\delta }_n\right\}}_{n=1}^N$，$\mathrm{\Gamma }$，$\mathrm{\Lambda }$

  步骤10. 根据式（23）更新$\mathit{\rho }$

  步骤11. 检查是否达到收敛条件$\frac{{‖\mathcal{O}-\mathcal{B}-\mathcal{S}‖}_F^2}{{‖\mathcal{O}‖}_F^2}\le \zeta$

  3.输出. 背景图像$\mathcal{B}$，目标图像$\mathcal{S}$