matlab中dfft,测控技术大作业 徐尧 02011417.doc

测试技术大作业 徐尧

指导教师 贾明平

课程：机械工程测试与控制技术Ⅱ

时间：2014年6 月1 日

摘要：本次实验通过对具体给定的实际题目来进行分析处理以获得信号处理方面处理结果和结论。项目设计围绕动态信号采集、分析与处理的基本原理，与方法进行，包含以下三个方面的内容：

信号的仿真、采集和处理：这个实验模拟不同频率、幅值、相位和白噪音的情况下，不同采样频率、采样长度与分辨率等之间的关系，通过matlab仿真信号并且进行DFFT处理以获得时域和频域的图像，从而进一步得到需要的结论；

基于计算机的声音信号的采集与分析：本实验利用计算机上的声卡和AD来获得若干人说同样话情况下的信号，然后用matlab进行频谱分析，并作出时域、频域图，进而分析其中的频率特性和采样方法；

机械运行数据分析与处理：实验依靠已经采集好的处理为ASCII码的TXT文件数据进行处理分析，在用matlab导入数据后对数据进行频谱分析，做出对应不同转速下的时域与频域图，从中进行观察并且分析发现电机可能出现的问题，进而提出问题所在。

创新方面：这次项目设计核心是利用matlab进行数据处理与分析，为后续进行GUI进行了基础搭建。

关键词：matlab；时域-频域分析；机械故障检测分析；

目 录

1、 信号仿真、采集与分析处理1

1.1 题目描述1

1.2 分析思路1

2、 基于计算机的声信号采集与分析9

2.1 题目描述9

2.2 实验过程描述9

2.3 题后讨论10

3、 机械运行数据分析与处理12

3.1 题目描述12

3.2 数据处理过程12

3.3 题后讨论12

4、心得

1、 信号仿真、采集与分析处理

1.1 题目描述

号采集过程中一般需要考虑以下几个参数：信号频率、采样频率、采样长度等，不同参数的数值设定对于信号采集的效果会产生直接影响， 为了掌握信号采集过程中这些参数对

采集过程及其效果产生的影响，可以通过 Matlab 或 C 语言对信号采集与分析处理的过程进行仿真分析，具体要求如下：

利用 Matlab或 C 语言产生信号，

其中：频率需要考虑低、中、高，典型的如：f1=30Hz、f2=400Hz、f3=2000Hz。每位同学可按照此思路自己取值，与他人不同；

n(t) 为白噪声，均值为零，方差为 0.7；频率、幅值、相位任意设定，要求每人不同；

对上述等式进行 DFFT 处理。

讨论：

(1)通过设置不同的采样频率，画出时域波形和傅里叶变换后的频谱图，用数据验证采样定理。讨论在采样点数一定(2的整数次方)的情况下，如 1024 点，2048点，4096点，采样频率对信号时域复现、频域分析的影响；

(2)采样频率、采样长度(采样点数)与频率分辨率的关系；

(3)通过设置不同幅值的信号与噪声，讨论噪声对信号时域分析和频域分析的影响；

(4)考虑矩形窗和Hanning窗对频谱的影响。

1.2 分析思路

1.2.1 背景知识介绍

采样定理：在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中最高频率fmax的2倍时(fs.max>2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。(以上摘自百度百科)为了避免发生混叠现象，采样频率fs必须大于信号最高频率fc的两倍，即fs>2fc，这就是采样定理。实际中，一般采样频率应选为被测信号中最高频率的2.56倍以上。(摘自《测试技术》)

Hanning窗：汉宁窗又称升余弦窗，汉宁窗可以看作是3个矩形时间窗的频谱之和，或者说是 3个 sinc(t)型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了 π/T，从而使旁瓣互相抵消，消去高频干扰和漏能。可以看出，汉宁窗主瓣加宽并降低，旁瓣则显著减小，从减小泄漏观点出发，汉宁窗优于矩形窗．但汉宁窗主瓣加宽，相当于分析带宽加宽，频率分辨力下降。

DFFT：快速离散傅里叶变换。

高斯白噪音：所谓高斯噪声是指它的概率密度函数服从高斯分布(即正态分布)的一类噪声。如果一个噪声，它的幅度分布服从高斯分布，而它的功率谱密度又是均匀分布的，则称它为高斯白噪声。(摘自百度百科)

1.2.2 实验参数的拟定和过程

本次实验首先设定的参数如下：

其中n(t)是白噪音。

白噪音的生成程序如下：

程序中，参数FC代表方差。最后获得的nt直接加到信号上就可以了。

采样频率对信号时域复现、频域分析的影响

在这里我们先设采样长度为L=2048.。然后在不同的采样频率下进行频谱分析。