发明名称 一种利用时延估计进行转子碰摩声发射源定位的方法
摘要 本发明公开了一种利用时延估计进行转子碰摩声发射源定位的方法,包括步骤:10):建立一维线性定位模型;20):设定衰减系数估计<img file="DDA0000227417671.GIF" wi="32" he="64" />的初始迭代值和时间延迟估计<img file="DDA0000227417672.GIF" wi="56" he="64" />的初始迭代值;30):测算<img file="DDA0000227417673.GIF" wi="96" he="64" />的真实值与估计值的估计误差<img file="DEST_PATH_IMAGE001.GIF" wi="24" he="19" />;40):根据估计误差,测算下一步迭代步长;50):测算出迭代后的时间延迟估计和衰减系数估计;60):测算出新的估计误差<img file="DEST_PATH_IMAGE002.GIF" wi="45" he="19" />;70):比较步骤30)和步骤60)测算的两次估计误差的相对误差,如果小于设定值,则由延迟时间估计<img file="DDA0000227417676.GIF" wi="56" he="64" />测得碰摩源的位置;否则,转向步骤40)继续迭代。该方法考虑到转子系统碰摩声发射信号传播过程中出现信号衰减及强噪声干扰问题,通过加入衰减估计和变步长步骤,精确实现碰摩故障源的定位。
申请公布号 CN102928817B 申请公布日期 2014.04.16
申请号 CN201210397465.6 申请日期 2012.10.18
申请人 东南大学 发明人 邓艾东;童航;秦康;曹浩
分类号 G01S5/20(2006.01)I 主分类号 G01S5/20(2006.01)I
代理机构 南京天翼专利代理有限责任公司 32112 代理人 黄明哲;朱芳雄
主权项 1.一种利用时延估计进行转子碰摩声发射源定位的方法,其特征在于,该方法包括以下步骤:10):建立一维线性定位模型:在转子试验台的碰摩导波板上设置第一传感器(1)和第二传感器(2)两个传感器,转子碰摩声发射源位于两个传感器之间,且转子碰摩声发射源和两个传感器位于同一直线上,转子碰摩后产生声发射信号,该声发射信号被两个传感器接收,声发射信号在传播中伴随着衰减和噪声干扰;20):设定衰减系数估计<img file="FDA0000388867720000011.GIF" wi="44" he="82" />的初始迭代值和时间延迟估计<img file="FDA0000388867720000012.GIF" wi="73" he="82" />的初始迭代值:通过碰摩声发射试验装置采集第一传感器(1)接收到的声发射信号x(t)和第二传感器(2)接收到的声发射信号y(t),然后以两路声发射信号的能量比作为衰减系数估计<img file="FDA0000388867720000013.GIF" wi="44" he="82" />的初始迭代值,最后设定一个时间延迟估计<img file="FDA0000388867720000014.GIF" wi="66" he="82" />的初始迭代值;30):测算y(t)的真实值与估计值的估计误差:假设两路声发射信号中,y(t)的能量比x(t)的能量小,转子碰摩声发射源信号先到达第一传感器(1),然后经一个延迟时间Δτ后到达第二传感器(2),则根据式(1)测算出y(t)的真实值与估计值的估计误差:<maths num="0001"><![CDATA[<math><mrow><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>&alpha;</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mover><mi>&tau;</mi><mo>^</mo></mover><mo>)</mo></mrow></mrow></math>]]></maths>式(1)其中,e(t)表示估计误差,<img file="FDA0000388867720000016.GIF" wi="49" he="81" />表示真实衰减系数α的衰减系数估计,<img file="FDA0000388867720000017.GIF" wi="78" he="81" />表示真实时间延迟Δτ的时间延迟估计,<img file="FDA0000388867720000018.GIF" wi="222" he="141" />表示第二传感器(2)接收的信号,此信号是第一传感器(1)接收到的信号x(t)在延时<img file="FDA0000388867720000019.GIF" wi="70" he="81" />之后的信号;40):根据估计误差e(t),测算下一步迭代步长;50):依据式(2),测算出迭代后的时间延迟估计和衰减系数估计;<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>&alpha;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mover><mi>&alpha;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&mu;</mi><mi>&alpha;</mi></msub><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mfrac><mrow><mi>d</mi><mrow><mo>(</mo><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mrow><mo>(</mo><mover><mi>&alpha;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mi>&Delta;</mi><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>&Delta;</mi><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&mu;</mi><mi>&Delta;&tau;</mi></msub><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mfrac><mrow><mi>d</mi><mrow><mo>(</mo><mi>e</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><mi>d</mi><mrow><mo>(</mo><mi>&Delta;</mi><mover><mi>&tau;</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced></math>]]></maths>式(2)其中,<img file="FDA0000388867720000022.GIF" wi="160" he="71" />表示迭代后的衰减系数估计,<img file="FDA0000388867720000023.GIF" wi="95" he="71" />表示迭代前的衰减系数估计,μ<sub>α</sub>表示<img file="FDA0000388867720000024.GIF" wi="49" he="81" />的迭代步长,<img file="FDA0000388867720000025.GIF" wi="169" he="142" />表示对迭代前的衰减系数估计的偏导,<img file="FDA0000388867720000026.GIF" wi="187" he="91" />表示迭代后的时间延迟估计,<img file="FDA0000388867720000027.GIF" wi="132" he="91" />表示迭代前的时间延迟估计,μ<sub>Δτ</sub>表示<img file="FDA0000388867720000028.GIF" wi="76" he="82" />的迭代步长,<img file="FDA0000388867720000029.GIF" wi="234" he="220" />表示对迭代前的时间延迟估计的偏导;60):将步骤50)测算出的迭代后的时间延迟估计和迭代后的衰减系数估计代入式(1),测算出新的估计误差e(t+1);70):比较步骤30)和步骤60)测算的两次估计误差的相对误差,相对误差为<img file="FDA00003888677200000210.GIF" wi="464" he="157" />如果相对误差小于1%,则停止测算,将步骤50)测得的时间延迟估计作为真实时间延迟Δτ,依据式(3),测得转子碰摩声发射源的位置;如果相对误差大于等于1%,则将步骤60)中测算出新的估计误差e(t+1)取代步骤30)中的e(t),然后重复步骤40)至步骤70);S<sub>1</sub>=(D-vΔτ)/2     式(3)其中,S<sub>1</sub>为转子碰摩声发射源到第一传感器(1)的距离,D表示第一传感器(1)和第二传感器(2)之间的距离,v为声发射信号的传播速度。
地址 211103 江苏省南京市江宁区润发路5号