釀酒葡萄與葡萄酒理化指標(biāo)及其質(zhì)量的分析與評(píng)估

1、釀酒葡萄與葡萄酒理化指標(biāo)及其質(zhì)量的分析與評(píng)估摘要 本文通過(guò)運(yùn)用spss首先對(duì)兩組品酒員對(duì)葡萄酒的評(píng)價(jià)進(jìn)行了差異顯著性分析以及評(píng)cronbach 系數(shù)分析對(duì)品酒員評(píng)價(jià)的可信度進(jìn)行了分析，得到了兩組品酒員評(píng)價(jià)之間存在顯著性差異，并且第一組品酒員的評(píng)價(jià)結(jié)果更為可信。接下來(lái)我們?nèi)〉谝唤M評(píng)酒員給每個(gè)樣品酒所給出的總分的平均分表示該樣品葡萄酒質(zhì)量的量化指標(biāo)進(jìn)行相關(guān)討論 我們通過(guò)對(duì)葡萄及葡萄酒的理化指標(biāo)歸一化處理，將它們之間的聯(lián)系現(xiàn)實(shí)的更為明顯，緊接著我們利用因子分析、主成分分析將決定葡萄酒質(zhì)量的因素維度較低，從而達(dá)到簡(jiǎn)化模型的目的，并通過(guò)分析將各成分因子重新命名，再將命名后的主成分和評(píng)論員對(duì)葡萄酒的評(píng)價(jià)作

2、為依據(jù)對(duì)釀酒葡萄等級(jí)評(píng)定，并引入層次分析法對(duì)評(píng)定體系進(jìn)行改進(jìn)構(gòu)想。最后利用matlab求解主成分因子與葡萄酒質(zhì)量之間的函數(shù)關(guān)系，并用bp人工神經(jīng)網(wǎng)絡(luò)進(jìn)行驗(yàn)證及討論。最后我們?cè)诿恳粋€(gè)問(wèn)題后面進(jìn)行了深入的反思與總結(jié)，得出了一些具體的改進(jìn)思路及方法，并得出在原有問(wèn)題上的修正。關(guān)鍵字 spss matlab cronbach 系數(shù)分析 主成分分析 層次分析法 bp人工神經(jīng)網(wǎng)絡(luò) 模型的改進(jìn)與修正 1. 問(wèn)題重述隨著經(jīng)濟(jì)的日益發(fā)展，人們生活水平的不斷提高，葡萄酒的受眾隨之增加，品酒行業(yè)逐漸專業(yè)化。品質(zhì)是現(xiàn)代葡萄酒生產(chǎn)追求的目標(biāo)之一，針對(duì)酒類的質(zhì)量檢測(cè)也成為食品行業(yè)工作的重中之重。根據(jù)已有數(shù)據(jù)，運(yùn)用數(shù)學(xué)建模

3、的方法，對(duì)葡萄酒進(jìn)行評(píng)價(jià)是一個(gè)重要問(wèn)題。確定葡萄酒質(zhì)量時(shí)一般是通過(guò)聘請(qǐng)一批有資質(zhì)的評(píng)酒員進(jìn)行品評(píng)。每個(gè)評(píng)酒員在對(duì)葡萄酒進(jìn)行品嘗后對(duì)其分類指標(biāo)打分，然后求和得到其總分，從而確定葡萄酒的質(zhì)量。釀酒葡萄的好壞與所釀葡萄酒的質(zhì)量有直接的關(guān)系，葡萄酒和釀酒葡萄檢測(cè)的理化指標(biāo)會(huì)在一定程度上反映葡萄酒和葡萄的質(zhì)量。從上述相關(guān)關(guān)系出發(fā)，參考相關(guān)數(shù)據(jù)，針對(duì)釀酒葡萄與葡萄酒的理化指標(biāo)和葡萄酒的質(zhì)量這幾個(gè)變量，進(jìn)行分析并建立數(shù)學(xué)模型，討論關(guān)于葡萄酒的評(píng)價(jià)問(wèn)題：特別注意指出模型中的優(yōu)點(diǎn)和不足之處，并做出改進(jìn)方向。2. 符號(hào)說(shuō)明符號(hào)說(shuō)明符號(hào)說(shuō)明第i個(gè)品酒員對(duì)第j款酒所給的總分s香氣指標(biāo)k整個(gè)品嘗的次數(shù)第i類芳香物質(zhì)（1

4、i8）總得分的方差q其他影響因素所有在第i位品酒師打分的方差k香氣指標(biāo)關(guān)于芳香物質(zhì)函數(shù)的常系數(shù)y葡萄酒的單寧m測(cè)得葡萄酒評(píng)分x1葡萄酒的總酚z葡萄酒標(biāo)準(zhǔn)分值（真實(shí)葡萄酒評(píng)分）x2葡萄酒的黃酮醇u外觀指標(biāo)第一組的總平均分v口感指標(biāo)理化指標(biāo)向量w整體評(píng)價(jià)指標(biāo)代表原來(lái)數(shù)據(jù)群體的因子香氣指標(biāo)占總分的權(quán)重兩類新的性質(zhì)因子外觀指標(biāo)占總分的權(quán)重口感類物質(zhì)成分因子口感指標(biāo)占總分的權(quán)重色澤類物質(zhì)成分因子整體評(píng)價(jià)指標(biāo)占總分的權(quán)重3. 模型建立與求解3.1. 數(shù)據(jù)處理3.1.2. 數(shù)據(jù)預(yù)處理 由于試題中給出的附件酒樣品未按照順序排列，我們首先將數(shù)據(jù)進(jìn)行了預(yù)處理，每組數(shù)據(jù)均按照樣品號(hào)和品酒員號(hào)順序排列，加和得到每一劑

