2015建模A題太陽(yáng)影子定位

1、精選文檔A題 太陽(yáng)影子定位 一， 摘要（宋體小四號(hào)，簡(jiǎn)明扼要的詳細(xì)敘述，字?jǐn)?shù)不可以超過(guò)一頁(yè)，不要譯成英文）本文針對(duì)太陽(yáng)影子定位技術(shù)，通過(guò)太陽(yáng)與地球相對(duì)運(yùn)動(dòng)的規(guī)律，建立桿長(zhǎng)、影長(zhǎng)、經(jīng)緯度、時(shí)間、日期的關(guān)系，建立模型。綜合分析了不同地點(diǎn)，不同的時(shí)間，不同的季節(jié)時(shí)影子長(zhǎng)度的形成規(guī)律及變化趨勢(shì)，運(yùn)用了軟件進(jìn)行分析，得出不同地區(qū)影子變化的模型。最后將具體情況運(yùn)用到建立的模型中，對(duì)實(shí)際問(wèn)題進(jìn)行可行性分析，根據(jù)條件的改變完善對(duì)模型的應(yīng)用和實(shí)用性檢驗(yàn)。第一問(wèn)中，我們通過(guò)兩種太陽(yáng)高度角的表示方法建立等式關(guān)系，根據(jù)控制變量法，分析出影子長(zhǎng)度分別與經(jīng)、緯度、桿長(zhǎng)、時(shí)間、日期的關(guān)系。然后，根據(jù)時(shí)差計(jì)算關(guān)系，當(dāng)北京時(shí)間

2、在9:00-15:00時(shí)，天安門廣場(chǎng)的時(shí)間，并應(yīng)用建立的模型。第二問(wèn)中，首先根據(jù)影子坐標(biāo)求出影子的長(zhǎng)度，擬合北京時(shí)間與影子長(zhǎng)度的函數(shù)，找出影子長(zhǎng)度的最低的點(diǎn)，從而根據(jù)時(shí)間求出當(dāng)?shù)亟?jīng)度，由于誤差的存在，我們將經(jīng)度、桿長(zhǎng)、緯度給定一定范圍，根據(jù)第一問(wèn)公式進(jìn)行搜索，從而確定可能的地點(diǎn)。關(guān)鍵字：（宋體小四號(hào) ） 真太陽(yáng)時(shí) 平太陽(yáng)時(shí) 赤緯角 太陽(yáng)高度角 熵值法 二，問(wèn)題提出如何確定視頻的拍攝地點(diǎn)和拍攝日期是視頻數(shù)據(jù)分析的重要方面，太陽(yáng)影子定位技術(shù)就是通過(guò)分析視頻中物體的太陽(yáng)影子變化，確定視頻拍攝的地點(diǎn)和日期的一種方法。1.建立影子長(zhǎng)度變化的數(shù)學(xué)模型，分析影子長(zhǎng)度關(guān)于各個(gè)參數(shù)的變化規(guī)律，并應(yīng)用你們建立的模

3、型畫出2015年10月22日北京時(shí)間9:00-15:00之間天安門廣場(chǎng)（北緯39度54分26秒,東經(jīng)116度23分29秒）3米高的直桿的太陽(yáng)影子長(zhǎng)度的變化曲線。2.根據(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)。3. 根據(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)與日期。4附件4為一根直桿在太陽(yáng)下的影子變化的視頻，并且已通過(guò)某種方式估計(jì)出直桿的高度為2米。請(qǐng)建立確定視頻拍攝地點(diǎn)的

4、數(shù)學(xué)模型，并應(yīng)用你們的模型給出若干個(gè)可能的拍攝地點(diǎn)。如果拍攝日期未知，你能否根據(jù)視頻確定出拍攝地點(diǎn)與日期？ 三，問(wèn)題分析第一問(wèn)：根據(jù)物體在太陽(yáng)光照射下將產(chǎn)生影子的自然現(xiàn)象,研究物體影子的形成原理, 通過(guò)分析太陽(yáng)光線照射物體的角度的日變化和年變化，引起物體影子的長(zhǎng)度和朝向有規(guī)律地變化來(lái)建立數(shù)學(xué)模型。利用Matlab軟件繪出影子長(zhǎng)短隨時(shí)間變化的圖像。將問(wèn)題中所給參數(shù)帶入，解決問(wèn)題。由于太陽(yáng)光線照射物體的角度的日變化和年變化，引起物體影子的長(zhǎng)度和朝向有規(guī)律地變化第二問(wèn)：通過(guò)對(duì)附件所給的影子坐標(biāo)的數(shù)據(jù)，求出影子的長(zhǎng)度，然后1通過(guò)第一問(wèn)的相關(guān)公式，對(duì)影長(zhǎng)和時(shí)間的關(guān)系進(jìn)行擬合，得到一個(gè)二次方程，得出影長(zhǎng)的

