第一篇：赣南师范学院数学建模训练题

赣南师范学院数学建模训练题

中国是一个严重缺水的国家。淡水资源总量为28124亿立方米，占全球淡水资源的6%，仅次于巴西、俄罗斯和加拿大，居世界第四位，但人均只有2100立方米，仅为世界平均水平的28%、美国的1/5，在世界上名列121位，是全球13个人均水资源最贫乏的国家之一。

中国水资源地区分布也很不平衡，长江流域及其以南地区，国土面积只占全国的36.5%，其水资源量占全国的81%；其以北地区，国土面积占全国的63.5%，其水资源量仅占全国的19%。我国水资源分布特点为：南多北少，东多西少，总量丰富，人均占有少，分布不均匀。

目前，中国年用水总量已突破6000亿立方米，全国年平均缺水量500多亿立方米，三分之二城市缺水。经水利部预测，202_年中国人口将达到16亿，届时人均水资源量仅有1750立方米。在充分考虑节水情况下，预计用水总量为7000亿至8000亿立方米，要求供水能力比现在增长1300亿至2300亿立方米，全国实际可利用水资源量接近合理利用水量上限，水资源开发难度极大。

根据202_ 年中国环境公报显示，全国地表水污染依然较重。长江、黄河、珠江、松花江、淮河、海河和辽河等七大水系总体为轻度污染。204 条河流 409 个地表水国控监测断面中，Ⅰ～Ⅲ类、Ⅳ～Ⅴ类和劣Ⅴ类水质的断面比例分别为 59.9%、23.7%和 16.4%。主要污染指标为高锰酸盐指数、五日生化需氧量和氨氮。26 个国控重点湖泊（水库）满足Ⅱ类水质的 1 个，占3.8%；中，Ⅲ类的 5 个，占19.2%；Ⅳ类的 4 个，占 15.4%；Ⅴ类的 6 个，占 23.1%；劣Ⅴ类的 10 个，占 38.5%。主要污染指标是总氮和总磷。从数据分析，我国湖泊（水库）的水质中 V 类占比高，且三湖中，除太湖环湖河流总体为轻度污染外，滇池巢湖环湖河流总体均为重度污染。

根据202_ 年中国环境公报显示，全国地下水质量状况不容乐观，水质为优良-良好-较好级的监测点总计为 1759 个，占全部监测点的 42.8%，2351 个监测点的水质为较差-极差级，占全部监测点的 57.2%。我国地下水水质较差的地区占比很高，且华北、东北、西北地区以及华东、中南、华南地区部分城市的水质状况继续变差。

随着经济、社会、人口的持续增长，水资源成为制约中国发展的重要因素。请在中国选择一个有代表性的地区，充分收集其有关水资源的相关数据，建立一个数学模型来制定一个有效的、可行的和具有成本效益的202_年水资源战略，以满足该地区直到202_年的用水需求。数学模型中必须解决水资源的存储、流动、环境污染（比如盐碱化、重金属污染等，考虑其中一种即可）和保护。如果可能的话，用你的模型，来探讨你的战略对经济和环境的影响。并解释它为什么是最好的水战略选择

第二篇：赣南师范学院自荐书格式

尊敬的各位领导、老师：

您们好！

本人是赣南师范学院数学与计算机科学学院数学专业毕业生，经慎重考虑欲加盟贵学校。

在大学三年里，时刻按照“宽专业、厚基础、强能力、高素质”标准去锻炼及发展自我，在不断的学习和实践中提高了自己的综合素质，已把自己塑造成为一个专业功底扎实、知识结构完善、适应能力强、富于协作精神的时代青年。

