北京科技大学概率论与数理统计上机报告3.doc
《北京科技大学概率论与数理统计上机报告3.doc》由会员分享,可在线阅读,更多相关《北京科技大学概率论与数理统计上机报告3.doc(24页珍藏版)》请在咨信网上搜索。
1、班级:信计1502组号:35 组员:何芝芝、吕瑞杰、陈彦睿专业: 信息与计算科学 班级: 信计1502(35组) 学生姓名: 吕瑞杰 陈炎睿 何芝芝 指导教师: 张志刚 完成时间: 2024年9月3日 概率论与数理统计第三次上机报告 - 1 -Matlab 概率论与数理统计上机练习(3)五、 假设检验【例】(离散型分布检验)某工厂近五年发生了63起事故,按星期几可以分为9 10 11 8 13 12,问该厂发生的事故数是有与星期几有关?clear allmi=9 10 11 8 13 12; % 周一到周六的事故数n=sum(mi); % 总的事故数r=0; % 总体中没有未知参数k=leng
2、th(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-
3、1);%右侧概率if kai2la2 xzx=*;elseif kai2la1 xzx=*;else xzx=-;endx=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 onif kai2=60);se6=length(find(se=60);%及格人数sy7=length(find(sy=80);se7=length(find(se=80);%优秀人数sy8=sy6/sy1;se8=se6/se1;%及格率sy9=sy7/sy1;se9=se7/se1;%优秀率fprintf(
4、t人数t 平均分t 最小值t 最大值 t极差tt标准差tt及格人数 及格率t优秀人数 优秀率n);fprintf( -n);fprintf(数学 %4dt%10.4ft%4dtt%4dtt%4dtt%10.4ft%4dt %10.4ft%4dt%10.4fn,sy1,sy4,sy2min,sy2max,sy3,sy5,sy6,sy8,sy7,sy9)fprintf(信计 %4dt%10.4ft%4dtt%4dtt%4dtt%10.4ft%4dt %10.4ft%4dt%10.4fn,se1,se4,se2min,se2max,se3,se5,sy6,sy8,se7,se9)fprintf(n)
5、;%方法一fprintf(检验数学和信计的方差是否相等n);h1,p1,varci1,stats1=vartest2(sy,se,alpha,both);if(h1=0) disp(结果:方差相等);else disp(结果:方差不相等);endfprintf(n);% %方法二% F=sy52/se52;%统计量F,满足F分布% alpha=0.05; %取显著水平为0.05% Fla1=finv(alpha/2,sy1-1,se1-1);Fla2=finv(1-alpha/2,sy1-1,se1-1);%求F的临界值% if (FFla1 & FFla2)% MM=数学分析1和数学分析2的
6、方差无显著差异;% else% MM=数学分析1和数学分析2的方差有显著差异;% end% fprintf(检验数学分析1和数学分析2的方差是否相等n);% fprintf(统计量F的值ttt显著性水平tt临界值ttttt检验结果n);% fprintf( %.4ftttt%.4fttt%.4fttt%15sn,F,alpha,Fla1,MM);% fprintf(nn);%方法一fprintf(检验数学和信计的平均分是否相等n);h2,p2,muci2,stats2=ttest2(sy,se,alpha,both);if(h2=0) disp(结果:平均分相等);else disp(结果:平
7、均分不相等);endfprintf(n);%方法二% %方法三% sw=(sy1-1)*sy52+(se1-1)*se52)/(sy1+se1-2);% T=(sy4-se4)/sw/sqrt(1/sy1+1/se1);%统计量T,满T分布% Tla1=tinv(alpha/2,sy1+se1-2);Tla2=tinv(1-alpha/2,sy1+se1-2);%求出T的临界值% if (abs(T)Ua) disp(优秀率无显著差异);else disp(优秀率有显著差异);endp=(sy6+se6)/(sy1+se1);U=(sy8-se8)/sqrt(sy8+se8)*p*(1-p);
8、Ua=norminv(1-alpha/2);if(abs(U)Ua) disp(及格率无显著差异);else disp(及格率有显著差异);endh3,p3,kstat3,critval3=lillietest(sy,alpha);if(h3=1) disp(数学不是正态分布)else disp(数学是正态分布)endh4,p4,kstat4,critval4=lillietest(se,alpha);if(h4=1) disp(信计不是正态分布)else disp(信计是正态分布)end%hist(sy)%直方图%h5,p5,stats5=chi2gof(sy)%可以检验分布%cdf=sy,
9、normcdf(sy,sy4,sy5)%h5,p5,ksstat,cv5=kstest(sy,cdf)% a=0:1:100;% a=a;% CDF=a,cdf(a,sy4,sy5);% h = kstest(sy,CDF,0.05);S=656568817476686982777466737277606281666876607480906960636867696260606067607767606060717260666186646060897374434061956970626663627874605076626584706983734371707371697460616070747848
10、936461797153606052636061656278606560858589696660;h,p,jbstat,critval=jbtest(S,alpha);if(h=0) disp(服从正态分布);else disp(不服从正态分布);endsavg=mean(S);svar=var(S);x=20:130;y=normpdf(x,savg,sqrt(svar);d=5;a=20:d:130;pdf=hist(S,a)./length(S)./d;plot(x,y,r);hold onscatter(a,pdf,filled);hold off输出: lx3_1_lrj_41521
11、335人数 平均分 最小值 最大值 极差标准差及格人数 及格率优秀人数 优秀率 -数学 54 70.6667 40 95 55 10.0885 53 0.9815 9 0.1667信计 60 65.5167 43 89 46 9.2313 53 0.9815 5 0.0833检验数学和信计的方差是否相等结果:方差相等检验数学和信计的平均分是否相等结果:平均分不相等优秀率有显著差异及格率有显著差异数学不是正态分布Warning: P is less than the smallest tabulated value, returning 0.001. In lillietest (line 20
12、6) In lx3_1_lrj_41521335 (line 99) 信计不是正态分布服从正态分布六、 方差分析【例1】(单因素方差分析)考虑温度对某化工产品得率的影响,选择五种不同温度进行试验,每一温度各做三次试验。方法一:自编程序clear allX=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,SET=0; ST=0; SA=0; SE=0;for i=1:a Ti
13、=0; for j=1:ni(i) T=T+X(i,j); Ti=Ti+X(i,j); ST=ST+X(i,j)2; end SA=SA+Ti2/ni(i);endST=ST-T2/n; % 总偏差平方和SA=SA-T2/n; % 效应平方和SE=ST-SA; % 误差平方和 F=(SA/(a-1)/(SE/(n-a); % F比alpha1=0.05; % 显著性水平la1=finv(1-alpha1,a-1,n-a); %由F分布的累积概率,求临界值,PFla=1-alphaalpha2=0.01; % 显著性水平la2=finv(1-alpha2,a-1,n-a); %由F分布的累积概率
14、,求临界值,PFla=1-alphap=1-fcdf(F,a-1,n-a); %计算F比值做为临界点的右侧概率 p=1-PXla2 xzx=*;elseif Fla1 xzx=*;else xzx=-;endfprintf( 来源tt平方和tt自由度tt均方和ttF比tt显著性n);fprintf( 效应Att%.2ftt%4dtt%.2ftt%.2ftt%.4fn,SA,a-1,SA/(a-1),F,p);fprintf( 误差tt%.2ftt%4dtt%.2ftttt%4sn,SE,n-a,SE/(n-a),xzx);fprintf( 总和tt%.2ftt%4dtt临界值=%.2f(%.2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 北京科技大学 概率论 数理统计 上机 报告
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。