5、量酒樣品的總評(píng)分?jǐn)?shù)，并求得各項(xiàng)明細(xì)指標(biāo)的平均值；其次將多組同類理化指標(biāo)取平均值，以便分析討論，進(jìn)行模型的建立與求解。3.1.3. 數(shù)據(jù)處理方法數(shù)據(jù)的預(yù)處理和部分圖示，采用辦公軟件microsoft excel 2010；數(shù)據(jù)的歸一化處理，采用辦公軟件microsoft excel 2010；數(shù)據(jù)結(jié)果的顯著差異分析，采用spss 19.0；數(shù)據(jù)結(jié)果的信度分析，采用spss 19.0；數(shù)據(jù)結(jié)果的主成分分析，采用spss 19.0；數(shù)據(jù)處理的擬合分析，采用matlab2012a3.2. 問(wèn)題1.1：附件1中兩組評(píng)酒員評(píng)價(jià)結(jié)果的顯著性差異分析要比較兩組評(píng)論員之間是否具有顯著性差異，而每組有10個(gè)評(píng)論員

6、，27或28個(gè)葡萄酒樣品，由于附件中的數(shù)據(jù)都是單項(xiàng)給出的，整體可比性不強(qiáng)，所以我們先用excel將數(shù)據(jù)進(jìn)行了初步的整理（求出每個(gè)品酒員對(duì)各樣品的總分，以及各個(gè)明細(xì)樣品的平均值），以第一組對(duì)紅葡萄酒樣品一的評(píng)價(jià)為例，處理結(jié)果如表1所示：表1. 紅葡萄酒樣品一的評(píng)價(jià)項(xiàng)目品酒員12345678910平均值葡萄酒樣品 1外觀分析澄清度544443432433.5色調(diào)1088666686666.6香氣分析純正度655545544544.6濃度867767774646.1質(zhì)量161412121414121412121212.8口感分析純正度655545644554.8濃度866766866666.3持久性

7、876767876666.6質(zhì)量221916191616191916161617.2平衡/整體評(píng)價(jià)11109109109999109.4總分1008478827579848169757277.93.2.1. 初步分析：針對(duì)處理后的數(shù)據(jù)，我們初步的想法是將每組中各個(gè)品酒員對(duì)每項(xiàng)樣品的總分作為一個(gè)樣本，因此每個(gè)樣本里將有270項(xiàng)數(shù)據(jù)。針對(duì)紅葡萄酒評(píng)價(jià)的兩組樣本，利用spss進(jìn)行一次相互獨(dú)立樣本之間的t檢驗(yàn)，探測(cè)兩個(gè)樣本是否具有顯著性差異。但隨后通過(guò)查閱統(tǒng)計(jì)分析中t檢驗(yàn)法的一些特性得到了相反的結(jié)論：鑒于本題樣本是每十個(gè)總分針對(duì)一個(gè)葡萄酒樣本，并不滿足獨(dú)立樣本t檢驗(yàn)方法的適用范圍，每組樣本中的數(shù)據(jù)是有

8、一定關(guān)聯(lián)的，因此我們采用了另外一種檢驗(yàn)方法配對(duì)樣本t檢驗(yàn)法。3.2.2. 二次分析 首先，我們對(duì)紅葡萄酒的評(píng)價(jià)進(jìn)行差異性分析（白葡萄酒的處理方法與其一致），將每項(xiàng)葡萄酒樣品各項(xiàng)明細(xì)評(píng)分的平均分作為一個(gè)樣本（包括總分的平均分），則每個(gè)樣本里共有1127=297項(xiàng)數(shù)據(jù)，兩組數(shù)據(jù)之間都是針對(duì)相同的樣品和指標(biāo)進(jìn)行評(píng)價(jià)，所以兩組中每項(xiàng)數(shù)據(jù)都是一一配對(duì)，互相關(guān)聯(lián)的。其次，我們對(duì)兩組葡萄酒的全部總分?jǐn)?shù)據(jù)共27*10*2=540組，進(jìn)行綜合處理并作圖1、圖2所示，可以得知兩組數(shù)據(jù)在數(shù)值上差異較大，初步判斷兩組數(shù)據(jù)具有顯著性差異，關(guān)于信度分析下文將給出分析討論。 圖1 圖2 再次，我們以兩組紅葡萄酒樣品1為例，

9、對(duì)其數(shù)據(jù)做了歸一化處理，利用雷達(dá)圖表示，如圖3所示。加強(qiáng)了上述關(guān)于顯著性差異的推論。圖3. 紅葡萄酒樣品歸一化雷達(dá)圖表示因此通過(guò)以上對(duì)數(shù)據(jù)和檢驗(yàn)方法的分析，我們采用配對(duì)t檢驗(yàn)法對(duì)兩組樣本進(jìn)行差異顯著性分析。3.2.3. 基本假設(shè)：1. 每組評(píng)論員對(duì)各項(xiàng)酒的評(píng)價(jià)分?jǐn)?shù)成正態(tài)分布。2. 評(píng)論員都具有一定的專業(yè)品酒技術(shù)，對(duì)每種酒的各個(gè)指標(biāo)的評(píng)價(jià)不會(huì)有大的偏差（保證評(píng)分的關(guān)聯(lián)度）。3.2.4. 模型的建立與求解：將處理后的數(shù)據(jù)導(dǎo)入到spss軟件中，針對(duì)兩組紅葡萄酒評(píng)價(jià)中的各項(xiàng)指標(biāo)評(píng)分的平均值進(jìn)行配對(duì)樣本t檢驗(yàn)法，設(shè)定置信度區(qū)間為95，得到結(jié)果如下：表2.成對(duì)樣本統(tǒng)計(jì)量均值n標(biāo)準(zhǔn)差均值的標(biāo)準(zhǔn)誤對(duì) 1紅葡萄

