内江师范学院数型实验报告册.docVIP

内江师范学院 数 学 模 型 实 验 报 告 册 编制 数学建模组 审定 牟廉明 专业： 数学与应用数学 班级： 级 班 学号： 姓名： 数学与信息科学学院 2014年月 实验名称： （实验） 指导教师： 实验时数： 4 实验设备：安装了VC++、mathematic、m的计算机实验日期：2014年 月 日 实验地点： 第五教学楼北02 实验目的： 掌握插值与拟合的原理，熟悉插值与拟合的软件实现。 实验准备： 在开始本实验之前，请回顾教科书的相关内容； 需要一台准备安装Windows XP Professional操作系统和装有VC++6.0的计算机。 实验内容及要求 下表给出了某工厂产品的生产批量与单位成本（元）的数据。从散点图可以明显地发现，生产批量在500以内时，单位成本对生产批量服从一种线性关系，生产批量超过500时服从另一种线性关系，此时单位成本明显下降。 生产批量 650 340 400 800 300 600 单位成本 2.48 4.45 4.52 1.38 4.65 2.96 生产批量 720 480 440 540 750 单位成本 2.18 4.04 4.20 3.10 1.50 要求：构造合适的模型全面地描述生产批量与单位成本的关系实验过程：x在500以内时，建立模型（1）；x超过500时，建立模型（2）。然后综合模型（1）和（2）建立回归模型。 2.符号说明： 符号 意义 符号 意义 y 单位成本 x1 生产批量大于500的部分 回归方程系数 x 生产批量 回归方程系数 x2 生产批量小于500的部分 回归方程系数 3.模型的建立 分段建立模型： 记生产批量时，单位成本为，生产批量时，单位成本为。为了大致地分析y与x的关系，首先利用表中表中数据分别作出对和对的散点图。 从图中可以发现对和对成线性关系。所以分别建立线性模型： 模型（1）： 模型（2）： 4.模型的求解 分段模型求解： 将x1和y1的数据分别输入MATLAB： X=[300 340 400 440 480]; Y=[4.65 4.45 4.52 4.20 4.04]; plot(X,Y,.) p=polyfit(X,Y,1) 得到模型（1）的回归系数估计值及其置信水平、检验统计量,F,p的结果见表1. 表一 参数 参数估计值 参数置信区间 5.5863 [4.5743 6.5983] -0.0031 [-0.0056 -0.0006] 然后，对数据进行残差分析： 从结果可以看出，应将第二个点去掉后再进行拟合： 去掉第二个点（340,4.45），再将剩下的点输入MATLAB: X=[300 400 440 480]; Y=[4.65 4.52 4.20 4.04]; plot(X,Y,.) p=polyfit(X,Y,1) 得到模型（1）的回归系数估计值及其置信水平、检验统计量,F, p的结果见表2. 表二 参数 参数估计值 参数置信区间 5.7293 [5.0902 ， 6.0596] -0.0034 [-0.0044 ， -0.0020] F=40.8967 p=0.0238682 可见R的平方非常接近1.说明模型较准确。 于是得到模型（1）： 将x2和y2的数据分别输入MATLAB： z=polyval(p,X); plot(X,Y,k+,X,z,b) hold on X=[540 600 650 720 750 800]; Y=[3.10 2.96 2.48 2.18 1.50 1.38]; plot(X,Y,.) p=polyfit(X,Y,1) 得到模型（1）的回归系数估计值及其置信水平、检验统计量,F,p的结果见表三. 表三 参数 参数估计值 参数置信区间 [5.4316 ， 8.8000] -0.0072 [-0.0096 ， -0.0047] 于是得到模型（2）： 对数据进行残差分析： 由图可知，数据无异常点。 综合模型（1）和（2）可得： 输入程序： z=polyval(p,X); plot(X,Y,k+,X,z,g); hold on X=[300 400 440 480 540 600 650 720 750 800]; Y=[4.65 4.52 4.20 4.04 3.10 2.96 2.48 2.18 1.50 1.38];