一种样点缺失条件下的跳频信号参数估计方法_王洪斌.pdf
《一种样点缺失条件下的跳频信号参数估计方法_王洪斌.pdf》由会员分享,可在线阅读,更多相关《一种样点缺失条件下的跳频信号参数估计方法_王洪斌.pdf(8页珍藏版)》请在咨信网上搜索。
1、 第2卷 第2期V o l.2 N o.2 2 0 2 3年4月 J o u r n a l o f A r m y E n g i n e e r i n g U n i v e r s i t y o f P L A A p r.2 0 2 3一种样点缺失条件下的跳频信号参数估计方法王洪斌,张邦宁,王 桁,吴彬彬,郭道省(陆军工程大学 通信工程学院,江苏 南京 2 1 0 0 0 7)摘要:在样点随机缺失条件下研究跳频信号参数估计问题具有现实意义。针对样点缺失条件下线性时频分析方法失效的问题,提出了一种基于正交匹配追踪(o r t h o g o n a l m a t c h i n g
2、 p u r s u i t,OMP)和卡尔曼滤波(K a l m a n f i l t e r,K F)的跳频信号参数实时估计方法。该方法对信号加滑动窗,将窗内样点随机缺失建模为一个信号稀疏表示问题,傅里叶正交基作为过完备字典,利用OMP直接估计窗内信号频率而不需要恢复信号。K F针对估计得到的信号频率进行平滑,当频率跳变时,K F的频率预测值将严重偏离历史值和频率估计值,偏离程度作为跳变时刻估计的依据。仿真结果和对比实验表明,在样点没有缺失时,该方法具有更优的跳频信号参数估计性能,在滑动窗长足够并满足信号稀疏度要求时,即使在样点缺失的条件下,依然可以获得有效的跳频信号参数估计结果。关键词
3、:跳频信号参数估计;样点缺失;正交匹配追踪;卡尔曼滤波 中图分类号:T N 9 7 5D O I:1 0.1 2 0 1 8/j.i s s n.2 0 9 7-0 7 3 0.2 0 2 2 0 9 1 4 0 0 1P a r a m e t e r E s t i m a t i o n M e t h o d o fF r e q u e n c y H o p p i n g S i g n a l s w i t h M i s s i n g O b s e r v a t i o n s WANG H o n g b i n,Z HANG B a n g n i n g,WAN
4、G H e n g,WU B i n b i n,GUO D a o x i n g(C o l l e g e o f C o mm u n i c a t i o n s E n g i n e e r i n g,A r m y E n g i n e e r i n g U n i v e r s i t y o f P L A,N a n j i n g 2 1 0 0 0 7,C h i n a)A b s t r a c t:I t i s o f p r a c t i c a l s i g n i f i c a n c e t o s t u d y t h e p a
5、r a m e t e r e s t i m a t i o n o f t h e f r e q u e n c y h o p p i n g(F H)s i g n a l s i n t h e c a s e o f r a n d o m m i s s i n g o b s e r v a t i o n s.I n v i e w o f t h e f a i l u r e o f t h e l i n e a r t i m e f r e q u e n c y a n a l y s i s i n t h e c a s e o f m i s s i n
6、g o b s e r v a t i o n s,a f r e q u e n c y h o p p i n g s i g n a l p a r a m e t e r e s t i m a t i o n m e t h o d i s p r o p o s e d i n t h i s p a p e r,b a s e d o n o r t h o g o n a l m a t c h i n g p u r s u i t(OMP)a n d K a l m a n f i l t e r(K F).B y t h i s m e t h o d,a s l i d
7、 i n g w i n d o w i s a d d e d t o t h e s i g n a l s,a n d t h e r a n d o m m i s s i n g o b s e r v a t i o n s i n t h e w i n d o w i s m o d e l e d a s a s i g n a l s p a r s e r e p r e s e n t a t i o n.A s a n o v e r c o m p l e t e d i c t i o n a r y,t h e F o u r i e r o r t h o n
8、 o r m a l b a s e s m a k e s u s e o f OMP t o e s t i m a t e t h e f r e q u e n c y o f t h e s i g n a l i n t h e w i n d o w w i t h o u t r e s t o r i n g t h e s i g n a l s.K F p e r f o r m s s m o o t h l y o n t h e e s t i m a t e d s i g n a l f r e q u e n c y.Wh e n t h e f r e q
9、u e n c y c h a n g e s,t h e f r e q u e n c y p r e d i c t i o n v a l u e o f K F w i l l s e r i o u s l y d e v i a t e f r o m t h e h i s t o r i c a l v a l u e a n d t h e f r e q u e n c y e s t i m a t e d v a l u e,a n d t h e d e g r e e o f d e-v i a t i o n w i l l o f f e r t h e s
10、u p p o r t f o r h o p p i n g t i m e e s t i m a t i o n.T h e s i m u l a t i o n r e s u l t s a n d c o m p a r a t i v e e x p e r i m e n t s s h o w t h a t t h e p a r a m e t e r e s t i m a t i o n o f t h e p r o p o s e d m e t h o d o u t p e r f o r m s t h a t o f o t h e r m e t h
11、o d s,w h e n t h e o b s e r v a-t i o n s a r e n o t m i s s i n g.W i t h t h e s u f f i c i e n t s l i d i n g w i n d o w l e n g t h a n d t h e e l i g i b l e s i g n a l s p a r s i t y r e q u i r e m e n t,e f f e c-t i v e F H s i g n a l p a r a m e t e r e s t i m a t i o n r e s u
12、l t s c a n a l s o b e o b t a i n e d e v e n i n t h e c a s e o f m i s s i n g o b s e r v a t i o n s.K e y w o r d s:p a r a m e t e r e s t i m a t i o n o f f r e q u e n c y h o p p i n g s i g n a l s;m i s s i n g o b s e r v a t i o n s;o r t h o g o n a l m a t c h i n g p u r s u i t;
13、K a l m a n f i l t e r 收稿日期:2 0 2 2-0 9-1 4第一作者:王洪斌,博士研究生,主要研究通信抗干扰、跳频信号分析和处理,w a n g_h o n g b i n_c o mm1 2 6.c o m。通信作者:郭道省,教授,博士生导师,主要研究网络空间安全,x y z g f g s i n a.c o m。跳频(f r e q u e n c y-h o p p i n g,F H)信号的载波频率以某种伪随机的方式变化,由于其在低截获概率、灵活的组网能力、抗干扰和多径衰落方面的固有优势,F H信 号 已 广 泛 应 用 于 卫 星 通 信1、无 线 通信
14、2、物理层安全3、智能电网4-5、水下通信6-7、物联网(i n t e r n e t o f t h i n g s,I o T)技 术8-9和 无 人 机(u n m a n n e d a e r i a l v e h i c l e,UAV)1 0等领域。在干扰分析、电子对抗和通信安全等领域中,估计F H信号的参数并对其进行跟踪是一项至关重要的且具有挑战性的任务,尤其在跳频模式未知的情况下。针对如F H信号这种时变非平稳信号的处理,经典的方法有短时傅里叶变换(s h o r t-t i m e F o u r i e r t r a n s f o r m,S T F T)1 1-
15、1 2、魏格纳分布(W i g n e r-V i l l e d i s t r i b u t i o n s,WV D)1 3及其改进算法等时频分析方法,还有利用F H信号稀疏性的压缩感知(c o m-p r e s s i v e s e n s i n g,C S)1 4和稀疏表示技术1 5等。但传统方法往往需要积累足够长的接收数据,难以应用于F H信号实时跟跳、参数快速估计等场景中。为了解决跳频信号实时处理的问题,出现了许多应对方法1 6-1 8,其中最具代表性的是利用自回归滑动平均(a u t o r e g r e s s i v e m o v i n g a v e r a
16、 g e,A RMA)模型对跳频信号的频率和跳变时刻实时跟踪。L i u等1 6首先将这一思想应用于F H信号的处理中,利用接收到的数据建立一个时域A RMA模型,根据模型系数估计信号频率,并预测下一时刻的接收数据,根据实际接收数据和预测数据检测跳变时刻。该方法虽然达到了实时跟踪F H信号并估计其参数的目的,但对接收信号的信噪比(s i g n a l-t o-n o i s e r a t i o,S N R)有较高的要求,而且频率估计精度不高。文献1 7 在时域A RMA模型的基础上提出了一种根据跳周期修正跳变时刻估计的方法,其核心思想就是对时域A RMA模型估计得到的跳变时刻进行差分聚类
17、,提高跳变时刻估计的正确率。但数据聚类仍是一种批处理的思想,降低了信号处理的实时性。文献1 8 在总结了前人工作的基础上认为,频率估计精度不高的原因在于采用时域信息建立的A R-MA模型对频域信息的表征能力不足。因此,文献1 8 提出改进算法,先将接收信号变换到频域,再随时间建立A RMA模型,这样的模型系数兼顾信号时频域特点,提高了频率估计能力。然而,A RMA模型虽然可以依据历史数据对接收信号做出预测,但受噪声的影响,其预测波动较大,低信噪比下的F H信号参数估计精度受到较大限制。此外,上述方法都没有考虑样点缺失条件下的跳频信号实时处理和参数估计问题。实际应用中,受冲击干扰、噪声、遮挡、信
18、道衰落、硬件问题和天线旁瓣接收等因素的影响,接收信号可能面临采样值异常,导致部分采样点不可用。文献1 9 在样点随机缺失的条件下研究了跳频信号的参数估计问题,提出了一种波形自适应时频核设计,结合了自动优化的预滤波窗和数据相关的自适应最优核(a d a p t i v e o p-t i m a l k e r n e l,A O K)。这种设计可以有效地抑制由于丢失观测样点引起的交叉项和伪影,同时保留跳频信号的自项。然而,该方法需要积累一定的数据量,不是一种实时处理方法,而且计算复杂度较高。本文 提 出 了 一 种 正 交 匹 配 追 踪(o r t h o g o n a l m a t c
19、 h i n g p u r s u i t,OMP)结合卡尔曼滤波(K a l m a n f i l t e r,K F)(OMP-K F)的跳频信号参数估计方法。利用OMP估计一个固定长度的滑动窗内的信号频率,估计结果作为K F的观测值,对下一时刻的频率进行预测,将预测值的突变作为跳变时刻估计的依据。不同于A RMA频域建模的是,利用K F对数据的平滑能力,降低了噪声的不利影响,提高了参数估计精度。1 信号模型在观测时间内,接收的F H信号在t时刻的表达式为s(t)=Hh=1AhTt-h T()ej 2 fh(t-h T)+v(t)(1)式中:j=-1,T为F H信号的跳周期,H表示观测
20、时间内的总跳数;Tt()表示归一化窗函数,当t-T/2,T/2(时,该窗函数等于1,否则等于0;Ah、fh分别表示第h跳的复振幅和瞬时频率,v(t)表示加性高斯白噪声。接收信号经过采样率为fs(t=1/fs表示采样间隔)的奈奎斯特采样后可以表示为sk=Hh=1Ahej 2 fhkt+vk 0kN-1(2)采用加窗方式对接收信号逐样点滑动处理,假设第n个滑动窗窗长为N(N小于F H信号的驻留时间),接收信号sn=s0,s1,sN-1TCN 1。在实际应用中,接收信号可能面临着采样点随机缺失的情况。假设sn的采样点随机缺失,缺失的位置满足均匀分布,可表示为yn=nsn(3)式中n为LN(L()&?
21、fn22(),a c c e p t H1i f?fn22()|?fn22(),a c c e p t H0(2 2)式中和表示检测阈值。将接受H1时的滑动窗的中间样点位置作为跳变时刻的估计点,如果存在连续多个跳变时刻检测点,则以多个点中的第一个点作为跳变时刻的估计点。依据估计得到的跳变时刻划分集合f中的频率估计值,对两次跳变时刻之间的所有频率估计值f按照数值大小排序,剔除异常值,然后计算剩余频率估计值f的均值作为该驻留时间内的频率估计值。当采样率fs=6 4 0 k H z,观测时间内有6跳,第1跳到第6跳的频率分别为1 3 0、2 8 0、6 0、1 8 0、3 0和2 3 0 k H z
22、。S N R=5 d B,滑动窗长分别为1 6、3 2和6 4时,在样点缺失和不缺失条件下的F H信号频率实时估计的时频表示如图3所示。如图3(g)、图3(h)和图3(i)所示,即使样点缺失4 0%时,预测值也可以很好地收敛于真实值。此外,随着窗长的增加,预测值更加接近真实值,但窗长的增加也带来了更大的跟踪延迟,如图3(a)、图3(b)和图3(c)所示,以真实的时频表示为标准,窗长为6 4的估计值和预测值比窗长为3 2或1 6的估计值和预测值具有更大的跟踪延时。但窗长是已知的,滑动延时也是固定的,只需要将这段延时补偿回去,就可以准确估计跳变时刻。36第2期 王洪斌,等:一种样点缺失条件下的跳频
23、信号参数估计方法图3 不同窗长样点缺失或不缺失时的时频表示(S N R=5 d B)2.3 算法复杂度分析对于OMP算法,在过完备字典与信号残差的相关计算中,需要计算一个NL型矩阵和L1型矩阵的内积,每一行做乘法L次,做加法L-1次,复杂度为O(N(2L-1)。在残差更新的过程中,迭代到第k步的支撑集的大小为Lk,残差更新的复杂度为2k L。由于该方法不对信号重构,故不考虑OMP重构信号复杂度。在迭代K步后,总计算复杂度为KN(2L-1)+2(1+2+K)L 2K LN+L K2(2 3)由于K6 d B时,即使样点随机缺失2 0%,也可以获得9 0%以上的估计正确率。当窗长取6 4,S N
24、R5 d B时,即使样点随机缺失4 0%,同样可以获得9 0%以上的估计正确率。从图7可以看出,样点缺失对频率估计的影响不容小觑,在高S N R下,样点缺失比例相差2 0%时,频率估计误差相差约一个数量级。56第2期 王洪斌,等:一种样点缺失条件下的跳频信号参数估计方法图6 样点缺失时不同窗长的正确跳变时刻估计比例图7 样点缺失时不同窗长的频率估计相对误差虽然如此,相较于现有文献中的大部分方法,OMP-K F方法在样点缺失条件下依然有效,而且即便样点随机缺失2 0%,当窗长取3 2时,在高S N R下,频率估计相对误差依然可以达到1 0-2,这充分体现了该方法在F H信号参数估计方面具有较强鲁
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 一种 缺失 条件下 信号 参数估计 方法 王洪斌
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。