MATLAB在化工中的应用_插值拟合与数值微分和积分重点.pptVIP

第四部分 插值、拟合与数值微分和积分;本章知识要点;插值、拟合、数值微分、数值积分在化工计算中的作用;插值简介;插值方法的选择;Matlab的插值（Interpolation）函数;一维插值interp1;已知x=[0:10]，y=[0 0.8415 0.9093 0.1411 -0.7568 -0.9589 -0.2794 0.6570 0.9894 0.4121 -0.5440];(y＝sinx)，比较一维线性、线性最近、立方和三次样条插值所得xi＝0,0.15,0.30,0.45,…,10处的值yi。如果初始数据点为x＝0,2,4,…,10，y＝sinx，以上方法插值效果。 ;初始数据对于插值的影响;spline与pchip;pchip与spline的区别;其它样条插值函数－csape;csape的使用;几个有用的函数;二维插值：interp2;二维插值函数的使用;拟合简介;最小二乘多项式拟合：polyfit;多项式次数对拟合效果的影响;最小二乘法拟合生成样条曲线;功能：平滑生成三次样条函数，即对于数据(xi,yi)，所求的三次样条函数y=f(x)满足 ：;功能：用最小二乘法拟合生成B样条曲线，即对于数据(xi,yi)，所求k次样条函数y=f(x)满足 ：;功能：平滑生成B样条函数 ;采用样条拟合函数csaps进行拟合，比较p的取值对于拟合效果的影响 ;数值微分;Matlab数值微分实现方法;函数diff;函数gradient;例题;数值积分;数值积分;梯形法数值积分：trapz();调用格式： q=quad(@fun,a,b) q=quad(@fun,a,b,tol) q=quad(@fun,a,b,tol,trace,p1,p2,……) 输入参数： fun 被积函数。在定义fun时，被积函数表达 式必须是向量形式，即表达式必须使用点 运算符(.*、./和.^)以支持向量 a,b 即积分限[a,b] tol 绝对误差限，默认值为1.e-6 p1,p2…… 直接传递给函数fun的已知参数 输出参数： q 积分结果 自适应Lobatto法数值积分：quadl() 调用格式同quad;不同积分函数的比较;广义积分;多重数值积分;样条函数在数值积分与微分中的应用;反应器停留时间分布的混合特性;;微分法进行动力学数据分析;function Cha4demo8 % 动力学数据 t = [0 20 40 60 120 180 300];CA = [10 8 6 5 3 2 1]; % 用最小二乘样条拟合法计算微分dCA/dt knots = 3;K = 3; % 三次B样条 sp = spap2(knots,K,t,CA); sp = spap2(newknt(sp),K,t,CA); pp = fnder(sp); % 计算B样条函数的导函数 dCAdt = fnval(pp,t); % 计算t处的导函数值 % ???制图形 ti = linspace(t(1),t(end),200); CAi = fnval(sp,ti); plot(t,CA,‘ro’,ti,CAi,‘b-’)，xlabel(‘t’)，ylabel(C_A) figure fnplt(pp); % dCAdti = fnval(pp,ti) % plot(ti,dCAdti,-) xlabel(t) ylabel(dC/dt) % 线性拟合 rA = dCAdt; y = log(-rA); x = log(CA); p = polyfit(x,y,1); k = exp(p(2)) n = p(1)