输入：高光谱图像观测矩阵X，参数 $\alpha$ 、 $\mu$ 、 $\lambda$ 、 $\sigma$ ，阈值 $\tau$ ，最大迭代次数 $T_{\max }$

输出：端元矩阵A和丰度矩阵S。

步骤1：通过VD算法预估端元数量P。

步骤2：通过VCA-FCLS算法进行初始化，获得A₀和S₀，并计算初始重构误差F₀。

步骤3：如式(12)~(15)采用乘性迭代法则交替更新A和S。

步骤4：计算每次迭代的图像重构误差F_K。

步骤5：根据式(16)判断迭代是否停止。如不满足迭代停止条件，重复步骤3~4。