数学建模——非线性回归（b站 数学建模加油站 学习笔记）

一、多项式回归求解

（一）一元多项式求解

直接上例题：

%导入数据
t = 1/30: 1/30: 14/30;
s =[11.86 15.67 20.6 26.69 33.71 41.93 51.13 61.49 72.9 85.44 99.08 113.77 129.54 146.48]
%使用求解多项式系数的命令
[p,S]= polyfit(t,s,2)
%做曲线观察拟合情况
Y = polyconf(p,t,S)
plot(t,s,'k+',t,Y,'r')

结果为：

s=489.29*t^2+65.89*t+9.13

（二）多元二项式求解

直接上例题：

%导入数据
x1 = [1000 600 1200 500 300 400 1300 1100 1300 300] ;
x2=[5 7 6 6 8 7 5 4 3 9];
y= [100 75 80 70 50 65 90 100 110 60]'; 
x= [x1' x2'];
%利用rstool命令求解
rstool(x ,y, 'purequadratic')

运行后会出现模型预测图，分别输入1000、6，点击“导出”即可获得拟合系数等数据：

---

二、非线性回归一般求解

 脚本代码为：

% 1)输入数据
x=2:16;
y=[6.42 8.20 9.58 9.5 9.7 10 9.93 9.99 10.49 10.59 10.60 10.80 10.60 10.90 10.76]
beta0=[8 2]'; %设定初始值，根据经验选择
% 2)求回归系数
[beta, r,J]=nlinfit(x' ,y','volum', beta0);
%使用非线性最小二乘法进行回归拟合，'volum' 为定义曲线函数的函数名
beta

%3)预测及作图
[YY, delta]=nlpredci( 'volum' ,x', beta,r,J); %预测因变量的值及置信区间
plot(x,y,'k+' ,x,YY,'r') % 绘制原始数据点和拟合曲线

volum自定义函数代码为：

function yhat=volum(beta,x)
yhat=beta(1)*exp(beta(2)./x);
end

运行结果为：

beta=11.6037

          -1.0641

所以回归方程为 

y=11.6037*