这样y的分量的个数比x的分量的个数多2。例如,在 例1中 x=[0123] y=[0.200.52.01.5-1] 周期边界条件时,无需指定边界条件值 在计算一个三次样条表达式的时候,必须将pp形中 不同域提取出来进行计算,这个过程可以由函数 unmkpp 完成,该函数的使用方法为: Breaks, coefs, npolys, ncoefs, dim]=unmkpp(pp)
这样y的分量的个数比x的分量的个数多2。例如,在 例1中, x=[0 1 2 3] y=[0.2 0 0.5 2.0 1.5 -1] 周期边界条件时,无需指定边界条件值。 在计算一个三次样条表达式的时候,必须将pp形中 不同域提取出来进行计算,这个过程可以由函数unmkpp 完成,该函数的使用方法为: [breaks,coefs,npolys,ncoefs,dim]=unmkpp(pp)
其中输入变量p是样条插值函数 cape的输出变量, unmkppl的输出变量有5个: breaks, coefs, polys, ncoefs, dim breaks:包含了插值节点 coefs:是一个矩阵,其第i行是第个多项式的系数:an,a1,a2a3即 (x)=a0(x-x)+a1(x-x1)2+a2(x-x1)+ 中的系数。 npoh:是多项式的个数; ncoefs:是每个多项式系数的个数; dim:是样条的维数
其中输入变量pp是样条插值函数csape的输出变量, unmkpp的输出变量有5个: breaks,coefs,npolys,ncoefs,dim。 0 1 2 3 3 2 0 1 1 1 2 1 3 , , , . ( ) ( ) ( ) ( ) : : dim : i i i i a a a a s x a x x a x x a x x a npolys ncoefs = − + − + − + − − − breaks:包含了插值节点; coefs:是一个矩阵,其第i行是第i个多项式的系数: 即 中的系数。 是多项式的个数; 是每个多项式系数的个数; 是样条的维数
其中输入变量pp是样条插值函数 cape的输出变量, unmkpp的输出变量有5个: breaks, coefs, polys, ncoefs, dim breaks:包含了插值节点; coefs:是一个矩阵,其第i行是第i个多项式的系数:ana,a2a3即 (x)=an(x-x21)2+a1(x-x1)2+a2(x-x1) 中的系数。 polys:是多项式的个数; ncoefs:是每个多项式系数的个数; dim:是样条的维数
其中输入变量pp是样条插值函数csape的输出变量, unmkpp的输出变量有5个: breaks,coefs,npolys,ncoefs,dim。 0 1 2 3 3 2 0 1 1 1 2 1 3 , , , . ( ) ( ) ( ) ( ) : : dim : i i i i a a a a s x a x x a x x a x x a npolys ncoefs = − + − + − + − − − breaks:包含了插值节点; coefs:是一个矩阵,其第i行是第i个多项式的系数: 即 中的系数。 是多项式的个数; 是每个多项式系数的个数; 是样条的维数
运行结果如下: pp- form: pp breaks:[0123 coefs:[3×4 double pleces order lm
运行结果如下: pp= form:'pp' breaks:[0 1 2 3] coefs:[3 4 double] pieces:3 order:4 dim:1
reaks- polys 0123 coefs- chefs 0.2800-0.18000.2000 4 1.04001.26001.28000.5000 dim= 0.6800-1.86000.68002.0000 因此所求的三次样条函数为 0.28x3-0.18x2+0.2x x∈[0,1] s()=-104x-1)+126(x-12+128(x-1)+05x∈2 0.68(x-2)-1.86(x-2)2+0.68(x-2)+20, x∈
3 2 3 2 3 2 0.18 0.2 , [0,1] 1.04( 1) 1.26( 1) 1.28( 1) 0.5, [1,2] ( 2) 1.86( 2) 0.68( 2) 2.0, [2,3] x x x x s x x x x x x x x x − + − − + − + − + − − − + − + 因此所求的三次样条函数为 0.28 ( )= 0.68 breaks= 0 1 2 3 coefs= 0.2800 -0.1800 0.2000 0 -1.0400 1.2600 1.2800 0.5000 0.6800 -1.8600 0.6800 2.0000 npolys= 3 ncoefs= 4 dim= 1