发明名称 一种去除振动观测信号中背景振动的有效方法
摘要 本发明涉及一种去除振动观测信号中背景振动的有效方法,属于信号技术领域。本发明从考虑本底振动与观测振动互相关性出发,提出了采用自互功率谱法去除本底振动的新方法,推导了计算公式。本发明中的自互功率谱法能够更有效地去除观测数据中的本底干扰,获得更贴近实际的交通环境振动的功率谱、时程、振动级以及加权振级。在本底振动占优势的低频段,本发明较现有方法计算精度有明显的提高。
申请公布号 CN103528671A 申请公布日期 2014.01.22
申请号 CN201310471586.5 申请日期 2013.09.30
申请人 黑龙江八一农垦大学 发明人 郑鑫;葛建锐;解恒燕;陶夏新;王福彤
分类号 G01H17/00(2006.01)I 主分类号 G01H17/00(2006.01)I
代理机构 代理人
主权项 1.一种去除振动观测信号中背景振动的有效方法,其特征在于:采用以下计算步骤:观测得到的环境振动由交通的环境振动和本底振动两部分组成,可以表达为h<sub>T</sub>(t)=h<sub>O</sub>(t)-h<sub>B</sub>(t)   (1);式中,h<sub>T</sub>(t)为交通系统引起的环境振动时程,h<sub>O</sub>(t)为观测记录总时程,h<sub>B</sub>(t)为本底时程。h<sub>O</sub>(t)是观测得到的,h<sub>B</sub>(t)也可以通过观测没有车辆通过时测点的振动获得;对(1)式两边做Fourier变换,得H<sub>T</sub>(f)=H<sub>O</sub>(f)-H<sub>B</sub>(f)   (2);其中,H(f)表示h(t)的Fourier谱。对(2)式两边求幅值的平方,<maths num="0001"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><mo>=</mo><mo>[</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><msup><mrow><mo>[</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mo>*</mo></msup></mtd></mtr><mtr><mtd><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msubsup><mi>H</mi><mi>B</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>H</mi><mi>O</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>式中,[·]*表示共轭,由于<img file="FSA0000096055620000012.GIF" wi="260" he="54" />与<img file="FSA0000096055620000013.GIF" wi="258" he="54" />互为共轭,则<img file="FSA0000096055620000014.GIF" wi="561" he="55" />可表示为<img file="FSA0000096055620000015.GIF" wi="407" he="54" />据此,(3)式可写成<maths num="0002"><![CDATA[<math><mrow><mrow><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>2</mn><mi>Re</mi><mo>[</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msubsup><mi>H</mi><mi>B</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>对(4)式求期望,<maths num="0003"><![CDATA[<math><mrow><mi>E</mi><mo>[</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>=</mo><mi>E</mi><mo>[</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>+</mo><mi>E</mi><mo>[</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>-</mo><mn>2</mn><mi>E</mi><mo>{</mo><mi>Re</mi><mo>[</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msubsup><mi>H</mi><mi>B</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>上式两边除以振动持时T得<maths num="0004"><![CDATA[<math><mrow><mfrac><mn>1</mn><mi>T</mi></mfrac><mi>E</mi><mo>[</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>T</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>=</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><mi>E</mi><mo>[</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>+</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><mi>E</mi><mo>[</mo><msup><mrow><mo>|</mo><msub><mi>H</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>-</mo><mfrac><mn>2</mn><mi>T</mi></mfrac><mi>E</mi><mo>{</mo><mi>Re</mi><mo>[</mo><msub><mi>H</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msubsup><mi>H</mi><mi>B</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>]</mo><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>据定义,可以写成P<sub>T</sub>(f)=P<sub>O</sub>(f)+P<sub>B</sub>(f)-2Re[P<sub>OB</sub>(f)]    (7);式中,<img file="FSA0000096055620000019.GIF" wi="437" he="94" />分别为三条时程的自功率谱密度函数,<img file="FSA00000960556200000110.GIF" wi="517" he="94" />为本底振动与观测振动的互功率谱密度函数;(7)式表明,交通环境振动的功率谱可以通过观测振动的功率谱加上本底振动的功率谱,再减去观测振动与本底振动的互功率谱实部的2倍获得。
地址 163319 黑龙江省大庆市高新区新阳路2号
您可能感兴趣的专利