看公式,可以推断出计算标准差分为几步:
- 计算平均值u=(x1+x2+...+xn)/n
- 计算方差s²=((x1-u)^2 +(x2-u)^2 +...+(xn-u)^2)/n
- 计算标准差σ=sqrt(s²)
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 sum = 0;
double avg = average(x, len);
for (int i = 0; i < len; i++) // 求和
sum += pow(x[i] - avg, 2);
return sum / len; // 得到方差
}
double standardDev(double *x, int len)
{
double var = variance(x, len);
return sqrt(var); // 得到标准差
}