incident_degree.m

用于mann-kendall趋势分析

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);      %得到邓氏面积关联度