转载 2015A国赛优秀论文

1
太阳影子定位
摘要
本题意在通过分析实际测量的竖直杆影子随时间变化的坐标或视频数据，确
认竖直杆所在的经纬度以及数据对应的记录日期。
第一问要求分析画出已知杆长、地理位置和日期的竖直杆的太阳影子变化曲
线，并分析影子长度关于各个参数的变化规律。本文利用天文学公式，建立了该
竖直杆太阳影子关于时间变化的计算求解模型，得出影长随时间的变化曲线如下。
此后，本文使用控制变量法，分别对经纬度、日期、杆长及当日时间这五个参数
对影子长度的影响展开讨论。经过误差分析并改进算法后，得到在11 点59 分24
秒，取到最短的影长为3.6638 米
第二问要求分析某竖直杆的太阳影子顶点坐标数据，确定其地点。本文沿用
第一问的天文学公式，以每一对时刻理论与实测的影长之比的差的平方和最小为
首要目标，以每一对时刻理论与实测方位角差值之差的平方和最小为次要目标，
建立多目标规划模型，采用分层求解法，遍历求解得到最能满足条件的位置为北
纬19 度，东经109 度，位于我国海南省。
第三问要求分析两组某竖直杆的太阳影子坐标数据，确定其地点及数据记录
日期。与第二问类似，在添加日期这一变量后，建立与第二问相同的多目标规划
模型，采用分层求解法遍历求解。得到附件二和附件三数据对应的最能满足条件
的位置分别为：
纬度 经度 日期 大致地点
附件二 40 79 5 月25 日 新疆图木舒克市
附件三 33 106 10 月31 日 陕西，汉中广元之间
第四问要求分析竖直杆的影子变化视频，分别在日期已知和未知的情况下求
出其所在地点与拍摄日期。本文首先利用动态追踪技术，找出影子顶点在每一时
刻的像素坐标，计算出每时刻的太阳高度角。进而利用天文学公式中太阳高度角
与日期、纬度、经度的确定关系，分别在日期已知和未知的情况下建立计算求解
模型，考虑到模型存在系统误差，难以得到完美解，故将问题转换为非线性规划
问题遍历求解。得到最能满足条件的位置与日期为：
纬度 经度 日期 大致地点
日期确定 43 115 7 月13 日 内蒙古锡林郭勒盟西南
日期未定 44 113 6 月25 日 内蒙古锡林郭勒盟西
本文亦对文中使用的天文学公式的误差及所求出的结果进行了讨论，在给出
最优解之外还讨论了其他解存在的可能性。
关键词：计算求解模型、控制变量法、多目标规划、分层求解法、误差分析
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
2
1. 问题重述
1.1 问题背景
现代科技的发展使得人们能够更为方便地记录高质量的视频文件。在分析视
频材料时，有时需要确定视频的拍摄地点及日期，而利用天文学知识，对视频物
体中的太阳影子变化进行分析是确定视频拍摄地点及日期的一种方法。
1.2 相关信息
本题给出了三组某处某固定直杆在水平地面上太阳影子的顶点坐标数据，每
组数据皆包含 21 个等间距时刻（北京时间）直杆太阳影子的x, y顶点坐标。其
中，第一组数据还额外包括了数据对应的日期。坐标系以直杆底端为原点，水平
地面为x, y平面，直杆垂直于地面。
除三组数据以外，本题亦给出了一根直杆在太阳下的影子变化视频，视频中
清晰的记录了2015 年7 月13 日8 点54 分06 秒到同日9 点34 分36 秒某地高
为2 米的直杆的太阳影子变化过程。以上信息中，各直杆的地理位置皆未知。
1.3 需要解决的问题
1) 建立太阳影子变化的数学模型，分析影子长度关于各个参数的变化规律并应
用这一模型画出2015 年10 月22 日北京时间9:00-15:00 之间天安门广场（北
纬39 度54 分26 秒，东经116 度23 分29 秒）3 米高的直杆的太阳影子长度
变化曲线。
2) 建立数学模型，根据题目提供的附件一中的影子顶点坐标数据，给出该直杆
所在的可能地点。
3) 建立数学模型，根据题目提供的附件二、三种影子顶点坐标数据，给出对应
直杆所在的可能的地点。
4) 分析题目提供的视频，确定视频拍摄地点的数学模型，分别讨论在拍摄日期
已知和未知两种情况下，能否确定视频的拍摄地点与日期，并给出该视频可
能的拍摄地点。
2. 模型假设
假设一：本文研究的所有对象所在的地面皆为平地。
假设二：本文研究的所有对象所处地海拔为0。
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
春哥给大家定期发资料啦，所有符合要求的，会被标记“定期发资料”，然后春
给会定期发布内部课讲义、各种数学建模资料、以及春哥独家整理的内部资料。
资料主要包括：优秀论文、优秀论文框架分析、排版格式问题汇总、各种论文
模版、各种数学建模教材PDF 版、各种算法源代码以及源代码使用说明以及历
年国赛赛题分析课件等资料。
（1）VIP 学员、传课学员（VIP 学员有获奖保证）
（2）关注微信公众号：shumohome，并回复“分享者”按照要求操作即可。
注：资料包括春哥分析整理的资料等，因为资料性价比非常高，所以请大家不要
分享给别人。
有其他疑问可以添加QQ 群：598397784 咨询
郑老师
郑老师被学员亲切的称为春哥，精通数学建模思路分析，数学建模思维培养，建
模论文写作等方面。所指导过的学员70%以上获得了省级二等奖及以上的奖项，
是
春哥带你学建模的所有课程按照数学建模学习四步走进行授课：
（1）排版与格式——在国赛评阅中，评阅用时非常短，阅卷者会首先看一篇论
文的排版与格式是否正确。虽然说评阅过程摘要最重要，但是本课程会提供摘要
模板，可以直接在比赛中掌握。
（2）论文与摘要——论文的重要程度，不言而喻，如何写论文，是非常重要的，
本课程从两个方面讲解如何写数学建模论文，一个是完全没有思路如何编凑论文，
一个是如何写出一篇完美的论文。
（3）解题思路——依照讲师整理的算法模型简介，以及优秀论文，对数学建模
赛题进行分析，让学员学会对于各种类型题目都能迅速分析出最佳算法。
（4）方法积累——讲师整理各类算法及其在数学建模中的需要掌握程度，让你
在最短的时间掌握一种算法，例如两分钟掌握层次分析在数学建模中的应用
3
3. 符号说明
符号 符号说明
i 
时刻i 的太阳高度角
 太阳赤纬
