信号分析与处理及其应用-虚拟测试系统.docVIP

PAGE PAGE 1 旧教材第九章 信号分析与处理及其应用 §9-1 概述 信号分析通常是指分析信号的类别、构成以及特征参数，主要涉及信号的表示和性质。本书的第三章（信号及其描述）即属于信号分析的内容。信号分析中的一个最基本的方法是：把频率作为信号的自变量，在频域里进行信号的频谱分析。信号的处理是指对信号进行某种加工变换或运算，如滤波、变换、调制/解调、增强、压缩、估计、识别等加工处理。广义的信号处理也把信号分析包括在内。 信号分析方法有两大类：模拟量分析法（计算机时代以前用的多） 数字量分析法（计算机时代以来广泛采用） 例：模拟法相关。 延时 延时 x(t) y(t) 乘法器 积分器 用R.L.C 20世纪信号分析领域最重大的进展：FFT(1965年)；小波分析。 FFT是1965年由美国人库利和图基提出的DFT的巧妙的快速算法。 §9-2 离散傅里叶变换（DFT） 卷积定理 卷积： 可证明： 所谓卷积定理，指的是： 设 IFTFT IFT FT IFT FT 则有 IFT IFT FT IFT IFT FT 离散傅里叶变换的推导 推导图解见P225图9-2. 参看周期函数的复指数形式的傅里叶级数展开， 有（DFT）： IDFT： Ts为时域采样间隔；f0为频率分辨力，f0=1/ T0 ；T0为时域截断长度，T0=NTs 。信号被截断后再延拓为周期函数，周期为T0。 x(k)表示在离散时间点t=kTs处的采样值，即x(k)=x(kTs), k=0,1,2,…(N-1)。k为时域离散值的序列号。 X(n)表示在离散频率点f=nf0处的采样值，即X(n)=X(nf0), n=0,1,2,…(N-1)。n为频域离散值的序列号。 DFT对（式（9-5）和（9-6））可写成如下形式： 式中 式（9-10）的含义是：N个时域点与N个频域点相对应，每一频域点都要有N个时域点来求得。 DFT对（式（9-10）和（9-11））也可写成矩阵表示式。 离散傅里叶变换的几个重要问题 1.采样定理：对于带宽有限(频谱变化范围为0~fm)的连续信号进行采样，只有当采样频率时，才不会在频率域产生混叠现象，才可在采样处理后仍有可能恢复其原信号波形。 2．信号经时域截断相当于乘以“矩形窗”，在频域上产生“泄漏”效应。采用适当的“窗函数”代替简单的“矩形窗”可使“泄漏”减少。 3．时域的采样间隔决定了频域分析的频带宽度，即 时域截断长度决定了频域分析的谱线间隔——频率分辨力，即 4．为了用FFT法计算DFT,一般N取2的整数次方。 5．保持其它参数不变，增加采样点数N可提高频率分辨力。 6.由于采样、截断等固有特点，一般信号经DFT得到的频谱只能是“估计值”。存在着“混叠”、“泄漏”、“栅栏”等效应引起的误差。对随机信号而言，还存在着用“样本”去估计“总体”而必然存在的统计误差。 *本章的9-3，9-4，9-5节对检测技术、测试技术、信号分析与处理等专业人员来说很重要，是用计算机作数字信号处理和编程的基础，请对本学科感兴趣的同学参照第三章（信号及其描述）的知识和“信号分析与处理”或“数字信号分析与处理”类参考书进一步自学。 §9-3 快速傅里叶变换（FFT） DFT方法计算量太大，限制了应用。直到1965年，美国的Cooly和Turkey提出了一种快速计算DFT的算法。例如：由N个计算出N个,当N=1024时，DFT的复数乘法次数约为105万次，Cooly和Turkey的复数乘法次数5120次，仅为DFT的1/200。人们称这种快速算法为快速傅里叶变换（FFT）。算法中，规定N取2的整数次幂，因此也称基2型FFT。 目前实现FFT主要有软件和硬件两种方法。FFT是自(功率)谱、互(功率)谱、频率响应函数、相干函数等经典频域分析和许多相关分析方法的基础。 FFT的参数选择: 对一个连续信号作FFT，一般按以下步骤选取参数： （1）估计的截止频率或按所需的最高频率对作低通滤波。 （2）估计所需的频率分辨率。由FFT得到的是离散频谱， 相邻两谱线间的频率间隔（即上节的频率分辨率）必须小于，才能分辨出中相邻的两频率峰值。 （3）由采样定理（或）确定采样频率或采样间隔。 （4）由确定的一个样本的最小采样点数（，最小采样长度）。计算式中，为采样周期，为频率分辨率。 （5）对基2型FFT（关于基2型FFT，请看参考书），按2的整数次幂(=N)圆整采样点数N，使N=。 （6）选取适当的窗函数。 为了减少或抑制泄漏，需要用各种不同的窗函数（如矩形窗、汉宁（Hanning）窗、汉明(Hamming)窗、三角窗等）对时域信号进行加权处理。 总结：1.若难以估计的截止频