当前位置:文档之家› 数学建模葡萄酒评价.docx

数学建模葡萄酒评价.docx

A题:葡萄酒的评价摘要本文主要进行了葡萄酒感官评价的可信度比较、酿酒葡萄评价分级、酿酒葡萄与葡萄酒的理化指标之间的联系、评价结果统计分析等方面的研究。

通过方差分析、层次分析等方法建立模型,解决了葡萄酒的评价问题。

问题一:利用方差分析法对评酒员评价数据进行分析,并用Excel 画出图表(见正文),直观地观察出两组评价数据范围接近,第二组评价数据波动不大,评价数据更可信。

问题二:要求根据酿酒葡萄的理化指标和葡萄酒的质量,对这些酿酒葡萄进行分级,我们认为影响酿酒葡萄品质的因素较多,酿酒葡萄各理化指标之间的关系又是极其复杂的,对其的评价是一个多指标、多属性的问题。

采用系统工程学的层次分析法(AHP )来确定影响葡萄品质的各因素的权重,应用综合评判法,对酿酒葡萄进行了评价和分级。

各等级下葡萄样品数如下表:等级优良中合格葡萄种类红葡萄54108白葡萄8892问题三:利用逐步回归法得到酿酒葡萄与葡萄酒的理化指标之间的关系,并用BP 神经网络进行比较验证。

问题四:通过聚类分析与神经网络相结合,分析酿酒葡萄与葡萄酒的理化指标和葡萄酒质量间的联系。

通过理化指标得到葡萄酒质量评价分数,并与第二组评酒员评价出的葡萄酒质量评价分数对比分析,可知现阶段还不能用酿酒葡萄与葡萄酒的理化指标来评价酒的质量。

本文的建模过程中,对于每个问题都充分考虑了影响因素,一定程度上体现了模型的可靠性,具有较强的适用性和普遍性。

关键词:方差分析 Excel 逐步回归分析Bp 神经网络聚类分析MatlabDPS 数据处理系统一、问题重述通过聘请一些有资质的评酒员品尝葡萄酒,根据他们反馈意见来确定葡萄酒的质量。

酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系,葡萄酒和酿酒葡萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。

已知某一年份一些葡萄酒的评价结果,及该年份这些葡萄酒的和酿酒葡萄的成分数据。

根据上述条件建立数学模型解决以下问题:1.分析两组评酒员的评价结果有无显着性差异,哪一组结果更可信。

2.根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。

3.分析酿酒葡萄与葡萄酒的理化指标之间的联系。

4.分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响,并论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量。

二、问题分析问题一:观察附表 1 中评酒员的评价结果数据,分析得出它们之间的差异。

根据评酒员对各组葡萄酒的评价结果数据,寻求结果数据更加稳定的一组,作为目标,利用求方差的数学模型,对各个数量指标进行分析比较,得出更有可信度组。

问题二:根据不同理化指标对于酿酒葡萄影响各不相同,用层次分析法构造比较矩阵。

计算得到各个因素所对应的权重,定一个分数指标,根据分数对葡萄进行分级。

问题三:题中葡萄与葡萄酒指标数分别为 60、17,考虑因变量太多,用逐步回归分析法建立求解模型,结合神经网络模型进行对比验证。

将多次测试值取平均数,获得可信数据。

问题四:考虑参数过多,为剔除微小影响因素,通过聚类分析法对影响指标进行归类,寻找主要因素,用神经网络建立模型,获得理化参数对葡萄酒的影响关系。

对理化指标仿真得到新的质量指标分数,与第二组评酒员评价数据比较分析,作为论证依据。

三、模型假设及符号说明3.1 模型假设(1)假设评酒员给出的评价数据不存在个人因素。

(2)假设一级指标只与一级指标相互影响,二级指标只与二级指标相互影响。

(3)假设葡萄分级时忽略二级指标对结果的影响。

3.2 符号说明i:表示第 i 个处理观测值总体平均数。

ij :表示试验误差。

i : 表示处理i 对试验结果产生的影响。

x ij:表示i ij 总和。

ss :表示误差平方和。

ess t表示处理间平方。

ss表示总变异的总平方和。

TW{W1,W2 ,....,W m} :表示权重系数集。

r ij (rij1, rij 2,,r ij4 ) :表示隶属度向量。

Vij (rij1, rij 2,rij 4) :表示评价等级。

P r i表示红葡萄的第i 个一级指标。

Pw j:表示白葡萄的第j 个一级指标。

Q r m:表示红葡萄酒的第m 个一级指标。

Q w n:表示白葡萄酒的第n 个一级指标。

pr a:表示红葡萄的第 a 个二级指标。

p w b:表示白葡萄的第 b 个二级指标。

qr c:表示红葡萄酒的第 c 个二级指标。

qw d:表示白葡萄酒的第 d 个二级指标。

四、模型的建立4.1 问题一 :通过建立方差分析模型对两组评酒员对葡萄酒的评分结果进行差异分析。