d 代表天数，当日期为 1 月 1 日时，d=1
 目标对象所在地的纬度。为正数时代表北纬，负数为南纬
l 目标对象所在位置经度。为正数时代表东经，负数为西经
i t 时刻i 下研究对象对应时角
L 直杆长度
yi L 直杆在时刻i 时的影长
’
yi L 直杆在时刻i 的理论影长
i 
直杆影子在时刻i 被测得的方位角
’
i 
直杆影子在时刻i 的理论方位角值
1 f 每一对理论的影长之比的差的平方和
2 f 方位角差值之差的平方和
ls t 本地太阳时间
lt t 本地时
i0 t i 时刻i 时子午线当地标准时间
E 时差
4. 问题分析
问题（一）
第一问要求建立影子长度变化的数学模型，分析影子长度关于各个参数的变
化规律，并画出题目指定的直杆的太阳影子长度变化曲线。由于本题要求研究太
阳影子以定位，结合实际生活体验可知，对于一个自身长度固定的竖直杆，其影
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
4
子长度不仅仅与其所在地理位置有关，还与地球自转情况、地球绕太阳公转的情
况有关。因地球自转与公转分别导致了昼夜交替及四季变化，故影子长度应该还
与对应的当天时刻和日期有关。
考虑到这些因素，可进一步搜索相关的文献，建立数学模型表现影长与地理
位置、时间的关系，通过研究关系式以得出题目要求求解的答案。
问题（二）
第二问要求根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数
学模型确定直杆所处的地点。本问相对于第一问，直杆长度未知，但因知道太阳
影子顶点在一个连续时间段内的数据，故相当于知道了此时间段内影子的影长以
及其角度变化情况。经过第一问的分析与解答可知，影长与直杆长度、当天日期、
时间、经纬度此5 个因素有关；本问中，日期、时间皆为已知量，而特殊地，当
其他因素恒定时，直杆长度与影长呈正比关系，故可以研究两个时刻的影长比值，
消除直杆长度这一未知量。这一比值仅仅与直杆所在的经纬度有关，因而能够达
到减少未知数的目的。
除了分析影长以外，角度也是一个不可忽略的因素。虽然影长信息已经足以
解出未知数，但再增加角度信息，能够提高方程组冗余度，使得模型的稳健性和
抗噪性得以提升。
用来表示角度的量应当能够体现影子的方位，参考资料可知，因方位角能够
表现影子在地平面上的方向，使用方位角[1]这一指标较为合适。在本题中，虽然
每个时刻的方位角的值难以得到，但是两个时刻的方位角差却与坐标系自身无关，
可以根据两个时刻的影子顶点在同一坐标系中的相对位置计算。
因而，综合考虑影长和角度两个指标，以直杆所在经纬度为x, y时，计算所
得两时刻的影子长度比与对应时刻测得的影子长度比的差值最小为第一个目标
函数，以计算所得两时刻方位角差值与测得的方位角差值的差值最小为第二个目
标，建立多目标规划模型，找出可行的直杆对应的经纬度值。
问题（三）
第三问与第二问类似，不同的是日期未知，要求根据某固定直杆在水平地面
上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点与日期。与第二
问相同，可综合考虑影长和角度两个指标，以直杆所在经纬度与日期分别为x, y, z
时，计算所得两时刻的影子长度比与对应时刻测得的影子长度比的差值最小为第
一个目标函数，以计算所得两时刻方位角差值与测得的方位角差值的差值最小为
第二个目标，建立多目标规划模型，找出可行的直杆对应的经纬度值与日期。
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
5
问题（四）
第四问要求分析题目提供的视频，确定视频拍摄地点的数学模型，分别讨论
在拍摄日期已知和未知两种情况下，能否确定视频的拍摄地点与日期。因此，需
要先对附件四的视频进行处理，找出正确的太阳高度角在拍摄时间段内各个时刻
的值。
因第一、二、三问建模过程中，确定经纬度并不一定需要直杆长度以及影长
这两个量同时已知，需要的仅仅是由之计算得到的太阳高度角在不同时刻的值，
因而得到此信息后可将问题转换为与前三问类似的问题。在前三问的建模过程中，
使用了关联太阳高度角、日期、纬度、经度、及时间的等式，故可利用该等式建
立计算求解模型求解。
5. 数据分析
5.1 附件四视频处理
本题附件四是一根直杆在太阳下的影子变化视频，视频中清晰的记录了2015
年7 月13 日8 点54 分06 秒到同日9 点34 分36 秒某地高为2 米的直杆的太阳
影子变化过程。直杆的地理位置皆未知。为了使用该视频中的信息解答第四问，
需要对视频做一定的处理，找出各个时刻直杆影子尖端部分的变化。
图 5.1.1 附件 4 视频截图
图 5.1.1 为附件4 的视频截图。因已知竖直杆垂直于地面，观察该视频发现：
竖直杆的中心位于画面中央，可认为摄像机的轴线是与直杆垂直的。随着时间的
推移，可以用肉眼观察到影子的长度变短，因此可以推断该视频拍摄在当地正午
时间前，则可以作为左上角的时间是该地所在时区的时间的佐证。
此外，由于在视频所记录的时间里，直杆影子转动的角度十分小，而在图示
画面中，可认为影子所在的直线是和摄像机的视线垂直的，因此可以认为影子长
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
6
度在透视矫正前后变化不大，因而不需要在处理视频的时候对视频做透视矫正，
即认为影子和直杆构成的平面垂直于摄像机的视线。
图 5.1.2 附件4 处理示意图
利用 Adobe After Effect 软件读取该视频，采用软件自带的动态追踪技术读取
视频每一帧影子顶部的像素坐标C、直杆顶端的坐标A 以及影子上的任一其他
点的像素坐标 B由此，即可得到同一坐标系内的两条向量BC, AC，进而解得两
向量的夹角大小，用于之后的各项运算。
图5.1.3 动态追踪操作示意图
动态追踪 C 点的具体操作如图所示，图中有一大一小两个白边方框，小方框
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
7
中是追踪的目标即影子顶点及周围的地面。由于影子顶点具有固定的形状与较深
的颜色，小白框内的像素的颜色信息有一种特征明显的分布，故可对每一帧图像
寻找最符合的像素区域来进行追踪。而另一方面，由于移动是连续的，为提高追
踪精度，大方框规定了下一帧目标点出现的范围。对追踪参数进行适当调试，可
得出图示的轨迹。
实际操作中，应将选择较为靠近直杆底座的点作为B 点，如此便不需要每次
都移动B 点，而可将之大致固定在一个位置，便于计算。在6.1.1 的分析中将体
现，ACB即为太阳高度角。
6. 模型的建立与求解
问题（一）的求解
6.1.1 模型的分析
为建立合适的数学模型，首先须对垂直直杆在太阳照耀下成影情况进行分析。
可构建示意图如下：
图 6.1.1.1 竖直直杆成影示意图
如图 6.1.1.1 所示，平面XOY 代表大地，此处可将大地看作是平坦的。OA为
竖直直杆，其中O点为直杆和大地的接触点，OB 为竖直杆OA的影子。黄色有向
线代表了由太阳射来的光线。由此构成了竖直直杆成影的示意图。由于直杆长度
OA已知为 3 米，故只需知道ABO的值，即可求出影子OB的长度。其中，ABO
为太阳高度角。
参考有关太阳高度角的资料[2]，可以得到太阳高度角的求解公式如下：
  arcsin sin sin  cos cos cos t 
