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）实际帐目数

3、（x2）同类商品竞争数（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.

4、3207.059911195.0提示：建立一个多元线性回归模型。三、实验环境Windows 操作系统；MATLAB 7.0.四、实验过程实验一：运用回归分析在MATLAB里实现输入：x=4.20 5.30 7.10 3.70 6.20 3.50 4.80 5.50 4.10 5.00 4.00 3.40 6.90'X=o nes(13,1) x;Y=1.40 1.20 1.00 1.90 1.30 2.40 1.40 1.60 2.00 1.00 1.60 1.80 1.40'plot(x,Y,'*');b,bi nt,r,ri nt,stats=regress

5、(Y,X,0.05);输出：b =2.6597-0.2288bint =1.88733.4322-0.3820-0.0757stats =0.495810.81680.00720.0903即 彳2.6597,勺 -0.2288,?）的置信区间为1.88733.4322,冷的置信区间为-0.3820 0.0757；r 2 =0.4958, F=10.8168, p=0.0072散点图33.544.55.566.577.5因P<0.05, 可知回归模型 y=2.6597-0.2288x 成立.估计某家公司商品流转费用率是6.5%的出口换汇成本。将x=6.5代入回归模型中，得到>>

6、 x=6.5;>> y=2.6597-0.2288*xy =1.1725实验二：在 MATLAB 里实现， 首先建立回归模型输出：x1=5.5 2.5 8.0 3.0 3.0 2.9 8.0 9.0 4.0 6.5 5.5 5.0 6.0 5.0 3.5 8.0 6.0 4.0 7.5 7.0'x2=31 55 67 50 38 71 30 56 42 73 60 44 50 39 55 70 40 50 62 59'x3=10 8 12 7 8 12 12 5 8 5 11 12 6 10 10 6 11 11 9 9'x4=8 6 9 16 15 17

7、8 10 4 16 7 12 6 4 4 14 6 8 13 11'Y=79.3 200.1 163.2 200.1 146.0 177.7 30.9 291.9 160.0 339.4 159.6 86.3 237.5 107.2 155.0 201.4100.2 135.8 223.3 195.0'X=o nes(20,1) x1 x2 x3 x4;b,bi nt,r,ri nt,stats=regress(Y,X,0.05);b,bi nt,stats输出：b =191.9158-0.77193.1725-19.6811-0.4501bint =103.1071280.7

8、245-7.14455.60072.06404.2809-25.1651-14.1972-3.72842.8283stats =0.903435.05090.0000644.6510即 ?)= 191.9158?严0.7719?2= 3.1725?3=-19.6811=-0.4501 ;?的置信区间为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<

9、;0.05, 可知回归模型 y=191.9158 -0.7719x1+3.1725*x2-19.6811*x3 -0.4501*x4成立. 分析哪些是主要的影响因素输入：x1=5.5 2.5 8.0 3.0 3.0 2.9 8.0 9.0 4.0 6.5 5.5 5.0 6.0 5.0 3.5 8.0 6.0 4.0 7.5 7.0'x2=31 55 67 50 38 71 30 56 42 73 60 44 50 39 55 70 40 50 62 59'x3=10 8 12 7 8 12 12 5 8 5 11 12 6 10 10 6 11 11 9 9'x4=8

10、 6 9 16 15 17 8 10 4 16 7 12 6 4 4 14 6 8 13 11'Y=79.3 200.1 163.2 200.1 146.0 177.7 30.9 291.9 160.0 339.4 159.6 86.3 237.5 107.2 155.0 201.4100.2 135.8 223.3 195.0'X=x1 x2 x3 x4;stepwise(X,Y);Coefficients with Error BarsCoeff. t-stat p-valX16.53444 0.7768 0.4473X24.02871 4.4192 0.0003X3-23

11、.935 -5.3877 0.0000X46.39571 1.7401 0.0989-30-20-101020Next step:Move X3 inNext StepesmrAl StepsMtercept = 169.4SSR-square - 0F - MaMRMSE = 72.5666Acj R-sq = -0.(J526315p= MaNModel HistoryExport .从表Stepwise Table中分析得出变量 x2和x3为主要的影响因素。-0.701486 -0.2424 0.81153.09067 7.0497 0.0000-19.514 -8.1596 0.000

12、0-0.418127 -0.2811 0.7823-rFirrr4rCoefficients with Error BarsCoeff. t-stat p-val-25-20-15-10-505X1x2x4Next StepI；Next step:Move no termsExporl .Mtercept = 186.046R-Square =0.902444F=7fi.6295RMSE = 23.9616Acj R-sq = O.S0522Sp=2.56246e-009806040esmr20Model History 移去非关键变量x1和x4后模型具有显著性.虽然剩余标准差(RMSE )都有了变化，统计量 F 的值明显增大，因此新的回归模型更好就得到最优模型。输入：X1=o nes(20,1) x2 x3;b,bi nt,r,ri nt,stats=regress(Y,X1);b,bi nt,stats输出：b =186.04843.0907-19.5140bint =110.4254261.67152.16574.0156-24.5597-14.4683