FFT捕获算法部分匹配滤波器长度选取研究

来源 :航空兵器 | 被引量 : 0次 | 上传用户:tofomy
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  摘要: 高动态环境下具有较大的多普勒频移, 不利于信号的快速捕获。 基于PMFFFT的快速捕获算法利用部分匹配滤波器结合FFT运算的方法可实现扩频信号的快速捕获, 大大减少捕获时间, 具有广泛的应用价值。 本文通过建立数学模型, 详细分析了PMFFFT算法的特性, 确定了部分匹配滤波器长度的影响因素和设计准则, 为算法的实现解决了首要问题。
  关键词: PMFFFT; 部分匹配滤波器; 高动态; 捕获性能
  中图分类号: TN914.42文献标识码: A文章编号: 1673-5048(2015)05-0046-06
  Abstract: It is not easy for signal fast acquisition in high dynamic environment with large Doppler frequency. The fast acquisition algorithm based PMFFFT makes partial matched filter combine with FFT to realize the fast acquisition of the spread spectrum signal, which greatly reduces the acquisition time and has the widespread application value. By building mathematical model, the characteristics of PMFFFT algorithm are analyzed in detail, and the influencing factors and design criteria of length selection for partial matched filter are confirmed, which solves the primary problem for the realization of algorithm.
  Key words: PMFFFT; partial matched filter; high dynamic; acquisition performance
  0引言
  卫星通信以其覆盖广、 全天候、 频带资源丰富等众多优势成为非常重要的通信手段。 许多卫星系统已经得到广泛应用, 例如GPS全球卫星定位系统, Gb(全球星)系统, 铱星(Iridium)系统都是非常成功的卫星通信系统。 卫星通信系统在传输过程中由于信号环境复杂恶劣, 不可避免会有较大的衰减导致信噪比较低; 而且在高动态环境下, 卫星与载体间的传输链路会存在较高的多普勒频移。 目前卫星通信中所采用的扩频通信技术以其抗干扰能力强、 可靠性高而广为认可。 扩频通信中最为重要的环节之一就是扩频信号的同步过程, 分为捕获和跟踪两个过程。 高动态环境下由于多普勒频移较大, 导致接收机相关运算产生较大的衰减, 并且为保证通信的可靠性, 卫星通信系统都使用了长度较长的伪随机码, 这些都不利于快速捕获。 卫星通信要求有较高的实时性, 因此研究如何能够在较大多普勒频移下实现长度较大的伪码快速捕获算法具有非常重要的意义。
  20世纪以来扩频通信技术不断发展, 扩频信号的快速捕获问题一直是研究的热点。 目前较为成熟的捕获算法主要有滑动相关捕获法、 匹配滤波器, 以及基于FFT的捕获算法等。 20世纪80年代, 文献[1]中首次将匹配滤波器引入到扩频捕获算法中。 20世纪90年代开始, 将FFT引入捕获算法成为新的热点。 文献[2]首次提出将匹配滤波器与FFT结合起来实现快速捕获的算法, 该算法在部分匹配滤波器的基础上加入了FFT运算以改良匹配滤波器在多普勒频移状况下的缺点, 使得系统在较大多普勒频移的条件下依然具有较高的检测概率和优秀的捕获时间。 文献[3]描述了大多普勒频偏下的一种非相干检测的快速捕获系统, 首次给出PMFFFT捕获模型, 并分析了高斯信道下的虚警概率和检测概率特性; 然后利用该模型对接收伪码进行匹配滤波处理, 进行FFT补偿, 可以增加频偏的搜索范围。 文献[4]中, Spangenberg等人详细分析了 PMFFFT 捕获系统的数学模型以及性能, 除了补零并加长 FFT 点数的处理方法外, 还提出对 FFT输入数据进行加窗处理来降低扇贝损失。 之后的研究大多是在此基础上如何改进算法的实现方法, 而鲜有实质的理论创新。 部分匹配滤波器的长度选择影响着整个算法的大多普勒频移下的捕获能力, 对于整个捕获算法的系统性能有着至关重要的影响, 但是对于整个系统的具体影响因素以及如何根据用户指标选取滤波器长度, 之前的研究鲜有给出综合的结论。
  本文研究了基于PMFFFT的快速捕获算法的原理性能, 全面考虑了匹配滤波器的长度选择对算法理论性能和捕获性能的影响, 可以解决捕获系统设计中匹配滤波器的长度选择问题, 从而为该算法在实际系统中的应用奠定坚实基础。
  1PMFFFT算法基本原理及数学模型
  PMFFFT算法是基于匹配滤波器算法与FFT算法相结合的捕获算法模型, 其结构示意图如图1所示[4]。
  该系统包含了P个部分匹配滤波器, 每个长度为X, XP=M; M为伪随机序列长度。 部分匹配滤波器的主要作用是对接收信号进行部分解扩, 得到一定的扩频增益, 从而提高信噪比。 第一个处理伪随机序列的前X个码片, 第二个处理后X个, 以此类推。 这P个部分匹配滤波器的输出为部分相关输出。 如果把这些输出进行相干累加就会与长度为M的标准匹配滤波器的响应相同。 对这些输出结果进行N点FFT运算, N≥P, 再进入一个最大信号选择器来进行多普勒估计。 FFT输出端的最大信号选择器的输出与多普勒频移最接近频率相同。 该方法将搜索伪码相位和估计多普勒频移同时进行, 即将原来的相位、 频率二维搜索过程变为伪码相位的一维搜索, 大大减少了捕获时间[5]。   航空兵器2015年第5期黄云青等: PMFFFT捕获算法部分匹配滤波器长度选取研究对其工作的数学原理进行分析, 假设接收的序列与本地伪随机序列已经同步, 设接收到的信号经过下变频和采样之后的量化输出信号为[6]
  式中:G1(fd)为PMF部分的频率响应, 是由部分相关累积引起的, 部分匹配滤波实质上是一个低通滤波过程, 如果令M=P, 则G1(fd)=1, 相当于除去了部分匹配滤波过程; G2(fd,k)是FFT部分的频率响应, 为FFT运算对归一化幅频响应的贡献。 由此也可看出PMFFFT算法是匹配滤波器算法和FFT算法的有机结合, 通过FFT方法来补偿匹配滤波器算法的响应。
  2PMFFFT算法的理论性能分析
  以GPS C/A码为例, M=1 023, 1/Tc=1.023 Mb/s。 令N=P=8, 此时X=128。 伪码同步时的归一化相关增益如图2所示。
  可以看出, 整个曲线的包络呈下降趋势, 这是由部分匹配滤波过程所引起的随多普勒频率的相关损失。 当多普勒频率落在两个FFT输出点之间时, 存在着明显的扇贝损失, 且当fd=12 MTc时扇贝损失最大[8], 最大扇贝损失L为
  L=Gfd,k=1Psin(π2)sin(π2P)≈2π=-1.9 dB(8)
  式(8)表明相关增益较多普勒频移为0时有最大约为1.9 dB的衰减。 因此门限设置的上限就会降低, 从而导致信号虚警概率的上升。
  2.1多普勒频率估计分辨率及范围
  当πfdTcM=πPNk时, G(fd,k)取最大值[9]。 通过设置合适的门限, 若N个FFT的输出值有超出门限的, 则可判断伪码已经捕获, 可估计多普勒频偏, 即fd=kNXTc。 FFT算法的频谱估计精度为Δf=1NXTc。
  由此可见, 该方法的捕获精度由分段匹配滤波器的长度与FFT的点数共同决定。 匹配滤波器长度越大, 频率估计精度越高,
  在伪码取得同步时, 部分匹配滤波器的响应G1(fd)相当于一个长度为X的矩形窗, 如图3所示, 它的主瓣宽度可表示为B=1XTc。
  PMFFFT系统幅频响应函数的第一个过零点对应的也是G1(fd)的第一个过零点, 即理论的单边分析带宽为fd=1XTc, 因此理论的捕获频率范围为fd≤1XTc。 为保证信号能量不会有太大损失, 一般取到1/4处作为实际分析带宽范围, 即fd≤14XTc。 可以看出匹配滤波器的长度越小, 分析带宽越大, 捕获的频率范围也就更大。
  2.2相关损失
  相关损失由部分匹配过程引起, 幅频响应包络可用G1(fd)表示。 则相关损失可表示为[10]
  K=20lg(G1)=20lgsin(πXfdTc)Xsin(πfdTc)(10)
  在码型确定的情况下, 相关损失取决于多普勒频移的大小和匹配滤波器长度。 C/A码在不同匹配滤波器长度下的相关损失如图4所示。
  匹配滤波器长度越小, 相关损失也越小, 但此时匹配滤波器规模与实现复杂度也就越大, 可根据需要来选择合适的匹配滤波器长度。
  例如, 若系统要求可接受的最大相关损失为-1 dB, 根据相关损失与匹配滤波器长度的曲线图可以判断, 滤波器长度为32和16时满足要求, 为降低硬件实现复杂度, 可选择匹配滤波器长度为32。 但如果已知最大多普勒频移达到10 kHz, 根据X≤14fdTc, 可知要满足X≤25.6, 此时匹配滤波器长度只能选择16。
  2.3扇贝损失
  扇贝损失是由FFT部分相位补偿不完全所造成的, 最大扇贝损失为
  由式(11)可知, 在码型确定的情况下, 最大扇贝损失取决于匹配滤波器长度X和FFT的点数N。 假定FFT点数为匹配滤波器的数目(不补零), 最大扇贝损失随着匹配滤波器长度的增大而减小, 但变化幅度较小, 如图5所示。
  3捕获性能分析
  实际工程中, 通常用检测概率、 虚警概率和平均捕获时间这三个指标来衡量一个捕获系统的工作性能[11]。 因此在选取匹配滤波器的长度时, 需要考虑它对这三个指标的影响。
  3.1虚警概率
  假设所需的信号并不存在, 虚警概率主要取决于噪声的统计特性和检测器的门限设置。 FFT输出信号的振幅用于信号检测器, 这些振幅值服从瑞利分布, 概率密度为[12]
  可以看出虚警概率取决于门限的设置和FFT的点数N。 在不补零的情况下, FFT的点数与匹配滤波器数目一致, 部分匹配滤波器长度越大, 虚警概率越小, 如图6所示。 3.3平均捕获时间
  本文采用的是单次逗留捕获判决方式, 该算法的平均捕获时间为[13]
  Tacq=(2-Pd)(1+KPfa)2Pdqτd(23)
  式中:q为搜索单元数, q等于码的长度或码长度的整数倍, 如果每次更新半个码片长度, 那么q等于两倍的码长, 即q=2M。 假设虚警惩罚时间为kτd, 其中:k为虚警惩罚次数; τd为单次驻留时间。 800 Hz多普勒频移下, 平均捕获时间与部分匹配滤波器长度之间的关系如图8所示。 部分匹配滤波器长度越大, 所需捕获时间越短。
  上述分析从算法的理论特性和捕获性能对部分匹配滤波器进行考察考。 理论特性方面, 部分匹配滤波器的长度影响着算法本身的多普勒频率估计精度和范围、 相关损失以及扇贝损失。 匹配滤波器长度越小, 频率估计精度越低, 频率估计范围越大, 相关损失越小, 最大扇贝损失越大; 然而匹配滤波器长度越小也意味着更大的硬件实现复杂度, 改变匹配滤波器长度无法有效减少扇贝损失的影响。 捕获性能方面, 工程中所常用的捕获评价标准为检测概率、 虚警概率和平均捕获时间, 部分匹配滤波器长度越长, 可以得到更加优秀的虚警概率和捕获时间, 检测概率却也随之降低。   因此在设计PMFFFT捕获系统时, 需要从理论特性和捕获性能双方面考虑如何选取系统的部分匹配滤波器长度。 理论特性方面根据系统环境的最大多普勒频移和系统所能接受的相关损失来确定匹配滤波器长度的限制要求, 再根据对捕获性能如检测概率及捕获时间的需求并考虑硬件实现成本来最终确定部分匹配滤波器的长度。
  5结论
  本文对基于PMFFFT快速捕获算法的理论原理进行了分析, 给出算法的结构, 建立了数学模型。 分别从多普勒频率、 相关损失和扇贝损失以及检测概率、 虚警概率和平均捕获时间几个方面给出部分匹配滤波器长度的改变对算法理论性能和捕获性能的影响, 并给出滤波器长度选择的指导思想。 实际中可以根据用户不同的指标需求选取不同的匹配滤波器长度, 解决了PMFFFT捕获系统设计中部分匹配滤波器长度的选择问题, 为PMFFFT算法的系统设计解决了首要的问题。
  参考文献:
  [1] Holmes J K. Coherent Spread Spectrum Systems[M]. New York: Wiley-Interscience, 1982.
  [2] Povey G J R, Talvitie J. Doppler Compensation and Code Acquisition Techniques for LEO Satellite Mobile Radio Communications[C]// Fifth International Conference on Satellite Systems for Mobile Communications and Navigation, London, UK,1996:16-19.
  [3] Stirling-Gallacher R A, Hulbert A P, Povey G J R. A Fast Acquisition Technique for a Direct Sequence Spread Spectrum Signal in the Presence of a Large Doppler Shift[C]// IEEE 4th International Symposium on Spread Spectrum Techniques and Applications Proceedings, 1996, 1: 156-160.
  [4] Spangenberg S M, Scott I, McLaughlin S, et al. An FFTBased Approach for Fast Acquisition in Spread Spectrum Communication Systems[J]. Wireless Personal Communications, 2000, 13(1-2): 27-55.
  [5] Ward R B. Acquisition of PseudoNoise Signals by Sequential Estimation[J]. Communication Technology, IEEE Transactions on, 1965, 13(4): 475-483.
  [6] Ward R B, Yiu K. Acquisition of PseudoNoise Signals by RecursionAided Sequential Estimation[J]. Communications, IEEE Transactions on, 1977, 25(8): 784-794.
  [7] Van Nee D J R, Coenen A. New Fast GPS CodeAcquisition Technique Using FFT[J]. Electronics Letters, 1991, 27(2): 158-160.
  [8] Lin D M, Tsui J B Y. Acquisition Schemes for Software GPS Receiver[C]//Proceedings of the 11th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS 1998), 1998, 11: 317-325.
  [9] Lin D M, Tsui J B Y. Comparison of Acquisition Methods for Software GPS Receiver[C]//Proceedings of the 13th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS 2000), 2000: 2385-2390.
  [10] Gardner W A, Spooner C M. Detection and Source Location of Weak Cyclostationary Signals: Simplifications of the MaximumLikelihood Receiver[J]. IEEE Transactions on Communications, 1993, 41(6): 905-916.
  [11] Sust M K, Kaufmann R F, Molitor F, et al. Rapid Acquisition Concept for Voice Activated CDMA Communication[C]// IEEE GLOBECOM 90, San Diego, CA, 1990: 1820-1826.
  [12] 唐学术.高动态GPS接收机的捕获与跟踪算法研究[D].北京:北京理工大学, 2008.
  [13] 刘艳华, 赵刚.基于PMFFFT的PN码捕获方法及性能[J].通信技术, 2009, 42(1):24-26.
