发明名称 限制条件下内河航道船舶大型化预测方法
摘要 本发明公开了一种限制条件下内河航道船舶大型化预测方法。本发明包括如下步骤:第一步,确定限制条件及影响因子;第二步,选择预测模型;第三步,样本选择及预处理;第四步,选取限制条件;第五步,选择向量机类型及参数;第六步,建立模型并确定评价指标;第七步,进行预测分析。本发明通过选取支持向量机模型,以往年过闸平均吨位为影响因子,以航道等级为限制条件下船舶大型化进行预测。选用往年平均吨位作为影响因子,预测值在数据上具有粘滞性,出现缓慢增长,快速增长,稳定增长,平稳四个阶段,为船舶制造业和交通基础设施建设部门提供了时间上的参考和依据。
申请公布号 CN102298706B 申请公布日期 2014.05.07
申请号 CN201110231032.9 申请日期 2011.08.12
申请人 河海大学 发明人 张玮;李俊星;夏莉敏;杨氾;阮桯
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 1.一种限制条件下内河航道船舶大型化预测方法,其特征在于,包括如下步骤:第一步,确定限制条件及影响因子;第二步,采用支持向量机SVM在高维特征空间中建立线性学习机;具体步骤包括:设高维空间中线性模型为:f(x,ω)=(ω·φ(x))+b          (1)其中<img file="FDA0000407614150000011.GIF" wi="94" he="67" />是将样本点映射到高维空间的非线性变换,回归机的凸规划问题表示为:<maths num="0001"><![CDATA[<math><mrow><mi>min</mi><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>|</mo><mo>|</mo><mi>&omega;</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>C</mi><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><msub><mi>&xi;</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>&xi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><img file="FDA0000407614150000013.GIF" wi="916" he="75" /><img file="FDA0000407614150000014.GIF" wi="816" he="78" /><maths num="0002"><![CDATA[<math><mrow><msub><mi>&xi;</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><msubsup><mi>&xi;</mi><mi>i</mi><mo>*</mo></msubsup><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>.</mo><mo>.</mo><mo>.</mo><mi>n</mi><mo>.</mo></mrow></math>]]></maths>其对偶最优化问题为:<img file="FDA0000407614150000016.GIF" wi="1662" he="124" /><maths num="0003"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mtd><mtd><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>&le;</mo><mi>C</mi><mo>,</mo><mn>0</mn><mo>&le;</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mo>*</mo></msubsup><mo>&le;</mo><mi>C</mi><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mrow></math>]]></maths>根据KKT条件,求解对偶问题得到:<img file="FDA0000407614150000019.GIF" wi="1258" he="72" /><maths num="0005"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>+</mo><mi>b</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><mi>b</mi><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><mi>NSV</mi></msub></mfrac><mo>{</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>&Element;</mo><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>C</mi><mo>)</mo></mrow></mrow></munder><mo>[</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><mi>&epsiv;</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>-</mo><msubsup><mi>&alpha;</mi><mi>j</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><mo>+</mo><munder><mi>&Sigma;</mi><mrow><msubsup><mi>&alpha;</mi><mi>j</mi><mo>*</mo></msubsup><mo>&Element;</mo><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mi>C</mi><mo>)</mo></mrow></mrow></munder><mo>[</mo><msub><mi>y</mi><mi>i</mi></msub><mo>+</mo><mi>&epsiv;</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><msub><mi>&alpha;</mi><mi>j</mi></msub><mo>-</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mo>*</mo></msubsup><mo>)</mo></mrow><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>}</mo></mtd></mtr></mtable></mfenced></math>]]></maths>其中,N<sub>NSV</sub>为标准支持向量NSV的数量,K(x<sub>i</sub>,x<sub>j</sub>)为核函数,包括:A,线性核函数K(x,x′)=(x·x′);B,多项式核函数K(x,x′)=((x·x′)+c)<sup>m</sup>c≥0,m为任意正整数;C,径向基函数(RBF)核函数<img file="FDA0000407614150000021.GIF" wi="570" he="140" />其中,||ω||<sup>2</sup>旧表示与模型复杂度相关的因素,C为惩罚函数,σ为参数,ξ为不敏感系数;第三步,样本选择及预处理:将历史年船舶平均吨位统计数据分为训练样本和预测样本,按照式(9)进行归一化处理:<maths num="0007"><![CDATA[<math><mrow><mi>y</mi><mo>=</mo><mfrac><mrow><mi>z</mi><mo>-</mo><msub><mi>z</mi><mi>min</mi></msub></mrow><mrow><msub><mi>z</mi><mi>max</mi></msub><mo>-</mo><msub><mi>z</mi><mi>min</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:y是归一化后样本取值;z是样本值;z<sub>min</sub>是样本训练集中最小样本值;z<sub>max</sub>是样本训练集中最大样本值;将广义船舶平均吨位预测视为一个多变量预测问题,样本集(x,y)的输入向量x中的每一个元素作为一个特征,在下一步从影响因子或者限制条件中选取并进行量化;第四步,选取限制条件或影响因子进行量化并体现在训练样本中;若选取的为限制条件,则直接根据限制条件中的极限值进行量化;若选取的为影响因子,则需根据对预测结果影响程度的大小进行取值;第五步,支持向量机类型及参数选择:选用ξ-SVR作为支持向量机类型,采用RBF核函数对于不同类型的样本拟合,其中RBF核函数的参数c,σ,ξ的选择是由训练样本的质量和容量来决定,具体是将样本分为n份,进行交叉验证,比较准确度和均方根误差,进而确定最优的参数组合;第六步,建立模型确定评价指标:根据向量机的基本结构和参数,构造求解向量机模型的凸最优化问题,得到对偶最优化问题函数后,对余下的数据进行预测能力和泛化能力的检验;评价模型的拟合程度和预测效果的验证指标为平均绝对误差MRE:<maths num="0008"><![CDATA[<math><mrow><mi>MRE</mi><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mfrac><mrow><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>|</mo></mrow><msub><mi>x</mi><mi>i</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,n为样本数量,x<sub>i</sub>为实际值,<img file="FDA0000407614150000032.GIF" wi="46" he="92" />为对应的预测值;第七步,预测分析:利用满足第六步所述评价指标误差限的支持向量机SVM模型,对目标值进行定量预测,并根据预测值定性分析结果的合理性,进行综合评价。
地址 211100 江苏省南京市江宁开发区佛城西路8号
您可能感兴趣的专利