我相信未来社会需要的是高素质的复合型人才，成功的学习者在充分的认识到书的价值的同时，也应认识到书的无价值。因而我在学习之外，积极参加了各种各样的课余活动，如“计算机协会”、“数学建模”等。所有这些活动都有利于我提高自身的计算机操作能力和团体协作能力。20xx下半学期在江西省会昌县第二中学教育实习一个多月，在教学和班主任工作两个方面得到较好的锻炼，掌握了教学各个环节的基本要求和方法，了解了班主任工作的主要内容及其意义，具备了独立从事教学工作和班主任工作的能力。在实习期间表现出色，受到当地教师、领导的认可，被学校评为“优秀实习生”。在班主任管理方面及和学生沟通方面获益匪浅，管理好一个班不仅仅靠能力，更重要的是用心，有激情，只要拥有朝气蓬勃的活力才能更好的把工作做好。

我坚信命运之神只垂青于有准备的人，目睹过去深感母校的培养恩深，注目将来惟有以热血、真诚、眼泪和汗水回报社会的赠与，让生命之烛高举奋斗之光！

我将贵校作为首选目标，是因为贵校拥有优美的教学环境和积极进取的拼搏精神。我相信在贵校领导的帮助指导下，我一定会学得更多，做得更好。

相信您的眼光，相信我的选择。

诚候佳音！

此致

敬礼！

自荐人：曾xx

20xx年xx月xx日

第三篇：202_数学建模A,B题

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

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

A题

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

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

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

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

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

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

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

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

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

B题

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

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

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

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

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

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

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

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

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

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

第四篇：202_数学建模A题

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

摘 要：

嫦娥三号卫星采用的是软着陆方式登陆月球，在卫星高速飞行的情况下，我们要精确地在月球预定区域内实现软着陆，需要对其运行轨道进行设计并制定相应控制策略。由于天体的运动均满足开普勒三大定律以及总能量守恒定律，我们据此建立一系列的方程，最终求得卫星在近月点处的速度大小Vab/(ac)*GM/a1.6922km/s 求得卫星在远月点处的速度大小Vbb/(ac)*GM/a1.6139km/s。其速度方向均为当前运动轨道的切方向。

