一种基于周期图法的冲床噪声功率谱估计改进方法

摘要

基于周期图法的冲床噪声功率谱估计改进方法，包括：初始化功率谱检测设备，设定传感器采样频率、离散傅里叶变换参数、窗函数形状、窗函数长度、迭代次数计数器i初始值、总迭代次数S、平均离散傅里叶变换近似值初始值；将冲床冲裁控制信号作为开始采样触发信号；将冲床收回液压锤控制信号作为结束采样触发信号；等待触发信号，触发传感器结束采集冲裁噪声x(t)的样本；对采集得到的噪声样本加窗处理；用迭代方法计算L点离散傅里叶变换平均值存储离散傅里叶变换平均值结果；重复以上步骤直到迭代计数器i=S；计算离散傅里叶变换平均值的平方除以数据长度，得到功率谱估计值

说明书

技术领域

本发明涉及一种用于噪声控制领域的功率谱估计方法，具体是一种基于周期图法的冲床噪声功率谱估计改进方法。

背景技术

在工业生产过程中，具体的在冲床的工作过程中，冲床撞击材料会产生一种冲裁噪声。这种噪声具有：重复性、短时性、高强度性的特征。这种重复撞击噪声会导致机器设备的声疲劳，长期作用将会缩短其使用寿命，甚至发生生产事故。强烈的噪声极易形成差拍型次声波，作用与人的躯体。人体各个部位都存在固有频率，身体为7-13HZ，内脏为4-6HZ，头部为8-12HZ，这些固有频率刚好在次声波频带内，所以冲压工人在强烈噪声环境中工作，常有头昏脑胀、恶心和心悸之感。降低冲床噪声已成为噪声控制工程中的当务之急。

不论是用传统的无源消噪技术还是新型的有源消噪技术都需要对噪声进行检测，为噪声控制提供噪声的先验信息。其中最主要的信息是噪声的功率谱信息。功率谱信息能反映出噪声所含有的主要频率成分，以及各个频率成分的大小。传统的无源消噪技术对功率谱信息的依赖不是很强，部分新型的有源消噪技术需要更多的噪声信息。所以功率谱估计方法的精度直接影响这类依赖噪声先验信息的新型有源消噪技术的性能。

近几十年，已有许多学者提出了各种经典的功率谱估计方法并对其进行了深入的研究，取得了一些重要的成果。经典功率谱估计方法中有一种方法名为周期图法，但是这种方法应用到具有重复、短时、高强度噪声背景下存在一些局限：

1.由于周期图功率谱估计方法先求噪声的傅里叶变换，而实际中测得的噪声是离散的有限长的，所以实际中计算的是离散有限长噪声的离散傅里叶变换，这个值是一个近似值。在对短时噪声进行功率谱估计时，采样数据太少，离散傅里叶变换的计算误差会很大，这个误差会直接影响功率谱估计的误差。

2.周期图功率谱估计方法应用在短时高强度（幅值变化剧烈）噪声中严重依赖噪声检测设备的采样频率。只有噪声检测设备达到足够高的采样频率才能有效测得这类噪声的功率谱。而提高设备采样频率的成本高昂。

3.周期图功率谱估计方法不能利用噪声的重复性这一重要的先验信息。

冲床作业中产生的冲裁噪声就是一类重复、短时、高强度噪声，用周期功率谱估计方法无法有效测得冲裁噪声的功率谱信息。如何将噪声的重复性信息利用起来，在保持噪声检测设备采样率不变的前提下，提高功率谱估计性能成为冲床噪声控制工程需要解决的一个问题。

发明内容

本发明要克服现有周期图功率谱估计方法在处理重复、短时、高强度噪声时的不足，提出一种基于周期图法的冲床噪声功率谱估计改进方法。

