2009_14_DSP的一些应用.pptVIP

复习：滤波器组合 滤波器是系统的一种。常见的滤波器组合可分为如下的三种形式： 1、串联：一个滤波器的输入为另一个滤波器的输出，即 * 复习：滤波器组合 滤波器是系统的一种。常见的滤波器组合可分为如下的三种形式： 2、并联：一个信号为多个滤波器的输入，系统的输出则为这些滤波器输出之和。 * 复习：滤波器组合 滤波器是系统的一种。常见的滤波器组合可分为如下的三种形式： 3、反馈（P75）： * IIR 相对于有限长脉冲响应（FIR）滤波器，递归滤波器（Recursive filter）又称为无限长脉冲响应滤波器（Infinite Impulse Response，简称IIR），简称IIR滤波器。在IIR滤波器的设计中，普遍采用的是离散序列的Z变换。 * IIR * DSP的一些具体应用 1、反褶 积； 2、最小平方滤波； 3、FIR与IIR； 4、二维或多维频谱计算。 * 二维信号的频谱计算 * 二维信号的频谱计算 在地震勘探中， １、垂直方向上，（每一道的）记录单位是时间（按向下递增方向排列）；频谱分析对应的单位是频率。 ２、水平方向上，每一道的间隔是长度；频谱分析对应的单位是波数。 * 二维信号的频谱计算 二维信号频谱计算全过程： １、对每道地震记录做频谱计算，即信号由时间－空间域变换到频率－空间域； 要点有二： （１）、每道都得做（一简单循环即可）； （２）、明白变换后高中低频所对应的部位。 * 二维信号的频谱计算 低频 低频 中频 中频 高频 * 二维信号的频谱计算 二维信号频谱计算全过程： ２、对每个频率所有各道做频谱计算，即信号由频率－空间域变换到频率－波数域； 要点有二： （１）、每个频率都得做； （２）、明白变换后高、中、低波数所对应的部位。 * 二维信号的频谱计算 低波数 中波数 高波数 低波数 中波数 * 二维信号的频谱计算 （不太准确：Spectrum2D.m） * 二维信号的频谱计算 　　上述两过程是可以调换的（即无先后之区别），所得到的结果是完全相同的。 重点是：明白其中高中低频率或波数所在的部位，以便于做滤波等处理。 同理可以对三维及更多维信号做频谱计算。 * DSP的一些具体应用 * DSP的授课内容 ）， CFST CFT SFT DFT (采样定理) 理论篇 重点: 1、计算； 2、性质的 推导。 应用篇 1、系统、Z变换（应用的理论基础）； 2、滤波、反褶积、多维信号的处理等。 前八周 后四周 * DSP的一些具体应用 ）， 褶积构成了DSP应用的理论基础： 1、系统是由固定的时间响应与输入系列的褶积来完成的； 2、滤波是系统的一种表现形式。 * DSP的一些具体应用 1、反褶 积； 2、最小平方滤波； 3、FIR与IIR； 4、二维或多维频谱计算。 * 反褶积 例1：不适定问题 “不适定问题”是个数学术语，通俗地讲，它是指病态、理论上具有多解的问题。 * 反褶积 例1：不适定问题 * 反褶积 例1：不适定问题 注意：可以利用计算机完整计算如此过程，得到唯一的、有意义的解。 * 反褶积 例2：反褶积 由于观察、记录误差等原因，导致这样的系统在实际工作中是存在的。 * 反褶积 例2：反褶积 * 反褶积 * 反褶积 * 反褶积 在所有可能的解当中，只有当满足 * 反褶积 也即意味着所得到的解最稳定（但不一定最正确）！ * 反褶积 数学上已经证明： 对于不适定系统而言，若其是对称（半）正定的，添加适当的正因素λ之后，系统所得到的解是最稳定的（某种程度上的最优解） * 反褶积 例3：循环反褶积 写成矩阵的形式，则有 * 反褶积 进一步地，可以得到 亦即系统对应的谱的表达式： * 反褶积 * 反褶积 * 反褶积 * 褶积运算的数学基础 从而导致（循环褶积的频谱表达式） * 褶积运算的数学基础 从而导致（循环褶积的频谱表达式） * 反褶积 进一步地， * 反褶积 上述三例中所添加的正数λ实际上是添加在系统的频谱上的；亦即添加了一个各种频率成分都存在且相等的信号，这样的信号就是白光，因此，我们称λ为白噪因子。 * 反褶积 例4：线性反褶积(有限长度) 对应的是一个Toeplitz矩阵系统。 对策： 2、对Toeplitz矩阵循环化之后， 求解对应的循环反褶积。 1、直接对Toeplitz矩阵求逆； * DSP的一些具体应用 1、反褶 积； 2、最小平方滤波； 3、FIR与IIR； 4、二维或多维频谱计算。 * 最小平方滤波 最小平方滤波，也称Wiener滤波或最佳滤波。 基本思想是：设计一个滤波器 使得待分析信号 经滤波处理后的输出 与期望输出