概率论与数理统计假设检验..ppt

% 定义甲机床对应的样本观测值向量 x = [20.1, 20.0, 19.3, 20.6, 20.2, 19.9, 20.0, 19.9, 19.1, 19.9]; % 定义乙机床对应的样本观测值向量 y = [18.6, 19.1, 20.0, 20.0, 20.0, 19.7, 19.9, 19.6, 20.2]; alpha = 0.05; % 显著性水平为0.05 tail = both; % 尾部类型为双侧 vartype = equal; % 方差类型为等方差 % 调用ttest2函数作两个正态总体均值的比较检验， % 返回变量h，检验的p值，均值差的置信区间muci，结构体变量stats [h,p,muci,stats] = ttest2(x,y,alpha,tail,vartype) 例6、首先用产生正态分布随机数命令生成两组均值μ分别为1 和2, 标准差σ均为4 的正态分布样本xx 和yy, 用双样本均值t 检验函数ttest2 来检验两个样本的均值是否相等. xx=normrnd(1,4,[1 100]); %生成μ=1, σ=4 的一组正态随机数. yy=normrnd(2,4,[1 100]); %生成μ=2, σ=4 的一组正态随机数. [h,sig,ci]=ttest2(xx,yy,0.05) 例7 设有甲、乙两种零件彼此可以代用，但乙零件比家零件制造简单，造价低，经过试验获得它们的抗压强度数据如下表(单位：kg/cm2) 甲种零件 88 87 92 90 91 乙种零件 89 89 90 84 88 87 已知甲、乙两种零件的抗压强度分别服从正态总体 和 ，问能否保证抗压强度质量下，用乙种零件代替甲种零件？ 解：根据题意构造假设： Matlab求解： x=[88 87 92 90 91]; y=[89 89 90 84 88 87]; [H,P,CI]=ttest2(x,y,0.05,-1) 输出： H = 0 P = 0.9000 CI = -Inf 4.1077 三、总体均值未知时的单个正态总体方差的卡方检验 调用格式： H = vartest(X,V) H = vartest(X,V,alpha) H = vartest(X,V,alpha,tail) [H,P] = vartest(...) [H,P,CI] = vartest(...) [H,P,CI,STATS] = vartest(...) [...] = vartest(X,V,alpha,tail,dim) vartest函数 % 定义样本观测值向量 x = [49.4 50.5 50.7 51.7 49.8 47.9 49.2 51.4 48.9]; var0 = 1.5; % 原假设中的常数 alpha = 0.05; % 显著性水平为0.05 tail = both; % 尾部类型为双侧 % 调用vartest函数作单个正态总体方差的双侧检验， % 返回变量h，检验的p值，方差的置信区间varci，结构体变量stats [h,p,varci,stats] = vartest(x,var0,alpha,tail) 四、总体均值未知时的两个正态总体方差的比较 F 检验 调用格式： H = vartest2(X,Y) H = vartest2(X,Y,alpha) H = vartest2(X,Y,alpha,tail) [H,P] = vartest2(...) [H,P,CI] = vartest2(...) [H,P,CI,STATS] = vartest2(...) [...] = vartest2(X,Y,alpha,tail,dim) vartest2函数 % 定义甲机床对应的样本观测值向量 x = [20.1, 20.0, 19.3, 20.6, 20.2, 19.9, 20.0, 19.9, 19.1, 19.9]; % 定义乙机床对应的样本观测值向量 y = [18.6, 19.1, 20.0, 20.0, 20.0, 19.7, 19.9, 19.6, 20.2]; alpha = 0.05; % 显著性水平为0.05 tail = both; % 尾部类型为双侧 % 调用vartest2函数作两个正态总体方差的比较检验， % 返回变量h，检验的p值，方差之比的置信区间varci，结构体变量stat