function r=incident_degree(x0,x1) 
%compute the incident degree for grey model. 

%Designed by NIXIUHUI,Dalian Fisher University. 
%17 August,2004,Last modified by NXH at 21 August,2004        

%數據初值化處理 
x0_initial=x0./x0(1);
temp=size(x1); 
b=repmat(x1(:,1),[1 temp(2)]); 
x1_initial=x1./b; 

%分辨系數選擇 
K=0.1; 

disp('The grey interconnect degree is: '); 


x0_ext=repmat(x0_initial,[temp(1) 1]); 
contrast_mat=abs(x0_ext-x1_initial); 

delta_min=min(min(contrast_mat));%delta_min在數據初值化后實際為零 
delta_max=max(max(contrast_mat)); 
a=delta_min+K*delta_max; 
incidence_coefficient=a./(contrast_mat+K*delta_max);%得到關聯系數 

r=(sum(incidence_coefficient'))'/temp(2);      %得到鄧氏面積關聯度