2011’中国电机工程学会杯全国大学生电工数学建模竞赛总结

第一篇：2011’中国电机工程学会杯全国大学生电工数学建模竞赛总结

2011’中国电机工程学会杯

全国大学生电工数学建模竞赛总结

“2011’中国电机工程学会杯全国大学生电工数学建模竞赛” 在全国众多高校、国内数学建模专家以及参赛学生们的大力支持、参与下，现已落下帷幕。在此，竞赛组委会向所有关心、支持竞赛的兄弟院校、专家、教师及参赛学生们表示衷心的感谢！

本届竞赛全国共有178所高校参赛，收到有效答卷2376份。经过论文评阅、复查、成绩公示、复核异议等环节，竞赛成绩现已确定，详细信息见附件1—3（只公布一、二、三等奖名单）。

组委会将对获得一、二、三等奖的学生颁发获奖证书，证书按照附件1中的地址邮寄。

证书制作环节咨询方式：

联系电话：0432－64806324

电子邮件： zhanglei@mail.nedu.edu.cn

联 系 人： 张老师

全国大学生电工数学建模竞赛组委会2012年2月22日

附件1：2011全国大学生电工数学建模竞赛一、二、三等奖成绩汇总表

附件2：2011全国大学生电工数学建模竞赛参赛学校及获奖情况统计表

附件3：2011全国大学生电工数学建模竞赛A、B题情况统计表

第二篇：全国大学生电工数学建模竞赛章程

一．总则

全国大学生电工数学建模竞赛（以下简称竞赛）是中国电机工程学会电工数学专委会主办的面向全国大学生的科技活动，目的是提高学生的综合素质、增强创新意识、培养学生应用数学知识解决实际工程问题的能力，激发学生学习数学的积极性，同时也将推动高校的教学改革与教育创新的进程。

二．竞赛内容

竞赛题目一般来源于电工、近代数学及经济管理等方面，经过适当的简化、加工的实际问题，主要包括：

1．信息处理问题；

2．控制理论及应用问题；

3．运筹与决策问题；

4．电路与电磁场理论相关问题。

参赛学生应学过普通高校的工科数学课程及相关专业的专业基础知识，不要求参赛者预先掌握深入的专门知识。竞赛题目比较灵活，能够使参赛学生充分发挥其创造能力。参赛者要根据题目要求，完成一篇包括模型的假设、建立和求解、算法的设计和计算机实现、结果的分析和检验、模型的改进等方面的论文（即答卷）。答卷的评定以假设的合理性、建模的创造性、结果的可行性和文字ss的清晰度为主要标准。

三．竞赛形式、规则和纪律

1．全国统一竞赛题目，采取网上竞赛方式，以学校为单位进行。

2．竞赛一般在每年11月末的三天内举行。

3．本科生以队为单位参赛，每队3人，专业不限。每队可设一名指导教师（或教师组），对参赛学生进行赛前的辅导及赛前准备工作。竞赛期间指导教师要回避参赛队员，禁止进行指导或参与讨论。

4．竞赛期间参赛学生可以使用各种图书资料、计算机和软件以及在网上浏览，不能与队外任何人讨论。

5．竞赛题目将按照规定时间准时在指定的网站公布，参赛队员在规定的时间内完成答卷，并准时在网上交卷。

6．各参赛院校应责成有关职能部门负责竞赛的组织和纪律监督工作，保证竞赛的规范性和公正性。

四．组织形式

竞赛由电工数学专委会主办，负责每年报名、拟定赛题、组织答卷的评定、评奖、印制获奖证书、举办颁奖仪式等。

五．奖励办法

1．电工数学专委会聘请专家组成评阅委员会，评选一、二、三等奖，获奖比例一般不超过参赛队数的二分之一，其余凡完成答卷的参赛队获得成功参赛奖。

2．一、二、三等奖均颁发获奖证书。

第三篇：2014全国大学生数学建模竞赛

嫦娥三号软着陆轨道设计与控制策略

摘要

随着月球探测任务的发展,未来月球探测考察目标将主要是 复杂地形特性的高科学价值区域。为了能够安全地在这些遍布岩石、的区域内完成高精度软着陆，这就要求导航和控制系统具有较强的自主性和实时性。本文针对最终着陆段安全、精确的需求，对月球软着陆导航与控制方法进行较深入研究，主要内容包括：

首先，提出一种基于单帧图像信息的障碍检测方法。该方法根据着陆区内障碍成像的特点，通过匹配相应的阴影区与光照区完成对岩石、弹坑的检测，利用图像灰度方差对粗糙区域进行提取：在检测出故障信息的基础上，选取安全着陆点以保证软着陆任务的成功。

其次，给出一种基于矢量观测信息的自主光学导航方法。该方法利用光学相机和激光测距仪测量值构建着陆点相对着陆器的矢量信息，结合着陆器的姿态信息确定着陆器的位置。为了消除测量噪声带来的干扰，利用扩展Kalman滤波理论设计了导航滤波器。

再次，提出一种李雅普诺夫函数障碍规避制导方法。该方法通过对状态函数、危险地形势函数的设计，以满足平移过程中减低障碍威胁与精确定点着陆器，设计PWPF（调频调宽）调节器实现定推理等效变推力控制效果。

最后，针对采用变推力主发动机的月球着陆器，提出一种垂直软着陆控制方法。该方法采用标称控制与闭环控制相结合的方式，规划标称轨迹以保证着陆器到达着陆点时其下降速度、加速度亦为零，设计闭环控制器产生附加控制量消除初始偏差、着陆器质量变化的干扰，以保证着陆器沿标称轨迹到达着陆点。

本文分别对所提出的最终着陆段导航与控制方法进行数学仿真以验证个方法的可行性。仿真结果表明，本文多给出导航方法能够达到较高的性能指标，满足在危险区域实现高精度软着陆的需要。

关键词： 月球软着陆；自主导航与控制；障碍检测；规避制导；适量测量

一、问题重述

嫦娥三号于2013年12月2日1时30分成功发射，12月6日抵达月球轨道。根据计划，嫦娥三号将在北京时间12月14号在月球表面实施软着陆。嫦娥三号如何实现软着陆以及能否成功成为外界关注焦点。嫦娥三号在着陆准备轨道上的运行质量为2.4t，其安装在下部的主减速发动机是目前中国航天器上最大推力的发动机，能够产生1500N到7500N的可调节推力，进而对嫦娥三号实现精准控制。其比冲（即单位质量的推进剂产生的推力）为2940m/s，可以满足调整速度的控制要求。在四周安装有姿态调整发动机，在给定主减速发动机的推力方向后，能够自动通过多个发动机的脉冲组合实现各种姿态的调整控制。嫦娥三号的预定着陆点为19.51W，44.12N，海拔为-2641m。嫦娥三号将在近月点15公里处以抛物线下降，相对速度从每秒1.7公里逐渐降为零。整个过程大概需要十几分钟的时间。在距月面100米处时，嫦娥三号要进行短暂的悬停，扫描月面地形，避开障碍物，寻找着陆点。之后，嫦娥三号在反推火箭的作用下继续慢慢下降，直到离月面4米高时再度悬停。此时，关掉反冲发动机，探测器自由下落。

嫦娥三号在高速飞行的情况下，要保证准确地在月球预定区域内实现软着陆，关键问题是着陆轨道与控制策略的设计。其着陆轨道设计的基本要求：着陆准备轨道为近月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软着陆过程共分为6个阶段，分别为着陆准备轨道、主减速段、快速调整段、粗避障段、精避障段、缓速下降阶段，要求满足每个阶段在关键点所处的状态；尽量减少软着陆过程的燃料消耗。

根据上述的基本要求，请你们建立数学模型解决下面的问题：

（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。

（3）对于你们设计的着陆轨道和控制策略做相应的误差分析和敏感性分析。

二、问题分析

对于问题一：

