使用MATLAB进行多元回归分析（自定义函数公式）——nlinfit函数的使用

具体参数设置参考官网：https://ww2.mathworks.cn/help/stats/nlinfit.html?s_tid=srchtitle#btldemp-2如下所示自定义公式为：y=b1*x1/(b2+x1)+b3*x2/(b4+x2)+b5*x3/(b6+x3)+b7*x4/(b8+x4);x为n行4列的矩阵，y为n行1列的矩阵；modelfun = @(beta,x)(beta(

潘聪明

6155人浏览 · 2020-11-06 12:24:08

潘聪明 · 2020-11-06 12:24:08 发布

具体参数设置参考官网：https://ww2.mathworks.cn/help/stats/nlinfit.html?s_tid=srchtitle#btldemp-2

如下所示**自定义公式为：y=b1x1/(b2+x1)+b3x2/(b4+x2)+b5x3/(b6+x3)+b7x4/(b8+x4);**

x为n行4列的矩阵，y为n行1列的矩阵；

modelfun = @(beta,x)(beta(1).*x(:,1)./(beta(2)+x(:,1))+beta(3).*x(:,2)./(beta(4)+x(:,2))+beta(5).*x(:,3)./(beta(6)+x(:,3))+beta(7).*x(:,4)./(beta(8)+x(:,4)));
rng('default') % for reproducibility
opts = statset('nlinfit');
opts.RobustWgtFun = 'bisquare';
beta0 = [2;2;2;2;2;2;2;2];初始beta值的设定
[beta,R,J,CovB,MSE,ErrorModelInfo]= nlinfit(x,y,modelfun,beta0,opts)
yuce_y=modelfun(beta,x)%依据自定义公式预测的y值