Matlab多项式回归实现

时间 2019-11-19

原文原文链接

　　多项式回归也称多元非线性回归，是指包含两个以上变量的非线性回归模型。对于多元非线性回归模型求解的传统解决方案，仍然是想办法把它转化成标准的线性形式的多元回归模型来处理。html

多元非线性回归分析方程

　　若是自变数 $X_1,X_2,\cdots,X_m$ 与依变数Y皆具非线性关系，或者有的为非线性有的为线性，则选用多元非线性回归方程是恰当的。例如，二元二次多项式回归方程为：ide

　　 $\widehat{y}=a+b_{11}x_1+b_{21}x_2+b_{12}x_1^2+b_{22}x_2^2+b_{11\times22}x_1x_2$

　　令 $b_1=b_{11},b_2=b_{21},b_3=b_{12},b_4=b_{22},b_5=b_{11\times22}$ ,及 $x_3=x_1^2,x_4=x_2^2,x_5=x_1\cdot x_2,$ 因而上式化为五元一次线性回归方程：函数

　　 $\widehat{y}=a+b_1x_1+b_2x_2+b_3x_3+b_4x_4+b_5x_5$

　　这样以来，即可按多元线性回归分析的方法，计算各偏回归系数，创建二元二次多项式回归方程。spa

-参考文献：智库百科，点击打开htm

多元二项式回归Matlab实现方法

一、多元二项式回归Matlab命令

rstool(x,y,'model',alpha)
输入参数说明：
x：n*m矩阵；
Y：n维列向量；
alpha：显著性水平(缺省时为0.05)；
mode：由下列4个模型中选择1个(用字符串输入,缺省时为线性模型)
　　'model'对应的字符串有：'linear'%线性常数和线性项（默认值）
　　'interaction'%常数、线性和相互做用项
　　'quadratic'%纯二次拟合函数
　　'purequadratic'%blog

二、实例演示说明

设某商品的需求量与消费者的平均收入、商品价格的统计数据以下,创建回归模型,预测平均收入为1000、价格为6时的商品需求量字符串

需求量(y) 100 75 80 70 50 65 90 100 110 60
收入(x1) 1000 600 1200 500 300 400 1300 1100 1300 300
价格(x2) 5 7 6 6 8 7 5 4 3 9get

解法一：选择纯二次模型'model'='quadratic'

y=β0+β1x1+β2x2+β11x1^2+β22x2^2
%直接用多元二项式回归以下
x1=[1000 600 1200 500 300 400 1300 1100 1300 300];%1*10矩阵
x2=[5 7 6 6 8 7 5 4 3 9];%1*10矩阵
y=[100 75 80 70 50 65 90 100 110 60]';%1*10矩阵
x=[x1' x2'];%10*2矩阵
rstool(x,y,'purequadratic')it

在x1对应的文本框中输入1000，X2中输入6，敲回车键，此时图形和相关数据会自动更新
此时在GUI左边的“Predicted Y1”下方的数据变为88.47981，表示平均收入为1000、价格为6时商品需求量为88.4791
点击左下角的Export按钮，将会导出回归的相关参数beta、rmse和residuals到工做空间(workspace)
在Export按钮下面能够选择回归类型
在Matlab命令窗口中输入
>>beta, rmse
beta =
110.5313%常数项
0.1464
-26.5709
-0.0001
1.8475
rmse =
4.5362%均方根偏差
由此得回归模型为：y=110.5351+0.1464x1-26.5709x2-0.0001x1^2+1.8475x2^2io

解法二：将上面的模型转换为多元线性回归

y=β0+β1x1+β2x2+β11x1^2+β22x2^2　　疑问：为何没有变量x1*x2呢？官方例子有：地址>>X=[ones(10,1) x1' x2' (x1.^2)' (x2.^2)'];>>[b,bint,r,rint,stats]=regress(y,X);>>b,statsb = 110.53130.1464 -26.5709-0.00011.8475stats =0.9702 40.6656 0.0005 20.5771%相关系数R^二、F值、与F值对应的几率（P值）、偏差方差估计值