CH10 DFT应用—信号谱分析.ppt

* 例(续)：语音信号的窄带谱图 采样率16kHz，R=16点=1mS，L=720点=45mS，Hamming窗 * 本章作业 10.2，10.8，10.15 CH10 DFT应用—信号谱分析 * 10.1 频谱仪工作原理 频谱仪用DFT完成模拟信号频谱分析过程 抗混叠低通滤波：非理想低通引入误差 信号采样离散化：量化和混迭引入误差 序列加窗有限长：时域加窗频域卷积畸变 DFT得离散频谱：频谱采样引起检测失误 * 信号频谱分析过程 信号频谱非带限 非理想低通 抗混叠引入误差 离散量化引入混叠 有限长引入窗频谱卷积 卷积“模糊化”频谱离散偏差 * 频谱分析注重点 要用V[k]作为Sc(jΩ)的有效估计，重点关心 单频信号中心频点的检测误差 可分辨相邻多频信号的最小频率间隔 单频信号频谱幅度的检测误差 多频信号频谱幅度的检测误差 由于V[k]获取过程中引入了多种偏差，必须了解如何减少误差，保证检测精度 * 频谱分析中的一些基本关系 采样率fs与低通截止频率fc：fs ≥ 2fc DFT点数N与窗长点数L：N≥L 频率f，角频率Ω，归一化角频率ω之间关系： ω = ΩT = 2πfT Ω = 2πf = ω/T f = Ω/2π = ω/2πT 采样频率 ωs=2π Ωs = 2π/T fs = 1/T 谱线频率间距 Δω=2π/N ΔΩ=Ωs/N Δf = fs/N 谱线对应频率 ωk = kΔω Ωk = kΔΩ fk = kΔf * 例：信号谱估计(习题10.6) xc(t)为带限5kHz的实信号，用10kHz采样后取1000点作DFT得X[k]。已知X[900]=1，X[420]=5，对尽可能多的Ω值求Xc(jΩ)。 解：∵ N=1000，T=1/fs=1/10000 ∴k=420， Ωk = 2πk/NT = 2π(4200) k=900， Ωk = Ωk-1000 = 2π(-1000) 上式用到DFT周期性，下式用对称性 ∴Xc(j2π(4200))= Xc(-j2π(4200))=X[420]T=5x10-4 Xc(j2π(1000))= Xc(-j2π(1000))=X[900]T=10-4 * 10.2 双余弦信号的频谱观察 工程上大量存在双频信号的实例，如积化和差 2cos(α)cos(β) = cos(α+β) + cos(α-β) 当 α β 就表征了通信调幅或载波调制， 当 α+β |α-β| 时代表了乐器调音的“差拍” 两个频率为ω0、ω1单频信号的合成信号，其傅里叶变换为δ(ω±ω0 )+δ(ω±ω1)，频谱为±ω0、±ω1处单一频线 满足采样定理离散化，频谱单周期(-π,π)内相同不变 * 加窗对频谱分析的影响 序列在时域通过加窗实现有限长，等效序列频谱与窗频谱在频域的卷积“模糊化”，单频时： V(ejω)=Aδ(ω-ω0)*W(ejω)=AW(ej(ω-ω0)) 单一频线与窗频谱卷积的结果，使单一频线改变为窗频谱的形状 可以选择窗频谱W(ejω)幅度为1并偶对称 单频谱峰值对应的频率、幅度仍精确的表征了原单一频线的频率、幅度！ * 双余弦序列加窗后频谱 频线相距较远相互影响很小 L=64的矩形窗频谱 3/8 δ(ω±2π/3) 1/2 δ(ω±2π/6) * 双频线靠拢后的效果 条件同前，但ω0不变，ω1减小逐渐向ω0靠近 频谱峰值的幅度偏离原频线的幅度 频谱峰值对应的频率也偏离原频线的频率 甚至双频线混淆不再可分辨 都是加窗惹的祸！ ω0=4π/28 ω1=4π/15 ω0=4π/28 ω1=4π/24 ω0=4π/28 ω1=4π/25 * 泄露—窗频谱旁瓣的影响 窗函数的选择与设计与FIR滤波器加窗法相同 注意窗长：L=M+1， M为FIR滤波器阶数 加窗引入的谱旁瓣对相邻谱峰 的幅度、中心都会产生误差 这种影响是相互的称为泄漏 关心频谱幅度测量、弱频谱 分量检测时泄漏影响很大 减小泄漏的方法是选择窗类 型、减小旁瓣 * 频线分辨—窗频谱主瓣的影响 加窗引入的谱主瓣对相邻谱峰的区分带来混淆 瑞利分辨率：二谱峰间峰谷幅度小于85%时的最小频线间隔Δω 窗谱主瓣Δωw越宽，分辨率越低 简化时，可取频率分辨率=Δωw Δωw与窗类型有关，旁瓣越低的 窗Δωw的相对宽度越大 Δωw与窗长L成反比 Δω 1 85% Δωw * 例：频谱分析时的加窗设计 50Hz工频电谐波测量，要求能确认比基波小50dB的谐波，测量精度达0.2% 解：一般周期梯形波，谐波10倍频程衰减40dB，故可取抗混叠低通1.6kHz，带外谐波小于60dB，采样频率3.2kHz。测量精度0.2%，并考虑量化及FFT误差，取16位ADC。 0.2%对应-54dB，查