MATLAB在机械振动信号中的应用

(完整版)MATLAB在机械振动信号中的应用申振（山东理工大学交通与车辆工程学院）摘要：综述现代信号分析处理理论方法如时域分（包括时域参数识别、相分等域分析（包括傅立叶变换功率谱分解)并合MATLAB中的相关函数来对所拟合的振动信号进时分析和频域分析并绘的频谱图进行说明关键词时域分析频域分析MATLAB信号是信息的载体，采用合适的信号分析处理方法以获取隐藏于传感观测信号中的重要信息括时域与频域信息等对于许多工程应用领域均具有重要意义。对获取振动噪声信号的分析处理，是进行状态监测、故障诊断、质量检查、源识别、机器产品的动态性能测试与优化设计等工作的重要环节，它可以预先发现机械部件的磨损和缺陷等故障，从而可以提高产品的质量，降低维护费用。随着测试技术的迅速发展,种信号分析方法也随之涌现,并广泛应用在各个领域［1]。时域描述简单直观，只能反映信号的幅值随时间的变化,而不能明确的揭示信号随时间的变化关系。为了研究信号的频率组成和各频率成分的幅值大小、相位关系，应对信号进行频谱分析把时域信号通过适当的数学方法处理变成频率

f

(或频率）为独立变量，相应的幅值或相位为因变量的频域描述。频域分析法将时域分析法中的微分或差分方程转换为代数方程，有利于问题的分析[2

.MATLAB是MathWorks公司于1982年出的一种功能强效高交互性好的数值计算和可视化计算机高级语言,它将数值分析、矩阵运算、信号处理和图形显示有机地融合为一体形成了一个极其方便、用户界面良好的操作环境。随着其自身版本的不断提高，MATLAB功能越来越强大应用范围也越来越广，如广泛应用于信号处理、数字图像处理、仿真、自动化控制、小波分析及神经网络等领域[3]

.本文主要运用了MATLAB对机械振信号进行分析分析过程包括时域分析和频域分析两大部分，时域分析的指标包括随机信号的均值、方差以及均方值。频域分析的性能指标包1

T12Txx(完整版)MATLAB在机械振动信号中的应用T12Txx括对功率谱分析、倒频谱分析。在进行上述分析之前先要对振动信号进行拟合.机械振动分为确定性振动和随机振动确定性振动又分为周期振动和非周期振动，周期振动又进一步分为简谐振动和复杂的周期振动。所以可以根据上述的分类来拟合振动信号[

.在计信号的处理程序时，运用MATLAB中的相关函数来对所拟合的振动信号进行时域分析和频域分析，并对绘出的频谱图进行说明。1

时域分析1.1均值

对于一个各态历经随机随机信号(t)

，其均值为x

x

T

1

0

x(t

（1）式中xt)

——样本函数T

-—测时间；—-常值分量。x1.2方差是述随机信号的波动分量定义为xlimt)]T0

（1）它表示信号(t)

偏离其均值平方的均值，方差的正平方根称为标准差。xx1均值

x

是随机信号(t)

平方的平均值，定义为

x

2

limT

1

T0

()

(3)它描述信号的能量或强度是(t)

平方的均值方值的正平方根值称为均根值

rms

参数、、之间的关系为xxx

2x

-

x

2

（4）1.4

时域统计分析

概率密度分析是以幅值大小为横坐标，以每个幅值间隔内出现的概率为纵坐标进行统计分析的方法。它反映了信号落在不同幅值强度区域内的概率情况。计算方法如下：(0

P[xx(

lim

/T

(0

i

)i

（5)2

n2nnn2nn概率密度函数()

(完整版)MATLAB在机械振动信号中的应用给出了信号取不同幅值大小的概率是随机信号的主要特征参数之一不同的随机信号有不同的概率密度函数图形可以此来识信号的性质，如正弦信号加随机噪声、窄带随机信号及宽带随机信号等。概率分布函数是信号幅值小于或等于某值R的概率定义为：F)

)dx

(6）概率分布函数又称为累积概率函数,表示了信号幅值落某一区间的概率［］

.22.1

频域分析傅里叶变换任何周期函数，均可展开成正交函数线性组合的无穷级数，如三角函数集的傅里叶级数。叶级数的表达形式如下：a(t)Asin(nt0n

（7)A022nnarctann

（n

（8）对于非周期信号或瞬变信号，利用如下的傅立叶变换进行频谱分析：

)t)ejdtxt)X)ej

