《投影寻踪方法及其应用幻灯片.ppt》由会员分享，可在线阅读，更多相关《投影寻踪方法及其应用幻灯片.ppt（29页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、投影寻踪方法及其应用第1页，共29页，编辑于2022年，星期六 报告内容报告内容1 投影寻踪方法的基本原理与建模步骤投影寻踪方法的基本原理与建模步骤 4 总结与讨论总结与讨论3 基于加速遗传算法的基于加速遗传算法的投影寻踪等级评价模型投影寻踪等级评价模型2 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型投影寻踪聚类模型第2页，共29页，编辑于2022年，星期六1 1 投影寻踪方法的基本原理与建模型步骤投影寻踪方法的基本原理与建模型步骤1.1 投影寻踪方法的基本原理投影寻踪方法的基本原理 投影寻踪投影寻踪(p(projectionpursuit，PP)方法属于方法属于直接由直接由样样本数据
2、本数据驱动驱动的探索性数据分析方法的探索性数据分析方法。它把高维数据。它把高维数据 x(i,j)通过某种组合投影到低维子空间上通过某种组合投影到低维子空间上 z(i)，对于投影到的，对于投影到的构形，采用投影指标函数构形，采用投影指标函数Q(z(i)来描述投影暴露原来描述投影暴露原系系统统某某种分类排序结构的可能性大小，寻找出使投影指标函数达到最种分类排序结构的可能性大小，寻找出使投影指标函数达到最优（即能反映高维数据结构或特征）的投影值优（即能反映高维数据结构或特征）的投影值z(i)，然后根据，然后根据该投影值来分析高维数据的分类结构特征（如投影寻踪聚类评该投影值来分析高维数据的分类结构特征
3、（如投影寻踪聚类评价模型），或根据该投影值与研究系统的实际输出值之间的散价模型），或根据该投影值与研究系统的实际输出值之间的散点图构造适当的数学模型以模拟系统输出（如投影寻踪等级评点图构造适当的数学模型以模拟系统输出（如投影寻踪等级评价模型）。价模型）。第3页，共29页，编辑于2022年，星期六1 1 投影寻踪方法的基本原理与建模步骤投影寻踪方法的基本原理与建模步骤1.2 投影寻踪方法的建模步骤投影寻踪方法的建模步骤 步骤步骤1 1：高维样本数据的预处理，：高维样本数据的预处理，确定系统输入。确定系统输入。步骤步骤2 2：构造投影指标函数。：构造投影指标函数。步骤步骤3 3：优化投影指标函数。
4、：优化投影指标函数。步骤步骤4 4：建立系统模型。：建立系统模型。第4页，共29页，编辑于2022年，星期六2 2 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型投影寻踪聚类模型AGA-PPCEAGA-PPCE 2.1 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型的建立投影寻踪聚类模型的建立 步骤步骤1 1：评价指标值的归一化处理。：评价指标值的归一化处理。设样本集为设样本集为x*(i,j)|i=1n,j=1p。其中。其中x*(i,j)为第为第i个样本第个样本第j个指标值。为个指标值。为消除各指标值的量纲和统一各指标值的变化范围，可采消除各指标值的量纲和统一各指标值的变化范围，可
5、采用下式进行极值归一化处理：用下式进行极值归一化处理：式式中中，xmin(j)、xmax(j)分分别别为为样样本本集集中中第第j个个指指标标值值的的最最小小值和最大值。值和最大值。第5页，共29页，编辑于2022年，星期六2 2 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型投影寻踪聚类模型AGA-PPCEAGA-PPCE 2.1 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型的建立投影寻踪聚类模型的建立 步骤步骤2：构造投影指标函数。：构造投影指标函数。PP方法就是把方法就是把p维数据维数据x(i,j)|j=1p综合综合成以成以a=(a(1),a(2),a(p)为投影方向的一维投
6、影值为投影方向的一维投影值z(i)然然后后根根据据z(i)|i=1n的的一一维维散散布布图图进进行行分分类类。式式(2.2)中中a为为单单位位长长度度向向量量。在在综综合合投投影影值值时时，要要求求投投影影值值z(i)的的散散布布特特征征应应为为：局局部部投投影影点点尽尽可可能能密密集集，最最好好凝凝聚聚成成若若干干个个点点团团；而而在在整整体体上上投投影影点点团团之之间间尽尽可可能能散散开开。据据此此投投影影指指标标函函数数可可构造为构造为Q(a)=Sz Dz(2.3)式中，式中，Sz为投影值为投影值z(i)的标准差，的标准差，Dz为投影值为投影值z(i)的局部密度，即的局部密度，即 第6页
7、，共29页，编辑于2022年，星期六2 2 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型投影寻踪聚类模型AGA-PPCEAGA-PPCE 2.1 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型的建立投影寻踪聚类模型的建立 步骤步骤3：优化投影指标函数。当各指标值的样本集给定时，投影指标函数：优化投影指标函数。当各指标值的样本集给定时，投影指标函数Q(a)只随投影方向只随投影方向a的变化而变化。可通过求解投影指标函数最大化问的变化而变化。可通过求解投影指标函数最大化问题来估计最大可能暴露高维数据某类特征结构的最佳投影方向题来估计最大可能暴露高维数据某类特征结构的最佳投影方向：max
8、Q(a)=Sz Dz(2.4)这这是是一一个个以以a(j)|j=1p为为优优化化变变量量的的复复杂杂非非线线性性优优化化问问题题，用用模模拟拟生生物物优优胜胜劣劣汰汰规规则则与与群群体体内内部部染染色色体体信信息息交交换换机机制制的的加加速速遗遗传传算算法法（AGA）来求解上述问题较为简便和有效。）来求解上述问题较为简便和有效。第7页，共29页，编辑于2022年，星期六加速遗传算法（加速遗传算法（AGA）步骤步骤1：模型参数的编码：模型参数的编码。设码长为。设码长为e，第，第j个参数的变化区个参数的变化区间为间为aj，bj，j=1p。把这些区间等分成。把这些区间等分成2e-1个子区间：个子区间
9、：cj=aj+Ijdj (2)式式中中子子区区间间长长度度dj=(bj-aj)/(2e-1)是是常常数数；搜搜索索步步数数Ij为为小小于于2e的的十十进进制制整整数数，是是变变数数；j=1p(下下同同)。把把Ij转转化化成成e位位二二进制数进制数ia(j,k)，k=1e(下同下同)，即式，即式(3)：第8页，共29页，编辑于2022年，星期六 步骤步骤2：初始父代群体的随机生成：初始父代群体的随机生成。生成生成n组各组各p个均匀随机数个均匀随机数(简称随机数简称随机数)u(j,i)，i=1n，经下式转换成十进制整数。，经下式转换成十进制整数。Ij(i)=INTu(j,i)2e(4)式式中中IN
10、T为为取取整整函函数数。由由式式(3)得得对对应应二二进进制制数数ia(j,k,i)，它它们们与与n组组模模型型参参数数cj(i)对对应应，并并把把它它们作为初始父代个体群。们作为初始父代个体群。编码与解码的逻辑过程：编码与解码的逻辑过程：cj(i)Ij(i)ia(j,k,i)第9页，共29页，编辑于2022年，星期六 步骤步骤3：父代个体适应能力评价：父代个体适应能力评价。把第把第i组参数代入式组参数代入式(1)得目标函数值得目标函数值fi，fi越小越小表示模型与观测值拟合得越好，适应能力越强，表示模型与观测值拟合得越好，适应能力越强，设第设第i个个体的适应能力与个个体的适应能力与fi成反比
11、。成反比。构造选择概率序列构造选择概率序列pi，把，把0，1区间分成区间分成n个子区间：个子区间：(0，p1)，(p1，p2)，pn-1，pn，它们与，它们与n个个体一一对应，个个体一一对应，fi越小的个体对应的越小的个体对应的子区间长度越大。子区间长度越大。第10页，共29页，编辑于2022年，星期六步骤步骤4：父代个体的概率选择：父代个体的概率选择。生成生成n个随机数个随机数ui|i=1n，若若uipi-1，pi，则选取第则选取第i个父代个体，其二进制数记为个父代个体，其二进制数记为ia1(j,k,i)。同理可得另外的。同理可得另外的n个父代个体个父代个体ia2(j,k,i)。选择是遗传算
12、法的关键，它体现了。选择是遗传算法的关键，它体现了优胜劣汰的思想。优胜劣汰的思想。第11页，共29页，编辑于2022年，星期六 步骤步骤5：父代个体的杂交：父代个体的杂交。由步由步4得到的父代个体配对成得到的父代个体配对成n对双亲。生对双亲。生成两随机数成两随机数u1和和u2，再转成十进制整数：，再转成十进制整数：IU1=INT(1+u1e)，IU2=INT(1+u2e)，设，设IU1IU2（否则互换其值）。杂交是指第（否则互换其值）。杂交是指第i对双亲对双亲ia1(j,k,i)和和ia2(j,k,i)随机变换一段二进制数，随机变换一段二进制数，从而生成第从而生成第i对子代个体：对子代个体：0
13、110100101(5)1010011100(6)(杂交前杂交前)(杂交后杂交后)第12页，共29页，编辑于2022年，星期六步骤步骤6：子代个体的变异：子代个体的变异。生成随机数生成随机数u1，u2，u3和和u4。当。当u10.5时子时子代个体取式代个体取式(5)，否则取式，否则取式(6)，记其二进制数，记其二进制数为为ia(j,k,i)。把。把u2、u3转换为十进制整数：转换为十进制整数：IU1=INT(1+eu2)，IU2=INT(1+eu3)设子代变异的概率设子代变异的概率(称为变异率称为变异率)为为pm。变异是。变异是当当u4pm时时对子代个体的对子代个体的IU1位和位和IU2位的值
14、进行位的值进行翻转操作翻转操作：0110100111(变异前变异前)(变异后变异后)第13页，共29页，编辑于2022年，星期六 步骤步骤7：进化迭代：进化迭代。第第i个子代个体经式个子代个体经式(3)、式、式(2)转化成第转化成第i组组模型参数。这模型参数。这n个子代个体作为新的父代，个子代个体作为新的父代，算法算法转入步转入步3，进入下一轮进化过程，重新评价、选，进入下一轮进化过程，重新评价、选择、杂交、变异择、杂交、变异，如此反复进化，使个体的适，如此反复进化，使个体的适应能力不断提高，直到最优个体的优化准则值应能力不断提高，直到最优个体的优化准则值小于某一指定值或最优个体的优化准则值不
15、再小于某一指定值或最优个体的优化准则值不再改善，则终止进化迭代，算法结束。改善，则终止进化迭代，算法结束。第14页，共29页，编辑于2022年，星期六 (a)初始分布初始分布(b)第第20次进化迭代后的分布次进化迭代后的分布 (c)第第60次进化迭代后的分布次进化迭代后的分布 第15页，共29页，编辑于2022年，星期六步骤步骤8：加速循环：加速循环。用第一次、第二次进化迭代所产生的优秀个体用第一次、第二次进化迭代所产生的优秀个体的变量变化区间作为变量新的初始变化区间，算法的变量变化区间作为变量新的初始变化区间，算法进入步骤进入步骤1，重新运行，重新运行SGA算法，如此加速循环，算法，如此加速
16、循环，优秀个体的变化区间将逐步调整和收缩，与最优优秀个体的变化区间将逐步调整和收缩，与最优点的距离将越来越近，直到最优个体的优化准则点的距离将越来越近，直到最优个体的优化准则函数值小于某一设定值或算法运行达到预定加速函数值小于某一设定值或算法运行达到预定加速(循环循环)次数，结束整个算法的运行。此时，就把次数，结束整个算法的运行。此时，就把当前群体中最佳个体或优秀个体的平均值指定为当前群体中最佳个体或优秀个体的平均值指定为AGA的结果。的结果。第16页，共29页，编辑于2022年，星期六AGA的主要特点的主要特点利用标准遗传算法演化过程中的利用标准遗传算法演化过程中的优秀个体子群优秀个体子群体
17、体来逐步调整、压缩算法的搜索空间。来逐步调整、压缩算法的搜索空间。控制参数设置：控制参数设置：二进制编码长度二进制编码长度e、杂交概率和变、杂交概率和变异率分别固定设置为异率分别固定设置为10、1.0和和1.0；群体规模；群体规模n与与优秀个体数目优秀个体数目s关系：关系：s/nn/(e2e)，(n,s)的常的常用配置有用配置有(300,10)、(400,20)和和(500,30)；每；每次加速循环中次加速循环中AGA只进行两次只进行两次SGA的进化迭代的进化迭代。p个变量、加速循环个变量、加速循环q次，优秀个体包围最优点的次，优秀个体包围最优点的概率为概率为(1-0.52s)pq第17页，共
18、29页，编辑于2022年，星期六 步骤步骤4：建立聚类模型。把由步骤：建立聚类模型。把由步骤3求得的最求得的最佳投影方向佳投影方向a*代入式代入式(2.2)后，得各样本点的后，得各样本点的投影值投影值z*(i)。投影值。投影值z*(i)与与z*(j)越接近，表越接近，表示样本示样本i与样本与样本j越倾向于分为同一类。按越倾向于分为同一类。按z*(i)值从大到小排序，据此可把各指标的样本集值从大到小排序，据此可把各指标的样本集进行分类。进行分类。第18页，共29页，编辑于2022年，星期六2 2 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型投影寻踪聚类模型AGA-PPCEAGA-PPCE
19、 2.2 基于加速遗传算法的基于加速遗传算法的投影寻踪聚类模型在气候区划中的应用投影寻踪聚类模型在气候区划中的应用 气候区划就是把研究区域划分成若干个分区，在同一分区内具有相似的气候条件，气候区划就是把研究区域划分成若干个分区，在同一分区内具有相似的气候条件，以便调整种植结构，因地制宜地发展农业生产。以便调整种植结构，因地制宜地发展农业生产。表表2.1 2.1 气候样本集及其投影值气候样本集及其投影值 气候因子气候因子j j=1j=2j=3j=4j=5j=6j=7j=8j=9j=10样样本本年平均年平均极端最高极端最高极端最低极端最低10年年年降水年降水年日照年日照年均相年均相对对无霜无霜海拔
20、海拔凌凌冻冻天天投影投影值值气温气温/CC气温气温/CC气温气温/CC积积温温/C C 量量/mm数数/h湿度湿度/%期期/d高度高度/m数数/dz*(i)1 1毕节毕节12.933.610.13672.0904.31236.0822501510.615.21.2952 2大方大方11.831.58.83332.81176.91265.9842561700.033.21.2833 3黔西黔西14.135.48.64047.4964.11263.6812741272.114.61.7224 4金沙金沙15.136.06.24703.31049.71091.681304920.07.52.1915
21、 5织织金金14.233.19.54264.21432.61165.6822801319.011.81.7286 6纳纳雍雍13.733.58.44005.61234.31447.7812681457.114.21.6257 7威宁威宁10.431.114.52572.8943.51960.3801902234.563.90.2718 8赫章赫章13.435.711.63948.9892.81400.8792441534.912.41.282 第19页，共29页，编辑于2022年，星期六图2.1 气候样本投影值z*(i)的散布图由表由表2.1和图和图2.1可知：可知：该样本集按投影值该样本集按
22、投影值z*(i)从大到小排序的样本序号依次为样从大到小排序的样本序号依次为样本点本点4、5、3、6、1、2、8和和7。其中，样本。其中，样本4可单独分为可单独分为A类，样本类，样本5、3、6、1、2和和8可分为可分为B类，样本类，样本7可单独分为可单独分为C类，该分类结果与文献类，该分类结果与文献“刘崇欣刘崇欣.黔西北林木气黔西北林木气候区划的聚类分析候区划的聚类分析.农业系统科学与综合研究农业系统科学与综合研究,1997,13(3)”的模糊聚类结果和根据表的模糊聚类结果和根据表2.1各样本的对比分析结果都是一致的。各样本的对比分析结果都是一致的。B类可进一步细分为由样本类可进一步细分为由样本
23、5、3和和6组成的组成的B-1类和由样本类和由样本1、2和和8组成的组成的B-2类。类。根据文献所示各样本的地理位置可知，以上分类具有明确的地理意义：根据文献所示各样本的地理位置可知，以上分类具有明确的地理意义：A类、类、B类和类和C类分别位于研究区域的东部、中部和西部，类分别位于研究区域的东部、中部和西部，B-1类和类和B-2类分别位于研究区域的中南部和中北类分别位于研究区域的中南部和中北部。部。第20页，共29页，编辑于2022年，星期六3基于加速遗传算法的投影寻踪等级评价模型基于加速遗传算法的投影寻踪等级评价模型AGA-PPGE 3.1 AGA-PPGE的建立的建立 步骤步骤1 1：评价
24、指标值的标准化处理。：评价指标值的标准化处理。设样本集为设样本集为x*(i,j)|i=1n,j=1p。其中。其中x*(i,j)为第为第i个样本个样本第第j个指标值。为消除各指标值的量纲和统一各个指标值。为消除各指标值的量纲和统一各指标值的变化范围，可采用下式进行指标值的变化范围，可采用下式进行标准化标准化处处理：理：式式 中中，Ex(j)、Sx(j)分分 别别 为为 原原 第第j个个 评评 价价 指指 标标x*(j,i)|i=1n的均值和标准差。的均值和标准差。第21页，共29页，编辑于2022年，星期六3基于加速遗传算法的投影寻踪等级评价模型基于加速遗传算法的投影寻踪等级评价模型AGA-PP
25、GE 3.1 AGA-PPGE的建立的建立 步骤步骤2：构造投影指标函数。：构造投影指标函数。PP方法就是把方法就是把p维数据维数据x(i,j)|j=1p综合成以综合成以a=(a(1),a(2),a(p)为投影方向的一维投影值为投影方向的一维投影值z(i)式式中中，a为为单单位位长长度度向向量量。在在综综合合投投影影值值时时，要要求求投投影影值值z(i)应应尽尽可可能能多多地地提提取取x(j,i)中中的的变变异异信信息息，即即z(i)的的标标准准差差Sz达达到到尽尽可可能能大大，同同时时要要求求z(i)与与已已知知标标准准等等级值级值y(i)的相关系数的绝对值的相关系数的绝对值|Rzy|达到尽
26、可能大。为此，投影指标函数可构造为达到尽可能大。为此，投影指标函数可构造为Q(a)=Sz|Rzy|式中，式中，Sz和和Rzy分别为分别为 第22页，共29页，编辑于2022年，星期六3基于加速遗传算法的投影寻踪等级评价模型基于加速遗传算法的投影寻踪等级评价模型AGA-PPGE 3.1 AGA-PPGE的建立的建立 步骤步骤3：优化投影指标函数。当给定标准等级：优化投影指标函数。当给定标准等级及其评价指标的样本数据及其评价指标的样本数据y(i)|i=1n和和x*(j,i)|j=1p,i=1n时，投影指标函数时，投影指标函数Q(a)只随投影只随投影方向方向a的变化而变化。不同的投影方向反映不同的数
27、据的变化而变化。不同的投影方向反映不同的数据结构特征，最佳投影方向就是最大可能暴露高维数据结构特征，最佳投影方向就是最大可能暴露高维数据某类特征结构的投影方向。可通过求解投影指标函数某类特征结构的投影方向。可通过求解投影指标函数最大化问题来估计最佳投影方向最大化问题来估计最佳投影方向：max Q(a)=Sz|Rzy|这是一个以这是一个以a(j)|j=1p为优化变量的复杂的非线性为优化变量的复杂的非线性优化问题，用优化问题，用AGA可方便地求解上述优化问题可方便地求解上述优化问题。第23页，共29页，编辑于2022年，星期六3基于加速遗传算法的投影寻踪等级评价模型基于加速遗传算法的投影寻踪等级评
28、价模型AGA-PPGE 3.1 AGA-PPGE的建立的建立 步骤步骤4：建立投影寻踪等级评价模型。把由步骤建立投影寻踪等级评价模型。把由步骤3求得的最佳投影方向的估求得的最佳投影方向的估计值计值a*代入投影值公式后，得第代入投影值公式后，得第i个样本投影值的计算值个样本投影值的计算值z*(i)，根据，根据z*(i)y(i)的散点图建立相应的等级评价模型。的散点图建立相应的等级评价模型。z*(i)与与y(i)之间一般呈单调非降关之间一般呈单调非降关系，当系，当z*(i)值超过某门限值时就判定为最高等级值超过某门限值时就判定为最高等级(N级级)，当指标值低于另，当指标值低于另门限值时就判定为最低
29、等级门限值时就判定为最低等级(1级级)，当，当z*(i)值介于这两门限值之间时则为中值介于这两门限值之间时则为中等等级，可用逻辑斯谛曲线来描述：等等级，可用逻辑斯谛曲线来描述：式式中中y*(i)为为第第i个个样样本本等等级级的的计计算算值值；模模型型参参数数c(1)、c(2)分分别别为为积积分分常常数和增长率，可通过求解如下优化问题来确定：数和增长率，可通过求解如下优化问题来确定：第24页，共29页，编辑于2022年，星期六3基于加速遗传算法的投影寻踪等级评价模型基于加速遗传算法的投影寻踪等级评价模型AGA-PPGE 3.2 AGA-PPGE在洪水灾情等级评价中的应用在洪水灾情等级评价中的应用
30、表表3.1河南省洪水灾情等级标准河南省洪水灾情等级标准灾情等级指标灾情等级指标一般灾一般灾较大灾较大灾大灾大灾特大灾特大灾成灾面积成灾面积(hm2)283.3直接经济损失直接经济损失(亿元亿元)85.0 根据表根据表3 3.1，可用如下方法随机产生各灾情指标值及其对应的标准灾情等级样本系列：，可用如下方法随机产生各灾情指标值及其对应的标准灾情等级样本系列：四个四个灾情等级值一般灾、较大灾、大灾、特大灾分别对应标准灾情等级值灾情等级值一般灾、较大灾、大灾、特大灾分别对应标准灾情等级值1、2、3、4。确定一般灾的左端点值和特大灾的右端点值，这里分别取为一般灾的右端点值的确定一般灾的左端点值和特大灾
31、的右端点值，这里分别取为一般灾的右端点值的0.5倍和特大灾的左端点倍和特大灾的左端点值的值的3倍，这样，所有灾级都有一个范围。倍，这样，所有灾级都有一个范围。利用均匀随机数在每个灾级范围内产生利用均匀随机数在每个灾级范围内产生5个值，考虑到直接经济损失一般与成灾面积具有正相关性，同个值，考虑到直接经济损失一般与成灾面积具有正相关性，同一样本点的直接经济损失的随机数应与成灾面积的随机数相同。一样本点的直接经济损失的随机数应与成灾面积的随机数相同。在灾情等级标准表中取边界值各一次，灾情等级值取与该边界值有关的两个灾级值的算术平在灾情等级标准表中取边界值各一次，灾情等级值取与该边界值有关的两个灾级值
32、的算术平均值。这样得到的样本点如表均值。这样得到的样本点如表3.2序号序号123所示。所示。第25页，共29页，编辑于2022年，星期六表表3.2洪水灾情等级的标准值和洪水灾情等级的标准值和PPGE模型的计算值的对比结果模型的计算值的对比结果序号序号灾情指标灾情指标投投影影洪水灾情等级洪水灾情等级序号序号灾情指标灾情指标投投影影洪水灾情等级洪水灾情等级i x*(1,i)x*(2,i)值值z*(i)标准值标准值计算值计算值i x*(1,i)x*(2,i)值值z*(i)标准值标准值计算值计算值138.707.900-1.1781.01.36917157.3038.600-0.4693.02.486
33、238.507.800-1.1801.01.36618283.3085.0000.4223.53.498332.106.500-1.2151.01.31519556.90167.1002.1714.03.967424.204.900-1.2571.01.25620649.50194.9002.7634.03.987536.407.400-1.1911.01.35021602.30180.7002.4614.03.979646.709.500-1.1351.51.43222446.50134.0001.4664.03.897797.6021.700-0.8402.01.89523694.9020
34、8.5003.0534.03.992860.4012.800-1.0562.01.552195072.929.900-1.0422.01.5749112.6025.200-0.7552.02.0331954148.1320.656-0.6792.02.1561056.2011.800-1.0802.01.5151956203.9227.521-0.4213.02.5591180.6017.600-0.9392.01.7361957179.1024.858-0.5323.02.38912136.7031.000-0.6152.52.2581963375.4694.9270.8344.03.726
35、13259.1076.1000.2513.03.3631964301.2447.8360.1123.03.23314200.1054.400-0.1663.02.9151975141.97116.4390.2573.03.36815280.1083.8000.3993.03.4811982279.84121.1270.7724.03.69916236.1067.6000.0883.03.2091984172.0651.619-0.2893.02.750第26页，共29页，编辑于2022年，星期六图图3.1河南省洪水灾情投影值河南省洪水灾情投影值z*(i)与标准灾情等级与标准灾情等级y(i)的散
36、点图的散点图图图3.1表明，可用表明，可用逻辑斯谛曲线逻辑斯谛曲线来描述来描述z*(i)与与y(i)间的关系，其中参数间的关系，其中参数N=4，c(1)、c(2)通过用通过用AGA优化估计，这样得到的河南省洪水灾情等级的优化估计，这样得到的河南省洪水灾情等级的PPGE模型为模型为用用PPGE模型和河南省模型和河南省1950至至1990年年41a中实际发生的中实际发生的9次大的洪灾损失资料进行灾情次大的洪灾损失资料进行灾情等级评估，结果参见表等级评估，结果参见表3.2序号序号1950至至1984这这9次洪水，它们的灾情等级标准值取神经网络次洪水，它们的灾情等级标准值取神经网络模型的评估结果，可见
37、两模型的评定结果相一致，模型的评估结果，可见两模型的评定结果相一致，PPGE模型的灾级分别率更高。如模型的灾级分别率更高。如1950年年的洪灾，成灾面积和直接经济损失的灾情指标值都在一般灾级与较大灾级的边界处附近，的洪灾，成灾面积和直接经济损失的灾情指标值都在一般灾级与较大灾级的边界处附近，因此因此PPGE模型评估该年的灾级为模型评估该年的灾级为1.574是合理的。是合理的。第27页，共29页，编辑于2022年，星期六 4总结与讨论总结与讨论 （1）应用投影寻踪方法过程中富有创新之）应用投影寻踪方法过程中富有创新之处是如何构造反映研究问题要求的投影指标处是如何构造反映研究问题要求的投影指标函数。函数。（2）用遗传算法等适应性强的优化方法优）用遗传算法等适应性强的优化方法优化投影指标函数和系统模型的参数。化投影指标函数和系统模型的参数。（3）通过应用投影寻踪方法解决各种实际）通过应用投影寻踪方法解决各种实际系统问题，不断提高对这些问题的认识。系统问题，不断提高对这些问题的认识。第28页，共29页，编辑于2022年，星期六谢谢大家！谢谢大家！第29页，共29页，编辑于2022年，星期六