c++读取数据处理皮尔逊

Pearson.h

#pragma once


#define M 10
#define N 2
#define K  2

long double calcCov(long double s[], long double t[])
{
	int i = 0, l = 0;
	long double mulE = 0, sourceE = 0, dataE = 0;

	for (i = 0, l = M; i<l; i++) {
		mulE += s[i] * t[i];
		sourceE += s[i];
		dataE += t[i];
	}
	return mulE - sourceE * dataE / l;

}


long double calcDenominator(long double source[], long double data[]) {
	int i, l;

	long double	sourceSquareAdd, sourceAdd, dataSquareAdd, dataAdd;

	sourceSquareAdd = 0;
	sourceAdd = 0;
	dataSquareAdd = 0;
	dataAdd = 0;

	for (i = 0, l = M; i<l; i++) {
		sourceSquareAdd += source[i] * source[i];
		sourceAdd += source[i];
		dataSquareAdd += data[i] * data[i];
		dataAdd += data[i];
	}
	return sqrt((sourceSquareAdd - sourceAdd * sourceAdd / l)*(dataSquareAdd - dataAdd * dataAdd / l));

}



l