5、最低值，從而可知正午時(shí)間，再算出經(jīng)度 四，建模過(guò)程第一問(wèn)1模型假設(shè)（1）：假設(shè)單一光源（太陽(yáng)光）照射（2）：直桿嚴(yán)格垂直于水平地面（3）：被照射直桿的形狀不會(huì)影響影子的長(zhǎng)度（4）：將整個(gè)天空視為一個(gè)天體圓（5）：不考慮大氣折射（6）：?jiǎn)栴}中給出的數(shù)據(jù)可靠2定義符號(hào)說(shuō)明 符號(hào)含義符號(hào)含義偏磁角（赤緯角）時(shí)差太陽(yáng)高度角入射角經(jīng)度緯度北京時(shí)間 太陽(yáng)方位角自1月1日算起的第幾天影子長(zhǎng)度時(shí)角桿長(zhǎng) 3模型建立：以桿影在陽(yáng)光下產(chǎn)生影子端點(diǎn)移動(dòng)的軌跡，代替太陽(yáng)運(yùn)行軌跡。運(yùn)用相對(duì)運(yùn)動(dòng)原理，將地球自轉(zhuǎn)及繞太陽(yáng)公轉(zhuǎn)的運(yùn)動(dòng)簡(jiǎn)化為地球不動(dòng)，太陽(yáng)繞地球轉(zhuǎn)動(dòng)。如圖一、二所示：圖一： 圖二：（1） 計(jì)算磁偏角（赤緯角）全年之

6、中，每一天太陽(yáng)和地球的運(yùn)轉(zhuǎn)與天體圓赤道之間所形成的夾角，也就是所謂的磁偏角都不同，會(huì)在與之間變化，其計(jì)算公式為： （2）由北京時(shí)間計(jì)算當(dāng)?shù)貢r(shí)間：按太陽(yáng)運(yùn)行位置，世界采取了時(shí)差制度并且遵循此制度，各國(guó)時(shí)間歷法都以此制度為基礎(chǔ)。按太陽(yáng)運(yùn)行位置，劃分時(shí)區(qū)，每個(gè)時(shí)區(qū)相差（每個(gè)時(shí)區(qū)相差1個(gè)小時(shí)）。當(dāng)?shù)貢r(shí)間的計(jì)算公式： 當(dāng)所得值為負(fù)數(shù)時(shí)，加上24小時(shí)。（3）計(jì)算時(shí)角因?yàn)榈厍蜃赞D(zhuǎn)一周約為24小時(shí)，所以，太陽(yáng)每小時(shí)大約自東向西移動(dòng)（即），故時(shí)角的計(jì)算公式為： 為正表示偏東，為負(fù)表示偏西。注意：計(jì)算中將其劃為弧度制。（4）計(jì)算太陽(yáng)高度角太陽(yáng)高度角簡(jiǎn)稱太陽(yáng)高度(其實(shí)是角度)。太陽(yáng)高度是決定地球表面獲得太陽(yáng)熱能數(shù)量

7、的最重要的因素，它在數(shù)值上等于太陽(yáng)在地球地平坐標(biāo)系中的地平高度。太陽(yáng)高度角的計(jì)算公式為： （5）利用太陽(yáng)高度角、桿長(zhǎng)及影長(zhǎng)列出函數(shù)式圖三：立竿測(cè)影模擬圖如圖三所示，由立竿見影的測(cè)量方式，得出影長(zhǎng)公式為： 4，模型求解：已知： 由（1）（2）（3）（4）（5）得：根據(jù)上式，圖四由圖四，在中由有幾何學(xué)原理，已知在時(shí)為遞增函數(shù)，故太陽(yáng)高度越小，影子越長(zhǎng)。影子有時(shí)比物體長(zhǎng)，有時(shí)比物體短。太陽(yáng)高度為45度時(shí)影子和物體一樣長(zhǎng)。由一天中太陽(yáng)位置的變化規(guī)律得出，早晚影子最長(zhǎng)，中午最短，早上到中午影子慢慢變短，中午到晚上影子慢慢又變長(zhǎng)。相似的，我們可以得出，早晚太陽(yáng)高度最小，中午最大，早上到中午太陽(yáng)高度慢慢變大

