我的代码现在
% Create some example points x and y
t = pi*[0:.05:1,1.1,1.2:.02:2]; a = 3/2*sqrt(2);
for i=1:size(t,2)
x(i) = a*sqrt(2)*cos(t(i))/(sin(t(i)).^2+1);
y(i) = a*sqrt(2)*cos(t(i))*sin(t(i))/(sin(t(i))^2+1);
end
请注意:点(x_i | y_i)不一定是等距的,这就是为什么t就是这样创建的。也不应该在进一步的代码中使用t,因为我的实际问题是未知的,我最后得到一堆x,y和z值。在这个例子中,我把它减少到2D。
现在我正在为x和y值创建ParametricSplines
% Spline
n=100; [x_t, y_t, tt] = ParametricSpline(x, y, n);
xref = ppval(x_t, tt); yref = ppval(y_t, tt);
具有功能
function [ x_t, y_t, t_t ] = ParametricSpline(x,y,n)
m = length(x);
t = zeros(m, 1);
for i=2:m
arc_length = sqrt((x(i)-x(i-1))^2 + (y(i)-y(i-1))^2);
t(i) = t(i-1) + arc_length;
end
t=t./t(length(t));
x_t = spline(t, x);
y_t = spline(t, y);
t_t = linspace(0,1,n);
end
生成的情节
plot(x,y,'ob',...
xref,yref,'xk',...
xref,yref,'-r'),...
axis equal;
如下所示:Plot Spline
问题:
如何更改代码,以便我总是直接在最初给定的点(xref_i|yref_i)
(显示为蓝色O)上得到一个结果点(x_j|y_j)
(在图中显示为黑色X)还有(x_j|y_j)
和(x_j+1|y_j+1)
之间的n个点?
E.g。如果n = 2,我想得到以下结果:
(xref_1|yref_1) = (x_1|y_1)
(xref_2|yref_2)
(xref_3|yref_3)
(xref_4|yref_4) = (x_2|y_2)
(xref_5|yref_5)
[...]
我想我唯一需要的是更改tt
的定义,但我无法弄清楚如何...感谢您的帮助!
答案 0 :(得分:0)
将此作为您的功能:
function [ x_t, y_t, tt ] = ParametricSpline(x,y,nt)
arc_length = 0;
n = length(x);
t = zeros(n, 1);
mul_p = linspace(0,1,nt+2)';
mul_p = mul_p(2:end);
tt = t(1);
for i=2:n
arc_length = sqrt((x(i)-x(i-1))^2 + (y(i)-y(i-1))^2);
t(i) = t(i-1) + arc_length;
add_points = mul_p * arc_length + t(i-1);
tt = [tt ; add_points];
end
t=t./t(end);
tt = tt./tt(end);
x_t = spline(t, x);
y_t = spline(t, y);
end
精华: 你必须以与距离向量t相同的方式构造tt,并在其间添加额外的nt点。