实用 MATLAB教学资料-21-22最小二乘续.pptVIP

例 某化学反应实验所得生成物的浓度随时间的变化数据，（1）拟合模型 参数 （2）绘出拟合曲线和数据点 t 1 2 3 4 5 6 7 8 y×103 4.00 6.40 8.00 8.80 9.22 9.50 9.70 9.86 t 9 10 11 12 13 14 15 16 y×103 10.00 10.20 10.32 10.42 10.50 10.55 10.58 10.60 [x,resnorm,residual,exitflag,output] = lsqnonlin(@fun,x0) fun:函数式m文件 原理：寻找满足如下条件的模型参数B 2、lsqcurvefit B-待定模型参数（向量） xdata-已知数据自变量（向量） ydata-已知数据因变量（向量） 格式 [B,resnorm,residual,exitflag,output]= lsqcurvefit(@fun,b0,xdata,ydata,lb,ub,options) 输入： fun-函数（模型）表达式 ! 不是残差平方和，也不是残差的表达式。 MATLAB自行计算残差平方和?。 匿名函数 函数式m文件 fun=@(B,xdata)B(1)*xdata.^2+B(2)*sin(xdata)...; [B,resnorm,residual,exitflag,output]= lsqcurvefit(fun,b0,xdata,ydata,lb,ub,options) 定义模型表达式——匿名函数 第一个，待定参数 第二个，自变量 function f = fun(B,xdata) %MATLAB根据位置读取 数据，第一个是待定参数，第二个是自变量 f =B(1)*xdata.^2+B(2)*sin(xdata)...; [B,resnorm,residual,exitflag,output]= lsqcurvefit(@fun,b0,xdata,ydata,lb,ub,options) 定义模型表达式——函数式m文件 格式 [B,resnorm,residual,exitflag,output]= lsqcurvefit(@fun,b0,xdata,ydata,lb,ub,options) 输入： xdata-实测数据自变量 ydata-实测数据因变量 例 测定活性为100%某催化剂的活性 y 随时间 t 变化情况，数据如下： ti 0 27 40 52 70 89 100 yi 100.0 82.2 76.3 71.8 66.4 63.3 61.3 用非线性最小二乘确定a,b,c,并绘制yi-ycal图 几种解法？ 例 测定活性为100%某催化剂的活性 y 随时间 t 变化情况，数据如下： ti 0 27 40 52 70 89 100 yi 100.0 82.2 76.3 71.8 66.4 63.3 61.3 用非线性最小二乘确定a,b,c,并绘制yi-ycal图 [x,resnorm,residual,exitflag,output] = lsqcurvefit(fun,x0,xdata,ydata) 例 测定活性为100%某催化剂的活性 y 随时间 t 变化情况，数据如下： ti 0 27 40 52 70 89 100 yi 100.0 82.2 76.3 71.8 66.4 63.3 61.3 用非线性最小二乘确定a,b,c,并绘制yi-ycal图 [B,resnorm,residual,exitflag,output] = lsqnonlin(fun,b0) nested function嵌套函数 实现参数传递 第9章 最优化技术 9.2 曲线拟合工具箱Curve Fitting Tool 启动 Curve Fitting Tool APPS点击Curve Fitting按钮 或 Command窗口输入cftool 数据区 从workspace中选择 给拟合任务命名 模型区 选择或自定义模型 拟合方法参数设置 操作区 自动或手动拟合 强制终止Stop 结果区 模型参数 拟合效果评价指标 残差平方和SSE 相关指数R2 … 绘图区 清单 拟合任务名称 拟合类型 拟合效果评价指标 残差平方和SSE 相关指数R2 …… 进入figure窗口编辑图形 创建新的拟合任务 显示残差图 Tools Data cursor-读取数据点坐标 Exclude out