2.2

功率谱分析2

经典功率谱估方法若xt)

为平稳随机信号自相关函数为绝对可时相关函数R)和功率谱密度S()xx为一个傅里叶变换对即3

K(完整版)MATLAB在机KS(1((e

（10)同理，在频域描述两个随机信号xt)

和(t)

相互关联程度的数字特征，可以定义为互谱功率密度简称互谱密度。而且互相关函数与互谱密度是一个傅里叶变换对.S(1((2

（11）2.2。2

改进的直接估法直接法和间接法的方差性能很差而且当数据长度太大时,曲线起伏加剧；若数据长度太小，则谱的分辨率又不好所以需要改进[提高的周期图法估计的另一种方式就是采用对采样数据分段使用非矩形窗,即Welch法。由于非矩形窗在边沿趋近于零，从而减少了分段对重叠的依赖。选择合适的窗函数，采用每段一半的重叠率能大大降低谱估计的方差。这种方法中，记录数据仍分成

K

段，即x(i)()x(N)

0K

（12每段M个取。窗函数w()

在计算周期图之前就与数据段相乘，于是可定义K

个修正周期图J(i))M

1

Mn

x(i)(n)(e

2

iK

（13)U

是窗口序列函数的平均能量U

1

Mn

w2(n)

(14）则定义谱估计为w

1K

i

(iM

（15）2.2.3AR模型功率谱估计法传统的功率谱估计方法是利用加窗的数据或加窗的相关函数估计值的傅里叶变换来计算4

q(完整版)MATLAB在机械振动信号中的应q的，具有一定的优势，如计算效率高，估计值正比于正弦波信号的功率等。但是同时也存在许多缺点，主要缺点就是方差性能较差、谱分辨率低。而参数模型法可大大提高功率谱估计的分辨率是现代谱估计的主要研究内容语音分析据压缩以及通信等领域有着广泛的应用[按照模型化进行功率谱估计主要思想如下：（1）选择模型；（2）从给出的数据样本估计假设的模型;（3）将估计的模型参数打入模型的理论功率谱公式中得出一个较好的谱估计值。假设产生随机序列x(n

的系统模型为一个线性差分方程即x(n((j)ijij

（16式中，(n)

表示白噪声序列，对上式进行Z变换可得

a()zj

bW(z)i

（17)j

i所以系统的传递函数为H(z)

X(z)W(z)

(z)(

（18式中，q(z)zjj

j

(19）q(z)bzii

（20）假定输入白噪声功率谱密度为(z),那么输出功率谱密度为ww(

(zB((z)A

))

(21）又根据zj，所以得5

(完整版)MATLAB在机械振动信号中的应用Px

2

(ej(j)

2

（22这样,当确定了系数a、和2后就可以求解得到随机信号的功率谱密度了通过上可jiw知，如果ii

时，则系统的差分方程变为qx(nx(nj)(njj上式即为自回归模型简称为AR（Auto—Regressive）模型，再将该式进行变换得

（23）H(

()z)Az)

zj

（24)j所以，AR型称为全极点模型AR型输出功率谱为

(j

)

qaj

j

2

（25)j显然，计算出2和后,就可以求解得到随机信号的功率谱密度p。wj本文采用型的一种Burg法进行功率谱估计。3

仿真研究仿真带噪声信号如下：(t)tt)ett)randn(该仿真带噪声信号由两个正弦信号e

t

ft1

e

t

ft2

和一个服从正态分布的高斯白噪声信号(t

叠加而成。f100Hzf3001

。其时域波形如图1所（程序详见附录16

(完整版)MATLAB在机械振动信号中的应用图1时域波形时域分析结果：序列的平均值为0。5050序列的最小值为-10.7448序列的最大值为12。0222序列的标准差为2。9153序列的方差为8。4992序列的均方值为2.9580图2

经典功率谱估图在功率谱中可以很明显的看到振动信号中有100Hz和两个主要的频率。表明信号中含有这两个频率的周期成分。如图27

(完整版)MATLAB在机械振动信号中的应用图3FFT频图上图3为频谱图从该谱中可以看到有三个主要峰值,即在0Hz，300Hz处。用Burg法进行PSD估计率谱图如图4，从中可以很明显的看到振动信号中有100Hz和300Hz个主要的频率.表明号中含有这两个频率的周期成分（程序详见附录图4Burg法进行PSD估计功率谱图在Welch法进行PSD功率谱估计当采用不同窗函数时的结果。从中可以很明显的看到振动信号中有100Hz和300Hz两个主要的频率。表明信号中含有这两个频率的周期成分。且海宁窗和布莱克曼窗较为明显（程序详见附录8

(完整版)MATLAB在机械振动信号中的应用图5Welch法行PSD功率谱估计功率谱图图5倒谱图理论上，傅立叶变换用于频谱分析，可以找出受噪声干扰的信号的频率成分，而这用时域分析是不能分辨的。对傅立叶变换做复共轭运算,可得到信号的功谱密度函数以显示各频率分量的能量分布。仿真带噪信号的傅立叶变换与功率谱分解结果如图和图、5、6所示。从图3和图4、6可以清楚看到，约在频率为100Hz、300Hz(即振动信号频率的倍频）处频谱幅值和能量出现局部极大值，对应机械振动的主振动源所在。4

结论信号是信息的载体，因此采用合适的信号分析处理方法以获取隐藏于传感观测信号中的重要信息包括时域与频域信息等于许工程应用领域均具有重要意义本在研究现代信号分析处理理论方法如时域分析包括时域参数识别相关析以及统计分析等域分析(包括傅立叶变换、功率谱分解等）的基础上，结合仿真数据对机械振动信号分析处理，具有一定9

(完整版)MATLAB在机械振动信号中的应用的参考价值。参考文献[1]凯工测试技术［M].西安:北业大学出版社，［2］许同乐。机械工程测试技术。北京：机械工业出版社,2010.［3]薛年喜MATLAB在数字信号处理中的应用（第二版）北京：清华大学出版社,2008［4]焦卫东。旋转机械振动信号分析。浙江.嘉兴学院学报.2007.附录附录一时域分析、频域分析程序A1=6；A2=8f1=100;f2=300；fs=1000t=0:1/fs:2N=length(t）;X1=A1*exp（-1.0*t).＊sin＊pi*f1X2=A2*exp（—2。0＊t)＊f2＊t)；R=randY=X1+X2+R；figure(1)；plot(t，Y);title（’振动信号的波形’xlabel('时间/秒')；ylabel（’幅度grid;on；％时域分析mi=min(Y);disp(mi）;%最小值mx=max（Y)；disp（mx)；％最大值st=std(Y)；disp标准差m=meandisp（m);％均值vr=var（Ydisp；％方差rm=rmsdisp(rm);％均方值10

(完整版)MATLAB在机械振动信号中的应用％频域分析l=length（Yr=fft(Y)/l;r=fftshift(rf=linspace(—fs/2,fs/2,l)figure（2，abs(r)）;gridholdon;figure(3）;psd（Y,2048，kaiser(512，5),0,0.95figure(4)；yc=rcepsplot（yc附录二Welch方法进行估计程序A1=6；A2=8f1=100；f2=300；fs=1000；nfft=1024；t=0：1/fs；N=length(t)X1=A1＊exp（—1.0。*sin(2*pi＊f1＊tX2=A2*exp（—1。5*t)＊sin（2＊piR=rand（1，NY=X1+X2+R；window1=boxcar(100)；window2=hamming（10