CN108512797A - 一种基于正交频分复用的雷达通信一体化信号设计方法 - Google Patents

一种基于正交频分复用的雷达通信一体化信号设计方法 Download PDF

Info

Publication number
CN108512797A
CN108512797A CN201810232683.1A CN201810232683A CN108512797A CN 108512797 A CN108512797 A CN 108512797A CN 201810232683 A CN201810232683 A CN 201810232683A CN 108512797 A CN108512797 A CN 108512797A
Authority
CN
China
Prior art keywords
signal
sequence
radar
indicate
frequency domain
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.)
Granted
Application number
CN201810232683.1A
Other languages
English (en)
Other versions
CN108512797B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810232683.1A priority Critical patent/CN108512797B/zh
Publication of CN108512797A publication Critical patent/CN108512797A/zh
Application granted granted Critical
Publication of CN108512797B publication Critical patent/CN108512797B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2614Peak power aspects
    • H04L27/2615Reduction thereof using coding
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2602Signal structure

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

该发明公开了一种基于正交频分复用的雷达通信一体化信号设计方法,属于雷达技术领域,涉及雷达通信一体化相关技术。本发明的解决方案是采用任意脉冲长度的OFDM波形作为雷达通信一体化信号设计的基础,首先利用时域限幅实现脉冲长度的任意化并限制信号的峰均值功率比,然后利用频域滤波与频域限幅保证信号的脉冲压缩性能,最后利用间隔式相位压缩法嵌入通信信息。该算法有效解决了传统雷达通信一体化设计方法中通信信息嵌入大幅降低OFDM信号的峰均值功率比性能与脉冲压缩效果的问题,从而实现了雷达系统与通信系统共用同一发射信号。

Description

