Question

我的代码现在

% 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的定义，但我无法弄清楚如何...感谢您的帮助！

Answer 1

将此作为您的功能：

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点。

在Matlab中更改参数样条曲线的网格点

1 个答案: