CN109085574B - Ofdm雷达通信一体化固定平台系统的信号处理方法 - Google Patents

Ofdm雷达通信一体化固定平台系统的信号处理方法 Download PDF

Info

Publication number
CN109085574B
CN109085574B CN201811218839.7A CN201811218839A CN109085574B CN 109085574 B CN109085574 B CN 109085574B CN 201811218839 A CN201811218839 A CN 201811218839A CN 109085574 B CN109085574 B CN 109085574B
Authority
CN
China
Prior art keywords
signal
pulse
baseband
frequency domain
cyclic prefix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811218839.7A
Other languages
English (en)
Other versions
CN109085574A (zh
Inventor
刘高高
王钊
蔡晶晶
秦国栋
武斌
鲍丹
李鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Shengxin Technology Co ltd
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201811218839.7A priority Critical patent/CN109085574B/zh
Publication of CN109085574A publication Critical patent/CN109085574A/zh
Application granted granted Critical
Publication of CN109085574B publication Critical patent/CN109085574B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种OFDM雷达通信一体化固定平台系统的信号处理方法,主要解决现有技术通信数据的误比特率高,恢复的一体化脉冲参考信号误差大,脉冲压缩处理效果差的问题。其实现步骤是:1、设置回波信号条件;2、将回波信号进行下变频处理和采样处理,得到处理后的基带信号,并移除该信号的循环前缀;3、对移除循环前缀后的信号进行傅里叶变换;4、对傅里叶变换后的信号进行解码和判决,得到通信信息;5、对通信信息进行编码和调制,得到参考信号;6、利用参考信号对移除循环前缀后的信号进行脉冲压缩处理。本发明具有通信数据误比特率低,脉冲压缩处理效果好的优点,可用于OFDM雷达通信一体化固定平台系统。

Description

OFDM雷达通信一体化固定平台系统的信号处理方法
技术领域
本发明属于雷达技术领域,特别涉及一种信号处理方法,可用于OFDM雷达通信一体化固定平台系统。
背景技术
传统的电子装备系统仅单独具有雷达系统的功能或通信系统的功能,所能实现的功能相对单一。当在要求同时实现两种功能的情况下,则造成系统体积庞大,且频率、时间等不同资源的占用,造成资源利用率降低。近些年出现了将雷达技术和通信技术相结合的电子装备系统,结合后的雷达通信一体化系统能够很好地解决上述问题,此种系统受到越来越多的关注。
雷达通信一体化,是指同一电子装备系统同时具备雷达功能以及通信功能的电子设备,这种设备在工作时能够实时完成雷达信号处理和通信信号处理两种功能。目前雷达通信一体化的实现方式包含以下三种:
一是分时使用雷达和通信功能,这种实现方式通过时间来分配两种功能,存在时间资源上利用率低的缺点。
二是将频段划分,分别用于雷达和通信功能,这种实现方式在频谱利用方面存在浪费。
三是将通信信号应用于雷达系统,这种实现方式是在接收端为未知的雷达参考信号的条件下,利用通信功能来恢复参考信号,进而完成后续的雷达信号处理,从而实现雷达通信一体化。这种实现方式提高了资源利用率。其中,一体化固定平台系统的信号是指一体化系统相对于检测目标保持相对静止的传输信号。在现有研究中,当雷达通信一体化系统处在相对于目标静止的条件下时,通信数据的误比特率较高,导致恢复的参考信号较原始参考信号之间存在误差。同时,在后续进行脉冲压缩的处理效果差,无法满足一体化系统的需要。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种OFDM雷达通信一体化固定平台系统的信号处理方法,以降低系统通信数据的误比特率,有效恢复参考信号,提高脉冲压缩的处理的效果。
本发明的技术思路是,将大带宽通信信号应用于雷达系统中,利用通信功能对传输的通信信息进行恢复,利用恢复的信号重建参考信号实现大带宽信号的雷达脉冲压缩处理。其实现步骤包括如下:
(1)雷达接收机通过天线接收调制方式为OFDM的时间回波脉冲信号s'(t),该回波脉冲信号的有效长度为τs
(2)对接收到的时间回波脉冲信号s'(t)依次进行下变频处理和采样率为fs的采样处理,得到基带脉冲信号s′r(n),其中,n为正整数,表示采样点数,基带脉冲的起始位置为No,脉冲总长度为Nsum,前M-1个脉冲单元为循环前缀,M为时间回波脉冲信号的距离单元数;
(3)对基带脉冲信号s′r(n)中的循环前缀进行移除,得到移除循环前缀后的信号sr(n);
(4)对移除循环前缀后的信号sr(n)进行傅里叶变换,得到基带频域信号Sr(k),其中k为正整数,表示脉冲索引;
(5)对基带频域信号Sr(k)进行通信处理,得到该基带频域信号中包含的通信信息
Figure BDA0001834244300000026
(5a)将基带频域信号Sr(k)向后延迟一个数据单元,并在空余位置处补零,得到延迟后的基带信号S′r(k);
(5b)对延迟后的基带信号S′r(k)进行共轭运算,得到延迟共轭信号S″r(k);
(5c)将基带信号Sr(k)与延迟共轭信号S″r(k)的对应元素逐个进行乘法运算,得到解码结果序列Sc′;
(5d)对解码结果序列Sc′中的每个元素进行判决:当解码结果序列Sc′中元素的模值大于π/2时判定为1,小于等于π/2时判定为0,得到全部元素的判定结果,即基带频域信号中包含的通信信息
Figure BDA0001834244300000021
(6)利用基带频域信号中包含的通信信息
Figure BDA0001834244300000022
构建参考信号sre(n):
(6a)对通信信息
Figure BDA0001834244300000023
按照下式进行编码运算,得到2DPSK信号
Figure BDA0001834244300000024
Figure BDA0001834244300000025
其中,
Figure BDA0001834244300000031
表示前一个通信信息数据点,
Figure BDA0001834244300000032
表示前一个2DPSK信号数据点,当n-1=0时,
Figure BDA0001834244300000033
为0,
Figure BDA0001834244300000034
表示异或运算;
(6b)对2DPSK信号
Figure BDA0001834244300000035
按照下式进行相位调制,得到雷达处理所需的参考信号sre(n):
Figure BDA0001834244300000036
其中,j表示虚数单位;
(7)利用雷达处理所需的参考信号sre(n)以及移除循环前缀后的信号sr(n)进行脉冲压缩处理,得到脉冲压缩后信号为src(n)。
本发明与现有技术相比具有以下优点:
第一,相对于传统的信号处理方法,本发明将雷达处理和通信处理进行了融合,实现信号处理机的雷达通信一体化;
第二,本发明利用通信处理获得原始通信信息并重建参考信号,达到了在接收端未知参考信号的条件下实现雷达的脉冲压缩处理。
附图说明
图1是本发明的实现流程图;
图2是用本发明对回波信号进行通信处理后的仿真误比特率图;
图3是用本发明对雷达信号的脉冲压缩仿真结果图。
具体实施方式
下面结合附图对本发明做进一步详细的描述。
参照图1,本发明的具体实现步骤如下:
步骤1、雷达接收机通过天线接收调制方式为OFDM的时间回波脉冲信号s'(t)。
(1a)设置发射信号条件:
雷达系统中发射信号形式多样和信号调制方式多变,本实例以雷达系统中常用的线性调频信号sLFM(t)为例,该发射信号的频率按照一次函数规律进行线性变化,其数学表达式为:
Figure BDA0001834244300000037
其中,rect(·)表示矩形函数,用于控制信号持续时间,τ为信号持续时间,f0表示信号的载波频率,A表示信号幅度,μ=B/τ表示调频斜率,B表示带宽;
(1b)设一体化雷达通信信号中通信信息采用2DPSK方式进行编码,雷达接收机通过天线接收到的时间回波脉冲信号为s'(t),其数学表达式为:
Figure BDA0001834244300000041
其中,
Figure BDA0001834244300000042
表示时间延迟,R为一体化系统与目标之间的距离,c为光速,该回波脉冲信号s'(t)的有效长度为τs
步骤2、获取基带脉冲信号。
(2a)对接收到的时间回波脉冲信号s'(t)进行下变频处理,即去除载波频率得到下变频处理后的基带信号;
(2b)对下变频处理后的基带信号进行采样率为fs的采样处理,得到基带脉冲信号s′r(n),该基带脉冲信号为离散的矩阵信号,其中,n为正整数,表示采样点数;基带脉冲s′r(n)的起始位置为No,脉冲总长度为Nsum;前M-1个脉冲单元为循环前缀,M为时间回波脉冲信号的距离单元数;采样率fs为根据信号处理需求不同的常量。
步骤3、对基带脉冲信号s′r(n)中的循环前缀进行移除,得到移除循环前缀后的信号为sr(n)。
(3a)计算回波脉冲信号中的提取位置:Ns=N0+M-1;
(3b)以提取位置Ns为起点,提取基带信号s′r(n)的后Nsum-Ns个单元,得到剩余回波脉冲信号sr(n);
步骤4、对剔除循环前缀的信号进行傅里叶变换。
对除去循环前缀的信号sr(n)按照傅里叶变换得到基带频域信号Sr(k),变换方式为:
Sr(k)=FFT(sr(n)),
其中,k为正整数,表示脉冲索引,FFT表示傅里叶变换。
步骤5、恢复通信信息
Figure BDA0001834244300000043
(5a)将基带频域信号Sr(k)向后延迟一个数据单元,并在空余位置处补零,得到延迟后的基带信号S′r(k);
(5b)对延迟后的基带信号S′r(k)进行共轭运算,得到延迟共轭信号
Figure BDA0001834244300000051
*表示共轭;
(5c)将基带信号Sr(k)与延迟共轭信号S″r(k)的对应元素逐个进行乘法运算,得到解码结果序列Sc′;
(5d)对解码结果序列Sc′中的每个元素进行判决:
当解码结果序列Sc′中元素的模值大于π/2时判定为1,小于等于π/2时判定为0,得到全部元素的判定结果,该判定结果即为基带频域信号中包含的通信信息
Figure BDA0001834244300000052
步骤6、构建雷达信号处理所需的参考信号sre(n)。
(6a)对通信信息
Figure BDA0001834244300000053
按照下式进行编码运算,得到2DPSK信号,记为
Figure BDA0001834244300000054
Figure BDA0001834244300000055
其中,
Figure BDA0001834244300000056
表示前一个通信信息数据点,
Figure BDA0001834244300000057
表示前一个2DPSK信号数据点,当n-1=0时,
Figure BDA0001834244300000058
为0,
Figure BDA0001834244300000059
表示异或运算;
(6b)对2DPSK信号
Figure BDA00018342443000000510
按照下式进行相位调制,得到雷达处理所需的参考信号sre(n):
Figure BDA00018342443000000511
其中,j表示虚数单位。
步骤7、利用参考信号sre(n),对除去循环前缀的信号sr(n)进行脉冲压缩处理,得到脉冲压缩后信号为src(n)。
(7a)对参考信号sre(n)进行共轭运算,得到参考信号sre(n)的频域共轭
Figure BDA00018342443000000512
(7b)将除去循环前缀的信号sr(n)与参考信号sre(n)的频域共轭
Figure BDA00018342443000000513
对应元素相乘,并进行快速傅里叶逆变换,完成信号的脉冲压缩。本实例通过计算脉冲压缩后信号src(n)的峰值位置对目标进行判断。
本发明的效果可通过以下仿真进一步阐述。
1.仿真条件:
仿真设定雷达通信一体化系统的工作载频为10GHz,信号带宽为20MHz,脉冲宽度为30μs,脉冲重频为10KHz,雷达通信一体化系统距离目标10km,通信信息采用随机生成的方式生成,通信编码方式采用2DPSK方式进行编码。
2.仿真内容:
仿真1:在仿真条件下,采用本发明在不同信噪比下,对比不同信噪比条件下的通信功能的误码率,结果如图2所示。
仿真2:在仿真条件下,固定信噪比为20dB条件下,本发明雷达处理脉冲压缩结果如图3所示。
3.仿真分析:
从图2中可以发现,随着信噪比的提高,系统的通信误码率逐步下降,通过和理论的2DPSK误码率进行对比发现,该方式能够保持较低的误码率,能够有效的对通信信息进行传输。
从图3中可以发现,利用通信恢复的信息重建参考信号能够实现雷达的脉冲压缩处理,在图中,空间场景中分布有多个目标,各个目标之间的距离差为166.7米,其中第二个目标对应系统的场景中心,在脉冲压缩后能够发现图中存在四个目标。通过系统采样率与脉冲宽度的乘积得到场景所对应的采样点数目,将其换算到本系统中能够和各个峰值点对应。此外进行脉冲压缩后系统的主副瓣之比普遍超过了13dB,相比于未加窗的线性调频信号具备相似或者更好的旁瓣性能,在雷达通信一体化系统相对于目标静止的条件下,系统的性能能够保持,不存在明显的损失。

Claims (2)

1.一种OFDM雷达通信一体化固定平台系统的信号处理方法,包括如下步骤:
(1)雷达接收机通过天线接收调制方式为OFDM的时间回波脉冲信号s'(t),该时间回波脉冲信号的有效长度为τs
(2)对接收到的时间回波脉冲信号s'(t)依次进行下变频处理和采样率为fs的采样处理,得到基带脉冲s′r(n),其中,n为正整数,表示采样点数,基带脉冲的起始位置为No,基带脉冲总长度为Nsum,前M-1个脉冲单元为循环前缀,M为时间回波脉冲信号的距离单元数;
(3)对基带脉冲s′r(n)中的循环前缀进行移除,得到移除循环前缀后的信号sr(n):
(3a)计算基带脉冲s′r(n)中的提取位置:Ns=N0+M-1;
(3b)以提取位置Ns为起点,提取基带脉冲s′r(n)的后Nsum-Ns个单元,得到移除循环前缀后的信号sr(n);
(4)对移除循环前缀后的信号sr(n)进行傅里叶变换,得到基带频域信号Sr(k),其中k为正整数,表示脉冲索引;
(5)对基带频域信号Sr(k)进行通信处理,得到该基带频域信号中包含的通信信息
Figure FDA0002914305870000012
(5a)将基带频域信号Sr(k)向后延迟一个数据单元,并在空余位置处补零,得到延迟后的基带信号Sr′(k);
(5b)对延迟后的基带信号S′r(k)进行共轭运算,得到延迟共轭信号S″r(k);
(5c)将基带频域信号Sr(k)与延迟共轭信号S″r(k)的对应元素逐个进行乘法运算,得到解码结果序列Sc′;
(5d)对解码结果序列Sc′中的每个元素进行判决:当解码结果序列Sc′中元素的模值大于π/2时判定为1,小于或等于π/2时判定为0,得到全部元素的判定结果,即基带频域信号中包含的通信信息
Figure FDA0002914305870000011
(6)利用基带频域信号中包含的通信信息
Figure FDA0002914305870000021
构建参考信号sre(n):
(6a)对通信信息
Figure FDA0002914305870000022
按照下式进行编码运算,得到2DPSK信号
Figure FDA0002914305870000023
Figure FDA0002914305870000024
其中,
Figure FDA0002914305870000025
表示前一个通信信息
Figure FDA0002914305870000026
数据点,
Figure FDA0002914305870000027
表示前一个2DPSK信号数据点,当n-1=0时,
Figure FDA0002914305870000028
为0,
Figure FDA0002914305870000029
表示异或运算;
(6b)对2DPSK信号
Figure FDA00029143058700000210
按照下式进行相位调制,得到雷达处理所需的参考信号sre(n):
Figure FDA00029143058700000211
其中,j表示虚数单位;
(7)利用雷达处理所需的参考信号sre(n)以及移除循环前缀后的信号sr(n)进行脉冲压缩处理,得到脉冲压缩后信号为src(n)。
2.根据权利要求1所述的方法,其中(7)中对参考信号sre(n)进行脉冲压缩,其实现步骤如下:
(7a)对参考信号sre(n)进行共轭运算,得到参考信号sre(n)的频域共轭
Figure FDA00029143058700000212
(7b)将移除循环前缀后的信号sr(n)与参考信号sre(n)的频域共轭
Figure FDA00029143058700000213
相乘,并进行快速傅里叶逆变换,完成脉冲压缩处理。
CN201811218839.7A 2018-10-19 2018-10-19 Ofdm雷达通信一体化固定平台系统的信号处理方法 Active CN109085574B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811218839.7A CN109085574B (zh) 2018-10-19 2018-10-19 Ofdm雷达通信一体化固定平台系统的信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811218839.7A CN109085574B (zh) 2018-10-19 2018-10-19 Ofdm雷达通信一体化固定平台系统的信号处理方法

Publications (2)

Publication Number Publication Date
CN109085574A CN109085574A (zh) 2018-12-25
CN109085574B true CN109085574B (zh) 2021-03-23

Family

ID=64843816

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811218839.7A Active CN109085574B (zh) 2018-10-19 2018-10-19 Ofdm雷达通信一体化固定平台系统的信号处理方法

Country Status (1)

Country Link
CN (1) CN109085574B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110361699A (zh) * 2019-07-23 2019-10-22 北京工业大学 一种适用于南极航空观测现场的冰雷达数据处理的方法
CN115442197B (zh) * 2022-08-30 2024-02-27 西安电子科技大学 采用无循环前缀ofdm的一体化信号设计及处理方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101323681B1 (ko) * 2012-03-23 2013-10-30 국방과학연구소 Ofdm 레이더 시스템의 속도추정장치 및 속도추정방법
CN102901839B (zh) * 2012-11-05 2014-08-13 杭州电子科技大学 一种低复杂度的声学多普勒测流系统速度估算方法
CN105306399B (zh) * 2015-07-24 2019-02-19 西安电子科技大学 一种雷达通信一体化信号的优化方法
JP6716435B2 (ja) * 2016-11-21 2020-07-01 株式会社東芝 レーダシステム及びそのレーダ信号処理方法
CN107086975B (zh) * 2017-03-21 2020-06-12 西安电子科技大学 一种实现ofdm雷达通信信号一体化的方法和多载波系统

Also Published As

Publication number Publication date
CN109085574A (zh) 2018-12-25

Similar Documents

Publication Publication Date Title
CN108768446B (zh) 低截获雷达通信一体化系统的信号波形设计方法
CN109061633B (zh) Ofdm雷达通信一体化机载平台系统的信号设计方法
CN105891789B (zh) 联合时频分布和压缩感知对抗雷达频谱弥散干扰的方法
CN102323575B (zh) Pd雷达在微弱信号检测过程中的距离走动校正方法
CN107247254B (zh) 一种非线性调频信号去斜处理方法
US8928524B1 (en) Method and system for enhancing data rates
CN109085574B (zh) Ofdm雷达通信一体化固定平台系统的信号处理方法
CN110824438A (zh) 用于处理ofdm雷达信号的方法和装置
CN107872279B (zh) 基于正交频率分解的雷达通信共享信号设计方法
CN107490785A (zh) 一种基于频域置零调制的ofdm雷达通信一体化波形设计方法
US20160018512A1 (en) Method for Generating and Compressing Multi-Sweep-Frequency Radar Signals
CN109061634B (zh) Ofdm雷达通信一体化小型无人机系统的信号设计方法
CN107037409B (zh) 基于压缩感知的mimo雷达波形分离方法
CN109061580B (zh) 一种调频间断连续波雷达缓和相互干扰方法
CN110045375A (zh) 一种一比特sar回波数据的处理方法及系统
WO2018014734A1 (zh) 一种适用于OvXDM系统的调制方法、装置及OvXDM系统
CN109870684A (zh) 一种基于cp-ofdm的碎片频谱背景下雷达距离像重构方法
CN104062641A (zh) 一种ofdm雷达距离像处理的新方法
CN103576151A (zh) 基于压缩感知的方位向多通道sar成像方法及系统
CN109343056A (zh) 一种针对非线性调频sar的rd成像方法和装置
CN112014807A (zh) 一种频率捷变雷达的自适应杂波抑制方法
CN109085575B (zh) Ofdm雷达通信一体化机载平台系统的信号处理方法
Aberman et al. Adaptive frequency allocation in radar imaging: Towards cognitive SAR
CN108983191B (zh) Ofdm雷达通信一体化系统的低速信号处理方法
CN104345297A (zh) 一种基于频域调相低旁瓣混沌雷达信号的产生方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Liu Gaogao

Inventor after: Wang Zhao

Inventor after: Cai Jingjing

Inventor after: Qin Guodong

Inventor after: Wu Bin

Inventor after: Bao Dan

Inventor after: Li Peng

Inventor before: Liu Gaogao

Inventor before: Cai Jingjing

Inventor before: Qin Guodong

Inventor before: Wu Bin

Inventor before: Bao Dan

Inventor before: Li Peng

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210903

Address after: 710000 5th floor, block B, productivity building, No.3, electronic West Street, electronic industrial park, high tech Zone, Xi'an City, Shaanxi Province

Patentee after: XI'AN SENSING TECHNOLOGY DEVELOPMENT Co.,Ltd.

Address before: 710071 Taibai South Road, Yanta District, Xi'an, Shaanxi Province, No. 2

Patentee before: XIDIAN University

TR01 Transfer of patent right
CP01 Change in the name or title of a patent holder

Address after: 710000 5th floor, block B, productivity building, No.3, electronic West Street, electronic industrial park, high tech Zone, Xi'an City, Shaanxi Province

Patentee after: Xi'an Shengxin Technology Co.,Ltd.

Address before: 710000 5th floor, block B, productivity building, No.3, electronic West Street, electronic industrial park, high tech Zone, Xi'an City, Shaanxi Province

Patentee before: XI'AN SENSING TECHNOLOGY DEVELOPMENT CO.,LTD.

CP01 Change in the name or title of a patent holder
CP02 Change in the address of a patent holder

Address after: Room 810, Building C, 8th Floor, Chuangye Building, No. 16 Gaoxin 1st Road, Xi'an City, Shaanxi Province, 710065

Patentee after: Xi'an Shengxin Technology Co.,Ltd.

Address before: 710000 5th floor, block B, productivity building, No.3, electronic West Street, electronic industrial park, high tech Zone, Xi'an City, Shaanxi Province

Patentee before: Xi'an Shengxin Technology Co.,Ltd.

CP02 Change in the address of a patent holder