发明名称 基于SCMP算法的欠定盲分离源信号恢复方法
摘要 本发明公开了一种基于SCMP算法的欠定盲分离源信号恢复方法。主要解决现有的欠定盲分离技术中复杂度高和恢复精度低的问题,其实现步骤是:1.对混合矩阵的转置A<sup>T</sup>进行QR分解,得到矩阵空间S<sub>A</sub>及对应的正交补空间<img file="DDA0000628443540000011.GIF" wi="66" he="46" />2.利用接收信号x(t)求出源信号在空间S<sub>A</sub>上的分量s<sub>1</sub>;3.根据源信号在空间S<sub>A</sub>以及空间<img file="DDA0000628443540000012.GIF" wi="43" he="45" />上的互补性质,利用近似l<sub>0</sub>范数方法搜索源信号s中N-M个零元素所在的位置;4.根据搜索出的N-M个零元素的位置,利用互补特性求出其他M个元素,得到恢复信号。本发明搜索速度快,恢复精度高,可用于通信及信号处理等领域。
申请公布号 CN104392146A 申请公布日期 2015.03.04
申请号 CN201410751676.4 申请日期 2014.12.09
申请人 西安电子科技大学;航天恒星科技有限公司 发明人 付卫红;陈杰虎;刘乃安;胡梅霞;黑永强;熊超;杨博
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于SCMP算法的欠定盲分离源信号恢复方法,包括如下步骤:(1)通过传感器接收待分离的混合信号x(t),利用混合矩阵估计方法估计出混合矩阵A,A是一个M×N维的混合矩阵;(2)对混合矩阵的转置A<sup>T</sup>进行QR分解:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>A</mi><mi>T</mi></msup><mo>=</mo><mi>QR</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>Q</mi><mn>1</mn></msub></mtd><mtd><msub><mi>Q</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>R</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000628443510000011.GIF" wi="529" he="160" /></maths>其中Q为列向量相互正交的N维方阵,Q<sub>1</sub>为方阵Q的前M列组成的子矩阵,Q<sub>2</sub>为方阵Q的后N‑M列组成的子矩阵,R<sub>1</sub>是一个M×M维的上三角矩阵;(3)利用接收到的混合信号x(t)求出源信号在矩阵空间R(A<sup>T</sup>)上的分量s<sub>1</sub>=A<sup>T</sup>(AA<sup>T</sup>)<sup>‑1</sup>x(t);(4)利用混合矩阵A的所有N个列对应的下标1,2,...N组成下标全集U={1,2,...N},从下标全集U={1,2,...N}中找出任意N‑M个元素构成下标子集Φ<sub>i</sub>,<img file="FDA0000628443510000012.GIF" wi="355" he="78" />对每个下标子集Φ<sub>i</sub>,求出误差向量e<sub>i</sub>=s<sub>1</sub>‑Q<sub>2</sub>(Q<sub>2</sub>[Φ<sub>i</sub>])<sup>‑1</sup>s<sub>1</sub>[Φ<sub>i</sub>]以及对应的误差函数<img file="FDA0000628443510000013.GIF" wi="473" he="150" />其中Q<sub>2</sub>[Φ<sub>i</sub>]表示Q<sub>2</sub>中下标子集Φ<sub>i</sub>所对应的N‑M个行向量组成的矩阵,s<sub>1</sub>[Φ<sub>i</sub>]表示下标子集Φ<sub>i</sub>所对应的N‑M个元素所组成的向量,e<sub>k</sub>表示误差向量e<sub>i</sub>中的第k个元素,δ为一小于0.5的常数;(5)从步骤(4)中求得的<img file="FDA0000628443510000014.GIF" wi="120" he="84" />个误差函数f<sub>i</sub>中,找出最小值对应的指标j以及对应的下标子集<img file="FDA0000628443510000015.GIF" wi="70" he="76" />即:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>j</mi><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mi>i</mi></munder><msub><mi>f</mi><mi>i</mi></msub><mo>,</mo><mover><mi>&Phi;</mi><mo>~</mo></mover><mo>=</mo><msub><mi>&Phi;</mi><mi>j</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000628443510000016.GIF" wi="497" he="111" /></maths>(6)根据步骤(5)中求得的下标子集<img file="FDA0000628443510000017.GIF" wi="98" he="71" />得到恢复信号<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>s</mi><mo>^</mo></mover><mo>=</mo><msub><mi>s</mi><mn>1</mn></msub><mo>-</mo><msub><mi>Q</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><msub><mi>Q</mi><mn>2</mn></msub><mo>[</mo><mover><mi>&Phi;</mi><mo>~</mo></mover><mo>]</mo><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>s</mi><mn>1</mn></msub><mo>[</mo><mover><mi>&Phi;</mi><mo>~</mo></mover><mo>]</mo><mo>.</mo></mrow>]]></math><img file="FDA0000628443510000018.GIF" wi="555" he="86" /></maths>
地址 710071 陕西省西安市太白南路2号