盲源分离算法在混合震动信号分离中的应用盲源分离算法在混合震动信号分离中的应用
提出了一种基于z变换域有理传递函数F的时间延迟正定盲源分离算法,并提出将其应用于人工勘探地震波时传
感器采集到的混合震动信号的信噪分离及横、纵波的分离。该算法既适用于分离以线性方式混合的信号,也适
用于以非线性方式混合信号。时间延迟长短的选择依赖于有待处理的震动信号数据的长度。仿真结果表明,该
算法能有效地对Matlab生成的人工模拟震动波进行信噪分离及横、纵波的分离,为震动信号数据的后期处理及
分析提供有利依据。
1 基于时间延迟的盲源分离算法原理基于时间延迟的盲源分离算法原理
在混合震动信号中,不同种类的波是
1.1 数据分析数据分析
通常地震勘探所采用传感器采集到的数据都是间隔一定的采样时间所对应的电压幅值。将采样时间作为横轴,电压幅值作
为纵轴,就得到了震动信号的波形图,即时间序列波形数据。
将传感器采集到的电压幅值数据看做是一个1行p列(采样点数)的矩阵。则将所得到的n组数据组成一个新的矩阵x(n×p阶的
矩阵),其中n为观测信号个数,p为采样点数。
1.2 算法原理算法原理
无论是线性或者非线性盲源分离[3-4]算法,最终就是求得解混矩阵[5]w从而达到实现信号分离目的。首先必须明确的是式
(1)中给出的约束最小化损失函数模型:
其中,G′是函数G的导数,F′为F的导数。式(2)中,对于权重w的更新算式(3)所示:
由式(3)可知,函数G的确切形式对求取解混矩阵w并不起决定性的作用,而功能函数F的选择则直接影响w的数值及最终
的分离结果。在此提出一种z变换域有理传递函数:
1.3 算法实现的具体步骤算法实现的具体步骤
(1)对观测到的混合震动信号x(t)n×p进行归一化处理,得到新矩阵 (t)n×p中的各分量互不相关且其每个元素均是归一化的
单位方差。
(2)求步骤(1)中得到的归一化矩阵 T(t)的自协方差矩阵Un×n。
(3)求步骤(2)中得到的自协方差矩阵的特征向量及特征值矩阵Fn×n和Dn×n(对角元素为特征值,其他元素均为零),使得矩
阵U、F和D满足式(5):
U·F=F·D (5)
(4)求矩阵D的对角矩阵 n×n(D矩阵各对角元素的二次方根取倒数并保留实部后,按照大小顺序降序排列)。
(5)由步骤(3)中的矩阵F和步骤(4)中的矩阵 求得白化矩阵vn×n,使得三者满足如下关系式:
评论0
最新资源