X0 = xlsread('boiler','1#炉','B2:B3600');%热效率
X1 = xlsread('boiler','1#炉','C2:E3600');%煤气流量
X2 = xlsread('boiler','1#炉','D2:G3600');%空气流量
X3 = xlsread('boiler','1#炉','E2:I3600');%热值
X4 = xlsread('boiler','1#炉','F2:K3600');%蒸汽流量
X5 = xlsread('boiler','1#炉','G2:M3600');%给水流量
y = [X0,X1,X2,X3,X4,X5]';
y = mapminmax(y);
y1=mean(y');
y1=y1';
for i=1:size(y,1)for j=1:size(y,2)y2(i,j)=y(i,j)/y1(i);end
end
for i=2:size(y,1)for j=1:size(y,2)y3(i-1,j)=abs(y2(i,j)-y2((i-1) ,j));end
end
a=1;b=0;
for i=1:size(y,1)-1for j=1:size(y,2)if (y3(i,j)<=a) a=y3(i,j); elseif (y3(i,j)>=b)b=y3(i,j);endend
end
for i=1:size(y,1)-1for j=1:size(y,2)y4(i,j)=(a+0.5*b)/(y3(i,j)+0.5*b);end
end
y5=sum(y4')/(size(y,2)-1);