精密测控与系统课程第一次作业.docVIP

采样窗长对于频域信息的影响 摘要：本文简单的讨论了DFT时采样窗长对于频域信息的影响，分别从整周期截取，泄露和频率分辨率的角度讨论了采样时窗长的选取原则。t)。那么，窗函数的选择对于DFT的效果会有很大的影响，为了避免栅栏效应引起的频谱的畸变，对周期信号作整周期截取是获得正确频谱的先决条件。 在满足整周期截取的条件下，窗长越长，则由于加窗所带来的泄露现象就越不明显，图是一个幅值为频率为Hz的的DFT结果。其中窗长分别为和可以看到为时，泄露更少，效果更好。图4周期截取的幅频特性 另外，当信号中有多个频率时，有关。DFT的变换对应的模拟频率分辩率为：，所以采样长度N 的选择为： 例如当一个=100Hz =250Hz =300Hz 采样频率=1000Hz 采用N=10，20， 图3 N=10时的幅频特性 图4 N=20使得幅频特性 图5 N=100时的幅频特性 （1）最小频率间隔=50Hz,=1000Hz, 为能区分250Hz、300Hz的信号频率，N20；由图（3）可以看出，当N=10时，只有两个峰值，由于频率分辩率大于50Hz,f2 和f3产生了混迭。 （2）N=20时刚好能够满足频率分辨的条件，在250Hz和300Hz位置出现了两个峰值。第一峰值的左侧，第一和第二峰值之间的虚假谱峰，主要由于旁瓣泄露引起。 （） 在选取窗函数的窗长时，应该要满足以下三个原则，得到的频域信息才比较准确。 截取的窗长要是信号周期的整数倍，保证整周期截取，避免栅栏效应； 在允许的范围内，应尽量去较大的窗长，减少加窗带来的泄露； 在处理多频率信号时，窗长的选择还要保证频率分辨率足够小，可以分辨出信号中不同频率的信息，即满足。 附录： clc; clear all; close all; %% 产生一个正弦信号 A = 1;%正弦信号幅值 f0=100;%正弦信号频率 f0_2=120;% fs=1000;%采样频率 Ts=1/fs;%采样时间间隔 phi =pi/4;%初始相位 w_s=10;%采样周期数初始值 while w_s0 % w_s=10; % 采样周期数 w_length=fix(w_s*(1/f0)*fs);%采样窗口长度（数据点数） t=(0:w_length)*Ts;%采样时间点序列 x_s= A*sin(2*pi*f0*t + phi)+A*sin(2*pi*2.5*f0*t + phi)A*sin(2*pi*3*f0*t + phi);%窗口中的采样样本 %采样结果绘图 figure(1); plot(fs*t,x_s); title(原始信号采样结果); xlabel(time/milliseconds); %% % %计算与显示信号的FFT，DFT窗长与窗口长度相同 X_DFT=fft(x_s,w_length)/w_length; half_N=fix(w_length/2);%只显示一半 f=fs/w_length*(0:half_N); figure(2); stem(f,2*abs(X_DFT(1:half_N+1))); title(DFT amplitude results); xlabel(Frequency [Hertz]); %% %% 计算与显示信号的FT,1024点固定长度DFT DFT_n=1024; X_DFT=fft(x_s,DFT_n)/DFT_n; half_N=DFT_n/2;%只显示一半 f=fs/DFT_n*(0:half_N); figure(3); stem(f,2*abs(X_DFT(1:half_N+1))); title(1024 DFT amplitude results); xlabel(Frequency [Hertz]); w_s = input(输入采样周期数(中值循环输入0) w_s= ); end