8、，中午到晚上太陽(yáng)高度慢慢又變小。（2）關(guān)于北京影長(zhǎng)問(wèn)題的探索利用Matlab繪出影子的變化規(guī)律圖。第一問(wèn)需要求解的題目中給出一下參數(shù)： 自1月1日算起的第295天。（2015年10月22日） 北京時(shí)間 9:00-15:00 東經(jīng)116度23分29秒 北緯39度54分26秒 3米注意：計(jì)算中將經(jīng)、緯度劃弧度制。將參數(shù)帶入Matlab中，繪出影子長(zhǎng)度變化曲線，如圖二所示：第二問(wèn)1. 模型假設(shè)：假設(shè)附件1所給頂點(diǎn)坐標(biāo)數(shù)據(jù)符合事實(shí)2. 定義符號(hào)說(shuō)明：方位角 ：高度角 ：赤緯角 ：物體所在地理緯度 ：當(dāng)?shù)貢r(shí)間 ：從1月1日起距當(dāng)?shù)厝掌诘奶鞌?shù) ：表示24小時(shí)制的時(shí)間數(shù)：太陽(yáng)某位置的方位時(shí)間：時(shí)角 ：影子的

9、長(zhǎng)度 ：影子頂點(diǎn)坐標(biāo)3,模型建立：因?yàn)橹窀拖鄬?duì)太陽(yáng)的位置所對(duì)應(yīng)于地球上所在點(diǎn)的相對(duì)位置, 由該點(diǎn)的地理緯度、日期和時(shí)間3個(gè)因素來(lái)決定。一般通過(guò)地平坐標(biāo)系及赤道坐標(biāo)系來(lái)同時(shí)表示太陽(yáng)的位置, 也就是以太陽(yáng)高度角 方位角 及 赤緯角、時(shí)角來(lái)表示。赤緯角是指地球赤道平面與太陽(yáng)和地球中心的連線之間的夾角，高度角指太陽(yáng)光的入射方向和地平面之間的夾角, 可知高度角的范圍是，方位角指經(jīng)過(guò)球心與太陽(yáng)位置點(diǎn)在地平圈上投影點(diǎn)的直線與地平圈正南向所夾的角。 定義方位角坐標(biāo)以正南向點(diǎn)為起始逆時(shí)針?lè)较驗(yàn)樨?fù), 分順時(shí)針?lè)较驗(yàn)檎? 亦分;正北向點(diǎn)為。因此根據(jù)太陽(yáng)位置的變化可以繪制出坐標(biāo)網(wǎng)圖, 在坐標(biāo)網(wǎng)圖中用同心圓來(lái)代表太陽(yáng)高

10、度圈, 用圓周上的刻度角來(lái)表示太陽(yáng)的方位角(自南向西為正值, 自南向東為負(fù)值)。故通過(guò)竿影軌跡點(diǎn)的坐標(biāo)可求出影長(zhǎng)，然后通過(guò)擬合影長(zhǎng)相關(guān)計(jì)算公式如下：影長(zhǎng)公式 (1)時(shí)角公式 (2) (3)太陽(yáng)高度角公式： (4) 太陽(yáng)方位角公式： (5) (6)由 可得(7) （8）根據(jù)(1)，求出胡各個(gè)時(shí)間對(duì)應(yīng)的影長(zhǎng)。如圖畫出的一段曲線圖像（如圖），并進(jìn)行二次擬合，求得p1= ,p2= , p3= ,根據(jù)韋達(dá)定理得：， ，為影子最短時(shí)的時(shí)間，可求出當(dāng)?shù)亟?jīng)度，即= ，由于二次擬合存在誤差， . 桿長(zhǎng)未知，假設(shè)1-3米之間，緯度范圍.4,模型求解：?jiǎn)栴}三1，模型假設(shè)：1，假設(shè)附件中有關(guān)葡萄和葡萄酒的數(shù)據(jù)符合事實(shí)

11、2，定義符號(hào)說(shuō)明： ： 表示葡萄酒的第i項(xiàng)指標(biāo)與釀酒葡萄的第j項(xiàng)指標(biāo)的相關(guān)系數(shù)（ 其中表示紅葡萄酒時(shí) ， ；表示白葡萄酒時(shí) ，） ；表示葡萄酒的第i項(xiàng)指標(biāo)與釀酒葡萄的第j項(xiàng)指標(biāo)顯著性實(shí)驗(yàn)的R-square值 ： 表示葡萄酒與釀酒葡萄第i項(xiàng)芳香物質(zhì)指標(biāo)的相關(guān)系數(shù)3，模型建立： 4，模型求解：第四問(wèn)1，模型假設(shè)：（1） 假設(shè)我們從前三問(wèn)中得到的數(shù)據(jù)準(zhǔn)確可靠（2） 假設(shè)附件1中品酒員所打分?jǐn)?shù)準(zhǔn)確可靠，可以反映葡萄酒的質(zhì)量2，定義符號(hào)說(shuō)明：（3） 表示第i瓶酒的第j個(gè)指標(biāo)無(wú)量綱化后的值（4） 表示第i種釀酒葡萄的第j個(gè)指標(biāo)無(wú)量綱化后的值（5） 表示第i瓶酒的綜合指標(biāo)3，模型建立：綜合指標(biāo)計(jì)算公式：每一

