灰色關聯(lián)分析

function output=grayrela(x0)

%參考因子與比較因子共同存儲在一個矩陣x0中,參考因子位于第一列

%斜率序列

for i=2:length(x0(:,1))

x1(i,:)=x0(i,:)-x0(i-1,:);

end

%標準化

m=length(x1(1,:));

for i=1:m

x2(:,i)=x1(:,i)/std(x1(:,i));

end

%排序

[y,pos]=sort(x2(:,1));

x2_sorted=x2(pos,:);

% 判定關聯(lián)性質

n=length(x1(:,1));

k=[1:n]';

for j=1:m

sig_j(j)=qiuhe(k.*x2_sorted(:,j))-qiuhe(x2_sorted(:,j))*qiuhe(k)/n;

end

%caculation of distantion

for j=2:m

dist_0i(:,j)=abs(sign(sig_j(:,j)./sig_j(:,1)).*x2_sorted(:,j)-x2_sorted(:,1));

end

%計算關聯(lián)系數(shù)

for i=1:n

for j=1:m

coef_rela(i,j)=(min(dist_0i)+0.5*max(dist_0i))/(dist_0i(i,j)+0.5*max(dist_0i));

end

end

for j=1:m

output(j)=qiuhe(coef_rela(:,j))/n;

end

其中：

function output=qiuhe(input)

output=0;

for i=1:length(input)

output=output+input(i);

end

計算方法參考文獻：王寧練：冰川平衡線變化的主導氣候因子灰色關聯(lián)分析 冰川凍土

總結

以上是生活随笔為你收集整理的matlab灰色关联代码,灰色关联分析matlab代码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯，歡迎將生活随笔推薦給好友。