數(shù)學(xué)建模論文-基于統(tǒng)計分析的葡萄酒評價模型

基于統(tǒng)計分析的葡萄酒評價模型摘要隨著經(jīng)濟(jì)的快速發(fā)展，人們對葡萄酒的品質(zhì)要求越來越高。目前，葡萄酒的評價都是建立在評酒員感官評價的基礎(chǔ)上。但是，由于評酒員的評價尺度、評價位置和評價方向的差異，導(dǎo)致了評價結(jié)果的不穩(wěn)定性。因此，一個迫切的任務(wù)是如何從葡萄的理化指標(biāo)和葡萄酒的理化指標(biāo)來評價葡萄酒的質(zhì)量。首先，利用秩和檢驗對兩組評酒員的評價結(jié)果進(jìn)行了顯著性差異分析，得出了兩組評酒員對紅葡萄酒的評價結(jié)果具有顯著性差異，對白葡萄酒的評價結(jié)果沒有顯著性差異，同時，使用方差分析得出第二組評酒員的評價結(jié)果更可信。其次，運(yùn)用聚類分析法分析了釀酒葡萄的理化指標(biāo)和葡萄酒質(zhì)量的關(guān)系，把釀酒葡萄分成五個等級。再次，采用相關(guān)性分析方法和回歸分析法分析了釀酒葡萄和葡萄酒的理化指標(biāo)之間的聯(lián)系，找出了一些重要因子，并對這些重要因子進(jìn)行了回歸分析，獲得了重要因子之間的確定性關(guān)系。最后，利用主成分分析和回歸分析法找出了釀酒葡萄和葡萄酒的理化指標(biāo)與葡萄酒質(zhì)量的線性關(guān)系：紅葡萄：；白葡萄：，通過該關(guān)系式可以大致的預(yù)測出葡萄酒的質(zhì)量。關(guān)鍵詞：秩和檢驗；主成分分析法；聚類分析法；相關(guān)性分析法；回歸分析法一、問題重述與分析1.1問題重述隨著經(jīng)濟(jì)的快速發(fā)展，人們對葡萄酒的品質(zhì)要求越來越高。目前，葡萄酒的評價都是建立在評酒員感官評價的基礎(chǔ)上。但是，由于評酒員的評價尺度、評價位置和評價方向的差異，導(dǎo)致了評價結(jié)果的不穩(wěn)定性。因此，一個迫切的任務(wù)是如何準(zhǔn)確地評價葡萄酒的質(zhì)量。為此，我們將解決以下問題：問題一：判斷兩組評酒員的的評價結(jié)果有無顯著性差異，并找出更可信的一組結(jié)果；問題二：根據(jù)釀酒葡萄的理化指標(biāo)和葡萄酒的質(zhì)量對這些釀酒葡萄進(jìn)行分級；問題三：分析釀酒葡萄與葡萄酒的理化指標(biāo)之間的聯(lián)系；問題四：分析釀酒葡萄和葡萄酒的理化指標(biāo)對葡萄酒質(zhì)量的影響，論證能否用葡萄和葡萄酒的理化指標(biāo)來評價葡萄酒的質(zhì)量。1.2問題分析本題是一個葡萄酒的評價問題。目前，葡萄酒的評價都是建立在評酒員的感官評價的基礎(chǔ)上。感官評價主要包括外觀、口感、香氣和平衡四個方面，它們必然與釀酒葡萄和葡萄酒的理化指標(biāo)存在重要的關(guān)聯(lián)。首先，數(shù)據(jù)的處理。在經(jīng)過分析和驗證后，適當(dāng)修正題中的個別有誤數(shù)據(jù)后，利用有效數(shù)據(jù)進(jìn)行建模求解。另外，由于數(shù)據(jù)之間的單位有區(qū)別，為了便于比較和分析，我們對所以的數(shù)據(jù)都進(jìn)行了標(biāo)準(zhǔn)化。其次，模型建立和求解。針對問題一，利用秩和檢驗對兩組評酒員的評價結(jié)果進(jìn)行顯著性差異分析，并使用方差選擇出比較可信的評價結(jié)果。針對問題二，運(yùn)用聚類分析法分析釀酒葡萄的理化指標(biāo)和葡萄酒質(zhì)量的關(guān)系，以此評判釀酒葡萄的等級。針對問題三，采用相關(guān)性分析方法分析釀酒葡萄和葡萄酒的理化指標(biāo)之間的聯(lián)系，找出了一些重要因子，以重要因子進(jìn)行了回歸分析，獲得了重要因子之間的確定性關(guān)系。針對問題四，利用主成分分析和回歸分析法找出了釀酒葡萄和葡萄酒的理化指標(biāo)與葡萄酒質(zhì)量的一些線性性關(guān)系。最后，模型檢驗，利用釀酒葡萄和葡萄酒的理化指標(biāo)與葡萄酒質(zhì)量的線性性關(guān)系檢驗葡萄酒的質(zhì)量。二、符號說明符號意義第組第個評酒員對第種紅葡萄酒的總評分第組第個評酒員對第種白葡萄酒的總評分第組評酒員對第種紅葡萄酒的平均總評分第組評酒員對第種白葡萄酒的平均總評分第組評酒員對第種紅葡萄酒的方差第組評酒員對第種白葡萄酒的方差三、模型假設(shè)1）評酒員都具有較好的資質(zhì)水平；2）所有芳香物質(zhì)對葡萄酒的質(zhì)量具有正影響；3）僅僅考慮釀酒葡萄與葡萄酒理化指標(biāo)中的一級指標(biāo)；4）葡萄酒的釀制過程對葡萄酒的質(zhì)量沒有影響。四、模型建立與求解4.1問題一4.1.1分析與建模問題一要求我們分析兩組評酒員的評價結(jié)果是否具有顯著性差異，并判斷哪組數(shù)據(jù)更可靠。評酒員對葡萄酒的評價是根據(jù)外觀、口感、香氣和平衡四個方面進(jìn)行感官評價的。為了研究的簡單，我們僅以該四個方面的評價分?jǐn)?shù)之和作為葡萄酒的質(zhì)量。顯著性差異的分析法主要包括檢驗、檢驗、秩和檢驗和符號檢驗等等。由于本題中的樣本量較小，總方差未知，我們采用秩和檢驗對兩組評酒員的評價結(jié)果進(jìn)行顯著性差異的分析。秩和檢驗的具體步驟：第一步：將兩個樣本數(shù)據(jù)混合并由小到大進(jìn)行等級排列（最小的數(shù)據(jù)秩次編為1，最大的數(shù)據(jù)秩次編為）;第二步：把容量較小的樣本中各數(shù)據(jù)的等級相加，即秩和，用表示;第三步：把值與秩和檢驗表中某α\o"顯著性水平"顯著性水平下的臨界值相比較，如果，則兩樣本差異不顯著；如果或，則表明兩樣本差異顯著。關(guān)于可信性結(jié)果的判斷，由于我們假設(shè)評酒員都具有較高的水平，因此，評價結(jié)果波動小的一組的可靠性強(qiáng)一些，故我們使用方差分析判斷兩組結(jié)果的可信性。4.1.2模型求解為了解決兩組紅葡萄酒的品酒員的評價結(jié)果有無顯著性差異及哪個結(jié)果更可信，我們運(yùn)用公式，求得樣品的平均分?jǐn)?shù)，再根據(jù)公式，求得各組樣品分?jǐn)?shù)的方差，結(jié)果如表1和圖1所示。從中可以看出，兩組紅葡萄酒的評價結(jié)果無顯著性差異，兩組白葡萄酒的評價結(jié)果有顯著性差異；第一組樣品的方差波動較大，表現(xiàn)出了結(jié)果的不穩(wěn)定性，所以第二組的結(jié)果更可信。表1：紅葡萄酒中的數(shù)據(jù)處理樣品編號第一組第二組秩和檢驗177.923.2968.173.690.2410275.844.167414.60.0368375.6128.2474.627.640.0749476.937.89048581.523.6572.112.290.9093675.520.4566.319.010.0574774.237.9665.356.410.3058872.328.016658.60.1115980.495.642431079.863.3668.832.560.01871171.479.0461.634.240.04061272.4126.0468.322.610.00171373.942.0968.813.760.01691477.114.2972.620.840.81881578.448.6465.737.210.05331667.374.0169.918.090.02511780.334.66091876.727.2165.445.240.17221976.423.4472.649.640.06352076.645.0475.835.160.38322179.257.9672.231.960.42652279.448.2471.621.840.08782377.410.4477.122.290.00352476.134.6971.59.650.13822579.595.8568.239.3612674.392.617237.40.569727773271.518.450.0012表2第二組白葡萄酒的數(shù)據(jù)處理樣品編號第一組第二組秩和檢驗1828377.923.290.1503274.2180.9675.844.160.7618378.359.8175.6128.240.8793479.440.2476.937.890.3834571113.881.523.650.0166668.4146.4475.520.450.3062777.535.2574.237.960.2109871.4165.2472.328.010.5443972.983.4980.495.640.06331074.3191.4179.863.360.40481172.3159.4171.479.040.76191263.3104.2172.4126.040.04481365.9153.6973.942.090.22531472102.877.114.290.42491572.4118.4478.448.640267.374.010.19841778.8129.7680.334.6111873.1140.8976.727.210.64901972.241.7676.423.440.12962077.857.9676.645.040.54472176.4155.4479.257.960.87922271124.879.448.240.08832375.939.2977.410.440.59352473.3100.0176.134.690.82032577.130.4979.595.850.13002681.365.6174.392.610.14952764.8129.9677320.01102881.372.4179.622.840.2234圖1紅葡萄酒質(zhì)量的方差圖圖2白葡萄酒質(zhì)量的方差圖4.2問題二4.2.1分析與建模由問題一的結(jié)果可知，第二組評價結(jié)果比較可信，因此我們在接下來的分析中都采用第二組評價的數(shù)據(jù)。針對問題二，由于葡萄的理化指標(biāo)種類很多，我們運(yùn)用運(yùn)用聚類分析法，對這些釀酒葡萄進(jìn)行分類，然后加入葡萄酒質(zhì)量評分，進(jìn)行分級。聚類分析指將物理或抽象對象的集合分組成為由類似的對象組成的多個類的分析過程。它是一種重要的人類行為。聚類分析的目標(biāo)就是在相似的基礎(chǔ)上收集數(shù)據(jù)來分類。聚類源于很多領(lǐng)域，包括數(shù)學(xué)，計算機(jī)科學(xué)，統(tǒng)計學(xué)，生物學(xué)和經(jīng)濟(jì)學(xué)。在不同的應(yīng)用領(lǐng)域，很多聚類技術(shù)都得到了發(fā)展，這些技術(shù)方法被用作描述數(shù)據(jù)，衡量不同數(shù)據(jù)源間的相似性，以及把數(shù)據(jù)源分類到不同的簇中。從統(tǒng)計學(xué)的觀點(diǎn)看，聚類分析是通過數(shù)據(jù)建模簡化數(shù)據(jù)的一種方法。傳統(tǒng)的統(tǒng)計聚類分析方法包括系統(tǒng)聚類法、分解法、加入法、動態(tài)聚類法、有序樣品聚類、有重疊聚類和模糊聚類等。采用k-均值、k-中心點(diǎn)等算法的聚類分析工具已被加入到許多著名的統(tǒng)計分析軟件包中，如SPSS、SAS等聚類分析是一種探索性的分析，在分類的過程中，人們不必事先給出一個分類的標(biāo)準(zhǔn)，聚類分析能夠從樣本數(shù)據(jù)出發(fā)，自動進(jìn)行分類。聚類分析所使用方法的不同，常常會得到不同的結(jié)論。不同研究者對于同一組數(shù)據(jù)進(jìn)行聚類分析，所得到的聚類數(shù)未必一致依據(jù)研究對象（樣品或指標(biāo)）的特征，對其進(jìn)行分類的方法，減少研究對象的數(shù)目。各類事物缺乏可靠的歷史資料，無法確定共有多少類別，目的是將性質(zhì)相近事物歸入一類。各指標(biāo)之間具有一定的相關(guān)關(guān)系。4.2.2模型的求解為了根據(jù)葡萄的理化指標(biāo)和葡萄酒的質(zhì)量對葡萄進(jìn)行分級，由附表2可以求得每個重要指標(biāo)下的的每種含量的平均值。用SPSS軟件進(jìn)行聚類分析得到有關(guān)釀酒葡萄分類的樹形圖，如圖3和圖4，因此，我們可以得到釀酒葡萄的分級，紅葡萄分為五級，白葡萄分為五級，見表3和表4.葡萄樣品4

