time=input('intervaltime='); %输入时间步,以 0.1s 为单位,每隔 time 时间画一条曲线
controllfactor=input('controllfactor='); %画曲线的条数
k=input('敏感系数 k=');
k1=input('k1=');
k2=input('k2=');
k3=input('k3=');
t=1;
L=1500;
x=0:15:1485;
x(1)=10;
temp=x;
f=linspace(0,0,100);
v=linspace(0,0,100);
v=v+tanh(15);
temp_v=v;
a=1:100;
t_1=[];
pos_1=[];
T=linspace(0,0,100);
tt=0.1;
pp=1;
while pp<=controllfactor
while t<=time*pp
for j=1:100
if j==100
xd(j)=temp(1)-temp(100);
else
xd(j)=temp(j+1)-temp(j);
end
if xd(j)<0
xd(j)=L+xd(j);
end
% a(j)=tanh(xd(j)-2)+tanh(2)-temp_v(j);
if j==100
a(j)=k*((6.75+7.91*tanh(0.13*(xd(j)-5)-1.57))-temp_v(j))+k1*(temp_v(1)-
temp_v(j))+k2*(temp_v(2)-temp_v(1))+k3*(temp_v(3)-temp_v(2));