李　超，姜娟娥，王俊瑞，郭秀娟

（河西學(xué)院 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，甘肅 張掖 734000）

數(shù)學(xué)建模之太陽影子定位模型

李超，姜娟娥，王俊瑞，郭秀娟

（河西學(xué)院 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，甘肅 張掖 734000）

本文針對(duì)太陽影子定位技術(shù)，通過太陽與地球相對(duì)運(yùn)動(dòng)的規(guī)律建立了太陽影子定位模型.利用最小二乘法原理分析了不同地點(diǎn)、不同季節(jié)影子長度的形成規(guī)律及變化趨勢，運(yùn)用M ATLAB軟件進(jìn)行了曲線擬合分析，得出了不同地區(qū)經(jīng)緯度與影子變化的模型，并運(yùn)用建立的模型對(duì)實(shí)際問題進(jìn)行了可行性分析，最后通過分析視頻中物體影長的變化，確定了視頻拍攝的地點(diǎn)和日期。

太陽高度角；赤緯角；M ATLAB；最小二乘法；曲線擬合

一、引言

現(xiàn)代技術(shù)的發(fā)展使得人們能夠方便地記錄高質(zhì)量的視頻文件。在分析視頻材料時(shí)，需要確定視頻的拍攝地點(diǎn)和日期，而利用太陽影子定位技術(shù)對(duì)視頻中的太陽影子變化進(jìn)行分析是確定視頻拍攝的地點(diǎn)和日期的一種重要方法。本文在已有研究的基礎(chǔ)上主要解決四個(gè)問題。問題一：根據(jù)已知條件和經(jīng)度、緯度、時(shí)差之間的內(nèi)在聯(lián)系建立數(shù)學(xué)模型，并用MATLAB制作出直桿的太陽影子的變化曲線圖；問題二：根據(jù)坐標(biāo)數(shù)據(jù)和數(shù)學(xué)模型，確定太陽的赤緯，得出所求地的緯度值，最后由Google地圖軟件搜尋出可能存在的地點(diǎn)；問題三：根據(jù)模型和影子頂點(diǎn)坐標(biāo)得出可能地點(diǎn)的經(jīng)緯度，建立三維立體緯度模型，同時(shí)通過MATLAB對(duì)所得的坐標(biāo)數(shù)據(jù)曲線擬合分析得出對(duì)應(yīng)的經(jīng)、緯度，得出可能的地點(diǎn)；問題四：對(duì)附件中的視頻進(jìn)行分析，抽出每四分鐘的圖像，用灰度處理圖像得出桿長和影長，再與實(shí)際桿長對(duì)比得出該段之間內(nèi)影長變化趨勢，根據(jù)模型確定可能的拍攝地點(diǎn)。

最后對(duì)模型進(jìn)行局部分析和檢驗(yàn)，通過改進(jìn)，該模型給我們提供了一個(gè)從視頻中讀取數(shù)據(jù)再進(jìn)行定位的方法。

二、模型假設(shè)與符號(hào)說明

（一）模型假設(shè)

1.北京時(shí)間為標(biāo)準(zhǔn)時(shí)間。

2.地球是一個(gè)光滑的理想球體。

3.太陽光直射地球途中無遮擋物且不受大氣層影響。

4.觀測影子變化途中天氣晴朗。

（二）符號(hào)說明（見表1）

表1

三、模型的建立與求解

（一）問題一：模型的建立與求解

1.模型的建立。通過對(duì)問題的分析建立天安門廣場(北緯39度54分26秒，東經(jīng)116度23分29秒)3米高直桿的太陽影子長度變化曲線的模型。一天時(shí)間內(nèi)太陽高度角的變化與地方時(shí)和太陽赤緯的變化密切相關(guān)，對(duì)于地球上的某個(gè)地點(diǎn)，太陽高度角隨著地方時(shí)和太陽赤緯角的變化而變化，如圖1所示。

急性缺血性腦卒中(acute ischemic stroke，AIS)具有發(fā)病急和癥狀重等特點(diǎn)[1]。AIS的發(fā)病機(jī)制與機(jī)體內(nèi)的炎癥反應(yīng)關(guān)聯(lián)較大[2]。趨化因子12(chemokine 12，CXCL12)對(duì)造血祖細(xì)胞產(chǎn)生調(diào)控，募集機(jī)體內(nèi)的神經(jīng)前體細(xì)胞，使骨髓源型祖細(xì)胞不斷介導(dǎo)病灶區(qū)的炎癥反應(yīng)，最終誘導(dǎo)缺血組織血管的再生[3]。白細(xì)胞介素-33(Interleukin-33，IL-33)能對(duì)基因轉(zhuǎn)錄實(shí)施調(diào)控，在機(jī)體的炎癥反應(yīng)中起重要作用[4]。本研究通過分析AIS患者血清CXCL12、IL-33水平變化及臨床意義，旨在為臨床AIS患者預(yù)后的評(píng)估提供依據(jù)?，F(xiàn)報(bào)道如下。

