主权项 |
一种基于粒子群优化的信号分离方法,其特征在于,包括以下步骤: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>Π</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>γ</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><munderover><mo>Σ</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>‾</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>‾</mo></mover><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><msqrt><mrow><munderover><mo>Σ</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>‾</mo></mover><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><msqrt><mrow><munderover><mo>Σ</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>‾</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>)>0.9,分离矩阵W取值范围[‑100,100];若0.9>min(γ<sub>ij</sub>)>0.7,分离矩阵W取值范围[‑10,10];若0.7>min(γ<sub>ij</sub>)>0.3,分离矩阵W取值范围[‑1,1];5)利用粒子群优化算法求解最优分离矩阵;6)优化结束后最后一代种群中性能最优的粒子为最优分离矩阵,将其与混合信号相乘,即可得到最优分离信号。 |