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的区别;二维插值：interp2;二维插值函数的使用;拟合简介;最小二乘多项式拟合：polyfit;多项式次数对拟合效果的影响;最小二乘法拟合生成样条曲线;功能：平滑生成三次样条函数，即对于数据(xi,yi)，所求的三次样条函数y=f(x)满足 ：;功能：用最小二乘法拟合生成B样条曲线，即对于数据(xi,yi)，所求k次样条函数y=f(x)满足 ：;功能：平滑生成B样条函数 ;采用样条拟合函数csaps进行拟合，比较p的取值对于拟合效果的影响 ;数值微分;Matlab数值微分实现方法;函数diff;函数gradient;例题;广义积分;样条函数在数值积分与微分中的应用;反应器停留时间分布的混合特性;;微分法进行动力学数据分析;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);样条函数：照一定的光滑性要求,连接起来的分段多项式。分段三次多项式(称为三次样条函数)是最常用的一种样条函数,它在每一段上,是(不同的)三次多项式,且满足整条曲线的斜率与曲率为连续等条件。在数学上,利用样条*放样,相当于进行插值*。它的基本思想是,把整个区间分段,各段分别用低次多项式逼近一个函数,...... 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 文献来源