发明名称 一种多通道压缩感知框架下的混合参数估计方法
摘要 一种多通道压缩感知框架下的混合参数估计方法,涉及多通道压缩感知技术领域,解决现有混合参数估计计算必须先完成重构混合信号,源信号重构效率低的问题。采集混合信号x<sub>i</sub>的压缩观测信号为y<sub>i</sub>,取非线性函数g(·),所述函数g(·)的输入为yW<sub>l</sub>,输出为Y,计算Y的熵,计算熵H(Y)的梯度<img file="DDA00002974544500011.GIF" wi="87" he="44" />沿着熵H(Y)的梯度方向更新反混合矩阵W<sub>l+1</sub>,使得熵H(Y)逐渐增大,所述更新反混合矩阵W的公式为:<img file="DDA00002974544500012.GIF" wi="429" he="62" />将迭代次数l的值加1,l=l+1,判断迭代次数l是否大于设定的总的迭代次数t,据经过t次迭代更新得到的反混合矩阵W<sub>t</sub>,计算混合矩阵A的估计值<img file="DDA00002974544500013.GIF" wi="72" he="71" />本发明可广泛应用于对混合参数估计的计算。
申请公布号 CN103219998B 申请公布日期 2016.01.20
申请号 CN201310101905.3 申请日期 2013.03.27
申请人 哈尔滨工业大学 发明人 付宁;徐红伟;乔立岩;于伟;殷聪如
分类号 H03M7/30(2006.01)I 主分类号 H03M7/30(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种多通道压缩感知框架下的混合参数估计方法,其特征在于它包括如下步骤:步骤一:采集混合信号x<sub>i</sub>的压缩观测信号为y<sub>i</sub>,1≤i≤m;其中,x<sub>i</sub>是第i个混合信号,m为混合信号的数目,混合信号x<sub>i</sub>的长度为N,观测信号y<sub>i</sub>的长度为M,即<img file="FDA000029745442000110.GIF" wi="435" he="70" />,且M<<N;设:反混合矩阵W是m行m列的实数阵列,即<img file="FDA000029745442000111.GIF" wi="250" he="88" />测量矩阵为Ф是M行N列的实数矩阵,即<img file="FDA000029745442000112.GIF" wi="232" he="70" />算法迭代次数l的初始值为1,总的迭代次数为L,反混合矩阵初始值为W<sub>0</sub>,更新步长为η;步骤二、在(0,1)内单调递增的函数中选取任意非线性函数g(·);步骤三、设定函数g(·)的输入为yW<sub>l‑1</sub>,输出为Y,即Y=g(yW<sub>l‑1</sub>),其中y=[y<sub>1</sub>,y<sub>2</sub>,...,y<sub>m</sub>],y<sub>i</sub>表示第i个混合信号的压缩观测值;步骤四、计算Y的熵<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><mi>H</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>+</mo><mi>E</mi><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mi>ln</mi><mi>g</mi><mo>&prime;</mo><mrow><mo>(</mo><mi>y</mi><msub><mi>W</mi><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>]</mo><mo>+</mo><mi>ln</mi><mo>|</mo><msub><mi>W</mi><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>|</mo><mo>,</mo></mrow>]]></math><img file="FDA00002974544200011.GIF" wi="902" he="133" /></maths>其中H(Y)表示Y的熵,H(y)表示混合信号观测值y的熵,函数ln是对数函数,g′表示g的一阶导数,<img file="FDA00002974544200012.GIF" wi="381" he="134" />表示ln g′(yW<sub>l‑1</sub>)的均值;步骤五、计算熵H(Y)的梯度<img file="FDA00002974544200013.GIF" wi="87" he="44" /><img file="FDA00002974544200014.GIF" wi="600" he="374" />其中<img file="FDA00002974544200015.GIF" wi="85" he="53" />是第l‑1次循环计算的反混合矩阵W<sub>l‑1</sub>的各个元素,<img file="FDA00002974544200016.GIF" wi="130" he="137" />是熵关于变量<img file="FDA00002974544200017.GIF" wi="87" he="53" />的偏导数,<img file="FDA00002974544200018.GIF" wi="65" he="41" />是熵的梯度;步骤六、沿着熵H(Y)的梯度方向更新反混合矩阵,使得熵H(Y)逐渐增大,所述更新反混合矩阵的公式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>W</mi><mi>l</mi></msub><mo>=</mo><msub><mi>W</mi><mrow><mi>l</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mi>&eta;</mi><mo>*</mo><mo>&dtri;</mo><mi>h</mi><mo>;</mo></mrow>]]></math><img file="FDA00002974544200019.GIF" wi="395" he="52" /></maths>步骤七、判断迭代次数l是否大于等于设定的总的迭代次数L,判断结果为是,则执行步骤八,判断结果为否,则将迭代次数l的值加1,l=l+1,返回步骤三;步骤八、根据经过L次迭代更新得到的反混合矩阵W<sub>L</sub>,计算混合参数<img file="FDA00002974544200021.GIF" wi="38" he="54" />为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><mi>A</mi><mo>^</mo></mover><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>W</mi><mi>L</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow>]]></math><img file="FDA00002974544200022.GIF" wi="195" he="69" /></maths>所述混合参数<img file="FDA00002974544200023.GIF" wi="578" he="215" />其中a<sub>ij</sub>表示产生混合信号x<sub>i</sub>时,第j个源信号的权重;上述混合信号x<sub>i</sub>是指第i个采集到的多路源信号的混合信号。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号