2015高教社杯全国大学生数学建模竞赛题目
（请先阅读“全国大学生数学建模竞赛论文格式规范”）
如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。 1.建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用你们建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。 2.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将你们的模型应用于附件1的影子顶点坐标数据，给出若干个可能的地点。
3. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地
点和日期。将你们的模型分别应用于附件2和附件3的影子顶点坐标数据，给出若干个可能的地点与日期。
4．附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高度为2米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。
如果拍摄日期未知，你能否根据视频确定出拍摄地点与日期？
太阳影子定位
摘要
本文通过分析物体的太阳影子变化，利用太阳影子定位技术建立确定视频拍摄的地点和日期的模型。
针对问题一，首先通过分析知影子长度的变化主要影响参数为：当地的经度、纬度、时刻、直杆长度、季节（日期)等，引入地理学参数：太阳赤纬、时角及太阳高度角，建立一个能够刻画影子长度变化和各个参数
转子动平衡
间关系的模型：；其次以实例对模型进行检验，在误差可允许的范围内，认为模型正确；进而对模型采用控制变量法分析影子长度关于各个参数的变化规律；然后求解出满足条件影子长度12时15分是最短，大约3.674米（表3）。影子长度的变化曲线（图5），9时至12时15分影子长度呈现下降趋势，12时15分之15时影子长度呈现上升趋势；最后考虑太阳照射中发生折射现象的推广。
针对问题二，
关键词
定厚机
1、问题重述：
如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。
1.建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用你们建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。
2.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将你们的模型应用于附件1的影子顶点坐标数据，给出若干个可能的地点。
3. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。将你们的模型分别应用于附件2和附件3的影子顶点坐标数据，给出若干个可能的地点与日期。
4．附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高度为2米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。
如果拍摄日期未知，你能否根据视频确定出拍摄地点与日期？
2、问题分析：
针对问题一，首先通过分析知影子长度的变化主要影响参数为：当地的经度、纬度、时刻、直杆长度、季节（日期)等，引入地理学参数：太阳赤纬、时角及太阳高度角，建立一个能够刻画影子长度变化和各个参数
间关系的模型；其次以实例对模型进行检验，在误差可允许的范围内，认为模型正确；进而根据所建模型分析影子长度关于各个参数的变化规律；然后做出影子长度的变化曲线；最后考虑太阳照射中发生折射现象的推广。
针对问题二，
三、模型假设：
1、求解此问题时忽略地球的自转
2、不考虑太阳光线在穿过大气层时的折射、太阳的视面角、高山阻挡、海拔高度等因素。
3、认为照射到地球上的太阳光可以看成是平行光线，地球上某地的水平地面是地球球面上过该地的切面。
四、符号说明：
：太阳赤纬
：太阳时角
：表示某地的地理纬度
燃油调压阀：表示某地的地理经度
：太阳高度角
五、问题一的模型建立与求解
5.1影响影子长度参数的确定
1、太阳赤纬
太阳赤纬是地球赤道平面与太阳和地球中心的连线之间的夹角且以年为周期，在周年运动中任何时刻的赤纬值都是严格已知的，可用下式计算：
式中称为日角，即，这里,式中为积日，就是日期在年内的顺序号（例如1月1日其积日为1，平年1212月31日的积日为365，闰年则为366等）。
（式中表示取整数部分）
故如果已知某日期的年、月、日，代入式（1），即可求得此日期的太阳赤纬值。
2、时角
时角表示一天体是否通过了当地的子午圈，其数值表示该天体与当地子午圈的角距离，并借用时间的单位以小时来计量，其中当地时间12点时的时角为零，令上午的时角为正，下午为负。
某地的时角计算式如下：
故若给定某地时刻的值，代入式（2），即可求出此时的时角值
3、太阳高度角
太阳高度角,是指太阳光线与地平面的夹角。应用球面三角形余弦公式，结合图形,可以推出任意时刻太阳高度角的计算公式为：
进一步可以得到：
式中，表示当地的地理纬度,表示太阳赤纬，表示太阳时角、表示太阳高度角。、的取值为北正南负。
故若已知某地的地理纬度、太阳赤纬、太阳时角，代入式（4），即可求得太阳高度角。
4、直杆的长度
直杆的影子始终在物体背着光源的一面，光从物体顶端照射到地面形成影子，直杆的长度直接影响着影子长度的变化。
5.2模型一的建立
首先分析影子长度变化直接受直杆长度和太阳高度角的影响，其关系式为：
其中表示影子的长度。
然后以太阳高度角与太阳赤纬、时角、地理纬度的关系为基础，即式（4），引进参数当地的经度，构建影子长度变化模型如下：
式中、、、、分别代表当地的经度、纬度、太阳赤纬、太阳高度角和时刻（北京时间）。编写程序计算时注意到，当时，反余弦函数值才存在，应采用判断，当其大于时则进行下一个的计算。
5.3模型一的检验
以实际某地情况为例，已知当地的经度为，纬度为，
5.4影子长度关于各个参数的变化规律
根据影子长度变化模型，以影子长度为因变量，依次选择各个参数为自变量，其余参数看为固定值，进而描述出影子长度关于各个参数变化规律。
1、影子长度关于直杆长度的变化规律
以直杆长度为自变量，以影子长度为因变量，太阳高度角为固定值，取直杆长度为米，间隔为0.5米，根据公式（5），得到与之间的关系，画出关系图像如图1：
由图1我们可以很直观的的看出，在其他参数不变时，直杆长度越长，影子长度越长，且两者的比值是不变的。
2、影子长度的日变化规律
以每日的时间作为自变量，影子长度为因变量，其它参数为固定值，取为2015年10月22日北京时间8:00-16:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆，根据公式（6），得到影子长度的日变化规律。
由于地球是个固体球,且自西向东自转,因此太阳高度也呈现出与之对应的日变化规律:地球上观察太阳为东升西落,早晨太阳从东方地平线上升起,晨昏线上太阳高度为，随着太阳的逐渐升高,太阳高度是逐渐增大的，影子的长度逐渐减小。当某地经线正对太阳光时,地方时为正午12点，即为北京时间12时15分,此时太阳高度角达到一天中的最大，影子长度最小。之后太阳逐渐西落,太阳高度也慢慢变小，影子逐渐变长,到西方地平线落下时,没有影子。
3、影子长度随纬度变化规律
以纬度作为自变量，影子长度为因变量。其它参数为固定值，由于正午太阳高度角最能反映太阳辐射的强弱变化,故取2015年10月22日北京时间12时东经116度23分29秒3米高的直杆，根据公式（6），得到影子长度的随纬度的变化规律。
压电陶瓷驱动器正午太阳高度由直射点向南北两侧递减，故影子长度由直射点向南北两侧递增。夏至:，影子长度从北回归线向南北两侧递增;同样冬至:，影子长度从南回归线向南北两侧递增;春秋分:，影子长度从赤道向南北两侧递增,且离直射点距离越近,与直射点纬度差越小,影子长度就越小。
4、影子长度随季节的变化规律
以季节作为自变量，影子长度为因变量。由于正午太阳高度角最能反映太阳辐射的强弱变化,故取2015年北京时间12时天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆，根据公式（6），得到影子长度的随季节的变化规律。
北回归线及以北地区,6月22日正午太阳高度达一年中最大值，此时影子长度最小,12月22
日达一年中最大值;南回归线以南地区,12月22日达一年中最小值,6月22日达一年中最大值;南北回归线之间,一年中因有两次太阳直射机会,赤道至北回归线之间12月22日的正午影子长度达一年中最大值,而赤道至南回归线之间6月22日达一年中最大值。
5.5问题一的求解
求解步骤：
1）将2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆符号化为，，；
2）计算太阳赤纬：
已知年份为2015年，积日为295，代入公式得：
得到太阳的赤纬角为
3）计算太阳时角：
需将当地时间转化为北京时间，对应北京时间，将代，得到2015年10月22日北京时间天安门广场时角值（如表1），其中。
表1：2015年10月冷轧酸洗22日北京时间天安门广场时角统计表
时刻 | 9:00 | 9:30 | 10:00 | 10:30 | 11:00 | 11:12 | 11:24 |
时角/度 | -48.61 | -41.11 | -33.61 | -26.11 | -18.61 | -15.61 | -12.61 |
时刻 | 11:36 | 11:48 | 12:00 | 12:12 | 12:15 | 12:24 | 12:36 |
亨润成型机炮筒图时角/度 | -9.61 | -6.61 | -3.61 | -0.61 | 0.14 | 2.39 | 5.39 |
时刻 | 12:48 | 13:00 | 13:30 | 14:00 | 14:30 | 15:00 | |
时角/度 | 8.39 | 11.39 | 18.89 | 26.39 | 33.89 | 41.39 | |
| | | | | | | |
4）计算太阳高度角：
将所得到的、、分别代入公式（4），得到2015年10月22日北京时间天安门广场太阳高度角值（如表2）
