连续时间信号频谱分析研究与MATLAB实现.pdfVIP

个计算频率点，这Jv个离散频率点以角频率表示为 决定了时域周期化的周期 ．对应于频率抽样的频率间 0 ∞=ko。=k；，k=0LJv一 (2) 隔 ，即频率分辨率 ：时域采样率决定了频域周期化的 周期 ，即频谱数值计算的范围 而在某时间段上对信 即得到离散频率点上的近似计算式： 号进行截取的方式 ．即不同窗函数的应用 决定了信 Iv一1 f(∞)=丁∑f(n)exp(-jon) (3) 号频谱估计的精度和有效范围。 n=O 设要分析连续时间非周期信号f(t)在频率范 围 对bLDFT计算式，显然有 [0，．]内的频谱，且要求分析的频谱分辨率f数值计算的 F(Q T·DFT[八(”)]=T·F(k) f41 频率间隔)为△／，则首先应根据信号频率范围确定采样 Ⅳr 率 ，在根据所要求的频率分辨率确定截取时间长度， 该式表明，利用DFT(FFT)计算连续时间傅里叶变换 从而计算出所需计算FFTfl~,序列长度 (点数) 最后根据 的频谱时，除了计算时域样点的离散傅里叶变化的频谱 信号时域波形特性选择使用不同的窗函数。 ( 还要将，(七)乘以取样时间间隔 才能得出结果。 (1)根据分析的信号频率范围确定采样率。要分 直接计算DFT的复杂度为序列长度的平方 即 析信号在频率范围[0，．，]内的频谱 ，则采样率必须满足 O(N)。1965年，库利 (TW ，Cooley)和图基 (JW ．Tuky) 采样定理，即，： ， 相应地 ，采样时间间隔71(也称为 在 《计算数学》杂志上发表了著名的 《机器计算傅里叶 时间分辨率)满足T一 ／，：1／2．／i 级数的一种算法》论文后 桑德一图基等快速算法相继 (2)根据频率分辨率要求确定分析信号／《，)的截取 出现，又经人们进行改进 ．很快形成了一套高效计算方 时间长度。要使所分析的频率分辨率达到 ，即每隔 法，这就是现在的快速傅里叶算法FFT。其思想是将长序 计算一个频率点，那么对信号的截取时间长度三必须满足 列分解为若干短序列进行DFT计算．然后通过若干旋转因 L≥i／△厂．根据截取时间长度 和采样时间间隔 就可以计 子的复数乘法和复数加法合成最终的结果。FFT算法中， 算出截取时间信号离散化之后的序列点数Ⅳ，也可以由计 如果序列长度是2的幂次，可将序列长为N的DFT分割为 算采样率 和频率间隔△，来等价计算出序列点数 即： 两个长为N／2的子序列的DFT，成为基2一FFT。快算傅里 叶变化算法只需要O(NlogⅣ)的计算复杂度，FFT算法为 N=J ／TJ+1=l ／△厂f+】， L- 一 L J _J ， 数字信号处理技术应用于各种信号的