泰勒系数
泰勒= fntlr(f,dorder,x)
泰勒= fntlr(f,dorder,x)
返回无通知的泰勒系数,直至给定的顺序dorder.
在给定的X
,描述的功能F
。
对于单变量函数和标量X
,这是向量
如果,更一般地,函数F
是D.
- 夸耀D> 1
甚至产品(d)> 1
和/或是M.
- 对于一些人来说M> 1
, 然后dorder.
预计是一个M.
- 正整数的vector,X
预计将成为一个矩阵M.
行,并且,在这种情况下,输出大小[PROD(d)* prod(dorder),size(x,2)]
,它的第j栏包含
为了I1 = 1:dorder(1)
,......,IM = 1:dorder(m)
。这里,D.一世F是部分衍生的F关于它的一世争论。
如果F
包含单变量功能和X
是标量还是1行矩阵fntlr(f,3,x)
生成与陈述相同的输出
df = fnder(f);[fnval(f,x);fnval(df,x);fnval(fnder(df),x)];
作为一个更复杂的例子,在21个相同间隔点的合理样条点看泰勒向量,其图是单位圈:
ci = rsmak('cirction');在= fnbrk(ci,'温度');t = linspace(在(1)中,(2),21);t(结束)= [];v = fntlr(ci,3,t);
我们情节CI.
和点一样v(1:2,:)
,确实验证这些是单位圈的点。
fnplt(ci),保持on,plot(v(1,:),v(2,:),'o')
接下来,验证v(3:4,j)
是一个向圆圈切相切的矢量v(1:2,j)
,我们使用matlab®颤动
命令将相应的箭头添加到我们的绘图中:
quivive(v(1,:),v(2,:),v(3,:),v(4,:))
最后,怎么样v(5:6,:)
还这些是第二衍生物,我们通过以下添加相应的箭头颤动
命令,从而完成一个合理的花键的第一和第二衍生物给出一个圆圈。
quivive(v(1,:),v(2,:),v(5,:),v(6,:)),轴相等,阻止
一个合理的花键的第一和第二衍生物给出一个圆圈
现在,我们的曲线是一个圆圈,你可能会期待第二衍生箭头直接向该圈子的中心点,这确实如此CI.
已经使用arclenth作为其独立变量。由于使用的参数不是arclength,因此我们使用的公式示例:B形状样条近似到圆圈,计算曲线的曲率CI.
在这些选定的点。为了便于比较,我们切换到那里使用的变量,然后只需使用那里的命令。
dspt = v(3:4,:);ddspt = v(5:6,:);kappa = abs(dspt(1,:)。* ddspt(2,:) - dspt(2,:)。* ddspt(1,:))。/ ...(和(dspt. ^ 2))。^(3/2);MAX(ABS(κ-1))ANS = 2.2204E-016
数值答案是放心:在所有测试的点处,曲率为1到圆形OFF。