发明名称 适用于高斯分布信号A/D转换器最佳输入功率设定方法
摘要 本发明属于信号处理技术领域。为在一般量化误差和过载误差之间取得平衡,使A/D转换器量化输出信噪比最优。为达到上述目的,本发明采取的技术方案是,适用于高斯分布信号A/D转换器最佳输入功率设定方法,包括下列步骤:步骤一:根据表1,查出A/D转换器量化位宽为n时的z值;步骤二:根据公式<img file="dda0000154533540000011.GIF" wi="254" he="62" />和已知条件A/D转换器最大量化幅值V<sub>max</sub>得到高斯分布输入信号的方差<img file="dda0000154533540000012.GIF" wi="80" he="80" />步骤三:根据<img file="dda0000154533540000013.GIF" wi="212" he="55" />计算A/D转换器高斯分布输入信号的最佳功率P<sub>x</sub>。本发明主要应用于信号处理。
申请公布号 CN102647190A 申请公布日期 2012.08.22
申请号 CN201210116678.7 申请日期 2012.04.18
申请人 天津大学 发明人 侯永宏;刘桂华;汪清;王洪宽;律会丽
分类号 H03M3/00(2006.01)I 主分类号 H03M3/00(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 刘国威
主权项 1.一种适用于高斯分布信号A/D转换器最佳输入功率设定方法,其特征是,包括如下步骤:假设待A/D转换器量化的高斯分布信号为x(n),其均值为0,方差为<img file="FDA0000154533510000011.GIF" wi="71" he="53" />则其概率密度函数为<img file="FDA0000154533510000012.GIF" wi="722" he="157" />A/D转换器的量化位宽为n,A/D转换器最大量化幅值为V<sub>max</sub>:步骤一:设定输入的高斯分布信号的功率P<sub>x</sub>近似等于输入的高斯分布信号的方差<img file="FDA0000154533510000013.GIF" wi="71" he="54" />步骤二:通过A/D转换器将高斯分布信号量化为n位,则量化级数M=2<sup>n</sup>,A/D转换器量化后的高斯分布信号的最大值和最小值分别为x<sub>q,max</sub>=(2<sup>n-1</sup>-1)Δx<sup>q,min</sup>=-2<sup>n-1</sup>Δ其中Δ为量化步长,且A/D转换器最大量化幅值V<sub>max</sub>=|x<sub>q,max</sub>|≈|x<sub>q,min</sub>|,则<maths num="0001"><![CDATA[<math><mrow><mi>&Delta;</mi><mo>&ap;</mo><mfrac><msub><mrow><mn>2</mn><mi>x</mi></mrow><mrow><mi>q</mi><mo>,</mo><mi>max</mi></mrow></msub><mi>M</mi></mfrac><mo>=</mo><mfrac><msub><mi>x</mi><mrow><mi>q</mi><mo>,</mo><mi>max</mi></mrow></msub><msup><mn>2</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup></mfrac><mo>=</mo><mfrac><msub><mi>V</mi><mi>max</mi></msub><msup><mn>2</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup></mfrac><mo>,</mo></mrow></math>]]></maths>此时,可得到限幅即过载误差的功率:<maths num="0002"><![CDATA[<math><mrow><msub><mi>N</mi><mi>s</mi></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mi>xq</mi><mo>,</mo><mi>max</mi></mrow><mo>&infin;</mo></msubsup><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mrow><mi>q</mi><mo>,</mo><mi>max</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>dx</mi><mo>+</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><msub><mi>x</mi><mrow><mi>q</mi><mo>,</mo><mi>min</mi></mrow></msub></msubsup><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mrow><mi>q</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>dx</mi></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mo>&ap;</mo><mfrac><mn>2</mn><mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt><msub><mi>&sigma;</mi><mi>x</mi></msub></mrow></mfrac><msubsup><mo>&Integral;</mo><mrow><mi>V</mi><mi>max</mi></mrow><mo>&infin;</mo></msubsup><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>V</mi><mi>max</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mi>x</mi><mn>2</mn></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup></mrow></mfrac></mrow></msup><mi>dx</mi></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mo>=</mo><mi>erfc</mi><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>max</mi></msub><mrow><msqrt><mn>2</mn></msqrt><msub><mi>&sigma;</mi><mi>x</mi></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>x</mi></msub><mo>&CenterDot;</mo><msub><mi>V</mi><mi>max</mi></msub></mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>x</mi><mn>2</mn></msubsup></mfrac></mrow></msup></mrow></math>]]></maths>其中x表示A/D转换器量化前的高斯分布信号,<img file="FDA0000154533510000018.GIF" wi="522" he="139" />是互补误差函数;截尾误差或舍入误差的功率为:<maths num="0005"><![CDATA[<math><mrow><msub><mi>N</mi><mi>q</mi></msub><mo>&ap;</mo><mfrac><msup><mi>&Delta;</mi><mn>2</mn></msup><mn>12</mn></mfrac><mo>=</mo><mfrac><mn>1</mn><mn>12</mn></mfrac><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>max</mi></msub><msup><mn>2</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><msup><mrow><mn>3</mn><mo>&CenterDot;</mo><mn>4</mn></mrow><mi>n</mi></msup></mfrac></mrow></math>]]></maths>总的量化误差功率即为:<maths num="0006"><![CDATA[<math><mrow><mi>N</mi><mo>=</mo><msub><mi>N</mi><mi>q</mi></msub><mo>+</mo><msub><mi>N</mi><mi>S</mi></msub></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mo>=</mo><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><msup><mrow><mn>3</mn><mo>&CenterDot;</mo><mn>4</mn></mrow><mi>n</mi></msup></mfrac><mo>+</mo><mi>erfc</mi><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>max</mi></msub><mrow><msqrt><mn>2</mn></msqrt><msub><mi>&sigma;</mi><mi>x</mi></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>x</mi></msub><mo>&CenterDot;</mo><msub><mi>V</mi><mi>max</mi></msub></mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>x</mi><mn>2</mn></msubsup></mfrac></mrow></msup></mrow></math>]]></maths>则量化信噪比为<maths num="0008"><![CDATA[<math><mrow><mi>snr</mi><mo>=</mo><mn>10</mn><mi>lg</mi><mrow><mo>(</mo><mfrac><msub><mi>P</mi><mi>x</mi></msub><mi>N</mi></mfrac><mo>)</mo></mrow><mo>=</mo><mn>10</mn><mi>lg</mi><mrow><mo>(</mo><mfrac><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mrow><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><msup><mrow><mn>3</mn><mo>&CenterDot;</mo><mn>4</mn></mrow><mi>n</mi></msup></mfrac><mo>+</mo><mi>erfc</mi><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>max</mi></msub><mrow><msqrt><mn>2</mn></msqrt><msub><mi>&sigma;</mi><mi>x</mi></msub></mrow></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>x</mi></msub><mo>&CenterDot;</mo><msub><mi>V</mi><mi>max</mi></msub></mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><msubsup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>x</mi><mn>2</mn></msubsup></mfrac></mrow></msup></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>令<img file="FDA0000154533510000022.GIF" wi="156" he="122" />对量化信噪比关于<img file="FDA0000154533510000023.GIF" wi="57" he="61" />求导相当于对y关于<img file="FDA0000154533510000024.GIF" wi="54" he="61" />求导<maths num="0009"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>y</mi></mrow><mrow><mo>&PartialD;</mo><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup></mrow></mfrac><mo>=</mo><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><msup><mrow><mn>3</mn><mo>&CenterDot;</mo><mn>4</mn></mrow><mi>n</mi></msup></mfrac><mo>+</mo><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><mo>&CenterDot;</mo><mi>erfc</mi><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>max</mi></msub><mrow><msqrt><mn>2</mn></msqrt><msub><mi>&sigma;</mi><mi>x</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>x</mi></msub><mo>&CenterDot;</mo><msub><mi>V</mi><mi>max</mi></msub></mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msubsup><mi>V</mi><mi>max</mi><mn>2</mn></msubsup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup></mrow></mfrac></mrow></msup><mo>=</mo><mn>0</mn></mrow></math>]]></maths>令<img file="FDA0000154533510000026.GIF" wi="322" he="84" />则上式可简化为<maths num="0010"><![CDATA[<math><mrow><mfrac><mi>z</mi><msup><mrow><mn>3</mn><mo>&CenterDot;</mo><mn>4</mn></mrow><mi>n</mi></msup></mfrac><mo>+</mo><mi>z</mi><mo>&CenterDot;</mo><mi>erfc</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mn>1</mn><msqrt><mi>&pi;</mi></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><msup><mi>z</mi><mn>2</mn></msup></mrow></msup><mo>=</mo><mn>0</mn></mrow></math>]]></maths>不同A/D转换器量化位宽n对应下的z值如表1所示,表1.不同n值下的z值<tables num="0001"><table><tgroup cols="6"><colspec colname="c001" colwidth="17%" /><colspec colname="c002" colwidth="17%" /><colspec colname="c003" colwidth="17%" /><colspec colname="c004" colwidth="17%" /><colspec colname="c005" colwidth="17%" /><colspec colname="c006" colwidth="17%" /><tbody><row><entry morerows="1">n</entry><entry morerows="1">z</entry><entry morerows="1">n</entry><entry morerows="1">z</entry><entry morerows="1">n</entry><entry morerows="1">z</entry></row><row><entry morerows="1">1</entry><entry morerows="1">0.8767</entry><entry morerows="1">11</entry><entry morerows="1">3.3653</entry><entry morerows="1">21</entry><entry morerows="1">4.9105</entry></row><row><entry morerows="1">2</entry><entry morerows="1">1.2096</entry><entry morerows="1">12</entry><entry morerows="1">3.5449</entry><entry morerows="1">22</entry><entry morerows="1">5.0421</entry></row><row><entry morerows="1">3</entry><entry morerows="1">1.5214</entry><entry morerows="1">13</entry><entry morerows="1">3.7173</entry><entry morerows="1">23</entry><entry morerows="1">5.1707</entry></row><row><entry morerows="1">4</entry><entry morerows="1">1.8096</entry><entry morerows="1">14</entry><entry morerows="1">3.8834</entry><entry morerows="1">24</entry><entry morerows="1">5.2965</entry></row><row><entry morerows="1">5</entry><entry morerows="1">2.0762</entry><entry morerows="1">15</entry><entry morerows="1">4.0438</entry><entry morerows="1">25</entry><entry morerows="1">5.4196</entry></row><row><entry morerows="1">6</entry><entry morerows="1">2.3242</entry><entry morerows="1">16</entry><entry morerows="1">4.1990</entry><entry morerows="1">26</entry><entry morerows="1">5.5402</entry></row><row><entry morerows="1">7</entry><entry morerows="1">2.5563</entry><entry morerows="1">17</entry><entry morerows="1">4.3494</entry><entry morerows="1">27</entry><entry morerows="1">5.6585</entry></row><row><entry morerows="1">8</entry><entry morerows="1">2.7747</entry><entry morerows="1">18</entry><entry morerows="1">4.4954</entry><entry morerows="1">28</entry><entry morerows="1">5.7746</entry></row><row><entry morerows="1">9</entry><entry morerows="1">2.9814</entry><entry morerows="1">19</entry><entry morerows="1">4.6374</entry><entry morerows="1">29</entry><entry morerows="1">5.8886</entry></row><row><entry morerows="1">10</entry><entry morerows="1">3.1779</entry><entry morerows="1">20</entry><entry morerows="1">4.7757</entry><entry morerows="1">30</entry><entry morerows="1">6.0005</entry></row></tbody></tgroup></table></tables>步骤三:将得到的z值和已知条件A/D转换器最大量化幅值V<sub>max</sub>代入<img file="FDA0000154533510000028.GIF" wi="288" he="68" />即得到使A/D转换器量化信噪比最大的高斯分布输入信号的方差<img file="FDA0000154533510000029.GIF" wi="55" he="62" />值;步骤四:利用公式<img file="FDA00001545335100000210.GIF" wi="211" he="55" />对高斯分布信号方差<img file="FDA00001545335100000211.GIF" wi="56" he="61" />进行修正得到A/D转换器高斯分布输入信号的最佳功率的精确值P<sub>x</sub>,R为A/D转换器的输入阻抗。
地址 300072 天津市南开区卫津路92号