概率统计上机练习(给学生3)2015.docVIP

Matlab 概率论与数理统计上机练习（3） 假设检验 【例】（离散型分布检验） 某工厂近五年发生了63起事故，按星期几可以分为[9 10 11 8 13 12]，问该厂发生的事故数是有与星期几有关？ clear all mi=[9 10 11 8 13 12]; % 周一到周六的事故数 n=sum(mi); % 总的事故数 r=0; % 总体中没有未知参数 k=length(mi); % 天数 pii=1/6; % 事故的概率 kai2=0; kai2=sum((mi-n*pii).^2)./(n*pii); % k2统计量的值 alpha1=0.05; % 显著性水平 alpha2=0.01; % 显著性水平 alpha3=0.001; % 显著性水平 la1=chi2inv(1-alpha1,k-r-1); % kai2分布的累计概率，即临界值 la2=chi2inv(1-alpha2,k-r-1); % kai2分布的累计概率，即临界值 la3=chi2inv(1-alpha3,k-r-1); % kai2分布的累计概率，即临界值 pz=1-chi2cdf(kai2,k-r-1);%右侧概率 if kai2la2 xzx=**; elseif kai2la1 xzx=*; else xzx=-; end x=0:0.1:la3; y=chi2pdf(x,k-r-1); plot(x,y); x1=kai2:0.1:la3; y1=chi2pdf(x1,k-r-1); hold on if kai2la3 fill([kai2,x1,la3],[0,y1,0],m) end fprintf( ---------------------------------------------------------------------------\n); fprintf( 样本数\t\t区间数\t\t未知参数\t\t自由度\t\t开方和\t\t 右侧概率\t\t显著性\n); fprintf( ---------------------------------------------------------------------------\n); fprintf( %4d\t\t%4d\t\t%4d\t\t%4d\t\t%.6f\t\t%.4f\t\t%4s\n,n,k,r,k-r-1,kai2,pz,xzx); fprintf( ---------------------------------------------------------------------------\n); fprintf(\n\n); hold off 离散型分布正态性检验 ------------------------------------------------------------------------------------- 样本数 区间数 未知参数 自由度 开方和 临界值 右侧概率 显著性 ------------------------------------------------------------------------------------- 63 6 0 5 1.6667 11.0705 0.8931 - ------------------------------------------------------------------------------------- 【练习3.1】（基本计算，两个正态总体的假设检验，检验水平），对数学分析I 求课程中“专业（数学、信计）””的考试人数、平均分、最小值、最大值、极差、标准差、及格人数、及格率、优良人数（大于等于80）、优良率；写出标准差的计算公式。 对“专业（数学、信计）”，检验方差、平均分是否相等。 对“专业（数学、信计）”，检验及格率、优秀率是否相等。 对“全体成绩”的分布进行检验，首先估计期望和方差，画出正态分布的密度函数曲线以及样本密度散点，对假设的正态分布进行检验。 方差分析 【例1】（单因素方差分析） 考虑温度对某化工产品得率的影响，选择五种不同温度进行试验，每一温度各做三次试验。 方法一：自编程序 clear all X=[90,92,88;97,93,92;96,96,93;84,83,88;84,86,82]; a=5; ni=[3,3,3,3,3]; %每个因素的样本数 n=sum(ni); %样本总数 %T=sum(sum(X)); %先求每列的和，再求总和 %求所有样本的和T，平方和，及ST,SA,SE T=0; S