4

─┬─────┐葡萄樣品27

27

─┘

├─┐葡萄樣品15

15

───────┘

├─┐葡萄樣品22

22

─────────┘

│葡萄樣品12

12

─┬─┐

├───┐葡萄樣品18

18

─┘

├─────┐

│

│葡萄樣品6

6

───┘

├─┘

├─┐葡萄樣品7

7

─────────┘

│

│葡萄樣品17

17

───┬─┐

│

├───────┐葡萄樣品24

24

───┘

├─────────┘

│

│葡萄樣品5

5

─────┘

│

│葡萄樣品13

13

─────┬───────────┘

├─┐葡萄樣品19

19

─────┘

│

│葡萄樣品20

20

───────────┬─────────┐

│

│葡萄樣品26

26

───────────┘

├───┘

├───────────┐葡萄樣品25

25

─────────────────────┘

│

│葡萄樣品11

11

───────────────┬───┐

│

│葡萄樣品16

16

───────────────┘

├───────┘

├───────┐葡萄樣品21

21

───────────────────┘

│

│葡萄樣品2

2

─────────┬─────────┐

│

│葡萄樣品9

9

─────────┘

├───────────────┐

│

│葡萄樣品23

23

───────────────────┘

├───┘

├─┐葡萄樣品8

8

