SCAU 统计学 实验4

#include <iostream>

using namespace std;

int main()
{
    double a[10]={6.5,6.6,6.7,6.8,7.1,7.3,7.4,7.7,7.7,7.7};    //方式1
    double b[10]={4.2,5.4,5.8,6.2,6.7,7.7,7.7,8.5,9.3,10.0};  //方式2
    //构建第一种排队方式等待时间标准差的95%的置信区间
    double n=10;
    double A=0.05;    //显著性水平
    double averageA=0;
    double averageB=0;
    double sumA=0;
    double sumB=0;

    for(int i=0;i<10;i++){
        sumA+=a[i];
    }
    averageA=sumA/n;

    double temp=0;
    for(int i=0;i<10;i++){
        temp+=(a[i]-averageA)*(a[i]-averageA);
    }
    double Sa2=temp/(n-1);
    cout<<" 第一种排队方式等待时间方差："<<Sa2<<endl;
    cout<<endl;

    //查表
    double Xleft=19.0228;
    double Xright=2.7004;

    double left=((n-1)*Sa2)/Xleft;
    double right=((n-1)*Sa2)/Xright;

    double compareA=right-left;

    cout<<" 第一种排队方式等待时间标准差的置信区间为："<<left<<"~"<<right<<endl;

    //构建第二种排队方式等待时间标准差的95%的置信区间
    cout<<"----------------------------------------------------------------"<<endl;
    for(int i=0;i<10;i++){
        sumB+=b[i];
    }
    averageB=sumB/n;

    temp=0;
    for(int i=0;i<10;i++){
        temp+=(b[i]-averageB)*(b[i]-averageB);
    }
    double Sb2=temp/(n-1);
    cout<<" 第二种排队方式等待时间方差："<<Sb2<<endl;
    cout<<endl;

    left=((n-1)*Sb2)/Xleft;
    right=((n-1)*Sb2)/Xright;
    double compareB=right-left;

    cout<<" 第二种排队方式等待时间标准差的置信区间为："<<left<<"~"<<right<<endl;


    //比较
    //cout<<endl;
    cout<<"----------------------------------------------------------------"<<endl;
    if(compareA>compareB){
        cout<<" 第二种排队方式更好，因为它的离散程度小于第一种排队方式。"<<endl;
    }else if(compareA<compareB){
        cout<<" 第一种排队方式更好，因为它的离散程度小于第二种排队方式。"<<endl;
    }
    return 0;
}

#include <iostream>
#include <math.h>

using namespace std;

int main()
{
    double n=10;
    double a[10]={78,63,72,89,91,49,68,76,85,55};
    double b[10]={71,44,61,84,74,51,55,60,77,39};
    double c[10];
    for(int i=0;i<10;i++){
        c[i]=a[i]-b[i];
    }
    double sum=0;
    for(int i=0;i<10;i++){
        sum+=c[i];
    }
    double d=sum/n;
    double temp=0;
    for(int i=0;i<10;i++){
        temp+=(c[i]-d)*(c[i]-d);
    }
    double Sd=sqrt(temp/(n-1));
    cout<<" 两种方法平均自信心得分之差的标准差："<<Sd<<endl;
    cout<<endl;

    //查表
    double t=2.2622;
    double Float=t*Sd/sqrt(n);
    cout<<" 两种方法平均自信心得分之差的95%的置信区间为:"<<d-Float<<"~"<<d+Float<<endl;

    return 0;
}