嫦娥三号的软着陆过程又分为6个阶段，对每个阶段都要进行一定的控制。通过对数值的分析，我们制定如下策略：在主减速阶段，首先卫星的发动机应全部用于水平方向的减速，竖直方向以自由落体状态加速下降，直至竖直方向速度达到56.8505m/s，水平方向速度为4.9522m/s，耗费燃料90kg，卫星下落高度约为14695m。之后发动机推力仍以最大输出工作，其推力在竖直方向上瞬间产生约3752N的向上阻力，在水平方向上产生6494N的阻力，使得卫星在主减速的后半阶段在竖直方向开始匀速运动，水平方向继续做减速运动水平速度最终降为0，直至卫星降落到距月3000m高度处；在快速调整阶段，姿态发动机调整卫星的运动方向为竖直向下，主发动机继续工作，发动机在水平方向上产生1097N推力，竖直方向上产生3752N推力，使得卫星水平方向的速度降为0m/s，竖直方向继续以57m/s的速度做匀速运动；在粗避障阶段，卫星对其正下方月面进行拍照获得数字高程图，利用matlab软件对图像进行灰度色差分析，进而初略确定了一个着陆点的像素点坐标(193，1169），过程中，姿态发动机根据图像进行卫星的位置调整，主发动机用于竖直方向上的减速工作，使得卫星在距月100m处达到悬停状态；精避障阶段，进行更精细的月面拍摄，采用同粗避障段类似方法得到更精确的像素点坐标为(318，651),同时利用姿态发动机调整运动方向的同时主发动机工作对水平方向进行速度控制，使其在距月30m处水平速度为0m/s；缓慢下降阶段，发动对卫星竖直方向进行减速控制，使其在距月4m处合速度变为0m/s；在最后阶段，关闭发动机，使其做自由落体运动，最终成功着陆。

关键词：

软着陆 开普勒定理 优化控制 基于C语言编程

一、问题重述

嫦娥三号于202_年12月2日1时30分成功发射，12月6日抵达月球轨道。此次嫦娥三号采用的是软着陆方法。

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

根据题目给出的已知量建立合适的模型求解分析下列几个问题：

（1）确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。

（2）确定嫦娥三号的着陆轨道和在6个阶段的最优控制策略。

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

二、问题分析

本问题的一个关键点在于求解出嫦娥三号卫星的绕月运行轨道。由于已知卫星的运行轨迹为近月点15km，远月点100km的椭圆形轨道。

对于第一问的关键就是求解出该椭圆轨道与月球赤平面的夹角。由于已知落点的经纬度坐标和软着陆过程中6个状态的数值，我们就可以使用逆推法求解出着陆准备轨道近月点和远月点的经纬度坐标位置。由于卫星的运动满足开普勒三大定律以及能量守恒定律，我们可以根据这几大定律确定几个等式联立方程组从而求得嫦娥三号在近月点和远月点的相应速度大小与方向。

对于第二问，要设计一种方案使得在软着陆的过程中耗能最少，并达到预设的各项指标。我们需要不断的去设计、计算、调整，在不断的尝试摸索中寻找出一个比较不错的软着陆方案。

三、模型假设及符号约定

3.1 模型假设（1）假设卫星在主减速阶段的运动轨迹为抛物线

（2）假设卫星在3000m处基本就处于目标地点的竖直上方（3）假设月球为球体

（4）假设在近月点即将开始软着陆状态的前后速度基本不变

3.2 符号约定 符号

含义 G

引力常量

M

月球的质量

m(t)

t时刻卫星的质量 月球的半径

月球表面的重力加速度 r

g'

Vi

ai

Si

Ei

Ti

Ft Ve i位置处卫星的速度大小 i方向的加速度

i位置处相对应的面积大小 i位置处卫星的能量 到i位置所用时间 t时刻发动机的推力 以m/s为单位的比冲 单位时间燃料消耗的质量 椭圆的半长轴 椭圆的半短轴 椭圆的焦点 m a

b

c

四、模型建立与求解

4.1 模型建立

4.1.1 问题1卫星轨道模型的建立

我们将嫦娥三号卫星绕月运行的椭圆轨道抽象出如下图所示的一个简单的几何图形，月球的月心位于椭圆的一个焦点F上，椭圆的半长轴为a，半短轴为b，半焦距为c。A点为近月点，速度为Va，B点为远月点，速度为Vb。易知A、B两点距月心的距离Laac，Lbac。在一个极小的时间段t内，卫星与月心连线扫过的面积分别为SaVa/2*t*La，SbVb/2*t*Lb。由开普勒第二定律可知，卫星与中心天体连线在单位时间内扫过的面积相等，所以SaSb，代入化简后可得公式(1)Vb(ac)/(ac)*Va；由于卫星运动的总机械能等于其动能和引力势能之和，所以在A点，卫星的总机械 能Ea1/2*mVa2(GMm/La)（公式2），同理B点的总机械能。卫星在运行过程中只有动能和引力势能Eb1/2*mVb2(GMm/Lb)（公式3）之间的转化，机械能守恒，所以EaEb（公式4）。[1] 由上面(1)(2)(3)(4)四个公式我们可以解得Vab/(ac)*GM/a，Vbb/(ac)*GM/a。

4.1.2 问题2卫星软着陆各阶段控制策略模型的建立

主减速段控制策略：

第一阶段，主减速发动机全力用于水平方向减速，竖直方向呈自由落体加速下落，至下落速度达到V1y(m/s)为止，用时记为T1，质量减少到M1；

V1yg'*T1(1)当0tT1时，由Fm*Ve可得出下式，tm(t)2400mdt024002.5510t(2)三号卫星运行速度V水平方向分量为：

(3)在0tT1内，卫星下落的高度ht0.5*g'*t^2, 水平方向位移

t s(t)(170002940ln(2400)2940ln(24002.2210t))dt

(4)第二阶段，竖直方向以V1y匀速下降，主减速发动机仍以最大推力工作，F(t)的竖直方向分力为m(t)*g'，水平方向分力为Fx(t)75002(m(t)*g')2。

