资源描述：
《数学建摸实验-有氧锻炼》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、实验报告课程名称姓名、！//:学建模叶美芳年级12级日期5.19学号料紳*料1235班级数学*班实验名称回归分析一.实验目的及要求：1.掌握回归分析的基本理论2.会运用回归分析相关理论进行编程和解决实际问题二.实验内容:1.背景：由于有氧锻炼中，人体的耗氧能力是衡量人身体健康状况的重要指标，而耗氧能力(人体单位重量单位时间内最大的耗氧量)在日常生活中难以直接测量，故本文建立相关数学模型,采用机理分析和回归分析以及数据拟合相结合的办法，以期望使普通大众能够方便快键地对自身身体状况有更为直接的了解，也能更加合理科学
2、地锻炼。2•题目：习题7:在有氧锻炼中人的耗氧能力y(mL/(min*kg))是衡量身体状况的重要指标，它可能与以下因素有关:年龄七，体重X2(kg)，1500m跑的时间七⑽⑴,静止时心跳速度义(次/min)，跑步后心速又5(次/min),对24名40至57岁的志愿者进行了测试，结果如下表1.1(节选)，试建立耗氧能力y与诸因素的之间的回归模型。表1.1序号1234♦•參21222324Y44.645.354.359.6••參39.446.145.454.7XI44404442•••57545250X289.5
3、75.185.868.2•••73.479.476.370.9X36.826.045.194.9蠢•參7.586.75.785.35X462624540♦•參58624848X5178185156166••參174156164146(1)若中许选择1个变量，最好的模型是什么？(2)若中只许选择2个变S，最好的模型是什么？(3)若不限制变量的个数，最好的模型是什么？(4)对最终模型观察残查，有无异常点，若有，剔除后如何?1.做法：本题不同小问需要建立不同模型，由于专业知识所限，并且提供的数据较少，难以做出精确符合
4、现实情况的模型，因此这里川最简单的线性回归法进行拟和模型基本形式如下:y=P（）+PIxI+...++'h%WC事实上，中的项（高次项和交互项）对于本题目来讲意义不大，因为所给定的5个自变量和因变量之间关系比较模糊，几个变量彼此之间的联系也很难说清，因此用自变量的一次线性拟和就足以适应本题的要求。但作为练:还是将每种回归方法都使用到Y,可以用于参考。具体采用的各个模型将在下面单独说明，这里不再重复。2.程序由于本题需要建立多组模型，并且要在不断的调试中发现最合理的，很多命令都要在这个过程中不断使用，这里仅仅给出
5、使用的最基本的命令。数据clearA二[...];%数据矩阵，略n二24;y=A（2,：）；%提取各个数据x1=A（3，:）;x2=A（4,:）;x3=A（5，:）;x4=A（6，:）;x5=A（7，:）;绘制散点阁（大致判断影响怙况）fori=l:5subplot（2,3，i），plot（A（i+2，:）,y，’+’），gridpauseendpause单参数hl归（第一M）X=[ones（n,l）,x4,J;%这里检验的是自变量x4，实际操作吋要分别检验xl~x5
6、b.bint，r，rini，s
7、=regr
8、ess（y’，X）;。/。回归分析程序（a=0.05）b,bint,s,%输出回归系数估计值、置信区间、以及统计ftrcoplot（r，rint）%残差图Polytool（x3y',2）%檢验一71I多项式回归的结果，输出交互式画而X5=[xl',x2x3,,x4,,x5'
9、;stepwise(X5，yr):XX=[x3xl*];rstool(XX，y./linear’)双参数回归（第二问）：用逐步回归法找出敁合理的两个变呈%利川输出的交互式画面，可以选fli最佳的两个变眾%当得到了最佳的两个变量后（这里
10、假没是x3xl）%检验二元情况下的交互项和高次项全部参数回归（第三问）：X5=[xl',x2,,x3,,x4x5,];%仍然川逐步回归法找出最合理的组合力‘式stepwise（X5，y’）笫五问要求对残差进行分析，并且剔除异常点，可以在该M得到最终模型后，采用regress得到的残差值和罝信区间并根据其绘制残差阁，然后再进行剔除操作重新检验。5.运行结果及分析散点图(2)656055鮝--•-5044<4404+♦67(3)•▲4*♦4▲±▼♦令404550556D6570(4)556C疇■5C‘54-44
11、-+▲44-140150160170180190(5)从左上到右下的顺序为xl~x5.可以由点的分布大致看出，除了x3自变量呈现比较明显的负相关趋势以外，对于其他的各个自变量都难以直接观测出其对于因变量的影响。根据这种结果，可以假设自变量x3（1500ni跑后心速）最直接的与锻炼耗氧能力相关，下面通过对各个自变量的单参数回归进行检验。单参数回归被检对象P0P1P1置信区间R"2FPs"
