发明名称 一种基于粒子群优化的信号分离方法
摘要 本发明公开了一种基于粒子群优化的信号分离方法,包括:(1)输入观测信号;(2)以最小化分离信号联合概率与边缘概率简单乘积间的差值为优化目标,建立独立分量分析算法优化模型;(3)根据奇异值分解法估计源信号个数,依此确定优化变量个数;(4)计算观测信号间的相关系数,确定优化变量取值范围;(5)利用粒子群优化算法对分离矩阵进行优化;(6)优化结束后最后一代种群中适应度最优的粒子为最优分离矩阵,将其与混合信号相乘,即可得到最优分离信号。该方法基于粒子群优化的独立分量分析,具有普遍适用性,对于各种盲源分离问题均具有良好的性能。
申请公布号 CN103427791B 申请公布日期 2016.09.28
申请号 CN201310344881.4 申请日期 2013.08.08
申请人 长安大学 发明人 李良敏;高强;房宏威;冯帆;杨本波;周劲草
分类号 H03H21/00(2006.01)I 主分类号 H03H21/00(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 汪人和
主权项 一种基于粒子群优化的信号分离方法,其特征在于,包括以下步骤:1)读取观测信号;2)以最小化分离信号联合概率与边缘概率简单乘积间的差值为优化目标,建立如下所示的独立分量分析算法优化模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>M</mi><mi>i</mi><mi>n</mi><mi>i</mi><mi>m</mi><mi>i</mi><mi>z</mi><mi>e</mi><mo>|</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>,</mo><mn>...</mn><mo>,</mo><msub><mi>y</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>-</mo><munderover><mo>&Pi;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mi>P</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>r</mi><mo>.</mo><mi>t</mi><mo>.</mo><mi>W</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000991214180000011.GIF" wi="813" he="119" /></maths>s.t.:y(t)=Wx(t)其中:x(t)为n维观测信号,y(t)为m维分离信号,W为分离矩阵,P(y<sub>1</sub>,y<sub>2</sub>,...,y<sub>m</sub>)为分离信号y<sub>1</sub>,y<sub>2</sub>,...,y<sub>m</sub>的联合概率,P(y<sub>i</sub>)为分离信号y<sub>i</sub>的边缘概率;3)根据奇异值分解法确定源信号个数,依此确定分离矩阵维数;4)计算观测信号间的相关系数γ<sub>ij</sub>,确定分离矩阵取值范围;分离矩阵W取值范围确定方法如下:4.1)根据下面的公式计算各观测信号间的相关系数γ<sub>ij</sub>:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>X</mi><mi>j</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>j</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>-</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow></mfrac></mrow>]]></math><img file="FDA0000991214180000012.GIF" wi="854" he="278" /></maths>其中X<sub>i</sub>(k),k=1,...,N为呈离散状态的第i组观测信号,N为观测信号长度,<img file="FDA0000991214180000013.GIF" wi="54" he="69" />为X<sub>i</sub>(k),k=1,...,N的均值;4.2)若min(γ<sub>ij</sub>)&gt;0.9,分离矩阵W取值范围[‑100,100];若0.9&gt;min(γ<sub>ij</sub>)&gt;0.7,分离矩阵W取值范围[‑10,10];若0.7&gt;min(γ<sub>ij</sub>)&gt;0.3,分离矩阵W取值范围[‑1,1];5)利用粒子群优化算法求解最优分离矩阵;6)优化结束后最后一代种群中性能最优的粒子为最优分离矩阵,将其与混合信号相乘,即可得到最优分离信号。
地址 710064 陕西省西安市雁塔区南二环中段33号