12、瓶酒對(duì)應(yīng)一個(gè)綜合指標(biāo)紅葡萄酒有27個(gè)綜合指標(biāo) （）白葡萄酒有28個(gè)綜合指標(biāo) （）用擬合方法找出綜合指標(biāo)的值與第二問(wèn)中葡萄酒的排名數(shù)的關(guān)系4，模型求解：將問(wèn)題三中的相關(guān)系數(shù)進(jìn)行篩選。將紅葡萄酒相關(guān)系數(shù)矩陣中絕對(duì)值小于0.6的記為-1表示沒(méi)有顯著線性關(guān)系；將矩陣中絕對(duì)值大于等于0.6且大于零的記為1表示正相關(guān)；將矩陣中絕對(duì)值大于等于0.6且小于零的記為0表示負(fù)相關(guān)。將白葡萄酒相關(guān)系數(shù)矩陣中絕對(duì)值小于0.4的記為-1表示沒(méi)有顯著線性關(guān)系；將矩陣中絕對(duì)值大于等于0.4且大于零的記為1表示正相關(guān)；將矩陣中絕對(duì)值大于等于0.4且小于零的記為0表示負(fù)相關(guān)。利用計(jì)算機(jī)編程求解出每瓶葡萄酒的綜合指標(biāo)（程序見附錄

13、）見下表：紅葡萄酒編號(hào)分?jǐn)?shù)綜合指標(biāo)白葡萄酒編號(hào)分?jǐn)?shù)綜合指標(biāo)168.62535178.25151273.62562276.755594375.1253337814471.6257092477.375577572.25181581.37532666.251489675.62511766.525407774.87530866.37560872.75735978.5132982.625630106881080.8752291162.3752131171.3755051268.75209581274.5301368.548981374.52841472.7591477.62563981566.25871

14、41579.125211669.625951666.875131774.75641780.753761864.875455321877231972.875881976.87574672076532762077.375152172.5162181.1253432271.8751322279.875662377.6251572377.25202471.6251272476.6251842567.25140122581.875292671.7515962675.87522572771.1259782777.87510212879.512利用matlab擬合綜合指標(biāo)的值與第二問(wèn)中葡萄酒的分?jǐn)?shù)得到下圖：

15、紅葡萄酒：去除一個(gè)奇點(diǎn)后用指數(shù)函數(shù)擬合得下圖：擬合結(jié)果：f(x) = a*exp(b*x) a = 6.06e+011 (-1.011e+013, 1.132e+013) b = -0.2746 (-0.5484, -0.0007818)R-square: 0.1055 白葡萄酒：用指數(shù)函數(shù)擬合后如下圖： 擬合結(jié)果：f(x) = a*exp(b*x) a = 1215 (-2.173e+004, 2.416e+004) b = -0.002948 (-0.2472, 0.2413) R-square: 0.000322 由R-square值可以看出兩組曲線擬合的結(jié)果不好，變換擬合函數(shù)嘗試數(shù)次后

16、所得擬合結(jié)果均不理想，因此我們認(rèn)為不能定量的用葡萄和葡萄酒的理化指標(biāo)來(lái)評(píng)價(jià)葡萄酒的質(zhì)量，只能根據(jù)圖像大致猜測(cè)綜合指標(biāo)與葡萄酒的質(zhì)量負(fù)相關(guān)五，模型評(píng)價(jià)與改進(jìn)第一問(wèn)中，我們所用的F檢驗(yàn)沒(méi)有T檢驗(yàn)合適，但是結(jié)果偏差不大第二問(wèn)中，topsis法靈活簡(jiǎn)便，操作方便。并用熵值法確定權(quán)重，具有合理性，但是缺少相應(yīng)的模型檢驗(yàn)。第四問(wèn)中，通過(guò)經(jīng)驗(yàn)設(shè)定綜合指標(biāo)進(jìn)行求解，簡(jiǎn)化了相應(yīng)的數(shù)學(xué)模型，只是缺少對(duì)綜合指標(biāo)設(shè)立的檢驗(yàn)，依據(jù)性不強(qiáng)。六，參考文獻(xiàn)1 Topsis綜合評(píng)價(jià)法，作者（無(wú)）2數(shù)學(xué)建模成績(jī)的評(píng)價(jià)與預(yù)測(cè)，作者（無(wú)）3陳光亭 裘哲勇 數(shù)學(xué)建模 高等教育出版社 2010年2月4王宏洲 數(shù)學(xué)建模優(yōu)秀論文 清華大學(xué)出

17、版社 2011年9月附錄：第二問(wèn)紅葡萄酒Ci（C語(yǔ)言）#include <stdio.h>#include <math.h>void main() int i,j,k;float s,t,p,q,z,e, b31,c31,d31,m27,n27,x272,g2731,l31,w31; float a2731=0.85 ,1.00 ,0.51 ,1.00 ,0.32 ,3.61 ,1.66 ,1.42 ,0.32 ,0.80 ,1.61 ,0.63 ,1.15 ,0.67 ,0.50 ,1.02 ,1.05 ,1.05 ,1.37 ,0.87 ,1.19 ,0.45 ,0

