• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 14
    • 下载费用:30 金币  

    重庆时时彩和值走势: 一种基于LWPT?DTW的间歇过程不等长时段同步化的方法.pdf

    关 键 词:
    一种 基于 LWPT DTW 间歇 过程 不等 时段 同步 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201610976111.5

    申请日:

    2016.11.07

    公开号:

    CN106354889A

    公开日:

    2017.01.25

    当前法律状态:

    实审

    有效性:

    审中

    法律详情: 实质审查的生效IPC(主分类):G06F 17/30申请日:20161107|||公开
    IPC分类号: G06F17/30 主分类号: G06F17/30
    申请人: 北京化工大学
    发明人: 王建林; 邱科鹏; 于涛; 赵利强
    地址: 100029 北京市朝阳区北三环东路15号
    优先权:
    专利代理机构: 北京思海天达知识产权代理有限公司 11203 代理人: 沈波
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201610976111.5

    授权公告号:

    |||

    法律状态公告日:

    2017.03.01|||2017.01.25

    法律状态类型:

    实质审查的生效|||公开

    摘要

    本发明公开了一种基于LWPT?DTW的间歇过程不等长时段同步化的方法,属于间歇过程建模及监测技术领域。该方法利用LWPT对使用K均值聚类算法后得到的间歇过程不等长时段进行高低频的多级分解,得到子时段不同频段的系数矩阵;再利用DTW对频段系数矩阵进行同步化,并通过ILWPT进行合成,获得同步化后的等长时段。本发明对间歇过程原始时段全部信息进行充分提取,降低了吉布斯现象对信号合成的影响,实现了间歇过程不等长时段的同步化,提高了同步化精度,运算快速且稳定,为间歇过程不等长子时段建模及监测提供了可靠的数据支撑。

    权利要求书

    1.一种基于LWPT-DTW的间歇过程不等长时段同步化的方法,其特征在于:该方法包括
    以下步骤,
    步骤一:采集间歇过程的批次过程数据,通过K均值聚类算法进行子时段的划分;
    步骤二:利用LWPT对各批次过程数据对应的子时段进行高低频的多级分解,得到子时
    段不同频段的系数矩阵;
    步骤三:利用DTW对不同子时段对应的频段系数矩阵进行同步化;
    步骤四:通过提升小波包逆变换对经过DTW同步化后的频段系数矩阵进行合成,得到间
    歇过程等长的子时段;
    步骤五:通过计算同步化精确率选择使同步化精确率达到最高的小波基函数和分解级
    数,获得同步化精度最高的间歇过程不等长子时段同步化结果。
    2.根据权利要求1所述的一种基于LWPT-DTW的间歇过程不等长时段同步化的方法,其
    特征在于:所述步骤一,具体包括:
    采集I个间歇过程批次的过程数据Xh(J×K),其中h(1≤h≤I)为批次的序号,J为测量变
    量的个数,K为样本个数;
    对每个间歇过程批次的过程数据通过K均值聚类算法进行子时段划分;从K个样本中任
    意选取H个样本为初始的聚类中心,分别计算每个样本到每个聚类中心的距离,根据距离大
    小把各样本划归最近的类;重新对每一类计算聚类中心,若聚类中心有变化,则重新计算每
    个样本到每个聚类中心的距离;若聚类中心不发生变化,则将重新计算的聚类中心作为聚
    类结果。
    3.根据权利要求1所述的一种基于LWPT-DTW的间歇过程不等长时段同步化的方法,其
    特征在于:所述步骤二,具体包括:
    选取小波基函数和分解级数,利用LWPT对各批次对应的间歇过程不等长子时段进行多
    级分解;对信号x[n]进行一级LWPT分解的过程如下:
    利用式(1)将x[n]根据奇偶性分解成采样的偶次项xe[n]和奇次项xo[n]
    x[n]=xo[n]+xe[n] (1)
    利用原始过程数据的相关性,用偶次项xe[n]预测奇次项xo[n],设预测算子为P,令预测
    误差为
    d1[n]=xo[n]-P(xe[n]) (2)
    通过预测误差d1[n]调整xe[n],设更新算子为U,令调整后的xe[n]为
    x1[n]=xe[n]+U(d1[n]) (3)
    经过分解、预测和更新三个步骤后,信号x[n]完成了一级分解,得到低频分量x1[n]和高
    频分量d1[n];分别对x1[n]和d1[n]重复以上的分解过程运算,获得原始信号的多级分解;设
    T和R代表待同步化子时段和参考子时段,分别为t×N维和r×N维的数据矩阵,t和r分别为
    待同步化子时段和参考子时段的时间长度,N是被测变量的个数,T和R经过LWPT多级分解后
    得到的频段系数矩阵分别为TCt和RCt,其中t=1,2,...,n,n=2L,L为分解的级数。
    4.根据权利要求1所述的一种基于LWPT-DTW的间歇过程不等长时段同步化的方法,其
    特征在于:所述步骤三,具体包括:
    利用DTW对不同子时段对应的频段系数矩阵进行同步化;设TC和RC分别代表T和R对应
    的某频段系数矩阵,TC和RC分别为a×N维和b×N维的数据矩阵,N是测量变量的个数,a和b
    分别为TC和RC的频段系数个数;TC和RC的局部距离为d(i,j),
    <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>c</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <mi>R</mi> <mi>C</mi> <mo>(</mo> <mrow> <mi>j</mi> <mo>,</mo> <mi>c</mi> </mrow> <mo>)</mo> <mo>-</mo> <mi>T</mi> <mi>C</mi> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> <mn>1</mn> <mo>&le;</mo> <mi>i</mi> <mo>&le;</mo> <mi>a</mi> <mo>,</mo> <mn>1</mn> <mo>&le;</mo> <mi>j</mi> <mo>&le;</mo> <mi>b</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
    基于局部距离采用递推和迭代的方式获得最短标准总体距离D*,
    <mrow> <msup> <mi>D</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mi>T</mi> <mi>C</mi> <mo>,</mo> <mi>R</mi> <mi>C</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <munder> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> <mi>F</mi> </munder> <mo>&lsqb;</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
    F是一条处于a×b网格中的路径,K是路径的长度,k是路径在网格中的序号,N(ω)是与
    ω(k)有关的标准因数,ω(k)为一非负的权函数;设D为正向累积距离矩阵,用动态规划算
    法由d(i,j)得到D(i,j),应用Sakoe-Chiba约束,到达点(i,j)有三条路径,即从(i-1,j)或
    (i-1,j-1)或(i,j-1)到(i,j),即D通过式(6)进行计算
    <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
    初始条件为D(1,1)=2d(1,1);最终得到的D(t,r),即两轨迹之间的最短总体距离;最
    优路径F*的起点为F*(1)=(t,r),根据Sakoe-Chiba约束,F*(2)从(t-1,r)、(t-1,r-1)和(t,
    r-1)中选??;依次类推,直到终点F*(K)=(1,1),K是最优路径F*的长度;从点(1,1)到(t,r),
    计算正向累积距离矩阵D,从点(t,r)到(1,1)寻找最优路径F*;根据最优路径F*对TC进行点
    对点的匹配,设TC'为TC同步化后的频段系数矩阵,当出现TC中的多个点与RC中的一点进行
    匹配时,采取多点平均法完成匹配,即当TC的第i点,i+1点,…i+n-1点与RC中第j个点的特
    征一致时,令
    <mrow> <msup> <mi>TC</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mi>j</mi> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>T</mi> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> <mo>+</mo> <mi>T</mi> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> <mo>+</mo> <mo>...</mo> <mo>+</mo> <mi>T</mi> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> </mrow> <mi>n</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
    同步化后,TC的数据长度和RC的数据长度一致。
    5.根据权利要求1所述的一种基于LWPT-DTW的间歇过程不等长时段同步化的方法,其
    特征在于:所述步骤四,具体包括:
    通过ILWPT对经过DTW同步化后的频段系数矩阵进行合成,得到间歇过程等长的子时
    段;ILWPT只需将步骤二所述的分解等式的左右变量互换和改变符号的正负即可
    xe[n]=x1[n]-U(d1[n]) (8)
    xo[n]=d[n]+P(xe[n]) (9)
    x[n]=xe[n]+xo[n] (10)
    根据上述的合成步骤将T的频段系数矩阵TCt进行合成,得到同步化后的时段T',时段T'
    和时段R具有相同的长度。
    6.根据权利要求1所述的一种基于LWPT-DTW的间歇过程不等长时段同步化的方法,其
    特征在于:所述步骤五,具体包括:
    设为使用不同的小波基函数ψi和分解级数Li得到的同步化时段,选取PCA相似因子
    作为多元时间序列相似性指标M,通过计算时段T与同步化后的时段T'的相似性M,选取最优
    的小波基函数ψ和分解级数L;利用式(11)计算M
    <mrow> <msub> <mi>M</mi> <mrow> <mi>i</mi> <mn>1</mn> <mo>,</mo> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <munderover> <mi>&Sigma;</mi> <mrow> <msub> <mi>a</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>1</mn> <mi>a</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>2</mn> <mi>a</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <msub> <mi>a</mi> <mn>2</mn> </msub> </mrow> </msub> </mrow> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>a</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>1</mn> <mi>a</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>2</mn> <mi>a</mi> <mn>2</mn> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
    其中,i1、i2表示两条多元时间序列,即T和T';λa1、λa2表示协方差矩阵的特征值;表
    示T的负载矩阵的第a1个负载向量与T'的负载矩阵的第a2个负载向量的夹角,A是主元个
    数,M的取值范围为[0,1];M越大,表明T和T'越相似,即同步化精度越高;分别计算基于
    MSDTW的不等长时段同步化和基于LWPT-DTW的不等长时段同步化的PCA相似因子M1和M2,
    设Ζ为设定的小波基函数和分解级数条件下间歇过程所有批次的同步化精确率,即
    选择使Ζ值达到最高的小波基函数和分解级数,并在此条
    件设定下,利用LWPT-DTW获得同步化精度最高的间歇过程不等长子时段同步化结果。

    说明书

    一种基于LWPT-DTW的间歇过程不等长时段同步化的方法

    技术领域

    本发明涉及一种间歇过程不等长时段同步化方法,属于间歇过程建模及监测技术
    领域,尤其涉及一种基于LWPT-DTW的间歇过程不等长时段同步化的方法。

    背景技术

    间歇过程是现代工业中重要的生产方式,在精细化工、制药及食品等领域有着广
    泛的应用。间歇过程具有复杂的过程特性,多时段是间歇过程的固有特征,若将单批次的多
    时段间歇过程当作一个整体进行过程建模及监测,将导致所建过程模型难以准确描述间歇
    过程,并造成较高的间歇过程监测误报率和漏报率;同时间歇过程的不同批次对应子时段
    的操作时间因生产原料的质量差异、初始条件的设置差异等出现不同步,也难以满足基于
    MPCA的间歇过程监测对过程数据长度一致性的要求,因此解决间歇过程不等长时段同步化
    问题,能够为基于MPCA的间歇过程监测提供可靠的过程数据。

    动态时间规整(Dynamic Time Warping,DTW)可根据两轨迹之间的局部特征对轨
    迹进行动态压缩、扩张,进而使两轨迹的长度一致,然而DTW算法计算量大,效率低,制约着
    其实际应用;DTW与离散小波变换构建的MSDTW(Multi-scale DTW)虽然能提高DTW算法的运
    行速度,但MSDTW算法仅仅是对信号的低频部分进行多级分解,没有挖掘高频部分的有用信
    息,且使用离散小波变换在处理具有突变点的信号时会出现吉布斯现象,增加了信号的合
    成误差,使得MSDTW算法具有一定的局限性;因此,DTW、MSDTW算法均难以应用于间歇过程,
    实现其不等长时段同步化。

    发明内容

    本发明的目的在于提高间歇过程不等长时段同步化精度,利用提升小波包变换
    (Lifting Wavelet Package Transform,LWPT)提出一种基于LWPT-DTW的间歇过程不等长
    时段同步化的方法,包括以下步骤:

    步骤一:采集间歇过程的批次过程数据,通过K均值聚类算法进行子时段的划分;

    步骤二:利用LWPT对各批次过程数据对应的子时段进行高低频的多级分解,得到
    子时段不同频段的系数矩阵;

    步骤三:利用DTW对不同子时段对应的频段系数矩阵进行同步化;

    步骤四:通过提升小波包逆变换(InverseLWPT,ILWPT)对经过DTW同步化后的频段
    系数矩阵进行合成,得到间歇过程等长的子时段。

    步骤五:通过计算同步化精确率选择使同步化精确率达到最高的小波基函数和分
    解级数,获得同步化精度最高的间歇过程不等长子时段同步化结果。

    所述步骤一,具体包括:

    采集I个间歇过程批次的过程数据Xh(J×K),其中h(1≤h≤I)为批次的序号,J为
    测量变量的个数,K为样本个数。

    对每个间歇过程批次的过程数据通过K均值聚类算法进行子时段划分。从K个样本
    中任意选取H个样本为初始的聚类中心,分别计算每个样本到每个聚类中心的距离,根据距
    离大小把各样本划归最近的类。重新对每一类计算聚类中心,若聚类中心有变化,则重新计
    算每个样本到每个聚类中心的距离;若聚类中心不发生变化,则将重新计算的聚类中心作
    为聚类结果。

    所述步骤二,具体包括:

    选取小波基函数和分解级数,利用LWPT对各批次对应的间歇过程不等长子时段进
    行多级分解。对信号x[n]进行一级LWPT分解的过程如下:

    利用式(1)将x[n]根据奇偶性分解成采样的偶次项xe[n]和奇次项xo[n]

    x[n]=xo[n]+xe[n] (1)

    利用原始过程数据的相关性,用偶次项xe[n]预测奇次项xo[n],设预测算子为P,令
    预测误差为

    d1[n]=xo[n]-P(xe[n]) (2)

    通过预测误差d1[n]调整xe[n],设更新算子为U,令调整后的xe[n]为

    x1[n]=xe[n]+U(d1[n]) (3)

    经过分解、预测和更新三个步骤后,信号x[n]完成了一级分解,得到低频分量x1
    [n]和高频分量d1[n]。分别对x1[n]和d1[n]重复以上的分解过程运算,获得原始信号的多级
    分解。设T和R代表待同步化子时段和参考子时段,分别为t×N维和r×N维的数据矩阵,t和r
    分别为待同步化子时段和参考子时段的时间长度,N是被测变量的个数,T和R经过LWPT多级
    分解后得到的频段系数矩阵分别为TCt和RCt,其中t=1,2,...,n,n=2L,L为分解的级数。

    所述步骤三,具体包括:

    利用DTW对不同子时段对应的频段系数矩阵进行同步化。设TC和RC分别代表T和R
    对应的某频段系数矩阵,TC和RC分别为a×N维和b×N维的数据矩阵,N是测量变量的个数,a
    和b分别为TC和RC的频段系数个数。TC和RC的局部距离为d(i,j),

    <mrow> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>c</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <mi>R</mi> <mi>C</mi> <mo>(</mo> <mrow> <mi>j</mi> <mo>,</mo> <mi>c</mi> </mrow> <mo>)</mo> <mo>-</mo> <mi>T</mi> <mi>C</mi> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>c</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> <mn>1</mn> <mo>&le;</mo> <mi>i</mi> <mo>&le;</mo> <mi>a</mi> <mo>,</mo> <mn>1</mn> <mo>&le;</mo> <mi>j</mi> <mo>&le;</mo> <mi>b</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>

    基于局部距离采用递推和迭代的方式获得最短标准总体距离D*,

    <mrow> <msup> <mi>D</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mi>T</mi> <mi>C</mi> <mo>,</mo> <mi>R</mi> <mi>C</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <munder> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> <mi>F</mi> </munder> <mo>&lsqb;</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>K</mi> </munderover> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&CenterDot;</mo> <mi>&omega;</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

    F是一条处于a×b网格中的路径,K是路径的长度,k是路径在网格中的序号,N(ω)
    是与ω(k)有关的标准因数,ω(k)为一非负的权函数。设D为正向累积距离矩阵,用动态规
    划算法由d(i,j)得到D(i,j),应用Sakoe-Chiba约束,到达点(i,j)有三条路径,即从(i-1,
    j)或(i-1,j-1)或(i,j-1)到(i,j),即D通过式(6)进行计算

    <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>min</mi> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mi>d</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>

    初始条件为D(1,1)=2d(1,1)。最终得到的D(t,r),即两轨迹之间的最短总体距
    离。最优路径F*的起点为F*(1)=(t,r),根据Sakoe-Chiba约束,F*(2)从(t-1,r)、(t-1,r-1)
    和(t,r-1)中选取。依次类推,直到终点F*(K)=(1,1),K是最优路径F*的长度。从点(1,1)到
    (t,r),计算正向累积距离矩阵D,从点(t,r)到(1,1)寻找最优路径F*。根据最优路径F*对TC
    进行点对点的匹配,设TC'为TC同步化后的频段系数矩阵,当出现TC中的多个点与RC中的一
    点进行匹配时,采取多点平均法完成匹配,即当TC的第i点,i+1点,…i+n-1点与RC中第j个
    点的特征一致时,令

    <mrow> <msup> <mi>TC</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mi>j</mi> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>T</mi> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> <mo>+</mo> <mi>T</mi> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mn>1</mn> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> <mo>+</mo> <mo>...</mo> <mo>+</mo> <mi>T</mi> <mi>C</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>,</mo> <mo>:</mo> <mo>)</mo> </mrow> </mrow> <mi>n</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>

    同步化后,TC的数据长度和RC的数据长度一致。

    所述步骤四,具体包括:

    通过ILWPT对经过DTW同步化后的频段系数矩阵进行合成,得到间歇过程等长的子
    时段。ILWPT只需将步骤二所述的分解等式的左右变量互换和改变符号的正负即可

    xe[n]=x1[n]-U(d1[n]) (8)

    xo[n]=d[n]+P(xe[n]) (9)

    x[n]=xe[n]+xo[n] (10)

    根据上述的合成步骤将T的频段系数矩阵TCt进行合成,得到同步化后的时段T',
    时段T'和时段R具有相同的长度。

    所述步骤五,具体包括:

    设为使用不同的小波基函数ψi和分解级数Li得到的同步化时段,选取PCA相似
    因子作为多元时间序列相似性指标M,通过计算时段T与同步化后的时段T'的相似性M,选取
    最优的小波基函数ψ和分解级数L。利用式(11)计算M

    <mrow> <msub> <mi>M</mi> <mrow> <mi>i</mi> <mn>1</mn> <mo>,</mo> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <msub> <mi>a</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>1</mn> <mi>a</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>2</mn> <mi>a</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&theta;</mi> <mrow> <msub> <mi>a</mi> <mn>1</mn> </msub> <msub> <mi>a</mi> <mn>2</mn> </msub> </mrow> </msub> </mrow> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>a</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>A</mi> </munderover> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>1</mn> <mi>a</mi> <mn>1</mn> </mrow> </msub> <msub> <mi>&lambda;</mi> <mrow> <mi>i</mi> <mn>2</mn> <mi>a</mi> <mn>2</mn> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>

    其中,i1、i2表示两条多元时间序列,即T和T';λa1、λa2表示协方差矩阵的特征值;
    表示T的负载矩阵的第a1个负载向量与T'的负载矩阵的第a2个负载向量的夹角,A是主元
    个数,M的取值范围为[0,1]。M越大,表明T和T'越相似,即同步化精度越高。分别计算基于
    MSDTW的不等长时段同步化和基于LWPT-DTW的不等长时段同步化的PCA相似因子M1和M2,设
    Ζ为设定的小波基函数和分解级数条件下间歇过程所有批次的同步化精确率,即
    选择使Ζ值达到最高的小波基函数和分解级数,并在此条
    件设定下,利用LWPT-DTW获得同步化精度最高的间歇过程不等长子时段同步化结果。

    本发明对间歇过程经K均值聚类算法获得的不等长时段进行高低频的多级分解,
    得到子时段不同频段的系数矩阵,利用DTW对频段系数矩阵进行同步化,并通过ILWPT进行
    合成,获得同步化后的等长时段。本发明能够对间歇过程原始时段全部信息进行充分提取,
    降低了吉布斯现象对信号合成的影响,实现了间歇过程不等长时段的同步化并提高了同步
    化精度,运算快速且稳定,为间歇过程不等长子时段建模及监测提供可靠的数据支撑。

    附图说明

    图1是本发明所述的一种基于LWPT-DTW的间歇过程不等长时段同步化方法流程
    图;

    图2是LWPT-DTW的具体实现图;

    图3是K均值聚类算法的流程图;

    图4是基于K均值聚类算法的第一子时段划分结果图;

    图5是MSDTW的同步化路径图;

    图6是LWPT-DTW的同步化路径图;

    图7是第一子时段的不同小波基函数的同步化精确率分析图;

    图8是本发明方法得到的同步化结果图;

    具体实施方式

    下面结合实例及附图对本发明作进一步的描述,需要说明的是,实例并不限定本
    发明要求?;さ姆段?。

    实施例

    利用青霉素发酵过程仿真平台(PenSim v2.0)生成40批次实验数据,反应时间设
    为400小时,采样间隔为1小时,各批次反应的初始条件设置在正常范围内,间歇过程变量如
    表1所示。

    表1青霉素发酵过程变量说明


    本发明的流程图及具体实现过程如图1及图2所示,对间歇过程各批次过程数据进
    行标准化处理,通过K均值聚类算法分别对每批次数据进行子时段划分,实现的流程图如图
    3所示。图4为其中5个批次的第一子时段划分结果,存在不同步的情况,因此,对不等长子时
    段进行同步化。

    选择其中一个批次数据为参考批次,将其余39个批次对应的子时段分别与参考批
    次数据对应的子时段进行同步化。对其中一个测试批次的第一个子时段与参考批次的第一
    个子时段分别进行基于MSDTW方法和基于LWPT-DTW方法的同步化,获得的同步化路径分别
    如图5和图6所示。经过2级分解后,基于MSDTW方法的同步化路径只有3条,而基于LWPT-DTW
    方法的同步化路径有4条;各个频段的路径不同反映出时段在各个频段中携带的信息不一
    样,因此基于LWPT-DTW方法的同步化可以更加充分提取原始时段的所有信息。

    在实现不等长时段同步化的情况下,通过选取合适的小波基函数和分解级数更多
    地保留待同步化时段的信息。分解级数设置为2级和3级,小波基函数如表2所示。

    表2小波基函数



    分别计算基于MSDTW的不等长时段同步化和基于LWPT-DTW的不等长时段同步化的
    PCA相似因子M1和M2,若M2>M1,表明在设定的小波基函数和分解级数条件下,基于LWPT-DTW
    的不等长时段同步化比基于MSDTW的不等长时段同步化具有更高的精度。选取表2中的小波
    基函数,第一子时段得到的同步化精确率Z如图6所示。因此,当第一子时段的分解级数为2,
    小波基函数为bior2.4时,基于LWPT-DTW的不等长时段同步化精确率最高,为74.36%。

    第一子时段的通风量变化曲线同步化结果如图7所示。其中R为参考时段的通风量
    变化曲线,T为待同步化时段的通风量变化曲线;T1和T2分别为MSDTW方法和LWPT-DTW方法
    的同步化结果,基于LWPT-DTW方法的PCA相似因子为0.9998,基于MSDTW方法的PCA相似因子
    为0.9948,因此基于LWPT-DTW方法的不等长时段同步化具有更高的同步化精度。

    关于本文
    本文标题:一种基于LWPT?DTW的间歇过程不等长时段同步化的方法.pdf
    链接地址://www.4mum.com.cn/p-6027155.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 百度七星彩下期预测 时时彩龙虎有方法吗 免费快乐飞艇计划软件 365彩票手机app下载 重庆时时全天计划单期 1.992高赔率奖金pk10app 广东时时11选五计划软件手机版 云南时时规则介绍 大乐透近100期走势图 快三人工精准计划软件 稳赚不赔的赚钱软件 多彩网极速快三稳赚方法 山东时时是什么 分分彩9码平刷稳赚方案 财神爷精选六肖已公开 时时彩 刷钱 稳赚