发明名称 基于深度优化粒子滤波器的桥梁振动频率检测方法
摘要 本发明公开了一种基于深度优化粒子滤波器的桥梁振动频率检测方法,在桥身两侧对称安装加速度传感器,对加速度传感器采集的样本进行分析,计算出k时刻的归一化权值、k+1时刻的归一化权值和k+1时刻统计样本的均值;筛选传感器的样本值,得到k+1时刻粒子滤波后的样本;将滤波后的k+1时刻样本,乘以k+1时刻归一化后的权值,得到此时唯一的样本值;重复上述步骤确定每时刻唯一的样本值,通过曲线拟合进行频谱计算,得到桥梁的振动频率。本发明节约了桥梁检测监控的硬件成本,提高了滤波的精度。
申请公布号 CN102564569B 申请公布日期 2014.04.02
申请号 CN201110456751.0 申请日期 2011.12.26
申请人 西北工业大学 发明人 贾蒙;樊养余;田维坚;马颖;王凤琴
分类号 G01H17/00(2006.01)I 主分类号 G01H17/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 1.基于深度优化粒子滤波器的桥梁振动频率检测方法,其特征在于包括下述步骤:1)在桥身两侧对称安装加速度传感器,加速度传感器安装的间距为5-10米,且保证桥梁两端的桥身两侧四个角处安装加速度传感器,加速度传感器的采样频率大于桥梁最大振动频率的10倍;2)对加速度传感器采集的样本进行分析,根据k时刻的采集样本<![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><mo>{</mo><msubsup><mover><mi>x</mi><mo>~</mo></mover><mi>k</mi><mn>1</mn></msubsup><mo>,</mo><msubsup><mover><mi>x</mi><mo>~</mo></mover><mi>k</mi><mn>2</mn></msubsup><mo>.</mo><mo>.</mo><mo>.</mo><msubsup><mover><mi>x</mi><mo>~</mo></mover><mi>k</mi><mi>N</mi></msubsup><mo>}</mo><mo>,</mo></mrow></math>]]></maths>计算出k时刻的归一化权值<![CDATA[<math><mrow><msub><mover><mi>w</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><mo>{</mo><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>k</mi><mn>1</mn></msubsup><mo>,</mo><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>k</mi><mn>2</mn></msubsup><mo>.</mo><mo>.</mo><mo>.</mo><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>k</mi><mi>N</mi></msubsup><mo>}</mo><mo>,</mo></mrow></math>]]></maths>其中N为传感器个数;计算k时刻的归一化权值的方法包括以下步骤:(1)从采样中得到k时刻的样本:<img file="FDA0000400193770000013.GIF" wi="346" he="88" />(2)分别计算样本均值<![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>,</mo></mrow></math>]]></maths>样本的方差<![CDATA[<math><mrow><mi>R</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mi>k</mi><mi>i</mi></msubsup><mo>-</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths>(3)利用先验概率公式计算样本的权值:<img file="FDA0000400193770000016.GIF" wi="566" he="153" />(4)归一化权值:<![CDATA[<math><mrow><msubsup><mover><mi>w</mi><mo>~</mo></mover><mi>k</mi><mi>i</mi></msubsup><mo>=</mo><mfrac><msubsup><mi>w</mi><mi>k</mi><mi>i</mi></msubsup><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>w</mi><mi>k</mi><mi>i</mi></msubsup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>3)采集k+1时刻的样本<img file="FDA0000400193770000018.GIF" wi="569" he="91" />计算出k+1时刻的归一化权值<![CDATA[<math><mrow><msub><mover><mi>w</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mo>{</mo><msubsup><mover><mi>w</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><msubsup><mover><mi>w</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mn>2</mn></msubsup><mo>.</mo><mo>.</mo><mo>.</mo><msubsup><mover><mi>w</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mo>}</mo></mrow></math>]]></maths>和k+1时刻统计样本的均值<![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mover><mi>x</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>i</mi></msubsup><mo>;</mo></mrow></math>]]></maths>4)筛选传感器的样本值,根据深度优化公式<img file="FDA00004001937700000111.GIF" wi="885" he="210" />判断样本的保留与去除,得到k+1时刻粒子滤波后的样本<img file="FDA00004001937700000112.GIF" wi="568" he="91" />5)将滤波后的k+1时刻样本,乘以k+1时刻归一化后的权值,得到此时唯一的样本值<![CDATA[<math><mrow><msub><mi>x</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mn>1</mn></msubsup><msubsup><mover><mi>w</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mn>1</mn></msubsup><mo>+</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mn>2</mn></msubsup><msubsup><mover><mi>w</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mn>2</mn></msubsup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mover><mi>w</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mo>;</mo></mrow></math>]]></maths>6)重复步骤2)~5),确定每时刻唯一的样本值,得到一组滤波后的样本{x<sub>1</sub>,...,x<sub>k+1</sub>,...x<sub>t</sub>},t表示时间长度,通过曲线拟合进行频谱计算,得到桥梁的振动频率。
地址 710072 陕西省西安市友谊西路127号