18、.76 ,0.70 ,1.14 ,1.17 ,0.57 ,0.92 ,0.24 ,0.24 ,0.97 ,0.89 ,1.13 ,0.13 ,1.13 ,1.55 ,0.94 ,0.70 ,1.30 ,0.47 ,0.74 ,1.83 ,0.59 ,1.68 ,1.02 ,0.77 ,1.00 ,1.01 ,1.06 ,0.68 ,0.77 ,1.36 ,0.77 ,0.34 ,0.55 ,0.97 ,1.15 ,0.85 ,0.99 ,0.87 ,0.24 ,1.04 ,3.52 ,1.05 ,0.64 ,1.05 ,1.26 ,0.59 ,0.95 ,0.81 ,1.35 ,0.84 ,1

19、.48 ,0.68 ,1.31 ,0.99 ,4.66 ,1.26 ,1.21 ,1.19 ,0.68 ,1.07 ,1.11 ,0.94 ,0.35 ,0.59 ,1.42 ,1.07 ,0.89 ,0.97 ,0.78 ,0.38 ,1.06 ,0.90 ,0.95 ,0.20 ,0.95 ,0.59 ,1.04 ,0.50 ,0.65 ,4.92 ,1.29 ,0.73 ,1.61 ,0.55 ,0.71 ,0.76 ,0.93 ,1.05 ,0.94 ,1.03 ,1.06 ,0.88 ,1.18 ,0.58 ,0.98 ,0.83 ,0.79 ,0.91 ,0.99 ,1.01 ,0

20、.35 ,1.01 ,0.77 ,1.05 ,0.08 ,1.05 ,1.48 ,0.75 ,1.31 ,1.32 ,2.50 ,0.87 ,1.20 ,0.96 ,1.25 ,0.13 ,0.19 ,1.03 ,0.86 ,0.99 ,1.03 ,0.99 ,0.99 ,1.06 ,2.15 ,1.43 ,0.76 ,0.98 ,1.41 ,1.00 ,0.55 ,0.38 ,1.02 ,1.44 ,0.97 ,0.15 ,0.97 ,0.44 ,0.44 ,0.00 ,1.54 ,3.09 ,1.25 ,0.73 ,0.92 ,0.83 ,0.46 ,0.22 ,1.20 ,0.99 ,1

21、.14 ,1.03 ,1.39 ,0.81 ,1.28 ,0.84 ,0.97 ,0.67 ,1.07 ,1.01 ,0.96 ,0.90 ,0.40 ,0.94 ,1.00 ,0.88 ,0.26 ,0.88 ,0.91 ,1.53 ,0.49 ,1.07 ,0.83 ,1.95 ,0.63 ,2.47 ,0.42 ,0.13 ,0.28 ,1.03 ,1.34 ,0.98 ,0.68 ,1.21 ,0.80 ,1.30 ,0.27 ,0.95 ,1.21 ,1.06 ,0.74 ,0.98 ,1.11 ,0.30 ,0.94 ,0.82 ,1.01 ,0.36 ,1.01 ,0.89 ,2

22、.68 ,2.28 ,2.12 ,0.27 ,0.83 ,1.04 ,0.62 ,1.03 ,1.24 ,3.26 ,0.97 ,0.87 ,1.05 ,0.34 ,0.96 ,1.08 ,0.97 ,0.89 ,1.02 ,1.63 ,0.89 ,1.35 ,1.02 ,0.60 ,0.41 ,0.94 ,0.95 ,1.26 ,1.03 ,1.26 ,2.06 ,0.82 ,1.00 ,0.71 ,0.85 ,0.51 ,2.05 ,0.57 ,2.50 ,1.02 ,1.65 ,0.95 ,0.86 ,0.94 ,1.03 ,0.88 ,1.07 ,0.98 ,0.78 ,0.78 ,1

23、.35 ,1.16 ,0.68 ,0.91 ,0.83 ,0.70 ,1.11 ,0.57 ,0.98 ,20.70 ,0.98 ,0.38 ,0.46 ,0.22 ,0.44 ,0.78 ,1.05 ,0.64 ,0.83 ,0.56 ,2.56 ,0.81 ,0.82 ,0.71 ,0.84 ,1.03 ,0.99 ,0.84 ,1.24 ,1.07 ,1.13 ,1.16 ,1.07 ,1.04 ,1.03 ,1.50 ,0.48 ,0.96 ,0.99 ,0.98 ,0.15 ,0.98 ,1.45 ,1.71 ,1.73 ,0.60 ,3.92 ,1.23 ,0.41 ,3.06 ,

