数模实验报告-计算机1班-20144138021-张建华-第10次

河北大学《数学模型》实验实验报告班级专业计算机1班姓名张建华学号20144138021实验地点C1-229指导老师司建辉成绩实验项目牙膏的销售量实验13-1牙膏的销售量一、实验目的掌握数据拟合及绘制散点图与拟合曲线组合图形方法；二、实验要求1.建立mydata.m;2.绘制y对x1的散点图。3．确定y对x1的的拟合，绘制散点图与拟合曲线组合图形。4.绘制y对x2的散点图。5．确定y对x2的的拟合，绘制散点图与拟合曲线组合图形。6.y对x1,x2的回归模型及其求解，销售量预测。三、实验内容1.建立mydata.m;2.绘制y对x1的散点图。程序如下：M=dlmread(‘mydata.m’);x1=M(:,5);y=M(:,6);plot(x1,y,’bo’);3．确定y对x1的的拟合，绘制散点图与拟合曲线组合图形。程序如下：clc;formatshortg;M=dlmread(‘mydata.m’);%读取ASCII码文件x1=M(:,5);y=M(:,6);plot(x1,y,’bo’);b=regress(y,[ones(size(x1)),x1]);%b=[β0β1]’，列向量xx=sort(x1);%按升序排序yy=[ones(size(xx)),xx]*b;holdon;plot(xx,yy,’-r’);holdoff;4.绘制y对x2的散点图。程序如下:clc;formatshortg;M=dlmread('mydata.m');%读取ASCII码文件x2=M(:,4);y=M(:,6);plot(x2,y,'bo');5．确定y对x2的的拟合，绘制散点图与拟合曲线组合图形。程序如下:clc;formatshortg;M=dlmread('mydata.m');%读取ASCII码文件x2=M(:,4);y=M(:,6);plot(x2,y,'bo');b=regress(y,[ones(size(x2)),x2,x2.^2]);%b=[β0β1β2]’，列向量xx=sort(x2);yy=[ones(size(xx)),xx,xx.^2]*b;holdon;plot(xx,yy,'-r');holdoff;5.y对x1,x2的回归模型及其求解，销售量预测。程序如下:clc;formatcompact;formatshortg;M=dlmread('mydata.m');%读取ASCII码文件x1=M(:,5);x2=M(:,4);y=M(:,6);[b,bint,r,rint,stats]=regress(y,[ones(size(x1)),x1,x2,x2.^2],0.05);disp('β0,β1,β2,β3估计值置信区间');[b,bint]R2=stats(1)F=stats(2)p=stats(3)disp('销售量预测');x1=0.2x2=6.5y=[1x1,x2,x2^2]*b四、实验结果及其分析1.2.3.4.5.6.实验13-2牙膏的销售量——模型改进（验证）一、实验目的1、掌握多元线性回归函数的使用方法及调用格式；2、了解完全二次多项式模型；二、实验要求1、y对x1,x2的回归模型的改进和求解，销售量预测。（编程）2、完全二次多项式模型（验证），运行结果并得出参数(β0,β1,β2,β3,β4,β5)’的值。三、实验内容1、改进的模型：y=β0+β1x1+β2x2+β3x22+β4x1x2+ε程序如下：formatcompact;formatshortg;M=dlmread('mydata.m');%读取ASCII码文件x1=M(:,5);x2=M(:,4);y=M(:,6);[b,bint,r,rint,stats]=regress(y,[ones(size(x1)),x1,x2,x2.^2，x1.*x2],0.05);disp('β0,β1,β2,β3,β4估计值置信区间');[b,bint]R2=stats(1)F=stats(2)p=stats(3)disp('销售量预测');x1=0.2x2=6.5y=[1x1,x2,x2^2]*b2、完全二次多项式模型（验证）y=β0+β1x1+β2x2++β3x1x2+β4x12+β5x22ε程序如下：clear;clc;formatcompact;formatshortg;M=dlmread('mydata.m');%读取ASCII码文件x1=M(:,5);x2=M(:,4);y=M(:,6);rstool([x1,x2],y,'quadratic')四、实验结果及分析1、：2、: