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

1、太阳影子定位摘要本文主要讨论分析并建立太阳影子定位的数学模型，利用所建模型及所给数据确定事件发生的位置，运用多种数学方法研究物体影子长度与当地经纬度、时间及日期之间的联系，具有一定的实际意义。主要利用MATLAB 对数据进行分析和处理。针对问题一：利用日期求出当天太阳直射点的纬度，然后利用太阳直射点的纬度，当地的纬度及时角求出太阳高度角，由三角函数关系得直杆影子长度，利用 MATLAB 对所求数据曲线拟合构建影子长度变化模型，分析影子长度关于各个参数的变化规律。利用所求模型求出问题一中直杆的影子长度及其变化曲线见图。针对问题二：建立时间与影子长度之间的模型，利用MATLAB 曲线拟合并对求出的

2、曲线函数分析，由时间差得出当地的经度。方法一根据相似三角形求出直杆高度，通过三角函数关系以及求太阳直射点纬度模型得直杆所在地的经纬度北纬 0 7'46'',东经107 24'18'' 或北纬 20 44'48'',东经107 24'18'' 。方法二利用太阳方位角和太阳高度角模型求出直杆所在地的经纬度南纬 8 2'16'',东经107 24'18'' 。针对问题三：利用相似三角形求出直杆的长度，通过曲线拟合求出正午最短的影子长度，根据三角函数关系求出正

3、午太阳高度角和太阳高度角的方程计算当地纬度和太阳直射点的纬度，再利用问题二构建的模型求经度，最后求出附件二直杆所在地的经纬度为（北 > ' " ' " . 、 一 、 . 一 一 、 、一、纬 39 52'48'',东经68 16'22''） ，日期为5月 24日，附件三直杆所在地的经纬度为（南纬51 10'15''，东经105 19'25''） ，日期为3月 24日。针对问题四：利用MATLAB 对视频进行处理，求出任意时刻影子长度，通过时间差求出视频拍摄

4、位置的经度。利用问题一求出太阳直射点的纬度和太阳高度角随时间的变化,.'rr .r 一 一一z一_，一一，一, r，一'"1继 而 求 出 直 杆 所 在 地 可 能 的 经 纬 度 （北纬 38.8626 ，东经 21 2329''） 或' " 、 . . . . ._一 _., .北纬 38.3626 ，东经 121 23'29'' 。 通过曲线拟合出当地正午时最短的影子长度，利用三角函数求出正午太阳高度角和求太阳高度角的方程，计算当地纬度、太阳直射点的经纬度（南纬 51 10'15'

5、9;，东经105 19'25''） ，日期为3月 24日。关键词： 太阳高度角太阳方位角太阳直射点曲线拟合相似三角形MATLAB一 问题重述如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。请分析题目，试建立数学模型讨论下列问题：1 .建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用你们建立的模型画出2015年10月 22日北京时间9:00-15:00之间天安门广场（北纬39度 54分26秒 ,东经116度 23分 29秒） 3米高的直杆的太阳影子长度

6、的变化曲线。2 .根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将你们的模型应用于附件1 的影子顶点坐标数据，给出若干个可能的地点。3 .根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。将你们的模型分别应用于附件2 和附件 3 的影子顶点坐标数据，给出若干个可能的地点与日期。4 .附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高 度为 2 米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。如果拍摄日期未知，你能否根据视频确定出拍摄地点与日期？二 问题分析问题一中，

7、由地理知识可知直杆高度、影子长度与太阳高度角存在三角函数关系，因此要求出太阳高度角随时刻的变化就能求出影子长度的变化。而太阳高度角的变化与太阳直射纬度、当地纬度和时角有关，因此要先求出该日期的太阳直射纬度。依次求解建立模型， 并利用 MATLAB 进行数据拟合求出直杆影子变化曲线。最后利用建立的模型求出直杆的影子长度及变化曲线。问题二中，法一利用附件1 所给的影子顶点坐标位置以及北京时间和日期，建立影子长度 L 与时间 M 之间的数学模型，运用MATLAB 对影子长度和北京时间进行曲线拟合并对拟合出来的曲线进行求导，得出当地正午太阳高度时直杆的最短影子长度及此时的北京时间。利用地理知识根据北京

