数字信号的处理地几个前沿课地训练题目.pdf

实用标准文案 第 10章 数字信号处理的几个前沿课题 前面介绍了数字信号处理的基本知识， 本章我们将介绍时谱分析、 小波变换、 地震观测 系统仿真与地面运动恢复等几个数字信号前沿课题，以便大家在实际工作中参考。 10.1 时谱（倒谱）分析 时谱分析 (Cepstrum analysis) 是一种非线性信号处理技术，它在语言、图像、和噪声 处理领域中都有广泛的应用。 时谱可分为两类：复时谱和功率时谱。 MATLAB信号处理工具箱提供复时谱分析的工具 函数。 复时谱（ Complex cepstrum ）的定义为： 1 j j n ? x( n) ln X (e ) e d （10-1 ） 2 由上式可见，复时谱实际上是序列 x(n) 的 Fourier 变换取自然对数，再取 Fourier 逆 变换，得到的复时谱仍然是一个序列。也就是说，复时谱是 x(n) 从时间域至频率域、频率 域至频率域、频率域至时间域的三次变换。 MATLAB信号处理工具箱函数 cceps 用于估计一个序列 x 的复时谱，调用格式为： xhat=cceps(x) 式中， x 为输入序列（实序列） ；xhat 为复时谱（复序列） 。 MATLAB信号处理工具箱还提供了序列实时（倒）谱的计算程序 rceps ，调用格式为 Y=rceps(x) ，其中 x 为实序列； y 为实时谱，执行的操作为： 1 j Cx ln X (e ) d 2 （10-2 ） 由此可知，我们不能从序列 x 的实时谱重构原始序列，因为实时谱是根据序列 Fourier 变换的幅值计算的， 丢失了相位方面的信息。 但如果需要， 可采用最小相位模式估计原始序 列。 由于复时谱从复频谱计算得到， 不损失相位信息， 因此复时谱是可逆的， 实时谱过程是 不可逆的。 时谱分析技术广泛地应用于语言信号分析、同态滤波技术中。 这里举一个说明复时谱在具有回声信号测量中的应用。 【例 10-1 】设原信号是一个 45Hz 的正弦波，在传播过程中遇到障碍产生回声，回声振幅衰 减为原信号的 0.5 ，并与原信号有 0.2s 的延迟。在某测点测到的信号是原信号和回声信号 的叠加。试使用复时谱分析该测点的信号。 %Samp10_1 t=0:0.01:1.49; % 时间信号，采样间隔为 0.01s sig=sin(2*pi*45*t); % 原始信号 echo=0.5*[zeros(1,20),sig(1:length(t)-20)]; %在信号前面补加 20 个零，并使其振幅衰减一半作为回声信号 sigecho=sig+echo; 精彩文档 实用标准文案 c=cceps(sigecho); % 求原始信号与回声信号叠加的复时谱