10、酒各項(xiàng)平均分113.285629719.397301.12555紅葡萄酒各項(xiàng)平均分212.820929718.634121.08126表3.成對(duì)樣本相關(guān)系數(shù)n相關(guān)系數(shù)sig.對(duì) 1紅葡萄酒各項(xiàng)平均分1 紅葡萄酒各項(xiàng)平均分2297.996.000表2為描述性統(tǒng)計(jì)表，給出了紅葡萄酒的各項(xiàng)平均分的平均值，標(biāo)準(zhǔn)差以及標(biāo)準(zhǔn)誤。表3為樣本間的配對(duì)性的測(cè)驗(yàn)結(jié)果，可以看到測(cè)出的p值（即sig值）等于00.001，因此在95%的置信水平上差異顯著，即兩組紅葡萄酒的各項(xiàng)平均分的平均值顯著相關(guān)，驗(yàn)證了我們之前的推論，因此也符合配對(duì)樣本t檢驗(yàn)的前提條件。表4.成對(duì)樣本檢驗(yàn)成對(duì)差分tdfsig.(雙側(cè))均值標(biāo)準(zhǔn)差均值的

11、標(biāo)準(zhǔn)誤差分的 95% 置信區(qū)間下限上限對(duì) 1紅葡萄酒各項(xiàng)平均分1 紅葡萄酒各項(xiàng)平均分2.464681.88444.10935.24949.679884.250296.000表4為最終的配對(duì)樣本t檢驗(yàn)結(jié)果報(bào)表。其中t表示用公式：所計(jì)算出的值，是最終計(jì)算出p值的依據(jù)。df為自由度，sig（雙側(cè)顯著性概率）則為我們所需要的p值?？梢钥吹阶罱K結(jié)果是p=0.0000.001，因此在95%的置信水平上差異顯著，所以兩組評(píng)論員評(píng)分無(wú)顯著性差異的假設(shè)不成立，即兩組評(píng)論員對(duì)紅酒的評(píng)分具有顯著性差異。圖四為以上結(jié)果的綜合圖形描繪。圖4. 結(jié)果的綜合圖形描繪對(duì)于白葡萄酒的組，利用同樣的方法進(jìn)行差異性分析，在spss

12、軟件里同樣也得到如下結(jié)果，如表5、表6、表7所示：表5.成對(duì)樣本統(tǒng)計(jì)量均值n標(biāo)準(zhǔn)差均值的標(biāo)準(zhǔn)誤對(duì) 1白葡萄酒各項(xiàng)平均分113.501930819.645421.11940白葡萄酒各項(xiàng)平均分213.914930820.210871.15162 表6.成對(duì)樣本相關(guān)系數(shù)n相關(guān)系數(shù)sig.對(duì) 1白葡萄酒各項(xiàng)平均分1 白葡萄酒各項(xiàng)平均分2308.996.000 表7.成對(duì)樣本檢驗(yàn)成對(duì)差分tdfsig.(雙側(cè))均值標(biāo)準(zhǔn)差均值的標(biāo)準(zhǔn)誤差分的 95% 置信區(qū)間下限上限對(duì) 1白葡萄酒各項(xiàng)平均分1 白葡萄酒各項(xiàng)平均分2-.412991.93277.11013-.62969-.19628-3.750307.000

13、圖5.白葡萄酒成對(duì)樣本檢驗(yàn) 如圖5所示，此時(shí)p值仍小于0.001，因此兩組評(píng)論員對(duì)白葡萄酒的評(píng)價(jià)在95%的置信水平上也具有差異顯著性，從而可得出結(jié)論：無(wú)論是針對(duì)紅葡萄酒還是白葡萄酒，兩組評(píng)酒員的評(píng)價(jià)結(jié)果都有顯著性差異。3.2.5. 問(wèn)題反思： 對(duì)于該問(wèn)題，我們將每項(xiàng)葡萄酒評(píng)分做了平均化處理，再把處理出來(lái)的平均分作為樣本，這樣每個(gè)樣本中的數(shù)據(jù)變成了與每個(gè)指標(biāo)得分情況對(duì)應(yīng)的一元數(shù)值，然后用配對(duì)樣本t檢驗(yàn)法對(duì)其進(jìn)行差異顯著性分析。通過(guò)這種方法處理出來(lái)的結(jié)果具有一定的科學(xué)性與準(zhǔn)確性。但是由于數(shù)據(jù)是進(jìn)行了平均化處理后的數(shù)據(jù)，因此模型所反映出的情況不具有完整性，對(duì)此我們便有了下面這種改進(jìn)思路：將每組中10

14、名品酒師對(duì)27項(xiàng)樣品的評(píng)價(jià)總分作為一個(gè)的矩陣（其中第i個(gè)品酒員對(duì)第j款酒所給的總分）。每個(gè)矩陣分別對(duì)應(yīng)該組品酒師的評(píng)分具體情況，再用這兩組數(shù)據(jù)構(gòu)成的矩陣進(jìn)行多元數(shù)據(jù)的t檢驗(yàn)，如此處理出的結(jié)果就應(yīng)該更具有完整性與可靠性。但由于本組成員所了解的統(tǒng)計(jì)學(xué)知識(shí)還不足以去運(yùn)用這樣多元性的、更為精確的處理手段，因此就將這作為一個(gè)以后改進(jìn)的思路。3.3. 問(wèn)題1.2：附件1中兩組評(píng)酒員評(píng)價(jià)結(jié)果的信度分析我們針對(duì)兩組中每位品酒師對(duì)每種酒所給出的總分得出的矩陣對(duì)其中的元素逐行進(jìn)行可信度分析。3.3.1模型討論首先我們根據(jù)品酒師打分為主觀作用，我們選擇評(píng)價(jià)評(píng)分者信度的kendall w協(xié)同檢驗(yàn)?zāi)Ｐ停ㄟ^(guò)spss對(duì)紅

15、白葡萄4組數(shù)據(jù)進(jìn)行模擬分析，得到的漸進(jìn)顯著性系數(shù)p均小于0.05（等于0.00），協(xié)同系數(shù)均為0.50左右，無(wú)法得出明確的結(jié)論。經(jīng)過(guò)思考，我們發(fā)現(xiàn)kendall模型對(duì)于采用等級(jí)評(píng)定方式評(píng)分的情況是可取的，然而當(dāng)評(píng)分者用其他非等級(jí)評(píng)定的方式（如百分制）評(píng)定時(shí)，其信度系數(shù)是不合適的，因?yàn)槠湫哦认禂?shù)僅能反映評(píng)分者之間的相對(duì)一致性，并沒(méi)有考慮評(píng)分者之間存在的評(píng)分絕對(duì)差別，在這種情況下評(píng)分者之間的系統(tǒng)誤差也看會(huì)被成是隨機(jī)誤差，此時(shí)計(jì)算出的信度系數(shù)不能正確反映評(píng)分者信度的高低。所以我們選擇內(nèi)部一致性信度分析。3.3.2. 問(wèn)題假設(shè)為簡(jiǎn)化模型，排除品酒員主觀因素和個(gè)人偏好不同而造成的影響，我們做出如下假設(shè)：

16、1、每位品酒員鑒定的同一種葡萄酒特質(zhì)完全相同。 2、對(duì)于每一種類葡萄酒被任何一品嘗相當(dāng)于對(duì)其質(zhì)量進(jìn)行一次客觀的測(cè)試，每一項(xiàng)指標(biāo)的評(píng)價(jià)，相當(dāng)于一道題目的得分。3、葡萄酒的特質(zhì)可以從品酒員的品嘗時(shí)考慮的方面（即客觀測(cè)試的一系列“題目”）測(cè)量出來(lái)。3.3.3. 模型的建立及求解對(duì)于抽象出來(lái)的測(cè)試模型，由于每道測(cè)試題目表示一個(gè)特質(zhì)指標(biāo)的評(píng)價(jià)，所以我們不能采用簡(jiǎn)單的二分法計(jì)分。因而我們選擇cronbach 系數(shù)分析。其中式中，k為整個(gè)品嘗的次數(shù)；為總得分的方差；為是所有在第i位品酒師打分的方差。我們先用excel進(jìn)行數(shù)據(jù)處理，得到的矩陣，以紅葡萄酒為例。 表8. 第一組紅葡萄酒 表9. 第二組紅葡萄酒

17、如表8、表9所示，不難看出對(duì)于紅葡萄酒，第一組標(biāo)準(zhǔn)化cronbachs alpha較高，即第一組品酒師對(duì)紅葡萄酒的評(píng)價(jià)更為可靠。同樣地對(duì)于白葡萄酒，第一組標(biāo)準(zhǔn)化cronbachs alpha比第二組更高，如表10、表11所示。表10. 第一組白葡萄酒 表11. 第二組白葡萄酒 即第一組更為可信。再對(duì)第一組進(jìn)行進(jìn)一步討論，對(duì)于第一組十位品酒師對(duì)紅白葡萄酒的評(píng)價(jià)的項(xiàng)已刪除的 cronbachs alpha 值（詳見(jiàn)附錄）進(jìn)行分析，沒(méi)有任何一位品酒師對(duì)于整體的標(biāo)準(zhǔn)化cronbachs alpha有明顯影響，所以十位品酒師的評(píng)價(jià)分?jǐn)?shù)也均可信。綜上所述，無(wú)論對(duì)于紅葡萄酒還是白葡萄酒，兩組品酒師的評(píng)價(jià)均存

18、在顯著差異，其中第一組品酒師的評(píng)價(jià)結(jié)果更可信，并且第一組十位品酒師的評(píng)價(jià)都屬于可信范圍。3.4. 問(wèn)題2：根據(jù)釀酒葡萄的理化指標(biāo)和葡萄酒的質(zhì)量對(duì)釀酒葡萄進(jìn)行分級(jí)3.4.1. 問(wèn)題分析如果單用葡萄酒為釀酒葡萄打分，則直接在葡萄酒的外觀、香氣、口感質(zhì)量權(quán)重中賦予權(quán)值，所得之和即可作為對(duì)釀酒葡萄的評(píng)分。而本題卻引入了另一類變量，就是葡萄的理化指標(biāo)。這是一項(xiàng)很龐大雜亂的數(shù)據(jù)，所以我們想到利用主成分分析法對(duì)這些數(shù)據(jù)進(jìn)行主要因子分析，找出這些指標(biāo)的公共因子，求的每一樣品在這些因子上的成份得分系數(shù)矩陣，給數(shù)據(jù)降低維度，再以此作為研究釀酒葡萄等級(jí)評(píng)定的依據(jù)。3.4.2. 模型建立（以紅葡萄為例）首先對(duì)數(shù)據(jù)進(jìn)行

19、歸一化處理，再用spss因子分析對(duì)處理后的數(shù)據(jù)進(jìn)行因子分析。并選擇具有 kaiser 標(biāo)準(zhǔn)化的正交旋轉(zhuǎn)法進(jìn)行處理。得到的結(jié)果如表12所示：表12.成份矩陣a成份12345678花色苷.852 總酚.851 單寧.759 dpph自由基.734 葡萄總黃酮.704 果梗比.585 蛋白質(zhì).584 l-.578 黃酮醇.562 -.524 百粒質(zhì)量-.542 出汁率.533 干物質(zhì)含量 .847 總糖 .792 還原糖 .777 可溶性固形物 .758 氨基酸總量 .559 白藜蘆醇 .812 a .724 可滴定酸 -.599 果皮質(zhì)量 -.597 .530 b .500.590 褐變度.618