8、时间与当地正午十二点之间的时间差值得出当地的经度。接下来求纬度，利用相似三角形求出直杆的高度为米，利用影长与直杆之间的三角函数关系求出正午太阳高度角H,根据问题一中的模型求太阳直射纬度，利用正午太阳高度角计算公式计算当地纬度,从而得出可能的经纬度位置。法二中利用太阳高度角和太阳方位角避开直杆高度，并利用问题一求太阳高度角的模型求出所在地的纬度。问题三中利用问题二中求经度的模型来对附件2、 3 中所给的影子长度和时间对所在位置的经度进行求解，利用相似三角形求出直杆的高度。将附件中的数据使用MATLAB 进行根据三角函数关系求出某一时刻的太阳高度角，利用问题一中建立的太阳高度角模型和正午太阳高度角

9、计算公式求出当地的纬度和太阳直射点的纬度。利用问题一求太阳直射点的纬度模型求当天的日期。问题四中，第一小问，首先利用视频中所给的日期，利用问题一建立的模型求出太阳直射点的纬度，再利用 MATLAB 对所给视频进行灰度处理，求出影长，利用三角函数关系，将影长数据代入求出太阳高度角随时间的变化，最后利用问题一中建立的太阳高度角模型的公式求得当地的纬度。第二小问，用 MATLAB 进行曲线拟合得出当地正午时最短的影子长度，根据三角函数关系求出某一时刻的太阳高度角，利用问题一中建立的太阳高度角模型和正午太阳高度角计算公式求出当地的纬度和太阳直射点的纬度，并通过问题一求太阳直射点的纬度模型求当天的日期。

10、另外本文以北京时间为标准，利用当地与北京的时间差确定视频拍摄的经度。三、模型假设1. 假设模型中一天中太阳直射纬度的变化忽略不计；2. 假设不考虑固定直杆的海拔高度；3. 假设地球自转时全球各地线速度一致；四、符号说明h太阳高度角；正午太阳高度角；l直杆杆长；L直杆影子的长度；地理纬度；t时角；太阳直射点的纬度；A太阳方位角。五、模型的建立与求解问题一：模型的建立与求解太阳直射点纬度1如图 O 为地球 , S 为太阳 , 为太阳直射点纬度, SOC ; 是黄经度（太阳在黄道上由春分点 自西向东运行到S 点所转过的角度,即SOA）;O 是黄道平面和天赤道平面的交线（棱）,OD、AS OS在黄道面

11、内，OD O ,AS O ;OE、AC、OC 分别是 OD、AS OS在天赤道面的射影,它们的垂足分别是E、 C、 C ; AOE 是黄赤交角（二面角的平面角），DOE SAC 23 26 21。在 Rt SAO中,AS OSsin ;在RtOCS中，SC OS sin在 Rt ACS 中，SC AS sin 23 26 21 OS sin sin 23 26 21''''''OS sin OS sin sin 23 2623 ， 所以， sin sin sin 23 2623 ， 得到太阳直射点纬度 的计算公式:图地球公转的周期是一个回归年（）

12、，现行公历的历年是历日的整数倍，它和回归年并 不精准相等；另外由于复杂的历史演变过程，以及一些人为的原因，上半年和下半年、冬 半年和夏半年以及各季、各月之间的天数也并不完全相同，因此用日期来推算太阳直射点 纬度（天文上称赤纬，一年中变动在 23 26'范围内）时，需按季节分段时行计算，以确保推算 的相对精确。因此，把全年的日期分成夏半年、冬半年以及自冬至日到春分日三个阶段来进行推导 太阳直射点纬度 。1 .一半年:从春分日（3月21日前后）到秋分日（9月23日前后）约186天，在此时段内视太阳在黄道上运转180。设从春分日开始，视太阳运行了 d天,则d大运行了 经度：把上式带入 arc

13、sin 0.397775sin 中得：180 arcsin 0.39775 sin d d 0,1,2 186 ;1862 .冬半年：自秋分日到冬至日（12月22日）前后总计90天，在此阶段上运行了 90 ,设从春分开始，视太阳运行了 d天,则d大运行了 经度：把上式带入arcsin 0.39775sin 中得：arcsin 0.39775 sin d 186d 186,187 276 ;3 .自冬至日到次年春分日总计 89天，在此阶段内也假设运转了 90。设从春分日开始， 假设太阳运转了 d天,则d大运行了 经度：把上式带入arcsin 0.39775sin中得： 90 一 "ar

14、csin 0.39775 cos n 276 n 276,277 365 ;89所以：太阳高度角h2我们知道中午太阳高度角最高；在北半球夏季比冬季的太阳高度角要高；低纬度地区要比中、高纬度地区的太阳高度角高。说明太阳高度角的变化是随时间、地理纬度和太阳直 射点而变化的。太阳高度角计算公式是一个多元函数方程。首先研究真太阳时（真太阳时二地方平均太阳时+时差）正午的太阳高度角计算，见图表示以M为测站中心的南北天子午 圈,ME为赤道截面，F为天顶。图在真太阳时正午时，太阳正位于当地子午面上（即O）, h 90，由此式可见，在固定的测站上，为常数真太阳时正午的太阳高度角只决定于太阳倾角的变化，在一年内

15、真太 阳时正午的太阳高度角变化在 23 26'之内。计算任意时刻的太阳高度角，见图因为是任意时刻，太阳不恰好处在当地子午面上,O为 太阳移动的位置，AGFED为天子午圈，BE为赤道截面，BF为大顶距，LG为赤道大顶距。赤道 截面与天顶的交角为纬度角，故有FE ，则GF 90。太阳与赤道截面的交角为太阳倾角则OL=，故OG= 90。太阳与地平面的交角为太阳高度角，则OB=h, FO=90-h。在GOF球面三角形中根据球面三角公式:则有:将此式简化得出：t是太阳位置与当地子午面的偏角，即时角（所谓时角，是指太阳所在的时圈与通过南点 的时圈构成的夹角，单位为度。自天球北极看，顺时针方向为正，

16、逆时针方向为负。时角表示太 阳的方位用为天球在一天24h内旋转360° ,所以每小时旋转15° ,当地时间12点时的时 角为零）。图求影子长度变化如图，利用三角函数，根据tanh,即可求出指定时刻直杆的影子长度。L图影子长度与太阳高度角的关系由上述求解过程分析得出影子长度的变化规律：1）与经纬度有关：就某一天来看，太阳直射纬度所在的位置，日影最短，为一圆点。在直射纬度以北、 以南的地区，正午日影随着正午太阳高度缩小而逐渐变长。2）与时间有关：就某个地点来看，一年中正午太阳高度增大时，日影逐渐缩短；正午太阳高度达最大时，日影最短；正午太阳高度减小时，日影逐渐增长；正午太阳高度

17、达最小时，日影最长。例如：6月22日，太阳直射北回归线，北回归线及其以北各地的正午太阳高度达到全年最大，其日影也达到全年最短；6月22日12月22日，在太阳直射点向南移动过程中， 北回归线及其以北各 地的正午太阳高度逐渐减小，那么其日影逐渐增长；12月22日，太阳直射南回归线，北回归线及其以北各地的正午太阳高度达到 全年最小，其日影也达到全年最长。12月22日6月22日，在太阳直射点向北移动过程中，北回归线及其以北各 地的正午太阳高度逐渐增大，那么其日影逐渐缩短。利用模型求解通过上述模型来对2015年10月22日北京时间9:00-15:0必问天安门广场（北纬39度54分26秒，东经116度23

