Matlab振动程序-代码作业.docxVIP

一、课题任务要求 随着机械工业不断向自动化、高精度、智能化等方向的发展，在机械设备运行及生产过程中进行参量测试、分析与诊断等处理过程已成为必要环节，许多信号处理方法如时域统计分析、相关分析、相干分析、频谱分析等已经被广泛被应用与机械工程测试领域。 本文为机械测试信号的时域和频域分析，其中时域分析包括对信号最大值、最小值、中值、方差的分析，频域分析包括FFT分析、频谱分析、倒谱分析。在进行上述分析之前先要对振动信号进行拟合。机械振动分为确定性振动和随机振动，确定性振动又分为周期振动和非周期振动，周期振动又进一步分为简谐振动和复杂的周期振动。所以可以根据上述的分类来拟合振动信号。在设计信号的处理程序时，用MATLAB中的相关函数来对所拟合的振动信号进行时域分析和频域分析，并对绘出的频谱图进行说明。 二、技术路线 对机械振动信号的时域和频域采集，根据所拟合的振动信号，选取所需要的时域性能指标和频域分析的性能指标对振动信号进行分析。其中时域分析包括对信号最大值、最小值、中值、方差的分析，频域分析包括FFT分析、频谱分析、倒谱分析。 现构造一个振动信号（在该程序中以两个衰减振动分量和一个随机数rand之和来拟合振动信号），再利用MATLAB中的函数mean（）、min（）、max（）、std（）对离散序列中的平均值、最大值、最小值、标准差等时域性能进行分析，通过调用函数fft（y）；psd（y）；rcep（y）对该振动信号进行频域内的性能分析。 在设计过程中的理论知识有离散傅立叶变换（DFT）、功率谱的概念和意义以及倒谱的概念和意义。 = 1 \* GB3 ①DFT的定义和意义： DFT的定义式为： X X(k)=DFT[x(n)]= 0≤k≤N-1 X(n)==x(n) 0≤n≤N-1 DFT的意义：DFT的意义在于它表示信号中的各个频率的分量的振动幅值的大小，亦即该分量对于振动信号影响的大小。通过DFT的快速算法FFT可以很方便的将振动信号的各个分量的幅值比重计算出来。即可以把信号的主频分量提取出来。 = 2 \* GB3 ②功率谱的概念和意义： 功率谱的定义式为:若X(Ω)=DFT[x(m)]，x(n)为N点序列。则 X（Ω） =DFT[xN (-m)] 从而有 DFT[R(M)]= DFT[x(m)] DFT[xN (-m)] 即 (Ω)= X(Ω)X(Ω)=|X(Ω)|^2 综上所述，先用FFT求出随机离散序列的DFT，再计算幅频特性的平方，再除以N，即得到该随机信号的功率谱估计。功率谱的意义在于它可以对信号中的周期成分进行分析。 = 3 \* GB3 ③倒谱的概念与意义： 倒谱可以分析复杂频谱图上的周期成分，分离和提取在密集泛频信号中的成分。倒频谱对于整个谱的形状不敏感。 三、实现程序 A1=3;A2=4; f1=100;f2=250;fs=1000; t=0:1/fs:2; N=length(t); X1=A1*exp(-0.9*t).*sin(2*pi*f1*t); X2=A2*exp(-1.5*t).*sin(2*pi*f2*t); R=rand(1,N); Y=X1+X2+R; figure(1); plot(t,Y); title(振动信号的波形); xlabel(时间/秒); ylabel(幅度); grid; hold on;%时域分析 m= mean(Y); disp (m); mi = min(Y); disp (mi); mx = max(Y); disp (mx); st = std(Y); disp (st);%频域分析 l=length(Y); r=fft(Y)/l;r=fftshift(r); f=linspace(-fs/2,fs/2,l); figure(2); plot(f,abs(r)); grid; hold on; figure(3); psd(Y,2048,1000,kaiser(512,5),0,0.95); figure(4); yc=rceps(Y); plot(yc); 四、运行结果 (1)时域分析结果： 序列的平均值为 0.5038 序列的最小值为 -6.7750 序列的最大值为 7.5083 序列的标准差为 1.6180 (2)频域分析结果: 这是所拟合得到的振动信号的图像。 上图为FFT频谱图，从该频谱中可以看到有三个主要高峰值，即在0Hz，100Hz，250Hz处。 在功率谱中可以很明显的看到振动信号中有100Hz和250Hz两个主要的频率。表明信号中含有这两个频率的周期成分。 五、总结 通过使用MATLAB中的相关函数编写的程序对这一所拟合的振动信号进行了时域分析和频域的分析