盲解卷积算法-盲信号实验报告.pdf

盲信号实验报告 盲解卷积算法 姓 名： 丁宪成 系 别： 电信学院 专 业： 电磁场与微波 学 号： 3110035012 指导教师： 陈文超 2011年07月13 日 盲解卷积算法 1. 原理： 几个重要概念： 1.1 褶积模型假设： 假设 ：地层是由具有常速的水平层组成；1 2 (P ) 假设 ：震源产生一个平面压缩波 波 ，法向入射到层边界上，在 (S ) 这种情况下，不产生剪切波 波 ； 3 假设 ：震源波形在地下传播过程中不变，即它是稳定的； 假设 ：噪音成分是零；4 假设 ：震源波形是已知的；5 假设 ：反射系数序列是一个随机过程。这意味着地震记录具有地；6 震子波的特征，即它们的自相关和振幅谱是相似的； 7 假设 ：地震子波是最小相位的，因此，它有一个最小相位的逆。 1.2 反滤波 f(t), f(t) x(t) 如果定义滤波算子为 则 与已知地震记录 的褶积得到一个 e(t) 对地层脉冲响应 的估计 e(t)=f(t) x(t)∗ ； （）1 x(t)=w(t)*f(t)*x(t)； （）2 δ(t)=w(t)*f(t)； （）3 1 f ()t ()*t （4） wt() 用流程图表示为： 1.3 震源反子波 计算震源反子波在数学上是利用 变换来实现的。例如，假设基本子z 1 0.5) 波为两点时间序列（ ，－ 1 1 1 1 2 wz( )1 zF(z) 1 z z ...； （5） 2 1 2 4 1 z 2 1 1 1 2 F(z) 1 z z ... （6） 1 2 4 1 z 2 1 1 F(z)的系数(1, , ,...)代表逆滤波算子f(t)有关的时间序列。可以看 2 4 出它有无限多个系数，然而它们递减的很快。如同任何滤波过程一样， 实际应用的算子都是被截断的。 1.4 最小平方反滤波 当输入子波良好，其 变换的逆可以用一收敛序列表示，则上面所描z 述的反滤波将得到一个很好的近似于尖脉冲的输出将下面问题列出.