其 中 ，  为太阳高度角，  为太阳赤纬，  的计算方式为[3] ：
 
360
=23.45 sin 81
365
 d
 
   
 
，d 代表天数，当日期为 1 月 1 日时，d=1。 为竖
直杆所在地的纬度，t为时角。同时，由于AO  BO，故
2 2
sin
y
L
L L
 

， , y L L
分别为直杆长度和影长。
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
8
通常，时角可由本地时换算得到。然而因本地时常常因地球自转及人为调整
影响，而与该地真实的本地太阳时不相同，故本文采用精度更高的调整后计算时
角的公式。参考相关文献可知[3]，=15  12 ls t  t  ，ls t 为本地太阳时间，
60 ls lt t t

  。
lt t 为本地时。本地时的计算方法为0 15 lt
l
t  t  ，其中，l为目标所在位置经度，当
l 为负时代表所在位置为西经，为正时即东经。0 t 为子午线当地标准时间。E 为
时差（EoT）。值得注意的是，此处的时差并非生活中使用的区时时差，而是定义
为：E  9.87sin 2B  7.53cos B 1.5sin B，  
360
81
365
B  d  。由此，即可得出
精度较高的时角数值，其计算方式为：
     
  0
9.87sin 2 7.53cos 1.5sin 360
15 12 , 81
15 60 365
l B B B
t t B d
   
       
 
