[发明专利]一种概率矩阵分解推荐方法

权利要求书

1.一种概率矩阵分解推荐方法，其特征在于该方法包括以下步骤：

步骤1)获得用户在线社交网络中的信任关系矩阵与用户商品评分矩阵；所述在线社交网络是在互联网上与其他人相联系的一个平台，用户对产品进行评分，同时分享给该用户的朋友，查询该用户的朋友的评分；

步骤2)随机生成U和Z，所述U∈R^d×m表示隐含用户特征矩阵，Z∈R^d×m表示隐含信任关系特征矩阵，表示d行m列的矩阵，d是用户根据经验指定的隐含特征数，m是用户的个数，将隐含用户特征矩阵和隐含信任关系特征矩阵的先验分布表示为：

p(U|σU2)=Πi=1mN(Ui|0,σU2)]]>

p(Z|σZ2)=Πk=1mN(Zk|0,σz2)]]>

所述i表示用户个数的变量，k表示用户个数的变量，U_i表示用户ui的特征列向量，Z_k表示第k个信任关系特征向量，表示是均值为0，方差为的高斯分布的概率密度函数，表示是均值为0，方差为的高斯分布的概率密度函数；

步骤3)对用户信任关系矩阵进行概率矩阵分解，用户信任关系矩阵的条件概率分布表示为：

p(T|U,Z,σT2)=Πi=1mΠk=1m[N(ti_k*|g(UiTZk),σT2)]Ii_k_T]]>

所述T表示m×m维的用户信任关系矩阵，表示U_i的转置，表示均值为方差为的高斯分布的概率密度函数，I_{i_k_T}是一个变量，表示用户u_i与用户u_k之间信任关系，当用户u_i信任用户u_k时，I_{i_k_T}＝1，否则I_{i_k_T}＝0；所述是将的值映射在[0,1]之间，g(UiTZk)=11+e(UiTZk),]]>是一个变量，ti_k*=d-(uk)d+(ui)+d-(uk)×ti_k,]]>t_{i_k}表示用户u_i与用户u_k之间的信任权值，d⁺(u_i)表示用户u_i信任的用户数量，d^-(u_k)表示用户u_k被信任的用户数量；

步骤4)随机生成V，所述V∈R^d×n表示隐含商品特征矩阵，表示d行n列的矩阵，n表示商品的个数，将隐含商品特征矩阵的先验分布表示为：

p(V|σV2)=Πj=1nN(Vj|0,σV2)]]>

所述j表示商品个数的变量，V_j表示商品i_j的特征向量，表示是均值为0，方差为的高斯分布的概率密度函数；

步骤5)对用户的商品评分矩阵进行概率矩阵分解，用户的商品评分矩阵的条件概率分布表示为：

p(R|U,V,σR2)=Πi=1mΠj=1m[N(ri_j|g(UiT(Vj+|Iu|-12Σj∈Iuyj)),σR2)]Ii_j_R]]>

所述R表示m×n维的用户的商品评分矩阵，r_{i_j}表示用户u_i对商品i_j的评分，N(ri_j|g(UiT(Vj+|Iu|-12Σj∈Iuyj)),σR2)]]>表示均值为g(UiT(Vj+|Iu|-12Σj∈Iuyj)),]]>方差为的高斯分布的概率密度函数，I_u表示被用户评过分的商品集合，|I_u|表示被用户评过分的商品的数量，I_{i_j_R}表示用户u_i是否给商品i_j评过分，若评过分，则I_{i_j_R}＝1，否则I_{i_j_R}＝0，y_j表示已被用户评过分的商品对即将要评分商品i_j的隐性影响的权重值，初始值为1.0；

步骤6)分别计算U_i、V_j、Z_k及y_j的更新值U_i′、V_j′、Z_k′与y_j′，Ui′=Ui-γ·∂L∂Ui,]]>Vj′=Vj-γ·∂L∂Vj,]]>Zk′=Zk-γ·∂L∂Zk,]]>yj′=yj-γ·∂L∂yj,]]>所述γ为预先定义的步长，γ足够小且γ＞0，γ的具体值根据经验确定，其中∂L∂Ui=Σj=1nIi_j_Rg′(UiT(Vj+|Iu|-12Σj∈Iuyj))×(g(UiT(Vj+|Iu|-12Σj∈Iuyj)-ri_j))×(Vj+|Iu|-12Σj∈Iuyj)+λCΣi=1mIi_k_Tg′(UiTZk)(g(UiTZk)-ti_j*))Zk+λUUi,]]>∂L∂Vj=Σi=1mIi_j_Rg′(UiT(Vj+|Iu|-12Σj∈Iuyj)×(g(UiT(Vj+|Iu|-12Σj∈Iuyj)-ri_j))Ui+λVVj,]]>

∂L∂Zk=λCΣi=1mIi_k_Tg′(UiTZk)g(UiTZk)-ti_j*)Ui+λZZk,]]>

∂L∂yj=Σi=1mIi_j_Rg′(UiT(Vj+|Iu|-12Σj∈Iuyj)×(g(UiT(Vj+|Iu|-12Σj∈Iuyj)-ri_j))Ui|Iu|-12+λyj]]>

所述λC=σR2/σT2,]]>λU=σR2/σU2,]]>λV=σR2/σV2,]]>λZ=σR2/σZ2,]]>λ为预先定义的规则化参数，λ足够小且λ＞0，λ的具体值根据经验确定；

步骤7)根据公式

L=12Σi=1mΣj=1nIi_j_R(ri_j-g(UiT(Vj+|Iu|-12Σj∈Iuyj))2)+λC2Σi=1mΣk=1nIi_k_T(ti_j*-g(UiTZk)2)+λU2||U||F2+λV2||V||F2+λZ2||Z||F2+λ2||yj||F2]]>

计算目标函数L的值，当目标函数L的值变化小于某个预先定义的很小的常数或在经过设定的迭代次数后终止迭代过程，否则令U_i＝U_i′，V_j＝V_j′，Z_k＝Z_k′，y_j＝y_j′，返回步骤6)，所述是欧几里得向量范数；

步骤8)当迭代运算终止后，得到U_i′、V_j′以及y_j′，使用r^i_j=Ui′T(Vj′+|Iu|-12Σj∈Iuyj′)]]>计算用户u_i对商品i_j的未知评分