clc,clear;
[rows,cols]=size();
p=0.5; %分辨系数
[m,n]=size(yangben);
R=[];
for irow=1:rows
yy=fangzhen(irow,:);
data=[yy;yangben];
data_gyh1=mean()
for i=1:m+1
for j=1:n
data_gyh(i,j)=data(i,j)/data_gyh1(j);
end
end
for i=2:m+1
for j=1:n
Dij(i-1,j)=abs(data_gyh(1,j)-data_gyh(i,j));
end
end
Dijmax=max(max(Dij));
Dijmin=min(min(Dij));
for i=1:m
for j=1:n
Lij(i,j)=(Dijmin+p*Dijmax)/(Dij(i,j)+p*Dijmax);
end
end
LijRowSum=sum(Lij');
for i=1:m
Rij(i)=LijRowSum(i)/n;
end
R=[R;Rij];
end
R