24、0.31 ,5.59 ,0.72 ,1.03 ,1.05 ,0.97 ,1.37 ,0.82 ,1.18 ,0.88 ,0.74 ,0.67 ,0.87 ,0.87 ,0.53 ,1.07 ,3.78 ,0.41 ,0.88 ,1.07 ,0.89 ,0.13 ,0.89 ,0.95 ,1.06 ,1.03 ,0.88 ,2.71 ,1.74 ,0.82 ,1.94 ,0.47 ,0.14 ,0.07 ,1.21 ,1.16 ,1.21 ,1.37 ,1.27 ,0.94 ,1.11 ,0.80 ,1.06 ,0.61 ,0.94 ,1.27 ,1.01 ,0.93 ,8.27 ,0.97 ,

25、0.59 ,1.09 ,0.03 ,1.09 ,0.67 ,0.16 ,1.04 ,1.18 ,3.59 ,0.78 ,0.98 ,1.41 ,0.89 ,2.26 ,1.15 ,0.97 ,0.94 ,0.94 ,0.68 ,0.65 ,0.73 ,1.42 ,0.67 ,0.83 ,1.19 ,1.01 ,0.83 ,1.05 ,1.19 ,3.38 ,0.97 ,0.52 ,1.08 ,0.12 ,1.08 ,0.89 ,0.82 ,1.48 ,1.74 ,0.36 ,0.95 ,1.00 ,1.00 ,0.95 ,1.31 ,3.80 ,0.94 ,1.13 ,0.90 ,1.03 ,

26、0.80 ,1.11 ,0.94 ,0.87 ,0.77 ,1.17 ,0.98 ,1.33 ,0.97 ,0.80 ,0.34 ,1.03 ,0.91 ,0.94 ,0.14 ,0.94 ,0.97 ,0.72 ,1.87 ,1.08 ,2.27 ,1.56 ,0.81 ,0.55 ,0.67 ,0.04 ,0.27 ,0.88 ,0.92 ,0.99 ,0.68 ,1.28 ,0.77 ,1.35 ,0.66 ,0.98 ,0.74 ,1.01 ,1.11 ,0.97 ,0.91 ,0.57 ,0.94 ,0.65 ,1.05 ,0.17 ,1.05 ,1.41 ,1.44 ,2.16 ,

27、0.57 ,0.68 ,1.45 ,0.76 ,1.38 ,1.12 ,0.95 ,0.23 ,1.00 ,0.84 ,0.95 ,1.03 ,0.73 ,1.29 ,0.81 ,0.50 ,0.62 ,0.96 ,1.07 ,0.70 ,1.00 ,0.77 ,0.27 ,0.99 ,0.72 ,0.99 ,0.11 ,0.99 ,0.92 ,1.01 ,0.80 ,0.72 ,2.77 ,0.96 ,1.04 ,0.88 ,1.06 ,0.15 ,1.24 ,1.04 ,1.20 ,1.10 ,1.37 ,1.29 ,0.85 ,1.23 ,1.86 ,1.49 ,0.71 ,1.06 ,

28、1.72 ,0.97 ,0.65 ,0.26 ,1.06 ,1.01 ,0.92 ,0.23 ,0.92 ,0.56 ,1.11 ,0.47 ,1.14 ,2.26 ,1.52 ,0.50 ,2.58 ,0.64 ,0.09 ,0.18 ,1.11 ,1.17 ,1.05 ,1.03 ,1.20 ,0.87 ,1.20 ,0.82 ,1.17 ,0.66 ,0.94 ,0.83 ,0.97 ,1.21 ,0.51 ,0.92 ,1.03 ,0.98 ,0.15 ,0.98 ,0.87 ,0.85 ,0.12 ,1.28 ,1.77 ,0.90 ,1.19 ,1.01 ,1.15 ,0.80 ,

29、0.88 ,1.01 ,0.97 ,0.99 ,1.37 ,1.01 ,0.97 ,1.07 ,0.72 ,0.95 ,1.60 ,1.00 ,0.85 ,1.04 ,1.05 ,0.40 ,1.03 ,0.95 ,1.01 ,0.05 ,1.01 ,0.55 ,0.18 ,0.40 ,0.52 ,4.00 ,1.22 ,0.86 ,1.71 ,0.99 ,0.32 ,0.27 ,0.95 ,1.00 ,0.97 ,0.68 ,0.77 ,1.25 ,0.84 ,1.28 ,1.88 ,0.80 ,0.89 ,1.21 ,1.03 ,1.09 ,0.47 ,1.08 ,2.66 ,1.01 ,

