机械故障诊断--研讨题1探究.docVIP

研讨1-信号采集、谱分析、窗函数、整周期采样、谱校正 信号采集过程中一般需要考虑以下几个参数：信号频率、采样频率、采样长度等，不同 参数的选择对于信号采集的效果会产生直接影响，为了掌握信号采集过程中这些参数对采集 过程及其效果产生的影响，可以通过Matlab 或C 语言对信号采集与分析处理的过程进行仿 真分析，具体要求如下： 利用Matlab 或C 语言产生信号： 其中：频率需要考虑低、中、高，典型的如：f1=30Hz、 f2=400Hz、f3=2000Hz。每位同 学可按照此思路自己取值，与他人不同； n(t) 为白噪声，均值为零，方差为0.7； 频率、幅值、相位任意设定。要求每人不相同； 对上述等式进行DFFT 处理。 讨论： 1）通过设置不同的采样频率，画出时域波形和傅里叶变换后的幅频谱图，用数据分析验证采样定理。讨论在采样点数一定（2 的整数次方）的情况下，如1024 点、2048 点、4096点，采样频率对信号时域复现、频域分析的影响； 2）采样频率、采样长度（采样点数）与频率分辨率的关系； 3）通过设置不同幅值的信号与噪声，讨论噪声对信号时域分析和频域分析的影响； 4）考虑矩形窗和汉宁窗对频谱的影响； 5）整周期采样及谱校正。  分别选择采样频率200、1000、4300、7000Hz，选择x=4*sin(2*pi*40*t+10)+7*sin(2*pi*500*t+50)+9*sin(2*pi*2100*t+70)+white_noise;(1~4问有小问题，因为上式x中相位应该除以180乘以pi变为角度，第五问纠正过来，而且第五问应该分析三个峰值，我只分析了一个，其余两个，可以谱校正取最大值（即峰值）时分别选取N/2，N/4,N/20从而取到第二第三峰值（即另外两个频率最大值）。) 1) MATLAB程序结果框图: 通过四个图的比较分析，当采样频率2时，即，波形无法显示在高频率段的状态，由图一、二可以看出波形出现混叠，而由图四、五可知，当，即时，波形在、、三个频率时的幅值高出其他频率，因此对于一个包含多个频率成分的信号进行处理时，要求采样频率必须高于信号频率成分中最高频率的两倍，即 这就是采样定理。 2） 采样点数一定的情况下，高频谱比低频的峰值频率带更窄，更可以清晰的看出峰值频率，并且两个峰值频率的混叠情况比低频的要少。 采样频率，采样点数N，采样间隔： 频率采样间隔 越小，分辨力越高。 三、设置=1，=1，=2 ，n(t)均值为零，方差分别为0.7， 1.5， 3.2 由图可知，白噪声1的主瓣幅值最高，白噪声3幅值最低，说明白噪声的方差越大，时域信号的相应越准确，噪声对信号的频域分析影响越小。 四、加窗对频谱的影响： 五、整周期采样，谱校正 只有信号的截断长度Tf0，获得准确的频谱。此即为整周期采样，而此信号包含多个频率成分，因此取、、的最小公倍数为T=21000Hz Matlab输出结果： MATLAB程序： 第一问： figure(1) %采样频率为200Hz fs=200; N=1024;%采样点数 v=0:N-1; t=v/fs;%时间序列 white_noise=0+sqrt(0.7)*randn(1,length(t)); x=4*sin(2*pi*40*t+10)+7*sin(2*pi*500*t+50)+9*sin(2*pi*2100*t+70)+white_noise; subplot(211); plot(t,x); title(采样频率为200Hz); xlabel(t); ylabel(y); number=1024; y=fft(x,number)/number; n=0:length(y)-1; f=fs*n/length(y); subplot(212); plot(f,abs(y)); xlabel(频率(Hz)); ylabel(幅值); set(gca,xlim,[0,100]); figure(2) %采样频率为1000Hz fs=1000; N=1024;%采样点数 v=0:N-1; t=v/fs;%时间序列 white_noise=0+sqrt(0.7)*randn(1,length(t)); x=4*sin(2*pi*40*t+10)+7*sin(2*pi*500*t+50)+9*sin(2*pi*2100*t+70)+white_noise; subplot(211); plot(t,x); title(采样频率为1000Hz); xlabel(t); ylabel(y); number=1024; y=fft(x,number)/number; n=0:length(y)-1; f=fs*n/length(y); subplot(21