c语言地震子波6,地震子波反褶积程序

《地震子波反褶积程序》由会员分享，可在线阅读，更多相关《地震子波反褶积程序(19页珍藏版)》请在人人文库网上搜索。

1、盲信号实验报告盲解卷积算法姓 名： 丁宪成 系 别： 电信学院 专 业： 电磁场与微波 学 号： 3110035012 指导教师： 陈文超 2011年07月13日盲解卷积算法 1. 原理：几个重要概念：1.1 褶积模型假设：假设1：地层是由具有常速的水平层组成；假设2：震源产生一个平面压缩波(P波)，法向入射到层边界上，在这种情况下，不产生剪切波(S波)；假设3：震源波形在地下传播过程中不变，即它是稳定的；假设4：噪音成分是零；假设5：震源波形是已知的；假设6：反射系数序列是一个随机过程。这意味着地震记录具有地；震子波的特征，即它们的自相关和振幅谱是相似的；假设7：地震子波是最小相位的，因此，。

2、它有一个最小相位的逆。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 最小平方反滤波。

3、当输入子波良好，其z变换的逆可以用一收敛序列表示，则上面所描述的反滤波将得到一个很好的近似于尖脉冲的输出.将下面问题列出方程式：给定输入子波(10.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)求解。

4、滤波系数，我们得到(a,b)为(0.95,0.38)。1.5 最佳维纳滤波器回顾研究反滤波和最小平方滤波时讨论的期