20、 .673 多酚氧化酶活力 .661 蘋果酸 .639 酒石酸 果穗質(zhì)量 .610 vc含量 -.546 固酸比 檸檬酸 -.571提取方法 :主成份。 a. 已提取了 8 個(gè)成份。 spss經(jīng)過(guò)因子分析和旋轉(zhuǎn)處理后一共得到了8個(gè)主成分，但在8個(gè)主成分之中，葡萄的理化指標(biāo)所載負(fù)荷主要是分配在前四個(gè)成分當(dāng)中。而觀察這四個(gè)成分的高負(fù)荷指標(biāo)，成分一中高負(fù)荷的指標(biāo)主要有單寧、酮類酚類物質(zhì)以及蛋白質(zhì)等，而這些物質(zhì)，特別是單寧，在很大程度上影響到了葡萄酒的口感，因此可將該成分命名為“口感類物質(zhì)”；成分二中高負(fù)荷的指標(biāo)有糖類、氨基酸類、vc含量等營(yíng)養(yǎng)成分類物質(zhì)的指標(biāo)，可命名為“營(yíng)養(yǎng)類物質(zhì)”；成分三主要有a、

21、b色澤，以及果皮含量指標(biāo)，可命名為“色澤類物質(zhì)”；成分四主要有蘋果酸，酒石酸等指標(biāo)，可命名為“酸類物質(zhì)”。這樣通過(guò)軟件做出的成份得分系數(shù)，就可求得對(duì)應(yīng)的因子變量解，如表13所示：表13.成份得分系數(shù)矩陣成份12345678氨基酸總量-.004.129-.173-.069.000.210-.059.052蛋白質(zhì).089-.041-.113.027-.041.053.139.140vc含量.039-.155-.002-.107.086.071-.212-.110花色苷.115-.007.178.004.002-.064.012-.014酒石酸.009-.005-.165-.015-.048-.00

22、9.017.464蘋果酸.024.000.319-.018.120-.140-.039-.043檸檬酸-.076-.065.142.067.022.021-.122.509多酚氧化酶活力-.108.061.229.111-.100.142.106-.040褐變度-.030-.051.285-.024-.023.065.006-.036dpph自由基.186-.013-.060.082.096.096.042-.065總酚.197.047-.041-.002.022-.021.070-.115單寧.171-.041.065.047.062.042-.220.074葡萄總黃酮.233.031-.0

23、77.064.075-.043.046-.095白藜蘆醇.045-.041-.019.014.357.157-.003-.103黃酮醇-.043.052.004.130.063.460-.058.020總糖.022.245-.041.020.046.017.125-.180還原糖-.047.178-.010-.001-.016.009.024.008可溶性固形物.014.234.005.061.018.030.065-.144可滴定酸.033.032.005.051.033.012-.391.058固酸比-.006.084.040.015-.004-.090.424-.084干物質(zhì)含量.005

24、.189.023.016.012.016-.050.056果穗質(zhì)量.029-.003-.004.342.044-.026.074.121百粒質(zhì)量.053.013-.015.296.065-.015-.025-.073果梗比-.049-.064.033-.120-.056.298-.026-.036出汁率.199-.014.002-.021.003-.184-.050-.021果皮質(zhì)量.017.047.057.411.059.172-.098.028l-.118-.006-.034.107.080.192.199-.207b.048.045.065.106.358-.088-.064.082a.

25、007.013.001.052.360.019-.010-.054提取方法 :主成份。 旋轉(zhuǎn)法 :具有 kaiser 標(biāo)準(zhǔn)化的正交旋轉(zhuǎn)法。 構(gòu)成得分。并據(jù)此算得每個(gè)樣品在每類成分上的得分，結(jié)果如表14所示：表14.基于四類主成分負(fù)荷的樣本得分表樣品口味類物質(zhì)營(yíng)養(yǎng)類物質(zhì)色澤類物質(zhì)酸類物質(zhì)1.00 1.31 0.17 3.05 -1.01 2.00 1.41 0.30 0.22 -0.79 3.00 0.68 2.21 -1.11 -0.66 4.00 -0.88 -0.23 -0.32 -0.60 5.00 0.43 -0.02 -0.23 1.79 6.00 0.02 1.36 -0.31 0

26、.08 7.00 -1.01 0.54 0.41 -1.31 8.00 -0.35 -0.13 2.94 0.99 9.00 2.78 -0.73 -0.89 -0.73 10.00 0.11 -2.66 -0.39 -0.82 11.00 -1.06 0.05 -0.13 -0.36 12.00 -0.84 1.74 -0.14 0.48 13.00 -0.27 -0.76 -0.74 -0.57 14.00 -0.38 -0.29 1.22 0.52 15.00 -0.18 -0.38 0.24 0.08 16.00 -0.31 -0.74 0.17 -1.27 17.00 0.53 0.

27、93 -0.25 2.06 18.00 -0.84 0.96 -0.04 -0.19 19.00 -0.06 -0.14 -0.27 -0.69 20.00 -0.40 -0.05 -0.99 1.02 21.00 -0.23 0.49 -1.18 -1.32 22.00 -0.48 0.75 0.26 -0.98 23.00 2.40 -0.09 -0.34 0.93 24.00 0.21 -0.07 -0.08 1.18 25.00 -0.48 -1.68 -0.61 0.42 26.00 -0.86 -1.01 -0.42 1.72 27.00 -1.27 -0.52 -0.06 0.0