嫦娥三号从15公里左右的高度下降到月球表面，在这一过程中不考虑月球表面太阳风的影响，忽略月球的自转速度引起的科氏力的影响，由于下降时间比较短也不考虑太阳、地球对嫦娥三号的摄动影响，嫦娥三号水平速度要从1.692km/s降为0m/s由于3000m处时嫦娥三号已经基本位于着陆点上方，所以此时假设在3000m处的速度只存在竖直向下的速度而不存在水平分速度，因为降落减速时间比较短只有垂直于月面的方向运动才能实现，所以在确定着陆点位置和着陆轨迹时应当考虑燃料最优情况下推力最大，方向自由的方法即取F7500N建立主减速段动力学模型。

三、符号说明

四、模型假设

对于问题一：

忽略月球的自传和太阳、地球对嫦娥三号卫星的引力摄动 月球近似为一个质量均匀的标准球体 将嫦娥三号是为一个质点

主减速忽略动作调整所产生的燃料消耗段不考虑太阳风的影响

五、模型建立与求解

5.1问题一的建模与求解 解法一： 假设嫦娥三号在t时刻在远月点开始缓慢下降，在n时刻到达近月点，整个过程遵循开普勒第三定律，即

v00

在t时刻有：v12R1 R0R0R1r0 R0r1r2 其中v1：远月点速度

v2：近月点速度

R0：远月点月心距

R1：近月点月心距（已知月球的半径为1738千米）

R017381001838km

R11738151753km 在t1时刻处v2 k2R1 R0R0R1R00.512k0.488 R0R1利用能量平衡式求得近地点速度为

20.51249012（）1.692km/s（沿切线方向）v2，比当地的环境速度17531.672km/s大vk0.0196km/s，径向速度vk0。

1同理解得v11.6139km/s（沿切线方向）

vri0

解得主减速段动力学模型的建立：

根据题意，在横向飞行的水平距离远远小于月球半径的平均值，所以可以将整个减速段过程简化为水平和竖直方向运动方程，根据牛顿第二定律、速度计算公式有：

axTx maytTymTxta

1.692km/s m0Qdt0Tyadt57m/s t0mQdt0tT22xTy27500N

v22atS

运用matlab编程解得S451810.4m； 其中 ax：水平方向加速度

ay：竖直方面加速度

a：月球表面重力加速度a Tx：推力的水平方向分力

Ty：推力的竖直方向分力

t:主减速段时间

S:嫦娥三号主减速段水平位移

Q:嫦娥三号发动机燃料秒消耗率

根据已知资料得到嫦娥三号着陆过程中纬度改变，经度基本不变，月球赤纬和地球纬度一样也分为南北各90个分度,又因为月球极区半径为1735.843km，所以每一个纬度的竖直高度差为19.2871

4g 6千米。即近月点位置坐标为19.0464W,28.9989N海拔15km,远月点位置坐标为160.9536E,28.9989S海拔100km。

解法2:轨迹方程法。

众所周知,太阳系中的八大行星都在按照各自的椭圆轨道绕太阳进行公转,太阳位于椭圆的一个焦点上,行星的运动遵循开普勒三定律,笔者发现,在各类物理竞赛中,常会涉及到天体运动速度的计算,本文拟从能量和行星运动的轨迹方程两个不同的角度来探索行星在近日点和远日点的速度。

该解法的指导思想是对椭圆的轨迹方程求导,并结合一般曲线的曲率半径通式求出近日点和远日点的曲率半径表达式,然后利用万有引力提供向心力列方程求解。如图1所示,椭圆的轨迹方程为

x2y221 5 2ba将5式变形为

a2x2b2y2a2b2 6

根据隐函数的求导法则将6式对x求导有

2a2x2b2yy0 7 即

a2xy2 8

by将7式再次对x求导得

2a22b2(yyyy)0 9 将8、9两式联立得

a2b2y2a4x2 10 y-43by根据曲率半径公式有 r(1y)11 y122 将8、10、11式联立并将A点坐标A(0，a)代入可得A点的曲率半径为

b2RA 12

a根据椭圆的对称性,远日点B的曲率半径为

b2RBRA 13

a 由于在A、B两点行星运行速度方向与万有引力方向垂直,万有引力只改变速度方向,并不改变速度大小,故分别根据万有引力提供向心力得

GMmmvA 14 (ac)2RAGMmmvB 15 2(ac)RB将13至15式联立可得 22vAbGMbGM，vB acaaca

5.2问题二的建模与求解 模型一：动力学模型

典型的月球软着陆任务中,探测器一般首先发射到100km的环月停泊轨道,然后根据所选定的着陆位置,在合适的时间给着陆器一个有限脉冲,使得着陆器转入近月点(在着落位置附近)为15km,远月点为100km的月球椭圆轨道,这一阶段称为霍曼转移段。当着陆器运行到近月点时,制动发动机开始工作,其主要任务是抵消着陆器的初始动能和势能,使着陆器接触地面时,相对月面速度为零,即实现所谓的软着陆,这一阶段称为动力下降段。着陆器的大部分燃料都是消耗在此阶段,所以月球软着陆轨迹优化主要是针对动力下降段这一阶段。由于月球表面附近没有大气,所以在飞行器的动力学模型中没有大气阻力项。而且从15km左右的轨道高度软着陆到月球表面的时间比较短,一般在几百秒的范围内,所以诸如月球引力非球项、日月引力摄动等影响因素均可忽略不计,所以这一过程可以在二体模型下描述。其示意图如图1所示,其中o为月球质心,x轴方向为由月心指向着陆器的初始位置,y轴方向为初始位置着陆器速度方向。

图 1 月球软着陆极坐标系

其动力学方程如下: rv 

v(F/m)sin/rr

22 ((F/m)cos2v)/r

mF/ISP

在上式中r为着陆器与月心距离,v为着陆器径向速度,为着陆器极角,为着陆器极角角速度,为月球引力常数,F着陆器制动发动机推力,m为着陆器质量,为制动发动机推力方向角,其定义为F与当地水平方向夹角,ISP为制动发动机比冲。根据动力下降段的起点位置可以确定动力学方程初始条件,由于起点处于霍曼转移轨道的近地点,故其初始条件为: r0rp

00

v00 01rprp(2ra)rarp其中rp和ra分别为霍曼转移段的近地点半径和远地点半径。

终端条件为实现软着陆, 即

rfR

vf0

f0

其中R为月球半径,终端条件中对终端极角f及终端时间tf无约束。

优化变量为制动发动机推力方向角(t)。

优化的性能指标为在满足上述初始条件和终端条件的前提下, 使着陆过程中燃料消耗最少,即

Jm(t)dt

t0f设计主减速段制导控制律 2动力下降段燃料最优精确着陆问题描述 2.1 燃料最优精确着陆问题

着陆器运动方程：考虑采用变推力发动机情况，有

rv

.vga

(1)

aTmmaT..其中r[rhrxry]T，v[vhvxvy]T分别表示着陆器相对期望着陆点的位置和速度矢量；T为推力器提供的推力矢量，幅值为 T，对应控制加速度矢量 a；g为火星的重力加速度矢量，此处认为是常值；m为着陆器质量，对应推力器质量排除系数。指标函数：考虑燃料消耗

min(m0mf)min0fTdt

(2)边界条件：即初始条件和终端条件

r(0)r0,v(0)v0,m(0)m0,r(tf)v(tf)[000]

(3)控制约束：考虑发动机一旦启动不能关闭，存在最大和最小推力约束

0T1TT

2(4)状态约束：为避免在着陆前撞击到火星地表，需确保整个下降段位于火星地平面以上，即

rh0

（5）进一步地，若着陆区域附近表面崎岖不平，仅仅确保地表约束不能满足需求时，可以考虑下降倾角约束，即将着陆器下降轨线约束到以着陆点为顶点的圆锥体内

2.2 等效后燃料最优精确着陆问题 定义等效变换变量

Ttrx2ry2rhtanalt

（6）

uaT

m

（7）

Tmzlnm等效着陆器运动方程： .r0I3..

yv00.00z其中p[uT0r0vI030z07*70ug0AcyBc(pg4)

（8）],g4[gTT0]T

t指标函数：

min0f(t)dt

（9）

边界条件：同式(3)。

控制约束：由文献[10]可知，控制约束(4)可等效表示为

u1T1ez0[1(zz0)(zz0)2]T2ez0[1(zz0)]

（10）（11）

2状态约束：地表约束同式(5)，倾角约束(6)可等效表示为

T

Sycy0

（12）

其中

