clc,clear,close all
y=[ 102.4 102.8 103.1 102.9 103.3 103.5 103.6 104.4 105.1 104.6 104.9 104.9 105.4
105.2 105.9 106.1 105.7 106.8 107.5 108 110.1 117.1 109.6 110.3 111 111.7
101.7 101.7 101.7 101.7 101.6 101.5 101.4 101.5 101.6 101.8 101.8 101.9 102.1
98.9 98.7 98.8 99 99.2 98.8 98.5 98.7 99.3 100.1 99.8 100.4 100.8
99.3 99.5 99.7 100 100.2 100.4 100.4 100.5 100.7 101.2 101.4 101.4 101.9
102.5 102.8 103.2 103.2 103.3 103.3 103.4 103.7 104 104 103.2 103 103.2
100 100 100.1 99.7 99.3 99.4 99.3 99.5 99.3 99.3 99.9 99.7 100.1
100.3 100.4 100.6 100.9 101.1 101.2 101.2 100.9 100.6 100.7 101 100.3 100.5
100.3 104.5 105 105 104.8 104.4 104.3 104.9 105.8 106 106.8 106.1 106.6];
y1=mean(y');
y1=y1';
for i=1:9
for j=1:13
y2(i,j)=y(i,j)/y1(i);
end
end
for i=2:9
for j=1:13
y3(i-1,j)=abs(y2(i,j)-y2((i-1) ,j));
end
end
a=1;b=0;
for i=1:8
for j=1:13
if(y3(i,j)<=a) a=y3(i,j);
elseif(y3(i,j)>=b) b=y3(i,j);
end
end
end
for i=1:8
for j=1:13
y4(i,j)=(a+0.5*b)/(y3(i,j)+0.5*b);
end
end
y5=sum(y4')/12