一种基于正交频分复用的雷达通信一体化信号设计方法
技术领域
本发明属于雷达技术领域,涉及雷达通信一体化相关技术。
背景技术
雷达通信一体化即最大程度地实现雷达系统和通信系统在硬件资源上的共享。雷达系统与通信系统的一体化能够有效地节省硬件资源,减小设备体积,是实现电子战系统一体化与装备小型化的关键技术之一。而能否实现雷达与通信信号的一体化关系到能否最大程度上实现两个系统在硬件资源上的共享。OFDM信号能够有效地抑制距离分辨单元间回波能量的串扰,极有可能被用作下一代高精度雷达的发射信号。另外,基于OFDM信号的无线通信技术已经很成熟且应用广泛,因而基于OFDM波形的雷达通信一体化信号设计目前已经成为国内外专家的重点研究课题。
对于OFDM信号,信号的峰均值功率比是一项重要的性能评价指标。而脉冲压缩性能是雷达信号的一项重要的性能评价指标。雷达通信一体化信号必须要兼顾以上两个性能评价指标,并在此基础上嵌入通信功能。但传统的雷达通信一体化设计方法中通信信息的嵌入往往会大幅降低OFDM信号的峰均值功率比性能与脉冲压缩性能,因此基于OFDM波形的雷达通信一体化信号设计难以实现。在文献“OFDM synthetic aperture radar imagingwith sufficient cyclic prefix,IEEE Transactions on Geoscience and RemoteSensing,2015,53(1):394-404”中,提出了一种循环前缀为零的OFDM信号,实现了发射信号脉冲长度的任意化,利用这种信号可以实现宽带高精度合成孔径成像。但是该文献未考虑该波形在通信方面的应用。在现代电子战中雷达间的通信必不可少,因此国内外研究者关于OFDM雷达通信一体化波形进行了大量研究。在文献“A dual-function MIMO radar-communications system using PSK modulation,Signal Processing Conference(EUSIPCO),2016 24th European.IEEE,2016:1613-1617”中,利用传统的相位调制在OFDM脉冲中嵌入通信信息,但是该文献没有实现OFDM脉冲长度的任意化,也没有考虑OFDM信号峰均值功率比。以上文献都没有考虑基于任意脉冲长度OFDM波形来进行雷达通信一体化波形的设计,且没有充分考虑嵌入通信信息对信号性能带来的影响,因而在实际的雷达通信一体化系统中难以应用。
发明内容
本发明的目的是针对背景技术存在的缺陷,研究设计一种基于任意脉冲长度OFDM波形的雷达通信一体化波形设计算法,解决现有雷达通信一体化波形设计算法通信信息嵌入大幅影响OFDM信号峰均值功率比与脉冲压缩效果的问题。
本发明的解决方案是采用任意脉冲长度的OFDM波形作为雷达通信一体化信号设计的基础,首先利用时域限幅实现脉冲长度的任意化并限制信号的峰均值功率比,然后利用频域滤波与频域限幅保证信号的脉冲压缩性能,最后利用间隔式相位压缩法嵌入通信信息。该算法有效解决了传统雷达通信一体化设计方法中通信信息嵌入大幅降低OFDM信号的峰均值功率比性能与脉冲压缩效果的问题,从而实现了雷达系统与通信系统共用同一发射信号。
为了方便描述本发明的内容,首先对以下术语进行解释:
术语1:恒模性
恒模性指一个序列的所有元素都具有相同的模值。
术语2:信噪比保持因子
在雷达应用中,如果OFDM信号的频域复加权系数不具有恒模性,那么在进行脉冲压缩的过程中噪声会被放大,接收信号的信噪比会下降。用当前信号的信噪比与恒模信号的信噪比的比值来评价信噪比的下降的程度,记为信噪比保持因子。
术语3:间隔
指的是对OFDM信号的频域复加权系数进行调制时,相邻的两个被调制的频域复加权系数之间的间隔。
术语4:压缩区间
一个压缩区间对应于传统PSK调制星座图中的一个星座点,调制后的相位位于对应压缩区间内。
术语5:压缩区间数目
对应于传统PSK调制星座图中星座点的数目,类比传统的PSK调制方法,拥有两个压缩区间的间隔式相位压缩法简称为BPC,拥有四个个压缩区间的间隔式相位压缩法简称为QPC,以此类推。
本发明提出了一种基于正交频分复用的雷达通信一体化信号设计算法,该方法包括:
步骤1:初始化
初始化步骤2到步骤4需要循环迭代的次数Q、迭代的初始值,假设OFDM信号含有N个子载波,不同子载波的频域复加权系数序列可以表示为:
S=[S0,S1,...,SN]T (1)
其中,Sk表示第k个子载波的频域复加权系数,[·]T为向量转置运算;则OFDM信号的时域表达式为:
其中,Δf=B/N=1/Tz表示相邻子载波的频率间隔,B表示信号的频宽,Tz表示OFDM信号一个周期内的时宽,TCP表示循环前缀的时间长度;式(2)中的信号进行采样,采样间隔Ts=1/B=Tz/N,得到时域离散序列:
其中,sn表示第n个采样点的采样值,M表示距离分辨单元的数目的最大值,TCP=(M-1)Ts为循环前缀的时间长度,且N≥M;
s=[s0,s1,...,sM-2,sM,...,sN+M-2]T可以由S=[S0,S1,...,SN]T经过N点FFT变换得到,所以频域复加权系数序列S与时域离散序列s等价,将交替对这两个序列进行约束,以作为算法迭代的初始值;
步骤2:过采样及时域滤波与限幅
在时域应该控制发射信号的峰均值功率比尽量小,为了尽量精确地估计信号的峰均值功率比,需要对时域信号进行L倍过采样,L为过采样倍数;时域信号的过采样可以通过对频域复加权系数序列补零实现;因此首先将复加权系数序列S(q)尾部扩展(L-1)N个零,得到序列q表示第q次迭代;对序列做L×N点IFFT变换得到时域过采样序列
于是,OFDM信号的峰均值功率比的估计值可以表示为:
为了实现发射脉冲长度的任意化,循环前缀的幅度要为零;因此需要对时域过采样序列进行时域滤波,得到新的序列为:
接着,给OFDM信号的峰均值功率比加入相应的约束,对时域离散序列进行限幅;将幅度过大的采样点的幅度限制在均值附近,限幅的控制因子为PAPRd,使得限幅后时域离散序列中的元素的幅度应该满足:
其中,为序列所有非零元素的平均能量;接着,对时域序列进行LN点FFT变换得到频域序列
步骤3:频域滤波与限幅
由于在时域进行了滤波及限幅会造成信号能量辐射到OFDM信号的频带外,为了抑制带外能量,需要对频域序列进行滤波,得到新的序列为:
截取序列的前N个元素,为了保证接收信号的信噪比,需要给这N个元素的幅度的抖动加入相应的约束;
信噪比保持因子可以表示为:
其中,SNRm表示脉冲压缩后第m个距离分辨单元的信噪比,SNRmax,m表示恒模信号脉冲压缩后第m个距离分辨单元的信噪比;
接着对序列的前N个元素进行限幅,限幅的控制因子为Gf,得到序列 中的元素应该满足:
其中,为序列中元素的平均能量;
步骤4:通信信息的嵌入
对序列进行间隔式线性相位压缩从而嵌入通信信息;首先,按照间隔γ选出待调制的元素:
其中,F(·)表示调制(线性相位压缩)函数,ε=0,1,...,(N-1)/γ;
接下来,进行线性相位压缩;线性相位压缩一种改进的相位调制方式,改进之处在于它放宽了PSK调制对相位的约束;当对应于第j个压缩区间的通信信息被发送时,线性相位压缩法可以表示为:
其中,的相位,的相位;[ej,ej+σ]表示第j个压缩区间,假设压缩区间的数目为J,则满足ej+1-ej=ej-ej-1=2π/J;σ为区间大小;
然后,返回步骤2反复迭代,直到达到预先设置的迭代次数;
步骤5:规范化处理,
经过Q次迭代,得到频域复加权系数S(Q);由于在步骤3中进行了频域滤波与限幅,时域序列受到一定影响,因此需要进行最后的时域滤波处理以保证信号满足循环前缀幅度为零;虽然最后的时域滤波会使得一小部分能量泄露到OFDM信号的频谱外,但通过多次仿真发现多次迭代后泄露出的能量非常少,可以忽略不计;然后将得到的频域复加权序列进行归一化处理;
接着,计算最后得到的OFDM信号的峰均值功率比及信噪比保持因子λ,并与预先设置的门限值PAPR0及λ0做比较,若峰均值功率比小于PAPR0且信噪比保持因子λ大于λ0则输出得到的OFDM信号作为发射信号;若不满足回到步骤1重新进行设计;
步骤6:雷达接收端处理
设雷达通信装置位于空中移动平台,在进行地面探测的同时,与地面站点进行通信;则雷达接收端接收到的信号为:
其中,M1表示M1个距离分辨单元,对应于无线通信中的M1条径;fc表示载波频率,gm表示第m个距离分辨单元的雷达反射截面积,为空中移动平台与第m个距离分辨单元之间的距离,c为光速,w(t)表示噪声,rm(t)表示第m个距离分辨单元的回波;
对式(12)中的信号进行采样,接着进行FFT变换,得到频域离散信号为:
其中,[α01,...,αN-1]T为RCS序列的傅里叶变换,[S0,S1,...,SN-1]T为发射信号的频域复加权系数,Wk表示噪声;由于雷达发射端和接收端位于同一平台[S0,S1,...,SN-1]T已知,直接计算出[α01,...,αN-1]T
步骤7:通信接收端处理,
位于地面站点的通信接收端接收到的信号为:
其中,M2表示无线信道包含M2条经;为第条径对应的信道参数,为地面站点与引起第条径的反射单元间的距离;
由式(12)、(14)可知,距离分辨单元数目的最大值M满足:
M=max{M1,M2} (15)
对式(14)中的信号进行采样,接着进行FFT变换,得到频域离散信号为:
其中,[β01,...,βN-1]T为序列的傅里叶变换,通过信道估计得到,再计算出[S0,S1,...,SN-1]T
接着,解调取出通信信息;如果:
则接收到的为信息为第j个压缩区域对应的通信信息;其中,的相位,为接收端得到的带有噪声的频域复加权系数。
本发明的有益效果:本发明采用任意长度的OFDM波形作为模板信号,利用时域限幅限制信号的峰均值功率比,然后利用频域滤波与频域限幅保证信号的脉冲压缩性能,最后利用间隔式相位压缩法嵌入通信信息。该算法有效解决了传统雷达通信一体化信号设计方法中通信信息嵌入大幅降低OFDM信号峰均值功率比性能与脉冲压缩效果的问题。本发明的优点兼顾了信号的峰均值功率比、脉冲压缩性能及通信性能,实现了雷达与通信波形的一体化。本发明可以应用于民用军事等领域。
附图说明
图1是本发明提供方法流程框图。
图2是本发明提供方法的典型应用场景。
图3是本发明具体实施方式中采用不同的间隔式线性相位压缩参数时OFDM信号的峰均值功率比的仿真结果。
图4是本发明具体实施方式中采用不同的间隔式线性相位压缩参数时OFDM信号的信噪比保持因子λ性能的仿真结果。
图5是本发明具体实施方式中采用不同的间隔式线性相位压缩参数时OFDM信号的误符号率性能的仿真结果。
图6是本发明具体实施方式中设计出的OFDM信号与同等条件下LFM信号抗距离分辨单元间回波能量串扰的性能的仿真结果。
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2016b上验证正确。下面就具体实施方式对本发明作进一步的详细描述。
步骤一:初始化,
本发明提出的设计算法中步骤二到步骤四需要循环迭代,因此需要初始化循环次数Q。另外,迭代的初始值也要进行初始化,假设OFDM信号含有N个子载波,不同子载波的频域复加权系数序列可以表示为:
S=[S0,S1,...,SN]T (18)
其中,Sk表示第k个子载波的频域复加权系数,[·]T为向量转置运算。则OFDM信号的时域表达式为:
其中,Δf=B/N=1/Tz表示相邻子载波的频率间隔,B表示信号的频宽,Tz表示OFDM信号一个周期内的时宽,TCP表示循环前缀的时间长度。式(18)中的信号进行采样,采样间隔Ts=1/B=Tz/N,得到时域离散序列:
其中,sn表示第n个采样点的采样值,M表示距离分辨单元的数目的最大值,TCP=(M-1)Ts为循环前缀的时间长度,,且N≥M。
s=[s0,s1,...,sM-2,sM,...,sN+M-2]T可以由S=[S0,S1,...,SN]T经过N点FFT变换得到,所以频域复加权系数序列S与时域离散序列s等价,我们将交替对这两个序列进行约束和设计。为了方便,我们以作为算法迭代的初始值。
步骤二:过采样及时域滤波与限幅,
在时域我们应该控制发射信号的峰均值功率比尽量小,为了尽量精确地估计信号的峰均值功率比,需要对时域信号进行L倍过采样,L表示过采样的倍数。时域信号的过采样可以通过对频域复加权系数序列补零实现。因此,我们首先将复加权系数序列S(q)尾部扩展(L-1)N个零,得到序列q表示第q次迭代。对序列做L×N点IFFT变换得到时域过采样序列
于是,OFDM信号的峰均值功率比的估计值可以表示为:
为了实现发射脉冲长度的任意化,循环前缀的幅度要为零。因此,我们需要对时域过采样序列进行时域滤波,得到新的序列为:
接着,给OFDM信号的峰均值功率比加入相应的约束,对时域离散序列进行限幅。将幅度过大的采样点的幅度限制在均值附近,限幅的控制因子为PAPRd,使得限幅后时域离散序列中的元素的幅度应该满足:
其中,为序列所有非零元素的平均能量。接着,对时域序列进行LN点FFT变换得到频域序列
步骤三:频域滤波与限幅,
由于在时域进行了滤波及限幅会造成信号能量辐射到OFDM信号的频带外,为了抑制带外能量,需要对频域序列进行滤波,得到新的序列为:
截取序列的前N个元素,为了保证接收信号的信噪比,需要给这N个元素的幅度的抖动加入相应的约束。
信噪比保持因子可以表示为:
其中,SNRm表示脉冲压缩后第m个距离分辨单元的信噪比,SNRmax,m表示恒模信号脉冲压缩后第m个距离分辨单元的信噪比。
接着对序列的前N个元素进行限幅,限幅的控制因子为Gf,得到序列 中的元素应该满足:
其中,为序列中元素的平均能量。
步骤四:通信信息的嵌入,
对序列进行间隔式线性相位压缩从而嵌入通信信息。首先,按照间隔γ选出待调制的元素:
其中,F(·)表示调制(线性相位压缩)函数,ε=0,1,...,(N-1)/γ。
接下来,进行线性相位压缩。线性相位压缩一种改进的相位调制方式,改进之处在于它放宽了PSK调制对相位的约束。当对应于第j个压缩区间的通信信息被发送时,线性相位压缩法可以表示为:
其中,的相位,的相位。[ej,ej+σ]表示第j个压缩区间,假设压缩区间的数目为J,则满足ej+1-ej=ej-ej-1=2π/J。σ为区间大小。
然后,返回步骤二反复迭代,直到达到预先设置的迭代次数。
步骤五:规范化处理,
经过Q次迭代,得到频域复加权系数S(Q)。由于在步骤三中进行了频域滤波与限幅,时域序列受到一定影响,因此需要进行最后的时域滤波处理以保证信号满足循环前缀幅度为零。虽然最后的时域滤波会使得一小部分能量泄露到OFDM信号的频谱外,但通过多次仿真发现多次迭代后泄露出的能量非常少,可以忽略不计。然后将得到的频域复加权序列进行归一化处理。
接着,计算最后得到的OFDM信号的峰均值功率比及信噪比保持因子λ,并与预先设置的门限值PAPR0及λ0做比较,若峰均值功率比小于PAPR0且信噪比保持因子λ大于λ0则输出得到的OFDM信号作为发射信号。若不满足回到步骤一重新进行设计。
步骤六:雷达接收端处理,
假设雷达通信装置位于空中移动平台,在进行地面探测的同时,与地面站点进行通信。则雷达接收端接收到的信号为:
其中,M1表示M1个距离分辨单元,对应于无线通信中的M1条径。fc表示载波频率,gm表示第m个距离分辨单元的雷达反射截面积,为空中移动平台与第m个距离分辨单元之间的距离,c为光速,w(t)表示噪声,rm(t)表示第m个距离分辨单元的回波。
对式(29)中的信号进行采样,接着进行FFT变换,得到频域离散信号为:
其中,[α01,...,αN-1]T为RCS序列的傅里叶变换,[S0,S1,...,SN-1]T为发射信号的频域复加权系数,Wk表示噪声。由于雷达发射端和接收端位于同一平台[S0,S1,...,SN-1]T已知,所以[α01,...,αN-1]T容易求出。
步骤七:通信接收端处理,
位于地面站点的通信接收端接收到的信号为:
其中,M2表示无线信道包含M2条经。为第条径对应的信道参数,为地面站点与引起第条径的反射单元间的距离。
由式(29)、(31)可知,距离分辨单元数目的最大值M满足:
M=max{M1,M2} (32)
对式(31)中的信号进行采样,接着进行FFT变换,得到频域离散信号为:
其中,[β01,...,βN-1]T为序列的傅里叶变换,可以通过信道估计得到,因此[S0,S1,...,SN-1]T容易求出。
接着,解调取出通信信息。如果:
则接收到的为信息为第j个压缩区域对应的通信信息。其中,的相位,为接收端得到的带有噪声的频域复加权系数。
本发明的效果通过以下仿真对比试验进一步说明:
仿真场景:假定距离分辨单元数量的最大值M=96,OFDM信号带宽为B=150MHz,子载波数量N=128,发射信号非零时域序列长度Nt=33,采样频率fs=150MHz。OFDM信号的过采样率取L=4,算法迭代次数Q=100,限幅控制因子PAPRd=1dB,Gf=0.05。通过进行5×105次独立的Monte Carlo实验,对不同条件下所设计出的OFDM波形的PAPR以及信噪比保持因子λ进行统计分析,给出相应的累计分布函数(CDF)。并和相同条件下(发射信号能量、信号带宽、信号脉宽)的LFM信号进行了抗距离分辨单元间回波能量串扰性能的对比。另外还对不同条件下所设计出的OFDM波形进行通信误符号率实验,给出它们在不同信噪比的高斯白噪声下的通信误符号率(在实验中,共有2.048×106个随机符号被传输)。每次试验中,算法迭代的初始值取其中θk是一个服从区间为[-π,π]的均匀分布的随机变量。
图3给出了采用不同的间隔式线性相位压缩参数时OFDM信号的峰均值功率比的仿真结果。由图3可以看出,OFDM信号的峰均值功率比随着间隔γ及压缩区间大小σ的增大而不断改善。当参数选择合适时,通信信息的嵌入对OFDM信号的峰均值功率比的影响很小。例如当参数γ=8时OFDM信号的峰均值功率比上升了不到0.5dB,能够满足发射系统对OFDM信号峰均值功率比的要求。
图4给出了采用不同的间隔式线性相位压缩参数时OFDM信号的信噪比保持因子λ性能仿真结果。由图4可以看出,OFDM信号的信噪比保持因子λ性能随着间隔γ及压缩区间大小σ的增大而不断改善。当参数选择合适时,通信信息的嵌入对OFDM信号的信噪比保持因子性能的影响很小。例如当参数γ=8时OFDM信号的信噪比保持因子下降了不到0.2dB,能够满足系统对于信噪比保持因子的要求。
图5给出了采用不同的间隔式线性相位压缩参数时OFDM信号的误符号率性能仿真结果。由图5可以看出,OFDM信号的误符号率随着压缩区间数目J的减小而不断改善。值得注意的是这里我们没有采用任何抗干扰编码及纠错技术,当参数选择合适时,误符号率能够满足无线通信系统的要求。另外,由于雷达发射信号的脉冲重复频率很大,且带宽较大,因此,通信数据率可以达到Mbps级别。
图6给出另外本发明设计出的OFDM信号与同等条件下LFM信号抗距离分辨单元间回波能量串扰的性能的仿真结果。由图6可以看出,与雷达系统中常用的LFM信号相比,利用本发明提出的算法设计出的OFDM信号的抗距离分辨单元间回波能量串扰的能力更强,也就是说该算法保留了OFDM信号在该方面的性能优势。
通过本发明的具体实施方式可以看出,该雷达通信算法能够满足雷达通信一体化系统的要求。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (1)