0100000S

0010000ctanalt

T000000

3.燃料最优精确着陆问题的离散化及变换 3.1 等效燃料最优精确着陆问题的离散化

首先将整个飞行时间均分成 n 段（对应 n +1 个点），每段步长为t，离散化后的着陆器运动方程为：yk1AykB(pkg4)

其中AR77,BR74分别为离散系统的系统矩阵和输入矩阵

12AetAcI3tActAc

2tt112BetsAcBcdsesAcdsBctBctBct2Bc

0026其中I3为三阶单位阵。

有系统性质可知，整个控制时域内系统状态满足 y3Ay2Bp2g4A3y0A2Bp0g4ABp1g4Bp2g4ynAyn1Bpn1g4Any0An1Bp0g4ABpn2g4Bpn1g4y1Ay0Bp0g4y2Ay1Bp1g4A2y0ABp0g4Bpn2g4Bp1g4

为表达方便，令

y0p00A0yp1111A ,pp2，2A2 Yy2nyn7n11pn4n11nA7n1700B1AB223ABn1An则(15)可等价于

000B012ABBB0002 ABB003AABB0n1AABBA2BABBn7n14n1000000Yy0pg4

分别定义如下常值矩阵：

最终可得离散化后的燃料最优化问题如下： 指标函数：式(9)可表示为

边界条件：式(3)可表示为

控制约束：式(10)和式(11)分别可表示为

状态约束：式(5)和式(12)分别可表示为

含有 p个线性约束和 q个二阶锥约束的最优化问题的标准形式为 指标函数

min(Tx)满足约束

DTxf0AxcibdinTiTi

（k=1,，n）

n*pp其中xR为待优化向量，R，线性约束参数DR,fR，二阶锥约束参数维数n(Ai,bi,ci,di)由相应约束确定

则式(17)～式(23)可最终转换为如下最优化问题： 指标函数：min(vpp)满足：

初值约束:MxΨ0pMx(Ψ0y0)A0g4r0末值约束:MxΨ0pMx(Ψ0y0)A0g4控制约束:Murkpvrkp 控制上限:(vzΨkTTTTv0T0

0

T1vr)p1vTz(Φky0Akg4)z0,z0 z0kT2e 控制下限：

4数值仿真结果与分析本节以某火星着陆器为例，计算了典型初始条件下满足各种约束的燃料最优精确着陆轨迹。其中探测器各参数分别取为：m02000kg,g[3.711400]ms2,c2kms,T11.3kN,T213kN.。着陆器初始位置矢量r0= [1500,-600, 800] m，初始速度矢量v0= [-30, 10, 40]m/s，倾角alt=86°。二阶锥优化问题可以通过大量免费的优化工具求解，如 CSDP、DSDP、OpenOpt、SeDuMi、SDPA、SDPLR等。本文选用 SDPT3 进行计算，通过执行线性搜索确定燃料最优下降时间tf为 43s，图 1 给出了相应的最优着陆轨迹、下降速度、加速度、控制推力、推力幅值以及探测器质量变化曲线。

由优化结果可以看出，探测器在给定时间飞行并软着陆到指定位置，且在整个下降过程始终与火星地表保持一定的安全距离，验证了下降倾角约束的有效性。其推力幅值曲线呈现“最大-最小-最大”的最优控制形式，不过为了保持发动机始终处于点火状态，在中间段对应最小推力约束，这与文献中的分析结论一致。此外，通过利用如 TOMLAB 等商业最优控制软件进行复核计算，也验证了此计算结果的燃料最优性能。

*

图 1 给定初始条件下火星着陆器动力下降段燃料最优计算结果

需要注意到，此燃料最优轨迹的获取对着陆器的实时在线计算性能提出了较高的要求，经测试，无论使用何种优化工具，计算给定飞行任务时间的最优轨迹均需数秒，而全局最优则需要数十秒甚至更长，这在实际任务中是不允许的。因此，可行的方案是通过在地面计算大量的燃料最优轨迹，并寻找规律，选取关键路径点状态存储到着陆器计算机中，通过在线查表或者在利用对计算量要求较小的反馈制导律完成安全着陆任务。

因此，为了研究探测器燃料最优轨迹特性，选取相同的探测器参数，暂不考虑推力器最小幅值约束和倾斜角约束（但考虑地表约束），固定初始高度为 1500m，初始位置水平方向从-8000m 到 8000m 内取值，分别选取各种不同的初始速度，可得燃料最优精确着陆轨迹簇如图 2 所示。

图 2 各种不同初始速度对应的火星着陆器动力下降段燃料最优轨迹簇

1)对任意探测器初始位置，特定初始速度对应的燃料最优着陆轨迹在末端必然收敛到一个固定的近似圆锥体内。

2)取决于探测器初始位置和速度的关系，燃料最优轨迹有两种形式：S 型和 C 型，其中 S 型主要对应于期望着陆点位置水平距离较大情况。3)当探测器初始水平速度为零时，圆锥体轴线垂直于火星地表，所有最优轨线关于该轴线中心对称。4)初始速度的大小也直接影响到任务的可靠性，因此需要在超声速进入段和降落伞减速段将着陆器速度下降到合理范围内。

上述结论对上注探测器关键点的选取有着较强的指导意义，比如基于最优轨线的斜率对路径点合并、基于最优轨线簇的对称性对上注轨线进行等效延伸、或者尝试仅将 S 型和 C 型的转折点作为路径点等，这样可以大大降低探测器自主存储与计算需求，进而有效提升任务的可靠性。重力转弯软着陆过程

对于最终着陆点，假设探测器的下降轨迹在一平面内，且月球引力场为垂直于月面XY的均匀引力场，引力加速度g沿-Z，如图1所示，制动推力方向沿探测器的本体轴z。重力转弯软着陆过程中探测器质心动力学方程可表示为

上式中各变量的物理意义如图1中所示，其中m>0为探测器质量；k>0为制动发动机比冲；u表示制动发动机的秒耗量

可通过一定的机构加以调节，故作为软着陆问题的控制变量。假定制动发动机的最大推力与初始质量比大于月面引力加速度，并且制动推进系统能够在一定的初始条件下将探测器停止月面上。

重力转弯过程中，探测器的高度、速度和姿态角度可由雷达高度表、多普勒雷达及惯性仪表测得。令软着陆初始条件探测器到达月面时速度减小到给定的值，故终端条件自由。软着陆燃耗最优问题的描述 对于最终着陆段，可假设

为一小角度。由此可将系统方程（1）化简为

要设计制导律实现软着陆，就是使

着陆时间

对于月球软着陆的燃耗最优控制问题，其性能指标可表示为

对于系统（2）的软着陆过程，燃耗最优问题等价于着陆时间最优问题，性能指标为

在月球重力转弯软着陆过程中，如果存在一个推力控制程序将探测器从初始条件转移到终端条件，并使性能指标（3）或（4）式最大，则称这个推力程序为软着陆燃耗最优或时间最优制导律。根据pontryagin极大值原理，系统的哈密顿函数及其对u的偏导数为

使哈密顿函数（5）式达到极大地控制输入u就是最优控制，科表示为。

如果存在一个有限区间

则最优控制u（t）取值不能由哈密顿函数确定。此时如果最优解存在，则称为奇异解，（8）式称为奇异条件。

最优制导问题的性质：1）对于自治系统（2）的时间最优控制问题，沿最优轨迹其哈密顿函数满足

将其对时间求导并将（2c）和（6c）式代入，得

另外，由于自由，根据横截条件有3）根据（6a）式。又由（9）式可得T（t）=0，4）根据极大值原理，系统的状态变量和共轭变量都是时间的连续可微函数，将切换函数对时间求导，利用（2），（6）式和性质2）得 软着陆最优控制中奇异条件的分析

对于月球重力转弯软着陆问题，最优制导律具有两个很好的性质。

定理一。月球重力转弯软着陆系统（2）的燃耗最优制导或时间最优制导问题不存在奇异条件。证明。用反证法，假设存在奇异条件，则在某个闭区间设，并由（5）式得

。根据反正假将（10）式两边对时间求导，并将（2）和（6）式代入化简得性质2），并考虑到或者情形1.得

