clc
clear
close all
step=1;
alphy=0.082;%变换角度
beta=0;%初始方向角
figure(1)
hold on
axis equal
title('欧拉螺线')
alphy_lin=0;
Q1=[0 0];%初始点坐标
k=0;
while true
Q2(1)=Q1(1)+step*cos(beta+alphy_lin);
Q2(2)=Q1(2)+step*sin(beta+alphy_lin);
plot([Q1(1),Q2(1)],[Q1(2),Q2(2)],'b')
drawnow
alphy_lin=alphy_lin+alphy;
beta=beta+alphy_lin;
if alphy_lin>=2*pi
k=k+1;
end
alphy_lin=mod(alphy_lin,2*pi);
beta=mod(beta,2*pi);
Q1=Q2;
if k>20
break
end
end
%% qq115075524
