摘要
本文通过对太阳光与地球影子的分析，建立了利用太阳影子的变化得出其定为和日期的数学模型。
问题一需要对确切时间地点的直杆的影长随时间的变化进行描述。由于太阳与地球的距离和杆长的长度数量级别相差太大，用长度与距离的比例关系进行计算会导致误差太大，所以采用太阳高度角进行求解。根据对实际情况的分析，可推导出太阳高度角与赤纬角、观测地点的纬度，一天之中的具体时间的函数关系。在给定地点日期的情况下，赤纬角、观测地点的纬度随之确定，可得到太阳高度角随时间的变化函数。直杆长的长度与直杆影长呈现直角三角形关系，二斜边与直杆影长的夹角巧好为太阳高度角。得到太阳高度角在一天之中随时间变化的函数关系，通过三角公式的变换，也就得到了一天之中直杆影长随时间变化的函数曲线。
针对问题二，先画出直杆影长与时间的散点图，用二次函数进行拟合，得出其最短影长所对应的北京时间。正午的影长最短，所以其最短影长所对应的时间即为观测点的正午时间。用观测点的正午时间与北京的正午时间作计算比较，即可得出观测点的经度值。同样，利用附件所给数据，画出时间与影长的散点图，利用问题一所得的模型对其进行拟合，得出可能的最优解。并求出每个最优解下的直杆地点。通过我们的计算拟合，得到其最优解为经度为110.7750，纬度为
18.50，通过地图的查询，其具体位置为海南省陵水县。
问题三与问题二类似，只是在给出数据时减少了具体的日期。在求解过程中依旧是先对其数据做二次函数的拟合，得出其最短影长所对应的北京时间，通过相应的模型求解出其经度。再利用问题一所给的模型，对其直杆影长随时间变化的散点图进行拟合，得出其可能的最优解及其情况下所对应的地点与日期。附件二可能的地点有两个，第一个为南纬16.470，东经880，时期为5月14号，具体位置为印度洋。第二个为北纬33.790，东经880，时期为4月15号，具体位置为西藏自治区。附件三可能的地点有两个，一个是南纬33.150，东经109.1250，日期为4月21日，地点在印度洋。第二个北纬46.10，东经109.1250，日期为6月11日，地点在蒙古国。
针对问题四，可利用比例从视频中得出相关的数据，再利用问题三建立的模型进行求解。
关键词：太阳高度角；时角；最小二乘法
一、问题重述
如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。
1.建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用你们建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。
2.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将你们的模型应用于附件1的影子顶点坐标数据，给出若干个可能的地点。
3. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。将你们的模型分别应用于附件2和附件3的影子顶点坐标数据，给出若干个可能的地点与日期。
4．附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高度为2米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。
如果拍摄日期未知，你能否根据视频确定出拍摄地点与日期？
二、问题分析
问题一需要在已知时间、地点的情况下，对直杆的影长随时间的变化的曲线进行描述。由于太阳与地球的距离和杆长的长度数量级别相差太大，用长度与距离的比例关系进行计算会导致误差太大。在计算球面弧长，球面高度等一般习惯用角度进行描述，且在计算影长时，用角度进行计算的准确率较高。在计算过程中，将太阳影子顶点与杆长顶点进行连线，此线段与水平地面的夹角可求，利用这个角度即可与杆长求出影长。在地理学中，这个角度被称为太阳高度角，它在一天内，随时间的变化而变化。所以在知道杆长的情况下，只要知道太阳高度角的变化情况即可对影长的变化进行描述。太阳高度角与赤纬角（太阳直射纬度）、地点纬度、时间三者有关。在确定地点为北京，时间为10月22日的情况下，太阳赤纬角只与时间有关，且为可测的函数关系。得到函数关系之后，即可对影长的变化曲线进行描述。
问题二给出各个时间下太阳影子的顶点坐标数据，并且给出了具体的日期。但是杆长，经纬度为未知数。在此条件下，需要建立模型求解出数据对应的直杆可能的地点。首先，可以直杆在地面的位置为原点，计算出各个时间下的杆影长。
再利用二次函数对时间与杆影长的关系进行拟合，从而求出最短杆影长所对应的北京时间。由于杆影长在正午时间最短，所以最短影长所对应的时间即为此地的正午时间，时间形式为北京时间。利用这个时间与北京正午时间的比较，即可得出此地的经度。用理想的杆影长与时间的曲线对实际的杆影长与时间的变化曲线进行拟合，即可得到可能的最优解，得出各个情况下的经纬度。
问题三同样是在告知各个时间下太阳影子的顶点坐标数据的情况下，要求建立模型求解出数据对应的直杆的位置与拍摄的日期。同问题二，利用二次函数的拟合可求出最短影长的所对应的北京时间。再用此时间与北京的正午时间作比较计算，得出此地的经度。最后，同样利用最小二乘法对实际的杆影长与时间的变化曲线进行拟合，得出其可能的最优解，并得出各个最优解情况下的直杆所对应的时间与地点。
三、模型假设
1.假设地球不动，太阳绕地球运动；
2.同时由于地球非常大，在某一点的周围可以看出是一个平面； 3.且在同一天里面忽略太阳的公转所带来的影响。 4.太阳光为平行光
四、符号说明
l：直杆影长； H：直杆长度；
h：太阳光与地面的夹角，即太阳高度角；
?：观测地点的纬度； ?：太阳赤纬角；
?：t时刻的太阳光与正午的太阳光的夹角，即时角；
t1：此地的正午时间所对应的北京时间。
五、模型建立与求解
5.1问题一的模型建立与求解
问题一需要在已知日期与地点的情况下，对太阳影子随时间变化的曲线进行描述。但是由于太阳离地球的距离、地球的半径、直杆的长度三者的数量级别相差太大，且地球在不断的自传与公转，利用长度比例进行计算时，会产生非常大的误差。所以在计算时，可以避开长度的计算，利用角度对直杆的影长进行计算。如下图：
图一，阳光直射图
图中，H为杆长，BC为太阳光，AC为地面,h为太阳光与地面所成的夹角，在图中为?ACB,l为直杆的影长。因为直杆与地面垂直，所以△ABC为直角三角形，影长的公式为：
l?H tanh式中，l为影长，H为直杆长度，h为太阳光与地面的夹角。
由上述公式可知，只要知道了太阳光与地面的夹角h，就可以知道影长与影长的变化。下面我们展开对角度h的推导与求解。 5.1.1太阳高度角
在实际过程中，地球是个圆形，在自转的同时绕太阳公转。但为了计算与建模的方便，假设地球不动，太阳绕地球运动；同时由于地球非常大，在某一点的周围可以看出是一个水平面；且在同一天里面忽略太阳的公转所带来的影响。
由此，可以通过计算得到一天之中太阳每一时刻相对于地球的位置，且可以求出每个时刻太阳光与地面的夹角h。
通过我们资料查询，上述太阳光与地面的夹角h在地理学中被称为太阳高度角。它与时间日期有关、经纬度、某一确定天数的时刻三者有关。现在通过相关的地理知识与数学知识进行推导。
为了推导的方便，需要引入两组描述太阳和地球相对位置的坐标系。 第一组给出太阳相对于固定在地球上的坐标系的位置，坐标系的一个轴指向天顶，另一个与之正交的轴指向地平线，即下图中的i, j坐标系，称为地平坐标系。
另一组坐标系也固定在地表的同一个位置，但是一个轴指向极点，即北天极，另一个与之正交的轴则指向与赤道平行的方向，即下图中的I, J坐标系，称为时角坐标系。
图二.坐标系
图中，?为两个坐标系之间的夹角,即观测的纬度。
两个坐标系之间有如下关系：
i?(sin?)I?(cos?)J
（cos?)?(sin?)J j? k?K
式中，?为两个坐标系之间的夹角,即观测的纬度。
下图给出了正午时分的太阳（在天空的最高点）在地平坐标系中的方位。
图三.地平坐标系
从而，指向太阳的单位矢量可以写成：
s?(cosz)i?(coszcosaz)j?(sinzsinaz)k
其中z称为天顶角而az称为太阳方位角。
第二个坐标系的一个轴始终指向北天极而另一轴位于赤道面上，方向指向南点。因而，指向正午太阳的单位矢量可以写成：
s?(sin?)I?(cos?cos?)J?(cos?sin?)K
图四.时角坐标系
?为太阳光与时角坐标系JOK平面的夹角，图中，?为太阳光在JOK的投影与J轴的夹角，正午时分为00。