下面证明这两种情形均与反证假设矛盾。根据式

及性质2）可知，由性质3）必有

根据

是时间t的斜率非零的线性函数，m和情形2.1）若定，根据横截条件有在区间内为常数。这与反证假设矛盾。

。下面再分三种情况进行分析。

又因为

不与此时由（6b）式有反证假设矛盾。2）若盾。3），与反证假设矛又

因

为

因此有成立，这与

此时（10）式在上根据定理一，重力转弯软着陆的最优制导律是一种开关（Bang-Bang）控制，只须控制发动机开关，不需要调节推力的大小。

定理2.对于月球重力转弯软着陆过程，其开关控制器的最优推力程序（7）最多进行一次切换。

证明。只要证明最多只在一个时间点成立即可。软着陆系统（2）在最优推力控制程序（7）的作用下，按最后轨迹降落。由性质3）知，为常数。根据性质4），若严格单调，因而在上至多有一个零点，即至多进行一次切换；若，则上为常数。由定理1，5 软着陆最优开关制导律

不可能在任何区间上成立，故必有既没有切换点。

对于最优推力控制程序（7），其切换函数中含有共轭变量，它是一个关于状态变量的稳式表达式。为实现实时制导，需求出关于状态变量的切换函数来。

根据定理一和定理二，重力转弯软着陆最优控制程序没有奇异值状态，并且在着陆过程中最多切换一次，其工作方式有4种：1）全开；2）全关；3）先开有关；4）先关后开。对于方式1）软着陆起始点即是开机点；方式2），3）不能实现软着陆；最后一种是通常情况下的最优着陆方式，即探测器先做无制动下降，然后打开发动机软着陆到月面。设开机时刻为到发动机工作时间为

式，在区间

内积分，并考虑

将（11）式中的对数按泰勒展开，忽略

并令

消掉T得到切换函数为

由切换函数（12）式可以看出，速度、位置的误差和制动发动机推动的将直接影响着陆的效果。一种方法是将终端高度从到达月面时实现软着陆设置为离月面还有几米时实现软着陆。另一种方法是考虑制动过程由一个主发动机和一组小推力发动机共同完成，通过调整开启的小发动机的数量，来实现变推力降落。具体地，令切换函数为

式中各符号的含义如图2所示

关机点可取为2m,可取为20m，可取为1m/s。为实现着陆的最优性，减速度

取为

其中T如（12）式中所示，m0为探测器的初始质量。

图三为最优着陆过程与其改进方法按图2降落的次优着陆过程的对比图。由此图中可看出，改进方法提高了着陆的安全性，当探测器的初始质量mo=350kg，发动机着陆过程多消耗燃料2.2kg。

时，改进方法比最优

（a）

（b）

问题三 协方差分析方法的基本原理 对于如下非线性函数关系

yfx1,x2xn（1）

可以使用一阶泰勒级数展开对其进行线性化，有

yyfffx1xnx1xn（2）x1xn其中，x1xn为x1xn的高阶项。从而得到线性化方程

yfxi（3）i1xin或表示为

YPX(4)

这里 P 是偏导数矩阵: Pif（5）xi若自变量x1xn是随机变量，则线性化方程的函数y的协方差矩阵为：

EYYTEPXXTPTPEXXTPT(6)即 CyPCXPT(7)式中Cx是自变量的协方差矩阵；Cy是函数Y的协方差矩阵。

协方差矩阵中对角线元素是方差，非对角线元素为协方差。显然，只要求出传递矩阵 P ,便可确定源误差与欲求量误差之间的关系。若给定各种源误差，如发动机安装误差、敏感器测量误差或发动机推力和点火时间等误差时，便可以分析其对目标轨道误差的影响以及对控制系统精度的影响，进一步对各系统及元部件提出适当的精度要求。计算向月飞行轨道误差的协方差迭代方程

考虑到轨道参数的误差之相对于轨道参数的标称值是小量，因此可以将轨道运动方程进行线性化，从而得到能够反映轨道参数偏差量的传播关系的误差方程。在应用双二体模型且在地球影响球范围内时，对轨道运动产生摄动影响的各项，如月球引力摄动、太阳引力摄动、大气阻力摄动和太阳光压摄动等对误差方程的影响很小，因此在误差方程中将它们忽略掉。反映轨道位置和速度误差的线性化方程如下：

vrg（8）vrrTur，其中u为地球引力常数。式中 gr3rrrx2ry2rz2（9）

写成状态方程形式：

0Irr(10)vG0vg式中 GT

r0Ir令FG0,Xv（11）

则式（9）变为

FX(12)X下面推导矩阵 F 的表达式：

guGTT3rrrruurrT33Trrrruuuur3333I3rrrrryzrxr（13）

式中 r x，r y 和 r z 是探测器在地心惯性坐标系里的轨道位置坐标。则Gu3T(Irr)(14)332rrrx2rxryrxT2rrryrxryrzryrxryrrzrxrzryzrxrzryrz（15）2rz

将式（15）、（14）代入（10），得： 0002-urx(132)Fr3r3urxryr5v3urxrzr5

积分式（11），得到： 0003urxryr520003urxrzr53urzryr5210000ry-u(13)32rr3urzryr5-urz(13)0r3r200100100（16）

0000

XteFtX0

（17）式中

(Ft)2(Ft)3(Ft)4(Ft)neIFt2!3!4!n!

（18）iNtFi.()i!i0Ft取前 6 阶截断，即：

eFttiFi!

（19）i06i

得到计算误差方程的迭代方程：

XtiteFtXti

（20）

eFt相当于式（4）中的 P 阵，由于误差方程是时变方程，因此每一步迭代都需要重新计算 P 阵，计算 P 阵需要利用标称轨道参数数据。

进一步根据式（7），得到协方差矩阵的迭代方程：

T

Ci1PCPiii

（21）向月飞行轨道误差的协方差分析

引起轨道误差的误差源主要是导航误差，包括位 置 误 差 和 速 度 误 差。其 中 ： 位 置 误 差 ：rrx,ry,rz,rx,ry,rz分别为在地心惯性坐标系中 X 轴、Y 轴、Z 轴的分量。速度误差：vvx,vy,vz,vx,vy,vz分别是在地心惯性坐标系 X 轴、Y 轴、Z 轴的分量。向月飞行轨道的初始轨道位置和速度误差由运载火箭的发射入轨精度决定，若探测器在飞行途中进行轨道修正，则经过轨道修正以后的轨道位置误差将由导航误差决定，速度误差将由姿态误差和制导误差决定。

上述误差决定了轨道误差协方差分析的计算初始条件，表 1 给出了在不进行中途轨道修正情况下，在地心惯性坐标系里，初始轨道位置误差和初始速度误差对轨道终点的位置和速度误差的影响。图 1 和图 2 给出了在算例三中探测器从近地轨道入轨点开始至进入月球轨道为止轨道位置的相应的轨道位置和速度总误差（3σ）的时间历程。

表 1 初始轨道位置和速度误差

对轨道终点误差的影响

图 1 轨道位置总误差时间历程（3σ）

图 2 速度总误差时间历程（3σ）基于敏感系数矩阵的制导误差分析

在月球软着陆主制动段，影响制导精度的误差源主要有偏离标准飞行轨迹的初始条件误差和导航与控制传感器误差。初始条件误差由主制动段以前的任务决定，传感器误差则由导航系统和传感器本身决定。此外，影响制导精度的因素还包括月球自转、月球不规则摄动等误差，对它们的研究可单独进行，这里暂不做介绍。2.1 误差模型建立

2.1.1 初始状态误差模型

记着陆器的实际初始状态为Xi，标准初始状态为Xn,则定义初始状态偏差xi为

xiXiXn

(7)对于主制动段这一特定的飞行过程，这些偏差都是确定的；而针对整个月球探测任务，这些偏差就变得具有随机性。在本文中，假定xi 的所有元素均服从零均值高斯分布，相互不独立，其相关性取决于前一阶段任务的特性。2.1.2 传感器误差模型

由于只研究误差对制导律的影响，所以这里假设需要测量的量均可由导航系统直接测得，误差大小

均考虑为典型误差值。由上一目设计的制导律可以看出，需要由导航与控制传感器测量的量主要为着陆器相对于着陆场坐标系的位置、速度和加速度。定义待测量量Q为