─────────────┬─────────────┐

│

│

│葡萄樣品14

14

─────────────┘

├───────┘

│

│葡萄樣品1

1

───────────────────────────┘

│

│葡萄樣品10

10

───────────────────────────────────────────────┘

│葡萄樣品3

3

─────────────────────────────────────────────────┘圖3紅葡萄的層次聚類分析的樹形圖葡萄樣品

5

─┬───────────┐

葡萄樣品

20

─┘

│

葡萄樣品

9

───────┬───┐

├─┐

葡萄樣品

28

───────┘

├─┤

│

葡萄樣品

23

───────┬───┘

│

│

葡萄樣品

26

───────┘

│

│

葡萄樣品

2

─────┬─────┐

│

│

葡萄樣品

19

─────┘

│

│

├───┐

葡萄樣品

10

─────┬─┐

├─┘

│

│

葡萄樣品

25

─────┘

├─┐

│

│

│

葡萄樣品

12

───────┘

├─┘

│

│

葡萄樣品

14

─────────┘

│

├───┐

葡萄樣品

4

─────┬─┐

│

│

│

葡萄樣品

22

─────┘

├───────┘

│

│

葡萄樣品

17

───────┘

│

├─┐

葡萄樣品

21

───────────────────┘

│

│

葡萄樣品

8

─────┬───────┐

│

├─┐

葡萄樣品

11

─────┘

├─────────┘

│

│

葡萄樣品

16

─────────────┘

│

│

葡萄樣品

3

─────────────────────────┘

├───────┐

葡萄樣品

6

───┬───┐

│

│

葡萄樣品

18

───┘

├─────┐

│

│

葡萄樣品

7

───────┘

├───────────┐

│

├─────────────┐

葡萄樣品

15

─────────────┘

├─┘

│

│

葡萄樣品

24

─────────────────────────┘

│

│

葡萄樣品

1

─────────────────────┬─────────────┘

│

葡萄樣品

13

─────────────────────┘

│

葡萄樣品

27

