2012年全国大学生数学建模竞赛A题一等奖论文
2012年全国大学生数学建模竞赛A题一等奖论文。高教社杯全国大学生数学建模竞赛编号专用页赛区评侧编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):全国统编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):基于数理分析的葡萄及葡萄酒评价体系摘要葡萄酒的质量评价是硏究葪萄酒的一个重要领域,目前葡萄酒的质量主要由评酒师感官评定。但感官评定存在人为因素,业界一自在尝试用葡萄的理化指标或者葡萄洏的理化指标定量评价葡萄洒的质量。本题要求我们根据葡萄以及葡萄酒的相关数据建模,并研究基」理化指标的葡萄酒评价体系的建立对于问题一,我们首先用配对样品t检验方法研究两组评酒员评价差异的显著性,将红葡萄酒与白葡萄酒进行分类处理,用SPSS软件对两组ⅳ酒员的评分的各个指标以及总评分进行了配对样本t检验。得到的部分结果显示:红葡萄酒外观色调、香气质量的评价存在显著性差异,其他单指标的评价不存在显著差异白葡萄、红葡萄以及整休的评价存在显著性差异接着我们建立了数掂可信度评价模型比较两组数据的可信性,将数据的可信度评价转化成对两组评酒员评分的稳定性评价。首先我们对单个评酒员评分与该组所有评酒员评分的均值的偏差进行了分析,偏差不稳定的点就成为噪声点,表明此次评分不稳定。然后我们用两组评酒员评分的偏差的方差衡量评酒员的稳定性。得到第2组的方差明显小于第1组的从而得出了第2组评价数据的可信度更高的结论。对于问题二,我们根据酿酒葡萄的理化指标和葡萄酒质量对葡萄进行了分级。方面,我们对酿酒葡萄的级理化指标的数据进行标准化,基于主成分分析法对其进行了因子分析,并且得到了27种葡萄理化指标的综合得分及其排序(见正文表5)。另一方面,我们又对附录给出的各单指标百分制评分的权重进行评价,并用信息熵法重新确定了权重,用新的权重计算出27种葡萄酒质量的综合得分并排序(见正文表6)。最后我们对两个排名次序用基于模糊数学评价方法将葡萄的等级划分为1-5级(见正文表8)。对于问一,首先我们将众多的葡萄理化指标用主成分分析法综合成6个主因子,并将葡萄等级也列为主因子之一。对葡萄的6个主因子,以及葡萄酒的10个指标用SPSS软件进行偏相关分析,得到酒黃酮与葡萄的等级正相关性较强等结论。之后对相关性较强的主因子和指标作多元线性回归。得到了葡萄酒10个单指标与主因了之间的多元回归方程,该回归方程定量表示两者之间的联系对于问题四,我们首先将葡萄酒的理化指标标准化处理,对葡萄酒的质量与荀萄的6个主因子和葡萄酒的10个单指标作偏相关分析,并求出多元线性回归方程。该方程就表示了葡萄和葡萄酒理化指标对葡萄酒质量的影响。之后,我们通过通径分析方法中的逐步回归分析得到葡萄与葡萄酒的理化指标只确定了葡萄酒质量信息的47%。从而得出了不能用葡萄和葡萄酒的理化指标评价葡萄酒的质量的结论。接着我们还采用通径分析屮的间接通径系数分析求出各自变量之间通过传递作用对应变量的影响,得到单宁与总酚传递性影响较强等结论最后,我们对模型的改进方向以及优缺点进行了讨论。关键词:配对样本t检验数据可信度评价主成分分析模糊数学评价综合评分信息熵偏相关分析多元线性回归1问题重述确定葡萄酒质量时一般是通过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分,然后求和得到其总分,从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系,葡萄酒和酿酒荀萄检测的理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件中给岀∫某年份一些葡萄酒的评价结果,并分別给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。我们需要建立数学模型并且讨论下列问题:1.分析附件1中两组评洒员的评价结果有无显著性差异,并确定哪一组的评价结果更可信。2.根据酿酒葡萄的理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。3.分析酿酒葡萄与葡萄酒的理化指标之间的联系。4.分析酿酒葡萄和葡萄酒的理化指标对葡萄酒质量的影响,并论证能否用荀萄和葡萄酒的理化指标来评价葡萄酒的质量2模型的假设与符号的约定2.1模型的假设与说明(1)评酒员的打分是按照加分制(不采用扣分制);(2)假设20名评酒员的评价八度在同一区间(数据合理,不需要标准化)(3)每位评酒员的系统误差较小,在本问题屮可以忽略不计(4)假设附件中给出的葡萄和荀萄酒理化指标都准确可靠。2.2符号的约定与说明符号符号的意义原假设显著性概率第1组评酒员对第号品种葡萄酒评分的平均值,第2组评洒员对第号品种葡萄酒评分的平均值第一组评酒员对指标评分的偏差的方差,第二组评酒员对指标评分的偏差的方差,=…,第1组10位评酒员对号酒样品第项指标评分的平均分第组第号评酒员对号酒样品第项指标评分与平均值的偏第1组第号评酒员对其项指标评分与平均值的偏差的平均第2组第个评酒员的总体指标偏差的方差重新确立的第项指标的权重第2组10个评酒员的总体指标偏差的方差评酒员指标的平均评分,=葡萄的第项指标,葡萄的第项因子,=葡萄酒的第项理化指标3问题一的分析与求解3.1问题一的分析题冂要求我们根据两组评酒员对27种红葡萄洒和28种白葡萄泙的10个指标相应的打分情况进行分析,并确定两组评酒员对葡萄酒的评价结果是否有显著性差异,然后判断哪组评酒员的评价结果更可信初步分析可知:由于评酒员对颜色、气味等感官指标的衡量人度不同,因此两组评酒员评价结果是否具有显著性差异应该与评价指标的类型有关,不同的评价指标的显著性差异可能会不同。同时,由于红葡萄酒和白葡萄酒的外观、口味竽指标羔异性较大,处理时需要将白葡萄酒和红葡萄酒的评价结果的显著性差昦分开讨论。基于以上分析,我们可以分别两组品尝同一种类酒样品的评酒员的评价结果进行两两配对,分析配对的数据是否满烂配对样品t检验的前提条件,而且根据常识可知评酒员对同一种酒的同一指标的评价在实际中是符合t检验的条件的。接着我们就可以对数据进行多组配对样品的t检验,从而对两组评洒员评价结果的显著性差异进行检验。由于对同一酒样品的评价数据只有两组,我们只能通过评价结果的稳定性来判定结果的可靠性。而每组结果的可靠性乂最终决定于每个评酒员的稳定性,因此将问题转化为对评酒员稳定性的评价。3.2配对样品的t检验简介统计知识指出:配对样本是指对冋一样本进行两次测试所获得的两组数据,或对两个完全相同的样本在不同条件下进行测试所得的两组数据。在本问中我们可以把配对样品理解为有27组两个完全相同的酒样品在两组不同评酒员的检测下得到的两组数据,两组屮各个指标的数据为每组评酒员对该指标打分的平均值配对样品的t检验可检测配对双方的结果是否具有显著性差异,因此就可以检验出配对的双方(第一组与第二组)对葡萄酒的评价结果是否冇差异性型对样品t检验具有的前提条件为:(1)两样品必须配对(2)两样品来源的总体应该满足正态性分布。配对样品t检验基本原理是:求出每对的差值如果两种处理实际上没有差异,则差值的总体均数应当为0,从该总体中抽出的样本其均数也应当在0附近波动;反之,如果两种处理有差异,差值的总体均数就应当远离0,其样本均数也应当远离0。这样,通过检验该差值总体均数是否为0,就可以得知两种处理有无差异。该检验相应的假设为:=,两种处理没有差別,4≠两和处理存在差别3.3葡萄酒配对样品的t检验问题一中配对样品为27组两个完全相同的酒样品在两组不同评酒员的检测下得到的两组数据,其中两组中各个指标的数据为各组10个评酒员对该指标打分的平均值。该问题中的10个指标分别为:外观澄清度、外观色调、香气纯正度、香气浓度、香气质量、口感纯正度、口感浓度、口感持久性、口感质量、平衡/总休评价。根据t检验的原理,对荀萄酒配对样品进行t检验之前我们要对样品进行正态性检验。首先我们根据附件一并处理表格中的数据,得到配对样品的两组数据,绘制红葡萄酒配对样品表格部分数据如表1表1红葡萄酒配对样品数据表澄清度澄清度平衡/整平衡/整(1组均值)(2组均值)体评价(1组体评价(2组均值)均值)2.3.18.4红29.6红263.63.78.8红273.73.78.8白葡萄酒配对样品表格部分数据如表2:表2白葡萄酒配对样品数据表澄清度澄清度平衡/整平衡/整(1组均值)(2组均值)体评价(1组体评价(2组均值)均值)白17.78.4白22.93.19.1日26白273.778.8从上表中我们能看出,将白葡萄酒和红葡萄酒中的每个指标分别进行样品的配对后,每一个指标的配对结果有27对,每一对的双方分别是1组和2组的评酒员对该指标的评分的平均值。3.3.1样本总体的K-S正态性检验配对样品的t检验要求两对应样品的总体满足正态分布,则总体中的样品应该满足正态性或者近似正态性,样本的正态性检验如卜以红葡萄酒的澄清度的27组数据为例分析:利用SPSS软作绘制两样品的直方图和趋势图如图1所示:图1红葡萄酒澄清度两组数据自方图我们假设两组总体数据都服从态分布,利用SPSS软件进行KS忙态性检验的具体结果见附录2.3。两组数据的近似相伴概率值P分别为0.239和0.329,大于我们一般的显著水平0.05则接受原来假设,即两组红葡萄酒的澄清度数据符合近似正态分布同理可用SPSS软件对其他指标的正态性进行检验,得到结果符合实际猜想,都服从近似正态分布。3.3.2葡萄酒配对样品t检验步骤两种葡萄酒的处理过程类似,这里我们以对红葡萄酒谜价结果的差异的显著性分析为例。step1:我们以第一组对葡萄酒的评价结果总体服从正态分布〃σ,以第二组对葡萄酒的评价结果总体服从正态分布μσ。我们已分别从两总体中获得了抽样样本和,并分别进行两样品相互配对。(具体数据见附录2.1)Step2:;引进一个新的随机变量,对应的样本为将配对样本的t检验转化为单样本t检验Step3:建立零假设4=,构造t统计量;Step4:利用SPSS进行配对样品t检验分析,并对结果做出推断3.4显著性差异结果分析3.3.1红葡萄酒各指标差异显著性分析由SPSS软件对红葡萄酒各指标的配对样品讠枍验后,得到各指标的显著性概率分布表。(结果如表3所示)表3红葡萄酒酒各指标显著性概率P指标外观澄清度外观色调香气纯正度香气浓度‖香气质量P0.6140.0020.1510.1000.010指标口感纯正度口感浓度口感持久性口感质量平衡/整体P0.4370.1580.2510.0550.674由统计学知识,如果显著性概率P显著水平α,则不能拒绝零假设,即认为两总体样本的均值不存在显著差异。则根据表3可得:两组评酒员对红葡萄酒各项指标的评价中除外观色调、香气质量存在显著性差异以外,其他8项指标都无显著性差异。3.3.2白葡萄酒各指标差异显著性分析代入白葡萄酒的评价数据,重复以上步骤,得到白荀萄酒各指标的显著性概率分布表。(结果如表4所示)表4白葡萄酒各指标显著性概率P分布表指标外观澄清度外观色调香气纯正度香气浓度香气质量P0,2990.0890.930.2380.714指标口感纯正度口感浓度口感持久性口感质量平衡/整体0,0000.0050.8630.0000.00l分析表4可得:两组评酒员对白葡萄酒各项指标的评价中只有凵感纯正度」感浓度、凵感质量、平衡/整体评价存在显著性差异,其他6项指标都无显著性差异3.3.3葡萄酒总体差异显著性分析(1)红葡萄酒总体差异显著性分析该问题的附件中已经给出了10项指标的杈重,因此将10项指标利用加权合并成总体评价。对于红葡萄酒两组评价结果构造两组配对t检验。得到显著性概率P=0.030
- 2020-12-04下载
- 积分:1
贝叶斯估计与跟踪实用指南
贝叶斯估计与跟踪实用指南pdfBayesianEstimationand Tracking园冒国后固目Bayesian Estimation and TrackingA practical Guide内容简介非高斯系统估计的文献一般都存在缺少实用方法的缺点,而高斯方法又缺乏严密、有组织的推导《贝叶斯估计与跟踪实用指南》一书设法解决这个领域中的不足,为读者提供全面的高斯噪声与非高斯噪声驱动的线性与非线性动态系统估计方法的介绍。本书特点在于提供了贝叶斯估计与跟踪的统一方法,重点强调了在贝叶斯框架下所有跟踪算法的推导过程,还给出了用于计算加权密度积分的高效的数值方法,包括加权高斯积分的线性与非线性卡尔曼滤波器,以及非高斯情况下的粒子滤波器。作者首先给出每种估计方法基本原理的详细推导过程,然后分步骤对每一种方法进行详细说明,使得跟踪滤波器的编码更简单易懂。本书利用了案例分析来展示所研究内容的应用。此外,该书还为每个算法提供了框图,让读者可以自己开发 MATLAB估计方法软件工具箱。《贝叶斯估计与跟踪实用指南》是一本优秀的估计与跟踪方法课程的研究生教材。本书还适用于估计与跟踪领域的科学研究人员、数学研究人员和工程技术人员,对他们具有重要的实际应用价值。对广大科技工作者来说,也是一本非常有价值的参考书。作者简介安东,J霍格,博土,美国约翰霍普金斯大学应用物理实验室科研人员,他为防空和导弹防御系(AMDD)提出了先进的目标跟踪方法。霍格博士的研究领域包括目标跟踪,信号和阵列处理与处理设计,主动、被动雷达和声纳设计,数字通信和编码理论以及时间频率分析。责任编辑:牛旭东xinwu@ndip.cn责任校对:苏向颖上架建议:计算机视觉封面设计:蒋秀芹httpt//www.ndip.enISBN978-7-118-09321-6WILEYWILEY eEsCopies of this book sold withouta Wiley sticker on the cover are9W78z118m0932161>unauthorized and illegal定价:79.00元装备科技译著出版基金贝叶斯估计与跟踪实用指南Bayesian Estimation and Tracking: A Practical Guide[美] Anton J.Haug著王欣于晓译图F荤版社北京著作权合同登记图字:军-2013-200号图书在版编目(CIP)数据贝叶斯估计与跟踪实用指南/(美)豪格(Hang,A.J)著;王欣,于晓译.一北京:国防工业出版社,2014.5(高新科技译丛)书名原文 Bayesian estimation and tracking:apractical guideISBN978-7-118-09321-6I.①贝…Ⅱ.①豪…②王…③于…Ⅲ.①贝叶斯估计-指南②贝叶斯估计一应用一目标跟踪一指南Ⅳ①0211.67-62②TN953-62中国版本图书馆CIP数据核字(2014)第068301号Authorized translation from the English language edition, entitled Bayesian Estimation andTracking, ISBN 978-0-470-62170-7, by Anton J. Haug, published by John Wiley SonsCopyright C 2012 by John Wiley Sons, Inc.All rights reserved. This translation published under license. No part of this book may be repro-duced in any form without the written permission of the original copyrights holderCopies of this book sold without a Wiley sticker on the cover are unauthorized and illegal本书简体中文版由 John Wiley&sons,lno.授权国防工业出版社独家出版。版权所有,侵权必究。阍前有:原社出版发行(北京市海淀区紫竹院南路23号邮政编码100048)北京嘉恒彩色印刷有限责任公司新华书店经售开本710×10001/16印张19%字数370千字2014年5月第1版第1次印刷印数1—2000册定价79.00元本书如有印装错误,我社负责调换)国防书店:(010)88540777发行邮购:(010)88540776发行传真:(010)88540755行业务:(010)88540717前言本书以基本原理为出发点,阐述了贝叶斯估计滤波器的完整的发展历程。书中硏究了由高斯噪声和非髙斯噪声驱动的线性、非线性动态系统。这里假设动态系统是连续的,因为只有在研究离散滤波器的时候,才会在离散时间对系统进行观测。本书的主要目标是给读者提供过去60年来所提出的大部分贝叶斯估计方法的一致性的全面描述,以表明每种以贝叶斯范式为基础的基本思想是如何与条件密度相关的。学习本书内容需要的预备知识包括线性代数、贝叶斯概率理论以及有限差分和插值的数值方法。书中第2章介绍了所有这些主题以及掌握书中其它内容所需的基础知识书中涵盖的部分内容来自于马里兰大学帕克分校数学系的一门研究生课程。该门课程的主要目的是让学生自己开发目标跟踪算法的 MATLAB工具箱。课后给学生们提出一些特定的跟踪问题,要求编写一个或多个 MATLAB子程序来实现某些跟踪方法。一般来说,学生开发的子程序是独立的,并且可以应用于多种难度不大的跟踪问题(例如,纯方位跟踪)。本书不包含任何习题集,使用本书的读者可以根据需要来调整跟踪问题的作业。此外,本书包含4个相当复杂的实用案例,使用本书的老师可以选择其中之一作为编码作业的框架。第一个案例贯穿本书的第Ⅱ和第Ⅲ部分;其余的3个案例则被分成了单独的章节,构成了本书的第Ⅳ部分。本书有两个鲜明特点:①运用大量的表格和图形,给出每种估计方法基本原理的详细推导过程;②每种方法都有非常详细的分步指导,使得跟踪滤波器的编码简单且易于理解。递推贝叶斯估计可以求解转移函数或变换函数的加权条件密度积分。转移函数反映动态状态向量随时间变化的路径,变换函数将状态向量变换成观测向量。有多种求解积分的数值方法,每种都有不同的估计方法。本书第Ⅱ部分和第Ⅲ部分的每个章节都涵盖求解这些积分的一个或多个数值逼近方法,第Ⅱ部分是高斯加权积分的卡尔曼滤波器方法,第Ⅲ部分是未知密度的加权密度积分的粒子滤波器方法。本书是作者在此领域多年研究的成果,希望能为贝叶斯估计和跟踪研究做出贡献。同时,也希望在研究改进这些方法并将其应用到新领域、解决新问题的基础上,本书能够开拓贝叶斯估计和跟踪研究的新方向。Anton J. haug目录刊言Ⅲ第Ⅰ部分预备知识第1章简介…………………………………………31.1贝叶斯推理…1.2贝叶斯层次估计法1.3本书研究范畴1.3.I目标1.3.2章节概述…1.4用 MATLAB⑧进行建模和仿真参考文献…第2章数学基础2.1矩阵线性代数简介9992.1.1矢量与矩阵的约定与记法…2.1.2和与乘积……………………………………………………102.1.3矩阵的逆…2.1.4分块矩阵的逆112.1.5矩阵的平方根…………………………………2.2矢量点发生器…………122.3多参数多维非线性函数估计………………………………152.3.1标量非线性函数估计152.3.2多维非线性函数估计182.4多变量统计概述……232.4.1一般定义32.4.2高斯密度………………………………………25参考文献………………………………………………………31第3章贝叶斯估计的基本概念……333.1贝叶斯估计3.2点估计式343.3基于概率密度函数的递推贝叶斯滤波器…363.4基于状态均值和协方差的递推贝叶斯估计,,,着.383.4.1状态矢量预测……393.4.2状态矢量更新…403.5一般估计方法讨论………………………………42参考文献43第4章实用案例:初步探讨4.1仿真、估计与评价过程,,,,,,,.,,,.非,,..,,鲁非·4.2利用DFAR浮标场进行匀速运动目标跟踪的场景仿真………454.2.1船舶动态模型………………………454.2.2多浮标观测模型464.2.3场景属性………………………………464.3 DIFAR浮标信号处理4.4DFAR似然函数……………………………………………53参考文献55第Ⅱ部分高斯假设:卡尔曼滤波估计器第5章高斯噪声:高斯加权分布的多维积分595.1第3章中重要结论总结5.2回顾卡尔曼滤波器校正(更新)方程的推导…………615.3贝叶斯点预测积分求解高斯密度5.3.1利用仿射变换来简化流程…645.3.2求解高斯加权积分的一般方法66参考文献…………………………68第6章线性卡尔曼滤波器706.1线性动态模型……………………………………………706.2线性观测模型…716.3线性卡尔曼滤波器…6.4IKF在 DIFAR浮标方位估计中的应用,,·,,,·.,,,..,.1.,·Ⅵ参考文献74第7章线性化卡尔曼滤波器:扩展卡尔曼滤波器…767.1一维情况………………767.1.1-维状态预测767.1.2一维状态估计误差方差预测777.1.3一维观测预测方程7.1.4一维预测方程的变换……787.1.5一维线性化EKF过程7.2多维情况………………………807.2.1状态预测方程817.2.2状态协方差预测方程………………………817.2.3观测预测方程……837.2.4多维预测方程的变换837.2.5线性化多维扩展卡尔曼滤波器过程………857.2.6二阶扩展卡尔曼滤波器…867.3多维协方差预测方程的另一种推导……7.4EKF在DFAR船舶跟踪的应用案例分析887.4.1船舶运动动态模型…………………………887.4.2 DIFAR浮标观测模型…………887.4.3初始化卡尔曼滤波器族中所有滤波器907.4.4选定加速度噪声参数………………………917.4.5EKF跟踪滤波器结果参考文献……93第8章 Sigma点类别:有限差分卡尔曼滤波器……948.1一维有限差分卡尔曼滤波器………………………………948.1.1一维有限差分状态预测…948.1.2一维有限差分状态预测……………………………968.1.3一维有限差分观测预测方程…968.1.4一维有限差分卡尔曼滤波器过程………………………968.1.5简化的一维有限差分预测方程978.2多维有限差分卡尔曼滤波器……………………………988.2.1多维有限差分状态预测…988.2.2多维有限差分状态协方差预测1008.2.3多维有限差分观测预测方程101Ⅶ
- 2020-12-07下载
- 积分:1