MATLAB与工程应用-第7章-动力学与振动.pptVIP

第七章 动力学与振动;7.1轨迹 举例说明:重力场中有两个物体,其中质量为m2的物体固定,而质量为m1的物体绕m2做平面圆周运动.做圆周运动的m1物体的轨道半径用变量r表示,角度用变量a表示. m1 r a m2 两物体系统;卫星绕地球转动时，m2等于地球的质量，m1等于卫星的质量，r为卫星球心与地球球心间的距离。其运动轨迹由下列方程组决定： 式中： ，其中t是时间变量，p为物体在地球表面做圆周运动的周期。在地球表面，r=6.373x106 m。 ;用龙格—库塔法可以实现求解： 引入新状态变量： 带入前面的微分方程组，可得四个一阶微分方程。;建立函数文件orbit.m function xd=orbit(t,x) xd=[x(2);x(1)*x(4)^2-4.0*pi^2/x(1)^2;