─────────────────────────────────────────────────┘圖4紅葡萄的層次聚類分析的樹形圖表3紅葡萄質(zhì)量分級表紅葡萄分級酒樣品第一級樣品2、9、23、8、14、1第二級樣品3、10第三級樣品13、14、20、26、25、11、16、21第四級樣品4、7、15、22第五級樣品12、18、6、7、17、24、5表4白葡萄質(zhì)量分級表白葡萄分級酒樣品第一級樣品9、20、5第二級樣品21、17、22、4、14、12、25、10、19、2、26、23第三級樣品27、13、1第四級樣品3、16、11、8第五級樣品24、15、7、18、64.3問題三4.3.1分析與建模我們將使用回歸分析法和相關(guān)性分析法分析了釀酒葡萄與葡萄酒的理化指標(biāo)之間的聯(lián)系。回歸分析（regressionanalysis)是確定兩種或兩種以上變數(shù)間相互依賴的定量關(guān)系的一種統(tǒng)計分析方法。運(yùn)用十分廣泛，回歸分析按照涉及的自變量的多少，可分為一元回歸分析和多元回歸分析；按照自變量和因變量之間的關(guān)系類型，可分為線性回歸分析和非線性回歸分析。如果在回歸分析中，只包括一個自變量和一個因變量，且二者的關(guān)系可用一條直線近似表示，這種回歸分析稱為一元線性回歸分析。如果回歸分析中包括兩個或兩個以上的自變量，且因變量和自變量之間是線性關(guān)系，則稱為多元線性回歸分析.它基于觀測數(shù)據(jù)建立變量間適當(dāng)?shù)囊蕾囮P(guān)系，以分析數(shù)據(jù)內(nèi)在規(guī)律，并可用于預(yù)報，控制等問題。相關(guān)分析（correlationanalysis），相關(guān)分析是研究現(xiàn)象之間是否存在某種依存關(guān)系，并對具體有依存關(guān)系的現(xiàn)象探討其相關(guān)方向以及相關(guān)程度，是研究隨機(jī)變量之間的相關(guān)關(guān)系的一種統(tǒng)計方法相關(guān)關(guān)系是一種非確定性的關(guān)系，例如，以和分別記一個人的身高和體重，或分別記每公頃施肥量與每公頃小麥產(chǎn)量，則與顯然有關(guān)系，而又沒有確切到可由其中的一個去精確地決定另一個的程度，這就是相關(guān)關(guān)系。4.3.2模型的求解我們先從附表2中找出它們的公共理化指標(biāo)。其中紅葡萄與葡萄酒中都含有花色苷J(rèn)，單寧J，總酚J，酒總黃酮J，白藜蘆醇J。白葡萄與葡萄酒中都含有單寧J，總酚J，酒總黃酮J，白藜蘆醇J。先考慮紅葡萄，運(yùn)用SPSS軟件分別在葡萄中公共指標(biāo)與葡萄酒中五個指標(biāo)之間進(jìn)行回歸分析，得到五個圖表及多元一次線性表達(dá)式。再通過Excel進(jìn)行數(shù)據(jù)分析得到紅葡萄中某個指標(biāo)與葡萄酒五個指標(biāo)的相關(guān)系數(shù)。綜合考慮用SPSS軟件及Excel得到的相關(guān)系數(shù)，找出指標(biāo)間相關(guān)系數(shù)最高的，即可得出葡萄與葡萄酒理化指標(biāo)之間的聯(lián)系。同理，得到白葡萄與葡萄酒理化指標(biāo)之間的聯(lián)系。1，紅葡萄與葡萄酒理化指標(biāo)間的聯(lián)系（1）花色苷鮮重與葡萄酒五指標(biāo)的關(guān)系由此我們可以得到，葡萄中的花色苷鮮重與葡萄酒中的花色苷有很大的相關(guān)性，與單寧負(fù)相關(guān)性。（2）總酚與葡萄酒五指標(biāo)的關(guān)系葡萄中的總酚與葡萄酒中的酒總黃酮有一定的正相關(guān)性，與酒中花色苷有很大的負(fù)相關(guān)性（3）單寧與葡萄酒五指標(biāo)的關(guān)系葡萄中單寧與葡萄酒中的各指標(biāo)相關(guān)性不是很大(4)葡萄總黃酮與葡萄酒五指標(biāo)的關(guān)系葡萄中的總黃酮與葡萄酒中的總酚有很大的相關(guān)性，與葡萄酒中的花色苷和單寧負(fù)相關(guān)性（5）白藜蘆醇與葡萄酒五指標(biāo)的關(guān)系葡萄中的白藜蘆醇與葡萄酒中的總酚有很大的相關(guān)性，與葡萄酒中的花色苷有負(fù)相關(guān)性。利用相關(guān)性分析對釀酒葡萄與葡萄酒理化指標(biāo)的關(guān)系分析，見表5.花色苷單寧總酚酒總黃酮白藜蘆醇花色苷鮮重總酚單寧葡萄總黃酮白藜蘆醇花色苷J(rèn)1單寧J0.7444371總酚J0.7648940.9210231酒總黃酮J0.6644940.8373150.9039291白藜蘆醇J0.1238640.3312490.4859850.3994391花色苷鮮重0.922630.7195710.7735210.7085330.2003771總酚0.6131520.8171820.8750570.8830720.4592460.7277631單寧0.6608630.7179120.7431620.701010.3146730.6881330.7546371葡萄總黃酮0.4413510.6836910.8153750.8228910.5663940.5664830.8950880.688381白藜蘆醇-0.03490.0490280.0762260.0469140.013475-0.06002-0.01942-0.049910.0211341表5紅葡萄主要指標(biāo)相關(guān)性比重從表5中數(shù)據(jù)可以得出，紅葡萄中的花色苷鮮重與葡萄酒的葡萄酒中的花色苷和總酚有很大的正相關(guān)性，葡萄中的總酚與葡萄酒中的總酚和酒總黃酮有很大的相關(guān)性，葡萄中的單寧與葡萄酒中的總酚與單寧有比較大的相關(guān)性，葡萄中的總黃酮與葡萄酒中的總黃酮和總酚有很大的相關(guān)性，葡萄中的白藜蘆醇與葡萄酒中的各指標(biāo)相關(guān)性不是很明顯。將表5的結(jié)果與用SPSS作出的結(jié)果綜合分析可以得到紅葡萄與葡萄酒的理化指標(biāo)之間的聯(lián)系：葡萄中的花色苷與葡萄酒中的花色苷和葡萄酒總黃酮有較大正相關(guān)性；葡萄中的總酚與葡萄酒中的葡萄酒總黃酮有較大正相關(guān)性；葡萄中的單寧與葡萄酒中的花色苷和總酚有較大正相關(guān)性；葡萄中的葡萄總黃酮與葡萄酒中的總酚和葡萄酒總黃酮有較大正相關(guān)性；葡萄中的白藜蘆醇與葡萄酒中的總酚較大正相關(guān)性，與花色苷有很大的負(fù)相關(guān)性。2，白葡萄與葡萄酒理化指標(biāo)間的聯(lián)系類似紅葡萄與葡萄酒理化指標(biāo)間的聯(lián)系，我們得到白葡萄的相關(guān)數(shù)據(jù)：（1）總酚與葡萄酒五指標(biāo)的關(guān)系白葡萄中的總酚與葡萄酒中的總酚和酒總黃酮有很大的相關(guān)性，（2）單寧與葡萄酒五指標(biāo)的關(guān)系白葡萄中的單寧與葡萄酒中的總酚與單寧有比較大的相關(guān)性（3）葡萄總黃酮與葡萄酒五指標(biāo)的關(guān)系白葡萄中的總黃酮與葡萄酒中的總黃酮和總酚有很大的相關(guān)性（4）白藜蘆醇與葡萄酒五指標(biāo)的關(guān)系白葡萄中的白藜蘆醇與葡萄酒中的總酚有很大的相關(guān)性，與葡萄酒中的單寧與白藜蘆醇有很大的負(fù)相關(guān)性白葡萄單寧總酚酒總黃酮白藜蘆醇總酚單寧葡萄總黃酮白藜蘆醇單寧1.000總酚0.8781.000酒總黃酮0.4230.5641.000白藜蘆醇-0.162-0.138-0.0311.000總酚0.4280.5470.744-0.1361.000單寧0.5740.5730.346-0.0020.5141.000葡萄總黃酮0.4950.5880.697-0.1010.9430.5481.000白藜蘆醇-0.0620.037-0.095-0.2130.056-0.0560.0651.000表6白葡萄主要指標(biāo)相關(guān)性比重白葡萄中的總酚與葡萄酒中的總酚和酒總黃酮有一定的相關(guān)性，與白藜蘆醇有不太明顯的負(fù)相關(guān)性。白葡萄中的單寧與葡萄酒中的總酚與單寧有比較大的相關(guān)性，白葡萄中的總黃酮與葡萄酒中的總黃酮和總酚有很大的相關(guān)性，白葡萄中的白藜蘆醇與葡萄酒中的各指標(biāo)相關(guān)性不是很明顯。綜合考慮以上結(jié)論，得到白葡萄與葡萄酒理化指標(biāo)間的聯(lián)系如下：葡萄中的總酚與葡萄酒中的總酚和葡萄酒總黃酮有較大正相關(guān)性；葡萄中的單寧與葡萄酒中的總酚和單寧大正相關(guān)性；葡萄中的葡萄總黃酮與葡萄酒中的總酚和葡萄酒總黃酮有較大正相關(guān)性；葡萄中的白藜蘆醇與葡萄酒中的總酚較大正相關(guān)性，與白藜蘆醇有很大的負(fù)相關(guān)性。4.4問題四4.4.1分析與建模為了分析釀酒葡萄和葡萄酒的理化指標(biāo)對葡萄酒質(zhì)量的影響，我們首先運(yùn)用主成分分析法得到釀酒葡萄與葡萄酒的所有主成分；其次將這些主成分看作自變量，把葡萄酒的質(zhì)量看作因變量，運(yùn)用matlab對其進(jìn)行回歸分析即可得到一個多元一次方程，得到結(jié)果。4.4.2模型求解主成分分析是設(shè)法將原來眾多具有一定相關(guān)性（比如P個指標(biāo)），重新組合成一組新的互相無關(guān)的綜合指標(biāo)來代替原來的指標(biāo)。主成分分析，是考察多個變量間相關(guān)性一種多元統(tǒng)計方法，研究如何通過少數(shù)幾個主成分來揭示多個變量間的內(nèi)部結(jié)構(gòu)，即從原始變量中導(dǎo)出少數(shù)幾個主成分，使它們盡可能多地保留原始變量的信息，且彼此間互不相關(guān).通常數(shù)學(xué)上的處理就是將原來P個指標(biāo)作線性組合，作為新的綜合指標(biāo)。最經(jīng)典的做法就是用F1（選取的第一個線性組合，即第一個綜合指標(biāo)）的方差來表達(dá)，即Var(F1)越大，表示F1包含的信息越多。因此在所有的線性組合中選取的F1應(yīng)該是方差最大的，故稱F1為第一主成分。如果第一主成分不足以代表原來P個指標(biāo)的信息，再考慮選取F2即選第二個線性組合，為了有效地反映原來信息，F(xiàn)1已有的信息就不需要再出現(xiàn)在F2中，用數(shù)學(xué)語言表達(dá)就是要求Cov(F1,F2)=0，則稱F2為第二主成分，依此類推可以構(gòu)造出第三、第四，……，第P個主成分。我們從累積貢獻(xiàn)率分析并簡化評價葡萄的主要指標(biāo)[3]，先運(yùn)用SPSS進(jìn)行主成分分析得出七個主成份，累計貢獻(xiàn)率為以上，如表7和表8。表7紅葡萄評價因子的特征根和累積貢獻(xiàn)率主成分特征值方差貢獻(xiàn)率（%）累計方差貢獻(xiàn)率（%）16.60224.45124.45124.64517.20541.65732.88710.69352.34942.83910.51662.86651.9747.31170.17661.525.6375.80771.2154.580.307表8白葡萄評價因子的特征根和累積貢獻(xiàn)率主成分特征值方差貢獻(xiàn)率（%）累計方差貢獻(xiàn)率（%）16.60624.45124.45124.64517.20541.65732.88710.69352.34942.83910.51662.86651.9747.31170.17661.525.6375.80771.2154.580.307運(yùn)用matlab軟件進(jìn)行多元線性回歸分析結(jié)果如下：1.紅葡萄酒與釀酒葡萄的理化指標(biāo)對葡萄酒質(zhì)量影響關(guān)系從圖5可以看出：除第一、第十一和第十七個數(shù)據(jù)外，其余數(shù)據(jù)的殘差離零點(diǎn)均較近，且殘差的置信區(qū)間均包含零點(diǎn)，這說明回歸模型能較好的符合原始數(shù)據(jù)，而第一、第十一和第十七個數(shù)據(jù)可視為異常點(diǎn)。2.白葡萄酒與釀酒葡萄的理化指標(biāo)對葡萄酒質(zhì)量影響關(guān)系從殘差圖可以看出：除第五和第十六個數(shù)據(jù)外，其余數(shù)據(jù)的殘差離零點(diǎn)均較近，且殘差的置信區(qū)間均包含零點(diǎn)，這說明回歸模型能較好的符合原始數(shù)據(jù)，而第五個和第十六個數(shù)據(jù)可視為異常點(diǎn)。圖5：紅葡萄殘差圖圖6：白葡萄殘差圖由此就論證了能用葡萄和葡萄酒的理化指標(biāo)來評價葡萄酒的質(zhì)量。五、模型評價與檢驗本文主要利用秩和檢驗、聚類分析、相關(guān)性分析、主成分分析和回歸分析等統(tǒng)計方法研究了兩組評酒員的評價結(jié)果的顯著性差異、釀酒葡萄的等級、釀酒葡萄和葡萄酒的理化指標(biāo)之間的聯(lián)系和釀酒葡萄和葡萄酒的理化指標(biāo)與葡萄酒質(zhì)量的關(guān)系等問題，獲得了釀酒葡萄和葡萄酒的理化指標(biāo)與葡萄酒質(zhì)量的相關(guān)性方程，基本上能預(yù)測出葡萄酒的質(zhì)量。當(dāng)然，該模型也存在著一些缺點(diǎn)：比如，預(yù)測的結(jié)果有一定的擬合，有待進(jìn)一步研究。參考文獻(xiàn)[1]鄔彤等，數(shù)據(jù)統(tǒng)計分析—SPSS原理及其應(yīng)用[M]，北京:北京師范大學(xué)出版社，2001。[2]謝輝等，統(tǒng)計方法在葡萄理化指標(biāo)簡化中的應(yīng)用[J]，新疆農(nóng)業(yè)科學(xué)，48（8）：1434-1437，2011。[3]薛薇，統(tǒng)計分析與SPSS的應(yīng)用（第三版）[M]，北京：中國人民大學(xué)出版社，2011。[4]趙靜,但琦，數(shù)學(xué)建模與數(shù)學(xué)實驗（第3版）[M]，北京:高等教育出版社，2008。[5]李運(yùn)等，統(tǒng)計分析在葡萄酒質(zhì)量評價中的應(yīng)用[J]，釀酒科技，178:79-82,2009。附錄一紅葡萄酒質(zhì)量的方差圖的matlab編程x=[1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21;22;23;24;25;26;27];y=[23.29;44.16;128.24;37.89;23.65;20.45;37.96;28.01;95.64;63.36;79.04;126.04;42.09;14.29;48.64;74.01;34.61;27.21;23.44;45.04;57.96;48.24;10.44;34.69;95.85;92.61;32];plot(x,y,'b')gridonholdonz=[73.69;14.6;27.64;37.16;12.29;19.01;56.41;58.6;23.16;32.56;34.24;22.61;13.76;20.84;37.21;18.09;8.25;45.24;49.64;35.16;31.96;21.84;22.29;9.65;39.36;37.4;18.45];plot(x,z,'r')gridonxlabel('編號')ylabel('葡萄酒質(zhì)量的方差')二白葡萄酒質(zhì)量的方差圖的matlab編程x=[1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21;22;23;24;25;26;27;28];y=[83;180.96;59.81;40.24;113.8;146.44;35.25;165.24;83.49;191.41;159.41;104.21;153.69;102.8;118.44;160.2;129.76;140.89;41.76;57.96;155.44;24.8;39.29;100.01;30.49;65.61;129.96;72.41];plot(x,y,'b')gridonholdonz=[23.29;44.16;128.24;37.89;23.65;20.45;37.96;28.01;95.64;63.36;79.04;126.04;42.09;14.29;48.64;74.01;34.61;27.21;23.44;45.04;57.96;48.24;10.44;34.69;95.85;92.61;32;22.84];plot(x,z,'r')gridonxlabel('編號')ylabel('葡萄酒質(zhì)量的方差')三白葡萄和葡萄酒的理化指標(biāo)對葡萄酒質(zhì)量的matlab編程：紅葡萄：x1=[0.61;0.90;2.04;0.21;0.84;0.44;1.44;0.10;1.35;0.01;0.08;1.26;1.14;1.87;0.84;0.75;1.31;1.56;0.93;1.51;0.52;0.74;0.42;0.19;0.19;0.26;0.43;0.70];x2=[0.21;0.36;2.58;0.10;1.21;0.83;1.46;0.75;0.02;0.20;0.77;1.69;0.12;0.03;0.01;0.72;0.91;1.63;0.04;0.74;0.70;1.22;0.70;0.52;0.51;1.69;1.06;0.55];x3=[0.72;0.80;0.16;0.75;0.38;1.07;0.98;0.79;0.56;0.10;0.15;0.36;1.49;0.72;0.50;0.83;0.29;1.56;0.93;0.75;1.05;0.35;0.79;1.15;0.08;0.55;3.28;0.29];x4=[0.57;1.14;0.01;1.04;0.66;1.19;1.40;0.44;0.47;0.14;0.28;0.24;1.51;0.65;0.65;0.31;0.68;1.14;0.76;0.67;1.42;0.50;0.60;1.92;0.54;0.45;2.66;0.05];x5=[0.81;0.64;0.57;0.60;0.43;0.31;1.41;1.58;0.72;1.10;1.16;0.71;1.75;0.72;1.04;1.01;0.04;0.42;0.26;1.18;0.71;1.16;0.27;1.60;1.32;1.40;0.12;1.19];x6=[1.48；0.05；0.36；0.45；1.31；0.14；0.76；1.10；0.12；0.79；2.00；1.09；1.85；1.21；0.44；0.27；1.65；1.00；0.54；0.47；0.20；0.08；1.21；0.15；1.28；1.47；0.06；0.64];x7=[0.35；0.25；0.52；0.18；0.79；0.55；1.84；1.81；1.33；0.78；1.00；0.31；1.09；0.28；1.93；0.81；0.68；1.14；0.22；0.85；0.32；0.18；0.92；1.46；1.06；1.19；0.38；1.20];x8=[1.26；0.54；1.51；0.34；0.73；1.19；1.33；0.21；1.52；0.52；1.58；0.75；0.80；0.42；1.46；1.09；0.03；1.39；1.78；1.11；0.70；0.02；0.60；0.33；0.54；0.35；1.20；0.09];x9=[1.56；0.67；1.21；0.16；0.91；1.30；1.05；0.94；1.05；0.59；1.53；0.81；0.36；0.10；1.32；1.26；0.49；1.52；1.46；1.72；0.71；0.20；0.55；0.54；0.48；0.28；1.06；0.16];x10=[1.51；0.80；0.43；0.60；0.94；0.12；0.30；1.06；1.01；0.27；0.95；0.01；2.70；0.68；0.50；0.96；0.98；0.71；0.14；0.23；0.60；0.20；0.99；0.64；1.70；1.68；0.61；1.05];x11=[0.32；0.86；0.22；0.23；0.36；0.78；0.66；0.47；0.01；0.29；0.60；0.63；0.47；0.74；0.94；0.79；0.42；0.72；0.16；1.15；0.90；0.06；0.72；3.64；0.48；0.39；2.12；0.25];；x12=[0.37；0.67；0.69；0.06；0.15；0.53；0.48；1.87；0.32；0.25；0.34；1.03；0.19；0.26；0.67；0.28；0.36；0.21；0.21；0.27；0.81；0.15；0.65；3.76；0.01；0.38；2.06；0.17];x13=[0.69；0.74；0.25；0.07；0.48；0.44；0.08；0.53；0.50；0.33；1.05；0.95；0.26；0.19；0.67；0.80；2.70；0.56；0.60；0.00；0.28；0.16；0.57；3.09；0.83；0.36；1.69；0.01];X=[ones(28,1),x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13];Y=[0.43；0.23；0.29；0.12；1.57；0.33；0.74；1.33；1.22；1.03；1.62；1.30；0.83；0.18；0.59；2.91；1.19；0.05；0.04；0.02；0.84；0.90；0.27；0.14；0.94；0.70；0.15；0.97];[b,bint,r,rint,stats]=regress(Y,X);b;bint;statsrcoplot(r,rint)白葡萄x1=[1.341.84;1.05;0.60;0.44;0.61;0.83;0.08;2.32;0.79;1.30;0.40;0.05;0.01;0.42;0.53;0.10;1.11;0.41;0.31;0.22;0.26;2.26;0.90;0.48;1.11;0.88];x2=[3.38;1.33;0.59;0.29;0.17;0.66;0.50;1.52;1.51;0.68;1.09;0.82;0.45;0.39;0.59;0.50;0.51;0.73;0.12;0.91;0.18;0.35;0.75;0.44;0.62;0.52;0.79;];x3=[0.78;1.05;0.60;0.66;0.51;0.57;1.47;0.60;2.93;0.12;0.57;1.29;0.87;0.15;1.11;0.93;0.15;1.02;0.33;0.57;0.33;0.57;2.03;0.57;0.06;0.21;0.75];x4=[1.23;1.43;0.98;0.79;0.09;0.19;1.25;1.30;1.58;0.42;1.41;1.01;0.61;0.01;1.74;0.57;0.28;1.28;0.03;0.87;0.04;0.26;1.57;0.08;0.69;1.53;0.54;];x5=[0.58;0.61;1.83;0.31;0.16;1.18;0.79;0.46;0.29;1.94;0.02;1.27;0.46;0.17;0.22;0.86;0.73;0.75;0.42;0.63;1.00;0.48;0.44;0.54;2.41;1.91;0.47];x6=[0.540.66;2.16;0.67;0.17;1.56;0.25;0.54;0.67;1.82;0.31;2.34;0.67;1.16;0.07;0.55;1.15;0.54;0.07;0.37;0.04;0.97;0.39;0.67;1.12;1.06;0.48;];x7=[0.18;0.04;2.26;0.62;0.24;1.75;0.25;0.23;0.45;1.60;0.24;1.89;0.27;0.54;1.08;0.00;0.38;0.95;0.07;0.47;0.07;0.87;0.16;0.10;2.33;1.33;0.32];x8=[0.24;1.83;1.67;0.87;0.56;0.87;1.32;2.38;0.97;0.61;0.12;0.29;1.46;0.46;1.27;0.82;0.29;0.95;0.28;1.26;0.24;0.61;0.46;0.44;0.50;0.73;0.54];x9=[1.47;1.30;0.90;0.11;0.09;0.16;0.95;1.53;2.00;0.76;0.890.570.530.500.870.880.040.760.440.561.060.032.480.340.730.940.61];x10=[1.06;1.86;1.35;0.37;0.13;0.10;0.87;0.96;1.27;0.68;0.93;0.65;0.45;0.48;1.22;0.83;0.68;0.67;0.45;0.52;1.05;0.06;2.69;0.46;0.46;1.16;0.57];x11=[1.30;1.31;2.06;0.27;0.49;0.03;1.12;1.64;1.95;0.59;0.92;0.28;0.30;0.41;1.13;0.84;0.66;0.97;0.44;0.48;0.97;0.06;1.25;0.52;0.64;1.26;0.45];X=[ones(27,1),x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11];Y=[0.61;0.88;1.03;0.17;0.40;1.06;1.31;1.13;1.93;0.43;2.24;0.56;0.43;0.52;1.21;0.15;1.00;1.29;0.52;1.33;0.42;0.27;1.66;0.25;0.58;0.37;0.25;];[b,bint,r,rint,stats]=regress(Y,X);rcoplot(r,rint)百科1．正文規(guī)范