1.一种基于正交频分复用的雷达通信一体化信号设计方法,该方法包括:
步骤1:初始化
初始化步骤2到步骤4需要循环迭代的次数Q、迭代的初始值,假设OFDM信号含有N个子载波,不同子载波的频域复加权系数序列可以表示为:
S=[S0,S1,...,SN]T (1)
其中,Sk表示第k个子载波的频域复加权系数,[·]T为向量转置运算;则OFDM信号的时域表达式为:
其中,Δf=B/N=1/Tz表示相邻子载波的频率间隔,B表示信号的频宽,Tz表示OFDM信号一个周期内的时宽,TCP表示循环前缀的时间长度;式(2)中的信号进行采样,采样间隔Ts=1/B=Tz/N,得到时域离散序列:
其中,sn表示第n个采样点的采样值,M表示距离分辨单元的数目的最大值,TCP=(M-1)Ts为循环前缀的时间长度,且N≥M;
s=[s0,s1,...,sM-2,sM,...,sN+M-2]T可以由S=[S0,S1,...,SN]T经过N点FFT变换得到,所以频域复加权系数序列S与时域离散序列s等价,将交替对这两个序列进行约束,以作为算法迭代的初始值;
步骤2:过采样及时域滤波与限幅
首先将复加权系数序列S(q)尾部扩展(L-1)N个零,得到序列q表示第q次迭代;对序列做L×N点IFFT变换得到时域过采样序列
于是,OFDM信号的峰均值功率比的估计值可以表示为:
对时域过采样序列进行时域滤波,得到新的序列为:
接着,给OFDM信号的峰均值功率比加入相应的约束,对时域离散序列进行限幅;将幅度过大的采样点的幅度限制在均值附近,限幅的控制因子为PAPRd,使得限幅后时域离散序列中的元素的幅度应该满足:
其中,为序列所有非零元素的平均能量;接着,对时域序列进行LN点FFT变换得到频域序列
步骤3:频域滤波与限幅
由于在时域进行了滤波及限幅会造成信号能量辐射到OFDM信号的频带外,为了抑制带外能量,需要对频域序列进行滤波,得到新的序列为:
截取序列的前N个元素,为了保证接收信号的信噪比,需要给这N个元素的幅度的抖动加入相应的约束;
信噪比保持因子可以表示为:
其中,SNRm表示脉冲压缩后第m个距离分辨单元的信噪比,SNRmax,m表示恒模信号脉冲压缩后第m个距离分辨单元的信噪比;
接着对序列的前N个元素进行限幅,限幅的控制因子为Gf,得到序列 中的元素应该满足:
其中,为序列中元素的平均能量;
步骤4:通信信息的嵌入
对序列进行间隔式线性相位压缩从而嵌入通信信息;首先,按照间隔γ选出待调制的元素:
其中,F(·)表示调制(线性相位压缩)函数,ε=0,1,...,(N-1)/γ;
接下来,进行线性相位压缩;线性相位压缩一种改进的相位调制方式,改进之处在于它放宽了PSK调制对相位的约束;当对应于第j个压缩区间的通信信息被发送时,线性相位压缩法可以表示为:
其中,的相位,的相位;[ej,ej+σ]表示第j个压缩区间,假设压缩区间的数目为J,则满足ej+1-ej=ej-ej-1=2π/J;σ为区间大小;
然后,返回步骤2反复迭代,直到达到预先设置的迭代次数;
步骤5:规范化处理,
经过Q次迭代,得到频域复加权系数S(Q);进行最后的时域滤波处理以保证信号满足循环前缀幅度为零;然后将得到的频域复加权序列进行归一化处理;
接着,计算最后得到的OFDM信号的峰均值功率比及信噪比保持因子λ,并与预先设置的门限值PAPR0及λ0做比较,若峰均值功率比小于PAPR0且信噪比保持因子λ大于λ0则输出得到的OFDM信号作为发射信号;若不满足回到步骤1重新进行设计;
步骤6:雷达接收端处理
设雷达通信装置位于空中移动平台,在进行地面探测的同时,与地面站点进行通信;则雷达接收端接收到的信号为:
其中,M1表示M1个距离分辨单元,对应于无线通信中的M1条径;fc表示载波频率,gm表示第m个距离分辨单元的雷达反射截面积,为空中移动平台与第m个距离分辨单元之间的距离,c为光速,w(t)表示噪声,rm(t)表示第m个距离分辨单元的回波;
对式(12)中的信号进行采样,接着进行FFT变换,得到频域离散信号为:
其中,[α01,...,αN-1]T为RCS序列的傅里叶变换,[S0,S1,...,SN-1]T为发射信号的频域复加权系数,Wk表示噪声;由于雷达发射端和接收端位于同一平台[S0,S1,...,SN-1]T已知,直接计算出[α01,...,αN-1]T
步骤7:通信接收端处理,
位于地面站点的通信接收端接收到的信号为:
其中,M2表示无线信道包含M2条经;为第条径对应的信道参数,为地面站点与引起第条径的反射单元间的距离;
由式(12)、(14)可知,距离分辨单元数目的最大值M满足:
M=max{M1,M2} (15)
对式(14)中的信号进行采样,接着进行FFT变换,得到频域离散信号为:
其中,[β01,...,βN-1]T为序列的傅里叶变换,通过信道估计得到,再计算出[S0,S1,...,SN-1]T
接着,解调取出通信信息;如果:
则接收到的为信息为第j个压缩区域对应的通信信息;其中,的相位,为接收端得到的带有噪声的频域复加权系数。
CN201810232683.1A 2018-03-21 2018-03-21 一种基于正交频分复用的雷达通信一体化信号设计方法 Active CN108512797B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810232683.1A CN108512797B (zh) 2018-03-21 2018-03-21 一种基于正交频分复用的雷达通信一体化信号设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810232683.1A CN108512797B (zh) 2018-03-21 2018-03-21 一种基于正交频分复用的雷达通信一体化信号设计方法

