没有平均值 C++ 的标准偏差

standard deviation without mean c++

本文关键字:标准 C++ 平均值      更新时间:2023-10-16

我有一个有效的程序来查找许多整数的标准偏差。但是,我要找到一种方法来获得没有平均值的标准差。

我理解公式是: std dev = sqrt [(B - A^2/N(/N]

哪里

A 是数据值的总和;

B 是数据值的平方和;

N 是数据值的数量。

但是我将如何在代码中编写它?这是我的偏差函数,但它使用平均值:

float calculateSD(int arr[])
{
float sum = 0.0, mean, standardDeviation = 0.0;
int i;
for(i = 0; i < SIZE; ++i)
{
    sum += arr[i];
}
mean = sum/SIZE;
for(i = 0; i < SIZE; ++i)
  //convert standardDeviation to float
    standardDeviation += static_cast<float>(pow(arr[i] - mean, 2));
//return standard deviation
return sqrt(standardDeviation / SIZE);
}    
#include <iostream>
#include <vector>
#include <numeric>
#include <math.h>
double stddev(std::vector<int> const& data)
{
    auto stats = std::make_pair(0.0,0.0);
    stats = std::accumulate(data.begin(), data.end(), stats,
                            [](std::pair<double,double> stats, double x) {
                                stats.first += x;
                                stats.second += x * x;
                                return stats;
                            });
    return sqrt((stats.second - pow(stats.first, 2.0) / data.size()) / data.size());
}
int main(int argc, const char *argv[])
{
    std::cout << stddev({1,1,1,1}) << std::endl;
    std::cout << stddev({1,2,1,2}) << std::endl;
    std::cout << stddev({1,10,1,10}) << std::endl;
}