為了做到層次分明、脈絡(luò)清晰，常常將正文部分分成幾個大的段落。這些段落即所謂邏輯段，一個邏輯段可包含幾個小邏輯段，一個小邏輯段可包含一個或幾個自然段，使正文形成若干層次。論文的層次不宜過多，一般不超過五級。1．參考文獻(xiàn)的規(guī)范及其作用為了反映文章的科學(xué)依據(jù)、作者尊重他人研究成果的嚴(yán)肅態(tài)度以及向讀者提供有關(guān)信息的出處，正文之后一般應(yīng)列出參考文獻(xiàn)表引文應(yīng)以原始文獻(xiàn)和第一手資料為原則。所有引用別人的觀點(diǎn)或文字，無論曾否發(fā)表，無論是紙質(zhì)或電子版，都必須注明出處或加以注釋。凡轉(zhuǎn)引文獻(xiàn)資料，應(yīng)如實說明。對已有學(xué)術(shù)成果的介紹、評論、引用和注釋，應(yīng)力求客觀、公允、準(zhǔn)確。偽注、偽造、篡改文獻(xiàn)和數(shù)據(jù)等，均屬學(xué)術(shù)不端行為致謝一項科研成果或技術(shù)創(chuàng)新,往往不是獨(dú)自一人可以完成的,還需要各方面的人力,財力,物力的支持和幫助.因此,在許多論文的末尾都列有"致謝"。主要對論文完成期間得到的幫助表示感謝，這是學(xué)術(shù)界謙遜和有禮貌的一種表現(xiàn)。(1)文參考文獻(xiàn)不編序號,僅在文末按其重要程度或參考的先后順序排列.

