2015数学建模A题太阳影子定位论文.doc

太阳影子定位摘 要本文主要讨论分析并建立太阳影子定位的数学模型，利用所建模型及所给数据确定事件发生的位置，运用多种数学方法研究物体影子长度与当地经纬度、时间及日期之间的联系，具有一定的实际意义。主要利用MATLAB对数据进行分析和处理。针对问题一：利用日期求出当天太阳直射点的纬度，然后利用太阳直射点的纬度，当地的纬度及时角求出太阳高度角，由三角函数关系得直杆影子长度，利用MATLAB对所求数据曲线拟合构建影子长度变化模型，分析影子长度关于各个参数的变化规律。利用所求模型求出问题一中直杆的影子长度及其变化曲线见图5.1.3.3。 针对问题二：建立时间与影子长度之间的模型，利用MATLAB曲线拟合并对求出的曲线函数分析，由时间差得出当地的经度。方法一根据相似三角形求出直杆高度，通过三角函数关系以及求太阳直射点纬度模型得直杆所在地的经纬度或。方法二利用太阳方位角和太阳高度角模型求出直杆所在地的经纬度。 针对问题三：利用相似三角形求出直杆的长度，通过曲线拟合求出正午最短的影子长度，根据三角函数关系求出正午太阳高度角和太阳高度角的方程计算当地纬度和太阳直射点的纬度，再利用问题二构建的模型求经度，最后求出附件二直杆所在地的经纬度为（北纬），日期为5月24日，附件三直杆所在地的经纬度为（南纬，东经），日期为3月24日。 针对问题四：利用MATLAB对视频进行处理，求出任意时刻影子长度，通过时间差求出视频拍摄位置的经度。利用问题一求出太阳直射点的纬度和太阳高度角随时间的变化，继而求出直杆所在地可能的经纬度或。通过曲线拟合出当地正午时最短的影子长度，利用三角函数求出正午太阳高度角和求太阳高度角的方程，计算当地纬度、太阳直射点的经纬度（南纬，东经），日期为3月24日。关键词：太阳高度角 太阳方位角 太阳直射点 曲线拟合 相似三角形 MATLAB一 问题重述如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。请分析题目，试建立数学模型讨论下列问题：1.建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用你们建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化曲线。2.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将你们的模型应用于附件1的影子顶点坐标数据，给出若干个可能的地点。3.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。将你们的模型分别应用于附件2和附件3的影子顶点坐标数据，给出若干个可能的地点与日期。4.附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高度为2米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。如果拍摄日期未知，你能否根据视频确定出拍摄地点与日期？二 问题分析问题一中，由地理知识可知直杆高度、影子长度与太阳高度角存在三角函数关系，因此要求出太阳高度角随时刻的变化就能求出影子长度的变化。而太阳高度角的变化与太阳直射纬度、当地纬度和时角有关，因此要先求出该日期的太阳直射纬度。依次求解建立模型，并利用MATLAB进行数据拟合求出直杆影子变化曲线。最后利用建立的模型求出直杆的影子长度及变化曲线。问题二中，法一利用附件1所给的影子顶点坐标位置以及北京时间和日期，建立影子长度L与时间M之间的数学模型，运用MATLAB对影子长度和北京时间进行曲线拟合并对拟合出来的曲线进行求导，得出当地正午太阳高度时直杆的最短影子长度及此时的北京时间。利用地理知识根据北京时间与当地正午十二点之间的时间差值得出当地的经度。接下来求纬度，利用相似三角形求出直杆的高度为2.76米，利用影长与直杆之间的三角函数关系求出正午太阳高度角H，根据问题一中的模型求太阳直射纬度，利用正午太阳高度角计算公式计算当地纬度,从而得出可能的经纬度位置。法二中利用太阳高度角和太阳方位角避开直杆高度，并利用问题一求太阳高度角的模型求出所在地的纬度。问题三中利用问题二中求经度的模型来对附件2、3中所给的影子长度和时间对所在位置的经度进行求解，利用相似三角形求出直杆的高度。将附件中的数据使用MATLAB进行曲线拟合得出当地正午时最短的影子长度，根据三角函数关系求出某一时刻的太阳高度角，利用问题一中建立的太阳高度角模型和正午太阳高度角计算公式求出当地的纬度和太阳直射点的纬度。利用问题一求太阳直射点的纬度模型求当天的日期。问题四中，第一小问，首先利用视频中所给的日期，利用问题一建立的模型求出太阳直射点的纬度，再利用MATLAB对所给视频进行灰度处理，求出影长，利用三角函数关系，将影长数据代入求出太阳高度角随时间的变化，最后利用问题一中建立的太阳高度角模型的公式求得当地的纬度。第二小问，用MATLAB进行曲线拟合得出当地正午时最短的影子长度，根据三角函数关系求出某一时刻的太阳高度角，利用问题一中建立的太阳高度角模型和正午太阳高度角计算公式求出当地的纬度和太阳直射点的纬度，并通过问题一求太阳直射点的纬度模型求当天的日期。另外本文以北京时间为标准，利用当地与北京的时间差确定视频拍摄的经度。三、模型假设 1. 假设模型中一天中太阳直射纬度的变化忽略不计；2. 假设不考虑固定直杆的海拔高度；3. 假设地球自转时全球各地线速度一致；四、符号说明 太阳高度角； 正午太阳高度角； 直杆杆长； 直杆影子的长度； 地理纬度； 时角； 太阳直射点的纬度； 太阳方位角。五、模型的建立与求解5.1 问题一：模型的建立与求解5.1.1太阳直射点纬度1如图5.1.1.1,为地球,为太阳,为太阳直射点纬度,;是黄经度(太阳在黄道上由春分点自西向东运行到点所转过的角度,即);是黄道平面和天赤道平面的交线(棱),在黄道面内,分别是在天赤道面的射影,它们的垂足分别是是黄赤交角(二面角的平面角),。在中,;在中，在中，所以，得到太阳直射点纬度的计算公式：图5.1.1.1地球公转的周期是一个回归年（365.2422），现行公历的历年是历日的整数倍，它和回归年并不精准相等；另外由于复杂的历史演变过程，以及一些人为的原因，上半年和下半年、冬半年和夏半年以及各季、各月之间的天数也并不完全相同，因此用日期来推算太阳直射点纬度(天文上称赤纬,一年中变动在范围内)时，需按季节分段时行计算，以确保推算的相对精确。因此，把全年的日期分成夏半年、冬半年以及自冬至日到春分日三个阶段来进行推导太阳直射点纬度。 .夏半年：从春分日（3月21日前后）到秋分日（9月23日前后）约186天，在此时段内视太阳在黄道上运转。设从春分日开始，视太阳运行了d天，则d天运行了经度：把上式带入中得： ；.冬半年： 自秋分日到冬至日（12月22日）前后总计90天，在此阶段上运行了，设从春分开始，视太阳运行了d天，则d天运行了经度：把上式带入中得： ；.自冬至日到次年春分日总计89天，在此阶段内也假设运转了。设从春分日开始，假设太阳运转了d天，则d天运行了经度：把上式带入中得： ；所以： 5.1.2太阳高度角h2我们知道中午太阳高度角最高;在北半球夏季比冬季的太阳高度角要高;低纬度地区要比中、高纬度地区的太阳高度角高。说明太阳高度角的变化是随时间、地理纬度和太阳直射点而变化的。太阳高度角计算公式是一个多元函数方程。首先研究真太阳时(真太阳时=地方平均太阳时+时差)正午的太阳高度角计算,见图5.1.2.1,AGFED表示以M为测站中心的南北天子午圈,ME为赤道截面,F为天顶。图5.1.2.1在真太阳时正午时,太阳正位于当地子午面上(即), ，由此式可见,在固定的测站上,为常数,真太阳时正午的太阳高度角只决定于太阳倾角的变化,在一年内真太阳时正午的太阳高度角变化在之内。计算任意时刻的太阳高度角,见图5.1.2.2,因为是任意时刻,太阳不恰好处在当地子午面上,O为太阳移动的位置,AGFED为天子午圈,BE为赤道截面,BF为天顶距,LG为赤道天顶距。赤道截面与天顶的交角为纬度角,故有FE,则GF。太阳与赤道截面的交角为太阳倾角则OL=,故OG=。太阳与地平面的交角为太阳高度角,则OB=h, FO=90-h。在GOF球面三角形中,根据球面三角公式：则有：将此式简化得出： 是太阳位置与当地子午面的偏角,即时角(所谓时角,是指太阳所在的时圈与通过南点的时圈构成的夹角,单位为度。自天球北极看,顺时针方向为正,逆时针方向为负。时角表示太阳的方位,因为天球在一天24内旋转360,所以每小时旋转15，当地时间12点时的时角为零）。图5.1.2.25.1.3求影子长度变化如图5.1.2.2，利用三角函数，根据即可求出指定时刻直杆的影子长度。图5.1.2.2影子长度与太阳高度角的关系由上述求解过程分析得出影子长度的变化规律：1) 与经纬度有关：就某一天来看，太阳直射纬度所在的位置，日影最短，为一圆点。在直射纬度以北、以南的地区，正午日影随着正午太阳高度缩小而逐渐变长。2) 与时间有关：就某个地点来看，一年中正午太阳高度增大时，日影逐渐缩短；正午太阳高度达最大时，日影最短；正午太阳高度减小时，日影逐渐增长；正午太阳高度达最小时，日影最长。例如：6月22日，太阳直射北回归线，北回归线及其以北各地的正午太阳高度达到全年最大，其日影也达到全年最短； 6月22日12月22日，在太阳直射点向南移动过程中，北回归线及其以北各地的正午太阳高度逐渐减小，那么其日影逐渐增长； 12月22日，太阳直射南回归线，北回归线及其以北各地的正午太阳高度达到全年最小，其日影也达到全年最长。 12月22日6月22日，在太阳直射点向北移动过程中，北回归线及其以北各地的正午太阳高度逐渐增大，那么其日影逐渐缩短。5.1.3利用模型求解通过上述模型来对2015年10月22日北京时间9:00-15:00之间天安门广场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长度的变化进行求解。每过15分钟计算一下本时刻的太阳高度角，数据见表5.1.3.1。时刻9:009:3010:0010:3011:0011:3012:00太阳高度角24.111828.257531.879534.867537.108738.501638.9476此时影子长度（米）6.70295.58154.82354.30563.96553.77133.7080时刻12:3013:0013:3014:0014：3015:00太阳高度角38.501637.108734.867531.879528.257524.1118此时影子长度（米）3.77133.96554.30564.82355.58156.7029表5.1.3.1太阳高度角随时刻的变化利用MATLAB对表5.1.3.1中的数据进行拟合求出直杆的影子长度的曲线变化见图5.1.3.3，详细代码见附录1。图5.1.3.3时刻与直杆影长之间的关系5.2问题二：模型的建立与求解5.2.1所在地经度的求解本文首先利用附件一中所给的影子顶点坐标求出影子长度，建立横坐标为时间M与纵坐标为影子长度L之间的函数关系，利用MATLAB整合出曲线方程以及曲线拟合图像，如图5.2.1.1：图5.2.1.1曲线方程为：正午太阳高度最大，直杆影子长度最短，由于图像横坐标为北京时间，因此曲线最低点对应的北京时间为当地正午十二时。对所求曲线方程进行求导，求出最低点的北京时间和最短的影子长度： 令,得出M=12.5991，。利用地理知识经度每相差在时间上相差小时得出当地的经度与北京的经度相差，由求导得出的结果可知北京时间比当地时间快，因此该地应在北京的西面，所以得出当地经度为：5.2.2所在地纬度的求解法一：由于直杆高度和地球最外层距离地面的距离都是固定的，利用直杆影长的变化和经过相同时间地球自转距离的差值，进而用相似三角形求出直杆的高度，原理如图5.2.2.1。图5.2.2.1相似三角形的比值公式如下：为确定时间内影子变化长度；为确定时间内地球自转距离的差值4；为地球最外层距离地球表面的距离5 。根据比值求出固定直杆的高度。利用直角三角形的角度关系求出当地正午太阳高度角，代入正午最短的直杆影子长度，由最终求得正午太阳高度角。通过问题一中求太阳高度角的计算公式，将正午的时角0代入，利用三角函数和角公式可推导出正午太阳高度角计算公式，纬差即为当地纬度与太阳直射点纬度之差，从而求出纬度位置，利用题目一所建立的日期与太阳直射点纬度之间的关系模型来求出4月18日太阳直射点纬度，最终求得当地纬度:最后得出直杆两个可能的位置: ，法二：利用太阳高度角求出所在地的纬度3 太阳方位角即太阳所在的方位，指太阳光线在地平面上的投影与当地子午线的夹角，可近似地看作是竖立在地面上的直线在阳光下的阴影与正南方的夹角。方位角以目标物正北方向为零，顺时针方向逐渐变大，其取值范围是0-360。因此太阳方位角一般是以目标物的北方向为起始方向，以太阳光的入射方向为终止方向，按顺时针方向所测量的角度，如图5.2.2.1：图5.2.2.1根据“天球理论”，可以导出太阳方位角与太阳直射点纬度、太阳高度角以及观测时刻的关系如下:由图5.2.2.1可以得出：由上述方程组可知：又有，将附件一中的北京时间转化为当地的时角代入方程，并利用题目一所建立的日期与太阳直射点纬度之间的关系模型来求出4月18日太阳直射点的纬度，最后由MATLAB解得： ，所以求出直杆所在位置为。5.3问题三：模型的建立与求解5.3.1所在地经度的求解本文首先利用附件2和附件3中所给的影子顶点坐标求出影子长度，建立横坐标为时间M与纵坐标为影子长度L之间的函数关系，利用MATLAB分别拟合出附件2、附件3数据的曲线方程以及图像，如图5.3.1.1和5.3.1.2，代码见附录2：图5.3.1.1图5.3.1.2附件2拟合出来的曲线方程为：附件3拟合出来的曲线方程为：正午太阳高度最大，直杆影子长度最短，由于图像横坐标为北京时间，因此曲线最低点对应的北京时间为当地正午十二时。对上述所求曲线方程分别进行求导，求出最低点的北京时间和最短的影子长度： 令,得出，。令,得出，。利用地理知识经度每相差在时间上相差小时得出当地的经度与北京的经度分别相差和，由求导得出的结果可知北京时间比当地时间快，因此两地都应在北京的西面，所以得出两地经度分别为：5.3.2所在地纬度与日期的求解对附件2，利用问题二的相似三角形模型求出固定直杆的长度，所以。利用直杆的高度与当地最短的影子长度的比值得出正午太阳高度角，求得。取附件二中某一时刻的影子顶点的坐标并求出影子长度，运用三角函数公式求出该时刻的太阳高度角，联立方程组：分两种情况讨论：若当地纬度与太阳直射点的纬度在同一半球，纬差为；若当地纬度与太阳直射点的纬度在不同半球，纬差为。用MATLAB求解得当地纬度，根据问题一建立的求太阳直射点的纬度的模型将代入三个阶段，筛选出符合条件的日期为5月24日。当地的经纬度为（北纬）。同理，对附件三进行求解得，太阳直射点的纬度为，当地经纬度为（南纬，东经），日期为3月24日。5.4问题四：模型的建立与求解5.4.1模型假设1) 假设直杆底座厚度忽略不计；2) 假设忽略风速对直杆高度、影子长度的影响；3) 假设拍摄过程中摄像头位置及角度的变化忽略不计。5.4.2第一小问模型的建立与求解1) 所在地纬度的求解本文首先利用视频中所给的日期，利用问题一中建立的求太阳直射点纬度的模型进行求解，根据公式求出太阳直射点的纬度。对用MATLAB进行灰度处理的21张图片进行影子长度的测量，测量数据见表5.4.2.1，利用影子的长度变化数据根据三角函数求出对应时刻的太阳高度角，根据问题一建立的求太阳高度角的数学公式，利用MATLAB6 对数据进行求解得出当地的纬度。时刻8:548:568:589:009:029:04影子长度2.39492.37202.33642.31442.37372.2423时刻9:069:089:109:129:149:16影子长度2.23302.18812.15672.14322.10332.0762时刻9:189:209:229:249:269:28影子长度2.05842.02711.99151.97791.94151.9101时刻9:309:329:34影子长度1.8811.87411.8254表5.4.2.12) 所在地经度的求解我们以北京时间为标准，求出北京与当地的时间差，从而根据两地的经度差求出当地的经度。假设将此2m的直杆放在北京，利用达到相同影子长度的时刻求时间差，选取某一时刻的影子长度，求出该时刻的太阳高度角，再利用视频中已知的日期求出当天太阳直射点的纬度，利用问题一建立的模型求出时角，由两地的时角差计算出两地的时间差，继而根据经度每隔时间相差一小时可求出分别位于北京东面和西面的两地经度为。因此，得出视频大体拍摄位置为和。5.4.3问题4.2模型的建立与求解利用MATLAB对上述影长与时间进行曲线拟合，图像如图5.4.2.2得出的表达式如下：；图5.4.2.2当地正午十二时时，直杆影子的长度最短，根据曲线求得最短影长L=0.3450。利用三角函数关系，求出影子最短时的正午太阳高度角。再联立方程组：筛选出合适