function [y a b]=guiyihua(r1)
[n N]=size(r1);
k1=0.9;%縮小系數(shù)
k2=1.1;%放大系數(shù)
a=max(r1).*k2;b=min(r1).*k1;s=a-b;
for i=1:n
    for j=1:N
        r1(i,j)=(r1(i,j)-b(1,j))/s(1,j);
    end
end
y=r1;