(2)文后參考文獻(xiàn)不注頁碼.本文來自

(3)文后參考文獻(xiàn)的著錄項目及次序與注釋基本相同

1)著錄參考文獻(xiàn)可以反映論文作者的科學(xué)態(tài)度和論文具有真實、廣泛的科學(xué)依據(jù)，也反映出該論文的起點(diǎn)和深度。

2)著錄參考文獻(xiàn)能方便地把論文作者的成果與前人的成果區(qū)別開來。

3)著錄參考文獻(xiàn)能起索引作用。

4)著錄參考文獻(xiàn)有利于節(jié)省論文篇幅。

5)著錄參考文獻(xiàn)有助于科技情報人員進(jìn)行情報研究和文摘計量學(xué)研究。[序號]主要責(zé)任者．文獻(xiàn)題名[文獻(xiàn)類型標(biāo)識]．出版地：出版者，出版年．

[1]周振甫．周易譯注[M]．北京:中華書局．1985．

[2]陳送．五四前后東西方文化問題論戰(zhàn)文選[C]．北京:中國社會科學(xué)出版社,1985．

[3]陳桐生．中國史官文化與《史記》[D]．西安:陜西師范大學(xué)文學(xué)研究所,1992年．

[4]白永秀,劉敢,任保平．西安金融、人才、技術(shù)三大要素市場培育與發(fā)展研究[R]．西安:陜西師范大學(xué)西北經(jīng)濟(jì)研究中心,1998．[序號]主要責(zé)任者．文獻(xiàn)題名[J]．刊名，年，卷(期)．

[5]何齡修．讀顧城《南明史》[J]．中國史研究，1998(3)．[序號]析出文獻(xiàn)主要責(zé)任者．析出文獻(xiàn)題名[A]．原文獻(xiàn)主要責(zé)任者(任選)．原文獻(xiàn)題名[C]．出版地：出版者，出版年．

[6]瞿秋白．現(xiàn)代文明的問題與社會主義[A]．羅榮渠．從西化到現(xiàn)代化[C]．北京：北京大學(xué)出版社，1990．[序號]主要責(zé)任者．文獻(xiàn)題名[N]．報紙名，出版日期(版次)．

[7]謝希德．創(chuàng)造學(xué)習(xí)的新思路[N]．人民日報，1998-12-25(10)．

e．國際、國家標(biāo)準(zhǔn)

[序號]標(biāo)準(zhǔn)編號，標(biāo)準(zhǔn)名稱[S]．

英文（例子）：本文來自

[01]Brown,H.D.TeachingbyPrinciples:AnInteractiveApproachtoLanguagePedagogy[M].PrenticeHallRegents,1994.

[02]Brown,JSetal.SituatedCognitionandtheCultureofLearning[J].EducationalReasercher,1,1989.

[