主权项 |
1.一种面向压缩感知的块稀疏信号重构方法,其特征是:所述方法的过程为:步骤一、采集块稀疏信号<b>x</b>的观测信号为<b>y,</b>观测信号<b>y</b>是长度为<i>m</i>的实数向量,即<img file="2010102406007100001DEST_PATH_IMAGE001.GIF" wi="78" he="34" />,设定块稀疏信号重构过程中各参数的初始状态值:其中,块稀疏信号<b>x</b>是长度为<i>N</i>、块稀疏度为<i>K</i>的实数向量,即<img file="2010102406007100001DEST_PATH_IMAGE002.GIF" wi="81" he="29" />,设定测量矩阵为<img file="2010102406007100001DEST_PATH_IMAGE003.GIF" wi="21" he="20" />是<i>m</i>行<i>N</i>列的实数矩阵,即<img file="2010102406007100001DEST_PATH_IMAGE004.GIF" wi="83" he="26" />,预先设定迭代误差err,分块向量<b>Group</b>的形式为:<img file="2010102406007100001DEST_PATH_IMAGE005.GIF" wi="277" he="48" />,其中,<i>N</i>=<i>M</i>×<i>d</i>,<i>M</i>为分块向量<b>Group</b>的分组数,<i>d</i>为分块向量<b>Group</b>的子块长度,设定残差的初始值<b>r</b><sub><b>0</b></sub>=<b>y</b>,恢复矩阵的初始值<img file="2010102406007100001DEST_PATH_IMAGE006.GIF" wi="57" he="26" />,迭代次数<i>l</i>的初始值为1,块稀疏信号<b>x</b>的重构向量<img file="2010102406007100001DEST_PATH_IMAGE007.GIF" wi="42" he="33" />;步骤二、根据第<i>l</i>-1次迭代后的残差<b>r</b><sub><i>l</i>-1</sub>,计算第<i>l</i>次迭代后的与残差<b>r</b><sub><i>l</i>-1</sub>最匹配的子空间<b>i</b><sub><i>l</i></sub>:<img file="DEST_PATH_IMAGE008.GIF" wi="236" he="37" />其中,测量矩阵<img file="712444DEST_PATH_IMAGE003.GIF" wi="21" he="20" />的<i>N</i>个列向量根据子块长度<i>d</i>将测量矩阵<img file="691902DEST_PATH_IMAGE003.GIF" wi="21" he="20" />分成<i>M</i>块,<img file="DEST_PATH_IMAGE009.GIF" wi="48" he="27" />为<img file="774434DEST_PATH_IMAGE003.GIF" wi="21" he="20" />的第<i>p</i>块对应的列向量组成的矩阵,<img file="DEST_PATH_IMAGE010.GIF" wi="102" he="25" />,<img file="DEST_PATH_IMAGE011.GIF" wi="108" he="25" />,<b>i</b><sub><i>l</i></sub>的值对应分块向量<b>Group</b>中的分组号,所述分组号为1,2,…,<i>M</i>,即<img file="356594DEST_PATH_IMAGE009.GIF" wi="48" he="27" />的转置与第<i>l</i>-1次迭代后的残差<b>r</b><sub><i>l</i>-1</sub>进行相乘运算后,得到一个长度为<i>d</i>的向量,将所述的向量中的多个元素取绝对值,再对所述的多个绝对值求其平均值,从而获得<i>M</i>个绝对平均值,从所获得的<i>M</i>个绝对平均值中选择最大的<i>K</i>个值的标号赋值给所述与残差<b>r</b><sub><i>l</i>-1</sub>最匹配的子空间<b>i</b><sub><i>l</i></sub>;步骤三、将步骤二中获得的子空间<b>i</b><sub><i>l</i></sub>所对应的测量矩阵<img file="DEST_PATH_IMAGE012.GIF" wi="23" he="21" />中的列向量组成的矩阵<img file="DEST_PATH_IMAGE013.GIF" wi="33" he="33" />与第<i>l</i>-1次迭代的恢复矩阵<img file="DEST_PATH_IMAGE014.GIF" wi="31" he="27" />的并集赋值给第<i>l</i>次迭代的过渡矩阵<img file="DEST_PATH_IMAGE015.GIF" wi="19" he="29" />,即:<img file="DEST_PATH_IMAGE016.GIF" wi="96" he="30" />;步骤四、根据观测信号<b>y</b>及步骤三获得的第<i>l</i>次迭代的过渡矩阵<img file="862312DEST_PATH_IMAGE015.GIF" wi="19" he="29" />,计算第<i>l</i>次迭代后的与观测信号<b>y</b>最匹配的子空间<b>t</b><sub><i>l</i></sub>:<img file="DEST_PATH_IMAGE017.GIF" wi="226" he="48" />其中,<i>j</i>为过渡矩阵<img file="DEST_PATH_IMAGE018.GIF" wi="21" he="30" />的分块标号,矩阵<img file="DEST_PATH_IMAGE019.GIF" wi="32" he="33" />为矩阵<img file="DEST_PATH_IMAGE020.GIF" wi="30" he="29" />的伪逆矩阵,<img file="DEST_PATH_IMAGE021.GIF" wi="150" he="32" />,矩阵<img file="186239DEST_PATH_IMAGE020.GIF" wi="30" he="29" />为过渡矩阵<img file="390956DEST_PATH_IMAGE015.GIF" wi="19" he="29" />所对应的测量矩阵<img file="776806DEST_PATH_IMAGE012.GIF" wi="23" he="21" />的列向量组成的矩阵,<b>t</b><sub><i>l</i></sub>的值对应分块向量<b>Group</b>中的分组号,即矩阵<img file="570319DEST_PATH_IMAGE019.GIF" wi="32" he="33" />的每一块分别与观测信号<b>y</b>进行乘积运算,将所述的多个乘积值分别取绝对值,再分别对所述的每一块块内的多个绝对值求其平均值,从中选择最大的<i>K</i>个值的标号赋值给所述与观测信号<b>y</b>最匹配的子空间<b>t</b><sub><i>l</i></sub>;步骤五、将步骤四获得的子空间<b>t</b><sub><i>l</i></sub>所对应的测量矩阵<img file="501366DEST_PATH_IMAGE012.GIF" wi="23" he="21" />中的列向量组成的矩阵<img file="DEST_PATH_IMAGE022.GIF" wi="33" he="33" />赋值给恢复矩阵<img file="DEST_PATH_IMAGE023.GIF" wi="23" he="31" />:<img file="DEST_PATH_IMAGE024.GIF" wi="72" he="33" />;步骤六、根据步骤五中获得的恢复矩阵<img file="DEST_PATH_IMAGE025.GIF" wi="20" he="27" />,计算第<i>l</i>次迭代后的残差<b>r</b><sub><i>l</i></sub>:<img file="DEST_PATH_IMAGE026.GIF" wi="167" he="38" />;其中,<img file="DEST_PATH_IMAGE027.GIF" wi="38" he="38" />为<img file="DEST_PATH_IMAGE028.GIF" wi="40" he="36" />的伪逆矩阵,<img file="DEST_PATH_IMAGE029.GIF" wi="150" he="31" />;步骤七、判断步骤六中所述的第<i>l</i>次迭代后的残差<b>r</b><sub><i>l</i></sub>是否小于预先设定迭代误差err,判断结果为是,则执行步骤十,判断结果为否,则执行步骤八;步骤八、判断步骤七中所述的迭代次数<i>l</i>的取值是否大于分组数<i>M</i>,判断结果为是,则执行步骤十,判断结果为否,则执行步骤九;步骤九、将迭代次数<i>l</i>的值加1,返回步骤二;步骤十、根据观测信号<b>y</b>及步骤六中获得的矩阵<img file="DEST_PATH_IMAGE030.GIF" wi="36" he="36" />,计算块稀疏信号<b>x</b>的重构向量<img file="DEST_PATH_IMAGE031.GIF" wi="18" he="30" />为:<img file="DEST_PATH_IMAGE032.GIF" wi="81" he="38" />,完成块稀疏信号<b>x</b>的重构。 |