Publications (2)

Publication Number Publication Date
CN108512797A true CN108512797A (zh) 2018-09-07
CN108512797B CN108512797B (zh) 2020-10-23

Family

ID=63377769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810232683.1A Active CN108512797B (zh) 2018-03-21 2018-03-21 一种基于正交频分复用的雷达通信一体化信号设计方法

Country Status (1)

Country Link
CN (1) CN108512797B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109507641A (zh) * 2018-11-07 2019-03-22 电子科技大学 一种合作的mimo雷达和mimo通信系统性能计算方法
CN109688082A (zh) * 2019-01-11 2019-04-26 电子科技大学 一种基于ofdm载波联合优化的通信雷达一体化系统
CN110324266A (zh) * 2019-05-16 2019-10-11 华南师范大学 一种自适应传输幅度位数压缩的ofdm可见光收发方法
CN111447158A (zh) * 2020-03-30 2020-07-24 电子科技大学 一种基于有限域傅里叶变换的信道编码特征提取方法
CN111836296A (zh) * 2020-06-29 2020-10-27 西南电子技术研究所(中国电子科技集团公司第十研究所) 一体化波形通信测量系统
CN112118013A (zh) * 2019-06-20 2020-12-22 南京大学 一种适用于多进制ldpc码的高性能、低复杂度的迭代大数逻辑译码方法
CN113067791A (zh) * 2021-03-22 2021-07-02 电子科技大学 一种基于相位与正交信号联合的信息调制和解调方法
CN113189548A (zh) * 2021-03-04 2021-07-30 西安电子科技大学 一种基于互信息的ofdm雷达通信一体化信号生成方法
CN113347132A (zh) * 2021-05-25 2021-09-03 西安空间无线电技术研究所 一种基于tds-ofdm的雷达通信一体化信号确定及传输方法
CN115242335A (zh) * 2022-07-19 2022-10-25 电子科技大学 一种ofdm雷达通信一体化信号联合优化设计与处理方法
CN115267686A (zh) * 2022-09-27 2022-11-01 烟台初心航空科技有限公司 具有低峰均比特性的雷达通信一体化调制信号产生方法
WO2023273781A1 (zh) * 2021-06-30 2023-01-05 华为技术有限公司 通信方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150244430A1 (en) * 2001-04-26 2015-08-27 Genghiscomm Holdings, LLC Cloud Radio Access Network
CN106911605A (zh) * 2017-03-01 2017-06-30 哈尔滨工业大学 一种基于fmcw的雷达通信一体化波形生成方法
CN107086975A (zh) * 2017-03-21 2017-08-22 西安电子科技大学 一种实现ofdm雷达通信信号一体化的方法
CN107490785A (zh) * 2017-09-19 2017-12-19 电子科技大学 一种基于频域置零调制的ofdm雷达通信一体化波形设计方法
CN107677995A (zh) * 2017-09-08 2018-02-09 电子科技大学 一种基于pmepr‑pslr联合优化的ofdm波形设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150244430A1 (en) * 2001-04-26 2015-08-27 Genghiscomm Holdings, LLC Cloud Radio Access Network
CN106911605A (zh) * 2017-03-01 2017-06-30 哈尔滨工业大学 一种基于fmcw的雷达通信一体化波形生成方法
CN107086975A (zh) * 2017-03-21 2017-08-22 西安电子科技大学 一种实现ofdm雷达通信信号一体化的方法
CN107677995A (zh) * 2017-09-08 2018-02-09 电子科技大学 一种基于pmepr‑pslr联合优化的ofdm波形设计方法
CN107490785A (zh) * 2017-09-19 2017-12-19 电子科技大学 一种基于频域置零调制的ofdm雷达通信一体化波形设计方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MUSSA A DIDA,ET AL: ""Constant envelope chirped OFDM for power-efficient radar communication"", 《2016 IEEE INFORMATION TECHNOLOGY, NETWORKING, ELECTRONIC AND AUTOMATION CONTROL CONFERENCE》 *
TIANXIAN ZHANG,ET AL: ""IRCI Free Range Reconstruction for SAR Imaging With Arbitrary Length OFDM Pulse"", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
TIANXIAN ZHANG,ET AL: ""OFDM Synthetic Aperture Radar Imaging With Sufficient Cyclic Prefix"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
李自琦等: ""基于分组格雷编码的OFDM雷达通信一体化系统峰均功率比抑制"", 《雷达学报》 *
王孟奇: ""基于OFDM的雷达通信一体化信号设计"", 《中国优秀硕士学位论文全文数据库(电子期刊)》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109507641A (zh) * 2018-11-07 2019-03-22 电子科技大学 一种合作的mimo雷达和mimo通信系统性能计算方法
CN109507641B (zh) * 2018-11-07 2023-05-02 电子科技大学 一种合作的mimo雷达和mimo通信系统性能计算方法
CN109688082A (zh) * 2019-01-11 2019-04-26 电子科技大学 一种基于ofdm载波联合优化的通信雷达一体化系统
CN109688082B (zh) * 2019-01-11 2021-02-12 电子科技大学 一种基于ofdm载波联合优化的通信雷达一体化系统
CN110324266A (zh) * 2019-05-16 2019-10-11 华南师范大学 一种自适应传输幅度位数压缩的ofdm可见光收发方法
CN110324266B (zh) * 2019-05-16 2022-03-11 华南师范大学 一种自适应传输幅度位数压缩的ofdm可见光收发方法及系统
CN112118013A (zh) * 2019-06-20 2020-12-22 南京大学 一种适用于多进制ldpc码的高性能、低复杂度的迭代大数逻辑译码方法
CN111447158A (zh) * 2020-03-30 2020-07-24 电子科技大学 一种基于有限域傅里叶变换的信道编码特征提取方法
CN111447158B (zh) * 2020-03-30 2021-06-29 电子科技大学 一种基于有限域傅里叶变换的信道编码特征提取方法
CN111836296A (zh) * 2020-06-29 2020-10-27 西南电子技术研究所(中国电子科技集团公司第十研究所) 一体化波形通信测量系统
CN111836296B (zh) * 2020-06-29 2022-07-29 西南电子技术研究所(中国电子科技集团公司第十研究所) 一体化波形通信测量系统
CN113189548A (zh) * 2021-03-04 2021-07-30 西安电子科技大学 一种基于互信息的ofdm雷达通信一体化信号生成方法
CN113189548B (zh) * 2021-03-04 2024-01-30 西安电子科技大学 一种基于互信息的ofdm雷达通信一体化信号生成方法
CN113067791A (zh) * 2021-03-22 2021-07-02 电子科技大学 一种基于相位与正交信号联合的信息调制和解调方法
CN113067791B (zh) * 2021-03-22 2022-07-15 电子科技大学 一种基于相位与正交信号联合的信息调制和解调方法
CN113347132A (zh) * 2021-05-25 2021-09-03 西安空间无线电技术研究所 一种基于tds-ofdm的雷达通信一体化信号确定及传输方法
WO2023273781A1 (zh) * 2021-06-30 2023-01-05 华为技术有限公司 通信方法及装置
CN115242335B (zh) * 2022-07-19 2023-09-22 电子科技大学 一种ofdm雷达通信一体化信号联合优化设计与处理方法
CN115242335A (zh) * 2022-07-19 2022-10-25 电子科技大学 一种ofdm雷达通信一体化信号联合优化设计与处理方法
CN115267686B (zh) * 2022-09-27 2023-01-17 烟台初心航空科技有限公司 具有低峰均比特性的雷达通信一体化调制信号产生方法
CN115267686A (zh) * 2022-09-27 2022-11-01 烟台初心航空科技有限公司 具有低峰均比特性的雷达通信一体化调制信号产生方法

