源代码+代码解释+代码应用案例+参考
本篇用来介绍非线性回归方程一般代码的样式,里面包含了简单非线性模型代码、简单预测的代码两个部分,欢迎大家补充、纠错、改错~
你敢改,俺的万万个论坛币就敢送~
1、非线性回归源代码
- [beta,r,J]=nlinfit(x,y,'modelfun',beta0) %非线性回归系数的命令
- nlintool(x,y,'modelfun',beta0,alpha) %非线性回归GUI界面
- [Y,DELTA]=nlpredci('modelfun', x,beta,r,J) %获取x处的预测值Y及预测值的显著性为1-alpha的置信区间Y±DELTA
2、非线性回归代码解释
beta:估计出的回归系数;
r:残差;
J:Jacobian矩阵;
x,y:输入数据x、y分别为矩阵和n维列向量,对一元非线性回归,x为n维列向量;
modelfun:M函数、匿名函数或inline函数,定义的非线性回归函数;
beta0:回归系数的初值;
3、应用案例
解:
(1)对将要拟合的非线性模型,建立M函数如下
- function yhat=modelfun(beta,x)
- %beta是需要回归的参数
- %x是提供的数据
- yhat=beta(1)*exp(beta(2)./x);
(2)输入数据
- 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]';
(3)求回归系数
- [beta,r ,J]=nlinfit(x',y',@modelfun,beta0);
- beta
- beta =
- 11.6036
- -1.0641
即得回归模型为
(4)预测及作图
- [YY,delta]=nlpredci('modelfun',x',beta,r ,J);
- plot(x,y,'k+',x,YY,'r')
4、参考
来源于MATLAB技术论坛。





雷达卡




的时候,雨也是晴;心雨
的时候,晴也是雨!
京公网安备 11010802022788号







