发明名称 基于信息熵的数字存储示波器垂直分辨率提高方法
摘要 本发明公开了一种基于信息熵的数字存储示波器垂直分辨率提高方法,对过采样的单次采集信号,采样最大熵原理估计各样本的概率和最大熵,根据样本概率计算样本的不确定度,进而得到置信区间,根据置信区间剔除样本中的粗大量化误差,根据信息量比率计算数据的融合权重系数,并根据融合权重系数对有效样本进行数据融合,将数据融合结果代替粗大量化误差的样本对样本数据进行重构,再采用基于平均的抽取算法,得到输出样本。本发明用于提高数字存储示波器的垂直分辨率,通过结合基于信息熵的数据融合和基于平均的抽取滤波算法,提高数字示波器的垂直分辨率,并提高输出信号的信噪比性能,并且不会影响示波器实际采样率下的模拟带宽。
申请公布号 CN104101751B 申请公布日期 2016.10.26
申请号 CN201410314235.8 申请日期 2014.07.03
申请人 电子科技大学 发明人 蒋俊;潘卉青;汪明;王翔辉
分类号 G01R13/02(2006.01)I 主分类号 G01R13/02(2006.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 一种基于信息熵的数字存储示波器垂直分辨率提高方法,其特征在于,包括以下步骤:S1:数字存储示波器的采集系统按照N倍于实际采样率的高采样率进行过采样,将采样样本以N个样本为一组,得到样本序列x<sub>1</sub>,x<sub>2</sub>,...x<sub>N</sub>;S2:从样本序列中剔除重复样本后得到N′个样本x<sub>1</sub>,x<sub>2</sub>,...x<sub>N'</sub>,采用最大离散熵估计得到去重后各个样本x<sub>d</sub>的概率p(x<sub>d</sub>),并得到最大离散熵H(x)<sub>max</sub>,d的取值范围为1≤d≤N′,其中离散熵的计算公式为:<maths num="0001"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>k</mi><munderover><mo>&Sigma;</mo><mrow><mi>d</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mo>&prime;</mo></msup></munderover><mi>p</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>d</mi></msub><mo>)</mo></mrow><mi>ln</mi><mi> </mi><mi>p</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>d</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000993344290000011.GIF" wi="563" he="131" /></maths>概率估计的约束条件为:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mo>&Sigma;</mo><mrow><mi>d</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>p</mi><mo>(</mo><msub><mi>x</mi><mi>d</mi></msub><mo>)</mo><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>p</mi><mo>(</mo><msub><mi>x</mi><mi>d</mi></msub><mo>)</mo><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mo>&Sigma;</mo><mrow><mi>d</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mo>&prime;</mo></msup></munderover><mi>p</mi><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>)</mo><msub><mi>g</mi><mi>m</mi></msub><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>)</mo><mo>=</mo><mi>E</mi><mo>(</mo><msub><mi>g</mi><mi>m</mi></msub><mo>)</mo><mo>,</mo><mo>(</mo><mi>m</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi><mo>)</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000993344290000012.GIF" wi="894" he="407" /></maths>其中,k为正常数,g<sub>m</sub>(x<sub>d</sub>)为各阶统计矩函数,m表示阶数,M表示最大阶数,E(g<sub>m</sub>)表示各阶统计矩的期望值;S3:计算N′个样本的测量不确定度u,计算公式为:<maths num="0003"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>d</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mo>&prime;</mo></msup></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>d</mi></msub><mo>-</mo><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mi>p</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>d</mi></msub><mo>)</mo></mrow></mrow></msqrt></mrow>]]></math><img file="FDA0000993344290000013.GIF" wi="478" he="143" /></maths>其中,<img file="FDA0000993344290000014.GIF" wi="37" he="39" />表示步骤S1中过采样得到的N个样本的平均值;得到样本的置信区间为<img file="FDA0000993344290000015.GIF" wi="283" he="59" />将N′个样本中在置信区间范围外的样本剔除,得到N″个样本;S4:计算N″个样本中各个样本的信息量比率ω<sub>v</sub>:<maths num="0004"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>v</mi></msub><mo>=</mo><mfrac><mrow><mi>I</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>v</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>H</mi><msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000993344290000016.GIF" wi="286" he="135" /></maths>其中,v的取值范围为v=1,2,...N″,I(x<sub>v</sub>)表示样本x<sub>v</sub>的自信息量,I(x<sub>v</sub>)=‑lnp(x<sub>v</sub>);对信息量比率归一化处理得到融合权重系数q<sub>v</sub>,计算公式为:<maths num="0005"><math><![CDATA[<mrow><msub><mi>q</mi><mi>v</mi></msub><mo>=</mo><mfrac><mrow><mn>1</mn><mo>/</mo><msub><mi>&omega;</mi><mi>v</mi></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>v</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup></munderover><mn>1</mn><mo>/</mo><msub><mi>&omega;</mi><mi>v</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000993344290000021.GIF" wi="267" he="191" /></maths>进行数据融合得到融合结果x<sub>f</sub>,计算公式为:<maths num="0006"><math><![CDATA[<mrow><msub><mi>x</mi><mi>f</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>v</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup></munderover><msub><mi>q</mi><mi>v</mi></msub><msub><mi>x</mi><mi>v</mi></msub></mrow>]]></math><img file="FDA0000993344290000022.GIF" wi="246" he="134" /></maths>S5:在原始N个样本x<sub>1</sub>,x<sub>2</sub>,...x<sub>N</sub>中,根据步骤S3得到的置信区间,将置信区间范围外的样本剔除,并加入同等数量的融合结果x<sub>f</sub>,得到N个样本x′<sub>1</sub>,x′<sub>2</sub>,...x′<sub>N</sub>;S6:将N个样本x′<sub>1</sub>,x′<sub>2</sub>,...x′<sub>N</sub>进行基于平均的抽取,得到平均结果<img file="FDA0000993344290000023.GIF" wi="277" he="135" />将平均结果A作为样本输出。
地址 611731 四川省成都市高新区(西区)西源大道2006号