30、0.10 ,1.01 ,2.42 ,0.58 ,2.16 ,0.76 ,1.85 ,0.90 ,1.10 ,1.02 ,0.92 ,1.63 ,1.33 ,1.01 ,1.14 ,1.00 ,1.37 ,1.01 ,0.99 ,1.06 ,0.62 ,0.60 ,1.20 ,0.90 ,0.56 ,0.96 ,0.56 ,4.55 ,1.03 ,1.08 ,0.88 ,0.15 ,0.88 ,1.01 ,1.53 ,0.70 ,0.92 ,0.86 ,1.21 ,1.12 ,1.14 ,0.96 ,0.89 ,0.39 ,1.10 ,1.00 ,1.09 ,1.03 ,0.89 ,1.21 ,

31、0.86 ,0.44 ,0.65 ,0.84 ,0.85 ,0.77 ,0.99 ,0.73 ,0.54 ,1.02 ,1.00 ,0.98 ,0.20 ,0.98 ,0.64 ,1.03 ,0.35 ,0.69 ,0.84 ,0.60 ,2.02 ,0.57 ,2.96 ,2.08 ,1.26 ,1.02 ,0.94 ,0.97 ,1.03 ,1.03 ,0.93 ,1.12 ,1.16 ,1.23 ,0.97 ,1.15 ,1.22 ,1.02 ,0.82 ,0.29 ,1.10 ,0.69 ,0.95 ,0.07 ,0.95 ,1.30 ,0.91 ,1.54 ,0.63 ,2.48 ,

32、1.21 ,0.60 ,0.96 ,1.00 ,0.61 ,0.41 ,0.99 ,1.08 ,0.94 ,1.03 ,1.08 ,0.86 ,1.21 ,2.16 ,1.34 ,0.76 ,1.14 ,1.28 ,0.99 ,0.54 ,0.25 ,1.01 ,0.59 ,0.97 ,0.13 ,0.97 ,0.45 ,0.49 ,0.15 ,0.60 ,2.55 ,0.98 ,0.78 ,1.49 ,0.65 ,0.44 ,0.85 ,0.74 ,0.69 ,0.90 ,1.03 ,1.27 ,0.70 ,1.48 ,1.20 ,1.42 ,1.04 ,0.87 ,1.15 ,1.03 ,

33、1.11 ,0.29 ,0.95 ,0.36 ,1.06 ,0.84 ,1.06 ,1.11 ,0.28 ,0.74 ,1.34 ,4.36 ,1.08 ,0.50 ,3.68 ,0.41 ,0.43 ,0.39 ,0.85 ,0.87 ,0.91 ,1.03 ,0.68 ,1.32 ,0.79 ,3.31 ,1.38 ,0.85 ,1.02 ,1.20 ,1.07 ,0.88 ,1.26 ,1.02 ,0.47 ,0.95 ,0.18 ,0.95 ,0.97 ,0.28 ,1.14 ,0.97 ,0.61 ,1.29 ,0.61 ,1.35 ,0.57 ,0.33 ,0.45 ,0.96 ,

34、0.94 ,0.96 ,1.03 ,1.04 ,0.91 ,1.14 ,1.18 ,0.84 ,0.89 ,0.89 ,1.04 ,1.10 ,1.63 ,1.05 ,1.01 ;for(i=0;i<27;i+) for(j=0;j<31;j+) gij=aij; for(j=0;j<31;j+) /原數(shù)據(jù)處理 s=0; for (i=0;i<27;i+) s+=aij*aij; bj=sqrt(s);for(j=0;j<31;j+) /歸一化后矩陣 for(i=0;i<27;i+) aij=aij/bj;for(j=0;j<31;j+) /每列元素按

35、冒泡法從小到大排序 for(i=1;i<27;i+) for(k=0;k<26;k+) if(akj>ak+1j) t=akj; akj=ak+1j; ak+1j=t; for(j=0;j<31;j+) /Z+ cj=a26j;for(j=0;j<31;j+) /Z- dj=a0j;for(i=0;i<27;i+)for(j=0;j<31;j+) s=0; for (i=0;i<27;i+) s+=aij*aij; s/=27; lj=s;p=0;s=0; for (j=0;j<30;j+) s+=a26j; s/=30; for (j=0

36、;j<30;j+) p+=lj; p/=30; e=sqrt(s*s+p*p)/2);printf("%fn",e);for(j=0;j<31;j+) wj=sqrt(lj*lj+a26j*a26j)/2);s=0;for(j=0;j<30;j+) s+=wj;printf("%f",s);w30=s*w30/e;for(j=0;j<31;j+) printf("%ft",wj);for(i=0;i<27;i+) /Di+ p=0; for(j=0;j<31;j+) p+=(gij-cj)*(gij

37、-cj); mi=sqrt(p);for(i=0;i<27;i+) /Di- q=0; for(j=0;j<31;j+) q+=(gij-dj)*(gij-dj); ni=sqrt(q);for(i=0;i<27;i+) /定義數(shù)組 xi0=i;for(i=0;i<27;i+) /Ci xi1=ni/(mi+ni);for(i=1;i<27;i+) /排序Ci for(j=0;j<26;j+)if(xj1<xj+11) t=xj1;xj1=xj+11;xj+11=t; z=xj0; xj0=xj+10; xj+10=z; for(i=0;i<27

38、;i+)printf("%.0ftc%.0f=%fn",xi0+1,xi0+1,xi1);/從優(yōu)到劣輸出樣品序號(hào)第二問(wèn)白葡萄酒Ci原理同上（C語(yǔ)言）第三問(wèn)相關(guān)系數(shù)矩陣：第四問(wèn)紅葡萄酒綜合指標(biāo)（C語(yǔ)言）#include <stdio.h>#include <math.h>void main()float a279=3.69 ,1.52 ,1.59 ,1.64 ,0.68 ,1.60 ,0.060359724,0.319634703,0.175660992,1.96 ,1.52 ,1.53 ,2.72 ,1.01 ,2.06 ,0.347068416,0