经过以上分析，可发现，当目标经纬度、日期、时间都已知时，可以求出相
应的太阳高度角。又因杆长已知，故可求出相应时刻的影长。因题目要求绘制出
2015 年10 月22 日北京时间9:00-15:00（即子午线标准时间1:00-7:00）之间天安
门广场（北纬39 度54 分26 秒,东经116 度23 分29 秒）3 米高的直杆的太阳影
子长度的变化曲线，故可写出太阳此地此日太阳影子长度与时间的关系式，从而
在给定的定义域中，画出此直杆影子长度的变化曲线。
6.1.2 模型的建立
根据 6.1.1 中的分析，可建立有关影长y L 与子午线当地标准时间0 t 的计算求
解模型如下：
其中 y L 为直杆的影长， L 是直杆长度，为3 米。 为太阳赤纬，计算方式
为：  
360
=23.45 sin 81
365
 d
 
   
 
，d 代表天数，当日期为 1 月 1 日时，d=1。
为竖直杆所在地的纬度，竖直杆在北半球时， 为正，否则为负。t 为时角，根
据6.1.1 中的分析，可以得出：
     
  0
9.87sin 2 7.53cos 1.5sin 360
15 12 , 81
15 60 365
l B B B
t t B d
   
       
 
 
 
2
0
1
1
sin sin cos cos cos
1,7
y L L
t
t
   
 
   
    

本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
9
0 t 为子午线当地标准时间，l 为目标所在位置经度，当竖直杆经度为东经
时，l 为正，否则为负。
可以看出，由于日期、杆长、经纬度皆为已知，所建立的模型实际上为y L
关于t0，定义域为1,7，解析式为
 2
1
1
sin sin cos cos cos
y L L
    t
 
   
    
的
函数关系式。
综上所述，第一问的模型为：
 
 
 
     
 
2 0
0
1
1 , 1,7
sin sin cos cos cos
360
=23.45 sin 81
365
9.87sin 2 7.53cos 1.5sin
15 12
15 60
360
81
365
y L L t
t
d
l B B B
t t
B d
   

 
    
    
  
    
  
    
      
  

  

6.1.3 模型的求解
采用 MATLAB 编程求解，可以得出直杆影子长度的变化曲线如下：
图 6.1.3.1 直杆的太阳影子长度变化曲线
图 6.1.3.1 表示了2015 年10 月22 日北京时间9:00-15:00（即子午线标准时
间1:00-7:00）之间天安门广场（北纬39 度54 分26 秒,东经116 度23 分29 秒）
3 米高的直杆的太阳影子长度的变化曲线，横轴北京时间（小时），纵轴为直杆影
本资料由：数模之家（春哥带你学建模）收集并分享，需要更多数学建模资料添加微信公众号：shumohome
数模之家
10
长（米）。
需要注意的是，由于本题中计算太阳赤纬、时角误差修正的公式并非完全精
确，因此本模型仍然存在一定的系统误差。在本问解答结束后，将会对本模型的
误差进行分析。
6.1.4 结果分析
根据本文的计算结果，直杆的太阳影子长度变化曲线如图6.1.3.1 所示，其在
11 点58 分30 秒时，直杆影长达到最小值3.8410 米；影长在9 点到11 点58 分
30 秒时随时间变短，此后随时间增长。
决定影子长度的除了当天的时间之外，还有杆长、日期、经度、纬度这四个
因素。由于杆长与影子长度呈正比例关系，故本文不多加讨论，而着重在日期、
经度、纬度三个方面给出讨论。为方便表现影子长度关于各个因素的关系，须控
制变量ÿ