文件

fntlr.

泰勒系数

句法

泰勒= fntlr(f,dorder,x)

描述

泰勒= fntlr(f,dorder,x)返回无通知的泰勒系数,直至给定的顺序dorder.在给定的X,描述的功能F

对于单变量函数和标量X,这是向量

T. F dorder. X 的) = [ F X 的) ; D. F X 的) ; ...... ; D. dorder. - 1 F X 的) ]

如果,更一般地,函数FD.- 夸耀D> 1甚至产品(d)> 1和/或是M.- 对于一些人来说M> 1, 然后dorder.预计是一个M.- 正整数的vector,X预计将成为一个矩阵M.行,并且,在这种情况下,输出大小[PROD(d)* prod(dorder),size(x,2)],它的第j栏包含

T. F dorder. X j 的) 的) 一世 1 ...... 一世 M. 的) = D. 1 一世 1 - 1 ...... D. M. 一世 M. - 1 F X 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。

也可以看看

|

这个主题有用吗?