39、.908675799,1.089279247,1.51 ,1.82 ,1.36 ,1.50 ,1.46 ,1.77 ,0.398909631,0.953742307,1.247736327,0.70 ,0.89 ,0.95 ,0.88 ,0.81 ,0.79 ,1.02952272,1.181854278,1.211064832,1.06 ,0.81 ,0.96 ,0.74 ,1.39 ,0.93 ,0.838708107,1.194361723,1.088826512,0.44 ,1.01 ,0.94 ,0.91 ,1.23 ,0.95 ,1.386083189,1.080603534,1.

40、067095255,0.34 ,0.55 ,0.62 ,0.56 ,0.51 ,0.50 ,1.435977316,0.969227715,1.451919594,3.48 ,1.66 ,1.62 ,1.58 ,0.28 ,1.55 ,0.209311948,0.771490967,0.66461427,1.47 ,1.78 ,1.81 ,2.02 ,1.07 ,1.73 ,0.344877942,0.915028787,1.095164795,0.53 ,0.77 ,0.69 ,0.64 ,0.90 ,0.61 ,1.389490593,1.15267024,0.362187613,0.04

41、 ,0.63 ,0.64 ,0.43 ,0.11 ,0.47 ,2.161024168,0.241016478,0.884643245,0.32 ,0.89 ,0.77 ,0.61 ,0.60 ,0.63 ,1.306495972,1.001588247,1.384914886,0.76 ,0.88 ,0.79 ,0.81 ,0.37 ,0.74 ,1.012242315,1.165971809,0.887359652,0.95 ,0.84 ,0.80 ,0.63 ,0.60 ,0.73 ,0.589480858,1.115147905,1.598152843,0.46 ,0.55 ,0.65

42、 ,0.37 ,0.25 ,0.31 ,1.288728795,1.148898154,0.864270192,0.65 ,0.66 ,0.65 ,0.54 ,0.32 ,0.52 ,1.228369071,1.180266031,0.82397682,0.89 ,1.26 ,0.98 ,1.00 ,0.46 ,1.39 ,1.002993648,1.112368473,1.137269105,0.27 ,0.61 ,0.69 ,0.72 ,0.48 ,0.62 ,1.416019666,1.086360929,1.020916335,0.75 ,0.82 ,0.82 ,0.79 ,2.51

43、,0.75 ,1.160951152,1.289060949,0.935802246,0.28 ,0.81 ,0.78 ,0.83 ,0.27 ,0.71 ,1.910093217,0.52392297,0.718489678,1.19 ,1.39 ,1.43 ,0.91 ,2.44 ,1.60 ,0.52327987,1.048243002,1.594078233,0.95 ,0.98 ,0.99 ,1.19 ,1.24 ,1.03 ,0.986930173,1.073059361,1.186164433,1.57 ,1.50 ,2.00 ,2.48 ,3.52 ,2.53 ,0.35534

44、3539,0.930315664,1.135005433,1.02 ,0.79 ,0.86 ,0.76 ,1.91 ,0.74 ,1.042665563,1.172523327,0.800434625,0.60 ,0.74 ,0.71 ,0.62 ,0.72 ,0.74 ,1.222771193,1.266229899,0.522002898,0.57 ,0.50 ,0.62 ,0.44 ,0.76 ,0.34 ,0.815343053,1.231884058,1.321079319,0.52 ,0.82 ,0.76 ,0.67 ,1.33 ,0.68 ,1.536739115,0.967440937,0.723469757,b2728=0.85 ,1.00 ,0.51 ,1.00 ,0.32 ,3.61 ,1.66 ,1.42 ,3.12 ,1.26 ,1.61 ,1.59 ,1.15 ,0.67 ,0.50 ,1.02 ,1.05 ,1.05 ,1.37 ,0.87 ,1.19 ,2.22 ,0.76 ,0.70 ,1.14 ,1.17 ,0.57 ,0.92 ,0.89 ,1.13 ,0.13