非线性回归分析.docxVIP

非线性回归分析 回归分析中，当研究的因果关系只涉及因变量和一个自变量时，叫做一元回归分析；当研究的因果关系涉及因变量和两个或两个以上自变量时，叫做多元回归分析。此外，回归分析中，又依据描述自变量与因变量之间因果关系的函数表达式是线性的还是非线性的，分为线性回归分析和非线性回归分析。通常线性回归分析法是最基本的分析方法，遇到非线性回归问题可以借助数学手段化为线性回归问题处理 两个现象变量之间的相关关系并非线性关系，而呈现某种非线性的曲线关系，如：双曲线、二次曲线、三次曲线、幂函数曲线、指数函数曲线(Gompertz)、S 型曲线(Logistic) 对数曲线、指数曲线等，以这些变量之间的曲线相关关系，拟合相应的 回归曲线，建立非线性回归方程，进行回归分析称为非线性回归分析 常见非线性规划曲线 1 双曲线 ? a ? b y x 二次曲线 三次曲线 幂函数曲线 指数函数曲线(Gompertz) 倒指数曲线 y=a eb / x 其中 a0， S 型曲线(Logistic) y ? 1 a ? be? x 对数曲线 y=a+blog x,x0 指数曲线 y=a ebx 其中参数 a0 1．回归： 确定回归系数的命令[beta，r，J]=nlinfit（x,y,’model’,beta0） 非线性回归命令：nlintool（x，y，’model’, beta0，alpha） 2．预测和预测误差估计： [Y，DELTA]=nlpredci（’model’, x，beta，r，J） 求nlinfit 或lintool 所得的回归函数在x 处的预测值Y 及预测值的显著性水平为 1-alpha 的置信区间Y， DELTA. t(s)1/302/303/304/305/306/307/30s(cm)11.8615.6720.6026.6933.7141.9351.13t(s)8/309/3010/3011/3012/30 t (s) 1/30 2/30 3/30 4/30 5/30 6/30 7/30 s (cm) 11.86 15.67 20.60 26.69 33.71 41.93 51.13 t (s) 8/30 9/30 10/30 11/30 12/30 13/30 14/30 s (cm) 61.49 72.90 85.44 99.08 113.77 129.54 146.48 解： 对将要拟合的非线性模型 y=a eb / x ，建立 M 文件 volum.m 如下： function yhat=volum(beta,x) yhat=beta(1)*exp(beta(2)./x); 输入数据： 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．求回归系数： x[beta,r ,J]=nlinfit(x,y,volum,beta0)； beta x 即得回归模型为： 4．预测及作图： y ? 11.6036e?1.0641 [YY,delta]=nlpredci(volum,x,beta,r ,J)； plot(x,y,k+,x,YY,r) 曲线方程曲线图形 曲线方程 曲线图形 变换公式 变换后的线性函数 y＝axb c＝ln a v＝ln x u＝ln y u＝c ? bv y＝aebx c＝ln a u＝ln y u＝c ? bv y＝a b ex c＝ln a v ? 1 u＝lxn y u＝c ? bv y＝a ? b ln x v＝ln x u＝y u＝a ? bv