其他文献
泰勒斯英国先进武器系统(ADW)打算对轻型多用途导弹(LMM)的便携式地对地型号进行改型,以作为可拆卸的步兵直接火力武器。  泰勒斯英国先进武器系统副经理David Beatty向本刊确认公司于12月初在英国国防部靶场从三脚架系统进行了地对地型号的首次验证点火,标志着“便携式轻型多用途导弹性能改进迈出了第一步”。  Beatty又说,在试验中,可拆卸型号使用激光驾束制导达到了“超过6千米”的射程。
期刊
为开发一条不依赖于传统的俄罗斯合作伙伴的新型生产线,“阿尔乔姆”(Artem)和Radionix 两家乌克兰公司已开始一项合作项目,研制温贝尔设计局的R-27空空导弹(北约代号“Alamo”)的防空改进型。  R-27导弹是俄罗斯设计的,苏联时期由莫斯科的温贝尔设计局研发了多种改型,但这些导弹及苏联设计的其他大多数空空导弹的生产线则位于基辅的“阿尔乔姆”工厂。  苏联解体后,“阿尔乔姆”公司继续全
期刊
摘要: 對适用于机载武器的组合导航算法展开研究, 给出了SINS/BDS组合导航状态方程, 推导了以伪距、 伪距率作为观测量的观测方程。 针对传统卡尔曼滤波算法在测量噪声不确定的情况下使用存在缺陷的問题, 研究了新的噪声自适应调整滤波与残差故障检测算法。 对新算法进行了软件设计, 并开展了软件仿真与产品车载试验。 试验结果表明, 所设计的算法能有效增强系统鲁棒性、 提高滤波算法精度, 有效隔离故障
期刊
MBDA 公司英国分公司的工程师设计了一种新颖的制导结构并获得了专利。该设计为导弹头部腾出许多内部空间,从而使更多的地方可用于战斗部;将导弹头部中的制导硬件减至最小,其余大部分硬件装在战斗部后面。  在许多战术导弹中,导引头都是安装在导弹头部的,该位置能够提供连续的视场。战斗部一般位于制导电子设备之后和推进系统之前,这种布局对战斗部来说远远达不到最优,因为导引头和制导组件会对战斗部性能产生非常不利
期刊
欧洲导弹集团(MBDA)致力于一项为期3年、由公司内部研究与发展部门资助的创新项目,通过充分利用活性材料技术来提高战斗部的毁伤效能,该项目的研究结果将对公司导弹的未来发展起到促进作用。  传统战斗部采用惰性或非活性金属外壳(如钢、铝、钨、铜,占战斗部重量的2/3),以沿初始速度矢量方向生成可能的最大毁伤。这种毁伤受破片的初始动能所限。  活性材料技术改为采用“过渡”金属族材料(如钛、钽、锆、铪、铝
期刊
英国已与MBDA公司签订了一份价值3亿英镑(4.62亿美元)的合同,为英国皇家空军(RAF)制造新型的“先进近距空空导弹”(ASRAAM),该消息宣布于2015年9月16日的2015伦敦防务系统与装备国际展上。  这份“能力保持”订单将以新制造的导弹来补充英国皇家空军ASRAAM导弹的整个库存。ASRAAM导弹是英国皇家空军“狂风”GR.4战斗机和“台风”战斗机最主要的近距空空武器,现有库存即将步
期刊
欧洲战斗机集团的“台风”战斗机继续分别与该集团的三家成员公司一起拓展其能力,宣布2015年7月中旬进行飞机新的能力试验。  7月15日,空客防务与航天公司宣布已完成“台风”的新型气动改进成套组件(AMK)的飞行试验,声称该组件显著提升战斗机性能,最大升力提高25%,最大迎角提高45%,战斗机横滚速率提高100%。  参与该项目的试飞员拉斐尔·贝尔特拉梅称其“令人印象深刻”。他补充道,“最初以为这不
期刊
摘要: 基于脉冲隔板总开孔面积一定的情况下, 研究了不同级间开孔形式对一脉冲壳体绝热层烧蚀和流动损失的影响, 对带硬质隔板的双脉冲发动机进行了内流场数值计算。 结果表明: 隔板背壁区涡流会加重绝热层烧蚀, 燃气再附着点附近烧蚀情况最严重; 改变开孔形状和布局对隔板流动损失影响不大, 但优化开孔倒角设计可改善隔板孔流通能力; 减小隔板外缘孔径并加开隔板中心孔能有效减轻绝热层烧蚀。  关键词: 固体火
期刊
摘要: 为设计高效记忆合金自调节式节流制冷器, 介绍了自调节式节流制冷器的工作原理, 给出利用记忆合金驱动阀针进行流量调节的机理。 对记忆合金通过控制节流阀针的位移调节制冷工质流量过程进行计算分析, 模拟了在不同初始高压和节流孔径的条件下, 节流制冷器流量随调节阀针位移的变化关系。 为实现螺旋翅片管换热器高效换热, 对其进行了系统的强化换热分析, 分别对内管管径、 芯管管径以及翅片尺寸等影响其换热
期刊
摘要: 基于一种半装药的高过载地面模拟试验发动机, 针对不同推进剂铝粉含量和绝热层材料开展了三次发动机过载试验, 获得了绝热层宏观烧蚀规律, 通过电镜分析了EPDM的烧蚀模型。 结果表明:(1)三元乙丙橡胶绝热层较丁羟基绝热层更抗烧蚀、 抗冲刷; (2)推进剂的铝粉含量越高, 燃气对三元乙丙绝热层的冲刷越严重, 绝热层的烧蚀率越大; (3)通过试验分析, 验证了三元乙丙绝热层的碳化、 热解、 基体
期刊