论文:Super-Resolution Delay-Doppler Estimation for OFDM Passive Radar

雷达与通信系统的对比

1.信息源不同。雷达系统的信息源来自于目标,可以说是来自于信道,追寻的是信息源的存在及其参数;通信系统的信息源来自于发射机的输入,目的是信息的传输。

2.在大多数情况下,雷达系统使用相同的天线来发送和接收信号,而通信系统具有不同的发送和接收天线。二个系统的射频部分相似,频段也可能重合,但是采用的信号处理有很大差异。

3.对于干扰的定义,二种系统也是不同的,甚至不同的雷达系统对干扰的定义也是不同的。

L0范数是指向量中非0的元素的个数如果我们用L0范数来规则化一个参数矩阵W的话,就是希望W的大部分元素都是0。换句话说,让参数W是稀疏的。

L1范数是指向量中各个元素绝对值之和,也有个美称叫"稀疏规则算子"(Lasso regularization)。

L1范数和L0范数都可以实现稀疏

无源雷达是一种不用发射机发射能量而靠接受温热物体或他源反射的微波能量探测目标的雷达,有天线和灵敏度极高的接收装置。无源雷达鉴别目标的能力,主要取决于目标之间的表面温差和目标的反射系数,天线波束与目标之间的入射余角,无线极化和波束宽度与接收机的最小可检测电子等。

无源雷达自身不发射信号,既带来优点也带来缺点。由于依赖于第三方发射机,操作员对辐射源无法主动控制,在被探测目标保持无线电静默、照射器又不工作的情况下,无源雷达就成了无源之水,不能发挥作用。此外,一些发射机的有效辐射功率较低,易受干扰和空射诱饵的影响而且要求发射机与目标、目标与接收机以及接收机与发射机之间信号不受阻挡,限制了无源雷达的使用。

论文思路:

  1. 由于发射信号是未知的,常规的基于匹配滤波器的接收机不容易实现,即直接信号和杂波对雷达信号处理造成干扰,给目标检测和参数估计带来困难。
  2. 一些无源雷达不解调接收信号,为了消除监督信道(SCs)中不需要的echo,引入一个附加的参考通道(RC),用于收集传输信号作为被动检测的参考。但是将RC中的信号视为噪声模板,实现近似匹配滤波器在实际应用中容易受到噪声的干扰,造成性能的下降。
  3. 由于调制类型通常是已知的,因此可以首先通过解调来估计数据信号,以便实现匹配的滤波器,解调提供了比直接使用RC信号更好的精度。此外,由于数据可以根据接收信号和估计的信道状态进行解调,因此不需要RC,降低了无源雷达的复杂性
  4. 使用OFDM进行调制。广播信号可以看作是一种特殊形式的脉冲多普勒雷达信号。目标可以使用相应的匹配滤波器组进行搜索。
  5. 解调误差导致信号模型中存在等效的脉冲噪声,从而增加杂波和强目标的旁瓣。因为目标的反射信号非常小,因此解调误差会对目标的检测和估计造成非常大的性能的下降。
  6. 目标和杂波的数量较少,因此目标和杂波反射信号在Delay-Doppler域是稀疏的;在正常工作条件下,通信系统的解调错误率通常很低,因此解调误差也是稀疏的。因此考虑了CS即压缩感知算法。
  7. CS算法关注的是可以用有限离散字典稀疏表示的信号,因此假设Interest频率位于固定频率网格上。然而,由于目标和杂波通常由连续域(时延和多普勒频率)的参数指定,离散化往往导致模型失配和恢复退化。
  8. 为了克服传统的基于稀疏的网格不匹配问题,论文引入了连续稀疏恢复超分辨率方法
  9. 因此使用L0范数表示目标和杂波反射信号的稀疏性,L1范数表示解调误差信号的稀疏性,目标和杂波的参数可以通过求解凸半定规划问题(SDP)来估计。论文同时提出了一种基于乘法器交替方向法(ADMM)的迭代算法,来解决这个SDP问题,其中每次迭代都涉及到封闭计算。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
