p1=0;%求分子项
rp2=0;%求分母项
for i=1:m
for j=1:n
for k=1:n
rp1=rp1+abs(r(i,k)*p(j,k));
end
rp2=rp2+rp1;
a=rp1/rp2;