数学建模中的主成分分析法

  • 格式:pdf
  • 大小:1.83 MB
  • 文档页数:2

科技视界

Science&TechnologyVision科技视界全国大学生数学建模竞赛袁已经成为许多高校学生课外科技活动的重要项目袁人们也越来越关注数学建模竞赛遥数学规划尧微分方程尧图论等是较为常见的建模方法遥而近年来袁越来越多的数据处理题目出现在数学建模竞赛当中袁数据处理的任务是降低数据的维数袁保留数据的有用信息遥主成分分析法作为一种主要的数据处理方法袁能够提取变量信息袁减少分析的维度袁使问题变得更简单尧直观遥因此袁尽快掌握主成分分析法的基本知识袁显得尤为迫切遥下面介绍主成分分析法的基本知识袁利用主成分分析法的思想方法建立数学模型遥1主成分分析的基本思想和数学模型1.1主成分分析简介主成分这个概念由美国统计学家KarlPearson在1901年提出袁当时只是进行了非随机变量的讨论遥是从多指标分析出发袁运用统计分析原理与方法提取少数几个彼此不想关的综合性指标而保持其原指标所提供的大量信息的一种统计方法遥1933年Hotelling则将此概念推广到了随机变量中遥主成分分析的原理袁是以较少数的综合变量取代原有的多维变量袁使数据结构简化袁把原指标综合成较少几个主成分袁再以这几个主成分的贡献率为权数进行加权平均袁构造出一个综合评价函数遥作为一种多指标分析方法袁在综合评价函数中袁各主成分的权术为其贡献率袁它反映了该主成分包含原数据的信息量占全部信息量的比重袁这样确定权术是客观尧合理的袁它克服了某些评价方法中人为确定权术的缺陷袁这种方法的计算比较规范袁便于在计算机上实现遥1.2主成分分析基本思想在许多实际问题中袁为了全面系统的反应问题袁我们通常用多个变量来刻画某一事物袁但由于这些变量间具有较强的相关关系袁变量间存在大量的重复信息袁直接用它们分析问题时袁往往会引起极大的误差遥因此人们希望用较少的新指标代替原来较多的旧变量袁同时要求这些新指标尽可能的反应原来的信息遥一般来说袁主成分与原始变量之间的关系院渊1冤各主成分都是原始变量的线性组合遥渊2冤主成分的个数远小于原始变量的个数遥渊3冤各主成分之间互不相关遥渊4冤主成分保留了原始变量的绝大部分信息遥1.3主成分分析的模型假设有n个样本袁有p个观测指标渊p约n冤袁得到原始数据矩阵X=渊X1袁X2袁噎袁Xp冤袁其相关系数矩阵为R遥数学上通常的做法是将原来p个指标做线性组合袁作为新的综合指标遥记这些新的综合指标为Z1袁Z2袁噎袁Zk遥最经典的方法就是用方差来表示遥Z1袁Z2袁噎袁Zk这些新指标之间互不相关袁且方差递减遥因此袁计算相关系数矩阵的特征值为姿1逸姿圆逸噎逸姿p袁向量l1袁l圆袁噎袁lp为相应的单位特征向量袁则第i个主成分为Zi=li'X渊i=1袁2袁噎p冤一般是按累计贡献量的大小取前k个袁多数情况下前几个主成分已代表了原来指标的大部分信息遥2主成分分析法的计算步骤主成分分析法做多指标评价的基本步骤如下院渊1冤对原来的p个指标进项标准化袁Zij=xij-xjsj遥i=1袁2袁噎袁n袁j=1袁2袁噎袁p曰渊其中袁n为样本个数袁p为原始指标的个数袁x为原始指标样本值袁s为样本标准差遥冤渊2冤根据标准化后的数据矩阵求出相关系数矩阵R遥渊3冤求出相关系数矩阵R的特征根姿和特征向量l袁以及贡献率姿i/pi=1移姿i渊4冤确定主成分F1袁F2袁噎袁Fk遥渊5冤计算综合评价值F=渊姿1/pi=1移姿i冤F1+渊姿2/pi=1移姿i冤F2+噎+渊姿k/pi=1移姿i冤Fk3利用Spss进行主成分分析的实例在进行多指标评价时袁由于要求评价结果客观尧全面袁就需要从各个方面用多个指标进行测量袁但这样就使得观测指标间存在信息重叠袁同时还会存在量纲尧累加时如何确定权重系数等问题遥为此袁就可以使用主成分分析法进行信息的浓缩袁并解决权重的确定问题遥本文以全国各市城镇单位就业人员工资水平这一问题来说明主成分评价的用法遥这里引用的是2012年山东省各市按行业分城镇单位就业人员平均工资这一数据袁希望对各地市工资水平给出分析与评价遥数据源自叶山东省统计年鉴2013版曳遥在Spss软件中打开文件之后袁操作步骤如下院渊员冤选择野分析冶寅野降维冶寅野因子分析冶选项遥渊2冤依次选中变量并点向右的箭头按扭遥渊3冤在野描述冶对话框中袁选中野相关系数冶选项组组中的野系数冶复选框遥渊4冤在野抽取冶对话框中袁选中野因子固定数量冶输入数字5遥渊5冤在野得分冶对话框中袁选中野显示因子得分系数矩阵冶遥渊6冤单击野确定冶按钮遥得到输出结果表1遥表员公因子方差

