发明名称 一种随机采样的模拟电路压缩传感测量和信号重构方法
摘要 一种随机采样的模拟电路压缩传感测量和信号重构方法属于电子系统测试和故障诊断领域。针对模拟电路输出响应中故障信号自身或在其正交空间里具有稀疏性分布特点,本发明根据电路拓扑结构选择测试节点,在分布式传感器测试网络下随机采样电路输出响应,将响应信号利用离散正交基在其变换域上进行稀疏化表示,稀疏信号在观测矩阵投影下完成压缩感知测量,当这些随机压缩采样点对信号重构的恢复率达到80%以上时,表明此时电路输出响应的压缩测量值有效,它们可构造特征集并用于模拟电路故障诊断。该方法解决了传统模拟信号采样占用大量硬件资源、信号重构计算量大等问题,利用随机采样的压缩感知测量方法提高了电子系统测试效率。
申请公布号 CN102435934B 申请公布日期 2013.08.21
申请号 CN201110280538.9 申请日期 2011.09.20
申请人 北京工业大学 发明人 袁海英;黎海涛;梅家平
分类号 G01R31/316(2006.01)I 主分类号 G01R31/316(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 刘萍
主权项 1.一种随机采样的模拟电路压缩传感测量和信号重构方法,其特征在于,包括以下步骤:1)构造待测模拟电路的分布式测试网络并获取响应信号选择待测模拟电路的测试节点并设置传感器获取响应信号,每个传感器对同一个节点的输出信号感知测量M次,传感器每次返回模拟电路的M个响应信号;模拟信号x(t)经过A/D随机采样离散化为数字信号x,N为离散化后的数字信号x即采样点个数,M&lt;N,x(t)∈R,R为实数,x∈R<sup>1×N</sup>;2)模拟信号x(t)的稀疏化表示模拟信号x(t)为一维时域连续信号,它通过离散正交基Ψ向变换域向量<img file="FDA00003263808100011.GIF" wi="48" he="95" />进行投影;<img file="FDA00003263808100012.GIF" wi="522" he="143" />,其中,<img file="FDA00003263808100013.GIF" wi="51" he="95" />中的非0元素为K个,0元素为N-K个,稀疏度K由模拟信号x(t)的傅立叶变换确定;<img file="FDA00003263808100014.GIF" wi="62" he="99" />为离散正交基展开系数;离散正交基Ψ取高斯白噪声矩阵,Ψ∈R<sup>N×N</sup>,Ψ<sup>H</sup>是Ψ的共轭转置矩阵,且ΨΨ<sup>H</sup>=Ψ<sup>H</sup>Ψ=I,I为单位矩阵;利用公式(1)将数字信号x由离散正交基Ψ<sup>H</sup>展开成级数和形式:<maths num="0001"><![CDATA[<math><mrow><mi>x</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>&Psi;</mi><mi>n</mi><mi>H</mi></msubsup><mover><msub><mi>y</mi><mi>n</mi></msub><mo>^</mo></mover><mo>=</mo><msubsup><mi>&Psi;</mi><mn>1</mn><mi>H</mi></msubsup><mover><msub><mi>y</mi><mn>1</mn></msub><mo>^</mo></mover><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msubsup><mi>&Psi;</mi><mi>N</mi><mi>H</mi></msubsup><mover><msub><mi>y</mi><mi>N</mi></msub><mo>^</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>数字信号x的稀疏化表示过程按如下步骤依次执行:离散正交基Ψ的寻找过程和矩阵中元素分布关系见公式(2)所示;其中旋转角度θ∈[0°,90°];空白处元素均为0;对角线除了两个元素为cosθ,其余元素均为1,两个cosθ之间元素1的个数从0个增加至N-2个;sinθ和-sinθ位于矩阵的反对角线上,sinθ与对角线出现的第一个cosθ同列,与对角线出现的第二个cosθ同行,-sinθ与对角线出现的第一个cosθ同行,与对角线出现的第二个cosθ同列;<img file="FDA00003263808100021.GIF" wi="1843" he="823" />首先取θ=0°,得到对角线元素为1的酉矩阵<img file="FDA00003263808100022.GIF" wi="673" he="366" />(1)当公式(1)展开式满足K/N&lt;10%时,表明寻找到具有K稀疏度的变换域向量<img file="FDA00003263808100023.GIF" wi="44" he="89" />将数字信号x在离散正交基Ψ上进行稀疏化表示,此时,模拟电路输出响应中故障信号自身具有分布稀疏性和可压缩特点,直接转到步骤3)处理;(2)当公式(1)展开式满足K/N&gt;10%时,表明还未寻找到具有K稀疏度的变换域向量<img file="FDA00003263808100024.GIF" wi="46" he="81" />将数字信号x在离散正交基Ψ上进行稀疏化表示,模拟电路故障响应信号在离散正交基Ψ<sup>H</sup>具有分布稀疏性特点;取θ=10°,重新寻找一个离散正交基Ψ,再按公式(1)利用离散正交基Ψ<sup>H</sup>将模拟信号x(t)转换为稀疏信号;分如下两种情况进行判定:①若φ<sub>i·</sub>满足K/N&lt;10%时,转到步骤3)处理;②若φ<sub>i·</sub>满足K/N&gt;10%时,在θ∈[0°,90°]范围内从0°到90°依次改变旋转角度θ,重新寻找一个新的离散正交基Ψ,继续对数字信号x进行稀疏化表示;直到φ<sub>i·</sub>满足K/N&lt;10%为止,再转到步骤3)处理;否则,不断重复寻找合适的离散正交基Ψ<sup>H</sup>来稀疏化表示数字信号x,当反复寻找离散正交基Ψ<sup>H</sup>稀疏化表示数字信号x的迭代次数超过10或者K/N&gt;50%,表明旋转角度θ增加至90°时,φ<sub>i·</sub>仍不能满足K/N&lt;10%,无法找到合适的离散正交基Ψ来稀疏化表示数字信号x,压缩感知测量过程结束;3)构造用于响应信号稀疏化表示的观测矩阵生成一个M行N列的观测矩阵<img file="FDA00003263808100031.GIF" wi="579" he="236" />,该观测矩阵为高斯白噪声测量矩阵或者贝努利测量矩阵;此时φ的每一行φ<sub>i·</sub>视为一个传感器的一次线性测量;每一列φ<sub>·j</sub>视为响应信号随机采样一次的压缩测量值;观测矩阵φ表示利用传感器完成M次线性测量,每个线性测量均包含N个响应信号的采样分量;4)数字信号x的压缩感知测量模拟信号x(t)经过离散化处理后得到数字信号x,通过公式(3)得到感知压缩信号y,矩阵φ与数字信号x相乘后拾取了模拟信号x(t)的部分特征信息;感知压缩信号y代表模拟信号x(t)在相应观测点上的离散采样值;<maths num="0002"><![CDATA[<math><mrow><mi>y</mi><mo>=</mo><mi>&phi;x</mi><mo>=</mo><msup><mi>&phi;&Psi;</mi><mi>H</mi></msup><mover><mi>y</mi><mo>^</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>进行信号重构的方法,步骤如下:利用M次线性测量和观测矩阵φ在概率意义上重构数字信号x或离散正交基Ψ下等价的变换域向量<img file="FDA00003263808100033.GIF" wi="42" he="94" />,在满足<img file="FDA00003263808100034.GIF" wi="325" he="96" />约束条件下,通过线性凸优化运算寻优目标函数如公式(4)所示:<maths num="0003"><![CDATA[<math><mrow><mi>min</mi><mo>|</mo><mo>|</mo><mover><mi>y</mi><mo>^</mo></mover><mo>|</mo><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<maths num="0004"><![CDATA[<math><mrow><mo>|</mo><mo>|</mo><mover><mi>y</mi><mo>^</mo></mover><mo>|</mo><mo>|</mo><mo>=</mo><msup><mrow><mo>(</mo><msup><mrow><mo>(</mo><mover><msub><mi>y</mi><mn>1</mn></msub><mo>^</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>+</mo><msup><mrow><mo>(</mo><mover><msub><mi>y</mi><mi>N</mi></msub><mo>^</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>;</mo></mrow></math>]]></maths>所述的模拟信号恢复重建过程,当变换域向量<img file="FDA00003263808100037.GIF" wi="37" he="94" />在观测矩阵φ和离散正交基Ψ<sup>H</sup>共同作用下,模拟信号恢复率达到80%以上时,这些离散正交基展开系数<img file="FDA00003263808100038.GIF" wi="61" he="102" />可视为随机采样模拟电路的压缩感知测量值。
地址 100124 北京市朝阳区平乐园100号