matlab多元线性回归模型

1、云南大学数学与统计学实验教学中心实验报告课程名称：数学实验学期：20102011学年下学期成绩：指导教师：学生姓名：学生学号：实验名称：回归分析实验编号：六实验日期：6月6日实验学时：2学院：专业：信息与计算科学年级：一、实验目的1,熟悉MATLAB勺运行环境.2.学会初步建立数学模型的方法3,运用回归分析方法来解决问题二、实验内容实验一：某公司出口换回成本分析对经营同一类产品出口业务的公司进行抽样调查，被调查的13家公司，其出口换汇成本与商品流转费用率资料如下表。试分析两个变量之间的关系，并估计某家公司商品流转费用率是6.5%的出口换汇成本.公司出口换汇成本（人民币元/美元）商品流转费用率（

2、%）公司出口换汇成本（人民币元/美元）商品流转费用率（%）11.404.2081.605.5021.205.3092.004.1031.007.10101.005.0041.903.70111.604.0051.306.20121.803.4062.403.50131.406.9071.404.80实验二：某建筑材料公司的销售量因素分析下表数据是某建筑材料公司去年20个地区的销售量（Y,千方），推销开支、实际帐目数、同类商品竞争数和地区销售潜力分别是影响建筑材料销售量的因素。1)试建立回归模型，且分析哪些是主要的影响因素。2)建立最优回归模型。地区i推销开支(x1)实际帐目数(x2)同类商品竞

3、争数(x3)地区销售潜力(x4)销售量Y15.53110879.322.55586200.138.067129163.243.050716200.153.038815146.062.9711217177.778.03012830.989.056510291.994.04284160.0106.573516339.4115.560117159.6125.044121286.3136.05066237.5145.039104107.2153.555104155.0168.070614201.4176.040116100.2184.050118135.8197.562913223.3207.0599

4、11195.0提示：建立一个多元线性回归模型。三、实验环境Windows操作系统；MATLAB7.0.四、实验过程实验一：运用回归分析在MATLAEB!实现输入：x=4.205.307.103.706.203.504.805.504.105.004.003.406.90'X=ones(13,1)x;Y=1.401.201.001.901.302.401.401.602.001.001.601.801.40'plot(x,Y,'*');b,bint,r,rint,stats=regress(Y,X,0.05);输出：b=2.6597-0.2288bint=1.88

5、733.4322-0.3820-0.0757stats=0.495810.81680.00720.0903即为=2.6597,?,1=-0.2288,用的置信区间为1.88733.4324区1的置信区间为-0.38202一_一，八r=0.4958,F=10.8168,p=0.0072因P<0.05,可知回归模型y=2.6597-0.2288x成立.-0.0757;2.5散点图21.533.544.555.566.577.5估计某家公司商品流转费用率是6.5%的出口换汇成本。将x=6.5代入回归模型中，得到>>x=6.5;>>y=2.6597-0.2288*xy=1

6、.1725实验二：在MATLAB!实现，首先建立回归模型输出：x1=5.52.58.03.03.02.98.09.04.06.55.55.06.05.03.58.06.04.07.57.0'x2=3155675038713056427360445039557040506259'x3=108127812125851112610106111199'x4=86916151781041671264414681311'Y=79.3200.1163.2200.1146.0177.730.9291.9160.0339.4159.686.3237.5107.2155.0201.

7、4100.2135.8223.3195.0'X=ones(20,1)x1x2x3x4;b,bint,r,rint,stats=regress(Y,X,0.05);b,bint,stats输出：b=191.9158-0.77193.1725-19.6811-0.4501bint=103.1071280.7245-7.14455.60072.06404.2809-25.1651-14.1972-3.72842.8283stats=0.903435.05090.0000644.6510即%=191.9158/产-0.7719居=3.172533=-19.6811%=-0.4501;陷的置信区

8、间为103.1071280.7245;肉的置信区间为-7.14455.6007;仅2的置信区间为2.06404.2809;眩的置信区间为-25.1651-14.1972;仅4的置信区间为-3.72842.8283;2r=0.9034,F=35.0509,p=0.0000因P<0.05,可知回归模型y=191.9158-0.7719x1+3.1725*x2-19.6811*x3-0.4501*x4成立.分析哪些是主要的影响因素输入：x1=5.52.58.03.03.02.98.09.04.06.55.55.06.05.03.58.06.04.07.57.0'x2=315567503

9、8713056427360445039557040506259'x3=108127812125851112610106111199'x4=86916151781041671264414681311'Y=79.3200.1163.2200.1146.0177.730.9291.9160.0339.4159.686.3237.5107.2155.0201.4100.2135.8223.3195.0'X=x1x2x3x4;stepwise(X,Y);X1esmr1234XXXXesmr6.534440.77680.4473-4.028714.41920.00039一-

10、23.935-5.38770.00006395711740100989rrfrrCoefficientswithErrorBarsCoeff.t-statp-val-30-20-1001020x2x4Nextstep:MoveX3inAHStepsbntercefrt=169.49SR-square-0F-MaMRMSE=72,5666Ac(jR-sq=-0,(?52631i5p=NaNModelHistory从表StepwiseTable中分析得出变量x2和x3为主要的影响因素。,M-0.701486-0.24240.8115-3.090677.04970.0000L-19514_81596

11、00000_n41ftd97_n9R11n7Q9QrirrrrCoefficientswithErrorBarsCoeff.t-statp-val-25-20-15-10-505Nextstep:I.MovenotermsExport.Utencept=186,046R-square=C.9O2444F=76,6295RMSE=23SEI6Ac|jR-sq0,305223p=256246e-W9ModelHistory806040移去非关键变量x1和x4后模型具有显著性.虽然剩余标准差(RMSE都有了变化，统计量F的值明显增大，因此新的回归模型更好.就得到最优模型。输入：X1=ones(20,1)x2x3;b,bint,r,rint,stats=regress(Y,X1);b,bint,stats输出：b=186.04843.0907-19.5140bint=110.4254261.67152.16574.0156-24.5597-14.4683stats=0.902478.62