% eta数据 % 输入:24个激光干涉数据 % 输出:6列eta数据 function [eta,etatstart,etaNpoints]=etadata(t0,sci,scib,tao,epsilon) global Fs M m1 m2 m3 m4 m5 m6; Npoints=length(sci(:,1)); prePoints=10*Fs+2*M; endPoints=2*M; etaNpoints=Npoints-prePoints-endPoints; etatstart=t0+prePoints/Fs; eta=zeros(etaNpoints,12); epsilon1=epsilon(:,1);epsilon2=epsilon(:,2);epsilon3=epsilon(:,3); epsilon4=epsilon(:,4);epsilon5=epsilon(:,5);epsilon6=epsilon(:,6); tao1=tao(:,1);tao2=tao(:,2);tao3=tao(:,3); tao4=tao(:,4);tao5=tao(:,5);tao6=tao(:,6); for ii=prePoints+1:prePoints+etaNpoints ti=t0+(ii-1)/Fs; Ltime=lasertravelTime(ti,0); delay0=Ltime(1:6)*Fs; delay0M=ceil(delay0); delay0e=delay0M-delay0; doppler0=1-Ltime(7:12); d3epsilon5tao5=epsilon5(ii-delay0M(3)+M:-1:ii-delay0M(3)-M+1)-tao5(ii-delay0M(3)+M:-1:ii-delay0M(3)-M+1); d1epsilon6tao6=epsilon6(ii-delay0M(1)+M:-1:ii-delay0M(1)-M+1)-tao6(ii-delay0M(1)+M:-1:ii-delay0M(1)-M+1); d2epsilon4tao4=epsilon4(ii-delay0M(2)+M:-1:ii-delay0M(2)-M+1)-tao4(ii-delay0M(2)+M:-1:ii-delay0M(2)-M+1); d5epsilon3tao3=epsilon3(ii-delay0M(5)+M:-1:ii-delay0M(5)-M+1)-tao3(ii-delay0M(5)+M:-1:ii-delay0M(5)-M+1); d6epsilon1tao1=epsilon1(ii-delay0M(6)+M:-1:ii-delay0M(6)-M+1)-tao1(ii-delay0M(6)+M:-1:ii-delay0M(6)-M+1); d4epsilon2tao2=epsilon2(ii-delay0M(4)+M:-1:ii-delay0M(4)-M+1)-tao2(ii-delay0M(4)+M:-1:ii-delay0M(4)-M+1); d3tao2tao5=tao2(ii-delay0M(3)+M:-1:ii-delay0M(3)-M+1)-tao5(ii-delay0M(3)+M:-1:ii-delay0M(3)-M+1); d1tao3tao6=tao3(ii-delay0M(1)+M:-1:ii-delay0M(1)-M+1)-tao6(ii-delay0M(1)+M:-1:ii-delay0M(1)-M+1); d2tao1tao4=tao1(ii-delay0M(2)+M:-1:ii-delay0M(2)-M+1)-tao4(ii-delay0M(2)+M:-1:ii-delay0M(2)-M+1); eta(ii-prePoints,1)=sci(ii,1)-(epsilon(ii,1)-tao(ii,1))/2-doppler0(3)*fracdelay(delay0e(3))*(d3epsilon5tao5+ d3tao2tao5)/2; eta(ii-prePoints,2)=sci(ii,2)-(epsilon(ii,2)-tao(ii,2))/2-doppler0(1)*fracdelay(delay0e(1))*(d1epsilon6tao6+ d1tao3tao6)/2; eta(ii-prePoints,3)=sci(ii,3)-(epsilon(ii,3)-tao(ii,3))/2-doppler0(2)*fracdelay(delay0e(2))*(d2epsilon4tao4+ d2tao1tao4)/2; eta(ii-prePoints,4)=sci(ii,4)-(epsilon(ii,4)-tao(ii,4))/2-doppler0(5)*fracdelay(delay0e(5))*d5epsilon3tao3/2+(tao(ii,1)-tao(ii,4))/2; eta(ii-prePoints,5)=sci(ii,5)-(epsilon(ii,5)-tao(ii,5))/2-doppler0(6)*fracdelay(delay0e(6))*d6epsilon1tao1/2+(tao(ii,2)-tao(ii,5))/2; eta(ii-prePoints,6)=sci(ii,6)-(epsilon(ii,6)-tao(ii,6))/2-doppler0(4)*fracdelay(delay0e(4))*d4epsilon2tao2/2+(tao(ii,3)-tao(ii,6))/2; eta(ii-prePoints,7)=(scib(ii,1)-sci(ii,1))/m5; eta(ii-prePoints,8)=(scib(ii,2)-sci(ii,2))/m6; eta(ii-prePoints,9)=(scib(ii,3)-sci(ii,3))/m4; eta(ii-prePoints,10)=(scib(ii,4)-sci(ii,4))/m3; eta(ii-prePoints,11)=(scib(ii,5)-sci(ii,5))/m1; eta(ii-prePoints,12)=(scib(ii,6)-sci(ii,6))/m2; end end
05-25

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

飞大圣

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值