数学建模作业.docx

1、数学建模作业 院 系： 数学学院 专 业： 信息与计算科学 年 级： 2014级 学生姓名： 王继禹 学 号： 201401050335 教师姓名： 徐霞 1、考察温度x对产量y的影响，测得下列10组数据：温度（）20253035404550556065产量（kg）13.215.116.417.117.918.719.621.222.524.3求y关于x的线性回归方程，检验回归效果是否显著，并预测x=42时产量的估值及预测区间（置信度95%）.解：(1)输入数据：x=20 25 30 35 40 45 50 55 60 65;X=ones(10,1) x;Y=13.2 15.1 16.4 17

2、.1 17.9 18.7 19.6 21.2 22.5 24.3;(2) 回归分析及检验：输入以下命令：b,bint,r,rint,stats=regress(Y,X)得结果：b = 9.1212 0.2230bint = 8.0211 10.2214 0.1985 0.2476stats = 0.9821 439.8311 0.0000 0.2333即 ， 的置信区间为8.0211，10.2214, 的置信区间为0.1985，0.2476， ，p z0=b(1)+b(2)*x得结果：z0 = 18.488所以，当x=42时，产量的估值为18.488kg及预测区间为16.3581，20.620

3、6 （置信度95%）。2、某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在曲线横坐标xi处测得纵坐标yi共11对数据如下：xi02468101214161820yi0.62.04.47.511.817.123.331.239.649.761.7求这段曲线的纵坐标y关于横坐标x的二次多项式回归方程.解：(1) 输入数据：x=0 2 4 6 8 10 12 14 16 18 20;y=0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7;(2)作二次多项式回归：p,s=polyfit(x,y,2)得结果： p =

4、0.1403 0.1971 1.0105S = R: 3x3 double df: 8 normr: 1.1097即这段曲线的纵坐标y关于横坐标x的二次多项式回归方程为(3) 预测及作图在matlab中输入的程序：x=0 2 4 6 8 10 12 14 16 18 20;y=0.6 2.0 4.4 7.5 11.8 17.1 23.3 31.2 39.6 49.7 61.7;p,s=polyfit(x,y,2)得出结果再输入：Y=polyconf(p,x,S);得出结果再输入：plot(x,y,k+,Y,r)得到试验点与回归曲线的图形（图3）。图33某校60名学生的一次考试成绩如下：9375

5、83939185848277767795948991888683968179977875676968848381756685709484838280787473767086769089716686738094797877635355(1) 计算均值，标准差，极差，偏度，峰度，画出直方图；(2) 检验分布的正态性；(3) 若检验符合正态分布，估计正态分布的参数并检验参数.解：在MATLAB中建立m文件：Untitled.m输入数据：x1=93 75 83 93 91 85 84 82 77 76 77 95 94 89 91;x2=88 86 83 96 81 79 97 78 75 67 69

6、 68 84 83 81;x3=75 66 85 70 94 84 83 82 80 78 74 73 76 70 86;x4=76 90 89 71 66 86 73 80 94 79 78 77 63 53 55;x=x1 x2 x3 x4;(1)计算均值，标准差，极差，偏度，峰度，画出直方图均值：j=mean(x)标准差：b=std(x)偏度：p=skewness(x)峰度：f=kurtosis(x)建立M文件：Untitled2.m：x1=93 75 83 93 91 85 84 82 77 76 77 95 94 89 91;x2=88 86 83 96 81 79 97 78 75

7、 67 69 68 84 83 81;x3=75 66 85 70 94 84 83 82 80 78 74 73 76 70 86;x4=76 90 89 71 66 86 73 80 94 79 78 77 63 53 55;x=x1 x2 x3 x4;j=mean(x) %b=std(x) %p=skewness(x) %f=kurtosis(x) %结果：Untitled2j = 80.1000b = 9.7106p = -0.4682f = 3.1529极差：用z表示极差。编写M文件：Untitled1.mx1=93 75 83 93 91 85 84 82 77 76 77 95

8、94 89 91;x2=88 86 83 96 81 79 97 78 75 67 69 68 84 83 81;x3=75 66 85 70 94 84 83 82 80 78 74 73 76 70 86;x4=76 90 89 71 66 86 73 80 94 79 78 77 63 53 55;X=min(x1);min(x2);min(x3);min(x4);Y=max(x1);max(x2);max(x3);max(x4);z=max(Y)-min(X)运行结果：z =44画出直方图：描绘直方图的命令：hist(data,k)；建立m文件：Untitled3.mx1=93 75

9、83 93 91 85 84 82 77 76 77 95 94 89 91;x2=88 86 83 96 81 79 97 78 75 67 69 68 84 83 81;x3=75 66 85 70 94 84 83 82 80 78 74 73 76 70 86;x4=76 90 89 71 66 86 73 80 94 79 78 77 63 53 55;x=x1 x2 x3 x4;hist(x,10)图4 频数直方图从图4可以知道，学生成绩可以大致看作近似服从正态分布。(2) 检验分布的正态性在Matlab中输入命令：x1=93 75 83 93 91 85 84 82 77 76

10、77 95 94 89 91;x2=88 86 83 96 81 79 97 78 75 67 69 68 84 83 81;x3=75 66 85 70 94 84 83 82 80 78 74 73 76 70 86;x4=76 90 89 71 66 86 73 80 94 79 78 77 63 53 55;x=x1 x2 x3 x4;normplot(x)运行结果：从图5可以看出，数据基本分布在一条直线上，故初步可以断定学生考试成绩为正态分布。图5 正态概率图(3) 若检验符合正态分布，估计正态分布的参数并检验参数在基本确定数据的分布后，就可以进行该数据的参数估计。muhat,sig

11、mahat,muci,sigmaci=normfit(x)在matlab中输入命令： x1=93 75 83 93 91 85 84 82 77 76 77 95 94 89 91;x2=88 86 83 96 81 79 97 78 75 67 69 68 84 83 81;x3=75 66 85 70 94 84 83 82 80 78 74 73 76 70 86;x4=76 90 89 71 66 86 73 80 94 79 78 77 63 53 55;x=x1 x2 x3 x4; muhat,sigmahat,muci,sigmaci=normfit(x)运行结果：muhat =

12、 80.1000sigmahat = 9.7106muci = 77.5915 82.6085sigmaci = 8.2310 11.8436 估计出学生成绩的均值为80，标准差为10，均值的0.95置信区间为77.6，82.6，标准差的0.95置信区间为8.2，11.8。已知60名学生的成绩服从正态分布，现在在方差未知的情况下，检验其均值m是否等于80.在matlab中的命令如下：h,sig,ci=ttest(x,80)程序： x1=93 75 83 93 91 85 84 82 77 76 77 95 94 89 91;x2=88 86 83 96 81 79 97 78 75 67 69

13、 68 84 83 81;x3=75 66 85 70 94 84 83 82 80 78 74 73 76 70 86;x4=76 90 89 71 66 86 73 80 94 79 78 77 63 53 55;x=x1 x2 x3 x4; h,sig,ci=ttest(x,80)结果：h = 0sig = 0.9367ci = 77.5915 82.6085说明：h =0，sig=0.9367，ci=77.5915 82.6085。检验结果（1）布尔变量h=0，表示不拒绝零假设，说明提出的假设学生成绩均值80是合理的。（2）95%的置信区间为77.6，82.6，它完全包括80，且精度很

14、高。（3）sig的值为0.9367，远超过0.5，不能拒绝零假设。所以，可以认为学生成绩的平均成绩为80.小学二（2）班班规一、 安全方面1、 每天课间不能追逐打闹。2、 中午和下午放学要结伴回家。3公路上走路要沿右边走，过马路要注意交通安全。4不能在上学路上玩耍、逗留。二、学习方面1、每天到校后，不允许在走廊玩耍打闹，要进教室读书。2、每节课铃声一响，要快速坐好，安静地等老师来上课。3、课堂上不做小动作，不与同桌说悄悄话， 认真思考，积极回答问题。4、养成学前预习、学后复习的好习惯。每天按时完成作业，保证字迹工整，卷面整洁。5、考试时做到认真审题，不交头接耳，不抄袭，独立完成答卷。三、升旗排

15、队和两操方面1、升旗时，要快速出教室排好队，做到快、静、齐，安静整齐地排队走出课室门，班长负责监督。2、上午第二节后，快速坐好，按要求做好眼保健操。3、下午预备铃声一响，在座位上做眼保健操。四、卫生方面1、每组值日生早晨7：35到校做值日。2、要求各负其责，打扫要迅速彻底，打扫完毕劳动工具要摆放整齐。3、卫生监督员（剑锋，锶妍，炜薪）要按时到岗，除负责自己的值日工作外，还要做好记录。五、 一日常规1每天学生到齐后，班长要检查红领巾。2劳动委员组织检查卫生。3、 每天负责领读的学生要督促学生学习。4、 上课前需唱一首歌，由文娱委员负责。5做好两操。6放学后，先做作业，然后帮助家长至少做一件家务事

16、。7如果有人违反班规，要到老师处说明原因。班训：坐如钟 站如松 快如风 静无声班规：课堂听讲坐如钟，精神集中认真听；排队升旗站如松，做操到位展雄风；做事迅速快如风，样样事情记得清；自习课上静无声，踏实学习不放松；个人努力进步快，团结向上集体荣；我为领巾添光彩，标兵集体记我功。加分标准序号考核项目加分值备注1单元考试满分+22单元考试85分以上+13课堂小测满分+14期中、期末考试满分+35在红领巾广播站投稿一次+26在校级活动中获奖+57作业十次全对得一颗星+38课堂上得到表扬+19班干部工作认真负责+110做好事、有利于班集体和学校的事+211进步比较明显+212连续一周该组值日卫生达标本组值日生每人加2分扣分标准序号考核项目扣分值备注1没交作业、不做晚作业-12忘带书本、学具-13迟到-14在课堂上被老师点名-25不穿校服，不戴红领巾-16吃零食、带钱、带玩具-27说脏话、打架-3请家长，写保证书8座位周围有垃圾-29课间操、眼保健操不认真做-110升旗时违反纪律-211来学校不进教室，在走廊聊天打闹-112体育课打闹说话、排队不整齐-2注：每人基本分60分起，学期末核算总分，作为学期评先依据。