《数字信号处理》课程设计项目3 振动信号处理.pptxVIP

《数字信号处理实验与课程设计教程--面向工程教育》 课程设计项目3 振动信号处理9.1 振动信号的时域及频域分析9.2 基于振动信号分析的电机轴承故障检测9.3 数字滤波器在地震信号分析中的应用 9.1 振动信号的时域及频域分析一、课题背景与研究意义 振动是指物体或结构随时间变化相对其平衡位置所做的往复运动，通常用位移、速度和加速度来描述，它是一种普遍存在的自然现象，如:地震、汽车驶过桥梁引起的振动、机器运转引起的振动等，对周围环境或机器结构造成了不利的影响。振动信号处理是通过数值计算方法对振动测试所得的信号进行分析、加工和识别，以提取其中的有用信息。二、课程设计目标要求 本课题研究的内容是：基于Matlab软件对振动信号进行基本的时域及频域分析，如在时域中对地震信号进行预处理，消除其中的趋势项及噪声干扰、振动位移信号测量；在频域中对地震信号进行FFT谱分析、计算振动测试信号的三分之一倍频程谱等，对实验结果进行分析并撰写4000~5000字的课题论文。课题论文包含以下内容：1)题目、摘要、关键词、引言2)内容：振动信号的时域及频域分析3)包括：理论、程序(及注释)、图形、结果分析4)结论、参考文献 三、课程设计内容与参考 设计MATLAB程序模块，实现以下功能：(一)振动信号的时域处理1.地震信号的预处理(1)基于最小二乘法的地震波趋势项消除设计要求：地震响应信号如图9.1 (上方图)所示，由于周围环境的干扰，其中混入了不需要的多项式趋势项(为4次函数:y=a0-a1k- a2k2- a3k3- a4k4，k=1~n)，该趋势项会直接影响信号的正确性，编程实现：采用多项式最小二乘法将其去除(下方图)，列出程序清单并显示结果波形。图9.1 消除地震信号中的趋势项 参考程序：见\课程设计项目参考源码\kcsj3\课程设计题3.1\kcsj311.m (在空格中填入注释) (2)基于五点滑动平均法的地震信号消噪设计要求：地震响应信号如图9.2(上方图)所示，其中含有随机干扰信号，该波形呈现出许多毛刺，为了削弱干扰信号的影响，对其采用五点滑动平均法进行消噪(下方图)，公式如下：其中x---采样数据；y---平滑后的结果；m---数据点数。 图9.2 五点加权平均法对地震响应信号进行平滑参考程序：见\课程设计项目参考源码\kcsj3\课程设计题3.1\kcsj312.m (在空格中填入注释) *思考:对原信号x及平滑后的信号y进行快速傅里叶变换(N=256)，观察幅频谱波形，由此得出什么结论？ 2.振动位移信号测量 在振动信号测试中，由于仪器设备或测试环境的限制，有的物理量需要通过对采集到的物理量进行转换处理才能得到，如：在振动台上进行的高层楼房模型试验中，要测试模型上的测点相对于振动台台面上的振动位移是非常困难的，可采用对速度信号进行积分得到位移信号。设计要求：基于matlab编程，利用辛普森数值积分公式对速度信号x(k)求积分,得振动位移信号y(k)(图9. 3):(9.4)其中：x(k)是x(t)=sin(2*t)+sin(0.6*t)+sin(3*t)经fs=300Hz采样所得的离散速度信号;△t—采样时间步长。 图9.3 振动位移信号测量参考程序：见\课程设计项目参考源码\kcsj3\课程设计题3.1\kcsj313.m (在空格中填入注释) (二)振动信号的频域处理1.基于FFT算法分析地震波中的频率成分设计要求：地震波是由各种地球表面运动与地壳运动等复杂运动混杂而产生的波，其中含有多种多样的频率成分，编程实现基于快速傅里叶变换(FFT)对从地震波信号截取的P波信号、S波信号和面波信号进行频谱分析，提取其中的频率成分，并绘出各自的原始信号及幅频谱。本课题地震数据来自从首都圈地震记录文件中截取的地震波信号文件：hns.dat、hns1.dat(P波数据）、hns2.dat(S波数据）及hns3.dat(面波数据)。（采样频率Fs=50Hz）参考程序：见\课程设计项目参考源码\kcsj3\课程设计题3.1\kcsj314.m (在空格中填入注释) 数据分析：地震波经过傅里叶变换之后，不同地震震相的频谱分布特征是不一样的，问：(1)P波的频谱分布高频成分居多,最高到（ ）Hz。(2)S波成分次之，频率范围是（ ）~（ ）HZ。(3)面波主要频率成分是低频，频率范围是（ ）~（ ） HZ。 图9.4 地震信号种各波形的频谱分析 2.振动信号的三分之一倍频程谱 三分之一倍频程谱是一种频域分析方法，具有谱线少频带宽的特点，常用于声学、人体振动、机械振动等测试分析中。由于人耳听觉频率范围为20~20000Hz，范围非常广，如果采用等宽频程，要表现出低频到高频的声音信息