《Matlab振动程序-代码作业》由会员分享，可在线阅读，更多相关《Matlab振动程序-代码作业(8页珍藏版)》请在人人文库网上搜索。

1、课题任务要求随着机械工业不断向自动化、高精度、智能化等方向的发展，在 机械设备运行及生产过程中进行参量测试、 分析与诊断等处理过程已 成为必要环节，许多信号处理方法如时域统计分析、相关分析、相干 分析、频谱分析等已经被广泛被应用与机械工程测试领域。本文为机械测试信号的时域和频域分析， 其中时域分析包括对信 号最大值、最小值、中值、方差的分析，频域分析包括FFT分析、频谱分析、倒谱分析。在进行上述分析之前先要对振动信号进行拟合。 机械振动分为确定性振动和随机振动， 确定性振动又分为周期振动和 非周期振动， 周期振动又进一步分为简谐振动和复杂的周期振动。 所 以可以根据上述的分类来拟合振动信号。在。

2、设计信号的处理程序时， 用MATLAB的相关函数来对所拟合的振动信号进行时域分析和频域 分析，并对绘出的频谱图进行说明。二、技术路线对机械振动信号的时域和频域采集， 根据所拟合的振动信号， 选 取所需要的时域性能指标和频域分析的性能指标对振动信号进行分 析。其中时域分析包括对信号最大值、最小值、中值、方差的分析， 频域分析包括FFT分析、频谱分析、倒谱分析。现构造一个振动信号 (在该程序中以两个衰减振动分量和一个随机数rand之和来拟合振动信号)，再利用MATLAB的函数mean()、 min ()、max()、std ()对离散序列中的平均值、最大值、最小值、 标准差等时域性能进行分析，通过。

3、调用函数fft (y)； psd (y)； rcep(y)对该振动信号进行频域内的性能分析。在设计过程中的理论知识有离散傅立叶变换(DFT、功率谱的概念和意义以及倒谱的概念和意义。 DFT的定义和意义:DFT的定义式为：X(k)=DFTx( n)二N 1knx(n)Wn 0N0 W k WN -11 N 10W nWN-1X( n)二X(k)WNkn =x (n)N k 0DFT的意义：DFT的意义在于它表示信号中的各个频率的分量的振动幅值的大小，亦即该分量对于振动信号影响的大小。 通过DFT的 快速算法FFT可以很方便的将振动信号的各个分量的幅值比重计算 出来。即可以把信号的主频分量提取出来。

4、。 功率谱的概念和意义: 功率谱的定义式为:若X( Q) =DFTx(m) , x(n)为N点序列。则X (Q) =DFTx N (-m)1从而有 DFTR(M)二亍 DFTx(m) DFTx N (-m)A丄丄即 SxN( Q )= N X( Q )X*( Q )= N |X( Q )|A2综上所述，先用FFT求出随机离散序列的DFT再计算幅频特性 的平方，再除以N,即得到该随机信号的功率谱估计。功率谱的意义 在于它可以对信号中的周期成分进行分析。 倒谱的概念与意义：倒谱可以分析复杂频谱图上的周期成分，分离和提取在密集泛频信号中的成分。倒频谱对于整个谱的形状不敏感。三、实现程序A1=3;A2。

5、=4;f1 = 100;f2=250;fs=1000;t=0:1/fs:2;N=le ngth(t);X1二A1*exp*t).*si n(2*pi*f1*t);X2二A2*exp*t).*si n(2*pi*f2*t);R=ra nd(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);d。

6、isp (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,;figure(4);yc=rceps(Y);plot(yc);四、运行结果(1) 时域分析结果：序列的平均值为序列的最小值为序列的最大值为序列的标准差为(2)频域分析结果:械动佶号的波形这是所拟合得到的振动信号的图像。O.G0.50-500400300 2K -1000100200。

7、300400500日 k rii 炉 w JjiArtR严打tip厘 ndotw h*lp一czwe-pwtUELM E 二匕：JJUHUAMJLLFure 2-LXfile Edl iew Insert Iools Dekiop Widow belp p已输I k叫轡區XTQ| 口图口?上图为FFT频谱图，从该频谱中可以看到有三个主要高峰值,即在 0Hz, 100Hz, 250Hz处* figure 33111111111050 1OQ 150200250300350 4W 45050QFrequency在功率谱中可以很明显的看到振动信号中有100Hz和250Hz两个主要的频率。表明信号中含有这两个频率的周期成分。五、总结通过使用MATLAB的相关函数编写的程序对这一所拟合的振动信号进行了时域分析和频域的分析，得到了关于该振动信号的时域分析结果并绘制出了频域分析图谱。通过这学期课程的学习，我认识到 MATLAB款非常强大的软 件，应用在各个行业，包括航空航天、人工智能、机械、控制、科学 运算等各个领域。目前我所学习和应用的只是它强大的功能非常小的 一部分，之后还要学习 MATLA更多的知识，还需更加努力。