航天器动力学09-数据转换.pptVIP

补充内容 数据转换 ;数据转换;r;x1;x1;特例;三、坐标转换矩阵的性质;;其分量形式为;根据相对运动微分方程求解：;X;坐标转换关系;*－－根据非惯性系的动力学方程积分结果 o－－根据惯性系积分结果，再进行转换;clear; %清除内存 global m g omega program=------------------start------------------- m=1; %物体质量 g=9.8; %重力加速度 H=5; %自由下落高度 omega=10; %转盘角速度 radius=1; %转盘半径 a=1;b=2;c=3; vx=0;vy=-1;vz=10; %%%%%%%%%初始条件 %非惯性系中的运动 y0(1)=0; y0(2)=radius;y0(3)=0;y0(4)=vx;y0(5)=vy;y0(6)=vz; %惯性系中的运动 y0(7)=a;y0(8)=b+radius;y0(9)=c;y0(10)=vx-radius*omiga;y0(11)=vy;y0(12)=vz; %%%%%%%%%%%%%%%积分工作 alltime=2*vz/g; %% 积分时间 step_time=0.01; %% 积分步长 options = odeset(RelTol,1e-10,AbsTol,1e-12); %% 积分的误差选项 [time,iy]=ode45(rg_kt_two,[0:step_time:alltime],y0,options); %% 积分结果放在iy中;计算主程序2;function dy=rg_kt_two(t,y) global m g omega dy=zeros(12,1); %非惯性系中的运动 dy(1)=y(4); dy(2)=y(5); dy(3)=y(6); dy(4)=omega*omega*+y(1)+2*omega*y(5); dy(5)=omega*omega*+y(2)-2*omega*y(4); dy(6)=-g; %惯性系中的运动 dy(7)=y(10); dy(8)=y(11); dy(9)=y(12); dy(10)=0; dy(11)=0; dy(12)=-g; ;卫星的运动微分方程是在惯性坐标系OXYZ中列写的;在卫星轨道平面内的轨道坐标系Ox’y’z’中; Ω;O;O;欧拉角的方向余弦矩阵;下面把欧拉角与轨道根数联系起来;欧拉角的来历：;地球在绕太阳公转，同时还有自转。公转比较简单，可以看成是质心的运动；自转就复杂的多了，相当于刚体绕质心的姿态运动。;首先地球存在自转。; 如果地球没有自转，则单摆的运动应该是平面运动。但实际上发现单摆相对于地球的运动明显不在一个平面内，摆的投影点作曲线运动。; 根据天文学知识，地球除自转外，地轴本身也在惯性空间中运动，如同一个大的陀螺。 ;自转轴; 公元前二世纪，古希腊天文学家喜帕恰斯在编制一本包含1022颗恒星的星表时，首次发现了岁差现象。; 牛顿第一个指出产生岁差的原因：由??地球的外型成微扁又上下凹陷的「梨」形。因此当受到太阳、月亮的引力时，会产生不对称的引力，这个引力对于地球的运动形成微扰动力矩。;地轴移动的另一个结果是：北极星在不断变化。; 地球自转轴在空间作岁差运动的同时，还伴随有许多短周期变化。英国天文学家詹姆斯·布拉德利（James Bradley）在1748年分析了20年恒星位置的观测资料后，发现了章动现象。;由于地球自转轴的章动，地轴不是正好指向北极星，从地球上看，北极星作一小幅度的摆动，周期为18.6年左右。由于恒星的相对位置几千年才有微小的移动，由此可知北极星的这种短周期的摆动是由于地球的自转轴运动所引起的。;为什么地轴的这种“点头”运动叫“章动”呢？;月亮平均盈亏周期为 29.5306日地球公转平均周期为365.2422日 　　　　　? 　1.因为一年平均是365.2422日，每月平均是29.5306日 　　所以　365.2422÷29.5306＝12.368266471(月) 　　一年实际上比月历上的日期多出了0.368个月 　 2.又因为　0.3682664×1９=6.9970629～7 　　也就是大约在１９年之后会出多出７个月。 　 3.因此每19年必须有7个闰