输入：稀疏张量 $X\in R^{I1\times I2\times I3}$ ，分解阶数R

输出：系数矩阵 $A\in R^{I1xR}$ ， $B\in R^{I2xR}$ ， $C\in R^{I3xR}$

将A，B，C初始化为大小为 $I_1\times R$ ， $I_2\times R$ ， $I_3\times R$ 的随机矩阵

重复以下步骤直到收敛：

1) 固定B，C，更新A：

计算 $A=X{\times }_2\text{khatri\_rao}\left(C,B\right)\times \text{pinv}\left(\text{khatri\_rao}\left(C,B\right)\right)\times A$

其中，X₂表示对应于X的第2个维度的张量乘法运算，khatri_rao表示Khatri-Rao乘积，pinv表示伪逆。

2) 固定A，C，更新B：

计算 $B=X{\times }_1\text{khatri\_rao}\left(C,A\right)\times \text{pinv}\left(\text{khatri\_rao}\left(C,A\right)\right)\times B$

其中，X₁表示对应于X的第1个维度的张量乘法运算，khatri_rao表示Khatri-Rao乘积，pinv表示伪逆。

3) 固定A，B，更新C：

计算 $C=X{\times }_3\text{khatri\_rao}\left(B,A\right)\times \text{pinv}\left(\text{khatri\_rao}\left(B,A\right)\right)\times C$

其中，X₃表示对应于X的第3个维度的张量乘法运算，khatri_rao表示Khatri-Rao乘积，pinv表示伪逆。

返回系数矩阵A，B，C