记第二阶段结束时刻，即距月面3000m时为T2，当T1tT2时，m(t)24007500*t/294024002.5510t(5)下落高度 h(t)0.5*g'*T12V1y*(tT1)(6)水平方向加速度

2ax(t)fx(t)/m(t)(75002(24002.5510*t)2*g')/(24002.5510t)(7)水平方向速度

ttVx(t)Vx(T1)T1222ax(s)dsVx(T1)(7500(24002.5510s)g')/(24002.5510s)dsT1 1=Vx(T1)2.5510gg'(24002.5510t)'(24002.5510T1)(75002u2/u)du

17500=Vx(T1)(75002u27500ln2.5510t75002u2u

(9)

'(24002.5510T1))gg'(24002.5510t)(8)水平方向位移 s(t)s(T1)快速调整段控制策略：

Vx(t)dt T1

在快速调整阶段，卫星仍然以Vyg'*T157m/s的速度下落，水平方向推力逐渐减小到0，这样卫星的方向就正对目标下方了。

在tT2时，发动机给与卫星水平方向一个力Fx，使得其水平方向的速度最终降为0m/s；同时主发动机在竖直方向上也对卫星施加一个力，使其在竖直方向上做匀速运动。水平方向位移Sx(T3)V0T30.5*a*T3^2。

水平方向推力为Fx(t)m(t)*ax(t)，合力为

gm1ve2tT2m(r)dr

两边对t求导得：

m'(t)gm1ve2m(t)

再由tT2时，m(t)m(T2)得，当T2tT24.96时：

m(t)m(T2)exp(gm1ve2(tT2)

水平方向速度

vx(t)vx(T2)ax(r)drvx(T2)tT2T2t

水平方向位移

s(t)s(T2)vx(r)drs(T2)(vx(T2)T2)(tT2)0.5(t2T2)T2t2

(s(T24.96)，15000)就是近月点坐标。

粗避障段控制策略：

以下步骤的实施均在matlab上实践，源代码见附录2 将抓拍到的图片利用matlab读入，求出图像平均灰度，以此来表示理想降落地点的高度。但通过比较发现平均灰度并不能代表理想地点高度，于是求出灰度众数average，作为理想降落高度。

然后新建图像，使新图像上的每点的灰度等于原图像对应点的灰度与灰度众数average的差的绝对值。于是新图像上灰度越高的点表示越不理想的点。（matlab运行效果图见附录3）

选取10 作为分界线，灰度值大于10的改为150，表示不理想的降落地点，灰度值小于10的改为0，表示理想的降落地点。将图片边界四周都像素点全部改为灰度150，表示不理想降落地点，因为越靠近图片边缘则越缺少信息。(matlab运行效果图见附录4)接下去通过将不理想的点向四周扩散来逐渐覆盖全图，留下0.001的空隙时停止扩散，则此时剩余的空隙则为相对理想的降落地点，求出这些降落地点中离中心飞船最近的地点，作为最终降落地点。

最后我们分别绘制了覆盖率为0.7时的图像、覆盖率为0.99的图像、覆盖率为0.999的图像(分别件附录5、6、7)。

精避障段控制策略：

此阶段只是在粗壁障阶段的基础上对着落点做了一个更精确的定位，所以选点的原理同上。

4.2 模型求解 4.2.1 问题1

由题目已知，卫星近月点为15km，远月点为100km，据此我们可以得出下面一个等式c15(cr)100，求解得c9r/2 km。据此又可以导出椭圆的半长轴的大小acr15243/2*r，半短轴ba2c2。然后将a、b、c等数值代入模型中，便可以直接求出卫星在近月点的速度大小Va1.6922km/s，卫星在远月点的速度大小Vb1.6139km/s。

4.2.2 问题2 主减速阶段：

由天体运动黄金代换式GMR2g'，得出g'GM/R21.6243m/s2 利用C语言编程使用迭代法搜索T1的值： 具体编程算法见附录1 Step 1 令t=T1, 下落高度 h(t)0.5*g'*t^2, deltt=1 竖直方向速度Vy(t)g'*t, 水平方向速度（初始速度有第一问可知为1692m/s）

Vx(t)16922940ln(2400)2940ln(24002.2210t)水平方向位移

S(t)(16922940ln(2400))tStep 2 令t=t+ deltt

29402.2210t(ulunuu)2400 24002.2210下落高度 h(t)0.5*g'*T1^2Vy*(tT1)竖直方向速度Vy(t)g'*T1 水平方向速度由（8）式计算可得 水平方向位移

Step 3 如果h(t)120002641 and Vx(t)2 如果h(t)120002641 and Vx(t)step 1 如果h(t)120002641, then end

最后利用C语言编程解得

572Vy(t)2, then go to step

572Vy(t)2, then T1T11, goto T135s,h14695.85m,Vx4.9522m/s,Vy56.8505m/s。

从而我们可以得出，在主减速阶段，首先卫星的发动机的全部动力要用于卫星水平方向上的减速，竖直方向以自由落体状态加速下降，直至竖直方向速度达到56.8505m/s，水平方向速度为4.9522m/s，此时卫星的质量约为2310kg，耗费燃料90kg，卫星下落高度约为14695m。之后发动机推力仍以最大输出工作，其推力在竖直方向上瞬间产生约3752N的向上阻力，在水平方向上产生6494N的阻力，使得卫星在主减速的后半阶段在 竖直方向开始匀速运动，水平方向继续做减速运动水平速度最终降为0，直至卫星降落到距月3000m高度处。

快速调整阶段：

根据分运动的等时性，可以利用卫星竖直方向上的运动方程hV0t，求得卫星在快速调整阶段的最大用时th/V010.52s，从而根据水平方向的运动公式VVoa*t求得水平方向的最小加速度a0.475m/s。所以可以得出该过程中发动机在水平方向上的最小推动力Fxmin1097N，竖直方向上的推力Fxymg'3752N，并求得水平位移S15234m。

五、模型评价与推广

5.1 模型评价

优点：本模型采用了化繁为简的方法，将复杂问题简单化，把卫星绕月运行抽象出一个椭圆的几何问题再结合物理模型利用数学方法求解相应的运动数值

缺点：本模型很多地方理想化，比如忽视了月球的扁率，卫星运行中姿态是变化的某些部分不应该用质点来处理

5.2 模型推广

利用本模型的算法可以求出任何天体在每个时刻的运行速度及其他物理量，方便研究天体运动，方便针对性的对需要探测的星球设计相应的卫星轨道。

六、参考文献

[1] 王健伟 李兴，近日点和远日点速度的两种典型求法，物理教师，第34卷第6期：58，202_

七、附录

附录1：

下面采用的是C语言程序来摸索求解T1的值（运行环境VC++ 6.0）文件夹内find_t1.cpp文件

#include #include #include #define g 1.6243

int main(){

int t,t1,deltt=1;

t=t1=4;

double h,vx,vy,compare,u1,u2;

h=0.5*g*t*t;

vx=1692-2940*log(2400)+2940*log(2400-2.2210*t);

vy=g*t;

while(1)

{

compare=sqrt(57*57-vy*vy);

if(h<12000+2641 && vx>compare)

{

t=t+deltt;

u1=g*(2400-2.5510*t1);

u2=g*(2400-2.5510*t);

vy=g*t1;

h=0.5*g*t1*t1+(t-t1)*vy;

vx=1692-2940*log(2400)+2940*log(2400-2.2210*t)-(1.0/2.5510)*(sqrt(7500*7500-u1*u1)+7500*log((7500-sqrt(7500*7500-u1*u1))/u1)-(sqrt(7500*7500-u2*u2)+7500*log((7500-sqrt(7500*7500-u2*u2))/u2)));

}

else if(h<12000+2641 && vx

{

t1++;

t=t1;

h=0.5*g*t*t;

vx=1692-2940*log(2400)+2940*log(2400-2.2210*t);

vy=g*t;

}

else if(h>12000+2641)

{

printf(“t1=%dn”,t1);

printf(“h=%lfn”,h);

printf(“vx=%lfn”,vx);

printf(“vy=%lfn”,vy);

break;

}

}

return 0;}

附录2：

big=imread('附件3 距2400m处的数字高程图');

[m,n]=size(big);

notenum=[1:256];

for i=1:m;for j=1:n;notenum(big(i,j)+1)=notenum(big(i,j)+1)+1;end end

max=0;maxnum=0;for i=1:256;if max

big1=big;for i=1:m for j=1:n big1(i,j)=abs(int8(big1(i,j))-average);end end

big2=big1;

for i=1:m for j=1:n if big2(i,j)>15 big2(i,j)=150;else big2(i,j)=0;end

end end

for i=1:m;big2(i,1)=150;end

for i=1:m;big2(i,m)=150;end

for j=2:n-1;big2(1,j)=150;big2(m,j)=150;end

for k=1:200 num=[m:n];

numsum=0;for i=1:m;for j=1:n;if big2(i,j)>20;num(i,j)=1;numsum=numsum+1;else num(i,j)=0;end

end end

if numsum/(m*n)>0.999 break end

for i=2:m-1;for j=2:n-1;sum=0;if num(i,j)==1;continue else

for i1=-1:1;for i2=-1:1;if num(i+i1,j+i2)==1;sum=sum+1;end

end

end

if rand(1)*3

end

end end end

distance=[m,n];for i=1:m;for j=1:n;if num(i,j)==0;distance(i,j)=(i-m/2)*(i-m/2)+(j-n/2)*(j-n/2);else distance(i,j)=m*m+n*n;end end end

min=m*m+n*n;mini=0;minj=0;

for i=1:m;for j=1:n;if num(i,j)== 0;if min > double(distance(i,j));min=double(distance(i,j))mini=i;minj=j;end end end end

附录3

附录4

附录5 覆盖率为0.7的图像

附录6 覆盖率为0.99的图像

附录7 覆盖率为0.999的图像

第五篇：赣南师范学院毕业生求职信

在求职中，不管是什么工作，都需要附带一份求职信，以下是一封求职信范文，请参考；

尊敬的各位领导：

您们好!

本人是赣南师范学院数学与计算机科学学院数学专业毕业生，经慎重考虑欲加盟贵学校。

在大学三年里，时刻按照“宽专业、厚基础、强能力、高素质”标准去锻炼及发展自我，在不断的学习和实践中提高了自己的综合素质，已把自己塑造成为一个专业功底扎实、知识结构完善、适应能力强、富于协作精神的时代青年。

我相信未来社会需要的是高素质的复合型人才，成功的学习者在充分的认识到书的价值的同时，也应认识到书的无价值。因而我在学习之外，积极参加了各种各样的课余活动，如“计算机协会”、“数学建模”等。所有这些活动都有利于我提高自身的计算机操作能力和团体协作能力。202_年下半学期在江西省会昌县第二中学教育实习一个多月，在教学和班主任工作两个方面得到较好的锻炼，掌握了教学各个环节的基本要求和方法，了解了班主任工作的主要内容及其意义，具备了独立从事教学工作和班主任工作的能力。在实习期间表现出色，受到当地教师、领导的认可，被学校评为“优秀实习生”。在班主任管理方面及和学生沟通方面获益匪浅，管理好一个班不仅仅靠能力，更重要的是用心，有激情，只要拥有朝气蓬勃的活力才能更好的把工作做好。

我坚信命运之神只垂青于有准备的人，目睹过去深感母校的培养恩深，注目将来惟有以热血、真诚、眼泪和汗水回报社会的赠与，让生命之烛高举奋斗之光!

我将贵校作为首选目标，是因为贵校拥有优美的教学环境和积极进取的拼搏精神。我相信在贵校领导的帮助指导下，我一定会学得更多，做得更好。

相信您的眼光，相信我的选择。

此致

敬礼

申请人：*** xx年**月**日