QX其估计值记为Q，则传感器误差定义为 YZUVWA

T

qQQ

(8)那么，单个测量量的估计误差模型可用误差向量 q的第j(j =1，2„7)个元素qj 来表示。由参考文献[5]可知，第 j个观测量的总估计误差qj 由以下四部分组成

~~-~qjbsqjnstqtqQtqtQjt

(9)jjbcjnc

j100100~~~~~针对主制动这一特定操作阶段，上述四部分误差具有如下特性：

qjbc—第 j 个观测量的测量误差，恒为常值，其分布服从零均值高斯分布； qjbs—第 j 个观测量的刻度因素误差系数，恒为常值，其分布服从零均值高斯分布； qjnc—第 j 个观测量的随机误差，其为一高斯白噪声；

qjns

—第 j 个观测量的刻度因素随机误差系数，其为一高斯白噪声。

2.2 制导误差分析

由于采用闭环制导，制导控制系统对随机误差具有一定鲁棒性，所以本文将着重对初始偏差和类似于qjbc和qjbs这样的传感器常值误差进行仿真研究，分析它们对制导精度的影响。2.2.1 误差分析系统建立

误差分析系统框图如图 1 所示，下面将对其结构进行分析。~~~~~~

图 1 误差分析系统结构图

图中所示初始状态偏差实际上是加在相应积分器中。

由前面的分析可知，观测量的实际输出值受到初始状态偏差、传感器测量误差以及传感器刻度因素误差的影响，故误差分析系统模拟程序的实际输入应包含以下几部分(以 X通道为例)：

XXxixbc~xbsX

(10)100~~

其中，X为观测量的实际输出值，X 为标准值，xi 为初始状态偏差(只在初始时刻存在)，xbc 为传感器测量偏差，xbs为传感器刻度因素误差系数。由图 1 可以看出，为了更准确地表示传感器误差模型，这里考虑了传感器的动态性能，其传递函数设为一阶惯性环节11Ts，其中，T 为传感器时间常数，因传感器的不同而取不同值。

由误差分析系统结构框图可以看出，其输入量主要包括：标准初始状态向量、初始状态偏差、传感器测量误差、传感器刻度因素误差系数、传感器时间常数、期望终端状态；输出量为加入误差前后的仿真终端状态向量。2.2.2 误差敏感系数矩阵求取

在有形如(7)式误差输入的情况下，首先根据图 1 生成一个模拟整个闭环制导控制系统的数字仿真程序，然后运行该程序，对比程序输出即可得到误差敏感系数矩阵。具体运行过程如下：

第一步：将传感器误差设置为零，初始状态设置为标准值，运行模拟程序。这一步称为标准运行。第二步： 将其中一个传感器误差设置为非零输入或者设置一个非标准初始状态，然后进行一系列运行。

第三步： 将第二步运行的系统输出和标准运行的系统输出进行比较即可确定各误差源的影响。如X 通道标准初始偏差为xi，输入该误差前后，X 通道终端状态分别为X0 和X1，则 X 通道对标准初始偏差xi的敏感性可用(X1X0)/xi来反映。

通过这种方法，可得到一组反映月球软着陆主制动段终端总误差向量pf和两个传感器误差向量~~~qbc、qbs以及初始状态偏差向量pi之间关系的误差敏感系数矩阵。由参考文献[6]可知，其相互关系可表示为

~~pfS1piS2qbcS3qbs（11）

其中，S1、S2和S3分别表示相对于pi、qbc和qbs的误差敏感系数矩阵。

终端误差向量能用这种形式表示的假设条件是动力学的线性化必须在标准轨迹区域内。验证该假设条件的方法有两种： 扩大输入误差仿真法和复合仿真法，这里略去其验证过程。2.2.3 误差分析

假设导航系统采用常规惯性测量单元，表 1 列出了其典型误差值，其中，位置误差能保持在10数量级，速度在10数量级，加速度为 10g 数量级。1-52~~

运用上述方法得到的敏感系数矩阵给出如下：

5.50210-3-4-3.850101.69210-3S1-38.36210-5.86010-4-3-2.57510-2.08010-4-1.05010-31.41810-11.40110-57.30110-5-1.00110-26.41110-53.24010-4-4.40710-2-2.57010-4-1.86210-3-5.58010-11.41010-57.90210-51.31210-55.71010-4-1.15710-38.10010-53.93610-21.73210-2-2.7431017.74610-1-4.02410-2-8.93910-23.21010-34.03010-31.23910-21.83310-2-2-18.742101.41410-1.19610-2-9.90110-3-2-2-2.69010-4.57710-6.81210-1-8.69510-2-5.2031002.11010-14.23510-16.17010-3-3.2811008.20210-2-5.76010-35.63310-1-3.4891022.4431014.401102-9.8331026.86410123.02010-9.85910-1-1.15410-3-40-3.13010-1.00010-1.37910-33.56010-4S2-2-3-5.402101.540101.04510-31.86410-3-34.77010-44.598109.99910-13.408100-7.21010-43.5041005.00010-55.64310-3-1.52710-19.36810-1-6.72110-1-1.30610-1-5.6314100-28.479103.73010-1S30-8.924104.61910-102.03310-5.49410-1-3.53310-1-2.8101001.60010-31.69210-16.75510-18.99610-1-209510-12.47310-21.66410-1-1.0271007.16510-23.344100-1.1121008.61310-17.8521003.246100-1.6181003.54010-14.98210-17.67010-1-1.122100-2.397100-2.38010-1-3.650100-2.5631002.55610-1-4.29110-23.401100-1.88810-1-5.103100-3.23010-13.56610-12.25610-10-1-7.005109.93010A1、A3:12.7592,30.1297j2.1329 A2:11.5522,30.6761j1.8978

由于数值仿真的起始点选为(1，0，-1)，靠近平衡点(1.5，0，-1.05)，仿真实验中混沌系统的基频w0=2.1329，基周期为为T0202.9443S。由前面的数值仿真实验知要使 Chua’s混沌系统保持其类随机性，仿真步长选在(0.0001，0.7)较为合适，用基周期来表达即为129940T015T0 ,15T0内，综观三个连续混沌系统仿真步长的理论计算，我们可以统一选取15000T0这样即可以提高仿真运算速度，又可以使混沌吸引子的形状和类随机性不发生变化，这个选择范围也与通常连续混沌系统数值仿真步长的经验取值相吻合六、模型结果及分析

七、结果分析

八、模型评价与改进方向

九、参考文献

第四篇：全国大学生数学建模竞赛

全国大学生数学建模竞赛

1、数模竞赛的起源与历史

数模竞赛是由美国工业与应用数学学会在1985年发起的一项大学生竞赛活动，目的在于激励学生学习数学的积极性，提高学生建立数学模型和运用计算机技术解决实际问题的综合能力，鼓励广大学生踊跃参加课外科技活动，开拓知识面，培养创精神及合作意识，推动大学数学教学体系、教学内容和方法的改革。我国大学生数学建模竞赛是由教育部高教司和中国工业与数学学会主办、面向全国高等院校的、每年一届的通讯竞赛。其宗旨是：创新意 识、团队精神、重在参与、公平竞争。1992载在中国创办，自从创办以来，得到了教育部高教司和中国工业与应用数学协会的得力支持和关心，呈现出迅速的发展发展势头，就2003年来说，报名阶段须然受到“非典”影响，但是全国30个省（市、自治区）及香港的637所院校就有5406队参赛，在职业技术学院增加更快，参赛高校由2002年的1067所上升到了2003年的1410所。可以说：数学建模已经成为全国高校规模最大课外科技活动。

2、什么是数学建模

数学建模（Mathematical Modelling）是一种数学的思考方法，是“对现实的现象通过心智活动构造出能抓住其重要且有

用的特征的表示，常常是形象化的或符号的表示。”从科学，工程，经济，管理等角度看数学建模就是用数学的语言和方法，通过抽象，简化建立能近似刻画并“解决”实际问题的一种强有力的数学工具。顾名思义，modelling一词在英文中有“塑造艺术”的意思，从而可以理解从不同的侧面，角度去考察问题就会有不尽的数学模型，从而数学建模 的创造又带有一定的艺术的特点。而数学建模最重要的特点是要接受实践的检验，多次修改模型渐趋完善的过程。

