python22版从一批保险丝中抽取10根，测试其融化时间，得到的数据如下（单位：min）:42,65...设这批保险丝的融化时间负荷正态分布，检验总体方差是否为12**2，答案。

目录

解题思路

一：题目

二：思路

开发环境

代码及效果

一：代码

二：效果

---

解题思路

一：题目

3从一批保险丝中抽取10根，测试其融化时间，得到的数据如下（单位：min）:
42,65,75,78,71,59,57,68,55,54.设这批保险丝的融化时间负荷正态分布，检验总体方差是否为12**2（a=0.05）

 

二：思路

问题属于数据统计分析。单个总体方差的假设检验问题。具体思路1·由样本数据计算出样本标准差
2做出假设，3计算检验统计量。4算出双侧临界值。

符号说明：

n为样本容量=10,
    S为保险丝的样本标准差值，
    ta/2(n-1)为自由度为（n-1）的t分布的上侧a/2的分位数。
    x^2检验统计量的值。

模型假设：

n=10,样本容量为10，
    a的值不变化，
    模型的xi是取值来自X-N(u,σ^2)的样本，始终满足题目的正态分布
    Ho:u=12*12;
    H0:u≠12*12
    ta/2(n-1)可以通过查询数据获得

目标方程：

 

开发环境

python 3.9  statsmodels 库

代码及效果

一：代码

代码示例：

from scipy.stats import t
import numpy as np
from scipy.stats import chi2
a=np.array([42,65,75,78,71,59,57,68,55,54])
v=np.var(a,ddof=1)
n=len(a)#样本容量
chiT=(n-1)*v/12**2#计算统计值
#小范围的值，大范围的值。
ca1=chi2.ppf(0.025,n-1);ca2=chi2.ppf(0.975,n-1)
print('统计量的值：',chiT)
print('s**2:',v)
if v<ca1 or v>ca2  :
 print('落入否定域,拒绝假设:')
else:
 print('接受假设:')

二：效果

效果

n=10 ,     

统计量的值： 7.61388888888889

s**2: 121.82222222222224

落入否定域,拒绝假设:

即显然x^2统计量的值，7.613888 > 2.7003894999803584，7.613888>19.02276779864163

所以拒绝原假设，

认为总体方差确实不为12*12