测试技术基础 教学课件 作者 李孟源第5章.pptVIP

图5.29 随机变量x和y的相关性 对于能量型变量x(t)和y(t)之间的相关程度常用相关系数ρxy表示： （5.51） 　　通常，|ρxy|≤1。当|ρxy|=1时，说明x(t)和y(t)两变量是理想的线性相关。|ρxy| =0表示x(t)和y(t)两变量之间完全无关。 　在MATLAB中， 利用corrcoef工具箱函数可以从输入矩阵中计算出相关系数矩阵。 　　采样定理要求fs≥2fm，但采样频率fs并非选得越高越好。 由N个时域采样点进行离散傅里叶变换，得到N个频域点，通常称为N条谱线，对应的频率范围为［-fs /2,fs /2］，因此相邻谱线的频率增量为 (5.34) 　　可见当采样点数N一定时，采样频率fs越高，频率增量大，频率分辨力越低。因此，在满足采样定理的前提下，采样频率不应选得过高， 一般取fs=(2~3)fm就够了。 　　3. 采样点数N的选择 　　离散傅里叶变换为使用计算机进行频谱分析提供了理论依据， 但还存在一个实际问题， 就是计算工作量太大， 即使利用计算机这个强有力的快速计算工具也要花费很长的时间。 于是人们力图寻找一种快而简便的算法， 使离散傅里叶变换真正具有实用价值。 1956年美国人库利(J.W.Cooley)、 图基(J.W.Tukey)提出了一种称为快速傅里叶变换(FFT)的算法， 问题才得到了圆满解决。 FFT算法将DFT算法的计算速度提高到原来的N/log2N倍，使傅里叶变换可以在很短的时间内完成。 目前已有很多关于离散傅里叶变换的硬件、软件和专用机可供使用。 　　4. 窗函数、截断和泄漏 　　信号的历程一般较长，在进行数字信号处理技术时要进行截断。 截断就是将无限长的信号乘以有限宽的窗函数。 “窗” 的意思是指透过窗口能够 “看到‘外景(信号)’” 的一部分。 最简单的窗是矩形窗， 如图5.22所示， 其函数为 (5.35) 其频谱函数为 (5.36) 　　对信号截取一段(-T, T)， 就相当于在时域中x(t)乘以矩形窗函数ω(t)， 于是有 图5.22 矩形窗 图5.23 常用窗函数 (a) 汉宁窗; (b) 哈明窗; (c) 高斯窗; (d) 三角窗 　　比较这五种窗，矩形窗旁瓣最高但主瓣最窄， 高斯窗旁瓣最低但主瓣却最宽，最理想的窗函数应该是主瓣窗窄而旁瓣低。 因此在处理数据时，要根据具体要求来选择窗函数。 一般来说应注意下述几点： 　　（1）如果要分析信号中那些幅值很小的频率成分(即次要的频率成分)，则不能用矩形窗，应该用泄漏最小的高斯窗， 因为那些幅度较小的谱密度将被矩形窗本身引起的皱波所淹没。　　（2） 如果仅仅分析信号的主要频率成分，而不考查频谱的细微结构， 则可用计算最为简单的矩形窗。 　　(3) 如果要两者兼顾，则可用汉宁窗或哈明窗，而哈明窗的应用最为广泛。  　　需要指出的是， 除了矩形窗外， 其它窗在对时域函数截断的同时， 还对时域函数的幅值有影响， 导致频域函数幅值下降， 因而要乘以一个修正系数进行修正，这点在计算时要特别注意。 　 　　5. 平均化处理 　　离散傅里叶变换是连续傅里叶变换的一种近似。 对信号进行截断分析， 用数学的语言来说就是抽出总体信号的一个样本进行分析。 如果多抽出一些样本进行离散傅里叶变换， 最后取其平均值， 必然会抵消一些随机误差而获得较高的精度， 这种方法称为平均化。 平均化处理在数据处理中得到了广泛的应用， 具体做法是先把足够多的点数采入计算机存储器， 然后一段接一段地进行分析， 最后取平均。 若总点数不够， 取用时可在各段之间交叉， 使同一数据能够多段重复使用。 5.4.4 快速傅里叶变换(FFT) 　　快速傅里叶变换(FFT)是离散傅立叶变换（DFT）的快速算法， 它在确定DFT的系数时使所要求的乘法及加法次数减少。 FFT的算法有很多种，其中大多数已编制了程序， 从而使应用于数字频谱分析、 滤波器模拟及相关领域的计算技术产生了较大的发展。 　　在按照式（5.31）的离散傅里叶变换直接计算X(k)值时， 对于N个X(k)中的每一个必须作N次x(n)乘以e-j2πkn/N，所以共有N2次复数乘法运算，而且还要作N(N-1)次复数相加的运算。 计算工作量直接与N的大小有关，当N较大时，计算工作量将急剧增加。  　FFT算法的实质就是把一个长数据序列x(n)经多次分选抽取， 最终分割成n/2个，每个有两个数据的序列作DFT计算，分别算出分割后比较短的子序列的频谱， 然后按一定的规则组合，即可得到整个序列x(n)的频谱。 　　例如有一数据序列{x(n)}， n=0， 1， 2， …，N-1，如图5.24所示(N=8)。 将序列{x(n)}按偶数项和奇数项