发明名称 基于加权主成分分析的码头功能调整综合测评方法
摘要 本发明属于交通运输规划与管理技术领域,为建立码头功能区调整综合测评模型,并对具体测评结果进行分级,为港口规划调整、码头规划设计调整和管理控制提供技术支持。为此,本发明采用的技术方案是,基于加权主成分分析的码头功能调整综合测评方法,包括下列步骤:码头功能调整综合测评指标体系的建立;权重的确定利用灰色关联度法能客观确定权重这一自然特性,确定码头功能调整综合测评各个测评指标权重;测评模型的建立1)采用灰色关联度法确定各测评指标的权重w<sub>j</sub>(j=1,2,...,26);2)测评指标原始数据的无量纲化处理阈值的选取与级别的划分。本发明主要应用于交通运输规划与管理。
申请公布号 CN103886432A 申请公布日期 2014.06.25
申请号 CN201410136376.5 申请日期 2014.04.04
申请人 天津市市政工程设计研究院 发明人 白子建;柯水平;刘凤松;郑利;张志学;徐汉清;陈强
分类号 G06Q10/06(2012.01)I;G06Q50/30(2012.01)I 主分类号 G06Q10/06(2012.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 刘国威
主权项 1.一种基于加权主成分分析的码头功能调整综合测评方法,其特征是,包括下列步骤:码头功能调整综合测评指标体系的建立采用系统分析法利用测评指标体系的建立及其综合测评实际上是更好的实现系统目标即使系统目标得到改进这一自然规律,通过码头功能调整的本质是码头各个功能子系统更加匹配这一总目标的分解,建立码头装卸能力、码头配套设施、财务分析、经济分析四个子目标,从这四个子目标出发,确定各自的测度指标从而建立码头功能区调整的综合测评指标体系;权重的确定利用灰色关联度法能客观确定权重这一自然特性,确定码头功能调整综合测评各个测评指标权重;①原始数据的处理采用平均值法对原始数据进行处理:即将各测评方案同一指标值相加除以测评方案数目求得平均值做为该指标的参考值,将各个指标的平均值依次排列做为参考数列,设n=26为测评指标的数目,则经过数据处理后的参考数列x<sub>0</sub>(t)表示为:{x<sub>0</sub>(t)}={x<sub>01</sub>,x<sub>02</sub>,...,x<sub>026</sub>}其中x<sub>01</sub>表示第一个测评指标的参考值,x<sub>026</sub>表示第26个测评指标的参考值;②计算关联系数设p为测评方案的数目,与参考数列作关联程度比较的p个数列即比较数列为:<maths num="0001"><![CDATA[<math><mrow><mo>{</mo><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>x</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>}</mo><mo>=</mo><mfenced open='[' close='}'><mtable><mtr><mtd><msub><mi>x</mi><mn>11</mn></msub></mtd><mtd><msub><mi>x</mi><mn>12</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>x</mi><mn>126</mn></msub></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>21</mn></msub></mtd><mtd><msub><mi>x</mi><mn>12</mn></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>x</mi><mn>226</mn></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><msub><mi>x</mi><mrow><mi>p</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>p</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>x</mi><mrow><mi>p</mi><mn>26</mn></mrow></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中x<sub>11</sub>表示第一个测评方案的第一个测评指标值,x<sub>126</sub>表示第一个测评方案的第26个测评指标值,x<sub>p26</sub>表示第<sub>p</sub>个测评方案的第26个测评指标值;将第k(k=1,2,...,p)个测评方案各指标的数值与参考数列x<sub>0</sub>(t)对应指标值的差值的绝对值记为:Δ<sub>ok</sub>(t)=|x<sub>0</sub>(t)-x<sub>k</sub>(t)|t=1,2,...,nx对于第k个测评方案,分别记26个Δ<sub>ok</sub>(t)中的最小数和最大数为Δ<sub>ok</sub>(min)和Δ<sub>ok</sub>(max),对p个比较数列,又记p个Δ<sub>ok</sub>(min)中的最小者为Δ(min),p个Δ<sub>ok</sub>(max)中的最大者为Δ(max),这样Δ(min)和Δ(max)分别是所有p个比较数列在各个测评指标中的绝对差值中的最小者和最大者,于是,第k个方案中第t个测评指标的关联系数通过下式计算:<maths num="0002"><![CDATA[<math><mrow><msub><mi>&zeta;</mi><mi>ok</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&Delta;</mi><mrow><mo>(</mo><mi>min</mi><mo>)</mo></mrow><mo>+</mo><mi>&rho;&Delta;</mi><mrow><mo>(</mo><mi>max</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>&Delta;</mi><mi>ok</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>&rho;&Delta;</mi><mrow><mo>(</mo><mi>max</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>式中ρ为分辩系数,用来削弱Δ(max)过大而使关联系数失真的影响,0<ρ<1;③求关联度计算第t个测评指标的关联度r<sub>t</sub>:<maths num="0003"><![CDATA[<math><mrow><msub><mi>r</mi><mi>t</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>p</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&zeta;</mi><mrow><mn>0</mn><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中ζ<sub>ok</sub>(t)为第k个方案中第t个测评指标的关联系数;④通过关联度确定权重对各个测评指标的关联度进行归一化处理,即每一个指标的关联度除以所有关联度的总和,然后可得到各个指标的权重矩阵:w=[w<sub>1</sub>,w<sub>2</sub>,....,w<sub>j</sub>,....,w<sub>26</sub>]其中w<sub>j</sub>为第j个测评指标的权重;测评模型的建立利用主成分分析法-这一客观测评方法能克服指标间信息重叠的这一自然规律,建立基于加权主成分分析法的码头功能调整综合测评模型,具体步骤为;1)采用灰色关联度法确定各测评指标的权重w<sub>j</sub>(j=1,2,...,26);2)测评指标原始数据的无量纲化处理I.对于成本型指标<maths num="0004"><![CDATA[<math><mrow><msub><mi>y</mi><mi>lj</mi></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>max</mi><mi>x</mi></mrow><mi>j</mi></msub><mo>-</mo><msub><mi>x</mi><mi>lj</mi></msub></mrow><mrow><msub><mrow><mi>max</mi><mi>x</mi></mrow><mi>j</mi></msub><mo>-</mo><mo>-</mo><msub><mrow><mi>min</mi><mi>x</mi></mrow><mi>lj</mi></msub></mrow></mfrac></mrow></math>]]></maths>II.对于效益型指标<maths num="0005"><![CDATA[<math><mrow><msub><mi>y</mi><mi>lj</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mi>lj</mi></msub><mo>-</mo><msub><mrow><mi>min</mi><mi>x</mi></mrow><mi>j</mi></msub></mrow><mrow><msub><mrow><mi>max</mi><mi>x</mi></mrow><mi>j</mi></msub><mo>-</mo><msub><mrow><mi>min</mi><mi>x</mi></mrow><mi>j</mi></msub></mrow></mfrac></mrow></math>]]></maths>III.对于区间型指标<img file="FDA0000487455370000023.GIF" wi="1224" he="378" />其中x<sub>lj</sub>表示第l个测评样本对象或方案对象的第j项指标值,maxx<sub>j</sub>、minx<sub>j</sub>分别是x<sub>lj</sub>中的最大值和最小值,[S<sub>1</sub>,S<sub>2</sub>]为指标的最优区间;将各指标通过无量纲化处理后再乘以100,转化成百分制计数,其中x<sub>lj</sub>表示第l个测评对象(样本,方案)的第j项指标值,共L个测评对象,26个测评指标;3)计算矩阵(y<sub>lj</sub>)<sub>L×26</sub>的相关矩阵R=(r<sub>ij</sub>)<sub>26×26</sub>,其中r<sub>lj</sub>为第l行与第j列的相关系数:(y<sub>lj</sub>)<sub>L×26</sub>是原始数据经过无量纲化处理后的L行26列的矩阵;4)计算R的特征值λ<sub>1</sub>≥λ<sub>2</sub>≥λ<sub>3</sub>≥......≥λ<sub>26</sub>及对应的标准正交化特征向量u<sub>1</sub>,u<sub>2</sub>,...,u<sub>26</sub>;5)求主成分:<maths num="0006"><![CDATA[<math><mrow><msub><mi>Z</mi><mi>k</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>26</mn></munderover><msub><mi>u</mi><mi>kj</mi></msub><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>26</mn></mrow></math>]]></maths>其中u<sub>kj</sub>表示第k个特征值对应的第j个标准正交化特征向量;其中y<sub>j</sub>是第j个标准化指标;6)计算累计方差贡献率E,使其满足<maths num="0007"><![CDATA[<math><mrow><mi>E</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&lambda;</mi><mi>j</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>26</mn></munderover><msub><mi>&lambda;</mi><mi>j</mi></msub><mo>&GreaterEqual;</mo><mn>85</mn><mo>%</mo></mrow></math>]]></maths>其中E的大小表示前m个主成分反映的信息量和的多少,85%表示提出的特征信息的精度,也可以取其它如90%、95%等,具体依据不同的要求确定;7)标准化指标y<sub>j</sub>的加权处理:<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>y</mi><mi>j</mi><mo>*</mo></msubsup><mo>=</mo><msub><mi>w</mi><mi>j</mi></msub><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mn>26</mn></mrow></math>]]></maths>其中w<sub>j</sub>为通过灰色关联度法求得的第j个测评指标的权重;8)计算每个测评对象(样本,方案)的综合测评指数值Z<maths num="0009"><![CDATA[<math><mrow><mi>Z</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>(</mo><msub><mi>&lambda;</mi><mi>k</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>26</mn></munderover><msub><mi>&lambda;</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>26</mn></munderover><msub><mi>u</mi><mi>kj</mi></msub><msubsup><mi>y</mi><mi>j</mi><mo>*</mo></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>以该值大小对包括样本,方案各个测评对象进行排序;阈值的选取与级别的划分根据各个指标的值以及最后的综合测评值确定级别,即级别的划分既要考虑各个分指标值的大小达到一定标准,又要考虑最后的综合测评值达到指定标准,进而保证整体方案的有效性。
地址 300051 天津市和平区营口道239号