发明名称 高精度延迟滤波器的多级插值设计方法
摘要 本发明涉及卫星导航信号模拟器的精密延迟信号产生技术。针对高倍率数字延迟滤波器的难以实现的问题,提出利用多个滤波器相互插值的设计方法得到精密时延滤波器。首先,根据滤波器的需求确定插值因子,分配滤波器参数和设计滤波器,再将滤波器进行相互插值得到所需设计的滤波器。
申请公布号 CN101915931A 申请公布日期 2010.12.15
申请号 CN201010221567.3 申请日期 2010.07.09
申请人 中国人民解放军国防科学技术大学 发明人 王跃科;杨俊;陈建云;张传胜;单庆晓;明德祥;钟小鹏
分类号 G01S19/23(2010.01)I 主分类号 G01S19/23(2010.01)I
代理机构 湖南省国防科学技术工业办公室专利中心 43102 代理人 冯青
主权项 1.高精度延迟滤波器的多级插值设计方法,其特征在于:采用多个滤波器相互插值的设计方法得到精密时延滤波器,该方法包括以下几个步骤:(1)确定两级插值因子对于插值倍数M<sub>0</sub>很大的情况,可以选取两级插值因子M<sub>1</sub>、M<sub>2</sub>,满足M<sub>0</sub>=M<sub>1</sub>M<sub>2</sub>。M<sub>1</sub>、M<sub>2</sub>的选取方法满足互质关系,大小相差不大;(2)分配滤波器参数根据两级插值因子分配情况和最终要求的超长阶滤波器特性确定,设计M<sub>1</sub>倍插值滤波器<img file="FSA00000179036000011.GIF" wi="138" he="66" />和M<sub>2</sub>倍插值滤波器<img file="FSA00000179036000012.GIF" wi="153" he="66" />的参数,设数字基带信号对应的采样频率为f<sub>s</sub>,信号带宽为B,两级插值因子分别为M<sub>1</sub>、M<sub>2</sub>,则M<sub>1</sub>倍插值滤波器设计参数设计如下:通带频率和阻带频率为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>f</mi><mrow><mi>PASS</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mi>B</mi><mrow><msub><mi>M</mi><mn>1</mn></msub><msub><mi>f</mi><mi>s</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>f</mi><mrow><mi>STOP</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>-</mo><mi>B</mi></mrow><mrow><msub><mi>M</mi><mn>1</mn></msub><msub><mi>f</mi><mi>s</mi></msub></mrow></mfrac></mrow></math>]]></maths>过渡带宽BW<sub>1</sub>为:BW<sub>1</sub>=f<sub>STOP1</sub>-f<sub>PASS1</sub>M<sub>2</sub>倍插值滤波器设计参数设计如下:通带频率和阻带频率为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>f</mi><mrow><mi>PASS</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mi>B</mi><mrow><msub><mi>M</mi><mn>2</mn></msub><msub><mi>M</mi><mn>1</mn></msub><msub><mi>f</mi><mi>s</mi></msub></mrow></mfrac></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>f</mi><mrow><mi>STOP</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>M</mi><mn>1</mn></msub><msub><mi>f</mi><mi>S</mi></msub><mo>-</mo><mi>B</mi></mrow><mrow><msub><mi>M</mi><mn>1</mn></msub><msub><mi>M</mi><mn>2</mn></msub><msub><mi>f</mi><mi>s</mi></msub></mrow></mfrac></mrow></math>]]></maths>过渡带宽BW<sub>2</sub>为:BW<sub>2</sub>=f<sub>STOP2</sub>-f<sub>PASS2</sub>(3)设计低倍率插值滤波器利用MATLAB滤波器设计工具,根据第二步得到的通带、阻带频率和过渡带宽,得到M<sub>1</sub>倍插值滤波器<img file="FSA00000179036000017.GIF" wi="138" he="65" />和M<sub>2</sub>倍插值滤波器<img file="FSA00000179036000018.GIF" wi="164" he="65" />(4)多级插值设计将<img file="FSA00000179036000021.GIF" wi="137" he="66" />看作信号,对其进行M<sub>2</sub>倍插值,插值的方式为在<img file="FSA00000179036000022.GIF" wi="138" he="67" />相邻两个数之间插入M<sub>2</sub>个零,插值后滤波器1变为h′<sub>M1</sub>(n),h′<sub>M1</sub>(n)=h<sub>M1</sub>(1),0,...0,h<sub>M1</sub>(2),0,...,0,h<sub>M1</sub>(3),......,h<sub>M1</sub>(M<sub>1</sub>)用h<sub>M2</sub>(n)对h′<sub>M1</sub>(n)进行滤波,得到最终的高阶滤波器h<sub>M0</sub>(n)h<sub>M0</sub>(n)=h′<sub>M1</sub>(n)×h<sub>M2</sub>(n),*表示卷积运算由于<img file="FSA00000179036000023.GIF" wi="137" he="66" />和<img file="FSA00000179036000024.GIF" wi="140" he="66" />均为有限时长序列,<img file="FSA00000179036000025.GIF" wi="139" he="66" />也必是有限时长序列,其长度为:<maths num="0005"><![CDATA[<math><mrow><msub><mi>N</mi><msub><mi>M</mi><mn>0</mn></msub></msub><mo>=</mo><msub><mi>M</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>N</mi><msub><mi>M</mi><mn>1</mn></msub></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>N</mi><msub><mi>M</mi><mn>2</mn></msub></msub></mrow></math>]]></maths>其中<img file="FSA00000179036000027.GIF" wi="226" he="61" />分别是两个滤波器的长度,上述设计结果,就是所要的M<sub>0</sub>倍插值滤波器。
地址 410073 湖南省长沙市砚瓦池正街47号