Matlab-质点的运动

1、源码

%质点运动的位置,速度和加速度
clear                                  %清除变量
tm=3;                                  %最大时间
dt=1;                                  %时间间隔
t1=0:tm;                               %以1秒为间隔的时间向量
x1=5*t1-t1.^3;                         %位置坐标
v1=diff(x1)/dt;                        %速度的平均值
v1=[5,v1];                             %补充初速度
a1=diff(v1)/dt;                        %加速度的平均值
a1=[0,a1];                             %补充初加速度

dt=0.1;                                %时间间隔
t2=0:dt:tm;                            %以0.1秒为间隔的时间向量
x2=5*t2-t2.^3;                         %位移
v2=diff(x2)/dt;                        %速度的平均值
v2=[5,v2];                             %补充初速度
a2=diff(v2)/dt;                        %加速度的平均值
a2=[0,a2];                             %补充初加速度

dt=0.01;                               %时间间隔
t3=0:dt:tm;                            %以0.01秒为间隔的时间向量
x3=5*t3-t3.^3;                         %位移
v3=diff(x3)/dt;                        %速度的平均值
v3=[5,v3];                             %补充初速度
a3=diff(v3)/dt;                        %加速度的平均值
a3=[0,a3];                             %补充初加速度

x=5*t2-t2.^3;                          %坐标的瞬时值
v=5-3*t2.^2;                           %速度的瞬时值
a=-6*t2;                               %加速度的瞬时值
figure                                 %创建图形窗口
plot(t1,x1,'-*',t2,x2,'-+',t3,x3,'.-',t2,x,'-o')%画位置坐标曲线族
legend('1','0.1','0.01','瞬时值')      %插入时间差图例
grid on                                %加网格
title('质点的位置','FontSize',16)      %标题
xlabel('\itt\rm/s','FontSize',16)      %横坐标标签
ylabel('\itx\rm/m','FontSize',16)      %纵坐标标签
figure                                 %创建图形窗口
plot(t1,v1,'-*',t2,v2,'-+',t3,v3,'.-',t2,v,'-o')%画速度曲线族
legend('1','0.1','0.01','瞬时值')      %插入时间差图例
grid on                                %加网格
title('质点的速度','FontSize',16)      %标题
xlabel('\itt\rm/s','FontSize',16)      %横坐标标签
ylabel('\itv\rm/m\cdots^-^1','FontSize',16)%纵坐标标签
figure                                 %创建图形窗口
plot(t1,a1,'-*',t2,a2,'-+',t3,a3,'.-',t2,a,'-o')%画加速度曲线族
legend('1','0.1','0.01','瞬时值')      %插入时间差图例
grid on                                %加网格
title('质点的加速度','FontSize',16)    %标题
xlabel('\itt\rm/s','FontSize',16)      %横坐标标签
ylabel('\ita\rm/m\cdots^-^2','FontSize',16)%纵坐标标签

2、运行结果

（1）图一

（2）图二

（3）图三