28、4由于要求對(duì)葡萄進(jìn)行評(píng)級(jí)，因此我們決定用百分制的得分來(lái)作為評(píng)級(jí)的依據(jù)，由于時(shí)間關(guān)系我們對(duì)評(píng)分的系統(tǒng)只能做一個(gè)簡(jiǎn)單化的處理?，F(xiàn)在我們掌握的數(shù)據(jù)一共有4個(gè)主成分因子得分，還有一個(gè)葡萄酒質(zhì)量得分（由于在一問(wèn)中已經(jīng)論證了第一組的評(píng)論員打分更可信，于是這里葡萄酒質(zhì)量得分以第一組為標(biāo)準(zhǔn)），所以我們將賦予主成分與葡萄酒所得分?jǐn)?shù)1:1:1:1:6的權(quán)重，以此求得釀酒葡萄的最終評(píng)分。如此，則每個(gè)成份所占分?jǐn)?shù)為10分，評(píng)論員評(píng)價(jià)所占分?jǐn)?shù)為60分。最后的等級(jí)評(píng)價(jià)我們將效仿目前權(quán)威的葡萄酒評(píng)分系統(tǒng)-美國(guó)著名的葡萄酒評(píng)論家羅伯特帕克推崇的是葡萄酒100分制評(píng)分體系，帕克的評(píng)分系統(tǒng)會(huì)給每一款酒一個(gè)基礎(chǔ)的分?jǐn)?shù)（50分）。在

29、50分的基礎(chǔ)上，按酒的質(zhì)量特點(diǎn)加分，將酒的品質(zhì)分成四類。而我們也效仿這一模式對(duì)釀酒葡萄進(jìn)行分級(jí)，將葡萄分成如下六個(gè)等級(jí)：96-100 extraordinary 頂級(jí)：各類指標(biāo)全面優(yōu)秀的釀酒葡萄。90-95 outstanding 優(yōu)秀：營(yíng)養(yǎng)成分飽滿，味道純正的釀酒葡萄。80-89 above average優(yōu)良：綜合指標(biāo)尚佳的釀酒葡萄。70-79 average 一般：略有瑕疵，尚無(wú)大礙的釀葡萄。60-69 below average低于一般：不值得推薦50-59 unacceptable 次品根據(jù)成份得分系數(shù)以及評(píng)論員評(píng)價(jià)數(shù)據(jù)處理后各樣品的釀酒葡萄得分明細(xì)以及對(duì)應(yīng)等級(jí)如表15所示：表15.

30、 釀酒紅葡萄得分明細(xì)及等級(jí)樣品口味類物質(zhì)營(yíng)養(yǎng)類物質(zhì)色澤類物質(zhì)酸類物質(zhì)瓶酒員評(píng)價(jià)最總得分等級(jí)1.00 6.35 5.82 10.25 0.94 37.62 60.99 低于一般2.00 6.61 6.09 3.39 1.57 48.18 65.84 低于一般3.00 4.82 10.00 0.18 1.96 48.24 65.19 低于一般4.00 0.95 4.99 2.08 2.14 41.16 51.33 次品5.00 4.18 5.43 2.30 9.21 43.98 65.11 低于一般6.00 3.18 8.27 2.12 4.15 43.32 61.04 低于一般7.00 0.64

31、 6.58 3.87 0.04 42.90 54.02 次品8.00 2.27 5.20 10.00 6.84 43.38 67.68 低于一般9.00 10.00 3.96 0.71 1.75 48.90 65.32 低于一般10.00 3.40 0.00 1.93 1.49 44.52 51.34 次品11.00 0.51 5.56 2.55 2.85 42.06 53.53 次品12.00 1.07 9.05 2.51 5.34 32.34 50.31 次品13.00 2.46 3.92 1.06 2.22 44.76 54.42 次品14.00 2.18 4.87 5.82 5.46

32、43.80 62.13 低于一般15.00 2.68 4.68 3.44 4.15 35.22 50.17 次品16.00 2.36 3.94 3.26 0.17 44.94 54.68 次品17.00 4.45 7.38 2.25 10.00 47.58 71.66 一般18.00 1.06 7.44 2.77 3.36 35.94 50.57 次品19.00 2.98 5.18 2.20 1.87 47.16 59.39 次品20.00 2.15 5.37 0.47 6.94 47.16 62.09 低于一般21.00 2.55 6.48 0.00 0.00 46.26 55.29 次品2

33、2.00 1.95 7.00 3.49 1.01 46.32 59.77 次品23.00 9.05 5.27 2.03 6.66 51.36 74.37 一般24.00 3.66 5.33 2.66 7.42 46.80 65.86 低于一般25.00 1.95 2.01 1.38 5.15 41.52 52.02 次品26.00 1.01 3.39 1.84 9.01 44.28 59.53 次品27.00 0.00 4.41 2.71 4.04 43.80 54.96 次品3.4.3. 針對(duì)白葡萄的模型建立依照對(duì)紅葡萄的等級(jí)評(píng)價(jià)，我們也首先對(duì)白葡萄進(jìn)行了主成分分析（結(jié)果見(jiàn)附件），同樣也將白

34、葡萄指標(biāo)分成了對(duì)應(yīng)的4個(gè)主成分，并求得了樣品在每個(gè)成分上的得分，最后進(jìn)行歸總評(píng)分，得出表16的結(jié)果：表16.釀酒白葡萄得分明細(xì)及等級(jí)白葡萄口味類物質(zhì)營(yíng)養(yǎng)類物質(zhì)色澤類物質(zhì)酸類物質(zhì)瓶酒員評(píng)價(jià)最總得分等級(jí)1.00 4.20 5.53 5.97 2.86 49.20 67.76 低于一般2.00 5.77 6.06 4.38 8.69 44.52 69.43 低于一般3.00 4.47 7.89 6.02 10.97 51.18 80.54 優(yōu)良4.00 5.32 10.00 5.89 7.61 47.64 76.47 一般5.00 7.26 8.43 4.54 6.52 42.60 69.35 低于

