发明名称 一种免疫随机共振的机床传动系统电流特征频率提取方法
摘要 一种免疫随机共振的机床传动系统电流特征频率提取方法,先利用电流传感器采集传动系统伺服驱动电机的电流信号,再计算传动系统中各级输出转轴的回转频率和啮合频率,采用巴特沃斯滤波型带通滤波器处理采集的电流信号,将滤波后的电流信号进行移频,再对移频的电流信号采用多模态免疫克隆选择方法优化二次采样随机共振系统参数,选取最大亲和力的参数作为最优二次采样尺度和结构系统参数,最后根据最优二次采样尺度实现输入信号的压缩,对应频谱中的最大谱峰分量变换后为机床运动系统的特征频率分量,本发明提高了大频率尺度伺服电机电流信号的二次采样随机共振系统优化质量,能够有效提取出传动系统微弱特征频率。
申请公布号 CN102175915B 申请公布日期 2013.01.02
申请号 CN201010594627.6 申请日期 2010.12.17
申请人 西安交通大学 发明人 梁霖;徐光华;刘飞;张熠卓
分类号 G01R23/02(2006.01)I;G06N3/12(2006.01)I 主分类号 G01R23/02(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贺建斌
主权项 1.一种免疫随机共振的机床传动系统电流特征频率提取方法,其特征在于:包括以下步骤:1)利用电流传感器采集传动系统伺服驱动电机的电流信号x(n),n为采样点数,n为大于零的整数;2)计算传动系统中各级输出转轴的回转频率f<sub>r</sub>和啮合频率f<sub>n</sub>,其中,f<sub>r</sub>=f<sub>i</sub>*Z     (1)f<sub>n</sub>=f<sub>r</sub>*m    (2)f<sub>i</sub>为传动系统的输入频率,Z为传动比,m为转轴齿数,则传动系统在电流信号中的特征频率分量为:f=f<sub>m</sub>+f<sub>n</sub>    (3)其中,f<sub>m</sub>为电流主频率;3)采用巴特沃斯滤波型带通滤波器处理采集的电流信号,其中频带选择为[0.8*f,1.2*f]Hz;4)滤波后的电流信号进行移频,将频带下限移至频率0Hz处;5)对移频的电流信号x(n)采用多模态免疫克隆选择方法优化二次采样随机共振系统参数,即二次采样尺度R及结构系统参数a和b,其中a、b为朗之万势阱函数<img file="FDA00001961873700011.GIF" wi="442" he="106" />的系统参数,二次采样尺度R为信号压缩比率R=f<sub>s</sub>/f<sub>sr</sub>,f<sub>sr</sub>为压缩后的频率值,亲和力函数为<img file="FDA00001961873700012.GIF" wi="271" he="132" /><img file="FDA00001961873700013.GIF" wi="570" he="119" />S<sub>p</sub>为频率分量f的功率,<maths num="0001"><![CDATA[<math><mrow><msub><mi>N</mi><mi>p</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>N</mi><mo>=</mo><mo>-</mo><mi>n</mi></mrow><mi>n</mi></munderover><msup><mrow><mo>|</mo><msub><mi>A</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>-</mo><mi>N&Delta;f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><msup><mrow><mo>|</mo><msub><mi>A</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>-</mo><mi>m&Delta;f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow></math>]]></maths>N<sub>p</sub>为输出信号的噪声功率,Δf为频率分辨率,A<sub>p</sub>(w)为第p个抗体表示参数的随机共振输出信号幅值谱;6)在优化出的参数集合中,选取最大亲和力的参数作为二次采样尺度R和结构系统参数a和b;7)根据二次采样尺度R实现输入信号的压缩,压缩采样频率为f<sub>sr</sub>=f<sub>s</sub>/R,f<sub>s</sub>为采样频率,其次,根据结构系统参数输出响应波形,对应频谱中的最大谱峰分量变换后为机床运动系统的特征频率分量,其中变换公式为f=f<sub>c</sub>+R*f<sub>0r</sub>,f<sub>0r</sub>为最大谱峰频率,f<sub>c</sub>是巴特沃斯带通滤波器的下限值。
地址 710049 陕西省西安市碑林区咸宁西路28号