请问这段代码是干什么的？请帮我讲解一下，越详细越好

function [t,w,p,sr]=tp(X,b)
w=b/norm(b);
t=X*w;
p=X'*t/(t'*t);

Xtp=t*p';
Xr=X-Xtp;

%+++ Compute selectivity ratio
for i=1:size(X,2)
  vartp(i)=sumsqr(Xtp(:,i));
  varr(i)=sumsqr(Xr(:,i));
end
sr=vartp./(varr+eps); % sr.
%+++ end

你说它干嘛就干嘛，具体情况具体使用