multBand?=?imread（‘C:\1111.tif‘）;
%?組合數
n?=?3;
%?波段數
band?=?size（multBand）;
combin?=?combntns（1:6n）;
sizeCombin?=?size（combin）;
OIF?=?zeros（sizeCombin（1）sizeCombin（2）+1）;
OIF（:1:end-1）?=?combin;
for?row?=?1:sizeCombin（1）
????bands?=?multBand（::combin（row:））;
????STDS?=?0;
????for?b?=?1:n
????????STDS?=?STDS?+?std2（bands（::b））;
????end
????CORRS?=?0;
????for?b1=