向量值函数

工具箱支持金宝app向量值样条函数。例如,如果你想要一个花键曲线通过给定的平面点,然后下面的代码定义了一些数据,然后创建和情节这样的样条曲线,使用弦长参数化和三次样条插值not-a-knot结束条件。

x = [19 43 62 88 114 120 130 129 113 76 135 182 232 298…348 386 420 456 471 485 463 444 414 348 275 192 106…30 48 83 107 110 109 92 66 45 23日22 30 40 55 55 52 34 20 16];y = [306 272 240 215 218 237 275 310 368 424 425 427 428…397 353 302 259 200 148 105 77 47 28 17 10 12 23 41 43……77 96 133 155 164 157 148 142 162 181 187 192 202 217 245 266 303);xy = (x, y);df = diff (xy, 1、2);t = cumsum ([0,√(1) * (df。* df))));简历= csapi (t, xy); fnplt(cv), hold on, plot(x,y,'o'), hold off

如果你想知道的面积曲线包围,你想评价积分,曲线上的点对应的参数值。样条曲线简历构造,可以完全在一个(有点复杂)命令:

面积= diff (fnval (fnint (…fncmb (fncmb(简历,[0 1]),‘*’,曾经(fncmb(简历,[1 0])))…),fnbrk(简历,“间隔”)));

来解释,y = fncmb(简历,[0 1])挑出曲线的第二个组件简历,Dx =曾经(fncmb(简历,[1 0]))提供第一个组件的导数,yDx = fncmb (y,‘*’, Dx)构造他们的逐点的产品。然后IyDx = fnint (yDx)构造的不定积分yDx最后,diff (fnval (IyDx fnbrk(简历,“间隔”)))计算不定积分的基本区间的端点,然后将第二个从第一个值的差别,因此得到的定积分yDx对其基本的间隔。取决于封闭区域向右或向左的曲线点的旅行增加参数,结果是积极的还是消极的。

此外,所有的值Y(如果有的话)的观点(X, Y)在样条曲线简历只是构造可以通过以下(有点复杂)命令:

X = 250;%定义值X Y = fnval (fncmb(简历,[0 1]),…意思是(fnzeros (fncmb (fncmb(简历,[1 0]),“-”,X))))

解释:1 x = fncmb(简历,[0])挑出曲线的第一个组件简历;”——“xmX = fncmb (x, x)组件的翻译X;t =意味着(fnzeros (xmX))提供了所有的参数值xmX是零,即。,for which the first component of the curve equalsX;y = fncmb(简历,[0,1])挑出曲线的第二个组件简历;最后,Y, Y = fnval (t)网站评估第二组件在这些参数曲线的第一个组件简历=X。

向量值函数的使用的另一个例子,假设您已经解决了粒子的运动方程在一些指定的力场在平面上,获得,在离散的时间,位置以及速度存储在4向量你如果,在标准的方式,你已经解决了数值相当于一阶系统。然后下面的语句,它使用立方埃尔米特插值,将产生一个阴谋的粒子路径:fnplt (spapi (augknt (t 4 2), t,重塑(z, 2 2 * n)))。