动态测试信号采集仿真与实例分析..docxVIP

动态测试信号采集仿真与实例分析课程：机械工程测试与控制技术Ⅱ 石卓指导老师:贾民平 2014 年 6 月 8 日摘要:测试技术项目设计主要包括三个部分：信号仿真、采集与分析处理，基于计算机的声信号采集与分析，机械运行数据分析与处理。通过对测试信号采集仿真，并运用Matlab数据编程处理软件，进行信号的DFFT处理，做信号的时域图、频谱图，以及研究不同的窗函数对信号的影响,并能够分析信号的频谱及其相应的特征。从而验证采样定理的正确性，以及采样定理在实际生活中的应用，针对不同的人如何分辨其声音的特征，以及掌握机械运行时通过频谱图进行故障分析等。关键词：Matlab；频谱分析；信号采样；故障诊断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）考虑矩形窗和汉宁窗对频谱的影响。1.2程序1、未加窗函数的程序： fs=5000; number=2000;N=2048;T=1/fs;t=(1:number)*T;n=(rand(1,number)-0.5)*sqrt(12*700);a1=pi/3;a2=pi/3;a3=pi/3;x=2*sin(2*pi*30*t+a1)+3*sin(2*pi*400*t+a2)+4*sin(2*pi*2000*t+a3)+n;subplot(2,1,1);plot(fs*(1:1000),x(1:1000));title(信号时域波形);xlabel(t)ylabel(x(t))X=fft(x,N)/number;f=fs/2*linspace(0,1,N/2+1);subplot(2,1,2);plot(f,2*abs(X(1:N/2+1)));title(信号的频谱图);xlabel(f(Hz))ylabel(幅值)2、加上矩形窗： fs=5000;number=2000;T=1/fs;N=2048;t=(1:number)*T;n=(rand(1,number)-0.5)*sqrt(12*0.7);a1=pi/3;a2=pi/3;a3=pi/3;x=2*sin(2*pi*30*t+a1)+3*sin(2*pi*400*t+a2)+4*sin(2*pi*2000*t+a3)+n;x=x/max(abs(x));wn=triang(2000);x1=x.*wn ;X=fft(x1,N);f=fs*(1:N);plot(f,abs(X)/max(abs(X)));xlabel(Hz);ylabel(幅值);title(矩形窗);3、加上汉宁窗后： fs=5000;number=2000;T=1/fs;N=2048;t=(1:number)*T;n=(rand(1,number)-0.5)*sqrt(12*0.7);a1=pi/3;a2=pi/3;a3=pi/3;x=2*sin(2*pi*30*t+a1)+3*sin(2*pi*400*t+a2)+4*sin(2*pi*2000*t+a3)+n;x=x/max(abs(x)); wn=hanning(2000);x1=x.*wn ;X2=fft(x1,N);f=fs*(1:N);subplot(2,1,2);plot(f,abs(X2)/max(abs(X2)));xlabel(Hz);ylabel(幅值);title(汉宁窗);1.3 实验数据1、在N=2048点，fs=4000Hz、5000Hz、6000Hz时：图一 fs=4000Hz 图二 fs=5000Hz图三 fs=6000Hz2、设定采样频率fs为5000Hz，分别取采样点数N为1024点,4096点时：图四 N=1024图五 N=40963、设定采样频率fs为5000Hz，采样