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

第四部分 插值、拟合与数值微分和积分 本章知识要点 插值简介 Matlab的插值（Interpolation）函数 调用格式： yi＝interp1(x,y,xi) 已知数据向量（x,y），计算并返回在插值向量xi处的函数值 yi=interp1(x,y,xi, ‘method’) yi=interp1(x,y,xi, ‘method’, ‘extrap’) ‘method’用于指定插值算法，其值可以是： ‘nearest’——最近插值 ‘linear’——线性插值（默认值） ‘spline’——分段三次样条插值 ‘pchip’——分段三次Hermite插值 ‘cubic’——与‘pchip’相同 初始数据对于插值的影响 二维插值：interp2 调用格式： zi=interp2(x,y,z,xi,yi,’method’) ‘method’算法属性值可以是； ‘nearest’——最近插值 ‘linear’——线性插值（默认） ‘spline’——三次样条插值(spline) ‘cubic’——立方插值 拟合简介 最小二乘多项式拟合：polyfit 多项式次数对拟合效果的影响 最小二乘法拟合生成样条曲线 拟合得到曲线函数sp以后，可利用fnval()计算任意自变量下的函数值。 数值微分 函数diff 函数gradient 例题 样条函数在数值积分与微分中的应用 反应器停留时间分布的混合特性 微分法进行动力学数据分析 样条函数：照一定的光滑性要求,连接起来的分段多项式。分段三次多项式(称为三次样条函数)是最常用的一种样条函数,它在每一段上,是(不同的)三次多项式,且满足整条曲线的斜率与曲率为连续等条件。在数学上,利用样条*放样,相当于进行插值*。它的基本思想是,把整个区间分段,各段分别用低次多项式逼近一个函数,...... 2、按照一定的光滑性要求,连接起来的分段多项式。分段三次多项式(称为三次样条函数)是最常用的一种样条函数,它在每一段上,是(不同的)三次多项式,且满足整条曲线的斜率与曲率为连续等条件。 查看全文 样条函数 在学术文献中的解释 1、样条插值约束属于强约束条件相应的拟合函数称为样条函数.与其它拟合方法相比样条函数更具一般性和真实性.下面是我们在该系统中所用的曲面拟会方法这两种方法在近年的应用中效果良好 文献来源 2、当函数为分段m次多项式,而且在分点上直至(m-1)阶导数连续时,通常称为样条函数,简称样条.采用三次样条函数替代边界元数值计算中常用的分段多项式插值作函数逼近,其优点非常突出 文献来源 3、对样条曲线进行数学模拟得到的函数便称为“样条函数”.”“样条函数”亦称“样条插值函数”实际上是一种改进的分段插值函数它不同一般分段插值函数在各分段点处存在间断而是能保持整个曲线的光滑并且保留分段插值函数的低次性特点 文献来源 4、样条函数就是指一类可以逼近δ函数的样条函数:δ(x)=0　　x≠0.∞x=0.(1)k阶δ样条函数Ωk(x)可以写成Ωk(x)=k+1j(-1)jCjk+1(x+k+12-j)k+k 文献来源 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) * 插值方法(interp,spline) 拟合方法(polyfit,csaps) 数值微分(polyder, fnder) 数值积分(quad, quadl, fnint) 插值、拟合、数值微分、数值积分在化工计算中的作用 表格式物性数据