测控技术大作业 徐尧 02011417

测试技术大作业02011417徐尧指导教师贾明平课程：机械工程测试与控制技术Ⅱ时间：2014年6 月1 日摘要：本次实验通过对具体给定的实际题目来进行分析处理以获得信号处理方面处理结果和结论。

项目设计围绕动态信号采集、分析与处理的基本原理，与方法进行，包含以下三个方面的内容：（1）信号的仿真、采集和处理：这个实验模拟不同频率、幅值、相位和白噪音的情况下，不同采样频率、采样长度与分辨率等之间的关系，通过matlab仿真信号并且进行DFFT处理以获得时域和频域的图像，从而进一步得到需要的结论；（2）基于计算机的声音信号的采集与分析：本实验利用计算机上的声卡和AD来获得若干人说同样话情况下的信号，然后用matlab进行频谱分析，并作出时域、频域图，进而分析其中的频率特性和采样方法；（3）机械运行数据分析与处理：实验依靠已经采集好的处理为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 语言产生信号，)()2s i n ()2s i n ()2s i n ()(333222111t n t f a t f a t f a t x ++++++=φπφπφπ其中：频率需要考虑低、中、高，典型的如：f 1=30Hz 、f 2=400Hz 、f 3=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 实验参数的拟定和过程本次实验首先设定的参数如下：)()302000*2sin()20200*2sin(*2)1030*2sin()(t n t t t t x ++++++=πππ其中n （t ）是白噪音。

白噪音的生成程序如下：程序中，参数FC 代表方差。

最后获得的nt 直接加到信号上就可以了。

（1）采样频率对信号时域复现、频域分析的影响在这里我们先设采样长度为L=2048.。

然后在不同的采样频率下进行频谱分析。

（以下图片figure1是时域图，figure2是频域图） fs=1000hz 时，图像如下：fs=2000hz=fmax时，图像如下：fs=4000hz时，图像如下：fs=5120hz=2.56fmax时，图像如下fs=10khz时，图像如下：分析与结论：通过观察图像，我们可以发现，随着采样频率fs的升高，时域中图像越来越近似于一个正弦信号，信号细节越好，这与输入信号是有关系的。

且可以看出，在fs<2fmax的时候，时域图像里的信号是杂乱无章的，看不到信号周期，而且fs较小时在频域图中也看不到高频分量，从而导致错误。

fs=4000hz的时候是一个分水岭，此时频域图正好在末尾出现了高频分量，从而可以证明采样定理的正确性。

但同时，随着采样频率上升，频域中的主瓣也变大了，信号效果也变得不好了。

（2）采样频率、采样长度（采样点数）与频率分辨率的关系采样频率与分辨率的关系可以在上图中看出：采样频率越大，图像的分辨性能越好，频率分辨率越高。

采样长度与频率分辨率的关系可以用以下的图表示，其中fs=5120hz：L=2048时的图像：L=4096时的图像：L=1024时的图像：由以上三幅图对比可以看见，采样频率一定时，采样长度越大，频域中频率的主瓣越小，频率分辨率越大。

（3）通过设置不同幅值的信号与噪声，讨论噪声对信号时域分析和频域分析的影响对于输入信号，我们改变其幅值，改变后为如下（信号幅值为原来10倍，白噪音不变）：)()302000*2sin(*10)20200*2sin(*20)1030*2sin(*10)(t n t t t t x ++++++=πππ我们在L=1024，fs=5120hz 的情况下作图：接着，我们将幅值改为初始的0.1倍，白噪音幅值不变，同样的，在L=1024，fs=5120hz 的条件下作图：从这两幅图中可以看出，信号幅值对于图像来说是比较重要的，当幅值远大于白噪音时，可以清楚的看到干扰被抑制了，分辨率相对比较高；同理，当信号幅值接近白噪音甚至比白噪音小的时候，在时域就基本上看不出信号了，但频域里面只能够勉强分辨出一到两个频率，但误差都很大。

改变噪音的幅值得到的结果也是相似的。

所以可以说，信号幅值和白噪音的幅值的比较对信号分辨率是有影响的。

将幅值改为原来的值后，噪音的方差改变为1，在L=2048，fs=5120hz时所得图如下：对比，方差为0.7时的图：观察两图，对照可以发现主频仍然可以看到，但是如果把坐标系统一后，能够发现方差为1时，频域图里的干扰噪音的高度相对会变大一些，0.7的相对矮（小）一些。

进而可以得出结论，白噪音对图像也是有影响的，主要作用于非主瓣、旁瓣的频率上。

（4）考虑矩形窗和Hanning窗对频谱的影响a矩形窗（fs=5120hz）：长度为300时图像：长度为900时图像：从中可见，矩形窗函数对时域图形影响较小，但对于频谱影响比较大。

可以明显发现，矩形窗函数越长的时候，主瓣越窄，频谱图效果更加好。

B.对于Hanning窗，同样的方法，长度为300时的图像：长度为900时的图像：通过观察可以发现，Hanning窗函数采样长度越长，主瓣越窄，分辨力越高。

将矩形窗函数与Hanning窗函数对比，能够发现Hanning窗的旁瓣小得多，泄漏也少得多，但是Hanning窗主瓣较宽。

2、基于计算机的声音信号的采集与分析2.1 题目描述现代计算机具有对声音、视频进行采样的功能，把模拟信号转换为数字信号。

通过计算机上的麦克风及声卡与AD，录制3人以上在不同环境噪声、不同发声状态下讲同一句话，如“机械工程测试与控制技术”语句。

先利用软件将录制语音转换为数据文件ASCII码（text文本），再利用1.1的软件进行频谱分析，画出时域、频域图形。

录音软件：如：蓝光影音Mp3录音机Matlab读取MP3文件：mp3readwrite/matlabcentral/fileexchange/13852-mp3read-and-mp3write?cont roller=file_infos&download=true讨论：（1）该设置至少为多少的采样频率？采样长度多长为合适？（2）不同人员讲话声音的时域、频域有什么区别？根据你的分析，该怎样区分不同人员的讲话声音？（3）要使他人不易识别你的讲话声音，该怎么处理？2.2 实验过程描述（1）数据采集：我采集了同宿舍三个人（均为男生）说的同一句话“测试技术”，然后转化为TXT文本，接着用matlab导入后进行频谱分析处理，最后做出时域和频域图。

（2）录音文件和转化后的TXT文件在附件里面。

（3）做出时域与频域图。

同学1同学2同学32.3 题目讨论2.3.1 该设置至少为多少的采样频率？采样长度多长为合适？查阅过资料后得知，人类的发声频率一般为65Hz 到1100Hz 之间，因而为了保证采样定理，我们取的采样频率为Hz Hz 281656.21100=⨯。

而采样长度应该与声音信号的长短来决定。

在采样频率一定的条件下采样点数（采样长度）越大，数据量越大，但同时采样较好。

在此处，我们不同信号采样长度不同，这也是因为录音时间不同，样本不同而产生的。

在这里大致约为4w 到10w 点.2.3.2 不同人员讲话声音的时域、频域有什么区别？根据你的分析，该怎样区分不同人员的讲话声音？在时域图中，我们虽然看到了每个人对应每个字的声音幅值不同，但是经过听取录音可以发现，这个是与个人发声习惯以及声音大小有关的，所以不可以用作区分的依据。

在频域图中，我们能发现每个人都有自己的发声频率，因而可以通过降噪滤波处理后利用各个峰值进行人员的区分。