靠在墙角的木棍

 高中的时候我们都做过靠在墙角的一根棍子，当一段以匀速运动时，求另一端的速度，此代码就是求解这一问题。

while 1
clc,clear;
AB=input('请输入AB的长度')
if AB==000  %如果AB等于000，则退出程序
break
end  
va=input('请输入A端的速度')
tmax=AB/va;%%最大时间
n=1;
t0=tmax/200;
for t=0:t0:tmax
x(n)=AB-t*va;%求得每刻的横坐标
y(n)=(AB^2-x(n).^2).^0.5;%求得每刻的纵坐标
fine(n)=acos((x(n))/AB);%求得每刻fine的值
vb(n)=va*cot(fine(n)); %求出vb的速度
vab(n)=va/sin(fine(n));%求出vab的速度
oumiga(n)=vab(n)/AB;%求出角速度
n=n+1;
end
figure
line([0,0],[0,AB+2],'Color','r','LineWidth',2)%A到端点的连线
line([0,AB+2],[0,0],'Color','r','LineWidth',2)%B到端点的连线
h1 = line('Color','b','Marker','.','MarkerSize',20);%画点A
h2 = line('Color','b','Marker','.','MarkerSize',20);%画点B
for i=1:length(x)
set(h1,'xdata',x(i),'ydata',0);%设置A的初始坐标
set(h2,'xdata',0,'ydata',y(i));%设置B的初始坐标
l(i)=line([x(i),0],[0,y(i)],'linewidth',1);
drawnow;                   %刷新屏幕
delete(l(i))          %删除生成的直线        
pause(t0)      %每t0 秒 展示图片
end
t=0:t0:tmax;
figure
subplot(2,1,1)
plot(t,vb);%绘制vb与t的图像
ylabel('速度（m/s）')
xlabel('时间（s）')
title('Vb的速度')
subplot(2,1,2)
plot(t,oumiga)%绘制角速度与t的图像
ylabel('角速度(rad/s)')
xlabel('时间（s）')
title('角速度')
end