xiangguangxing.m

来自「稀土测量 matlab 源程序 粒子群优化算法的改进 lS-SVM等源程序」· M 代码 · 共 22 行

%相关性分析，n为数据长度 

clc
clear
load input
y=input(:,5);
x=input(:,4);
n=150;
sxx1=0;
sxy1=0;
syy1=0;
syy2=mean(y);
for i=1:150
   sxx1= sxx1+x(i)^2;
   sxy1=sxy1+x(i)*y(i);
   syy1=syy1+y(i)^2;
end
sxx=sxx1-((n*mean(x))^2)/n;%sxx
sxy=sxy1-mean(x)*mean(y)*n;
syy=syy1-((n*mean(y))^2)/n;
r=sxy/(sqrt(sxx*syy));%相关系数