35、一般6.00 3.39 8.12 7.91 5.96 41.04 66.43 低于一般7.00 3.72 7.49 9.37 2.61 46.50 69.70 低于一般8.00 0.72 5.29 4.50 9.13 42.84 62.49 低于一般9.00 7.47 6.76 7.21 9.63 43.74 74.81 一般10.00 7.77 7.24 7.23 6.19 44.58 73.01 一般11.00 0.99 8.51 6.39 12.61 43.38 71.89 一般12.00 5.91 7.45 7.12 10.00 37.98 68.47 低于一般13.00 0.00 7

36、.87 8.79 7.91 39.54 64.11 低于一般14.00 5.13 8.63 4.69 5.28 43.20 66.92 低于一般15.00 2.41 4.88 7.55 5.46 43.44 63.74 低于一般16.00 3.61 0.00 5.71 5.21 44.40 58.94 次品17.00 2.43 8.58 6.35 9.27 47.28 73.91 一般18.00 1.67 7.74 7.79 4.54 43.86 65.59 低于一般19.00 5.27 2.23 4.81 10.76 43.32 66.40 低于一般20.00 6.98 7.19 4.77

37、4.86 46.68 70.48 一般21.00 5.49 8.93 3.62 6.93 45.84 70.81 一般22.00 6.66 4.49 8.39 8.79 42.60 70.92 一般23.00 6.81 7.22 5.53 8.04 45.54 73.14 一般24.00 8.51 7.64 10.54 7.88 43.98 78.55 一般25.00 7.85 4.33 4.95 8.07 46.26 71.47 一般26.00 10.00 6.48 7.36 8.09 48.78 80.72 優(yōu)良27.00 5.68 5.21 13.62 9.44 38.88 72.83

38、一般28.00 8.14 8.25 6.87 6.55 48.78 78.60 一般3.4.4. 模型的改進(jìn)構(gòu)想我們對(duì)于釀酒葡萄的評(píng)級(jí)把標(biāo)準(zhǔn)主要是建立的數(shù)據(jù)的歸一化以及主成分提取的基礎(chǔ)之上，由于時(shí)間關(guān)系，對(duì)影響釀酒葡萄等級(jí)的給個(gè)因素只是進(jìn)行了一個(gè)簡(jiǎn)單的權(quán)重化處理，所以我們對(duì)此的改進(jìn)構(gòu)想是運(yùn)用層次分析法將各個(gè)因子構(gòu)建一個(gè)多層次，決定權(quán)重因素更多，更為復(fù)雜的等級(jí)評(píng)價(jià)系統(tǒng)，將 作為各因素的比較矩陣；作為權(quán)重向量；利用求得其最大特征向量，并以此作為排序求得最終的權(quán)重向量，最后再進(jìn)行評(píng)級(jí)處理，這樣將使得評(píng)定結(jié)果更具可靠性。3.5. 問(wèn)題3：分析釀酒葡萄與葡萄酒的理化指標(biāo)之間的聯(lián)系先將葡萄酒以及釀酒葡萄的

39、指標(biāo)整合在一起，總共包括35組指標(biāo)。經(jīng)過(guò)對(duì)數(shù)學(xué)建模中相關(guān)性的研究，得知在spss軟件中共有三種研究相關(guān)性的方法，分別為雙變量相關(guān)性分析、偏變量相關(guān)性分析以及距離相關(guān)性分析。而前兩種方法主要是針對(duì)兩組或三組數(shù)據(jù)進(jìn)行關(guān)聯(lián)度測(cè)試，而本題中是要求對(duì)多項(xiàng)指標(biāo)進(jìn)行對(duì)比，因此我們先采用距離相關(guān)分析法對(duì)這35組指標(biāo)進(jìn)行初步的關(guān)聯(lián)度探測(cè)。通過(guò)spss的處理，得到了各指標(biāo)間的相似性距離，也就是pearson相關(guān)系數(shù)?？梢钥吹狡咸丫婆c釀酒葡萄各項(xiàng)指標(biāo)之間存在較多pearson相關(guān)較大的組合，比如酒的花色苷與葡萄的花色苷，酒的總酚、酒總黃酮與葡萄的dpph自由基，酒的單寧與果梗比由于具有相似性的指標(biāo)太多，所以我們抽取

40、幾組具有代表性的指標(biāo)進(jìn)行更進(jìn)一步的分析。3.5.1. 色酒的色澤與釀酒葡萄的顏色比較，由于花色苷是影響酒和酒的釀酒葡萄的重要因素，所以酒中花色苷與葡萄中花色苷的理化指標(biāo)的聯(lián)系反映了葡萄的顏色對(duì)葡萄酒色澤的影響。于是我們對(duì)以上兩組理化指標(biāo)做了一次雙因素相關(guān)性分析，結(jié)果如表17所示： 表17.相關(guān)性 花色苷w花色苷花色苷wpearson 相關(guān)性1.923*顯著性（雙側(cè)）.000n2727花色苷pearson 相關(guān)性.923*1顯著性（雙側(cè)）.000n2727*. 在 .01 水平（雙側(cè)）上顯著相關(guān)。兩者在0.01的置信水平上顯著相關(guān)，且相關(guān)系數(shù)為0.923，說(shuō)明酒的花色苷含量跟葡萄花色苷有很密切的

41、關(guān)系，或者說(shuō)，酒的色澤從很大程度上取決與葡萄的顏色。下面用matlab對(duì)它們進(jìn)行二次擬合，處理數(shù)據(jù)如下：linear model poly2: f(x) = p1*x2 + p2*x + p3coefficients (with 95% confidence bounds): p1 = 0.000412 (-0.002904, 0.003728) p2 = 2.22 (0.9547, 3.485) p3 = 22.21 (-63.17, 107.6)goodness of fit: sse: 2.041e+005 r-square: 0.8516 adjusted r-square: 0.83

