matlab3.ppt

应用数学系分析教研室 五、样条工具箱（splines toolbox） 三次插值样条函数 PP形式的样条函数的构造及操作 B形式样条函数的构造和应用 其他函数 举例 Matlab5中的样条工具箱为学习样条和使用样条函数提供了一个理想的软件环境.它提供了函数最常用的两种表示形式:B形式和PP形式(Piecewise Polynomial). B形式,即以B样条函数基底的线性组合来表示样条函数; PP形式,即以分段多项式来表示样条函数. 工具箱中的操作函数具体分为以下几类: 对三次插值样条函数的操作; 对pp形式样条函数的操作; 对B形式样条函数的操作; 对张量样条函数的操作; 其他. 3.1 三次插值样条函数 3.1.1 三次插值样条函数的定义 给定区间 [a,b] 的一个划分 和区间[a,b]上的一个函数f(x) 定义1 若函数s(x) 满足下列条件: 1)一致通过n+1个型值点 即 (i=0,1,2,…,n) 2)二阶连续,即 3)三次分段,即在每一个小区间 上均为三次多项式. 这样的 s(x)称为[a,b]上以 为节点的三次插值样条函数,其几何图形称为三次样条曲线. 由定义可知,为确定 s(x),还需附加2个边界条件: 3.1.2 三次插值样条函数的构造 记节点处的一阶导数为 记节点处的二阶导数为 1.用节点处的一阶导数构造三次插值样条函数。由插值条件可得,在 上 可表示为 记 对于给定端点斜率： 可解方程组 求得 对于给定端点的二阶导数，即 可以解方程组 2 用节点处的二阶导数构造三次插值样条函数 同样对于给定的三种不同条件:端点斜率.二阶导数,周期特性,有相应的方程组,可解得 3.1.3 工具箱中关于三次插值样条的函数 样条工具箱提供以下5种函数来构造不同边界条件和要求的三次插值样条函数. (一) csape 功能:构造各种边界条件下的三次插值样条函数 格式:pp=scape(x,y,[,conds[,valconds]]) 说明:(x,y)是插值点的序列,pp为指定conds条件下以(x,y)为插值点所返回的pp形式的三次样条函数. (二) csapi 功能:构造‘not-a-knot’边界条件下的三次样条函数. 格式:values=csapi(x,y,xx) pp=csapi(x,y) 说明:values返回以(x,y)为插值点序列,函数在向量xx处的值, pp返回函数. 举例:多变元 图形: (三) csaps 功能:构造光滑的三次插值样条函数. 格式:values=csaps(x,y,p,xx) pp=csaps(x,y,p) 说明:p=0时,以最小方差直线拟合数据,光滑程度最低; p=1时,以‘variational’为边界构造函数,光滑程度最高. 举例: 一般情况 fnplt(cscvn([0 1 1 0 –1 –1 0 0; 0 0 1 2 1 0 –1 –2])); 周期情况 fnplt(cscvn([0 .8 .9 0 0 -.9 -.8 0;.5 1 0 –1 0 1 .5])); 图形: 一般情况 周期情况 (五) getcurve 功能:交互式的创造三次样条函数. 格式:[xy,spcv]=getcurve 3.2 PP形式的样条函数的构造及操作 3.2.1 分段多项式形式的样条函数 定义2 给定区间 的一个划分 当 满足下列两个条件： 1）在每个子区间 上， 是次数不高于n的多项式，但至少有一个子区间上为n次多项式； 2） 在区间 上存在n-1阶的连续导数，即 则称 为关于划分的n次样条函数. 由定义可知,零次样条是阶梯函数,一次样条是折线函数. 3.2.