3、竞赛的内容：

竞赛题目一般来源于工程技术和管理科学等方面经过适当简化加工的实际问题，不要求参赛者预先掌握深入的专门知识，只需要学过普通高校的数学课程。题目有较大的灵活性供参赛者发挥其创造能力。参赛者应根据题目要求，完成一篇包括模型假设、建立和求解、计算方法的设计和计算机实现、结果的分析和检验、模型的改进等方面的论文（即答卷）。竞赛评奖以假设的合理性、建模的创造性、结果的正确性和文字表述的清晰程度为主要标准。

4、竞赛的步骤

建模是一种十分复杂的创造性劳动，现实世界中的事物形形色色，五花八门，不可能用一些条条框 框规定出各种模型如何具体建立，这里只是大致归纳一下建模的一般步骤和原则：

1）模型准备：首先要了解问题的实际背景，明确题目的要求，收集各种必要的信息．

2）模型假设：为了利用数学方法，通常要对问题做必要的、合理的假设，使问题的主要特征凸现出来，忽略问题的次要方面。

3）模型构成：根据所做的假设以及事物之间的联系，构造各种量之间的关系,把它问题化

4）模型求解：利用已知的数学方法来求解上一步所得到的数学问题，此时往往还要作出进一步的简化或假设。注意要尽量采用简单的数学工具。

5）模型分析：对所得到的解答进行分析，特别要注意当数据变化时所得结果是否稳定。

6）模型检验：分析所得结果的实际意义，与实际情况进行比较，看是否符合实际，如果不够理想，应该修改、补充假设，或重新建模，不断完善。

7）模型应用：所建立的模型必须在实际应用中才能产生效益，在应用中不断改进和完善。

5、模型的分类

按模型的应用领域分类： 生物数学模型、医学数学模型、地质数学模型、数量经济学模型、数学社会学模型

按是否考虑随机因素分类 ：确定性模型、随机性模型按是否考虑模型的变化分类 ：静态模型、动态模型按应用离散方法或连续方法 ：离散模型、连续模型

按建立模型的数学方法分类 ：几何模型、微分方程模型、图

论模型、规划论模型、马氏链模型 按人们对事物发展过程的了解程度分类 ：

白箱模型： 指那些内部规律比较清楚的模型。如力学、热学、电学以及相关的工程技术问题。

灰箱模型： 指那些内部规律尚不十分清楚，在建立和改善模型方面都还不同程度地有许多工作要做的问题。如气象学、生态学经济学等领域的模型。

黑箱模型：

指一些其内部规律还很少为人们所知的现象。如生命科学、社会科学等方面的问题。但由于因素众多、关系复杂，也可简化为灰箱模型来研究。

6、数学建模应用

今天，在国民经济和社会活动的以下诸多方面，数学建模都有着非常具体的应用。

1分析与设计： 例如描述药物浓度在人体内的变化规律以分析药物的疗效；建立跨音速空气流和激波的数学模型，用数值模拟设计新的飞机翼型。预报与决策： 生产过程中产品质量指标的预报、气象预报、人口预报、经济增长预报等等，都要有预报模型。使经济效益最大的价格策略、使费用最少的设备维修方案，是决策模型的例子。3 控制与优化： 电力、化工生产过程的最优控制、零件设计中的参数优化，要以数学模型为前提。建立大系统控制与优化的数

学模型，是迫切需要和十分棘手的课题。规划与管理 生产计划、资源配置、运输网络规划、水库优化调度，以及排队策略、物资管理等，都可以用运筹学模型解决 报名时间：从大赛的通知文稿发出后，就可以报名了，报名截止时间一般在开始比赛的前7到10天。

竞赛时间：每年的9月的第三个星期五上午8时至下一个星期一上午8时。

报名方式：如果有分赛区（每个赛区应至少有6所院校的20个队参加），就联系分赛区报名，没有分赛区，则直接向主委会报名。

大学生以队为单位参赛，每队3人（须属于同一所学校），专业不限。竞赛分本科、专科两组进行，本科生参加本科组竞赛，专科生参加专科组竞赛（也可参加本科组竞赛），研究生不得参加。每队可设一名指导教师（或教师组）。

第五篇：2011高教社杯全国大学生数学建模竞赛题目

2011高教社杯全国大学生数学建模竞赛题目

（请先阅读“全国大学生数学建模竞赛论文格式规范”）

A题

城市表层土壤重金属污染分析

随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。

按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、„„、5类区，不同的区域环境受人类活动影响的程度不同。

现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10 厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。

附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。

现要求你们通过数学建模来完成以下任务：

(1)给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。

(2)通过数据分析，说明重金属污染的主要原因。

(3)分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4)分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？

附件1：A区和全市六区交通网络与平台设置的示意图。

B题

交巡警服务平台的设置与调度

“有困难找警察”，是家喻户晓的一句流行语。警察肩负着刑事执法、治安管理、交通管理、服务群众四大职能。为了更有效地贯彻实施这些职能，需要在市区的一些交通要道和重要部位设置交巡警服务平台。每个交巡警服务平台的职能和警力配备基本相同。由于警务资源是有限的，如何根据城市的实际情况与需求合理地设置交巡警服务平台、分配各平台的管辖范围、调度警务资源是警务部门面临的一个实际课题。

试就某市设置交巡警服务平台的相关情况，建立数学模型分析研究下面的问题：

（1）附件1中的附图1给出了该市中心城区A的交通网络和现有的20个交巡警服务平台的设置情况示意图，相关的数据信息见附件2。请为各交巡警服务平台分配管辖范围，使其在所管辖的范围内出现突发事件时，尽量能在3分钟内有交巡警（警车的时速为60km/h）到达事发地。

对于重大突发事件，需要调度全区20个交巡警服务平台的警力资源，对进出该区的13条交通要道实现快速全封锁。实际中一个平台的警力最多封锁一个路口，请给出该区交巡警服务平台警力合理的调度方案。

根据现有交巡警服务平台的工作量不均衡和有些地方出警时间过长的实际情况，拟在该区内再增加2至5个平台，请确定需要增加平台的具体个数和位置。

（2）针对全市（主城六区A，B，C，D，E，F）的具体情况，按照设置交巡警服务平台的原则和任务，分析研究该市现有交巡警服务平台设置方案（参见附件）的合理性。如果有明显不合理，请给出解决方案。

如果该市地点P（第32个节点）处发生了重大刑事案件，在案发3分钟后接到报警，犯罪嫌疑人已驾车逃跑。为了快速搜捕嫌疑犯，请给出调度全市交巡警服务平台警力资源的最佳围堵方案。

附件2：全市六区交通网络与平台设置的相关数据表（共5个工作表）。

全市路口节点标号 路口的横坐标X

路口的纵坐标Y

路口所属区

域

发案率(次

数)

说明：

15

413 403 383.5 381 339 335 317 334.5 333 282 247 219 225 280 290

359 343 351 377.5 376 383 362 353.5 342 325 301 316 270 292 335

A A A A A A A A A A A A A A A

1.7 2.1 2.2 1.7

A列：是全市交通网络中路口节点

B列：路口节点的横坐标X，是在交D列：路口节点所属的区

米

坐标的长度单位为毫米

C列：路口节点的纵坐标Y，是在交

2.1 2.5

2.4 2.4 2.1 1.6 2.6 2.4 2.2 2.5 2.1

E列：各路口节点的发案率是每个

地图距离和实际距离的比例是1:1

337 17 415 18 432 19 418 20 444 21 251 22 234 23 225 24 212 25 227 26 256 27 250.5 28 243 29 246 30 314 31 315 32 326 33 327 34 328 35 336 36 336 37 331 38 371 39 371 40 388.5 41 411 42 419 43 411 44 394 45 342 46 342 47 325 48 315 49 342 50 345 51 348.5 52 351 53 348 54 370 55 371 56 354 57 363 58 357 59 351 328 A 335 A 371 A 374 A 394 A 277 A 271 A 265 A 290 A 300 A 301 A 306 A 328 A 337 A 367 A 351 A 355 A 350 A 342.5 A 339 A 334 A 335 A 330 A 333 A 330.5 A 327.5 A 344 A 343 A 346 A 342 A 348 A 372 A 374 A 372 A 382 A 380.5 A 377 A 369 A 363 A 353 A 374 A 382.5 A 387 A 382 A 2.6

