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）图三