发明名称 确定水力机械振动信号功率谱高低能量区域的方法
摘要 本发明涉及一种确定水力机械振动信号功率谱高低能量区域的方法,在水力机械振动信号的功率谱中,低频高能量区域与高频低能量区域之间存在着一个明显的交界点,将水力机械振动信号功率谱中的低频高能量区域和高频低能量区域能量随振动信号频率的变化趋势采用最小二乘法拟合的一次分段函数来拟合,则偏差平方和最小的一次分段函数的交点即为水力机械振动信号功率谱中的低频高能量区域和高频低能量区域能量分界点。
申请公布号 CN103162815A 申请公布日期 2013.06.19
申请号 CN201210532879.5 申请日期 2012.12.11
申请人 哈尔滨电机厂有限责任公司;哈尔滨大电机研究所;国家水力发电设备工程技术研究中心 发明人 赵越;李正;明君;李屹;王从寅;赵明晗;赵英男;李振龙;辛艳东;郭全宝;刘智良;孙德友;童文波;覃咏麟;魏方睿;吴可君;邵国辉;徐用良;刘登峰;张金伟;赵伟;张颖
分类号 G01H17/00(2006.01)I 主分类号 G01H17/00(2006.01)I
代理机构 代理人
主权项 1.一种确定水力机械振动信号功率谱高低能量区域的方法,其特征是:1、启动计算机系统;2、对采集的振动信号进行功率谱分析;3、在频率分析区间(x<sub>0</sub>,x<sub>n-1</sub>)内,有x<sub>i</sub>=x<sub>0</sub>+i·h,其中:i=0,1,2…,n-1,且h为频率分辨率。令x<sub>m</sub>=x<sub>i+1</sub>,分别建立一次分段函数<img file="FDA00002563996000011.GIF" wi="53" he="43" /><maths num="0001"><![CDATA[<math><mrow><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>k</mi><mn>0</mn></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mn>0</mn></msub><mo>,</mo><mi>x</mi><mo>&le;</mo><msub><mi>x</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><msub><mi>k</mi><mn>1</mn></msub><mi>x</mi><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><mi>x</mi><mo>></mo><msub><mi>x</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中:<img file="FDA00002563996000013.GIF" wi="28" he="43" />表示一次分段函数;x<sub>m</sub>表示一次分段函数<img file="FDA00002563996000014.GIF" wi="27" he="43" />的交点;k<sub>0</sub>、k<sub>1</sub>表示一次分段函数<img file="FDA00002563996000015.GIF" wi="27" he="42" />中一次项的系数;b<sub>0</sub>、b<sub>1</sub>表示一次分段函数<img file="FDA00002563996000016.GIF" wi="27" he="42" />中常数项;4、令振动信号的每一测量数据(x<sub>i</sub>,y<sub>i</sub>)对上述一次分段函数的偏差为V<sub>i</sub>,则有<img file="FDA00002563996000017.GIF" wi="252" he="59" />即:<maths num="0002"><![CDATA[<math><mrow><msub><mi>V</mi><mi>i</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>k</mi><mn>0</mn></msub><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>b</mi><mn>0</mn></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>x</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>k</mi><mn>1</mn></msub><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>b</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>></mo><msub><mi>x</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>若x<sub>i</sub>≤x<sub>m</sub>时有n<sub>1</sub>个测量数据,x<sub>i</sub>&gt;x<sub>m</sub>时有n<sub>2</sub>个测量数据,即n<sub>1</sub>+n<sub>2</sub>=n,则偏差的平方和Q<sub>i</sub>为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>Q</mi><mi>i</mi></msub><mo>=</mo><mi>&Sigma;</mi><msubsup><mi>v</mi><mi>i</mi><mn>2</mn></msubsup><mo>=</mo><mover><mi>&Sigma;</mi><msub><mi>n</mi><mn>1</mn></msub></mover><msubsup><mi>v</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><mover><mi>&Sigma;</mi><msub><mi>n</mi><mn>2</mn></msub></mover><msubsup><mi>v</mi><mi>i</mi><mn>2</mn></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>[</mo><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>k</mi><mn>0</mn></msub><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msub><mi>b</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><msub><mi>n</mi><mn>1</mn></msub></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>[</mo><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>k</mi><mn>1</mn></msub><msub><mi>x</mi><mi>i</mi></msub><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>令<maths num="0004"><![CDATA[<math><mrow><mfrac><msub><mrow><mo>&PartialD;</mo><mi>Q</mi></mrow><mi>i</mi></msub><mrow><mo>&PartialD;</mo><msub><mi>k</mi><mn>0</mn></msub></mrow></mfrac><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mfrac><msub><mrow><mo>&PartialD;</mo><mi>Q</mi></mrow><mi>i</mi></msub><mrow><mo>&PartialD;</mo><msub><mi>k</mi><mn>1</mn></msub></mrow></mfrac><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mfrac><msub><mrow><mo>&PartialD;</mo><mi>Q</mi></mrow><mi>i</mi></msub><mrow><mo>&PartialD;</mo><msub><mi>b</mi><mn>0</mn></msub></mrow></mfrac><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mfrac><msub><mrow><mo>&PartialD;</mo><mi>Q</mi></mrow><mi>i</mi></msub><mrow><mo>&PartialD;</mo><msub><mi>b</mi><mn>1</mn></msub></mrow></mfrac><mo>=</mo><mn>0</mn></mrow></math>]]></maths>以确定相应的k<sub>0</sub>、b<sub>0</sub>、k<sub>1</sub>和b<sub>1</sub>,式中:<img file="FDA00002563996000026.GIF" wi="106" he="118" />Q<sub>i</sub>对k<sub>0</sub>的偏导数;<img file="FDA00002563996000027.GIF" wi="106" he="118" />Q<sub>i</sub>对k<sub>1</sub>的偏导数;<img file="FDA00002563996000028.GIF" wi="105" he="118" />Q<sub>i</sub>对b<sub>0</sub>的偏导数;<img file="FDA00002563996000029.GIF" wi="106" he="119" />Q<sub>i</sub>对b<sub>1</sub>的偏导数;分别计算x<sub>m</sub>=x<sub>i+1</sub>时的Q<sub>i</sub>值;5、对所有Q<sub>i</sub>值进行排序,则偏差平方和Q<sub>i</sub>的最小值Q<sub>imin</sub>所对应的振动信号功率谱低频高能量区域和高频低能量区域一次分段函数<img file="FDA000025639960000210.GIF" wi="28" he="43" />的交点频率值x<sub>B</sub>即为水力机械振动信号功率谱中的低频高能量区域和高频低能量区域能量分界点;6、输出计算结果。
地址 150040 黑龙江省哈尔滨市香坊区三大动力路99号哈尔滨电机厂有限责任公司技术管理部