張印　李偉東　楊洋

【摘 要】本文針對(duì)2015年全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽A題[1]中的第2問(wèn)，依據(jù)最小二乘原理，并利用實(shí)際測(cè)定的參照物的影子長(zhǎng)度隨時(shí)間變化的數(shù)據(jù)，建立了確定圖片拍攝地點(diǎn)的優(yōu)化模型。在Matlab環(huán)境下編程求解，得出參照物的經(jīng)緯度。將此經(jīng)緯度作為已知數(shù)據(jù)，重新計(jì)算參照物的影子長(zhǎng)度，并與實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，發(fā)現(xiàn)二者基本吻合。說(shuō)明模型是合理有效的。

【關(guān)鍵詞】最小二乘法；太陽(yáng)高度角；太陽(yáng)影子定位模型

1 問(wèn)題背景及提出

如何確定視頻的拍攝地點(diǎn)和拍攝日期是視頻數(shù)據(jù)分析的重要方面。其中一種簡(jiǎn)單的情形是：根據(jù)相關(guān)影長(zhǎng)數(shù)據(jù)確定參照物位置，這一技術(shù)也稱為太陽(yáng)影子定位技術(shù)。在本文中，根據(jù)2015年全國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽A題附件1的數(shù)據(jù)[1]，確定出該地點(diǎn)的位置。

2 影長(zhǎng)與時(shí)間的函數(shù)關(guān)系

2.1 影響影長(zhǎng)的各個(gè)參數(shù)

太陽(yáng)高度角是指太陽(yáng)光的入射方向和地平面之間的夾角，可以用太陽(yáng)高度角來(lái)估計(jì)影子場(chǎng)電影。太陽(yáng)高度角h與實(shí)測(cè)地的地理緯度θ，太陽(yáng)直射角緯度φ，時(shí)角Ω滿足以下函數(shù)關(guān)系式[2]：

太陽(yáng)直射點(diǎn)每時(shí)都在向西移動(dòng)，每小時(shí)移過(guò)15度經(jīng)度，在地理題的計(jì)算中可粗略取每天移動(dòng)0.25度緯度[3]。據(jù)此求解出太陽(yáng)直射角緯度與日期之間的關(guān)系，進(jìn)而通過(guò)MATLAB編程，得到了每天太陽(yáng)直射點(diǎn)的緯度，部分表見表1，在 4 月 18日的太陽(yáng)直射點(diǎn)緯度為φ=7。

將H，φ中的數(shù)據(jù)代入（5）中，用 MATLAB 編程計(jì)算，我們得到 21 個(gè)時(shí)間點(diǎn)所對(duì)應(yīng)的經(jīng)度位置，我們發(fā)現(xiàn)，這 21 個(gè)經(jīng)度變化都十分小，故我們經(jīng)度該位置的經(jīng)度就選取求出的 21 個(gè)經(jīng)度的平均值，得到 4 個(gè)經(jīng)度，進(jìn)而得到四個(gè)固定直桿可能所處的地理位置：（18.24N，105.9528E），（18.24N，37.4972E），（7.3969N，107.1988E），（7.3969N，36.2512E）。

4 地理位置的檢驗(yàn)

根據(jù)測(cè)定的關(guān)于x，y 的坐標(biāo)，我們發(fā)現(xiàn)隨著時(shí)間的增加，影子長(zhǎng)度也在不 斷增加，說(shuō)明當(dāng)?shù)貢r(shí)間已經(jīng)過(guò)了正午時(shí)間 12：00，而此時(shí)的北京時(shí)間為 14：42-15：42，結(jié)合北京時(shí)間以東八區(qū) 120 度為標(biāo)準(zhǔn)，每向西 15 度，當(dāng)?shù)貢r(shí)間往前 1 小時(shí)的規(guī)律，我們計(jì)算出在79.5 E 時(shí)，北京時(shí)間為 14：52，當(dāng)?shù)貢r(shí)間恰好為12：00，因此我們所確定的經(jīng)度必須大于 79.5E。進(jìn)而確定可能的地理位置為：（18.24N，105.9528E）和（7.3969N，107.1988 E）。

根據(jù)地圖確定地理位置為（18.24N，105.9528E）的地點(diǎn)為越南，地理位置為（7.3969N，107.1988E）的地點(diǎn)為泰國(guó)灣（為靠近馬來(lái)西亞的一個(gè)海峽），故這個(gè)位置不符合我們的標(biāo)準(zhǔn)，故我們確定的直桿所處位置大致在越南。

5 模型的檢驗(yàn)

我們將所確定的地理位置代入關(guān)于太陽(yáng)高度角的式子（5）中，畫出 14：42-15：42的影子長(zhǎng)度的變化曲線l1，同時(shí)擬合出附件 1 中影子長(zhǎng)度隨時(shí)間變化曲線l2 ，兩者之間的誤差e=||，可以這21個(gè)時(shí)刻的誤差分別為：0.0004，0.0005，0.0005，0.0006，0.0006，0.0006，0.0006， 0.0006，0.0006，0.0006，0.0005，0.0005，0.0005，0.0005，（下轉(zhuǎn)第130頁(yè)）（上接第136頁(yè)）0.0004，0.0004，0.0004，0.0004，0.0004，0.0004，0.0004。兩條影子長(zhǎng)度隨時(shí)間變化的曲線（見圖1）。

【參考文獻(xiàn)】

[1]http：//www.shumo.com/2015cumcm.html[OL].

[2]http：//baike.baidu.com/view/86609.htm[OL].

[3]http：//baike.baidu.com/link？url=_aUuZtIGSo3DBZgjWiSGC2NqH9NNBjEcXX37 JYu5MWsF6M22qioNdjuY81eNjeP31LlpvJKTBhf_ePWF9g5CXa#1[OL].

[4]王國(guó)安，米鴻濤，鄧天宏，李亞男，李蘭霞.太陽(yáng)高度角和日出日落時(shí)刻太陽(yáng)方位角 一年變化范圍的計(jì)算[J].氣象與環(huán)境科學(xué)，2007，30：161-163.

[5]http：//wenku.baidu.com/view/f8fb67f481c758f5f61f6749.html？from=search[OL].

[責(zé)任編輯：楊玉潔]