16.2　运动学方程求解过程的MATLAB实现

运动学方程（16-1）和方程（16-2）为二阶常微分方程，可通过ode45函数求解。对应倾倒过程的第一阶段（无滑动的定轴转动）编制自定义函数function k=qingdao_fun1（t,y,flag,j0,oc,v,r,m,g,1,p,H1,H2,c），并存为文件qingdao_fun1 .m，以方便GUI主程序调用，代码如下：

对应倾倒第二阶段（又滑动又旋转的平面运动）编制自定义函数function k=qingdao_fun2（t,y,flag,j0,jc,oc,v,r,m,g,1,p,u,H1,H2,c），代码如下：