发明名称 未知数目稀疏源盲分离的混叠矩阵估计方法
摘要 本发明是一种未知数目稀疏源盲分离的混叠矩阵估计方法。它属于工程领域,具体涉及到盲源分离技术领域。它的目的是解决目前基于经典聚类算法的稀疏源盲分离的混叠矩阵估计方法都要求源信号数目已知,并且估计精度较差的问题。本发明根据稀疏源混叠信号呈线性聚类的几何特点,基于聚类中心与每类数据致密点的距离关系,提出了一个新的聚类有效性准则,并根据此准则估计出源信号数目。同时利用霍夫变换寻找每一类数据的致密点,以代替聚类中心来估计混叠矩阵,提高了混叠矩阵的估计精度。本发明适用于源信号数目未知情况下的稀疏源盲分离的混叠矩阵的估计,提高了混叠矩阵的估计精度,广泛适用于语音识别、医学信号处理、无线通讯等领域。
申请公布号 CN101339620B 申请公布日期 2011.11.02
申请号 CN200810136921.5 申请日期 2008.08.13
申请人 哈尔滨工业大学 发明人 彭喜元;付宁;乔立岩;彭宇
分类号 G06K9/62(2006.01)I;G06K9/00(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 王吉东
主权项 1.未知数目稀疏源盲分离的混叠矩阵估计方法,其特征在于它的具体过程为:步骤A、将观测信号中的观测数据点中绝对值小于预定阈值的点去掉;所述观测信号在时域就具有较好的稀疏性;步骤B、将剩余的观测数据点投影到上半单位超球面上,获得所述观测数据点投影到上半单位超球面之后的数据集为Z={z<sub>j</sub>|j=1,2,…,L},其中z<sub>j</sub>是M维向量,L是数据点数;步骤C、估计源信号数目,并获得混叠矩阵,具体过程为:步骤C1、源信号数目为c,其取值的上限为c<sub>max</sub>,选取源信号数目c的初始值为2;步骤C2、以源信号数目c为聚类簇数目,对观测数据点运行20次K-means聚类算法,并取最优划分,获得第i类数据的聚类中心m<sub>i</sub>、第i类数据的第j个数据<img file="FSB00000496613000011.GIF" wi="78" he="98" />和第i类数据的点数n<sub>i</sub>;步骤C3、用霍夫变换求第i类数据的致密点p<sub>i</sub>;步骤C4、根据公式<maths num="0001"><![CDATA[<math><mrow><mi>FP</mi><mo>-</mo><mi>index</mi><mrow><mo>(</mo><mi>c</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>c</mi></munderover><mfrac><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>p</mi><mi>i</mi></msub><mo>-</mo><msub><mi>m</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mrow><mfrac><mn>1</mn><msub><mi>n</mi><mi>i</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>i</mi></msub></munderover><msup><mrow><mo>|</mo><mo>|</mo><msubsup><mi>z</mi><mi>j</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></msubsup><mo>-</mo><msub><mi>m</mi><mi>i</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>c</mi><mrow><mi>r</mi><mo>-</mo><mn>1</mn></mrow></msup></mrow></math>]]></maths>计算聚类有效性准则的函数值FP-index(c),其中,r称作惩罚因子;步骤C5、如果所述源信号数目c<c<sub>max</sub>,则令c=c+1,返回执行步骤C2;否则,执行步骤C6;步骤C6、比较源信号数目c从2至c<sub>max</sub>所对应的每个函数值FP-index获得最小函数值所对应源信号数目c,将其作为源信号数目的估计,然后将其对应的每类数据致密点按列向量组合起来作为混叠矩阵的估计;在上述过程中所述的i=1,2,…,c,所述j=1,2,…,n<sub>i</sub>。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号