发明名称 一种射频直接带通采样数字接收系统中时钟抖动消除方法
摘要 本发明涉及一种射频直接带通采样数字接收系统中时钟抖动消除方法。所述方法通过采用射频直接带通采样数字接收系统的非均匀采样输出值,利用带通采样定理,数字下变频原理和射频直接带通采样数字接收系统的特点,实现射频直接带通采样数字接收系统的时钟抖动的消除。该发明方法可根据射频直接带通采样数字接收实际需求选择不同的级联级数以满足所要求的时钟抖动消除性能。本发明所述方法实际实现复杂度低,易于硬件实现,消除效果好,对射频直接带通采样数字接收系统输入信号类型无限制,非常适用于射频直接带通采样数字接收系统在现代通讯,雷达,高精密数字仪器中的应用。
申请公布号 CN105447218A 申请公布日期 2016.03.30
申请号 CN201410808895.1 申请日期 2014.12.22
申请人 北京大学深圳研究生院 发明人 张晓宇;邹月娴;徐祥俊
分类号 G06F17/50(2006.01)I;H03K5/01(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 一种射频直接带通采样数字接收系统中时钟抖动消除方法,其特征在于针对射频直接带通采样数字接收系统(Direct Radio Frequency Sampling Receive,DRFSR),存在的时钟电路抖动,实现射频直接带通采样数字接收系统的时钟抖动误差的估计和消除,具体原理表述如下,对于射频直接带通采样数字接收系统,不失一般性,设其输入信号为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>r</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>Re</mi><mo>{</mo><mrow><mo>(</mo><msub><mi>s</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>j</mi><msub><mi>s</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mi>t</mi></mrow></msup><mo>}</mo></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>s</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>cos</mi><mo>[</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mi>t</mi><mo>]</mo><mo>-</mo><msub><mi>s</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>sin</mi><mo>[</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mi>t</mi><mo>]</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000640411510000011.GIF" wi="1348" he="160" /></maths>其中,s<sub>I</sub>(t)+js<sub>Q</sub>(t)为接收的带通信号对应的基带等效信号,f<sub>c</sub>为载波频率。设ADC的采样抖动RMS(均方根)值为ζ<sub>n</sub>,射频直接带通采样数字接收系统采用带通采样,则ADC输出信号可表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>r</mi><mrow><mo>(</mo><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>Re</mi><mo>{</mo><mrow><mo>(</mo><msub><mi>s</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>j</mi><msub><mi>s</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></msup><mo>}</mo></mtd></mtr><mtr><mtd><mo>=</mo><msub><mi>s</mi><mi>I</mi></msub><mrow><mo>(</mo><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>cos</mi><mo>[</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>]</mo></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>s</mi><mi>Q</mi></msub><mrow><mo>(</mo><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>sin</mi><mo>[</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mi>n</mi><msup><msub><mi>T</mi><mi>s</mi></msub><mo>&prime;</mo></msup><mo>+</mo><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><msub><mi>&zeta;</mi><mi>n</mi></msub><mo>]</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000640411510000012.GIF" wi="1552" he="246" /></maths>其中,T<sub>s</sub>′为ADC采样周期,利用带通采样定理和数字下变频原理可得,采样信号变换到基带的近似表达式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>z</mi><mi>n</mi></msub><mo>&ap;</mo><mi>s</mi><mrow><mo>(</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><msub><mi>&zeta;</mi><mi>n</mi></msub></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000640411510000013.GIF" wi="1149" he="75" /></maths>其中,T<sub>s</sub>为基带信号采样周期,z<sub>n</sub>是射频直接带通采样数字接收系统经过带通采样和数字下变频后获得的基带信号(为含有抖动的非均匀离散信号),s(nT<sub>s</sub>)为理想的均匀基带离散信号,并由下式定义s(nT<sub>s</sub>)=s<sub>I</sub>(nT<sub>s</sub>)+js<sub>Q</sub>(nT<sub>s</sub>)                       (4)由(3)式可得,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>z</mi><mi>n</mi></msub><mo>&ap;</mo><mi>s</mi><mrow><mo>(</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>s</mi><mrow><mo>(</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mo>&infin;</mo></munderover><mfrac><mrow><msup><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mo>)</mo></mrow><mi>k</mi></msup><msup><msub><mi>&zeta;</mi><mi>n</mi></msub><mi>k</mi></msup></mrow><mrow><mi>k</mi><mo>!</mo></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000640411510000014.GIF" wi="1334" he="147" /></maths>由(5)式可得,时钟抖动消除由下式计算:s(nT<sub>s</sub>)=z<sub>n</sub>‑e(n)                       (6)e(n)为重构得到的误差信号,且e(n)由下式计算:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>s</mi><mrow><mo>(</mo><mi>n</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mo>&infin;</mo></munderover><mfrac><mrow><msup><mrow><mo>(</mo><mi>j</mi><mn>2</mn><mi>&pi;</mi><msub><mi>f</mi><mi>c</mi></msub><mo>)</mo></mrow><mi>k</mi></msup><msup><msub><mi>&zeta;</mi><mi>n</mi></msub><mi>k</mi></msup></mrow><mrow><mi>k</mi><mo>!</mo></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000640411510000015.GIF" wi="1278" he="162" /></maths>利用(6)式,即可求得公式(5);综上,本算法实现了对射频直接带通采样数字接收系统中时钟抖动的消除。
地址 518055 广东省深圳市南山区西丽大学城北大园区