42、93 rmse: 92.22該擬合結(jié)果在95%的置信區(qū)間內(nèi)，準(zhǔn)確性較高。最終對(duì)應(yīng)的公式為：圖6.擬合對(duì)比圖3.5.2. 味單寧是決定酒味道的主要理化指標(biāo)，而在通過(guò)上面的處理方法我們得知酒的單寧指標(biāo)與釀酒葡萄的總酚，總黃酮指標(biāo)具有很高的相似性，pearson相關(guān)系數(shù)分別為0.817、0.684.于是這里我們將研究酒的風(fēng)味與釀酒葡萄的總酚，總黃酮的關(guān)系。用spss中多元線性回歸分析對(duì)酒的單寧，葡萄的單寧、總分以及黃酮醇進(jìn)行處理，得到結(jié)果如表18所示：表18. 相關(guān)性酒的單寧葡萄總酚黃酮醇pearson 相關(guān)性酒的單寧1.000.817.579葡萄總酚.8171.000.346黃酮醇.579.346

43、1.000sig. （單側(cè)）酒的單寧.000.001葡萄總酚.000.039黃酮醇.001.039.n酒的單寧272727葡萄總酚272727黃酮醇272727這里三組數(shù)據(jù)之間處理后的sig值都為0，它們?cè)谥眯艆^(qū)間內(nèi)顯著相關(guān)，因此對(duì)它們做回歸擬合的可靠性很高。下面是線性回歸分析對(duì)這三組數(shù)據(jù)進(jìn)行的回歸分析，所得系數(shù)如表19所示：表19.系數(shù)a模型非標(biāo)準(zhǔn)化系數(shù)標(biāo)準(zhǔn)系數(shù)tsig.共線性統(tǒng)計(jì)量b標(biāo)準(zhǔn) 誤差試用版容差vif1(常量)1.893.6942.727.012葡萄總酚.307.046.7016.683.000.8811.136黃酮醇.024.008.3363.206.004.8811.136a.

44、 因變量: 酒的單寧 對(duì)應(yīng)的酒的單寧（y）與葡萄總酚（x1）、黃酮醇（x2）的線性關(guān)系式如下：擬合結(jié)果與原數(shù)據(jù)對(duì)比圖形如下：圖7. 擬合結(jié)果與原數(shù)據(jù)對(duì)比同樣對(duì)白葡萄中酒的單寧（y）與葡萄總酚（x1）、黃酮醇（x2）進(jìn)行多元化回歸分析，所得結(jié)果如表20：表20.系數(shù)a模型非標(biāo)準(zhǔn)化系數(shù)標(biāo)準(zhǔn)系數(shù)tsig.b標(biāo)準(zhǔn) 誤差試用版1(常量)1.222.3943.100.005白葡萄總酚.063.066.254.953.350白葡萄總黃酮.102.116.235.881.387a. 因變量: 白葡萄酒單寧 回歸后的殘差統(tǒng)計(jì)如表21：表21.殘差統(tǒng)計(jì)量a極小值極大值均值標(biāo)準(zhǔn) 偏差n預(yù)測(cè)值1.55183.0721

45、1.8506.3286028殘差-.938142.23665.00000.6410628標(biāo)準(zhǔn) 預(yù)測(cè)值-.9093.717.0001.00028標(biāo)準(zhǔn) 殘差-1.4083.357.000.96228a. 因變量: 白葡萄酒單寧 3.6. 問(wèn)題4.1：因子分析簡(jiǎn)化模型尋求函數(shù)關(guān)系的擬合3.6.1. 問(wèn)題分析 從葡萄和葡萄酒的理化指標(biāo)來(lái)看，其成分較為復(fù)雜，如果直接對(duì)所有理化指標(biāo)和葡萄酒品質(zhì)的關(guān)系進(jìn)行分析，其對(duì)應(yīng)關(guān)系過(guò)于復(fù)雜，難以得出結(jié)論。為簡(jiǎn)化模型，我們對(duì)其進(jìn)行如下處理：1.我們以品酒師信度較高的第一組的總平均分作為代表每一個(gè)樣品酒品質(zhì)的標(biāo)準(zhǔn)。2.對(duì)于各種理化性質(zhì)，我們只考慮一級(jí)指標(biāo)，暫時(shí)不考慮芳香類

46、物質(zhì)（后文會(huì)提出所求函數(shù)關(guān)系關(guān)于芳香類物質(zhì)的修正）3.為減少變量個(gè)數(shù)，我們通過(guò)因子分析、主成分分析及其嵌套將原來(lái)與得分相關(guān)聯(lián)的理化指標(biāo)向量簡(jiǎn)化為兩個(gè)可以代表原來(lái)數(shù)據(jù)群體的因子4.我們通過(guò)matlab 對(duì)總平均分與因子的關(guān)系進(jìn)行三維的擬合，最終得到近似的關(guān)系=3.6.2. 數(shù)據(jù)處理及模型建立：首先，因?yàn)橛绊懫咸丫瀑|(zhì)量的應(yīng)該是各理化成分的相對(duì)含量，為了結(jié)果的準(zhǔn)確性，所以我們有必要對(duì)所有葡萄酒所得總平均分進(jìn)行歸一化處理。另一方面為了減小變量向量的維度，我們利用spss對(duì)紅白葡萄理化性質(zhì)指標(biāo)進(jìn)行因子分析分析，所得旋轉(zhuǎn)成分矩陣見(jiàn)附錄。對(duì)于紅葡萄我們得到8個(gè)因子即可以看成一個(gè)8維的特征向量，對(duì)于白葡萄我們的到了9個(gè)因子即