CN108282424B - 用于四个数据集联合盲源分离的四阶张量联合对角化算法 - Google Patents
用于四个数据集联合盲源分离的四阶张量联合对角化算法 Download PDFInfo
- Publication number
- CN108282424B CN108282424B CN201810079041.2A CN201810079041A CN108282424B CN 108282424 B CN108282424 B CN 108282424B CN 201810079041 A CN201810079041 A CN 201810079041A CN 108282424 B CN108282424 B CN 108282424B
- Authority
- CN
- China
- Prior art keywords
- matrix
- tensor
- algorithm
- joint
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 85
- 238000000926 separation method Methods 0.000 title claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 133
- 238000000034 method Methods 0.000 claims abstract description 24
- 230000002087 whitening effect Effects 0.000 claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000009466 transformation Effects 0.000 claims abstract description 9
- 238000002156 mixing Methods 0.000 claims description 28
- 238000005457 optimization Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 10
- 238000009795 derivation Methods 0.000 claims description 5
- 241000287196 Asthenes Species 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 239000004615 ingredient Substances 0.000 claims description 3
- 230000006870 function Effects 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 5
- 101001056160 Homo sapiens Methylcrotonoyl-CoA carboxylase subunit alpha, mitochondrial Proteins 0.000 description 3
- 102100026552 Methylcrotonoyl-CoA carboxylase subunit alpha, mitochondrial Human genes 0.000 description 3
- 238000010219 correlation analysis Methods 0.000 description 3
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 3
- 239000010977 jade Substances 0.000 description 3
- 101100126955 Arabidopsis thaliana KCS2 gene Proteins 0.000 description 2
- 241000132023 Bellis perennis Species 0.000 description 2
- 235000005633 Chrysanthemum balsamita Nutrition 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000001605 fetal effect Effects 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 238000007499 fusion processing Methods 0.000 description 2
- 238000012880 independent component analysis Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0272—Voice signal separating
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
- H04L25/03082—Theoretical aspects of adaptive time domain methods
- H04L25/03089—Theory of blind algorithms, recursive or not
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L2025/03592—Adaptation methods
- H04L2025/03598—Algorithms
- H04L2025/03611—Iterative algorithms
- H04L2025/03636—Algorithms using least mean square [LMS]
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L2025/03592—Adaptation methods
- H04L2025/03598—Algorithms
- H04L2025/03611—Iterative algorithms
- H04L2025/03649—Algorithms using recursive least square [RLS]
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种用于四个数据集联合盲源分离的四阶张量联合对角化算法,包括以下步骤:S1、观测信号:对四个数据集观测信号分别进行预白化;S2、目标张量:预白化之后,构造一组互四阶累积量张量;S3、初始化因子矩阵;S4、代价函数收敛计算:若一次扫描过后,算法收敛,则计算结束,若算法仍未收敛,则以此次更新所得的因子矩阵作为初始值,进行下一次扫描,遍历更新雅克比旋转矩阵,更新因子矩阵,直至收敛为止。本发明所述的用于四个数据集联合盲源分离的四阶张量联合对角化算法,该算法基于正交旋转变换,因而能够获得最小二乘意义上的最优解,是一种针对不多于四个数据集信号的J‑BSS方法。
Description
技术领域
本发明涉及生物医学、通信、语音、阵列等信号处理领域,一种可用于四个数据集联合盲源分离的四阶张量联合对角化算法。
背景技术
联合盲源分离(Joint Blind Source Separation,J-BSS)作为一种新兴的、基于数据驱动的多数据集融合处理技术,在生物医学、通信、语音、阵列等信号处理领域获得了大量关注[1]–[14]。
同传统的盲源分离(BSS)不同[14]–[24],J-BSS通常假设如下多数据集信号模型:
x(m)(t)=A(m)s(m)(t),m=1,2,...,M, (1A)
其中,x(t)(m)∈CN和s(t)(m)∈CR分别表示第m个数据集的观测信号和源信号, A(m)∈CN×R表示第m个数据集的混合矩阵。
参数M、N、R分别表示数据集个数、观测信号通道数及源信号个数。
J-BSS旨在仅已知多数据集观测信号x(m)(t)的前提下,对多数据集混合矩阵 A(m)及源信号s(m)(t)进行联合辨识。
在上述模型中,多数据集信号往往是同一目标或物理现象,经由不同的观测方式所获得的多组数据。
在这里,“不同的观测方式”是不同的设备或实验方案,如脑电(EEG)与功能磁共振(fMRI)[2];或是同一种设备和方案下,不同的被试个体,如多被试fMRI[3],[7];或者是相同的设备、方案及被试个体,在不同变换域的观测 (如频域、时域、空域、统计域等[3],[10])。
在处理上述多数据集信号时,由于J-BSS利用了不同集合间的相关性及相异性,经融合处理后的信息能够更为全面地描述观测对象在多层次、多剖面、多模态意义上的特征,因此在辨识能力、分离精度等方面之BSS具有性能优势。
关于J-BSS的研究最早可追溯至上世纪30年代,文献[1]所提出的典范相关分析方法(Canonical Correlation Analysis:CCA)。在后来的较长一段时间内, CCA被广泛用于两个数据集信号的融合预处理之中。
近些年来,CCA被拓展至多于两个数据集的情形,即多集合CCA(Multiset CCA:M-CCA)[6],并在多被试fMRI数据融合中得以应用。随后,M-CCA 的耦合矩阵分解模型进一步演变为矩阵广义联合对角化(Generalized joint diagonalization:GJD)模型[8]。
此外,由经典的独立成分分析(Independent Component Analysis:ICA)发展而来的独立矢量分析(Independent Vector Analysis:IVA)也是J-BSS的主要方法之一[26]–[28]。
最近一两年,面向多集合、多维源信号分离问题的联合独立子空间分析 (jointindependent subspace analysis:J-ISA)逐步被人们所关注[29]。
在上述研究中,GJD是J-BSS的一类重要方法。这类方法通过计算互累积量,将多数据集信号转化为一组相互耦合,且具有可联合对角化结构的矩阵切片,从而通过对上述结构进行代数拟合实现J-BSS。
目前,有关GJD的研究工作大多基于信号二阶互协方差。而仅有的基于四阶互累积量的方法[8],也只是考虑了四阶互累积量在矩阵层面的低阶可联合对角化结构,未能充分挖掘其高阶可联合对角化结构。此外,该方法并未对信号的时域特性(如非平稳性)加以利用,进而导致算法的精准度不足。
发明内容
根据上述提出的技术问题,而提供一种用于四个数据集联合盲源分离的四阶张量联合对角化算法,用以弥补现有技术未利用信号联合高阶统计特性的缺陷,进而导致算法的精准度不足的缺点。本发明采用的技术手段如下:
一种用于四个数据集联合盲源分离的四阶张量联合对角化算法,包括以下步骤:
S1、观测信号:
对四个数据集观测信号分别进行预白化:
记第m个数据集的白化矩阵为W(m)∈CR×N,m=1,…,4;
则预白化处理之后,第m个数据集观测信号为:
y(m)(t)=W(m)x(m)(t)=U(m)s(m)(t)∈CR, (1)
其中,酉矩阵U(m)@W(m)A(m)∈CR×R是经预白化之后,第m个数据集信号的混合矩阵;
S2、目标张量:
预白化之后,构造一组互四阶累积量张量如下:
Tk==cum(y(1)(tk),y(2)(tk),y(3)(tk),y(4)(tk)),k=1,...,K, (2)
其中,“cum(·)”用于计算互四阶累积量,其定义为:
其中,“E(·)”用于计算数学期望,上标“*”表示共轭,根据定义可知,每一个张量Tk∈CR×R×R×R表示四组观测信号在时刻tk的互四阶累积量张量,k= 1,…,K;
其公式表示如下:
Tk=Dk×1U(1)×2U(2)*×3U(3)*×4U(4), (4)
公式(4)中Dk=cum(s(1)(tk),s(2)(tk),s(3)(tk),s(4)(tk))∈CR×R×R×R为四组数据集的源信号在时刻tk的互四阶累积量张量;
S3、初始化因子矩阵;
S4、代价函数收敛计算:
首先,优化准则:
定义目标函数λ为数据张量Tk与拟合张量Dk×1U(1)×2U(2)*×3U(3)*×4U(4)间的均方误差:
其中||·||F表示弗罗贝尼乌斯范数,通过最小化λ,计算混合矩阵在最小二乘意义下的最优估计值;
注意到U(m)为酉矩阵,m=1,2,3,4,根据酉变换的保范性,公式(7)改写为:
其中,off(·)用于将其输入张量的超对角元置0,diag(·)用于对其输入张量的非对角元置0,为常数。因此,对λ进行最小化等价于对进行最大化,因此本发明通过如下准则估计混合矩阵U(m),m=1,2,3,4:
公式(7)–(9)表明,通过最大化Tk×1U(1)H×2U(2)T×3U(3)T×4U(4)H的超对角元范数平方和,即对T1,...,TK进行联合对角化,获得混合矩阵U(1),U(2),U(3),U(4)在最小二乘意义下的最优估计值;
其次,雅克比迭代:
雅克比迭代将待求解的酉矩阵写成一系列雅克比旋转矩阵的乘积,进而通过分别优化每一个雅克比旋转矩阵,最终实现目标函数的最大化;
若一次扫描过后,算法收敛,则计算结束,若算法仍未收敛,则以此次更新所得的因子矩阵作为初始值,进行下一次扫描,遍历更新雅克比旋转矩阵,更新因子矩阵,直至收敛为止。
作为优选在步骤S1之前,做如下假设:
(A3)源信号为非高斯、非平稳信号;
(A4)所考虑的多数据集模型为正定的,即N≥R;
(A5)数据集的数目M=4。
作为优选在实际中,在假设条件(A1)和(A4)下,白化矩阵W(m)∈CR×N可通过观测信号x(m)(t)的二阶协方差矩阵的奇异值分解而获得。
作为优选由假设(A1),(A2)及(A3)可知,每一个张量Dk均为超对角张量,即:
定义U′(m)@U(m)D(m)P,,m=1,2,3,4,其中D(1),D(2),D(3),D(4)∈CR×R为对角矩阵,且满足D(1)D(2)*D(3)*D(4)=I,I为单位矩阵,P∈CR×R为排序矩阵;
不难得知,若将公式(4)中的U(m)替换为U′(m),等式依然成立;
因此,U′(m)在不考虑幅度/相位模糊以及顺序模糊时,和混合矩阵U(m)等价;
通过对方程(4)进行求解,可对混合矩阵U(m)进行估计。
雅克比旋转矩阵的闭式最优解
因此,优化准则(12)可改写为:
其中
首先根据定义可知:
其中,
根据公式(15),有:
其中Mi,j=[m1,i,j,...,mK,i,j,q1,i,j,...,qK,i,j]∈C3×2K;
作为优选迭代停止条件的判定为:
计算相邻两次扫描所获得张量的对角元范数和的相对变化值:
当其小于阈值ε时,则认为算法收敛,迭代终止;
或,判断当前扫描所获得的雅克比旋转矩阵是否接近于单位矩阵,具体而言,定义当ξ(1)、ξ(2)、ξ(3)和ξ(4)同时小于阈值ε时,则说明和已经近似于单位阵,不能更新目标张量,可认为算法收敛,迭代终止。
与现有技术相比较,本发明所述的用于四个数据集联合盲源分离的四阶张量联合对角化算法,是一种针对不多于四个数据集信号的J-BSS方法。该方法对多数据集信号计算其在不同时间段的互四阶累积量张量,并通过对其进行张量联合对角化,实现多数据集信号的J-BSS。
此外,我们提出了一种基于Givens连续旋转的张量联合对角化算法。同现有技术[30]中提出的仅适用于单数据集BSS的张量对角化算法相比,该算法适用于张量的加载因子矩阵不相同的情况,从而适用于J-BSS。同现有技术[8]中提出的非正交张量对角化算法相比,该算法基于正交旋转变换,因而能够获得最小二乘意义上的最优解。
本发明所述的用于四个数据集联合盲源分离的四阶张量联合对角化算法,应用了信号的时域特性,本发明对多数据集信号计算其在不同时间段的互四阶累积量张量,提高了算法的精准度。
本发明所述的用于四个数据集联合盲源分离的四阶张量联合对角化算法,利用四个因子矩阵构建了张量,在算法的过程中,同时交替更新四个因子矩阵,算法精准,同现有技术中的算法只是交替更新前三个因子矩阵相比,第四个因子矩阵在得出前三个因子矩阵后,自动求出。
本发明的优化准则直接对四个因子矩阵进行优化,现有技术的算法只优化前三个因子矩阵。由于优化准则不同,且被优化的张量也不同,所以后续的优化过程都是不同的,现有技术中由于前三个因子矩阵已被求出,故而第四个因子矩阵通过还原目标张量而获得,但是这样获得的因子矩阵精度低。本发明直接优化四个因子矩阵,提高了第四个因子矩阵的精度。
本发明在计算张量时,将信号分为k段,共计算出k个张量,同现有技术中算法信号不分段,计算出n个张量,也就是同因子矩阵的维度个张量相比,利用了更多的信号信息,从而提高算法精度。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明算法流程图。
图2-1是本发明实施例1中无噪声条件下JTD和GOJD算法的A-PI值随扫描次数的变化曲线(K=3,N=3)。
图2-2是本发明实施例1中无噪声条件下JTD和GOJD算法的A-PI值随扫描次数的变化曲线(K=3,N=50)。
图2-3是本发明实施例1中无噪声条件下JTD和GOJD算法的A-PI值随扫描次数的变化曲线(K=50,N=3)。
图2-4是本发明实施例1中无噪声条件下JTD和GOJD算法的A-PI值随扫描次数的变化曲线(K=50,N=50)。
图3-1是本发明实施例2中A-PI随SNR变化曲线示意图(快拍数=100000)。
图3-2是本发明实施例2中A-PI随快拍数变化曲线示意图(SNR=5dB)。
图4本发明实施例3中八路ECG数据图。
图5是本发明实施例3中实际FECG信号分离结果图。
具体实施方式
如图1所示,一种用于四个数据集联合盲源分离的四阶张量联合对角化算法,首先吗,做如下假设:
(A3)源信号为非高斯、非平稳信号;
(A4)所考虑的多数据集模型为正定的,即N≥R;
(A5)数据集的数目M=4。
在上述假设中,(A1)和(A2)是现有的联合盲源分离方法中普遍用到的假设条件;条件(A3)和(A4)则表明本发明所述方法适用于正定情况下,非高斯及非平稳信号的联合盲源分离;条件(A5)则限定本发明方法所能处理的数据集数目不超过四个,本发明以四个数据集为例阐述所提出的联合盲源分离算法,低于四个数据集的情况通过类似的推导获得。
具体包括以下步骤:
S1、观测信号:
对四个数据集观测信号分别进行预白化:
记第m个数据集的白化矩阵为W(m)∈CR×N,m=1,…,4;
则预白化处理之后,第m个数据集观测信号为:
y(m)(t)=W(m)x(m)(t)=U(m)s(m)(t)∈CR, (1)
其中,酉矩阵U(m)@W(m)A(m)∈CR×R是经预白化之后,第m个数据集信号的混合矩阵;在实际中,在假设条件(A1)和(A4)下,白化矩阵W(m)∈CR×N可通过观测信号x(m)(t)的二阶协方差矩阵的奇异值分解而获得[8]。
S2、目标张量:
预白化之后,构造一组互四阶累积量张量如下:
Tk==cum(y(1)(tk),y(2)(tk),y(3)(tk),y(4)(tk)),k=1,...,K, (2)
其中,“cum(·)”用于计算互四阶累积量,其定义为:
其中,“E(·)”用于计算数学期望,上标“*”表示共轭,根据定义可知,每一个张量Tk∈CR×R×R×R表示四组观测信号在时刻tk的互四阶累积量张量,k= 1,…,K;
其公式表示如下:
Tk=Dk×1U(1)×2U(2)*×3U(3)*×4U(4), (4)
公式(4)中Dk=cum(s(1)(tk),s(2)(tk),s(3)(tk),s(4)(tk))∈CR×R×R×R为四组数据集的源信号在时刻tk的互四阶累积量张量;
由假设(A1),(A2)及(A3)可知,每一个张量Dk均为超对角张量,即:
定义U′(m)@U(m)D(m)P,,m=1,2,3,4,其中D(1),D(2),D(3),D(4)∈CR×R为对角矩阵,且满足D(1)D(2)*D(3)*D(4)=I,I为单位矩阵,P∈CR×R为排序矩阵;
不难得知,若将公式(4)中的U(m)替换为U′(m),等式依然成立;
因此,U′(m)在不考虑幅度/相位模糊(由对角矩阵D(m)表征)以及顺序模糊 (由排序矩阵P表征)时,和混合矩阵U(m)等价。通过对方程(4)进行求解,可对混合矩阵U(m)进行估计。
联合盲源分离的目的,正是在不考虑上述两类平凡模糊的前提条件下,对混合矩阵进行辨识,进而对源信号进行估计。
同传统的BSS有所不同的是,当m取不同索引值时,不同的U′(m)相对于真实混合矩阵U(m)的排序矩阵P是完全相同的,这意味着经由J-BSS所估计的混合矩阵及源信号成分,在顺序上是自然对齐的,这也是J-BSS相较于BSS的优势之一。
本发明中利用四个因子矩阵构建了张量,在算法的过程中,同时交替更新四个因子矩阵,同现有技术中的算法只是交替更新前三个因子矩阵相比,第四个因子矩阵在得出前三个因子矩阵后,自动求出。
本发明在计算张量时,将信号分为k段,共计算出k个张量,同现有技术中算法信号不分段,计算出n个张量,也就是因子矩阵的维度个张量相比,利用了更多的信号信息,从而提高算法精度。
S3、初始化因子矩阵;
S4、代价函数收敛计算:
首先,优化准则:
定义目标函数λ为数据张量Tk与拟合张量Dk×1U(1)×2U(2)*×3U(3)*×4U(4)间的均方误差:
其中||·||F表示弗罗贝尼乌斯范数(Frobenius norm),通过最小化λ,计算混合矩阵在最小二乘意义下的最优估计值;
注意到U(m)为酉矩阵,m=1,2,3,4,根据酉变换的保范性,公式(7)改写为:
其中,off(·)用于将其输入张量的超对角元置0,diag(·)用于对其输入张量的非对角元置0,为常数。因此,对λ进行最小化等价于对进行最大化,因此本发明通过如下准则估计混合矩阵U(m),m=1,2,3,4:
公式(7)–(9)表明,通过最大化Tk×1U(1)H×2U(2)T×3U(3)T×4U(4)H的超对角元范数平方和,即对T1,...,TK进行联合对角化,获得混合矩阵U(1),U(2),U(3),U(4)在最小二乘意义下的最优估计值;接下来,我们提出一种基于雅克比迭代的张量联合对角化算法。
本发明的优化准则直接对四个因子矩阵进行优化,现有技术中算法只优化前三个因子矩阵。由于优化准则不同,且被优化的张量也不同,所以后续的优化过程都是不同的,现有技术中由于前三个因子矩阵已被求出,故而第四个因子矩阵通过还原目标张量而获得,但是这样获得的因子矩阵精度低。本发明直接优化四个因子矩阵,提高了第四个因子矩阵的精度。
其次,雅克比迭代:
雅克比迭代将待求解的酉矩阵写成一系列雅克比旋转矩阵的乘积,进而通过分别优化每一个雅克比旋转矩阵,最终实现目标函数的最大化;
若一次扫描过后,算法收敛,则计算结束,若算法仍未收敛,则以此次更新所得的因子矩阵作为初始值,进行下一次扫描,遍历更新雅克比旋转矩阵,更新因子矩阵,直至收敛为止。
迭代停止条件的判定为:
计算相邻两次扫描所获得张量的对角元范数和的相对变化值:
当其小于阈值ε时,则认为算法收敛,迭代终止;
或,判断当前扫描所获得的雅克比旋转矩阵是否接近于单位矩阵,具体而言,定义当ξ(1)、ξ(2)、ξ(3)和ξ(4)同时小于阈值ε时,则说明和已经近似于单位阵,不能更新目标张量,可认为算法收敛,迭代终止。
从上述描述可知,雅克比迭代将优化问题(9),化解为一系列子优化问题 (12)。由于每一个雅克比旋转矩阵仅有两个未知参数,其最优解有望具有简单的闭式解形式。接下来,我们将推导每一步迭代中所涉及的闭式最优解。
雅克比旋转矩阵的闭式最优解
因此,优化准则(12)可改写为:
其中
首先根据定义可知:
其中
根据公式(15),有:
其中Mi,j=[m1,i,j,...,mK,i,j,q1,i,j,...,qK,i,j]∈C3×2K;
上述推导假设J-BSS问题为复值,当J-BSS问题为实值时,亦通过类似的计算将其转化为实值四阶张量的联合对角化问题,并通过实值雅克比旋转对之进行求解,计算步骤与复值情况类似。
需要注意的是,此时实值雅克比旋转矩阵定义为:
算法性能通过如下式定义的性能指标(Performance Index:PI)进行评价:
通过计算每一次扫描结束时算法所对应的PI值,可反映算法的收敛性。图 2-1至图2-4绘制了在K和N取不同值时,JTD算法在10次独立实验中,其 A-PI值随扫描次数的变化曲线。
作为参考,我们绘制在相同条件下,广义正交矩阵对角化算法(GOJD,[8]) 的A-PI曲线(其数据矩阵由Tk=固定后三维索引获得)。看到,所提出的JTD 算法呈现出单调下降、线性收敛的收敛模式。同GOJD算法相比,JTD算法平均5次扫描即能达到收敛,具有更优的收敛性能。
图2-1至图2-4,为无噪声条件下JTD和GOJD算法的A-PI值随扫描次数的变化曲线。这里展示四种条件下的实验结果:(a)如图2-1所示,张量的尺寸N和数量K均取较小值;(b)如图2-2所示,张量的尺寸N取较大值,数量 K取较小值;(c)如图2-3所示,张量的尺寸N取较小值,数量K取较大值; (d)如图2-4所示,张量的尺寸N和数量K均取较大值。
实施例2,本实验中,根据公式(1A)产生四组观测信号x(m)(t),m=1,2,3,4:
其中混合矩阵A(m)∈CN×R由程序随机产生,噪声项n(m)(t)∈CN×T为零均值单位方差的高斯白噪声。
源信号由下式构造:
sr(t)=Qrs′r(t), (22)
其中s′r(t)∈C4为分段调幅的随机相位信号,其相位及调幅系数由程序随机产生。矩阵Qr∈R4×4由程序随机产生,用于引入集合间的相关性。不难得知,所构造的源信号为非平稳非高斯信号(假设 A3),且不同数据集的源信号满足组间相关和组内独立的假设条件(假设A1, A2)。
在公式(21)中,σs和σn分别代表信号强度和噪声强度。进一步可定义信噪比(SNR)如下:SNR@20lgσs/σn。
我们用所提出的基于JTD的J-BSS算法对上述多数据集信号进行分离。并将其同下述同类型的BSS及J-BSS算法进行对比:
(a)特征矩阵近似联合对角化算法(Joint Approximate Diagonalization ofEigenmatrices:JADE)[16];
(b)三阶张量联合对角化算法(Simultaneous Third-Order TensorDiaognalization:STOTD)[30];
(c)多集合典范相关分析算法(Multiset Canonical Correlation Analysis:MCCA)[6];
(d)广义正交矩阵联合对角化算法(Generalized Orthogonal JointDiagonalization:GOJD)[8]。
其中JADE和STOTD为基于四阶累积量的BSS算法,MCCA和GOJD为 J-BSS算法。在GOJD算法的实验中,我们分别采用了基于四阶累积量和二阶协方差的算法,分别记为GOJD4和GOJD2。
作我们用采样互四阶累积量为(2)所定义的互四阶累积量的估计值(2),其计算公式如下:
在本实验中,我们固定各数据集的观测通道数N=10,源信号数为R=5, A-PI曲线中的每个点由200次蒙特卡洛独立实验的结果平均所获得。首先,固定时间采样点数T=100000,L=T/10,α=0.3。
令SNR从0dB变化至10dB。接下来,固定SNR=5dB,令时间采样点数由10000变化至100000。在上述两种情况下,所比较算法的A-PI曲线如图3-1 和图3-2所示。看到,所提出的JTD算法性能明显优于除GOJD2之外的其它算法。
同GOJD2相比,由图3-1看出,在SNR为2–8dB时,JTD算法的分离精度更高。图3-2则表明,当快拍数较高时(80000–100000)所提出的JTD算法精度更高,反之则是GOJD2更具优势。
这主要是因为GOJD2所采用的二阶协方差较之JTD所采用的的四阶累积量,在短快拍的情况下,具有更小的有限采样误差。
如图3-1和图3-2所示,JTD,GJD,MCCA,JADE,STOTD分离计算机合成信号时的性能对比。
实施例3,实验旨在从实际采集的,混合了孕妇心电信号(Mother ECG, MECG),胎儿心电信号(Fetal ECG,FECG),及其它干扰信号的观测数据中,提取FECG信号。
采用的八路ECG数据取自鲁汶大学的DAISY数据库0,信号波形如图4 所示。采样频率为250Hz,采样点数为2500,采样时间是10s,如图4所示,是实际采集的八路ECG数据,数据取自鲁汶大学的DAISY数据库。
按照下式,将八路观测信号划分为四个数据集,每个数据集中含有五路观测信号:
将四个数据集的观测信号预白化后分段,每段长度为L=500,重叠率为α=0.5,不难得知此时K=10。对每一段信号计算互四阶累积量以构建K个目标张量。用JTD算法对其进行联合对角化,获得解混矩阵后,由公式(1)可获得各源信号成分,结果如图4所示,表示第m个数据集的源信号的估计值,如图5所示,能够看到MECG信号和FECG信号被成功分离。
从图中看出,每个数据集第1路信号和第2路信号为MECG,前三个数据集中第3路信号为FECG。结果表明,基于JTD的J-BSS算法能够成功分离实际采集的FECG信号。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
参考文献
[1]Hotelling H.Relations between two sets of variates[J].Biometrika,1936,28(3/4):321-377.
[2]Sui J,Adali T,Yu Q,et al.A review of multivariate methods formultimodal fusion of brain imaging data[J].Journal of Neuroscience Methods,2012,204(1):68-81.
[3]Kim T,Attias H T,Lee S Y,et al.Blind source separation exploitinghigher-order frequency dependencies[J].IEEE Trans on Audio,Speech andLanguage Processing,2007,15(1): 70-79.
[4]Lee J H,Lee T W,Jolesz F A,et al.Independent vector analysis(IVA):multivariate approach for fMRI group study[J].Neuroimage,2008,40(1):86-109.
[5]Anderson M,Adali T,Li X L.Joint blind source separation withmultivariate Gaussian model:Algorithms and performance analysis[J].IEEE Transon Signal Processing,2012, 60(4):1672-1683.
[6]Li Y O,Adali T,Wang W,et al.Joint blind source separation bymultiset canonical correlation analysis[J].IEEE Trans on Signal Processing,2009,57(10):3918-3929.
[7]Correa N M,Adali T,Li Y O,et al.Canonical correlation analysis fordata fusion and group inferences[J].IEEE Signal Processing Magazine,2010,27(4):39-50.
[8]Li X L,Adal1 T,Anderson M.Joint blind source separation bygeneralized joint diagonalization of cumulant matrices[J].Signal Processing,2011,91(10):2314-2322.
[9]Congedo M,Phlypo R,Chatel-Goldman J.Orthogonal and non-orthogonaljoint blind source separation in the least-squares sense[A].Signal ProcessingConference(EUSIPCO)[C]. Bucharest,Romania,2012:1885-1889.
[10]Gong X F,Wang X L,Lin Q H.Generalized Non-Orthogonal JointDiagonalization With LU Decomposition and Successive Rotations[J].IEEE Transon Signal Processing,2015, 63(5):1322-1334.
[11]Lahat D,Jutten C.Joint independent subspace analysis:a quasi-Newton algorithm[A]. International Conference on Latent Variable Analysis andSignal Separation[C].Liberec, Czech Republic,2015:111-118.
[12]Lahat D,Adali T,Jutten C.Multimodal data fusion:an overview ofmethods,challenges,and prospects[J].Proceedings of the IEEE,2015,103(9):1449-1477.
[13]Adali T,Levin-Schwartz Y,Calhoun V D.Multimodal data fusion usingsource separation: Two effective models based on ICA and IVA and theirproperties[J].Proceedings of the IEEE, 2015,103(9):1478-1493.
[14]Hérault J,Jutten C,Ans B.Détection de grandeurs primitives dansun message composite par une architecture de calcul neuromimétique enapprentissage non supervisé[A]. Proceedings of the 10th Workshop Traitementdu signal et ses applications(GRETSI)[C]. Nice,France,1985:1017–1022.
[15]Comon P.Independent component analysis,a new concept?[J].SignalProcessing,1994, 36(3):287-314.
[16]Cardoso J F,Souloumiac A.Blind beamforming for non-Gaussiansignals[A].IEE Proceedings F(Radar and Signal Processing)[C].1993,140(6):362-370.
[17]Oja E,Hyvarinen A.A fast fixed-point algorithm for independentcomponent analysis[J]. Neural Computation,1997,9(7):1483-1492.
[18]Belouchrani A,Abed-Meraim K,Cardoso J F,et al.A blind sourceseparation technique using second-order statistics[J].IEEE Trans on SignalProcessing,1997,45(2):434-444.
[19]Novey M,Adali T.Complex ICA by negentropy maximization[J].IEEETrans on Neural Networks,2008,19(4):596-609.
[20]Sawada H,Araki S,Makino S.Underdetermined convolutive blindsource separation via frequency bin-wise clustering and permutation alignment[J].IEEE Trans on Audio,Speech and Language Processing,2011,19(3):516-527.
[21]Yu Y,Petropulu A P.PARAFAC-based blind estimation of possiblyunderdetermined convolutive MIMO systems[J].IEEE Trans on Signal Processing,2008,56(1):111-124.
[22]Lin Q H,Liu J,Zheng Y R,et al.Semiblind spatial ICA of fMRI usingspatial constraints[J].Human Brain Mapping,2010,31(7):1076-1088.
[23]Cichocki A,Zdunek R,Phan A H,et al.Nonnegative matrix and tensorfactorizations: applications to exploratory multi-way data analysis and blindsource separation[M].John Wiley&Sons,2009.
[24]Handbook of Blind Source Separation:Independent componentanalysis and applications[M].Academic Press,2010.
[25]Li Y O,Wang W,Adali T,et al.CCA for joint blind source separationof multiple datasets with application to group fMRI analysis[A].IEEEInternational Conference on Acoustics, Speech and Signal Processing[C].2008:1837-1840.
[26]Kim T,Eltoft T,Lee T W.Independent vector analysis:An extensionof ICA to multivariate components[A].International Conference on IndependentComponent Analysis and Signal Separation[C].Springer Berlin Heidelberg,2006:165-172.
[27]Anderson M,Li X L,Adal1 T.Nonorthogonal independent vectoranalysis using multivariate Gaussian model[A].International Conference onLatent Variable Analysis and Signal Separation[C].Springer Berlin Heidelberg,2010:354-361.
[28]Anderson M,Adali T,Li X L.Joint blind source separation withmultivariate Gaussian model:Algorithms and performance analysis[J].IEEE Transon Signal Processing,2012, 60(4):1672-1683.
[29]Lahat D,Jutten C.Joint independent subspace analysis:a quasi-Newton algorithm[C]. International Conference on Latent Variable Analysis andSignal Separation[C].Liberec, Czech Republic,2015:111-118.
[30]De Lathauwer L,De Moor B,Vandewalle J.Independent componentanalysis and (simultaneous)third-order tensor diagonalization[J].IEEE Transon Signal Processing,2001, 49(10):2262-2271.
Moor D.Database for the identification of systems(DaISy)[Z].http://www.esat.kuleuven.ac. be/sista/daisy,2010。
Claims (6)
1.一种用于四个数据集联合盲源分离的四阶张量联合对角化算法,其特征在于包括以下步骤:
S1、观测信号:
对四个数据集观测信号分别进行预白化:
则预白化处理之后,第m个数据集观测信号为:
其中,酉矩阵U(m)@W(m)A(m)∈CR×R是经预白化之后,第m个数据集信号的混合矩阵,x(m)(t)为第m个数据集的观测信号,s(m)(t)为第m个数据集的源信号,A(m)为第m个数据集的混合矩阵;
S2、目标张量:
预白化之后,构造一组互四阶累积量张量如下:
其中,“cum(·)”用于计算互四阶累积量,其定义为:
其公式表示如下:
S3、初始化因子矩阵;
S4、代价函数收敛计算:
首先,优化准则:
其中||·||F表示弗罗贝尼乌斯范数,通过最小化λ,计算混合矩阵在最小二乘意义下的最优估计值;
注意到U(m)为酉矩阵,m=1,2,3,4,根据酉变换的保范性,公式(7)改写为:
其次,雅克比迭代:
雅克比迭代将待求解的酉矩阵写成一系列雅克比旋转矩阵的乘积,进而通过分别优化每一个雅克比旋转矩阵,最终实现目标函数的最大化;
若一次扫描过后,算法收敛,则计算结束,若算法仍未收敛,则以此次更新所得的因子矩阵作为初始值,进行下一次扫描,遍历更新雅克比旋转矩阵,更新因子矩阵,直至收敛为止。
5.根据权利要求1所述的用于四个数据集联合盲源分离的四阶张量联合对角化算法,其特征在于:
雅克比旋转矩阵的闭式最优解求解如下:
因此,优化准则(12)可改写为:
首先根据定义可知:
其中,
根据公式(15),有:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810079041.2A CN108282424B (zh) | 2018-01-26 | 2018-01-26 | 用于四个数据集联合盲源分离的四阶张量联合对角化算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810079041.2A CN108282424B (zh) | 2018-01-26 | 2018-01-26 | 用于四个数据集联合盲源分离的四阶张量联合对角化算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108282424A CN108282424A (zh) | 2018-07-13 |
CN108282424B true CN108282424B (zh) | 2021-02-12 |
Family
ID=62805422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810079041.2A Active CN108282424B (zh) | 2018-01-26 | 2018-01-26 | 用于四个数据集联合盲源分离的四阶张量联合对角化算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108282424B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111366905B (zh) * | 2020-04-12 | 2023-09-01 | 南京理工大学 | 空间微动群目标多通道盲源分离方法 |
CN113114597B (zh) * | 2021-03-25 | 2022-04-29 | 电子科技大学 | 一种基于jade算法的四输入信号分离方法及*** |
CN116776108A (zh) * | 2023-06-14 | 2023-09-19 | 中国人民解放军空军预警学院 | 一种基于三阶累积量和张量分解的欠定联合盲源分离方法与*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103780522A (zh) * | 2014-01-08 | 2014-05-07 | 西安电子科技大学 | 基于双重迭代的非正交联合对角化瞬时盲源分离方法 |
CN105282067A (zh) * | 2015-09-16 | 2016-01-27 | 长安大学 | 一种复数域盲源分离方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102903368B (zh) * | 2011-07-29 | 2017-04-12 | 杜比实验室特许公司 | 用于卷积盲源分离的方法和设备 |
-
2018
- 2018-01-26 CN CN201810079041.2A patent/CN108282424B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103780522A (zh) * | 2014-01-08 | 2014-05-07 | 西安电子科技大学 | 基于双重迭代的非正交联合对角化瞬时盲源分离方法 |
CN105282067A (zh) * | 2015-09-16 | 2016-01-27 | 长安大学 | 一种复数域盲源分离方法 |
Non-Patent Citations (1)
Title |
---|
Jacobi-type Joint Diagonalization for Complex Symmetric Matrices;WANG Ke;《新型工业化》;20131231;第3卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108282424A (zh) | 2018-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Adali et al. | Diversity in independent component and vector analyses: Identifiability, algorithms, and applications in medical imaging | |
Li et al. | Joint blind source separation by generalized joint diagonalization of cumulant matrices | |
Aminghafari et al. | Multivariate denoising using wavelets and principal component analysis | |
Zhou et al. | Linked component analysis from matrices to high-order tensors: Applications to biomedical data | |
Bobin et al. | 5 Blind Source Separation: The Sparsity Revolution | |
CN108282424B (zh) | 用于四个数据集联合盲源分离的四阶张量联合对角化算法 | |
Zhang | Morphologically constrained ICA for extracting weak temporally correlated signals | |
JP5323860B2 (ja) | 混合された信号を複数の成分信号に分離する方法 | |
Staicu et al. | Significance tests for functional data with complex dependence structure | |
Soldati et al. | ICA analysis of fMRI with real-time constraints: an evaluation of fast detection performance as function of algorithms, parameters and a priori conditions | |
CN103530505B (zh) | 一种人脑语言认知建模方法 | |
Spurek et al. | ICA based on asymmetry | |
Hashemi et al. | Towards accelerated greedy sampling and reconstruction of bandlimited graph signals | |
Karhunen et al. | Finding dependent and independent components from related data sets: A generalized canonical correlation analysis based method | |
CN105760700B (zh) | 一种适于多被试复数fMRI数据分析的自适应定点IVA算法 | |
Li et al. | Flexible complex ICA of fMRI data | |
Boukouvalas | Development of ICA and IVA algorithms with application to medical image analysis | |
Fukumoto et al. | Node clustering of time-varying graphs based on temporal label smoothness | |
Malik et al. | FastICA based blind source separation for CT imaging under noise conditions | |
Weylandt et al. | Simultaneous grouping and denoising via sparse convex wavelet clustering | |
CN112399366A (zh) | 基于Hankel矩阵及WKNN方差提取的室内定位法 | |
Prasad | Independent component analysis | |
Gong et al. | A gridless method for direction finding with sparse arrays in nonuniform noise | |
Zhang et al. | Point cloud segmentation based on hypergraph spectral clustering | |
Kawanabe et al. | Joint low-rank approximation for extracting non-Gaussian subspaces |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |