建模教程正文

69.17W，42.00N，1752.013km，远月点与之相对，可求出其经纬度为110.83E，42.00S，距离月心1837.013km。与之相反，为西偏南2.44°。针对问题（2）我们分阶段进行研究。主减速阶段轨道是一条平缓的抛物线，建立630.75sg；快速调整阶段，运行轨道近似为直线下降，用能量守恒和平衡定理即可求出10.32sATLAB和VsualBasic26.95m2.3最后匀减速2.63,竖直软着陆时，根据平衡条件和运动学方程可求出应先自由下落20.34s,56.23s,0.45m，消耗燃料量和57.32s，嫦娥三号剩余质量1121kg2.22s742.491279kg。关键字二维动力学模型蚁群算法轨道设 三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机能够产生的预定着陆点为19.51°W，44.12°N，海拔为-2641m。MATLABVisualBasicm即GM

tmnmn月心连线所在直线为z轴，方向竖直向上。在此处研究的范围内，假设月面是绝对的平面，则可以规定此处平均半径所在平面为xoyzn((xx0)2yy0)2h 其中，x0、y0为近月点对月心投影的经纬度，hx0=69.17°W，与近月点相对，则远月点的位置为110.83°E，42.00°S，距离月心1837.013km。 v2=m2-1 已知月球平均半径为1737.013km，近月点距月球表面15km，远月点距月球表面即可求得：嫦娥三号近月点速度大小为1.69km/s，远月点速度大小为1.61km/s。69.17°W42.00°NW，44.12°N，则运用三角函数关系可得，tanη=0.04269,得η=2.44°，即嫦娥三号近方向为西偏南2.44°。1xoz平面直线为y轴。1v=(F/m)siny-m/r2+w=-((F/m)cosy+2vw)/

ωψ为制动发动机推力方向角，其定义为：与当地水平方向w0

r0=v0=

r+r+ prp

r2=R+ J= i0°<y<it1=(1-exp(v1-v0/ve)vem0)/位秒)：500<t1<700。1DACAnum_num_KQr55表 DACA算法结最优燃料消耗着陆径向速度-着陆极角角速度5.12E-图2m1=1289.24kg。1mv2+mgh-(m-mt)gh=mvv 21 12 12 快速调整阶段t2=10.32s，m2=1264.47kg。轨迹较陡，近似为直线。我们先用MATLAB软件读出附图3中各位置的高度，再用VisualBasic程序来求各点减速过程。而单位时间的秒耗量m又可分为竖直方向的秒耗量m^和水平方方向的秒耗mvt mvt 1mvt

C e31 e3132 e t31=

J=m(t31+t33)+m^(t31+t32+t33J=m3-gg(mg+mv) m31

g eve

gg(m31g+mv) 1gmvegmvet2+mt31l2t+ 求解该线性规划，可得出该平移过程燃料消耗最小时m33=1277kg，t31t32=4.37st33=2.06s Dh1=2v3=h-h=Dh

-mt)g=(m-mt 3 3t35t35 35adt=,,和VisualBasic程序来设计计算，计算公式也相同。得出嫦娥三号向西北移动l=0.41m,其

得t5=2.22s图3表31234567向下向左向左向左向右向下向右距离向上向上向上向下向下向下向上表 研究各点与预定着陆点偏离距离1234567在问题二中我们依据燃料最省原则确定了嫦娥3号的着陆轨道，但总耗时为段的参数为:第一阶段t1=600s,m1=1200.08kg;第二阶段t2=12.06s，m2=1100.00kg;第3阶段m33=1096.07kg，t31=1.00s,t32=3.24s,t33=1.00s;第4阶段t44=4.28s,t45灵敏度分析2：初始速度不变阶段一参数不变阶段二推力变为6500N是，t2m2=1189kg;阶段三m33=1156kg，t31=1.45s,t32=3.89s,t33第4阶段t44=4.68s,t45=6.34s,m4=1030kg;第3阶段t5=2.22sJ=1370kg；消耗总时间t=664.7s。同理比较得误差较大，灵敏度比较高。离散的优化问题，从而限制了其在连续优化问题中的应用。【4】针对PID控制器参数【1Proceedingsofthe25thChineseControlConference，7-11August，2006.【2】赵吉松，谷良贤,高原，月球软着陆轨道的时间逼近法快速优化设计，宇航学报，第29卷第5期：1531-1535，2008.28卷第2期：409-413，2007.天大学学报，第38卷第11期：1537-1541，2012.【6】张圣勤，MATLAB7.0实用教程，北京：机械工业出版社，2009781731 2424244 895 1144.49holdon;formatlong;f=@(x)1/1.62*(1.62-2940/(1500-x.^1));[y,n]=quad(f,-314.8,0)formatlong;f=@(x)1/1.62*(2*pi)*exp(-x.^2/2);[