文档

mkpp

做分段多项式

描述

例子

= mkpp (休息时间系数建立一个分段多项式从它的断点和系数。使用ppval求分段多项式在特定点上的值,或者unmkpp提取分段多项式的细节。

= mkpp (休息时间系数d指定分段多项式是向量值的,这样它的每个系数的值都是长度的向量d

例子

全部折叠

创建一个分段多项式,在区间[0,4]中有一个三次多项式,在区间[4,10]中有一个二次多项式,在区间[10,15]中有一个四次多项式。

break = [0 4 10 15];Coefs = [0 1 -1 1 1;0 0 1 -2 53;1 6 1 4 77];Pp = mkpp(break,coefs)
页=带字段的结构:形式:'pp'打破:[0 4 10 15]coefs: [3x5 double]件:3顺序:5 dim: 1

在区间[0,15]的多个点处计算分段多项式,并绘制结果。在多项式相交的断点处画出垂直虚线。

Xq = 0:0.01:15;Plot (xq,ppval(pp,xq)) line([4 4],ylim,“线型”“——”“颜色”“k”) line([10 10],ylim,“线型”“——”“颜色”“k”

创建并绘制一个分段多项式,四个区间在两个二次多项式之间交替。

前两个子图显示了一个二次多项式,其负移到区间[-8,-4]和[-4,0]。多项式是

1 - x 2 - 1 2 - x 2 4 + x

第三个子图显示了一个分段多项式,由这两个二次段在四个区间交替构造。添加垂直线来表示多项式相交的点。

Subplot (2,2,1) cc = [-1/4 1 0];Pp1 = mkpp([-8 -4],cc);Xx1 = -8:0.1:-4;情节(民ppval (pp1民),“k -”) subplot(2,2,2) pp2 = mkpp([-4 0],-cc);Xx2 = -4:0.1:0;情节(xx2 ppval (xx2 pp2),“k -”次要情节(2,1,2)页= mkpp ([8 4 0 4 8], [cc、cc、cc; cc]);Xx = -8:0.1:8;情节(xx ppval (pp、xx),“k -”)举行线(4 [4],ylim,“线型”“——”) line([0 0],ylim,“线型”“——”) line([4 4],ylim,“线型”“——”)举行

输入参数

全部折叠

断点,指定为长度向量L + 1的开始和结束的严格递增元素l间隔。

数据类型:|

多项式系数,用an表示l——- - - - - -k第i行矩阵系数(我,:)包含某一阶局部系数的k第i个区间的多项式,[休息(我),优惠(i + 1).换句话说,多项式是系数(我,1)* (X-breaks(我))^ (k - 1) +系数(我,2)* (X-breaks(我))^ (k-2) +…+ coefs(i,k-1)*(X-breaks(i)) + coefs(i,k)

数据类型:|

维度,指定为整数的标量或向量。指定d表示分段多项式具有大小的系数值d

数据类型:|

输出参数

全部折叠

分段多项式,以结构形式返回。元素使用此结构ppval函数在一个或多个查询点处计算分段多项式。该结构具有这些字段。

描述
形式

“页”分段多项式

休息时间

长度向量L + 1的开始和结束的严格递增元素l时间间隔

系数

l——- - - - - -k每一行的矩阵系数(我,:)包含某一阶局部系数的k第i个区间的多项式,[休息(我),优惠(i + 1)

件数,l

订单

多项式的阶数

昏暗的

目标的维度

因为多项式系数系数为每个区间的局部系数,则必须减去对应结区间的下端点才能在常规多项式方程中使用这些系数。换句话说,对于系数(a, b, c, d)在间隔中(x1, x2),对应的多项式为

f x 一个 x x 1 3. + b x x 1 2 + c x x 1 + d

扩展功能

另请参阅

|||

R2006a之前介绍