算法总结.doc

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
算法总结

单因素方差分析: 1)均衡数据: p=anoval(x) 返回值p 是一个概率,当 p α 时接受H0,x为m×r的数据矩阵,x 的每一列是一个水平的数据(这里各个水平上的样本容量ni= m)。另外,还输出一个方差表和一个Box 图。 2)非均衡数据 p=anova1(x,group) x 为向量,从第 1 组到第 r 组数据依次排列;group 为与 x 同长度的向量,标志 x 中数据的组别(在与 x 第i 组数据相对应的位置处输入整数i(i=1,2,……,r) 双因素方差分析 p=anova2(x,reps) 其中 x 不同列的数据表示单一因素的变化情况,不同行中的数据表示另一因素的变化情 况。如果每种行—列对(“单元”)有不止一个的观测值,则用参数 reps 来表明每个“单 元”多个观测值的不同标号,即 reps 给出重复试验的次数t。 例 3 一种火箭使用了四种燃料、三种推进器,进行射程试验,对于每种燃料与每 种推进器的组合作一次试验,得到试验数据如表 8。问各种燃料之间及各种推进器之间 有无显著差异? 解 记燃料为因素A,它有4个水平,水平效应为 αi ,i=1,2,3,4。推进器为因素B, 它有3个水平,水平效应为βj , j=1,2,3 。我们在显著性水平 α=0.05 下检验 编写如下的Matlab程序: x=[58.2 56.2 65.3 49.1 54.1 51.6 60.1 70.9 39.2 75.8 58.2 48.7]; [p,t,st]=anova2(x) 求得p=0.4491 0.7387,表明各种燃料和各种推进器之间的差异对于火箭射 程无显著影响。 例 4 一火箭使用了 4 种燃料,3 种推进器作射程试验,每种燃料与每种推进器的 组合各发射火箭 2 次,得到如表 9 结果。 试在水平 0.05 下,检验不同燃料(因素A)、不同推进器(因素B)下的射程是 否有显著差异?交互作用是否显著? 解 编写程序如下: clc,clear x0=[58.2,52.6 56.2,41.2 65.3,60.8 49.1,42.8 54.1,50.5 51.6,48.4 60.1,58.3 70.9,73.2 39.2,40.7 75.8,71.5 58.2,51.0 48.7,41.4]; x1=x0(:,1:2:5);x2=x0(:,2:2:6); for i=1:4 x(2*i-1,:)=x1(i,:); x(2*i,:)=x2(i,:); end [p,t,st]=anova2(x,2) 求得 p=0.0035 0.0260 0.0001,表明各试验均值相等的概率都为小概率,故 可拒绝均值相等假设。即认为不同燃料(因素A)、不同推进器(因素B)下的射程有 显著差异,交互作用也是显著的。 回归分析: 多元线性回归: [b,bint,r,rint,stats]=regress(Y,X,alpha) alpha为显著性水平(缺省时设定为0.05),b,bint为回归系数估计值和它们的置信区间,r,rint 为残差(向量)及其置信区间,stats 是用于检验回归模型的统计量,有四个数值,第一个是R^2,第二个是F,第三个是与F对应的概率p, pα 拒绝 H0,回归模型成立,第四个是残差的方差s ^2 一元多项式回归:(类似插值) [p,s]=polyfit(x0,y0,n); x0=17:2:29;x0=[x0,x0]; y0=[20.48 25.13 26.15 30.0 26.1 20.3 19.35 24.35 28.11 26.3 31.4 26.92 25.7 21.3]; [p,s]=polyfit(x0,y0,2); p 得到 p =-0.2003 8.9782 -72.2150 即a2=-0.2003 a1=8.9782 a0= -72.2150 polytool(x0,y0,2) 可以得到一个交互式画面 多元二项式回归: rstool(x,y,model,alpha) 也会产生交互式画面 linear(线性) purequadratic(纯二次): interaction(交叉): quadratic(完全二次): 编程如下: x1=[120 140 190 130 155 175 125 145 180 150]; x2=[100 110 90 150 210 150 250 270 300 250]; y=[102 100 120 77 46 93 26 69 65 85]; x=[x1 x2]; rstool(x,y,purequadratic) 图的左下方有两个下拉式菜单,一个菜单Export用

文档评论(0)

ranfand + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档