matlab实现CP分解 代码

使用tensor-box，在matlab中实现cp分解

1、X = sptenrand([5 4 3], 10)
//随机生成一个稀疏三维张量5*4* 3，其中有10条数据

X is a sparse tensor of size 5 x 4 x 3 with 10 nonzeros
(1,4,1) 0.4966
(2,2,3) 0.8998
(3,2,3) 0.8216
(3,3,1) 0.6449
(3,3,3) 0.8180
(3,4,1) 0.6602
(4,1,2) 0.3420
(4,1,3) 0.2897
(5,2,2) 0.3412
(5,3,2) 0.5341

2、P = parafac_als(X,2)
//用ALS交替最小二乘来分解张量，分解参数R=2，即分解为两个秩一张量；其中分解参数最多可以设置为min（5*4；5* 3；4*3）

CP_ALS:
Iter 1: fit = 3.219563e-001 fitdelta = 3.2e-001
Iter 2: fit = 3.645517e-001 fitdelta = 4.3e-002
Iter 3: fit = 3.732887e-001 fitdelta = 8.7e-003
Iter 4: fit = 3.809608e-001 fitdelta = 7.7e-003
Iter 5: fit = 4.021826e-001 fitdelta = 2.1e-002
Iter 6: fit = 4.427524e-001 fitdelta = 4.1e-002
Iter 7: fit = 4.734919e-0