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

连续时间信号频谱分析研究及MATLAB实现 723001 【摘 要】 【关键词】 【中图分类号】 【文献标识码】 【文章编号】 （一）引言 机不能直接加以处理。为了实现数值计算，还需要对F (Ω) 进 在信号处理过程中，频域分析方法往往比时域分析方法 行离散化处理。将频率段 0 ~ 1/ T Hz 划分为 N 个计算频率 更方便和有效。对于确知连续时间信号，其频域分析可以通 点，这N 个离散频率点以角频率表示为： 过连续时间傅里叶变换来进行，但是，这样计算出来的结果 2π k k 仍然是连续函数，计算机不能直接加以处理。为了实现数值 ω ω0 N ,k 0,1, LN =−1 （2） 计算，还需要对其进行离散化处理，即采用离散傅里叶变换 即得到离散频率点上的近似计算式： （DFT）进行分析。DFT 的快速算法的出现，使 DFT 在数字通 N −1 信、图像处理、功率谱估计、系统分析与仿真、雷达信号处 F (ω) T ∑f (n)exp(=−jωn) （3） 理、光学、医学等各个领域都得到广泛应用。本文以正弦信 n 0 号为例，介绍用DFT的快速算法即快速傅里叶变化（FFT）实 对比DFT计算式，显然有 现确知连续时间信号的频谱分析，给出了MATLAB语言工具下 F (Ω) 2π ≈T •DFT[ f (n)] T =•F(k) （4） 的分析程序。 Ω k NT （二）信号频谱分析的原理 该式表明，利用 DFT(FFT)计算连续时间傅里叶变换的频 对于时间连续信号 f (t) ，其频谱分析可以通过连续时间 谱时，除了计算时域样点的离散傅里叶变化的频谱F (k ) ，还 傅里叶变换（CTFT）来进行。连续时间傅里叶变化特别适合 要将F (k ) 乘以取样时间间隔T ，才能得出结果。 于对时间连续信号的理论分析，但是，由于函数f (t) 和其频 直接计算DFT的复杂度为序列长度的平方，即O(N2 ) 。1965 年，库利（T.W.Cooley）和图基（J.W.Tuky）在《计算数学》 谱函数都是连续函数，不能够直接用计算机来处理，因此在 杂志上发表了著名的《机器计算傅里叶级数的一种算法》论 进行数值计算时必须将其离散化，然后利用离散傅里叶变换 文后，桑德—图基等快速算法相继出现，又经人们进行改进， （DFT）实现近似计算。 很快形成了一套高效计算方法，这就是现在的快速傅里叶算 设对连续时间信号f (t) 的截取时间段长度为