协方差及计算代码

 1 void cox(int x[],int y[],int n,double *p)
 2 {
 3    int i; 
 4    double averx,avery,sumx=0,sumy=0;
 5    double sum=0;
 6    for(i=0;i<n;i++)
 7    {
 8       sumx+=x[i];
 9       sumy+=y[i];
10    }
11    averx=sumx/n;
12    avery=sumy/n;
13    for(i=0;i<n;i++)
14      sum+=(x[i]-averx)*(y[i]-avery);
15    *p=sum/n;//原文是除以(n-1)
16 }

原文在最后除以(n-1)，把数组的下标当作数组的长度了。

 

参考： 协方差计算的c语言代码