% wavenum = csvread('wavenumber.csv');
% wavenum =wavenum';
zerod = csvread('1old-peiru-Absorbance.csv');
zerod = zerod(1:338,2:3);
zerod = zerod';
N = 338;
L = 10;
E=zeros(1,338);
for i = 1 : N - L
window = zerod(:, i: i + L);
cor = corrcoef(window');
S=0;
for j= 1:2
for k =1:2
S = S + cor(j,k);
end
end
E(i) = (S-2)/2;
end
plot(E);
title('移动窗口特征波段选取');
xlabel('波段(wavenumber:0.2-80.85cm-1)');
ylabel('Correlation coefficient')