实用 MATLAB教学资料-5-6.pptVIP

for i=1:n%构造Hermite插值多项式 la=1; lp=0.0; for j=1:n if(j~=i) la=la*(x-xe(j))/(xe(i)-xe(j)); lp=lp+1/(xe(i)-xe(j)); end end temp1=1-2*(x-xe(i))*lp; temp2=ye(i)*temp1*la^2; temp3=dye(i)*(x-xe(i))*la^2; f=f+temp2+temp3; end %至此，多项式构造完毕 f=simplify(f)%多项式表达形式化简 pretty(f)%以习惯的表达形式输出 if(nargin==4) %如果输入向量为4个，那么除了输出构造的多项式之外，还要计算插值 f=subs(f,x,xi);% 用xi替换f 表达式中的x end 参考：宋叶志、贾东永 编著. MATLAB数值分析与应用，机械工业出版社，2009年7月第1版. 已知数据如下，请分别用Lagrange插值、分段线性插值、分段三次Hermite插值和三次样条插值构造插值函数，并绘出节点以及插值函数曲线。 作 业 x -5 -4 -3 -2 -1 0 f(x) 0.038462 0.058824 0.1 0.2 0.5 1 x 1 2 3 4 5 ? f(x) 0.5 0.2 0.1 0.058824 0.038462 ? 如果不添加边界条件，则边界为not-a-knot 如果添加边界条件，只能是类型I ; 边界条件书写格式[ b1 ye b2 ] 　三、the cubic spline interpolation 　　　三次样条插值 2、spline 函数 格式1 yi=spline(xe, ye, xi) spline函数 求xi处的函数值， 不能给出插值多项式的信息 pp=spline(xe, ye) yi=ppval( pp, xi ) spline函数 格式2 格式3 yi=ppval( spline(xe, ye), xi ) spline函数 不能给出插值多项式的信息 ppval 也可以用 fnval替换 例4 已知某液相反应 浓度随时间变化的实验数据，用spline计算x=0.1和x=0.4min的浓度。 x / min y / g·L-1 0 5.19 0.2 3.77 0.6 2.3 1 1.57 2 0.8 5 0.25 10 0.094 x / min y / g.L- x / min y / g·L-1 y/ g·L-1min-1 0 5.19 -9.45 0.2 3.77 0.6 2.3 1 1.57 2 0.8 5 0.25 10 0.094 0 例4 已知某液相反应 浓度随时间变化的实验数据，用spline计算x=0.1和x=0.4min的浓度。 第3章 插值 3.3 MATLAB插值函数曲线的绘制 一、fnplt 函数 先构造插值函数 y= interp1(xe, ye, ‘linear’,‘pp’) y= pchip(xe, ye) y=csape(xe, ye, conds , [b1,b2]) y=spline(xe, [b1 ye b2]) fnplt(y) 再绘制函数曲线 cscvn构造自然（variational, or natural）边界条件的三次样条函数 插值节点可以重复（绘制封闭曲线） y = cscvn(points) fnplt(y) 其中points是插值节点 cscvn 与 fnplt 联用 例5 画圆圈 cscvn 与 fnplt 联用 二、getcurve 功能：实时画图 原理：用三次样条曲线（cscvn函数）连接散点 使用方法： 在命令窗口输入getcurve，进入绘图窗口。 在网格区，点击鼠标一下，形成一个点。绘图过程中，这些点由折线连接。在网格外区域点击一下结束绘图。 二、getcurve 例6 画圆圈 getcurve 第3章 插值 3.4 MATLAB样条GUI工具箱 Graphical User Interface（GUI） 图形用户界面进行构造插值函数 在command窗口输入splinetool 点击import data，输入节点数据 样条曲线信息： 增加曲线Replicate 删除曲线Delete 更改曲线名称Rename 方法信息： 选择样条方法（插值和拟合两类） 设置边界条件类型和数值 节点信息： 增减Edit 改变值的大小 显示构造的函数曲线： 更改轴名 Edit 显示