圖1

圖2

通過圖形和各個(gè)角度的關(guān)系得到太陽高度角的計(jì)算公式，從而建立了影子長度變化的三角函數(shù)模型：，根據(jù)太陽赤緯角計(jì)算公式

其中

圖2為固定直桿在太陽光照射下的影子，其中H=3 (單位:米)，可得固定直桿影長公式為聯(lián)立以上式子可得影長公式為：

2.模型的求解。利用Matlab軟件編寫程序分析影子長度關(guān)于太陽高度角、太陽赤緯、地方時(shí)日期時(shí)間的變化規(guī)律，可得2015年10月22日北京時(shí)間9:00-15:00之間天安門廣場(北緯39度54分26秒，東經(jīng)116度23分29秒)3米高的直桿的太陽影子長度的變化曲線，如圖3所示。

圖3

圖4

根據(jù)圖像分析，在北京9:00到12:00之間太陽照射固定直桿的影子從大約6米逐漸遞減到大約3.6米，正午12:00以后太陽的影子又開始逐漸增長。

（1）在不同的地方直桿影子的變化趨勢基本相同，但在同一時(shí)刻直桿的影子隨著經(jīng)緯度的變化而變長。

（2）在每個(gè)地區(qū)直桿影子最短的時(shí)刻是在中午12點(diǎn)左右，但在12:00之前直桿的影子由長變短，12點(diǎn)以后直桿的影子逐漸變長。利用Matlab軟件編寫程序繪制不同日期、相同地方、相同時(shí)間段內(nèi)的直桿影子變化情況，如圖5所示。

圖5

圖6

根據(jù)圖像可知，相同長度直桿的影子的變化趨勢基本相同，但是在不同的日期、相同的時(shí)間點(diǎn)上影子的長度不相同。在中午12時(shí)，7月的影子最短，1月的影子最長。由此可知，影子的長度與季節(jié)有關(guān)，而季節(jié)的交替變更與日期有關(guān)。這說明影子的長度受日期的影響。

（二）問題二：模型的建立與求解

為了解決問題二，我們需運(yùn)用最小二乘法基本原理③求出誤差平方和最小的擬合。在實(shí)驗(yàn)過程中實(shí)驗(yàn)數(shù)據(jù)不可能精確獲得，往往有一定的誤差。為了盡量減少誤差，需要用一個(gè)擬合函數(shù) f（ x），使擬合函數(shù) f（ x）在xi處的函數(shù)值f（ xi）與對(duì)應(yīng)的實(shí)驗(yàn)實(shí)測數(shù)據(jù)值yi（ i=1，2，3…N）形成誤差平方和最小值這就是最小二乘法。

通過分析上述問題，根據(jù)固定直桿的太陽影子頂點(diǎn)坐標(biāo)數(shù)據(jù)，建立了三角函數(shù)數(shù)學(xué)模型，并且分析了直桿影子的變化趨勢。在此基礎(chǔ)上，通過反解倒推確定直桿所處的經(jīng)緯度，建立了直桿影長與經(jīng)緯度相關(guān)的數(shù)學(xué)表達(dá)式，并且應(yīng)用最小二乘法曲線擬合的方法確定直桿所處的地點(diǎn)進(jìn)行擬合分析得。因赤緯角日變化很小，一年內(nèi)任何一天的赤緯角用下式計(jì)算：。其中N為日數(shù)，自每年1月1日開始計(jì)算P為真太陽時(shí)，a為北京標(biāo)準(zhǔn)差，Q為時(shí)差，（x,y）表示影子在平面內(nèi)的坐標(biāo)點(diǎn)p=a+q，，222lxy=+， 聯(lián)立這四個(gè)式子建立緯度與經(jīng)度的數(shù)學(xué)模型為:

應(yīng)用Matlab軟件編寫的程序?qū)ι鲜鼋⒌臄?shù)學(xué)模型進(jìn)行擬合分析，如圖6所示。該圖是緯度117.78度，經(jīng)度25.245度時(shí)的擬合情況。根據(jù)圖像的變化趨勢和最小二乘法的原理可以估計(jì)出實(shí)際的直桿位置，實(shí)際影子的長度與模型影子的長度之間的吻合程度比較接近。在14:42到15:20的時(shí)間段內(nèi)，實(shí)際影長的曲線與模擬影長的曲線之間的距離漸漸接近，到15:20左右以后實(shí)際影長與模型影長之間曲線的接近程度減弱。

