数学建模-数据处理分解.ppt

* * * * * * * * * * * 例4 财政收入预测问题：财政收入与国民收入、工业总产值、农业总产值、总人口、就业人口、固定资产投资等因素有关。下表列出了1952-1981年的原始数据，试构造预测模型。 解 设国民收入、工业总产值、农业总产值、总人口、就业人口、固定资产投资分别为x1、x2、x3、x4、x5、x6，财政收入为y，设变量之间的关系为： y= ax1+bx2+cx3+dx4+ex5+fx6 使用非线性回归方法求解。 技能节竞赛第二题 * * 1．? 对回归模型建立M文件model.m如下: function yy=model(beta0,X) a=beta0(1); b=beta0(2); c=beta0(3); d=beta0(4); e=beta0(5); f=beta0(6); x1=X(:,1); x2=X(:,2); x3=X(:,3); x4=X(:,4); x5=X(:,5); x6=X(:,6); yy=a*x1+b*x2+c*x3+d*x4+e*x5+f*x6; * * 2. 主程序liti6.m如下: X=[598.00 349.00 461.00 57482.00 20729.00 44.00 ………………………………………………………….. 2927.00 6862.00 1273.00 100072.0 43280.00 496.00]; y=[184.00 216.00 248.00 254.00 268.00 286.00 357.00 444.00 506.00 ... 271.00 230.00 266.00 323.00 393.00 466.00 352.00 303.00 447.00 ... 564.00 638.00 658.00 691.00 655.00 692.00 657.00 723.00 922.00 ... 890.00 826.00 810.0]; beta0=[0.50 -0.03 -0.60 0.01 -0.02 0.35]; betafit = nlinfit(X,y,model,beta0) To MATLAB(liti6） * * betafit = 0.5243 -0.0294 -0.6304 0.0112 -0.0230 0.3658 即y= 0.5243x1-0.0294x2-0.6304x3+0.0112x4-0.0230x5+0.3658x6 结果为: 返 回 * * 1、考察温度x对产量y的影响，测得下列10组数据： 求y关于x的线性回归方程，检验回归效果是否显著，并预测x=42℃时产量的估值及预测区间（置信度95%）. 2、某零件上有一段曲线，为了在程序控制机床上加工这一零件，需要求这段曲线的解析表达式，在曲线横坐标xi处测得纵坐标yi共11对数据如下： 求这段曲线的纵坐标y关于横坐标x的二次多项式回归方程. * * * * 4、混凝土的抗压强度随养护时间的延长而增加，现将一批混凝土作成12个试块，记录了养护日期x（日）及抗压强度y（kg/cm2）的数据： * * * * * * * * * * * * * * * * * * * * * * * * 教学目的 教学内容 2、掌握用数学软件求解回归分析问题。 1、直观了解回归分析基本内容。 1、回归分析的基本理论。 2、用数学软件求解回归分析问题。 * * 一元线性回归 多元线性回归 回归分析 数学模型及定义 *模型参数估计 *检验、预测与控制 可线性化的一元非线 性回归（曲线回归） 数学模型及定义 *模型参数估计 *多元线性回归中的 检验与预测 逐步回归分析（了解） * * 一、一元线性回归数学模型 例1 测16名成年女子的身高与腿长所得数据如下： 以身高x为横坐标，以腿长y为纵坐标将这些数据点（xI，yi）在平面直角坐标系上标出. 散点图 解答 * * 一元线性回归分析的主要任务是： 返回 * * 模型参数估计 1、回归系数的最小二乘估计 * * * * 返回 * * 检验、预测与控制 1、回归方程的显著性检验 * * （Ⅰ）F检验法 （Ⅱ）t检验法 * * （Ⅲ）r检验法 * * 2、回归系数的置信区间 * * 3、预测与控制 （1）预测 * * （2）控制 返回 * * 一、多元线性回归