改进方法首先利用加窗法截取多段有效噪声，然后求取每段噪声的离散傅里叶变换，再计算离散傅里叶变换的平均值，结合周期图法求得噪声功率谱。该方法提高了离散傅里叶变换的估计精度，间接提高了功率谱估计的精度。该方法主要针对具有重复、短时、高强度噪声的功率谱估计，除了能有效应用于冲床的冲裁噪声功率谱估计，还适用于其它具有重复、短时、高强度噪声的功率谱估计，如打桩机、锻造机、射击靶场的噪声的功率谱估计。用该方法改进传统功率谱估计设备不需要改变硬件设备，只需要更新软件中的计算方法，成本低。

本发明是通过以下技术方案实现的，本发明在周期图功率谱估计方法的基础上，根据冲床的冲裁噪声的特性改进噪声的离散傅里叶变换值的计算方法，提高检测设备对冲裁噪声的功率谱估计精度。冲床的冲裁噪声具有重复性和短时性，所以本发明用以下数学公式描述冲裁噪声信号：

x(t)=s(t)+u(t),t∈[0,∞)

T₁=ξT,(ξ≤1)

其中x(t)表示含有高斯白噪声干扰的冲裁噪声信号，s(t)表示冲裁噪声，u(t)表示高斯白噪声干扰，t表示时间，T₁表示一次冲裁噪声有效时长，T表示冲裁周期，ξ表示冲裁噪声占空比。

本发明方法的具体步骤如下：

(1)初始化功率谱检测设备。设定传感器采样频率、离散傅里叶变换参数、窗函数形状、窗函数长度、迭代次数计数器i初始值、总迭代次数S、平均离散傅里叶变换近似值初始值。

(2)将冲床冲裁控制信号作为开始采样触发信号。等待触发信号，触发传感器开始采集冲裁噪声样本序列。

(3)将冲床收回液压锤控制信号作为结束采样触发信号。等待触发信号，触发传感器结束采集冲裁噪声x(t)的样本。

(4)对采集得到的噪声样本加窗处理。默认选择长度为N的矩形窗，也可选择改变窗的长度和形状。比较有效的窗还有汉明窗和布莱克曼窗。加窗的具体做法是对步骤(3)中得到的噪声样本进行截取或补零，样本长度大于N则截取长度为N的样本序列，样本长度小于N则在样本序列末尾补零；然后将样本序列与窗函数序列做点乘，得到一次冲裁噪声样本$>\{{\hat{x}}_{\left(i\right)}\left(n\right),n=\mathrm{0,1},...,N-1\};$>

(5)用迭代方法计算L点离散傅里叶变换平均值迭代更新公式如下：

$>\overline{X_N\left(k\right)}=\left(\begin{array}{cc}{X}_{\left(1\right)}\left(k\right)& i=1\\ \frac{(i-1)\overline{X_N\left(k\right)}+X_{\left(i\right)}\left(k\right)}{i}& i>1\end{array}\right)$>

$>X_{\left(i\right)}\left(k\right)=\underset{n=0}{\overset{N-1}{\Sigma }}{\hat{x}}_{\left(i\right)}{e}^{-j2\mathrm{\pi kn}/L}$>

其中，i表示当前迭代计数器次数，k表示频率，X_(i)(k)表示对第i次冲裁噪声的离散傅里叶变换值。

(6)存储离散傅里叶变换平均值结果，迭代计数器加1，i=i+1。

(7)重复步骤(2)到(6)直到迭代计数器i=S。

(8)计算离散傅里叶变换平均值的平方除以数据长度，得到功率谱估计值公式如下：

$>{\hat{P}}_{\mathrm{xx}}\left(k\right)=\frac{1}{N}{\left|\overline{X_N\left(k\right)}\right|}^2$>

在冲床噪声检测控制工程中采用本发明提出的方法能够获得足够的功率谱估计精度和分辨率，能抑制白噪声的干扰。

本发明的优点是：用加窗法和平均法计算冲裁噪声离散傅里叶变换近似值，解决了传统方法因离散傅里叶变换近似值方差大导致的功率谱估计精度低的缺陷，且方法简单，易于实现。

