matlab - 非重叠圆的螺旋

Question

我想创建一个从不相互重叠的圆形标记螺旋。这是我到目前为止得到的，但它与第一​​个标记重叠，最后一个标记彼此相距太远。

t = pi : pi/20 : 20*pi;
t = asind(1./t);
r = t;
x = r .* cos(t);
y = r .* sin(t);
plot(x,y,'o-');
axis equal; hold on 

第二个图中显示了未将 t 重新定义为 asinf(1/t) 的绘图，如下所示。

t = pi : pi/20 : 20*pi;
r = t;
x = r .* cos(t);
y = r .* sin(t);
plot(x,y,'o-');

关于角度 t 的间距必须如何实现标记不重叠的任何想法？

标签： matlabplotspiral

Answer 1

您可以近似弧长，大大简化Gilles-Phillipe 的解。这是一种简化，这意味着标记之间的距离并非处处相同。然而，距离相当一致，尤其是更远的地方。

这里的近似是假设螺旋在局部是一个圆。然后弧长r*dt在螺旋中的位置r与原点有一段距离，以改变dt弧度角。

我们现在不再需要求解符号方程。我在一个循环中编写了代码。我确信可以将它矢量化，使整个事情变成两行代码，但我将把它作为练习留给读者。

这是代码：

d = 1; % step size
q = 1/(2*pi); % spiral constant -- radius grows by q every 1 radian turn
N = 300; % number of points

t = 0; % initial angle
r = d; % initial radius
p = zeros(100,2);
p(1,:) = [r*cos(t),r*sin(t)]; % first point
for ii=2:N
   dt = d/r;
   t = t+dt;
   r = r+dt*q;
   p(ii,:) = [r*cos(t),r*sin(t)];
end

clf
plot(p(:,1),p(:,2),'o-') 
axis equal