圖7

圖8

圖7為緯度116.34度、經(jīng)度25.248度時(shí)的擬合情況，表示的是實(shí)際影子長度和模型影子長度的擬合情況，由圖像的變化趨勢和最小二乘法的原理可以觀測出在14.6到15的范圍內(nèi)有重合點(diǎn)，在范圍以外也很接近，預(yù)測模型和實(shí)際模型的擬合度好，則預(yù)測模型的緯度116.34度、經(jīng)度為25.248度的可以近似代替的經(jīng)緯度。

（三）問題三：模型的建立與求解

根據(jù)某固定直桿在水平地面上的太陽影子頂點(diǎn)坐標(biāo)數(shù)據(jù)，應(yīng)用上述所建立的模型一與模型二采用了最小二乘法非線性擬合的方法變量經(jīng)緯度與日期進(jìn)行預(yù)測，用MATLAB軟件編程繪制出的影子長度與時(shí)間的變化曲線圖，如圖8所示。

為檢驗(yàn)?zāi)Ｐ偷臏?zhǔn)確性，在地理位置緯度為160.51度，經(jīng)度為5.5869度，積日為329.54秒內(nèi)太陽光照射下模型影長與實(shí)際影長之間的長度較為接近，而且由圖8可以進(jìn)一步看出在13.10這一時(shí)間左右直桿的實(shí)際影長與模型影長大體一致，從而可以說明該模型可實(shí)施。

（四）問題四：模型的建立與求解

表2　太陽影子長度匯總表

根據(jù)視頻中已知直桿的影子變化情況，抽取出每兩分鐘的圖像，再通過灰度處理圖像得出視頻中直桿影子的長度隨時(shí)間變化而變化，如表2所示。所得數(shù)據(jù)應(yīng)用所計(jì)算的數(shù)據(jù)建立有關(guān)經(jīng)度與維度的模型，通過該模型的建立可以估計(jì)出測量時(shí)的地點(diǎn)，如表3所示。

表3　估測拍攝地點(diǎn)經(jīng)緯度匯總表

如果拍攝日期未知，可以利用以上問題二的模型來確定拍攝的日期。

四、模型的評(píng)價(jià)

（一）模型優(yōu)點(diǎn)

本文經(jīng)過數(shù)據(jù)的合理處理，采用三角函數(shù)模型、最小二乘法模型并對(duì)模型，進(jìn)行局部分析和檢驗(yàn)。該模型給我們提供了一個(gè)從視頻中讀取數(shù)據(jù)并進(jìn)行定位的方法。它在實(shí)際生活中具有較強(qiáng)的適用性，應(yīng)用前景非常廣闊。

（二）模型缺點(diǎn)

由于所采用的數(shù)據(jù)樣本、讀取數(shù)據(jù)和測量時(shí)的誤差，對(duì)模型造成了數(shù)據(jù)對(duì)比方面的不平衡，加上考慮問題時(shí)忽視了天氣變化對(duì)現(xiàn)實(shí)中太陽直射地面時(shí)桿影的影響，從而導(dǎo)致桿影的測量值與真實(shí)值之間存在誤差。另外，在視頻的截取技術(shù)方面還存在一定的誤差，導(dǎo)致測量視頻中直桿影子的長度時(shí)存在誤差。

[1] 陳東彥,李冬梅,王樹忠.數(shù)學(xué)建模(第一版)[M].北京:科學(xué)出版社,2007.

[2] 蔡山,張浩,陳洪輝,沙基昌.基于最小二乘法的分段三次曲線擬合方法研究[J].科學(xué)技術(shù)與工程,2007(3).

[3] 李植華,王百眾,王元知.關(guān)聯(lián)度分析[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究,1988(2).

[4] 郭仁忠,張克權(quán).Q型聚類分析中變量相關(guān)性的處理方法分析[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,1987(3).

［責(zé)任編輯 房曉偉］

G64

A

1673-9132(2016)32-0229-03

10.16657/j.cnki.issn1673-9132.2016.32.146

郭秀娟（1981- ），女，漢族，碩士研究生，講師，研究方向：從事數(shù)學(xué)課程與教學(xué)論研究。

郭秀娟

河西學(xué)院大學(xué)生科技創(chuàng)新項(xiàng)目——數(shù)學(xué)建模之太陽影子定位模型（2015-114）。