2015年全国大学生数学建模竞赛A题分析报告.pdf

太阳影子定位技术问题的数学模型 摘要 本文涉及的是太阳影子定位技术问题。在已知视频中物体的太阳影子变化的情况下，要 确定视频的拍摄地点和拍摄日期。首先，分析了文中四个问题的关系，发现前三个问题的已 知条件逐步减少，问题难度依次递进。第四问则给出一个实际问题，该问题需要转化成数学 模型利用前三问的方法求解；随后，建立了 L-G 模型、MinZ-模型等，并应用非线性最小二 乘法、遗传算法等算法对模型求解。得到基于模型的合理结果。最后，将第四问的实际问题 转化数学模型并求解，进而解决问题。 对于问题一，要解决的问题是杆长与影子长度的关系，根据天文、几何知识，我们建立 了模型来刻画问题给出的参数之间联系，如赤纬角模型、时角模型、太阳高度角模型、影子 长度模型（L-G 模型）等；分析了各参数对影子长度的影响；最后运用MATLAB 绘制出具体 给定参数下的3 米高直杆的影子变化曲线；从曲线可以看出在9:00 到 15:00 这段时间里， 影子长度先变短后变长，最短为3.627 米，最长为7.182 米。 问题二提供了一个关于时间、影子坐标的附件1，杆长未知，为了确定直杆所处的地点， 本问建立了MinZ-模型，首先将经度、纬度、杆长离散化，搜索出大概的可行解，然后运用 非线性最小二乘算法，选取 matlab 中的lsqcurvefit 命令，以可行解为初值，再运用非线性 最小二乘算法，选取MATLAB 中的lsqcurvefit 命令，在控制残差在10−8之内范围的情况下得 到了三个可能地点皆在海南省昌江县内，最小误差的地点为海南省江黎族自治县，北纬 19.3025°，东经108.6988°，此时对应直杆高度为2.0219m。同时，将结果代入问题一的模 型进行检验，验证了模型的稳定性和算法的合理性。 问题三沿用问题一的模型和问题二的算法，由于一个已知量变成一个变量，根据算法特 点，在增加一个变量的情况下，算法搜索影长差时只需要增加一重循环。关于附件2 数据， 残差最小对应的位置为北纬39.8926°，东经79.7438°，具体地点在新疆维吾尔自治区喀什地 区巴楚县。日期为2015 年4 月24 日或6 月19 日；关于附件3 数据则是残差最小对应的位 置为北纬31.3625°，东经110.1602°，具体地点在湖北省神农架林区。此时对应直杆高度为 2.9871m。对应日期为1 月25 日或10 月5 日。 问题四是将实际问题转化成数学模型并求解的过程，用aviread 函数将视频读入MATLAB 软件，并用size 函数读取视频帧数Vf=61000 ，利用mat2gray 函数可以对图像灰度化处理。 取合理阈值使影子与背景的分界清晰可见，记录直杆顶点及底端坐标，然后依次测录各图片 中影子顶点的坐标（单位为像素），最后应用问题二和问题三的方法求解，结合运用MATLAB 遗传算法工具箱，在给出日期条件下，残差最小对应的位置为北纬40.3344°，东经113.2556°， 具体地点为内蒙古自治区鄂尔多斯市达拉特旗县；若日期未知，也得到了比较合理的视频拍 摄地点和日期。 论文最后做了误差分析并给出了模型改进意见；论文的特色在于思路清晰，方法简洁， 将数学模型与计算机算法、图形很好的结合起来，用图形图表体现结果。每个问题都进行了 结果分析和模型验证。 关键词：影子定位 L-G 模型 minZ-模型 非线性最小二乘法 遗传算法 一、问题的重述 太阳影子定位技术是通过太阳影子变化来确定视频拍摄地点和拍摄时间，现以该项技术 为基础建立数学模型，完成下列任务： 1、分析影子长度关于各个参数的变化规律，并画出2015 年 10 月22 日北京时间 9:00-15:00 之间天安门广场（北纬39 度54 分26 秒,东经116度23 分29 秒）3 米高的直杆 的太阳影子长度的变化曲线。 2、根据附件 1 所给的影子顶点坐标数据，建立数学模型确定直杆所处地点，给出若干 个可能的地点。 3、根据附件2、附件3 所给的影子顶点坐标数据，建立数学模型确立若干个可能的地 点与日期。 4、附件4 为一根直杆在太阳的影子变化下的视频，已知所给直杆高度为2 米，请根据 视频中影子的变化，确定若干个可能的拍摄地点。如果拍摄日期未知，请根据视频确定日期 和拍摄地点。