我有一个矩阵值函数,当x到1时,我试图求它的极限。
因此,在这个例子中,我有三个矩阵v1-3,分别表示在0.85,0.9,0.99的采样值。我现在所做的工作效率很低,主要有以下几点:
for i=1:101
for j = 1:160
v_splined = spline([0.85,0.9,0.99], [v1(i,j), v2(i,j), v3(i,j)], [1]);
end
end
必须有更好、更有效的办法来做到这一点。特别是当足够快的时候,我将面对v将是4-5维向量的情况。
谢谢!
发布于 2014-04-24 17:41:12
免责声明:天真地推断出有风险的业务,自己承担风险。
我要说的是
如果你想使用拉格朗日多项式来形成抛物线,你的结果是:
v_splined = v1*(1-.9)*(1-.99)/( (.85-.9)*(.85-.99) ) ...
+v2*(1-.85)*(1-.99)/( (.9-.85)*(.9-.99) ) ...
+v3*(1-.85)*(1-.9)/( (.99-.85)*(.99-.9) );
我把这个不简化的,所以你可以看到它是如何来自拉格朗日多项式,但显然简化是容易的。还请注意,这消除了循环的需要。
https://stackoverflow.com/questions/23279752
复制相似问题