2.5

1.9

1.8

1.9

1.4

1.4

2.4

1.1

1.6

1.2

0.8

1.3

1.4

2.1

1.6

1.5 案发地P点的标号：32 1.4

1.7

1.4

1.1

0.1

1.2

1.4

1.7

1.4

1.4

1.7

1.1

1.4

1.2

1.6

1.4

1.2

1.1

0.8

0.6

1.4

0.9

0.5

0.8

1.1

0.9

369 61 335 62 381 63 391 64 392 65 395 66 398 67 401 68 405 69 410 70 408 71 415 72 418 73 422 74 418.5 75 405.5 76 405 77 409 78 417 79 420 80 424 81 438 82 438.5 83 434 84 438 85 440 86 447 87 448 88 444.5 89 441 90 440.5 91 445 92 444 93 140 94 145 95 160 96 142.5 97 150 98 186 99 158 121 157 158 159 388 A 395 A 381 A 375 A 366 A 361 A 362 A 359 A 360 A 355 A 350 A 351 A 347 A 354 A 356 A 364.5 A 368 A 370 A 364 A 370 A 372 A 368 A 373 A 376 A 385 A 392 A 392 A 381 A 383 A 385 A 381.5 A 380 A 360 A 130 B 118 B 96 B 71 B 70 B 145 B 73.5 B 68 B 145 B 138.5 B 135 B 0.7

0.6

1.2

1.4

0.8

0.7

0.8

0.8

0.9

1.1

0.9

1.1

0.8

0.9

1.1

0.8

1.1

0.8

0.8

0.8

0.8

1.4

1.1

0.9

1.2

1.4

1.1

0.9

1.4

0.9

0.9

0.8

1.6

1.6

1.6

2.1

1.8

1.6

2.6

2.6

1.1

0.9

0.5

105 137.5 106 144 107 139 108 144.5 109 151 110 151.5 111 150 112 158 113 159 114 164 115 163 116 149 117 143 118 137 119 131 120 130 121 127 122 125 123 129 124 130 125 124 126 136 127 136 128 142 129 148 130 142 131 147 132 128 133 136.5 134 142 135 147 136 154 137 148.5 138 140 139 140 140 137.5 141 138 142 143 143 151 144 153 145 143 146 143 147 143 114 B 113 B 112 B 117 B 115 B 113 B 118 B 111 B 118 B 109 B B 105 B 99.5 B 102 B 103 B 103 B 100 B 102 B 98 B 96 B 90 B 90 B 96 B 90 B 96 B 96 B 91 B 91 B 71 B 76 B 79 B 81 B 86 B 74.5 B 70 B 63 B 63 B 59 B 63 B 69 B 63 B 60 B 57 B 51.5 B 0.7

0.4

0.8

0.2

0.8

0.6

0.9

0.8

1.1

0.8

0.4

0.7

1.2

0.8

0.9

0.5

0.6

0.6

0.8

0.9

0.4

0.7

1.1

0.8

0.8

0.7

0.6

0.7

1.2

0.8

1.1

0.8

0.9

1.1

0.6

0.7

0.8

0.4

1.1

0.8

1.1

0.7

0.6

0.8

108.5 148 160 149 162 150 141 151 143 152 151 153 150 154 164 155 171 156 165.5 157 181 158 176 159 170 160 168 161 166 162 176 163 180 164 183 165 202 166 137.5 167 167 168 376 169 210 170 263 171 284 172 278.5 173 295 174 299 175 362 176 410 177 395 178 277 179 235 180 200 181 167 182 225 183 400 184 414 185 424 186 411 187 420 188 403 189 376 190 380 191 377 65 B 59 B 49 B 40 B 44 B 33 B 124 B 125 B 139 B 131 B 141 B 140 B 145 B 150 B 145 B 149 B 145 B 131 B 462 C 399 C 400 C 390 C 445 C 409 C 425 C 382 C 444 C 443 C 408.5 C 520 C 496 C 465 C 466.5 C 462 C 443 C 447 C 422 C 400 C 396 C 401 C 404 C 406 C 404 C 424 C 1.1

0.6

0.4

0.8

0.5

0.1

0.6

0.7

1.1

1.4

1.6

0.8

0.6

0.8

0.6

0.7

1.1

1.1

2.6

2.2

1.4

2.6

2.2

1.9

2.2

2.6

2.2

2.1

2.2

1.7

2.2

1.9

2.4

2.4

1.2

1.4

1.2

1.4

0.8

1.2

0.9

0.8

0.8

192 374 193 370 194 368 195 374 196 365 197 356 198 358 199 354 200 357 201 359 202 347 203 261 204 270 205 313 206 324 207 333 208 334 209 323 210 312 211 317 212 316 213 315.5 214 316 215 318 216 291.5 217 284 218 281 219 299 220 302 221 305 222 281 223 274 224 273.5 225 267 226 270 227 275 228 276 229 270 230 276 231 288 232 293.5 233 296 234 303 235 298.5 424 C 423 C 427.5 C 431 C 448 C 450 C 459 C 495 C 513 C 528 C 553 C 537.5 C 514 C 511.5 C 511 C 511 C 497 C 497 C 498 C 451 C 448 C 444 C 434 C 412 C 415 C 425 C 421 C 434 C 451 C 457 C 458.5 C 448 C 444 C 446 C 440 C 422 C 419 C 415 C 405 C 403 C 392.5 C 387 C 386 C 378 C 0.8

0.4

0.9

1.2

1.4

1.4

1.2

1.1

1.2

0.4

0.5

0.8

1.4

0.4

0.8

0.7

0.8

0.7

1.1

1.1

0.8

0.7

0.9

1.2

1.4

1.4

1.6

1.4

1.4

1.2

1.1

0.8

0.9

1.1

0.9

1.1

1.2

0.8

1.4

1.4

1.4

1.1

1.4

1.6

236 293 237 296 238 276 239 250 240 247 241 262 242 269 243 276 244 282 245 282 246 273 247 276 248 138.5 249 155 250 163 251 173 252 183 253 238 254 213 255 189 256 210 257 201 258 150 259 135 260 143 261 142 262 140 263 121 264 109 265 138.5 266 167 267 168 268 184 269 194 270 200 271 212 272 220 273 246 274 246 275 252 276 257 277 255.5 278 249 279 247 376 C 372 C 352 C 350 C 384 C 399 C 397 C 402 C 398.5 C 386.5 C 389 C 361 C 378 C 396 C 390 C 364 C 370 C 382 C 412 C 413 C 433 C 434 C 400 C 395 C 407 C 414 C 430 C 432 C 441 C 442 C 442 C 435 C 440 C 442 C 442 C 443 C 443 C 444 C 455 C 458 C 460.5 C 466 C 464 C 469 C 1.2

1.7

1.4

1.2

1.4

1.2

1.3

1.2

1.1

1.2

1.1

0.8

1.2

1.5

1.2

1.2

0.7

0.8

1.1

0.9

0.9

1.2

1.2

1.1

0.8

1.1

0.9

0.5

1.2

1.6

1.4

1.2

0.9

1.4

1.6

1.7

2.1

1.4

1.2

1.5

1.2

1.1

0.8

280 254 281 251.5 282 259 283 261 284 255 285 240 286 241 287 236 288 235 289 232 290 235.5 291 245 292 225 293 225 294 219 295 219 296 228.5 297 213 298 211 299 208.5 300 206 301 206 302 200 303 200 304 200 305 200 306 206 307 194 308 184 309 184 310 193.5 311 193 312 184 313 184 314 192.5 315 192 316 192 317 170 318 168 319 167 320 101 321 91 322 70 323 46 472 C 477 C 478 C 470 C 494 C 495 C 514 C 514 C 496 C 487 C 486.5 C 474 C 457.5 C 451 C 451 C 462 C 472 C 481 C 487 C 496 C 507 C 515 C 514 C 507 C 497 C 484 C 466 C 466 C 463.5 C 475 C 475 C 484 C 484 C 496.5 C 496.5 C 507 C 514 C 516.5 C 507 C 495.5 C 343 D 355 D 377 D 371 D 0.7

