Pearson线性相关系数计算公式

为了计算股市的一些数值，我这两天一直在钻研线性代数，准备写一个小工具来管理自己的基金。

但是自己的线代忘记得差不多了，没办法只有再拣起来了。先复习一下关于Pearson线性相关系数的计算公式，呵呵，google了好久才找到的。

经过一定简化后得到公式：

根据简化后公式写的PHP函数如下：

function pearson($x,$y)
{
if(!is_array($x) || !is_array($y))
return 0;
if(count($x)!=count($y))
return 0;

$n=count($x);

$x2=0;$y2=0;$xy=0;$s_x=0;$s_y=0;

for($i=0;$i<$n;$i++) {
$x2 +=$x[$i]*$x[$i];
$y2 +=$y[$i]*$y[$i];
$xy +=$x[$i]*$y[$i];
$s_x+=$x[$i];
$s_y+=$y[$i];
}

return ($n*$xy-$s_x*$s_y)/(sqrt($n*$x2-$s_x*$s_x)*sqrt($n*$y2-$s_y*$s_y));
}

?>