盲信号实验报告
盲解卷积算法
姓 名：系 别：专 业：学 号：指导教师： 2011年07月13日
丁宪成
电信学院
电磁场与微波
3110035012
陈文超
盲解卷积算法
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）
1f(t)??(t)* （4）
w(t)用流程图表示为：
1.3 震源反子波
计算震源反子波在数学上是利用z变换来实现的。例如，假设基本子波为两点时间序列（1，－0.5)
11121?1?z?z?...； （5） w(z)?1?zF(z)?12421?z2F(z)?111?1?z?z2?... （6） 1241?z211F(z)的系数(1,,,...)代表逆滤波算子f(t)有关的时间序列。可以看
24出它有无限多个系数，然而它们递减的很快。如同任何滤波过程一样，
地震子波反褶积程序