1.1

0.8

0.4

1.4

1.4

0.8

0.7

0.7

0.8

0.8

1.2

1.4

1.6

1.4

1.2

1.6

1.4

1.2

0.8

1.2

0.7

1.2

1.3

1.4

1.4

1.4

1.5

0.8

0.7

0.9

0.6

0.8

0.7

0.9

0.8

0.6

1.1

1.4

2.4

1.7

2.5

2.4

324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 56 424 20 442 74 326 76 302 15 240 28 161 34.5 164.5 30 181 27 206 42 242 30 246 31 254 39 254 50 289 72 288 60 246 95 299 81 297 80 287 67 314 21 330 36 360 76 344 97 339 103 337 104 341 97 345 89 345 81 344.5 81 350 89 350 92.5 351.5 88 353 81.5 353 87 359 84 361 76 355 58.5 370 34 306 38 418.5 61 425 57 429 60 433 85 369 D 2.1

D 2.2

D 2.6

D 2.1

D 2.6

D 0.4

D 0.1

D 0.6

D 0.2

D 1.4

D 1.6

D 1.1

D 1.2

D 0.7

D 1.1

D 0.7

D 1.4

D 1.6

D 1.4

D 1.7

D 1.1

D 1.2

D 0.8

D 2.4

D 1.2

D 1.1

D 1.6

D 0.8

D 0.8

D 0.4

D 0.7

D 1.1

D 1.4

D 0.9

D 1.1

D 0.9

D 0.8

D 0.6

D 0.1

D 1.4

D 1.4

D 1.6

D 1.4

D 1.9

368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 107.5 362 131 366.5 170 342 174 340 232.5 264 202 223 241 210 235 197.5 228 173 214 164 278 196 267 168 90 167 123 177.5 143 153 192 264 145 285 133 255 90 198 21 15 60 68 70 84 27 149 62 143 58 176 58 160 72 163 70 176.5 90 178 115 168 115 177.5 123 168 123 164 123 155 143 164 144 168 149 177 128 178 128 188 164 194 156 177 168 177 156 169 167 168 D 1.4

D 1.2

D 1.2

D 1.5

E 2.4

E 1.9

E 2.4

E 2.6

E 2.6

E 2.6

E 2.6

E 2.4

E 1

E 1.1

E 1.9

E 2.6

E 2.4

E 2.4

E 1.7

E 1.1

E 0.8

E 0.2

E 1.6

E 0.9

E 1.4

E 0.6

E 0.7

E 0.7

E 0.8

E 0.6

E 0.8

E 0.7

E 0.6

E 0.7

E 0.9

E 1.2

E 0.9

E 0.9

E 1.5

E 1.7

E 0.8

E 1.1

E 0.8

E 0.8

412 172 413 167 414 160 415 163 416 186 417 269 418 295 419 302 420 316 421 278 422 284 423 257.5 424 239 425 241 426 246 427 246.5 428 240 429 236 430 231 431 232 432 239.5 433 242 434 235 435 237.5 436 246 437 200 438 170 439 59 440 72 441 90 442 74 443 60 444 90 445 151 446 160 447 90 448 129 449 142 450 152 451 155 452 163 453 171 454 171 455 171 167 E 164 E 164 E 153.5 E 168 E 133 E 112 E 112 E 141 E 143 E 173 E 170 E 198 E 198 E 199 E 202 E 202 E 201 E 199 E 206.5 E 207.5 E 206 E 209.5 E 212 E 208 E 194 E 222 E 189 E 189 E 187.5 E 198 E 196 E 211 E 236 E 244 E 222 E 248 E 265 E 255.5 E 258 E 258 E 258 E 252.5 E 247 E 0.9

0.2

0.7

1.2

1.6

1.6

1.1

1.4

1.6

1.7

1.4

1.9

0.4

0.3

0.6

0.4

0.4

1.1

0.1

0.6

0.5

0.2

0.4

0.1

0.4

1.1

1.6

0.8

0.9

0.6

0.7

0.4

1.6

1.4

1.5

0.8

1.7

2.1

1.1

0.6

0.8

1.1

0.4

1.2

456 214 457 244 458 268 459 259 460 188 461 184 462 171 463 171 464 163 465 154 466 151 467 148 468 162 469 177 470 187 471 155 472 159 473 125 474 107 475 382.5 476 373 477 330 478 400 479 441 480 417 481 332 482 321 483 403 484 420 485 455 486 295.5 487 294 488 316 489 308 490 327 491 316 492 314 493 313 494 317 495 318.5 496 320 497 326.5 498 325 499 323 235 E 238 E 237 E 255 E 261 E 253 E 263 E 268 E 268 E 268.5 E 275 E 274 E 277.5 E 281 E 284 E 316 E 292 E 267 E 285 E 267 F 250 F 219 F 247 F 442 F 312 F 246 F 275 F 140 F 269 F 335 F 238 F 244 F 300 F 257.5 F 255 F 236 F 230 F 223 F 215 F 222 F 229 F 227.5 F 220 F 213 F 1.1

1.1

1.1

1.1

1.4

1.2

1.1

0.8

0.9

0.7

0.4

1.5

1.5

0.7

1.4

1.6

1.8

1.8

1.6

2.4

1.9

0.8

2.3

1.7

1.5

1.9

1.7

2.1

2.4

1.9

1.4

1.1

1.5

1.2

0.8

1.4

0.9

0.6

0.2

0.3

0.5

0.7

0.6

0.4

500 329 501 332 502 334 503 346 504 342 505 356 506 358 507 345 508 348 509 357.5 510 359 511 404 512 403 513 379 514 386 515 373 516 363 517 362 518 354 519 348 520 349 521 353 522 371 523 371 524 375 525 388.5 526 405 527 389 528 388 529 353 530 334 531 336 532 352 533 353 534 362.5 535 370 536 388 537 395.5 538 395 539 408.5 540 430 541 450 542 394 543 387 212 F 226 F 210.5 F 209 F 200 F 202 F 195 F 194 F 188 F 188 F 159 F 161 F 202 F 202 F 213 F 213 F 212 F 218 F 216.5 F 215 F 222 F 223 F 224 F 218.5 F 219 F 218 F 213.5 F 224.5 F 233 F 229.5 F 232 F 239 F 247 F 236 F 236 F 236 F 237 F 237.5 F 233 F 227 F 237 F 268 F 254 F 250 F 0.7

0.7

0.6

0.7

0.8

0.7

0.6

0.4

0.4

0.6

1.1

1.2

0.8

0.7

0.8

0.6

0.4

0.8

0.6

0.9

0.7

0.8

0.8

0.6

0.4

1.1

0.8

0.9

0.6

0.8

0.7

1.1

1.2

0.6

0.8

1.1

1.2

1.4

1.1

1.5

1.4

0.1

1.4

0.9

544 383 545 369 546 367.5 547 362 548 350 549 348 550 355 551 367 552 367 553 375 554 376 555 381 556 378 557 380 558 371 559 356.5 560 338 561 372 562 398 563 392 564 382.5 565 396 566 411 567 424 568 435 569 434 570 430 571 441 572 470 573 468 574 455 575 453 576 425 577 462 578 481 579 462 580 440 581 423 582 435 250 F 249.5 F 249 F 249 F 251 F 255 F 265 F 265 F 257.5 F 258 F 260 F 260 F 266 F 270.5 F 284 F 281 F 297 F 307 F 308 F 277 F 276 F 270 F 291 F 297 F 319 F 307 F 295 F 309 F 342 F 432 F 361 F 400 F 433 F 437 F 457 F 447 F 449 F 448 F 507.5 F 1.1

0.8

0.7

0.8

0.6

1.4

1.1

0.8

1.2

1.4

1.1

1.7

1.4

1.2

1.1

1.4

1.2

1.4

1.5

1.1

0.9

1.4

1.2

0.8

0.9

0.7

1.4

1.2

0.2

1.2

0.6

0.6

0.8

1.4

0.6

1.2

1.4

0.4