Also Published As

Publication number Publication date
CN108512797B (zh) 2020-10-23

Similar Documents

Publication Publication Date Title
CN108512797A (zh) 一种基于正交频分复用的雷达通信一体化信号设计方法
CN104052691B (zh) 基于压缩感知的mimo‑ofdm系统信道估计方法
Nissel et al. Enabling low-complexity MIMO in FBMC-OQAM
CN109039960A (zh) 一种水声稀疏信道估计变步长稀疏度自适应匹配追踪方法
CN106453162B (zh) 针对多输入多输出正交频分复用系统的信道估计方法
Like et al. Signal classification in fading channels using cyclic spectral analysis
CN106534002B (zh) 一种基于压缩感知的电力线信道估计方法
US20100002816A1 (en) Method and apparatus for multiple signal identification and finding the basis functions of the received signal
CN109787932A (zh) 一种基于ofdm的雷达通信一体化信号设计方法
CN107359899B (zh) 稀疏信道条件下正交频分复用扩频水声通信无导频判决反馈信道估计方法
CN106936743B (zh) 一种电力线通信系统脉冲噪声抑制方法
CN108768903A (zh) 基于crdsa类协议的低轨道卫星随机接入系统多用户检测方法
CN105116384A (zh) 基于混沌序列的调频调相的雷达正交波形设计方法
CN109600327A (zh) 一种基于虚部干扰利用的信道估计方法
CN105187352A (zh) 一种基于ofdm前导的整数频偏估计方法
CN104794264A (zh) 一种基于稀疏频率的雷达通信波形设计方法
CN105791194A (zh) 一种可抗窄带干扰的同步方法及其系统
CN105959004B (zh) 一种基于大规模mimo的单精度adc自适应门限量化方法
CN109870684A (zh) 一种基于cp-ofdm的碎片频谱背景下雷达距离像重构方法
CN104539562A (zh) 基于多输入多输出正交频分复用的宽带短波信道估计方法
Kellett et al. Radar communications via random sequence encoding
Bai et al. Impulsive noise mitigation with interleaving based on MUSIC in power line communication
Anoh et al. Wavelet packet transform modulation for multiple input multiple output applications
Sejdic et al. Channel estimation using DPSS based frames
CN116125457A (zh) 基于otfs调制的雷达嵌入式通信方法

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
GR01 Patent grant
GR01 Patent grant