clc;clear;close
n=1:-0.01:0;
r=0.5;
k=0.20826;t0=3.27216;y0=0.82581;
t1=t0;y1=-r*y0;t=t0+n*pi;
%t=[t,t+pi/4];
c1=1+k*sin(t1);c2=y1+k*cos(t1);
x=c1+c2*(t-t0)-k*sin(t);
y=c2-k*cos(t);
plot(x,y)
xlabel('Displacement')
ylabel('velocity')
axis([-1 1 -1 1])
上面那条曲线是不可能的,因为不管t为何值y=c2-k*cos(t)= -0.6194- 0.2083cos(t)<0的.怎么可能大于0呢