随机振动多窗口谱分析法模型建立及统计分析.docVIP

随机振动多窗口谱分析法模型建立及统计分析 针对随机振动中谱估计算法，分析了周期图估计算法对噪声 系列的功率谱相关系数及加窗造成的信息损失，在此基础上分析了由 直接谱估计算法构建的多窗口谱分析法，建立了其噪声模型，推导了 估计方差减小与参与运算的窗口数的关系。针对白噪声序列，得到 MTSA估计器避免加窗而造成信息损失时窗口数与信号序列的关系，最 后推出了 MTSA估计器的统计特征。 关键词：谱估计：周期图估计：多窗谱分析 :ATN911.7 引言 随着现代工业的发展，尤其是国防工业，对设备的可靠性要求越 来越高，如由于飞机的颠簸、爬升、俯冲等会对机载设备造成强烈的 振动、冲击，严重地考验着这些设备的可靠性。振动试验的目的在于 确定所设计的设备在运输、工作过程中能承受外来或者自身产生的振 动而不被破坏。在可靠性振动试验中，随机振动是通过谱均衡控制的 方法来达到的，即通过谱均衡控制使测试件承受的振动功率谱与设定 的参考功率谱一致［1，2］。在随机振动中，由于不同试件的特性参数存 在很大差异，且系统参数不稳定的影响导致整个系统的动态特性发生 改变，因而对谱估计算法提出很高的要求，同时谱估计在电力、化工、 冶金、矿山、船运等行业应用广泛［3］。目前振动信号频谱分析最常 用的方法是基于FFT的周期图法，这种方法的分辨力性能虽然较好， 但方差性能很差，不满足一致性估计条件。为改进周期图法的这一问题2 提出Y很多改进算法，如Welch功率谱估计算法。Welch算法可有效 降低功率谱估计的方差，但分辨力性能同时也降低。为了提高Welch 算法的性能，提出了增加分段数的方法使估计方差减小，但分数段增 加时，段长减小带来分辨力下降，因此，Welch算法的估计方差与分 辨力是一对不可调和的矛盾。从众多文献的研究［4, 5, 6〕知周期图法及 其各种改进算法存在分辨力性能和方差性能之间的固有矛盾。本文针 对随机振动谱估计分析采用多窗口谱分析法(Multitaper Spectrum Analysis, MTSA),并对MTSA谱估计器和传统算法进行对比分析，在此 基础上对其统计特性进行分析以确立其噪声模型，为振动谱均衡控制 提供良好的估计算法基础。 MTSA谱估计器的统计分析及其噪声模型的建立 由于MTSA谱估计器相当于K个直接谱估计器的平均，因此它的 统计特性与直接谱估计器的特性有很大的关联。以下通过白噪声序列 的周期图估计算法的特性分析，建立MTSA谱估计器的噪声模型，并对 其统计特性进行分析。 2.1内噪声序列功率谱周期图估计算法分析 对于零均值的平稳高斯过程，加w(t)窗的Fourier变换为： (1) 直接谱估计器可表示为： (2) 因此，直接谱估计器的协方差为： (3) 其中，相关序列。由此可以得到直接谱估计器在不同频率上的相 关系数力： (4) 当窗函数为矩形窗，即，直接谱估计器演变为周期图谱估计器， 它在Fourier频率点处不相关。对于0噪声序列式(3)变为： (5) 令得： (6) 根据式(4)计算周期图估计器针对白噪声序列的功率谱相关系数。 取，16时，图1为计算得到的结果。由图中可以看出，周期图估计器 的不相关频率问隔与Fourier频率是一致的。从上面的分析可知直接 谱估计器进行离散频率取样时，它的不相关离散频率不再是Fourier 频率，下面针对白噪声序列来分析直接谱估计器系数不相关时的离散取样点。取DPSS窗与MS窗，，，计算结果如图2和图3。从图中 可以看出，使用零阶窗函数时，在谱窗带宽外直接谱估计的结果不相 关；随着谱窗阶数的增加，相关的频率点扩展到谱窗带宽之外；然而 由图2 (c，d, e，f)、图3 (c, d，e, f)可以看到在谱窗带宽之内，某些频率 点仍然不相关；通过比较0频率点、0. 25频率点与f的相关系数，表 明相关系数的分布与频率的测量点无关。 3.2直接谱估计信息损失分析及MTSA谱估计器噪声模型 由于MTSA谱估计器相当于K个直接谱估计器的平均，因此它的统 计特性与直接谱估计器的特性有很大的关联。设零均值的白噪声过 程，其方差为，因此其功率谱密度为，方差的估计值为，并可求 得的方差为 (7) 从理论上来说，一个最好的谱估计器得到的功率谱密度为。对于 白噪声过程，周期图谱估计器是的无偏估计，因此对进行平滑可得 到： (8) 对于使用归一化窗函数w(t)的直接谱估计器，对它进行频域平 滑： (9) 并存方差为：由cauchy不等式可以得到，且只有当w（t）为矩形窗时等号成立， 因此对于任何非矩形窗，周期图估计器与直接谱估计器关系如下： （11） 从式（11）可以看出虽然直接谱估计器使估计的结果得到平滑， 但使用窗函数对信号序列进行加权，导致信息量的减小，从而使估计 的方差增大。下面分析MTSA谱