[发明专利]基于灰色神经网络的爆管预测方法有效

专利信息
申请号: 201310151844.1 申请日: 2013-04-27
公开(公告)号: CN103258243B 公开(公告)日: 2016-11-30
发明(设计)人: 徐哲;杨洁;车栩龙;孔亚广;薛安克 申请(专利权)人: 杭州电子科技大学
主分类号: G06Q10/04 分类号: G06Q10/04
代理公司: 杭州君度专利代理事务所(特殊普通合伙) 33240 代理人: 杜军
地址: 310018 浙*** 国省代码: 浙江;33
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种基于灰色神经网络的爆管预测方法。本发明首先,对于给定的爆管因素和爆管率数据序列,通过静态灰色建模,对爆管率序列进行预测。预测结果和原爆管率序列进行比较,得到残差。然后,利用神经网络在这些残差和爆管因素间建立神经网络逼近模型。经过反复训练的神经网络就是残差和所选的灰色模型数据间的映射关系。最后预测时,再将灰色模型的预测值用神经网络的补偿值进行补偿。本发明结合灰色建模方法和神经网络模型建立灰色神经网络模型,克服传统爆管模型需要大量的数据的缺点,能较好的解决小样本预测问题,提高预测精度。
搜索关键词: 基于 灰色 神经网络 预测 方法
【主权项】:
基于灰色神经网络的爆管预测方法,其特征在于该方法包括以下步骤:步骤(1)收集爆管因素、整理统计爆管数据并计算爆管率,所述的爆管数据为可量化的爆管因素,包括管径、埋深、管网运行压力和管长;步骤(2)建立GM(0,N)模型;基于步骤(1)收集统计的N‑1个爆管因素和1个爆管率λ,建立以N‑1个爆管因素为因子变量,以爆管率为行为变量的GM(0,N)模型,具体步骤如下:1)设为系统特征数据序列,此序列表示行为变量的n个观测值;<mrow><msubsup><mi>X</mi><mn>2</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>=</mo><mo>{</mo><msubsup><mi>x</mi><mn>2</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msubsup><mi>x</mi><mn>2</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>x</mi><mn>2</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>}</mo></mrow><mrow><msubsup><mi>X</mi><mn>3</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>=</mo><mo>{</mo><msubsup><mi>x</mi><mn>3</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msubsup><mi>x</mi><mn>3</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>x</mi><mn>3</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>}</mo></mrow>...... ...... ...... ...... ......<mrow><msubsup><mi>X</mi><mi>N</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>=</mo><mo>{</mo><msubsup><mi>x</mi><mi>N</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msubsup><mi>x</mi><mi>N</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>x</mi><mi>N</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>}</mo></mrow>这N‑1个序列称为相关因素序列,表示N‑1个因子变量各自的n个观测值;上述各数据序列的一次累加生成1‑AGO序列记为i=1,2,…,N,所谓一次累加生成即:设为原始序列,D为序列算子其中称D为的一次累加生成算子,而生成的新的序列为一次累加生成序列;2)基于序列建立模型<mrow><msubsup><mi>x</mi><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>2</mn></mrow><mi>N</mi></munderover><msub><mi>b</mi><mi>i</mi></msub><msubsup><mi>x</mi><mi>i</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mi>a</mi></mrow>其中参数列可用最小二乘估计而得到为其中:<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>B</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>x</mi><mn>2</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>x</mi><mn>3</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><msubsup><mi>x</mi><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mn>2</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>x</mi><mn>3</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><msubsup><mi>x</mi><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>...</mo></mtd><mtd><mo>...</mo></mtd><mtd><mo>...</mo></mtd><mtd><mo>...</mo></mtd></mtr><mtr><mtd><mrow><msubsup><mi>X</mi><mn>2</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><msubsup><mi>x</mi><mn>3</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><msubsup><mi>x</mi><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow></mtd><mtd><mrow><mi>Y</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msubsup><mi>x</mi><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced>3)应用步骤2)中第2)步模型进行检测,得到预测的行为变量序列为:<mrow><msubsup><mover><mi>X</mi><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>=</mo><mo>{</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>}</mo></mrow>通过以上步骤,就建立了以N‑1个爆管因素为因子变量,以爆管率为行为变量的GM(0,N)模型;步骤(3)建立神经网络模型;将一次累加的相关因素序列作为BP神经网络的输入,由GM(0,N)模型预测得到的特征序列与一次累加的特征数据序列的残差序列E={e(1),e(2),…e(n)}作为网络的输出,其中建立BP神经网络模型;首先,为了避免隐含层某神经元处于饱和状态以及数据间数量级的差异,保证网络对样本具有足够的输入敏感性和良好的拟合性,在对BP神经网络进行训练之前,要对学习祥本数据进行预处理;即对所有数据进行归一化处理,将样本数据转化为[‑1,1]区间上的值;在使用经过学习后的网络时,对网络的输出数据还应进行反归一化,恢复最终预测值;归一化的具体算法是:<mrow><mi>p</mi><mi>n</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mrow><mo>(</mo><mi>p</mi><mo>-</mo><msub><mi>min</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow><mrow><mo>(</mo><msub><mi>max</mi><mi>p</mi></msub><mo>-</mo><msub><mi>min</mi><mi>p</mi></msub><mo>)</mo></mrow></mfrac><mo>-</mo><mn>1</mn></mrow>式中p表示所收集的一组数据;minp表示这组数据中的最小值;maxp表示这组数据中的最大值;pn表示映射后的数据;反归一化具体算法是:<mrow><mi>p</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>&lsqb;</mo><mrow><mo>(</mo><msub><mi>max</mi><mi>p</mi></msub><mo>-</mo><msub><mi>min</mi><mi>p</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>&rsqb;</mo><mo>+</mo><msub><mi>min</mi><mi>p</mi></msub></mrow>然后,应用Matlab中的工具箱用基本的反向传播算法来训练网络,以得到隐含层和输出层相应的权值;这样,经过反复训练的神经网络就是残差序列和一次累加爆管相关因素序列的映射关系;步骤(4)预测爆管率;预测时,先将GM(0,N)模型的预测值用神经网络的补偿值进行误差补偿,以得到预测值然后进行一次累减生成得到<mrow><msubsup><mover><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>=</mo><mo>{</mo><msubsup><mover><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msubsup><mover><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mover><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>^</mo></mover><mn>1</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>}</mo></mrow>所谓一次累减生成即:设为原始序列D为序列算子其中<mrow><msubsup><mi>x</mi><mi>i</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mi>d</mi><mo>=</mo><msubsup><mi>x</mi><mi>i</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>x</mi><mi>i</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo><mi>k</mi><mo>=</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>,</mo><mo>...</mo><mi>n</mi></mrow>称D为的一次累加生成算子,记为1‑IAGO;至此,经过步骤(1),(2),(3),(4),就建立了灰色神经网络的供水管网爆管预测模型。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于杭州电子科技大学,未经杭州电子科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201310151844.1/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top