4.1.1数学模型反应全部观测值总变异的总平方和是个观测值X ij与总平均数x的离均差平方和,记为:SS Tkx )2用n( x i反映重复 n 次的处理间变异,称为处理间平方和记为SS ti 1k nx i )2(x ij为各处理内离均差平方和之和,反映了各处理内的变异即误差,称为处理i 1 j 1内平方和或误差平方和,记为SS e处理内自由度为观测值的总个数减k 处理内自由度记为df e由于:因此:各项平方和除以各自的自由度便得到总均方、处理间均方和处理内均方,分别记为MS T MS t MS e4.1.3F检验通过 MS t与 MS e2的比较来推断是否为零i 是否相等即12k在20的条件下,MSt服从自由度df1k 1与df2k n 1 的F分布。

MS e若实际计算的 F F即 p0.0520。

df ,df,不能否定 :21若F F F即 0.01p0.05 ,否定 H O:20 。

接受 H A:20 df, df2112若F F0.01 df1 , df2,即 p0.01,否定 H O:220 。

接受 H A:差异小的一组评酒员的评价可信。

4.2 问题二4.2.1 建立层析结构模型建立层次模型之前,应对酿酒葡萄进行分析。

通过分析出影响目标相关因素,将评估酿酒葡萄的等级作为目标层的元素。

对葡萄进行评级时可以从酿酒葡萄的理化指标和葡萄酒的质量二个方面分别考虑,将葡萄的理化指标、葡萄酒质量作为第二层的元素。

从中提取相关的类作为第三层的元素,例如从葡萄酒质量中的外观分析、香气分析、口感分析等,和葡萄理化指标中的各指标的含量。

严格对应葡萄与两个评价因素的映射,将第三层的某些类细化为族。

确定的层次模型示例如图4-2 所示。

酿酒葡萄等级葡萄酒质量葡萄理化指标外 香 口 氨 蛋花 基VC 色观 气 感 酸白含 苷分 分 分 总 质析析析量量4.2.2 构造判断矩阵判断矩阵表示针对上一层次某因素而言,本层次与之有关的各因素之间的相对重要性。

其形式如下:b ij 表示对 A 层而言, B 层中因素 b i 对 b j 的相对重要程度,通常取1、3、5、7、9 及其他们的倒数, 2、 4、 6、 8 表示第 i 个因素相对于第j 个因素的影响介于上述两个相邻等级之间。

判断矩阵 B 具有如下特征: b ii 1 、 b ji1 、 b ijbik,其中 (i , j , k 1,2, , n) 。

b ijbjk判断矩阵中的 b ij 是根据经验经过反复研究验证后确定。

用层次分析法应保持判断矩阵的一致性,矩阵中的 b ij 满足上述三条关系式时,说明判断矩阵具有完全的一致性。

层次单排序是根据判断矩阵计算对于上一层某因素而言本层次与之有联系的因素的重要性次序的权重(的排序)。

在矩阵运算中表现为求最大特征值对应的特征向量。

采用方根法计算准则层的因素相对于目标层的层次单排序。

各因素的权重向量为 ( 1 ,2 ,..., n ) T其中对进行归一化处理,得到(1,2 ,...,n )Tii, 其中 njj 1计算矩阵的最大特征根其中 ( A ) 表示向量 A 的第 i 个元素。

计算一致性指标:C. R.= R.C. I.max n,其中 C. I .n1 R.I.RI 为平均随机一致性指标,其变化情况如表4-2 所示,当CR<0.1,认为判断矩阵的一致性是可以接受的,否则应对判断矩阵作适当调整。

4.2.4 层次总排序为了得到层次结构中每一层次的所有因素相对于总目标的层次总排序,需要将计算出来的层次单排序再次进行适当计算。

假设层次结构模型是由目标层 (A) 、准则层 (B) 和指标层 (C) 所组成,准则层有 m个因素,指标层有 n 个因素。

已知 B 层对 A 层的层次单排序为:( 1 , 2 ,..., m )TC 层对 B 层的准则Bj 的层次单排序为:j( 1 j , 2 j ,...,nj)T则 C 层各指标对 A 层的层次总排序的方法为:C 层各指标对目标层的层次总排序为:( 1 ,2 ,..., n ) T评语集是对各个安全因素可能做出的总的评价结果的集合。

根据实际评估需要进行设计。

把评语集定义为以下几个等级:V V1,V2,V3,V4 ={优,良,中等,合格} 。

设 F 是各种等级因素的集合,把 F 中的葡萄等级按照某一准则分类,一般将相近或相似的等级因素分为一组,设 F 中的等级因素有m组,即FF1 , F2 ,..., F m每个F i又有:F i F1i , F2 i ,..., F ni, n 为组成F i的子因素的个数。

依次对各个子因素继续划分下层。

权重系数集是根据各子因素对上一层父因素的重要程度, 对每个子因素分配的权重系数的集合。

每个Fi映射一个函数值Wi ,Wi组成的集合W W1,W2,...,Wm 即为权重系数集,其中Wi满足归一性和非负性条件,即:同理,可以得到下一层的权重系数:W ij (i 1,2,..., m; j 1,2,..., n) 。

权重系数集受主观因素影响较大,特别是当某一因素出现所有的专家都一致认为其中一个因素是重要的,而其它因素为 0 时,则在评估过程重会夸大该因素而忽略其它因素。

前面步骤通过AHP法可以弱化该影响。

通过层次分析法得到酿酒葡萄的分数并进行分级。

4.3 问题三分析酿酒葡萄与葡萄酒的理化指标之间的联系。

由于给出的数据过多,我们采用逐步回归法进行求解。

根据题目建立模型:对回归模型进行标准化得到得出:标准化的回归模型的矩阵:标准化前后的关系得到标准化后的矩阵:采用逆紧凑变换法来对给出的多种指标进行分析。

由( R( 0) E ) 经高斯消元法变换为( E R 1 ) ,既可求出解。

在逐步回归分析中,每引进一个变量或者剔除一个变量,都要对R 进行一次求解求逆紧奏变换法变换。

对给出的因素进行分析,引入方差贡献最大者。

回归平方和越大,回归方程的效果就越好。

得到葡萄与葡萄酒的理化指标之间个各个关系式。

采用 bp 神经网络对给出的葡萄与葡萄酒的理化指标进行建模。

节点输出模型:隐节点输出模型:Q i f (W ij X i Qj )输出节点输出模型:Y k f (TikQj Qk)其中, f 为非线形作用函数;q 为神经单元阈值。

相关主题