18、分29秒）3米高的直杆的太阳影子长度的变化进行求解。每过15分钟计算一下本时刻的太阳高度角，数据见表。时亥IJ9:009:3010:0010:3011:0011:3012:00太阳局度角此时影子长度（米）时亥IJ12:3013:0013:3014:0014: 3015:00太阳局度角此时影子长度（米）表太阳高度角随时刻的变化利用MATLAB对表中的数据进行拟合求出直杆的影子长度的曲线变化见图，详细代码见附录1。图时刻与直杆影长之间的关系问题二：模型的建立与求解所在地经度的求解本文首先利用附件一中所给的影子顶点坐标求出影子长度，建立横坐标为时间M与纵坐标为影子长度L之间的函数关系，利用MATLA

19、B整合出曲线方程以及曲线拟合图像， 如 图：图曲线方程为：正午太阳高度最大，直杆影子长度最短，由于图像横坐标为北京时间，因此曲线最低 点对应的北京时间为当地正午十二时。对所求曲线方程进行求导，求出最低点的北京时间 和最短的影子长度：令 L' 0,得出 M= , l 0.4952。利用地理知识经度每相差15在时间上相差1小时得出当地的经度与北京的经度相差8.9865 ,由求导得出的结果可知北京时间比当地时间快，因此该地应在北京的西面，所以得出当地经度为:所在地纬度的求解法一：由于直杆高度和地球最外层距离地面的距离都是固定的，利用直杆影长的变化和经过相同时间地球自转距离的差值，进而用相似三

