标准偏差计算输出不正确

时间:2014-04-14 17:41:53

标签: c++

我正在尝试计算数组的标准偏差但是我的答案返回为0.我认为问题源于“计数”搞砸了。我收到数据的数组只是四个数字1,4,6,7。代码输出的答案为0,但这是不正确的。任何帮助将非常感激。

#include <iostream> 
#include <iomanip>
#include <fstream>
#include <cmath>
#include <string>

using namespace std;
double mean(double *mydata, double N);
double standard_dev(double *mydata, double m, int N);

int main()
{
    int N(0);
    char filename [100];
    double m, stdev;
    double next;

    int count=0;
    cout<<"Enter name of file: ";
    cin>>filename;

    ifstream myfile;
    myfile.open(filename);
    while(myfile>>next)
    {
        count++;
    }

    N=count;
    double *mydata;
    mydata=new double[N];

    for(int i=0; i<N; i++)
    {
        myfile>>mydata[i];
    }
    m =  mean(mydata, N);
    stdev = standard_dev(mydata, m, N);
    cout<<"The standard deviation is:" <<stdev<<endl;

    myfile.close();
    delete[] mydata;
    return 0;
}

double mean(double *mydata, double N)  
{
    double sum(0), m;
    for(int i=0; i<N; i++)
    {
        sum +=mydata[i];
    }
    m=(sum/(double)N);
    return (m);
}
double standard_dev(double *mydata, double m, int N)
{
    double *mydata2= new double [N];
    for(int i=0; i<N; i++)
    {
        mydata2[i]= pow((mydata[i]-m),2);
    }
    double sum(0), S, X;
    for(int i=0; i<N; i++)
    {
        sum+=mydata2[i];
    }
    X=sum/N;
    S=sqrt(X);
    return (S);
}

包含4个数字的数组的代码只是: 1 4 6 7 它们在编码方面是垂直的。

1 个答案:

答案 0 :(得分:2)

您正在读错数据。首先,您将所有数据都读入next以增加count。然后,您尝试使用for循环读取不存在的数据。我建议您使用std::vector来避免动态内存分配问题;

vector<double> mydata;
while(myfile>>next)
{
    mydata.push_back(next);
}

另一个最简单的解决方案就是将元素数量放在文件中。例如:

文件:

  

4 1 4 6 7

从文件中读取N并使用for循环