盲信号实验报告..doc

盲信号实验报告 盲解卷积算法 姓 名： 丁宪成 系 别： 电信学院 专 业： 电磁场与微波 学 号： 3110035012 指导教师： 陈文超 2011年07月13日 盲解卷积算法 原理： 几个重要概念： 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） （4） 用流程图表示为： 1.3 震源反子波 计算震源反子波在数学上是利用z变换来实现的。例如，假设基本子波为两点时间序列（1，－0.5) ； （5） （6） 的系数代表逆滤波算子f(t)有关的时间序列。可以看出它有无限多个系数，然而它们递减的很快。如同任何滤波过程一样，实际应用的算子都是被截断的。 1.4 最小平方反滤波 当输入子波良好，其z变换的逆可以用一收敛序列表示，则上面所描述的反滤波将得到一个很好的近似于尖脉冲的输出.将下面问题列出方程式：给定输入子波(1－0.5)，寻找一个二项滤波器(a,b)，使实际输出与期望输出(1,0,0)之间的误差在最小二乘的意义上最小。将滤波器(a,b)与输入子波(1,-0.5)褶积以计算实际输出。误差的累积能量L定义为实际的与期望的输出的系数之差的平方和： （7） 目的是寻找系数(a,b)使L最小，这要求变量L随系数(a,b)而变并使之为零对上式进行简化，取L对a和b的偏导数，并使其为零，得到： （8） 有两个方程和两个未知数即滤波器系数(a,b)，可变成下列普通矩阵形式： （9） 求解滤波系数，我们得到(a,b)为(0.95,0.38)。 1.5 最佳维纳滤波器 回顾研究反滤波和最小平方滤波时讨论的期望输出，即零延迟尖脉冲(1,0,0)。重写最小平方反滤波的求解方程如下： （10） 两边同时除以2得到： 上述结论被维纳普遍化以推导出将输入转换为任意期望输出的滤波器(Robinson和Treltel，1980)。一个滤波器长度为n的矩阵方程的普遍形式是： （11） 式中、和(i＝0,l,2,3,…,n-1)分别为输入子波自相关、维纳滤波系数和期望输出与输入子波的互相关。其计算过程的流程图可以用下面图形表示为： 与最小平方滤波的联系：最佳维纳滤波器是最佳的，是指它的实际输出与期望输出之间的最小平方误差最小。当期望输出是零延迟尖脉冲(1,0,0,…,0)时，维纳滤波与最小平方滤波相同，即后者是前者的特例。 1.6 脉冲反褶积 第一类期望输出(零延迟尖脉冲)的处理称为脉冲反褶积。期望脉冲 (1,0,0,…,0)与输入子波的互相关得到序列(1,0,0,…,0)。正则方程的一般形式改写成特殊形式： (12) 方程已用作了归一化。在前面讨论过的最小平方反滤波器有与矩阵方程相同的形式。因此，脉冲反褶积在数学上是与最小平方反滤波相同的。而在实际上两类滤波器有一些区别。在脉冲反褶积(统计反褶积)的情况，方程在左边的自相关矩阵由输入地震记录(假设6)求得，而在最小平方反滤波(确定性反褶积)的情况则直接由已知的震源子波求得。 结论：如果输入子波不是最小相位，则脉冲反褶积不能将它