20、角形求出直杆的高度，原理如图。图相似三角形的比值公式如下：Li L2为确定时间内影子变化长度；Pi P2为确定时间内地球自转距离的差值4；X为地球最外层距离地球表面的距离5。6根据比值求出固定直杆的高度l 0.06 3 102.7226m。465 3 60利用直角三角形的角度关系tanH L求出当地正午太阳高度角H,代入正午最短的直 L杆影子长度，由H arctan(最终求得正午太阳高度角H 79.6914。通过问题一中求太阳高度角的计算公式 sinh sin sin cos cos cost ,将正午的时角0代入，利用三角函数和角公式 cos(A B) cos A cosB sin A si

21、n B可推导出正午太阳高度角计算公式H 90纬差| ,纬差即为当地纬度与太阳直射点纬度之差，从而求出纬度位置，利用题目一所建立的日期与太阳直射点纬度之间的关系模型来求出4月18日太阳直射点纬度，最终求得当地纬度：最后得出直杆两个可能的位置一 ' "，一 ' 一 ' "，一 '北纬 0 746 ,东经107 2418 ， 北纬 20 4448 ,东经 107 2418法二：利用太阳高度角求出所在地的纬度3太阳方位角即太阳所在的方位，指太阳光线在地平面上的投影与当地子午线的夹角，可近似地看作是竖立在地面上的直线在阳光下的阴影与正南方的夹角。方位角

22、以目标物正北方向为零，顺时针方向逐渐变大，其取值范围是0-360 °。因此太阳方位角一般是以目标物的北方向为起始方向，以太阳光的入射方向为终止方向，按顺时针方向所的角度，如图：图根据“天球理论”，可以导出太阳方位角与太阳直射点纬度、太阳高度角以及观测时刻的关系如下:由图可以得出：由上述方程组可知：又有 sinh sin sin cos cos cost ， 将附件一中的北京时间转化为当地的时角代入方程，并利用题目一所建立的日期与太阳直射点纬度之间的关系模型来求出4月 18日太阳直射点的纬度，最后由MATLAB 解得：8.0378 ，所以求出直杆所在位置为南纬8 2'16

23、9;',东经107 24'18'' 。问题三：模型的建立与求解所在地经度的求解本文首先利用附件2和附件3中所给的影子顶点坐标求出影子长度，建立横坐标为时问M与纵坐标为影子长度L之间的函数关系，利用 MATLAB分别拟合出附件2、附件3数据的曲线方程以及图像，如图和，代码见附录2:图图附件2拟合出来的曲线方程为：附件3拟合出来的曲线方程为：正午太阳高度最大，直杆影子长度最短，由于图像横坐标为北京时间，因此曲线最低 点对应的北京时间为当地正午十二时。对上述所求曲线方程分别进行求导，求出最低点的 北京时间和最短的影子长度：'令 L1 0,得出 M1 15.20

24、79, L1 0.6223。令 L2 0 彳导出 M2 12.7379, L2 3.4809。利用地理知识经度每相差15在时间上相差1小时得出当地的经度与北京的经度分别相 差48 77和11 44 ,由求导得出的结果可知北京时间比当地时间快，因此两地都应在北京 的西面，所以得出两地经度分别为：所在地纬度与日期的求解对附件2,利用问题二的相似三角形模型求出固定直杆的长度，所以0.06 3 106465 3 602.5m 。利用直杆的高度与当地最短的影子长度的比值得出正午太阳高度角H arctan ,求L得H 76.02。取附件二中某一时刻的影子顶点的坐标并求出影子长度，运用三角函数公式求出该时刻

25、的太阳高度角，联立方程组：分两种情况讨论：若当地纬度与太阳直射点的纬度在同一半球，纬差为 -;若当地 纬度与太阳直射点的纬度在不同半球，纬差为 。用MATLAB求解得当地纬度39.87 ,根据问题一建立的求太阳直射点的纬度的模型将代入三个阶段，筛选出符合条件的日期为5月24日。当地的经纬度为（北纬39 52 48，东经68 16 22 ）。同理，对附件三进行求解得，太阳直射点的纬度为3.1430 ,当地经纬度为（南纬51 1015 ,东经105 1925 ）,日期为3月24日。问题四：模型的建立与求解模型假设1）假设直杆底座厚度忽略不计；2）假设忽略风速对直杆高度、影子长度的影响；3）假设拍摄

26、过程中摄像头位置及角度的变化忽略不计。第一小问模型的建立与求解1）所在地纬度的求解本文首先利用视频中所给的日期，利用问题一中建立的求太阳直射点纬度的模型进行180求解，根据公式 arcsin 0.39775 sin d 求出太阳直射点的纬度43 49 36 。186对用MATLAB进行灰度处理的21张图片进行影子长度的测量，测量数据见表，利用 影子的长度变化数据根据三角函数求出对应时刻的太阳高度角，根据问题一建立的求太阳 高度角的数学公式sinh sin sin cos cos cost ,利用MATLAB6对数据进行求解得出 当地的纬度38.8626 。时亥IJ8:548:568:589:0

27、09:029:04影子长度时亥IJ9:069:089:109:129:149:16影子长度时亥IJ9:189:209:229:249:269:28影子长度时亥IJ9:309:329:34影子长度2） 表所在地经度的求解我们以北京时间为标准，求出北京与当地的时间差，从而根据两地的经度差求出当地的经度。假设将此2m 的直杆放在北京，利用达到相同影子长度的时刻求时间差，选取某一时刻的影子长度，求出该时刻的太阳高度角，再利用视频中已知的日期求出当天太阳直射点的纬度，利用问题一建立的模型sinh sin sin cos cos cost 求出时角t 50.2 ，由两地的时角差计算出两地的时间差，继而根据

28、经度每隔15 时间相差一小时可求出分别位r r -*/ -t一一n-t' V. / .-i-r t . r rr 一、 、r'_'''于北京东面和西面的两地经度为21 23'29''和 121 23'29''。._»，ir、一，一rr- »、.'”一因 此 ， 得 出 视 频 大 体 拍 摄 位 置 为 （北纬 38.8626 ，东经 21 2329''）和一 ，一 , "北纬 38.3626 ，东经 121 23'29'' 。问题模型的建立与求解利用 MATLAB 对上述影长与时间进行曲线拟合，图像如图得出的表达式如下：2L 0.07359M 2 2.201M16.16；图当地正午十二时时，直杆影子的长度最短，根据曲线求得最短影长L= 。利用三角函数关系，求出影子最短时的正午太