matlab数理统计最新完整版本.pptVIP

方差分析方差分析例9一火箭使用了4种燃料，3种推进器作射程试验，每种燃料与每种推进器的组合各发射火箭2次，得到结果如下：B1B2B3A158.2,52.656.2,41.265.3,60.8A249.1,42.854.1,50.551.6,48.4A360.1,58.370.9,73.239.2,40.7A475.8,71.558.2,51.048.7,41.4试在水平0.05下，检验不同燃料（因素A）、不同推进器（因素B）下的射程是否有显著差异？交互作用是否显著？方差分析解编写程序如下：clc,clearx0=[58.2,52.656.2,41.265.3,60.849.1,42.854.1,50.551.6,48.460.1,58.370.9,73.239.2,40.775.8,71.558.2,51.048.7,41.4];x1=x0(:,1:2:5);x2=x0(:,2:2:6);fori=1:4x(2*i-1,:)=x1(i,:);x(2*i,:)=x2(i,:);endp=anova2(x,2)求得p=0.00350.02600.001，表明各试验均值相等的概率都为小概率，故可拒绝均值相等假设。即认为不同燃料（因素A）、不同推进器（因素B）下的射程有显著差异，交互作用也是显著的。第五讲

MATLAB数理统计引言数理统计研究的对象是受随机因素影响的数据数据样本少则几个，多则成千上万，人们希望能用少数几个包含其最多相关信息的数值来体现数据样本总体的规律。面对一批数据如何进行描述与分析，需要掌握参数估计和假设检验这两个数理统计的最基本方法。我们将用MATLAB的统计工具箱(StatisticsToolbox)来实现数据的统计描述和分析。频数表和直方图一组数据（样本）往往是杂乱无章的，作出它的频数表和直方图，可以看作是对这组数据的一个初步整理和直观描述。将数据的取值范围划分为若干个区间，然后统计这组数据在每个区间中出现的次数，称为频数，由此得到一个频数表。以数据的取值为横坐标，频数为纵坐标，画出一个阶梯形的图，称为直方图，或频数分布图。若样本容量不大，能够手工作出频数表和直方图，当样本容量较大时则可以借助MATLAB这样的软件了。例1学生的身高和体重学校随机抽取100名学生，测量他们的身高和体重，所得数据如表身高体重17275171621666216055155571735816655170631675317360身高体重16955168671686517567176641685016149169631716117864身高体重16964165521645917374172691695217357173611667016357身高体重17165169621705817264169581677217576164591666316954身高体重16747168651656416857176571705715851165621725316966数据输入两种方法一种是在交互环境中直接输入如果在统计中数据量比较大，这样作不太方便；另一种办法是先把数据写入一个纯文本数据文件data.txt中，数据列之间用空格键或Tab键分割，该数据文件data.txt存放在matlab\work子目录下，在MATLAB中用load命令读入数据具体作法是：loaddata.txt这样在内存中建立了一个变量data，它是一个数据矩阵。作频数表及直方图用hist命令实现，其用法是：[N,X]=hist(Y,M)数组（行、列均可）Y的频数表它将区间[min(Y),max(Y)]等分为M份（缺省时M设定为10），N返回M个小区间的频数，X返回M个小区间的中点。hist(Y,M)数组Y的直方图。例1的程序loaddata.txt;high=data(:,1:2