太陽(yáng)影子定位技術(shù)研究

太陽(yáng)影子定位技術(shù)的研究 摘要 影子在生活中非常常見(jiàn)，而太陽(yáng)影子與物體本身、物體所在的地理位置、當(dāng)?shù)氐娜掌诤彤?dāng)天時(shí)間都有密切的關(guān)系，本文主要研究影長(zhǎng)隨物體本身的高度（）、物體所在的地理位置（與有關(guān)）、當(dāng)?shù)氐娜掌冢ㄅc有關(guān)）和當(dāng)天時(shí)間（與t有關(guān)）的變化規(guī)律，從而根據(jù)影子變化的相關(guān)信息推測(cè)出物體所處的位置。首先，我們查找到了包含物體本身的高度、當(dāng)?shù)貢r(shí)間、季節(jié)和地理位置這四個(gè)因素的太陽(yáng)高度角公式：，然后采用控制變量法分別做出影長(zhǎng)隨單個(gè)因素的變化圖像，進(jìn)而分析影長(zhǎng)與各個(gè)因素之間的聯(lián)系。當(dāng)物體本身的高度、季節(jié)、地理位置都已知時(shí)，則可以直接作出影子長(zhǎng)度隨時(shí)間的變化關(guān)系圖。由圖可知：2015年10月22日天安門(mén)廣場(chǎng)3米高的直桿的太陽(yáng)影子在9:00到15:00之間的變化大致呈對(duì)稱(chēng)分布，最小影長(zhǎng)為：3.663米、最大影長(zhǎng)為：6.616米。其次，考慮到影子頂點(diǎn)的位置隨太陽(yáng)方位角會(huì)變化，所以需要考慮太陽(yáng)方位角（），對(duì)問(wèn)題一的模型進(jìn)行修正，消去桿長(zhǎng)之后，利用附件1中每隔3分鐘的影子長(zhǎng)度變化差列出方程組，解得附件1的地點(diǎn)可能為：（北緯18度46分2秒,東經(jīng)111度23分29秒），即海南附近。然后，我們根據(jù)附錄二，繪制影長(zhǎng)隨時(shí)間的變化曲線，再對(duì)其進(jìn)行擬合分析，得到當(dāng)?shù)亟?jīng)度為74度26分。由于緯度和赤緯不隨時(shí)間變化，所以我們采用最小二乘法進(jìn)行擬合，最終解得、、N=275，即在新疆地區(qū)（北緯39度20分，東經(jīng)74度26分），日期為10月2日。同時(shí)對(duì)于附錄三，采用同樣的處理辦法，得到、、N=325，即在湖北地區(qū)（北緯31度27分，東經(jīng)108度44分），日期為11月23日。最后，以視屏中旗桿所在地為圓心建立坐標(biāo)系，每隔3分鐘對(duì)視屏進(jìn)行截屏，得到一組和附件1類(lèi)似的桿影的頂點(diǎn)坐標(biāo)隨時(shí)間變化的數(shù)據(jù)，把模型二對(duì)附件一的處理方法運(yùn)用到這組數(shù)據(jù)中，確定出視頻的拍攝位置在陜西境內(nèi)。當(dāng)拍攝日期未知時(shí)，采用模型三處理這組數(shù)據(jù)，也得到視頻的拍攝地點(diǎn)在陜西境內(nèi)。關(guān)鍵詞：太陽(yáng)高度角太陽(yáng)方位角控制變量最小二乘法擬合一問(wèn)題重述視頻數(shù)據(jù)分析在現(xiàn)代生活非常常見(jiàn)，而確定視頻的拍攝地點(diǎn)和拍攝日期更是其重中之重，太陽(yáng)影子定位技術(shù)就是通過(guò)分析視頻中物體的太陽(yáng)影子變化，確定視頻拍攝的地點(diǎn)和日期的一種方法。太陽(yáng)影子的變化與物體本身的高度及所處的地理位置、時(shí)間和季節(jié)等因素有密切的關(guān)系，知道其中的幾個(gè)因素便可以對(duì)物體所處地點(diǎn)進(jìn)行定位。因而我們需要建立數(shù)學(xué)模型來(lái)確定這些因素之間的關(guān)系，對(duì)此進(jìn)行以下幾步工作：一、建立影子長(zhǎng)度變化的數(shù)學(xué)模型，分析影子長(zhǎng)度關(guān)于各個(gè)參數(shù)的變化規(guī)律，并應(yīng)用所建立的模型畫(huà)出2015年10月22日北京時(shí)間9:00-15:00之間天安門(mén)廣場(chǎng)（北緯39度54分26秒,東經(jīng)116度23分29秒）3米高的直桿的太陽(yáng)影子長(zhǎng)度的變化曲線。二、根據(jù)某固定直桿在水平地面上的太陽(yáng)影子頂點(diǎn)坐標(biāo)數(shù)據(jù)，建立數(shù)學(xué)模型確定直桿所處的地點(diǎn)。并將所建的模型應(yīng)用于附件1的影子頂點(diǎn)坐標(biāo)數(shù)據(jù)，得到若干個(gè)可能的地點(diǎn)。三、根據(jù)某固定直桿在水平地面上的太陽(yáng)影子頂點(diǎn)坐標(biāo)數(shù)據(jù)，建立數(shù)學(xué)模型確定直桿所處的地點(diǎn)和日期。并將所建的模型分別應(yīng)用于附件2和附件3的影子頂點(diǎn)坐標(biāo)數(shù)據(jù)，得到若干個(gè)可能的地點(diǎn)與日期。四、建立確定視頻拍攝地點(diǎn)的數(shù)學(xué)模型，附件4為一根直桿在太陽(yáng)下的影子變化的視頻，并且已通過(guò)某種方式估計(jì)出了直桿的高度為2米。應(yīng)用所建立的模型給出若干個(gè)可能的拍攝地點(diǎn)。二問(wèn)題分析太陽(yáng)影子的變化與物體本身的高度、時(shí)間、季節(jié)和地理位置等因素都有重要關(guān)系，因此研究影子的變化規(guī)律有重要的意義。針對(duì)問(wèn)題一，我們要查找包含物體本身的高度、當(dāng)?shù)貢r(shí)間、季節(jié)和地理位置的函數(shù)表達(dá)式，再對(duì)影長(zhǎng)隨單個(gè)因素的變化趨勢(shì)進(jìn)行分析。當(dāng)物體本身的高度、季節(jié)、地理位置都是已知時(shí)，可以直接作出影子長(zhǎng)度隨時(shí)間的變化關(guān)系圖。針對(duì)問(wèn)題二，在問(wèn)題一的基礎(chǔ)上，考慮到影子頂點(diǎn)的位置隨太陽(yáng)方位角會(huì)變化，所以需要對(duì)問(wèn)題一的模型進(jìn)行修正，由于h和桿高、影長(zhǎng)聯(lián)系在一起的，而附件1中的桿高未知，所以要先把它代入之前的太陽(yáng)高度角公式消去h,得到矯正以后的太陽(yáng)高度角公式，然會(huì)利用附件1中每隔3分鐘的影子長(zhǎng)度變化差，列出方程組，解得當(dāng)?shù)貢r(shí)間t,用它和北京時(shí)間的時(shí)間差算出它與北京的經(jīng)度差，進(jìn)而算出當(dāng)?shù)氐慕?jīng)度。同時(shí)，利用附件1中的影長(zhǎng)、日期和當(dāng)?shù)氐臅r(shí)間，可以算出當(dāng)?shù)氐木暥取?duì)于問(wèn)題三：我們根據(jù)附錄，繪制影長(zhǎng)隨時(shí)間的變化曲線，再對(duì)其進(jìn)行擬合分析，得到函數(shù)關(guān)系，根據(jù)函數(shù)關(guān)系，我們預(yù)測(cè)太陽(yáng)影長(zhǎng)最短時(shí)的北京時(shí)間、計(jì)算時(shí)差，又由于緯度和赤緯不隨時(shí)間的變化，所以我們采用最小二乘法進(jìn)行擬合，確定參數(shù)和。對(duì)于問(wèn)題四:由于附件4給了一個(gè)視頻，所以我們先每隔3分鐘對(duì)視頻截屏，以旗桿與地面的交點(diǎn)為圓心建立坐標(biāo)系，得到影的頂點(diǎn)坐標(biāo)，把視屏數(shù)據(jù)轉(zhuǎn)化為類(lèi)似于附件1的數(shù)字?jǐn)?shù)據(jù)，再運(yùn)用模型二對(duì)這組數(shù)據(jù)列方程組，求解出視頻的拍攝地點(diǎn)。三模型假設(shè)1.假設(shè)太陽(yáng)是固定不動(dòng)的。2.假設(shè)地球離太陽(yáng)的距離無(wú)窮遠(yuǎn)，太陽(yáng)光線是平行光。3.假設(shè)地球?yàn)榍蝮w。4.假設(shè)地球的運(yùn)行軌道為圓。5.假設(shè)大氣層均勻，太陽(yáng)光線穿過(guò)大氣層時(shí)沒(méi)有折射。6.不考慮高山遮擋、海拔高度等因素的影響。四符號(hào)說(shuō)明1.：太陽(yáng)高度角。2.：地球緯度。3.：太陽(yáng)赤緯。4.：時(shí)角。5.：太陽(yáng)方位角。6.：影長(zhǎng)。7.：物體高度。8.T：當(dāng)?shù)貢r(shí)間。9.x：影子的橫坐標(biāo)。10.y：影子的縱坐標(biāo)。五模型的建立與求解5.1建立模型一太陽(yáng)影子的變化與物體的本身的高度及所處的地理位置、時(shí)間和季節(jié)等參數(shù)有密切的關(guān)系。根據(jù)地理學(xué)知識(shí)，我們知道，物體高度越高，太陽(yáng)影子越長(zhǎng)，物體所處緯度越遠(yuǎn)離太陽(yáng)直射點(diǎn)的緯度，太陽(yáng)影子越長(zhǎng)，查閱資料得：（1）其中，表示太陽(yáng)高度角，表示地理緯度，表示太陽(yáng)赤緯，表示時(shí)角。時(shí)角t是以正午十二點(diǎn)為零度開(kāi)始算，每小時(shí)為15度。即十四點(diǎn)和十點(diǎn)分別為30度和-30度。太陽(yáng)赤緯是地球赤道平面與太陽(yáng)和地球中心的連線之間的夾角。因赤緯值日變化很小，一年內(nèi)任何一天的赤緯角δ可用下式計(jì)算：（2）（3）式中N為日數(shù)，自每年1月1日開(kāi)始計(jì)算。由于幾何關(guān)系，太陽(yáng)高度角與影長(zhǎng)有如下關(guān)系（4）其中為物高，為影長(zhǎng)，聯(lián)立（1）、（4）可得影長(zhǎng)與各個(gè)參數(shù)的關(guān)系，即為最終所建數(shù)學(xué)模型。（5）5.1.1模型一的求解通過(guò)對(duì)模型的分析，我們知道太陽(yáng)影子的變化與物體的本身的高度及所處的地理位置、時(shí)間和季節(jié)等參數(shù)有密切的關(guān)系，因此我們采用控制變量法來(lái)分別分析影子長(zhǎng)度與各個(gè)因素之間的關(guān)系：(1)太陽(yáng)影子和緯度的變化關(guān)系：我們保持物體的本身的高度及時(shí)間、季節(jié)等參數(shù)不變，單獨(dú)分析太陽(yáng)影子和緯度之間的變化關(guān)系。對(duì)此我們作了以下圖像：圖1太陽(yáng)影子隨緯度變化通過(guò)對(duì)圖1分析可知，對(duì)于同一日期，同一時(shí)刻，緯度由直射點(diǎn)向北或向南變化時(shí)，太陽(yáng)影子的長(zhǎng)度總是越來(lái)越長(zhǎng)。(2)太陽(yáng)影子和物高的變化關(guān)系我們保持物體所處的緯度及時(shí)間、季節(jié)等參數(shù)不變，單獨(dú)分析太陽(yáng)影子和物體高度的變化關(guān)系。對(duì)此我們作了以下圖像：圖2太陽(yáng)影長(zhǎng)隨物高的變化通過(guò)圖2分析可知，同一日期，同一時(shí)刻，同一緯度不同物體的太陽(yáng)影子與物體的高度成正比。(3)太陽(yáng)影子長(zhǎng)度和季節(jié)的變化關(guān)系我們保持物體本身高度及所處的緯度和時(shí)間不變，單獨(dú)分析太陽(yáng)影子和季節(jié)的變化關(guān)系。太陽(yáng)影子的變化與物體所處緯度的不同呈不同的規(guī)律:a.當(dāng)緯度大于時(shí)，保持物體本身高度及所處的緯度和時(shí)間不變，太陽(yáng)影子長(zhǎng)度隨日期的變化如下圖所示：圖3南北回歸線外影長(zhǎng)隨日期的變化由圖知，物體所處緯度大于北回歸線度數(shù)時(shí)，太陽(yáng)影子長(zhǎng)度自1月1日起隨日期的增加先減小后增加。b.當(dāng)緯度介于和之間時(shí)，保持物體本身高度及所處的緯度和時(shí)間不變，太陽(yáng)影子長(zhǎng)度隨日期的變化如下圖所示：圖4南北回歸線內(nèi)影長(zhǎng)隨日期的變化由圖可知，物體所處緯度小于北回歸線度數(shù)時(shí)，由于太陽(yáng)直射點(diǎn)的移動(dòng)，太陽(yáng)影子長(zhǎng)度自1月1日起隨日期的增加先減小后增大再減小最后又增大。對(duì)于問(wèn)題一，應(yīng)用我們所建立的模型，對(duì)2015年10月22日北京9:00-15:00之間天安門(mén)廣場(chǎng)（北緯39度54分26秒,東經(jīng)116度23分29秒）3米高的直桿的太陽(yáng)影子長(zhǎng)度的變化進(jìn)行了如下求解：聯(lián)立(2)、（3）式可知當(dāng)N=295，把和(39.97度)代入(5)得到太陽(yáng)影長(zhǎng)y隨地方時(shí)t的變化關(guān)系：圖5影長(zhǎng)在9點(diǎn)到12點(diǎn)的變化規(guī)律由圖5可知，2015年10月22日天安門(mén)廣場(chǎng)3米高的直桿的太陽(yáng)影子在9:00到15:00之間的變化規(guī)律為：先減小后增大，呈對(duì)稱(chēng)分布，太陽(yáng)影子長(zhǎng)度分別在9：00和15：00點(diǎn)時(shí)出現(xiàn)最大值，最大影長(zhǎng)約為：6.616米，在12：00左右時(shí)出現(xiàn)最小值，最小影長(zhǎng)為：3.663米。5.2建立模型二經(jīng)過(guò)對(duì)模型一的分析，我們知道，太陽(yáng)影子的長(zhǎng)短與太陽(yáng)高度角有關(guān)，而太陽(yáng)影子的朝向卻與太陽(yáng)方位角有關(guān)，因此，我們通過(guò)查閱資料，得到太陽(yáng)方位角有如下公式(6)其中為太陽(yáng)方位角陰影長(zhǎng)度的南北分量陰影長(zhǎng)度的南北分量遮擋物高度太陽(yáng)方位角陰影長(zhǎng)度太陽(yáng)光線h陰影長(zhǎng)度的東西分量太陽(yáng)高度角圖6太陽(yáng)方位角示意圖又根據(jù)圖6的幾何關(guān)系知：(7)其中，x為影子的橫坐標(biāo)。通過(guò)聯(lián)立（1）（6）（7）可得模型二：（8）5.2.1模型二的求解附件一提供了影子在一段時(shí)間內(nèi)的坐標(biāo)，利用上述模型，根據(jù)相同時(shí)間間隔內(nèi)影長(zhǎng)的增長(zhǎng)量得到下列方程組：（9）其中由（7）式給出，由（2）式給出。運(yùn)用mathematica8.0解方程組（9）得出（單位：弧度），則由公式：（10）得當(dāng)?shù)貢r(shí)間T=14.368時(shí)=14時(shí)22分。它與北京時(shí)間的時(shí)間差為20分鐘，由時(shí)差可得該地的經(jīng)度與北京的經(jīng)度相差5度，所以該地的經(jīng)度為111度23分29秒。聯(lián)立（1）、（6）式得：（11）將問(wèn)題二中已知的、值及求解出的值代入（11）式中，解得:=20度46分2秒，綜上所述，附件1可能為：（北緯20度46分2秒,東經(jīng)111度23分29秒），即海南附近。5.3模型三的建立首先，我們根據(jù)附錄二，繪制影長(zhǎng)隨時(shí)間的變化曲線，再對(duì)其進(jìn)行擬合分析得到函數(shù)關(guān)系，預(yù)測(cè)到當(dāng)太陽(yáng)影長(zhǎng)最短時(shí)的北京時(shí)間，得到當(dāng)?shù)貢r(shí)間與北京時(shí)間的差值，再根據(jù)殘差進(jìn)行修正，得到當(dāng)?shù)氐慕?jīng)度。而緯度和赤緯不隨時(shí)間的變化，所以我們采用最小二乘法對(duì)公式（8）進(jìn)行擬合，確定參數(shù)和，并運(yùn)用mathematical8.0進(jìn)行求解，確定的值。5.3.1模型三的求解對(duì)附件2進(jìn)行曲線擬合得到圖7：圖7附件2的曲線圖由圖7得到函數(shù)關(guān)系：，預(yù)測(cè)當(dāng)太陽(yáng)影子長(zhǎng)度最短時(shí)北京時(shí)間是15:12，即當(dāng)?shù)貢r(shí)間與北京時(shí)間的差值為3.2小時(shí)，經(jīng)過(guò)修正后，時(shí)間差值為2.73小時(shí)，即當(dāng)?shù)亟?jīng)度與北京經(jīng)度的差值為40.95°，得到當(dāng)?shù)氐慕?jīng)度為74度26分。由最小二乘法對(duì)公式（8）進(jìn)行擬合，最終確定、、N=275。即在新疆地區(qū)（北緯39度20分，東經(jīng)74度26分），日期為10月2日。同上，我們對(duì)于附錄三采同樣的處理方法，擬合曲線得到圖8：圖8附件3的曲線圖由圖8得到函數(shù)關(guān)系:,預(yù)測(cè)當(dāng)太陽(yáng)影子長(zhǎng)度最短時(shí)北京時(shí)間是12:44，即當(dāng)?shù)貢r(shí)間與北京時(shí)間的差值為0.73個(gè)小時(shí)，經(jīng)過(guò)修正后，時(shí)間差值為0.51個(gè)小時(shí)，即當(dāng)?shù)亟?jīng)度與北京經(jīng)度的差值為7.65°，得到當(dāng)?shù)氐慕?jīng)度為108度44分。由最小二乘法對(duì)公式（8）進(jìn)行擬合，最終確定，，N=325，即在湖北地區(qū)（北緯31度27分，東經(jīng)108度44分），日期為11月23日。5.4模型四的建立與求解：附件4只給了一個(gè)視頻，所以我們先每隔3分鐘對(duì)視頻截屏，以旗桿與地面的交點(diǎn)為圓心建立坐標(biāo)系，得到桿影的頂點(diǎn)坐標(biāo)，把視屏數(shù)據(jù)轉(zhuǎn)化為類(lèi)似于附件1的數(shù)字?jǐn)?shù)據(jù)：表1視頻中影長(zhǎng)隨時(shí)間的變化關(guān)系時(shí)間：7月13號(hào)北京時(shí)間X坐標(biāo)Y坐標(biāo)影長(zhǎng)9:0300280.0523639:0619.760.186.352310.0523639:0919.420.1584.867930.0523619:1218.420.1582.246010.0523619：1518.410.1680.454740.0523619：1818.20.2179.539290.0523589：2117.830.2177.92250.0523579：2417.460.3276.313010.0523479：2717.040.3274.477930.0523469：3016.770.3273.298240.0523459：3316.560.3772.385260.052338要確定視頻拍攝地的地理位置，就必須知道桿的高度、太陽(yáng)影子的變化規(guī)律及當(dāng)?shù)氐臅r(shí)間。而題中桿的高度、當(dāng)?shù)貢r(shí)間都是未知的，因此我們首先要消去桿高（）這個(gè)未知量，進(jìn)而求得當(dāng)?shù)貢r(shí)間（t），對(duì)此我們進(jìn)行如下處理：聯(lián)立（1）、（5）、（6）并利用相同時(shí)間間隔內(nèi)影長(zhǎng)的增長(zhǎng)量得到下列方程組：（12）N=194（13）聯(lián)立（2）、（3）、（13）式解得：0.383（14）把、、、代入（12）式解得：t=-0.883386（單位：弧度）則由公式（9）得:當(dāng)?shù)貢r(shí)間T=8：38,它與北京時(shí)間的時(shí)間差為25分鐘，由時(shí)差可得該地的經(jīng)度與北京的經(jīng)度相差6.25度，所以該地的經(jīng)度為110度8分29秒。把t=-0.5（單位：弧度）、=0.383、=0.999975代入（10）式，解得:=36度40分9秒綜上所述，視頻拍攝地大概位于（北緯36度40分9秒,東經(jīng)110度8分29秒），在陜西省境內(nèi)。當(dāng)拍攝日期未知時(shí)，經(jīng)分析發(fā)現(xiàn)，和問(wèn)題三相似，所以我們采用模型三進(jìn)行求解：首先，我們根據(jù)表1，繪制影長(zhǎng)隨時(shí)間的變化曲線，再對(duì)其進(jìn)行擬合分析得到函數(shù)關(guān)系，預(yù)測(cè)到當(dāng)太陽(yáng)影長(zhǎng)最短時(shí)的北京時(shí)間，得到當(dāng)?shù)貢r(shí)間與北京時(shí)間的差值，再根據(jù)殘差進(jìn)行修正，得到當(dāng)?shù)氐慕?jīng)度。而緯度和赤緯不隨時(shí)間的變化，所以我們采用最小二乘法對(duì)公式（8）進(jìn)行擬合，確定參數(shù)和，并運(yùn)用mathematical8.0進(jìn)行求解，確定的值。對(duì)表2中的數(shù)據(jù)進(jìn)行曲線擬合得到圖9：圖9視頻中影長(zhǎng)隨時(shí)間的變化關(guān)系圖由圖9得到函數(shù)關(guān)系：，預(yù)測(cè)當(dāng)太陽(yáng)影子長(zhǎng)度最短時(shí)北京時(shí)間是12：25，即當(dāng)?shù)貢r(shí)間與北京時(shí)間的差值為0.47個(gè)小時(shí)，經(jīng)過(guò)修正后，時(shí)間差值為0.41個(gè)小時(shí)，即當(dāng)?shù)亟?jīng)度與北京經(jīng)度的差值為6.25°，得到當(dāng)?shù)氐慕?jīng)度為110度15分。由最小二乘法對(duì)公式（8）進(jìn)行擬合，最終確定36度40分，0.385，N=193，即在西安地區(qū)（北緯36度40分，東經(jīng)110度15分），日期為7月12日。利用模型三所求解的日期與視頻所給出的日期在誤差允許的范圍內(nèi)，說(shuō)明我們的模型合理。參考文獻(xiàn)[1]百度百科，http//.www.baike.baidu,2015年9月11日星期一[2]何銀濤，張梅，黃華，垂直單軸跟蹤光伏支架方針間距計(jì)算界面設(shè)計(jì)，太陽(yáng)能，2015（7）[3]武琳，基于太陽(yáng)陰影軌跡的經(jīng)緯度估計(jì)技術(shù)研究，天津大學(xué)，2010。附錄研究參數(shù)規(guī)律PrivateSubCommand1_Click()N=Text1.Textb=(2*3.1415926*(N-1)/365)delta=(0.006918-0.399912*Cos(b)+0.070257*Sin(b)-0.006758*Cos(2*b)+0.000907*Sin(2*b)-0.002697*Cos(3*b)+0.00148*Sin(3*b))x=3t=0phi=10*3.1415926/180y=Sqr((x/(Sin(phi)*Sin(delta)+Cos(phi)*Cos(delta)*Cos(t)))^2-x^2)PrintyEndSub求解方位角F=(Cos[t1]-Cos[t1+3])/(Cos[t1+3]-Cos[t1+6])*(Sqrt[1-1.217*(Sin[t1+3])2]-Sqrt[1-1.208*(Sin[6+t1])2])+Sqrt[1-1.217(Sin[t1+3])2]-Sqrt[1-1.226*(Sin[t1])2];Plot[F,{t1,0,/2}]Export["E:/a.txt",Table[{t1,F},{t1,0,/2,0.001}],"Table"]求解時(shí)間與影長(zhǎng)PrivateSubCommand1_Click()DimPiAsDoublePi=4*Atn(1)phi=39.906*Pi/180delta=-0.1880374t=(Text3.Text)*Pi/180y=Sqr((3/(Sin(phi)*Sin(delta)+Cos(phi)*Cos(delta)*Cos(t)))^2-9)PrintyEndSub問(wèn)題二代碼Solve[Sqrt[1-((Cos[0.1841]*Sin[1.01])/(1.0365/Sqrt[1.0365^2+0.4973^2]))^2]==Sin[\[Phi]]*Sin[0.1841]+Cos[\[Phi]]*Cos[0.1841]*Cos[1.01]]問(wèn)題三代碼fp=Table[Prime[x],{x,20}]{1.109142,1.104962,1.101012,1.097295,1.093757,1.090422,\1.087214,1.084188,1.08132,1.078567,1.075945,1.073407,1.071014,\1.068697,1.066489,1.064357,1.062298,1.060331,1.058438,1.056612,\1.054833}FindFit[fp,(Tan[d]*Sin[f])/(x-Cos[f]),{d,f},x]目錄TOC\o"1-2"\h\z\u第一章總論 11.1項(xiàng)目名稱(chēng)與承辦單位 11.2研究工作的依據(jù)、內(nèi)容及范圍 11.3編制原則 31.4項(xiàng)目概況 31.5技術(shù)經(jīng)濟(jì)指標(biāo) 51.6結(jié)論 6第二章項(xiàng)目背景及建設(shè)必要性 82.1項(xiàng)目背景 82.2建設(shè)的必要性 9第三章建設(shè)條件 113.1項(xiàng)目區(qū)概況 113.2建設(shè)地點(diǎn)選擇 錯(cuò)誤！未定義書(shū)簽。3.3項(xiàng)目建設(shè)條件優(yōu)劣勢(shì)分析 錯(cuò)誤！未定義書(shū)簽。第四章市場(chǎng)分析與銷(xiāo)售方案 134.1市場(chǎng)分析 134.2營(yíng)銷(xiāo)策略、方案、模式 14第五章建設(shè)方案 155.1建設(shè)規(guī)模和產(chǎn)品方案 155.2建設(shè)規(guī)劃和布局 155.3運(yùn)輸 185.4建設(shè)標(biāo)準(zhǔn) 185.5公用工程 205.6工藝技術(shù)方案 215.7設(shè)備方案 215.8節(jié)能減排措施 24第六章環(huán)境影響評(píng)價(jià) 256.1環(huán)境影響 256.2環(huán)境保護(hù)與治理措施 266.3評(píng)價(jià)與審批 28第七章項(xiàng)目組織與管理 297.1組織機(jī)構(gòu)與職能劃分 297.2勞動(dòng)定員 297.3經(jīng)營(yíng)管理措施 307.4技術(shù)培訓(xùn) 30第八章勞動(dòng)、安全、衛(wèi)生與消防 318.1編制依據(jù)及采用的標(biāo)準(zhǔn) 318.2安全衛(wèi)生防護(hù)原則 318.3自然災(zāi)害危害因素分析及防范措施 328.4生產(chǎn)過(guò)程中產(chǎn)生的危害因素分析及防范措施 328.5消防編制依據(jù)及采用的標(biāo)準(zhǔn) 348.6消防設(shè)計(jì)原則 358.7火災(zāi)隱患分析 358.8總平面消防設(shè)計(jì) 358.9消防給水設(shè)計(jì) 368.10建筑防火 368.11火災(zāi)檢測(cè)報(bào)警系統(tǒng) 378.12預(yù)期效果 37第九章項(xiàng)目實(shí)施進(jìn)度 389.1實(shí)施進(jìn)度計(jì)劃 389.2項(xiàng)目實(shí)施建議 38第十章項(xiàng)目招投標(biāo)方案 4010.1招標(biāo)原則 4010.2項(xiàng)目招標(biāo)范圍