清华大学数学实验报告.docVIP

PAGE  PAGE 15 清华大学数学实验报告 篇一：清华大学数学实验报告10 数学实验报告10回归分析 一、实验目的 1.掌握数据的参数估计、假设检验的基本原理、算法，及用MATLAB实现的方法； 2.练习用这些方法解决实际问题。 二、实验内容 题目1 【问题描述】 电影院调查电视广告费用和报纸广告费用对每周收入的影响，得到下面的数据，建立回归模型并进行检验，诊断异常点的存在并进行处理。 【分析与解】 初步判断可以发现电视台每周的收入与电视、电影广告费用具有线性关系，因而可以进行线性拟合。设电视广告费用为x1，电影广告费用为x2，每周收入为y，则有y=??0+ ??1x1+ ??2x2。 为绘制残差和置信区间，在Matlab中编程如下： y = [96 90 95 92 95 95 94 94]; x1 = [1.5 2 1.5 2.5 3.3 2.3 4.2 2.5]; x2 = [5 2 4 2.5 3 3.5 2.5 3]; n = 8; m = 2; X = [ones(n,1),x1#39;,x2#39;]; [b,bint,r,rint,s] = regress(y#39;,X); b,bint,s, rcoplot(r,rint) 运行以上代码，得到如下图像： 以及相应的输出结果： b = 83.2116 1.2985 2.3372 bint = 78.805887.6174 0.4007 2.1962 1.4860 3.1883 s = 0.908924.9408 0.0025 0.4897 通过分析图像可知，数据1为异常点， 将其舍去再进行上述分析， 得到如下结果： 以及相应的输出结果： b = 81.4881 1.2877 2.9766 bint = 78.787884.1883 0.7964 1.7790 2.3281 3.6250 s = 0.976884.3842 0.0005 0.1257 通过分析上述结果可知，R2=0.9768，因而所得到的模型精度比较高。同时，所得到的线性拟合的结果为y=81.4881+1.2877x1+2.9766x2。 若采用二次模型，则有y=β0+β1x1+β2x2+β3x12+β4x1x2+β5x22。 通过除 去异常点，再次进行回归分析可以得到： b = 8.5413 -3.0821 3.8868 9.3397 2.8304 -4.7488 同时，可以得到剩余方差为s2 = 0.0200。该结果说明二次回归模型的准确度要高于线性规划的模型。 【结论】电视台的收入与电视广告广告费用和电影广告费用可以用 线性规划模型来描述，其关系为y=81.4881+1.2877x1+2.9766x2，也可以采用多项式模型来描述，其关系为y=8.5413-3.0821x1+3.8868x2 +9.3397x12+2.8304x1x2-4.7488x22并且二次回归模型的准确度要高于线性规划。 题目2 【问题描述】 汽车销售商认为汽车销售量与汽油价格、贷款利率有关，两种类型的汽车(普通型和豪华型)18个月的调查资料如下表所示，其中y1时普通型汽车销量(千辆)，其中y2时豪华型汽车销量(千辆)， x1是汽油价格(元/gal)，x2是贷款利率(%)。 (1)对普通型和豪华型汽车分别建立如下模型： 111222 ??1=??0+??1??1+??2??2，??2=??0+??1??1+??2??2 给出β的估计值和置信区间，决定系数R2，F值及剩余方差等。 (2)用??3=0，1表示汽车类型，建立统一模型： ??=??0+??1??1+??2??2+??3??3 给出??的估计值和置信区间，决定系数R2，F值和剩余方差等，以x3=0，1代入统一模型，将结果与(1)的两个模型的结果比较，解释两者的区别。 (3)对统一模型就每种类型汽车分别作??1和??2与残差的散点图，有什么现象，说明模型有何缺陷? (4)对统一模型增加二次项和交互项，考察结果有什么改进。 【分析与解】 (1)为分析模型中各个变量的相互关系，画出y1与x1、x2，y2与x1、x2的关系图像，在Matlab中编程如下： clear clc x1 = [1.89,1.94,1.95,1.82,1.85,1.78,1.76,1.76,1.75,1.74, 1.70,1.70,1.68,1.60,1.61,1.64,1.67,1.68]; 篇二：清华大学数学实验报告1 数学实验报告1 插值与数值积分 电13 苗键强 2011010645 一、实验目的 1. 掌握用 MATLAB 计算拉格朗日、分段线性、三次样条三种插值的方法，改变结点的数目，对三种插值结果进行初步分析。 2. 掌握用 MATL