怎么样在真实推荐系统中应用矩阵分解算法

怎么样在真实推荐系统中应用矩阵分解算法

里面有自带的LU分解函数，你可在帮助文档搜索 LU给你一个我自己写的LU分解函数function [myl,myu,x]=MYLU(A,b)%依据《现代电力系统分析》编制的LU分解程序。%matlab自带的lu函数与书上所讲略有不同，不方便参照课本步骤进行后续计算。%A为待分解的矩阵，myl为分解后下三角矩阵，myu为分解后的上三角矩阵,未考虑单独的对角阵D%clc%A=[5,0,0,5,1;2,1,1,0,0;0,1,2,0,10;1,0,5,1,0;1,1,0,3,10] %两组测试用矩阵方程%b=[1,3,3,4,1]';%M=[1,2,1,1;2,1,0,0;1,0,1,0;1,0,0,1]%b=[-1,1,2,0]'%A=[2,3,1;3,7,-1;5,-4,2]%b=[12,13,5]'n=length(A);myl=speye(n); %可以先不分配内存，由系统自动分配内存myu=speye(n);%==============LU分解================for i=1:n for j=1:n if (i=j) x=0; for k=1:j-1 x=x+myl(i,k)*myu(k,j); end myl(i,j)=A(i,j)-x; end endend%B=myl*myu %测试LU分解矩阵，B=A表明结果正确%===============解方程Ax=b================y=(zeros(1,n))';x=(zeros(1,n))';y(1)=b(1)/myl(1,1);%解Ly=bfor i=2:n p=0; for k=1:i-1 p=p+myl(i,k)*y(k); end y(i)=(b(i)-p)/myl(i,i);end%C=(myl*y) %测试Ly=b，如果相等表明结果正确x(n)=y(n);for i=1:n-1 p=0; for k=n-i+1:n p=p+myu(n-i,k)*x(k); end x(n-i)=y(n-i)-p;end%D=myl*myu*x里面注释部分可以删除