作者简介院刘佳渊1986要冤袁女袁山东淄博人袁从事高等数学教学尧数学建模竞赛指导等工作遥数学建模中的主成分分析法ThePrincipalComponentsAnalysisofMathematicalModeling刘佳(淄博职业学院袁山东淄博255314)揖摘要铱本文通过介绍主成分分析法的基本知识袁使学生在短时间内了解在数学建模中怎样利用主成分分析的思想尧方法建立数学模型作一个初步的探讨遥揖关键词铱主成分分析曰数据分析曰累计贡献率揖Abstract铱Byintroducingthebasicknowledgeofthprincipalcomponentanalysismethod,sothatstudentsinashortperiodoftimetotounderstandtheuseofthought,themethodofprincipalcomponentanalysistoestablishthemathematicalmodelforapreliminarydiscussionofhowinmathematicalmodeling.揖Keywords铱Principalcomponentsanalysis曰Dataanalysis曰Thecumulativecontributionrate初始提取农林牧渔1.000.838采矿业1.000.842

制造业1.000.796电力尧燃气及水的生产供应业1.000.834建筑业1.000.830批发和零售业1.000.890223.com.cn. All Rights Reserved.科技视界

Science&TechnologyVision科技视界提取方法院主成分分析表2解的总方差

提取方法院主成分分析表3成分矩阵

提取方法院主成分已提取5个主成分解的总方差这一表格显示了各主成分解释原始变量总方差的情况袁这里选取5个主成分时累计方差贡献率达到了85.713%袁因此选取前5个主成分可以代表各地市工资水平遥表3给出了主成分系数矩阵袁可以说明主成分在各变量上的载荷袁从而得到各主成分的表达式院F1=0.442x1-0.517x2+0.609x3+噎+0.883x17+0.820x18+0.938x19F2=-0.115x1+0.211x2-0.309x3+噎+0.157x17-0.147x18+0.055x19F3=-0.519x1+0.586x2+0.097x3+噎+0.183x17+0.183x18+0.002x19F4=0.599x1+0.413x2-0.482x3+噎+0.095x17-0.032x18-0.178x19F5=-0.036x1+0.129x2+0.426x3+噎+0.099x17-0.174x18-0.060x194各地市职工平均工资水平综合评价与得分将17个地市的数据带入主成分表达式可得各地市的5的主成分得分遥再利用各因子的方差贡献率作为相应因子的权术可得17个地市的职工平均工资水平的综合得分公式F=52.096/85.713F1+11.289/85.713F2+9.056/85.713F3+8.130/85.713F4+5.142/85.713F5各地市平均工资水平得分最高的前3个地区如表4所示遥表源本文选择的5个主成分集中了原始变量的85.713%的信息袁效果较好遥2012年

袁山东省各地市职工平均工资水平可以用这5个主成分来代替袁利用这5个主成分来综合评测各地市平均工资水平袁得到排名前三位的是青岛尧济南尧东营遥揖参考文献铱咱员暂司守奎.数学建模算法与应用[M].北京院国防工业出版社袁2011院595-601咱2暂林明海.对主成分分析法运用中的十个问题的解析[J].统计与决策,2007渊16冤院16-18.咱责任编辑院曹明明暂交通运输尧仓储和邮政业1.000.732住宿和餐饮业1.000.891信息传输尧软件和信息技术服务业1.000.884金融业1.000.954房地产业1.000.887租赁和商务服务业1.000.838科学研究和技术服务业务1.000.951水利尧环境和公共设施管理业1.000.855居民服务尧修理和其他服务业1.000.788教育1.000.946卫生和社会工作1.000.856文化尧体育和娱乐业1.000.758公共管理尧社会保障和社会管理1.000.917成份初始特征值提取平方和载入合计方差的%累积%合计方差的%累积%19.89852.09652.0969.89852.09652.09622.14511.28963.3852.14511.28963.38531.7219.05672.4421.7219.05672.44241.5458.13080.5711.5458.13080.5715.9775.14285.713.9775.14285.7136.6873.61489.3277.6183.25292.5808.4772.51195.0909.2871.51296.60210.2421.27197.87311.2001.05098.92312.105.55399.47613.059.30999.78614.026.13999.92515.014.075100.000161.106E-0165.822E-016100.00017-5.999E-018-3.157E-017100.00018-2.090E-017-1.100E-016100.00019-1.729E-016-9.101E-016100.000成份12345农林牧渔.442-.115-.519.599-.036采矿业-.517.211.586.413.129制造业.609-.039.097-.482.426电力尧燃气及水的生产供应业.492-.306.596.057.374建筑业.804.336-.186.160-.105批发和零售业.693.517-.229.041.297交通运输尧仓储和邮政业.568-.555-.037-.311.051住宿和餐饮业.802.229-.283.319.117信息传输尧软件和信息技术服务业.243.875.212-.074.094金融业房地产业.916.201.053-.021.060租赁和商务服务业.837-.316-.094-.014-.172科学研究和技术服务业务.900-.102.249.162-.206水利尧环境和公共设施管理业.703-.382-.248.038.390居民服务尧修理和其他服务业.453-.209.551.459-.153教育.946-.101.002.081-.186卫生和社会工作.883.157.183.095.099文化尧体育和娱乐业.820-.147.183-.032-.174公共管理尧社会保障和社会管理.938.055.002-.178-.060.631.294.158-.510-.430