附图说明

图1为本发明方法的程序流程图。

图2为本发明实施例中冲裁噪声十个周期内的检测图。

图3为本发明实施例中第一次冲裁噪声的检测图。

图4为本发明实施例中未改进方法和改进方法得到的功率谱比较图。

具体实施方式

以下结合附图和实施例对本发明的技术方案作进一步描述。

如图1所示，功率谱检测设备首先初始化设备参数，然后等待冲床冲裁控制信号，触发传感器采集噪声样本序列，当冲床收回液压锤控制信号发出时，触发传感器停止采集噪声样本，这样就完成了冲裁噪声截取和加矩形窗的步骤。加窗程序不做任何处理就是选择矩形窗，也可根据噪声的特点选择非矩形窗，进一步提高算法性能。这里功率谱检测设备使用了迭代的方法计算离散傅里叶变换平均值，不需要存储多次冲裁噪声，节省内存空间。最后计算离散傅里叶变换平均值的平方除以数据长度得到冲裁噪声的功率谱估计值。

如图2所示，冲裁噪声周期为1秒，一次冲裁噪声持续时间为0.1秒，信噪比为10dB。以此为实施例，本发明的冲裁噪声功率谱估计流程如下：

(1)初始化功率谱检测设备。设定传感器采样频率为40KHz，离散傅里叶变换数据点数为4096，窗函数选择长度为4000的矩形窗，迭代次数计数器i为1，总迭代次数为10，平均离散傅里叶变换近似值为0。

(2)将冲床冲裁控制信号作为开始采样触发信号。等待触发信号，触发传感器开始采集冲裁噪声样本序列。

(3)将冲床收回液压锤控制信号作为结束采样触发信号。等待触发信号，触发传感器结束采集冲裁噪声样本。

(4)对采集得到的噪声样本加窗处理。默认选择长度为4000的矩形窗，也可选择改变窗的长度和形状。比较有效的窗还有汉明窗和布莱克曼窗。加窗的具体做法是对步骤(3)中得到的噪声样本进行截取或补零，样本长度大于4000则截取长度为4000的样本序列，样本长度小于4000则在样本序列末尾补零。然后将样本序列与窗函数序列做点乘，得到一次冲裁噪声样本如图3所示为第一次冲裁噪声的检测图。

(5)用迭代方法计算离散傅里叶变换平均值迭代更新公式如下：

$>\overline{X_N\left(k\right)}=\left(\begin{array}{cc}{X}_{\left(1\right)}\left(k\right)& i=1\\ \frac{(i-1)\overline{X_N\left(k\right)}+X_{\left(i\right)}\left(k\right)}{i}& i>1\end{array}\right)$>

$>X_{\left(i\right)}\left(k\right)=\underset{n=0}{\overset{3999}{\Sigma }}{\hat{x}}_{\left(i\right)}{e}^{-j2\mathrm{\pi kn}/4096}$>

其中，i表示当前迭代计数器次数，k表示频率，X_(i)(k)表示对第i次冲裁噪声的离散傅里叶变换值。

(6)存储离散傅里叶变换平均值结果，迭代计数器加1，i=i+1。

(7)重复步骤(2)到(6)直到迭代计数器i=10。

(8)计算离散傅里叶变换平均值的平方除以数据长度，得到功率谱估计值计算公式如下：

$>{\hat{P}}_{\mathrm{xx}}\left(k\right)=\frac{1}{4000}{\left|\overline{X_N\left(k\right)}\right|}^2$>

结果显示在图4中，其中实线是未改进算法对冲裁噪声的功率谱估计结果，虚线是本发明的改进算法对冲裁噪声的功率谱估计结果。虚线的波峰明显，通过虚线可以测出噪声中含有9个主要频率分量：130Hz,290Hz,400Hz,500Hz,611Hz,772Hz,810Hz,881Hz,1000Hz。