function?VV=MSR（S）
????if?length（size（S））>2
????????S=rgb2gray（S）;
????end
????[mn]=size（S）;
????S=double（S）;

????sizes=m*n;
????mean=sum（sum（S）2）/sizes;
????R=zeros（mn）;
????R=repmat（meanmn）;

????%N（ij）高斯函數
????c=15;%尺度一
????for?i=1:m
????????for?j=1:n
????????????N（ij）=exp（-（i^2+j^2）/c^2）;
????????end
????end
????K=1/sum（N（:））;
????N1=K*N;
????G1=conv2（SN1‘same‘）;??

????c=80;%尺度二
????for?i=1:m
????????for?j=1:n
????????????N（ij）=exp（-（i^2+j^2）/c^2）;
????????end
????end
????K=1/sum（N（:））;
????N1=K*N;
????G2=conv2（SN1‘same‘）;??

????c=200;%尺度三
????for?i=1:m
????????for?j=1:n
????????????N（ij）=exp（-（i^2+j^2）/c^2）;
????????end
????end
????K=1/sum（N（:））;
????N1=K*N;
????G3=conv2（SN1‘sam