发明名称 一种机电设备温度预警方法
摘要 本发明公开了一种机电设备温度预警方法,包括步骤:一、机电设备运行温度的检测;二、温度时间序列的传输;三、温度时间序列的分析处理及温度预警,其具体过程如下:步骤301、预测下一采样时刻机电设备温度时间序列的温度预测值,步骤302、确定温度预警识别框架{1,…,k,θ},步骤303、温度预警计算机根据预先建立并训练好的三层BP神经网络,确定出DS证据理论中证据源的基本概率分配;步骤304、温度预警计算机根据DS证据理论的证据组合规则对n个证据源的基本概率分配进行时域融合,得到预警结果M(A)。本发明的方法实现简单,温度预警可靠性高,降低了预警的不确定性,实时性高,适用范围广,便于推广使用。
申请公布号 CN103712702B 申请公布日期 2016.01.06
申请号 CN201410012721.4 申请日期 2014.01.11
申请人 西安科技大学 发明人 黄梦涛;高杏梅;颜一鸣
分类号 G01K1/02(2006.01)I;G06N3/02(2006.01)I 主分类号 G01K1/02(2006.01)I
代理机构 西安创知专利事务所 61213 代理人 谭文琰
主权项 一种机电设备温度预警方法,其特征在于该方法包括以下步骤:步骤一、机电设备运行温度的检测:温度采集装置(1)对机电设备的运行温度进行周期性检测,并将检测得到的温度采样值按照时间的先后顺序进行排列,形成温度时间序列;步骤二、温度时间序列的传输:温度采集装置(1)将所述温度时间序列通过数据通信模块(2)实时传输给温度预警计算机(3);步骤三、温度时间序列的分析处理及温度预警,其具体过程如下:步骤301、温度预警计算机(3)对其接收到的温度时间序列,应用自回归移动平均模型建立温度预测模型ARIMA(p,d,q),并根据温度预测模型ARIMA(p,d,q)得到下一采样时刻的机电设备温度时间序列的温度预测值;其中,温度预测模型ARIMA(p,d,q)的表达式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close='' separators=''><mtable><mtr><mtd><mi>&Phi;</mi><mrow><mo>(</mo><mi>B</mi><mo>)</mo></mrow><msup><mo>&dtri;</mo><mi>d</mi></msup><msub><mi>X</mi><mi>t</mi></msub><mo>=</mo><mi>c</mi><mo>+</mo><mi>&Theta;</mi><mrow><mo>(</mo><mi>B</mi><mo>)</mo></mrow><msub><mi>&epsiv;</mi><mi>t</mi></msub></mtd><mtd><mi>t</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>T</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000455848520000011.GIF" wi="881" he="75" /></maths>其中,X<sub>t</sub>为t时刻的温度预测值,T为样本容量,B为后移算子且BX<sub>t</sub>=X<sub>t‑1</sub>,X<sub>t‑1</sub>为t‑1时刻的温度采样值,B<sup>k</sup>为k步后移算子且B<sup>k</sup>X<sub>t</sub>=X<sub>t‑k</sub>,X<sub>t‑k</sub>为t‑k时刻的温度采样值,<img file="FDA0000455848520000012.GIF" wi="37" he="43" />为后向差分算子且<img file="FDA0000455848520000013.GIF" wi="206" he="52" />d为差分次数,c为常数,ε<sub>t</sub>为t时刻的白噪声;Φ(B)为自回归算子且Φ(B)=1‑φ<sub>1</sub>B‑φ<sub>2</sub>B<sup>2</sup>‑…‑φ<sub>p</sub>B<sup>p</sup>,p为自回归阶数,φ<sub>1</sub>,φ<sub>2</sub>,…,φ<sub>p</sub>为自回归系数;Θ(B)为移动平均算子且Θ(B)=1‑θ<sub>1</sub>B‑θ<sub>2</sub>B<sup>2</sup>‑…‑θ<sub>q</sub>B<sup>q</sup>,q为移动平均阶数,θ<sub>1</sub>,θ<sub>2</sub>,…,θ<sub>q</sub>为移动平均系数;步骤302、温度预警计算机(3)根据机电设备运行的温度范围、温度预警等级和DS证据理论确定温度预警识别框架{1,…,k,θ},预警等级编号j=1~k,其中,k为预警等级数,θ为预警结果不确定;步骤303、温度预警计算机(3)根据预先建立并训练好的三层BP神经网络,确定出DS证据理论中证据源的基本概率分配:首先,将步骤一中得到的温度采样值或步骤301中得到的温度预测值输入预先建立并训练好的三层BP神经网络中,得出三层BP神经网络的输出;然后,以三层BP神经网络的输出作为DS证据理论的证据源并根据基本概率分配计算公式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>M</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>B</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mi>&alpha;</mi></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>B</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>o</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000455848520000021.GIF" wi="393" he="295" /></maths>确定出DS证据理论中证据源的基本概率分配;其中,i为机电设备的运行温度的采样时刻顺序且i=1,2,…,n,n为将要进行时域融合的连续采样时刻的个数;j为预警等级编号,k为预警等级数,α为对证据源的信任程度;B<sub>i</sub>(j)为第i采样时刻的温度采样值或温度预测值对应的三层BP神经网络判断结果为j的输出;M<sub>i</sub>(j)为对第i采样时刻的证据源判断结果为j的基本概率分配;M<sub>i</sub>(θ)为对第i采样时刻的证据源判断结果为θ的基本概率分配;其中,预先建立并训练三层BP神经网络的具体过程为:步骤3031、建立三层BP神经网络:以步骤一中得到的温度采样值或步骤301中得到的温度预测值作为BP神经网络的输入,输入节点数为1,以预警等级数k作为BP神经网络的输出节点数,以预警等级对应的期望输出矩阵作为BP神经网络的输出,隐层节点数为输入节点数与输出节点数之和,建立三层BP神经网络;其中,所述期望输出矩阵为k行1列矩阵,预警等级为j时第j行的矩阵元素为1,其余k‑1行的矩阵元素均为0;步骤3032、训练三层BP神经网络:在机电设备运行温度的可能取值范围内,随机提取m个温度值作为所述三层BP神经网络的输入,并以m个温度值分别对应的m个期望输出矩阵作为所述三层BP神经网络的输出,构建训练样本,对所述三层BP神经网络进行训练并得到训练好的三层BP神经网络;步骤304、温度预警计算机(3)根据DS证据理论的证据组合规则对n个证据源的基本概率分配进行时域融合,得到预警结果M(A);所述证据组合规则为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>-</mo><mi>K</mi></mrow></mfrac><mo>&times;</mo><munder><mi>&Sigma;</mi><mrow><mi>a</mi><mo>&cap;</mo><mi>b</mi><mo>=</mo><mi>A</mi></mrow></munder><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>&times;</mo><msub><mi>M</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000455848520000022.GIF" wi="695" he="135" /></maths>其中,K为M<sub>1</sub>(a)对应的证据源与M<sub>2</sub>(b)对应的证据源的不一致因子,A为温度预警识别框架{1,…,k,θ}的任意一个元素且A=1,…,k或θ,a为M<sub>1</sub>对应的焦元且a=1,…,k或θ,b为M<sub>2</sub>对应的焦元且b=1,…,k或θ;M<sub>1</sub>(a)和M<sub>2</sub>(b)为两个独立的证据源对应的基本概率分配。
地址 710054 陕西省西安市雁塔路中段58号