当前位置:
文档之家› 连续时间信号频谱分析研究及MATLAB实现
连续时间信号频谱分析研究及MATLAB实现
周期 ,即频 谱数值 计算 的范 围 ;而 在某 时间段 上对信
即得到离散频率点上 的近似计算式 :
N
. 一 .
-
号进 行截取 的方 式 .即不 同窗函数 的应用 .决定 了信 号频谱估计 的精度和有效范 围。 设 要 分 析 连 续 时 间 非 周 期 信 号 厂 f在 频 率 范 围 ( )
f()x (j n f) n ep一 0 ) 3 3
n= 0
对 ID T : F 计算式 显然有 I c
Ff J . T D T ( ] T Fk 1l () 2 ・ F[ n = ・ ( 4 f ) )
一
l
该式表明 ,利用D TF T计算连续 时间傅 ̄ n 变换 F (F ) - I - 的频谱时 ,除 了计算时域样点的离散傅里叶变化的频谱 , 七 ,还要gF k乘 以取 样时间间隔 ,才能得出结果 。 () q () 直 接 计 算 DF 的 复 杂 度 为 序 列 长 度 的 平 方 ,即 T
取信 号 的时 间段 、如何 选择 时域采 样率 ,以及在 时 间
Hz
,
采 样 时 间间 隔 ( 间分 辨率 ) =l , 5 ,而 时 为T / = ms f
段 上对信 号进 行截取 的方式 。截取 信号 的 时间段 长度 根 据要 求的频 率分辨 率可 以得 出信 号 时域截 断长 度为
2 、频 谱 分 析 的若 干 问题 讨 论
( )根据频率分辨率要 求确定分析信号厂( 的截取 2 f ) 时间长度 。要使所分析的频率分辨率达到△ ,即每隔 厂 计算一个频率点 ,那么对信号 的截取时间长度三 必须满足 L / ≥1A f,根据截取时间长度 和采样时间间隔 可以计
算出截取 时间信号离散化之后的序列点数Ⅳ ,也可 以由计
算采样率 和频率间隔 来等价计算出序列点数,即:
N= / + = /厂+ ( { l1 I I 5 )
( )根 据信 号 时域 波形 特 性来 应 用不 同 的窗 函 3 数 。使 用 窗 函数 可 以控制频 谱主 瓣宽度 、旁瓣 抑制 度
等参 数 ,从 而更 好地进 行波 形频 谱分析 和滤波 器参 数
3 、分 析 实例
试对 f: i 2 t+ i 2 8t ( s (t0 s ( 0 进行频谱分析 .要求 ) n7 ) n ̄ ) 5
~1 H ,频率分辨率 为1 z 0 。根据 H 频域 离散 化将 对应 于 时域 信号 的周期化 。因此 ,对于 分析 的频率范 围为0 0 z 20 连 续非 周期信 号频 谱进 行数值 计算 时 ,要确 定如 何截 所要 求的分析频率范围可 以确定信 号的采样率 为 = 0
◆ 建设维 护◇
个计算频率点 ,这Ⅳ个离散频率点 以角频率表示为 :
’
决定 了时域周 期化 的周期 ,对应 于频率 抽样 的频率 间
隔 ,即频 率分辨 率 ;时域采 样率决 定 了频域 周期化 的
0 k =k 3= 0 3
、 பைடு நூலகம்N
, = 0 1L k ,, N 一1
…
l
时间分 辨率 ) 满足T= / 12 1 L≤ / 。
在 《 计算 数学》杂志上发表了著名的 《 机器计 算傅里 叶 级数 的一种算 法》论 文后 ,桑德一 图基等快速算 法相继 出现 ,又经人们进 行改进 .很快形成 了一套高效计算 方 法 ,这就是现在 的快速傅里叶算法FT F 。其思想是将长序 列分解为若干短序列进行D T F 计算 ,然后通过若干旋转 因 子的复数乘法和复数加法合成最终 的结果 。FT F 算法 中, 如果序列长 度是2 的幂次 ,可将序列长 为N F 分割 为 的D T 两个长为N 2 / 的子序列 的D T F ,成 为基2 FT - F 。快算傅里 叶变化算法 只需要O Nlg 的计算复 杂度 ,FT 法为 ( o Ⅳ) F算 数字信号处理技术应用于各种信号的实时处理创造 了条 件 .大大推动了数字信号处理技术 的发展[】 1 。
D J ) 5 ,库利 ( . . o l ) 图基 (. .uy ( 2。16 年 v 9 TW C o y和 e JW Tk )
信号时域波形特性选 择使 用不同的窗函数 。 ( )根 据分析的信号频率 范围确定采样率 。要 分 1 析信号 在频率范 围[ ] 0 内的频谱 ,则采样率 必须满 足 , 采样定理 ,lf≥ L,相应地 ,采样 时间间隔 ( 为  ̄ .2 l, J 也称
对该 频谱 进 行频域 采样 。根据 时域 频域 的对 应关 系 ,
的设计 。将 窗 函数与信 号 的时域波 形或频 谱进 行相 乘
的过程 ,就称 为对信号做 时域 加 窗和频域 加窗 。不 同 窗函数 与信 号时域 波形相 乘就是 以不 同 窗函数对 时 间
无限长 的连续信号厂 ( 进行 时间段截取 。 )
通过 以上 原理 可知 ,连续 非周期 信 号频谱 的数值
计 算 必须首 先对 信号 时域 采样 ,得 到时 间离散 化 的信 号 ,时域采 样必 须满 足或近 似 满足采 样定 理 。根 据 时 域 频域 的对应 关 系 ,时域采样 将 导致所 得 的抽样 信号 频 谱周期 化 。然 而 ,为 了使 周期 化后 的抽样 信号 频谱 便 于计算 机 处理 .还必 须再 将其 频域离 散化 ,方 法是
[ 】 0 内的频 谱 ,且 要求分 析的频谱分辨率 ( , 数值计算 的 频率间隔 ) 为 .则首先应 根据信号频率 范围确定采样 率 ,在根 据所 要求 的频率分 辨率确 定截 取 时间长度 .
从而计 算出所需计算F T F 的序列长度 ( 点数 ) ,最后根 据
’
1
F( ) e = o
田