非线性最小二乘数据拟合高斯牛顿法.pptVIP

1 可选性曲线绘制 首先解决计算问题: 煤炭浮沉试验综合综合表的计算是选煤工作者最基本最常用的一种计算，表中除浮沉试验基础数据外均为空白栏，需要填入各计算结果。 试按表中的原煤基元灰分与产率计算累计浮物产率、沉物产率、灰分及邻近密度物产率。 50~0.5mm级原煤浮沉试验综合表 密度级 /kg·L-1 产率 /% 灰分 1% 累 计 分选密度±0.1 /kg·L-1 浮 物 沉 物 产率 /% 灰分 /% 产率 /% 灰分 /% 密度 /kg·L-1 产率 /% 1 2 3 4 5 6 7 8 9 -1.3 48.00 5.30 48.00 5.30 100.00 28.41 1.30 58.30 1.30~1.40 10.30 9.26 58.30 6.00 52.00 49.75 1.40 18.00 1.40~1.50 7.70 23.19 66.00 7.99 41.70 59.75 1.50 11.70 1.50~1.60 4.00 40.00 70.00 9.82 34.00 68.05 1.60 5.75 1.60~1.80 3.50 47.50 73.50 11.62 30.00 71.79 1.70 3.50 +1.80 26.50 75.00 100.00 28.41 26.50 75.00 合计 100.00 28.41 当变量之间为非线性相关时，可用非线性最小二乘数据拟合（高斯—牛顿法）。 [beta,r,j]=nlinfit(x,y, ′fun′,beta0) [ypred,delta]=nlpredci(FUN,inputs,beta,r,j) ci=nlparci(beta,r,j) nlintool(x,y, ′fun′,beta0) nlinfit 非线性拟合函数。beta是以x,y为数据返回的系数值。fun是系数向量和数组x的函数，返回拟合y值的向量。beta0为选取的初始值向量，r为拟合残差，j为Jacobian矩阵值。 nlpredci 非线性最小二乘预测置信区间。nlparci 非线性模型参数置信区间。ypred为预测值，delta为置信区间的半长值，inputs为矩阵。ci为b的误差估计。 nlintool 非线性拟合交互式图形工具。显示95%置信区间上下的两条红线和其间的拟合曲线。移动纵向虚线可显示不同的自变量及其对应的预测值。还可有其他参数。 非线性拟合相关命令 例题1: 某筛分试验测得筛孔尺寸与累计产率间的关系如表所示，试建立模型。 产率与筛孔关系 筛孔/mm 50 25 13 6 3 0.5 产率/% 1.5 11 30 55 70 88 此为一近似指数关系，即 型。 用MATLAB编程 首先建立一个′ep′函数文件。 function f=ep(b,x)； f=(b(1) * exp(－b(2) * x) 输入 x=[50 25 13 6 3 0.5]’; y=[1.5 11 30 55 70 88]; b0=[1 1]; %初始值 [b,r,j]=nlinfit(x,y, ′ep′,b0);b [yp,d]=nlperdci(′ep′,x,b,r,j); ci=nlparci(b,r,j) %输出预测值、残差及置信区间 nlintool(x,y, ′ep′,b0); %绘制非线性拟合曲线图3.14 plot(x,r, ′――′) %绘残差曲 结果 得 b= 91.4158 0.0855 ci= 90.2085 92.6231 0.0828 0.0883 移动垂直于x轴的虚线可以显示y随x的变化。数学模型为 y=91.4158e-0.0855x 例2 已知浮煤累计灰分及产率为x和y，试进行非线性拟合。 x=[3.46 7.33 9.47 10.48 11.19 11.79 12.78 20.50]； y=[10.69 56.84 76.89 82.15 84.70 86.32 88.45 100]； 现采用含有双曲正切的函数拟合。 f=100 * (b(1)+b(2) * x +b(2) * b(3)*tanh(b(4)*(x-b(5))))； 输入 x=x’； b0=[0.2 0 -0.4 -0.1 5]；%初始值 [b,r,j]=nlinfit(x,y, ′mf′,b0)； [yp,d]=nlpredci (′mf′,x,b,r,j); nlintool(x,y, ′mf′,b0) 结果 b= 0.3414 0.0136 -0.3797 -0.3365 6.2948 即 y=100*(0.3414+0.0136