最近写的一段c++算法中需要计算标准差,于是就造了个轮子,顺便记录下来...
看如下公式,可以推断出计算标准差分为几步:
/**
* 求平均值
*/
double average(double *x, int len)
{
double sum = 0;
for (int i = 0; i < len; i++) // 求和
sum += x[i];
return sum/len; // 得到平均值
}
/**
* 求方差
*/
double variance(double *x, int len)
{
double average = average(x, len);
for (int i = 0; i < len; i++) // 求和
sum += pow(x[i] - average, 2);
return sum/len; // 得到平均值
}
/**
* 求标准差
*/
double average(double *x, int len)
{
double variance = variance(x, len);
return sqrt(variance); // 得到标准差
}
最后,注意double *x在使用后的释放问题,遵循"谁申请谁释放"(new[]/delete[])原则
if(NULL != x)
{
delete[] x;
x = NULL;
}