CN114124238B - 基于时分体制的otfs通信雷达一体化波形设计方法 - Google Patents

基于时分体制的otfs通信雷达一体化波形设计方法 Download PDF

Info

Publication number
CN114124238B
CN114124238B CN202111219969.4A CN202111219969A CN114124238B CN 114124238 B CN114124238 B CN 114124238B CN 202111219969 A CN202111219969 A CN 202111219969A CN 114124238 B CN114124238 B CN 114124238B
Authority
CN
China
Prior art keywords
radar
communication
signal
transmitter
target
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
CN202111219969.4A
Other languages
English (en)
Other versions
CN114124238A (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.)
Beijing University of Posts and Telecommunications
Original Assignee
Beijing University of Posts and Telecommunications
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 Beijing University of Posts and Telecommunications filed Critical Beijing University of Posts and Telecommunications
Priority to CN202111219969.4A priority Critical patent/CN114124238B/zh
Publication of CN114124238A publication Critical patent/CN114124238A/zh
Application granted granted Critical
Publication of CN114124238B publication Critical patent/CN114124238B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/50Transmitters
    • H04B10/516Details of coding or modulation
    • H04B10/548Phase or frequency modulation
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/282Transmitters
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/411Identification of targets based on measurements of radar reflectivity
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/50Transmitters
    • H04B10/516Details of coding or modulation
    • H04B10/548Phase or frequency modulation
    • H04B10/556Digital modulation, e.g. differential phase shift keying [DPSK] or frequency shift keying [FSK]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

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

Abstract

本发明涉及无线通信和雷达信号处理技术领域,提供一种基于时分体制的OTFS通信雷达一体化波形设计方法,具体包括:产生基于OTFS调制的通信雷达信号;对通信雷达信号进行功率控制,并根据功能需求分配合适的发射功率;为雷达探测和通信传输分配工作时隙,插入保护带;发射机处于通信传输工作时段,对通信传输信号OTFS解调还原;发射机处于雷达探测工作时段,对雷达回波OTFS解调判读与匹配滤波接收,进而获取探测信息;本发明使用OTFS信号能支持高速移动通信不受高多普勒环境的影响,满足通信雷达一体化系统多样化的功能需求,有效提升时间资源利用率。

Description

基于时分体制的OTFS通信雷达一体化波形设计方法
技术领域
本发明属于无线通信和雷达信号处理技术领域,具体涉及一种基于时分体制的OTFS通信雷达一体化波形设计方法。
背景技术
传统的雷达传感系统与无线通信系统分别按照各自的功能用途,以不同的工作频段独立设计,呈现为并行存在、纵向发展的状态,但随着5G通信、车联网和智能电子作战等新型民用和军事领域需求的井喷式增长,与雷达、通信相关的无线设备呈指数级增加,高速数据速率对带宽的需求爆发式增长,导致Sub-6GHz等电磁频谱受到严重的相互干扰与拥塞,同时频谱资源的利用率还很不充分,为缓解无线频谱资源紧缺的现状,将通信设备和雷达设备实现频段共享逐渐成为发展趋势,诸如将5G毫米波通信频段与车载毫米波雷达工作频段在系统上从硬件架构到信号处理进行方法趋同,就是一个很好的尝试;随着通信技术和雷达技术共享工作频段发展趋势的推进,横向上将两个相互独立的系统通过融合,打造能够满足通信和雷达要求的兼容复合性新系统,即通信雷达一体化已成为可能,其一体化技术的关键是设计能够同时满足通信功能和雷达需求的信号波形。
目前,基于正交频分复用(Orthogonal Frequency Division Multiplexing,OFDM)的通信雷达一体化波形设计是研究热点和主流体制,OFDM信号作为雷达信号具有近似图钉状的模糊函数,能实现距离域和多普勒域的解耦,OFDM信号作为通信波形具有频谱效率高、可对抗多径效应的多载波调制优势,但是,OFDM一体化信号的性能依赖于子载波间的正交性,当探测高速运动目标的多普勒频移对信号正交性影响较大时,将导致严重的载波间干扰(Inter-Carrier Interference,ICI),难以保证可靠的通信。为解决OFDM的上述问题,近年来提出了一种在时延-多普勒域调制信息的新型二维调制技术,即OTFS调制;基于时分体制的正交时间频率空间(Orthogonal Time Frequency Space,OTFS)调制信号支持高多普勒环境下的高速移动通信,允许在探测目标时存在更大的多普勒频移估计;但目前,基于OTFS的通信雷达一体化相比于基于OFDM的通信雷达一体化在对高速运动目标的通信和雷达探测方面,其比较优势的技术解决途径都集中在OTFS雷达探测目标的参数估计上,并没有形成完整的基于OTFS调制的通信雷达一体化方案。
发明内容
针对现有技术存在的问题,本发明提供一种基于时分体制的OTFS通信雷达一体化波形设计方法,具体包括以下步骤:
步骤1.对通信数据与雷达数据进行数字调制和OTFS调制,得到通信发射信号与雷达发射信号:
步骤1.1.通信信源根据待通信传输的内容产生串行比特通信信号流,雷达数据产生器按照伪随机PN序列生成雷达比特数据流;
步骤1.2.将串行比特通信信号流与雷达比特数据流通过星座映射进行数字调制,得到相对应的信息符号x[k,l],数字调制的方式包括正交相移键控(Quadrature PhaseShift Keying,QPSK)和正交振幅键控(Quadrature Amplitude Modulation,QAM);
步骤1.3.将信息符号x[k,l]放置于时延-多普勒域平面网格,进行OTFS调制;
步骤2.根据通信传输的覆盖能力和雷达探测的作用距离需求,发射机对雷达信号和通信信号进行功率控制,为雷达信号和通信信号分配合适的发射功率:
步骤2.1.由雷达方程确定雷达探测目标的作用距离与各雷达参数间的关系;
步骤2.2.由弗里斯传输公式确定通信传输的覆盖范围与各通信参数间的关系;
步骤2.3.发射机进行发射信号的功率控制,为通信信号分配满足通信覆盖能力的发射功率,为雷达信号分配满足雷达探测距离需求的发射功率;
步骤3.通过选通转换开关控制通信与雷达的分时工作,为雷达探测功能和通信传输功能固定工作时隙,并插入防止雷达信号与通信信号互相干扰的保护带,实现雷达探测与通信传输以时分复用的方式共享天线、发射机和接收机平台;
步骤4.发射机处于通信传输的工作时段,则执行本步骤,发送OTFS通信信号与目标接收机进行信息传输,否则执行步骤5:
步骤4.1.发送OTFS连续时间通信信号,并经过信道传输后到达接收端,生成接收信号y(t);
步骤4.2.接收端对接收信号y(t)进行OTFS解调;
步骤4.3.对时延-多普勒域符号y[k,l]进行检测,得到估计信号
Figure BDA0003312208580000022
并通过星座解映射获取对应星座点上的估计比特,进而恢复出原始发送的串行比特通信信号流;
步骤5.发射机处于雷达探测功能的工作时段,发送雷达信号进行目标探测,探测到目标后反射雷达回波至发射机,由发射机进行雷达信号处理:
步骤5.1.建立时延-多普勒域雷达信道h(τ,ν);
步骤5.2.发射的雷达信号经过时延-多普勒域雷达信道h(τ,ν),在探测到目标后反射雷达回波至发射机,生成接收信号y(t),发射机对y(t)进行OTFS解调,应用魏格纳变换与SFFT得到时延-多普勒域上的接收数据符号y[k,l];
步骤5.3.对得到时延-多普勒域上的接收数据符号y[k,l]通过雷达匹配滤波接收处理,进而获取探测目标对应于相对距离的时延特性和对应于相对运动速度的多普勒特性;
进一步的,所述步骤1.3.中时延-多普勒域平面网格定义为Г,并满足下式(1):
Figure BDA0003312208580000021
上式(1)中,k为多普勒域的索引,l为时延域的索引,1/(NT)表示多普勒轴的分辨率,1/(MΔf)表示时延轴的分辨率,(k)N表示对k取平移运算,并满足下式(2):
Figure BDA0003312208580000031
上式(2)中,目标运动方向将导致多普勒频移取值可正可负,多普勒域索引k经过平移得到的(k)N的数值为正表示目标朝着靠近发射机的方向运动,数值为负则表示目标朝着远离发射机的方向运动。
进一步的,所述步骤1.3.将信息符号x[k,l]放置于时延-多普勒域平面网格,进行OTFS调制,包括以下步骤:
步骤1.3.1.使用辛有限傅里叶逆变换(Inverse Symplectic Finite FourierTransform,ISFFT)将时延-多普勒域中的信息符号x[k,l]映射为时间-频率域上的符号X[n,m],映射过程满足下式(3):
Figure BDA0003312208580000032
上式(3)中,X[n,m]表示时间-频率域网格上的信号,x[k,l]表示时延-多普勒域网格上的信号,N代表一帧中时域上OTFS符号的数量,M代表频率域上子载波的数量;
步骤1.3.2.对时间-频率域平面网格上的符号X[m,n]应用海森堡变换(Heisenberg Transform),使用发射脉冲成形函数获得时间域上准备发射的连续时间的通信信号或者雷达信号s(t),连续时间的通信信号或者雷达信号s(t)变换公式为下式(4):
Figure BDA0003312208580000033
上式(4)中,gtx(t)为发射矩阵脉冲成形滤波器的冲激响应,t代表当前时刻。
进一步的,所述步骤1.3.2.中时间-频率域平面网格定义为Λ,并满足下式(5):
Λ={(nT,mΔf),n=0,1,...,N-1,m=0,1,...,M-1}…(5),
上式(5)中,n为时间的索引,m为频率的索引,T是一帧中一个OTFS符号的持续时间,Δf是相邻子载波间的频率间隔,且Δf=1/T。
进一步的,所述步骤2.1.中,由雷达方程得到下式(6):
Figure BDA0003312208580000041
上式(6)中,Pr,rad表示发射机接收到的目标回波功率,Pt,rad表示雷达信号的发射功率,Gt表示发射天线的增益,λ表示工作的波长,且探测信号所用波长λ乘以工作的载波频率fc等于光速,σ表示目标的散射截面积,Rrad为雷达探测目标的作用距离;
由上式(6)得知,当雷达发射机接收到的回波功率等于发射机的最小可检测雷达信号功率时,即得到雷达检测目标的最大作用距离Rmax,rad,如下式(7):
Figure BDA0003312208580000042
上式(7)中,Smin,rad表示发射机的最小可检测雷达信号功率。
进一步的,所述步骤2.2.中,由弗里斯传输公式得到下式(8):
Figure BDA0003312208580000043
上式(8)中,Pr,com表示目标接收机接收到的通信信号的功率,Pt,com表示通信信号的发射功率,Gr表示目标接收机的接收天线的增益,Gt表示发射天线的增益,Rcom为发射机与目标通信的通信覆盖半径;
由公式(8)得知,当目标接收机接收到的通信信号的功率等于目标接收机的通信信号接收功率阈值时,获得发射机与目标接收机正常通信的最大通信覆盖半径Rmax,com,如下式(9):
Figure BDA0003312208580000044
上式(9)中,Smin,com表示目标接收机的通信信号接收功率阈值。
进一步的,所述步骤2.3.中,发射机进行发射信号的功率控制,具体包括以下三种情况:
当雷达探测距离和通信传输覆盖的需求为雷达探测的最大作用距离大于通信传输的最大通信覆盖半径时,则发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com满足下式(10):
Figure BDA0003312208580000045
当雷达探测距离和通信传输覆盖的需求为雷达最大作用距离等于最大通信覆盖半径时,发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com满足下式(11):
Figure BDA0003312208580000051
当雷达探测距离和通信传输覆盖的需求为雷达最大作用距离小于最大通信覆盖半径时,发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com满足下式(12):
Figure BDA0003312208580000052
进一步的,所述步骤3.中为通信信号和雷达信号分配固定的工作时隙,并插入防止雷达信号与通信信号互相干扰的保护带,具体方法如下:
设定发射机在T1时间段内发射OTFS雷达信号执行探测功能,设定发射机在T2时间段内发射OTFS通信信号执行通信功能,则T1和T2之间插入的保护带Tg满足下式(13):
Figure BDA0003312208580000053
上式(13)中,τmax表示雷达探测时发射机与待测目标的最大往返时延,Rmax,rad表示雷达探测目标的最大作用距离,c表示光速。
进一步的,所述步骤4.1.中,OTFS连续时间通信信号经过信道传输后到达接收端并生成接收信号y(t)的过程满足下式(14):
y(t)=∫∫h(τ,ν)s(t-τ)ej2πν(t-τ)dνdτ…(14),
上式(14)中,s(t)表示发送的通信信号,τ表示时延,ν表示多普勒频移,h(τ,ν)为时延-多普勒信道。
进一步的,所述步骤4.2.接收端对接收信号y(t)进行OTFS解调,包括以下步骤:
步骤4.2.1.应用魏格纳变换(Wigner Transform)将接收信号y(t)转换为时间-频率域信号Y[n,m];
魏格纳变换的公式为下式(15):
Figure BDA0003312208580000054
上式(15)中,(·)*表示共轭运算,grx(t)表示接收脉冲成形函数,Y(t,f)为grx(t)与y(t)的互模糊函数;
以间隔t=nT和f=mΔf进行采样,得到时间-频率域信号Y[n,m]如下式(16):
Figure BDA0003312208580000066
步骤4.2.2.对时间-频率域信号Y[n,m]使用辛有限傅里叶变换(SymplecticFinite Fourier Transform,SFFT),将时间-频率域信号Y[n,m]映射为时延-多普勒域上符号y[k,l],映射过程满足下式(17):
Figure BDA0003312208580000061
上式(17)中,N代表时域上OTFS符号的数量,M代表频率域上子载波的数量。
进一步的,所述步骤5.1.中,时延-多普勒域雷达信道h(τ,ν)的建模表达式为下式(18):
Figure BDA0003312208580000062
上式(18)中,hp表示第p(p=1,...,P)个目标的复信道增益,τp表示发射机与第p个目标的双程往返时延,νp表示第p个目标的双程多普勒频移;
其中,τp满足下式(19):
Figure BDA0003312208580000063
上式(19)中,Rp(m)表示第p个目标与发射机的相对距离;
νp满足下式(20):
Figure BDA0003312208580000064
上式(20)中,Vp(m/s)表示第p个目标与发射机的相对运动速度,fc为工作的载波频率;
由式(18)~(20)得到离散形式的时延-多普勒域雷达信道h(τ,ν),如下式(21):
Figure BDA0003312208580000065
上式(21)中,h[k,l]表示多普勒抽头k和时延抽头l处的复信道增益,δ(·)为狄利克雷函数。
进一步的,所述步骤5.2.中,时延-多普勒域上的接收数据符号y[k,l]与发送端信息符号x[k,l]间的输入输出关系满足下式(22):
Figure BDA0003312208580000071
上式(22)中,<·>N表示模N运算,<·>M表示模M运算,w[k,l]为噪声矩阵,并且满足w[k,l]~CN(0,σ2),表示w[k,l]为服从均值为零、方差为σ2的复高斯分布的加性高斯白噪声,α[k,l]表示相移;
α[k,l]的表达式如下式(23):
Figure BDA0003312208580000072
进一步的,所述步骤5.3.中,对得到的接收数据符号y[k,l]进行雷达处理,具体包括如下步骤:
步骤5.3.1.分别将N×M维接收数据信号符号矩阵y[k,l]、发送信号符号矩阵x[k,l]、噪声矩阵w[k,l]和雷达信道增益矩阵h[k,l](0≤k≤N-1,0≤l≤M-1)按列优先的顺序进行排列,得到MN×1维接收符号列矢量
Figure BDA0003312208580000073
发送符号列矢量/>
Figure BDA0003312208580000074
噪声列矢量/>
Figure BDA0003312208580000075
和雷达信道列矢量/>
Figure BDA0003312208580000076
其中,y[k,l]、x[k,l]、w[k,l]和h[k,l]分别为
Figure BDA0003312208580000077
的第k+Nl个元素(0≤k+Nl≤MN-1),将各矢量之间关系用矩阵与矢量形式表达为下式(24):
Figure BDA0003312208580000078
上式(24)中,
Figure BDA0003312208580000079
的第一列为/>
Figure BDA00033122085800000710
其余列为/>
Figure BDA00033122085800000711
的循环移位;将公式(24)中MN×MN维矩阵/>
Figure BDA00033122085800000712
的第i行第j列元素表示为下式(25):
Figure BDA00033122085800000713
上式(25)中,i满足0≤i=k′+Nl′≤MN-1,j满足0≤j=k″+Nl″≤MN-1;
步骤5.3.2.对得到的接收数据符号y[k,l]进行匹配滤波处理,处理过程为下式(26):
Figure BDA0003312208580000081
上式(26)中,(·)H表示共轭转置,
Figure BDA0003312208580000082
为匹配滤波处理估计的雷达信道矢量,第k+Nl个元素为多普勒抽头k和时延抽头l处的复信道增益/>
Figure BDA0003312208580000083
非零元素所在位置代表探测到的目标;
若当
Figure BDA0003312208580000084
时,/>
Figure BDA0003312208580000085
不为0,表示在多普勒抽头/>
Figure BDA0003312208580000086
和时延抽头/>
Figure BDA0003312208580000087
位置处探测到第p个目标,目标相应的多普勒频移为/>
Figure BDA0003312208580000088
对应的目标与发射机的相对运动速度/>
Figure BDA0003312208580000089
为下式(27):
Figure BDA00033122085800000810
目标相应的时延为
Figure BDA00033122085800000811
对应的目标与发射机的相对距离/>
Figure BDA00033122085800000812
为下式(28):
Figure BDA00033122085800000813
与现有技术相比,本发明具有的显著的有益效果在于:
1.本发明所述的基于时分体制的OTFS通信雷达一体化波形设计方法,使用OTFS信号支持高速移动通信不受高多普勒环境的影响,而在探测目标时允许更大的多普勒频移估计,并以适当的复杂度为代价实现了对多个目标的探测。
2.本发明所述的基于时分体制的OTFS通信雷达一体化波形设计方法,引入了功率控制功能,根据特定的通信传输的覆盖能力和雷达探测的作用距离需求,为通信信号和雷达信号分配合适的发射功率,从而满足通信雷达一体化系统多样化的功能需求。
3.本发明所述的基于时分体制的OTFS通信雷达一体化波形设计方法,基于时分复用的方式控制通信传输与雷达探测功能共享硬件和分时工作,并插入保护带防止雷达信号与通信信号相互干扰,操作简便实用,同时灵活调整时隙的分配,有效提升时间资源利用率,为通信时段内加入跳频图案以增加保密性和抗干扰性提供了技术支撑。
附图说明
图1是本发明所述基于时分体制的OTFS通信雷达一体化波形设计方法的流程图;
图2是本发明所述基于时分体制的OTFS通信雷达一体化波形设计方法的原理框图;
图3是本发明所述步骤3时隙分配的方法示意图。
具体实施方式
下面结合说明书附图,进一步详细描述本发明的具体实施方式。
如图1所示,一种基于时分体制的OTFS通信雷达一体化波形设计方法,包括如下步骤:
步骤1.对通信数据或雷达数据进行数字调制和OTFS调制,得到通信发射信号或雷达发射信号;
步骤1.1.通信信源根据待通信传输的内容产生串行比特通信信号流,雷达数据产生器按照伪随机PN序列生成雷达比特数据流;
步骤1.2.将串行比特通信信号流或者雷达比特数据流通过星座映射使用正交相移键控(Quadrature Phase Shift Keying,QPSK)或者正交振幅键控(QuadratureAmplitude Modulation,QAM)进行数字调制,得到相对应的信息符号x[k,l];
步骤1.3.将信息符号x[k,l]放置于时延-多普勒域平面网格,进行OTFS调制;
其中,时延-多普勒域平面网格定义为Г:
Figure BDA0003312208580000091
其中k为多普勒域的索引,l为时延域的索引,1/(NT)表示多普勒轴的分辨率,1/(MΔf)表示时延轴的分辨率,(k)N表示对k取平移运算:
Figure BDA0003312208580000092
目标运动方向将导致多普勒频移取值可正可负,多普勒域索引k经过平移得到的(k)N的数值为正表示目标朝着靠近发射机的方向运动,数值为负则表示目标朝着远离发射机的方向运动,本发明方法所能探测的目标的时延范围是(0,1/Δf],多普勒频移范围是(-1/2T,1/2T];
步骤1.3.1.使用辛有限傅里叶逆变换(InverseSymplecticFiniteFourierTransform,ISFFT)将时延-多普勒域中的信息符号x[k,l]映射为时间-频率域上的符号X[n,m]:
Figure BDA0003312208580000093
其中,N代表一帧中时域上OTFS符号的数量,M代表频率域上子载波的数量;
所述时间-频率域平面网格定义为Λ:
Λ={(nT,mΔf),n=0,1,...,N-1,m=0,1,...,M-1},
其中n为时间的索引,m为频率的索引,T是一帧中一个OTFS符号的持续时间,Δf是相邻子载波间的频率间隔,且Δf=1/T;
步骤1.3.2.对时间-频率域平面网格上的符号X[n,m]应用海森堡变换(Heisenberg Transform),使用发射脉冲成形函数获得时间域上准备发射的连续时间的通信信号或者雷达信号s(t):
Figure BDA0003312208580000101
其中,gtx(t)为发射矩阵脉冲成形滤波器的冲激响应,t代表当前时刻;
步骤2.根据通信传输的覆盖能力和雷达探测的作用距离需求,发射机对雷达信号和通信信号进行功率控制,为雷达信号和通信信号分配合适的发射功率;
步骤2.1.由雷达方程确定雷达探测目标的作用距离与各雷达参数间的关系:
Figure BDA0003312208580000102
其中,Pr,rad表示发射机接收到的目标回波功率,Pt,rad表示雷达信号的发射功率,Gt表示发射天线的增益,λ表示工作的波长,探测信号所用波长λ乘以工作的载波频率fc等于光速,σ表示目标的散射截面积,Rrad为雷达探测目标的作用距离;
当雷达发射机接收到的回波功率等于发射机的最小可检测雷达信号功率时,得到雷达检测目标的最大作用距离:
Figure BDA0003312208580000103
其中,Rmax,rad表示雷达探测目标的最大作用距离,Smin,rad表示发射机的最小可检测雷达信号功率;并且,当接收到的目标回波功率Pr,rad大于最小可检测信号功率Smin,rad时,雷达探测才能可靠地检测目标,当Pr,rad等于Smin,rad时,雷达探测达到雷达检测目标的最大作用距离Rmax,rad,当实际目标距离大于Rmax,rad后,Pr,rad将小于Smin,rad,导致雷达探测无法有效可靠地检测目标;
步骤2.2.由弗里斯传输公式确定通信传输的覆盖范围与各通信参数间的关系:
Figure BDA0003312208580000104
其中,Pr,com表示目标接收机接收到的通信信号的功率,Pt,com表示通信信号的发射功率,Gr表示目标接收机的接收天线的增益,Rcom为发射机与目标通信的通信覆盖半径;
当目标接收机接收到的通信信号的功率等于目标接收机的通信信号接收功率阈值时,获得发射机与目标接收机能正常通信的最大通信覆盖半径:
Figure BDA0003312208580000111
其中,Rmax,com表示发射机与目标接收机能进行通信的最大通信覆盖半径,Smin,com表示目标接收机的通信信号接收功率阈值,并且,当接收到的通信信号功率Pr,com大于目标接收机的通信信号接收功率阈值Smin,com时,发射机才能与目标正常通信,当Pr,com与Smin,com相等时,对应能够通信的极限距离即最大通信覆盖半径,当发射机与目标的距离大于最大通信覆盖半径后,Pr,com将小于接收阈值Smin,com,导致无法正常通信传输;
步骤2.3.发射机进行发射信号的功率控制,为通信信号分配满足通信覆盖能力的发射功率,为雷达信号分配满足雷达探测距离需求的发射功率;
当雷达探测距离和通信传输覆盖的需求为雷达探测的最大作用距离大于通信传输的最大通信覆盖半径,实现更远距离的目标探测,但在最大通信覆盖半径之外只能探测目标而无法与目标通信时,则发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com满足如下关系:
Figure BDA0003312208580000112
当雷达探测距离和通信传输覆盖的需求为雷达最大作用距离等于最大通信覆盖半径,在最远距离处发射机既能和目标通信也能进行探测时,则发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com需满足如下关系:
Figure BDA0003312208580000113
当雷达探测距离和通信传输覆盖的需求为雷达最大作用距离小于最大通信覆盖半径,实现在更远距离与目标通信,但在雷达探测目标的最大作用距离之外只能与目标通信而无法进行有效探测时,则发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com需满足如下关系:
Figure BDA0003312208580000114
步骤3.通过选通转换开关控制通信与雷达的分时工作,为雷达探测功能和通信传输功能固定工作时隙,并插入防止雷达信号与通信信号互相干扰的保护带,实现雷达探测与通信传输以时分复用的方式共享天线、发射机和接收机平台;
对通信信源数据与雷达数据生成的通信信号和雷达信号分配固定的工作时隙,设定雷达探测时隙的长度为T1,即发射机在T1时间段内发射OTFS雷达信号执行探测功能,通信时隙的长度为T2,即发射机在T2时间段内发射OTFS通信信号执行通信功能,在T1和T2之间插入一个保护带Tg,防止在T1里发射的雷达信号经目标反射后,反射的雷达回波在T2里与通信信号产生冲突干扰,Tg的时间长度需满足:
Figure BDA0003312208580000121
其中,τmax表示雷达探测时发射机与待测目标的最大双程往返时延,取决于雷达探测目标的最大作用距离Rmax,rad,c表示光速;
通过选通转换开关实现通信功能与雷达功能的分时工作,控制发射机按如下时序循环工作,T1时间段发射雷达信号执行探测功能,Tg时间段作为保护带不发送任何信号,T2时间段发射通信信号执行通信功能;
步骤4.发射机处于通信传输的工作时段,则执行本步骤,发送OTFS通信信号与目标接收机进行信息传输,否则执行步骤5;
步骤4.1.发送OTFS连续时间通信信号,并经过信道传输后到达接收端,生成接收信号y(t):
y(t)=∫∫h(τ,ν)s(t-τ)ej2πν(t-τ)dνdτ,
其中,s(t)表示发送的通信信号,τ表示时延,ν表示多普勒频移,h(τ,ν)为时延-多普勒信道;
步骤4.2.接收端对接收信号y(t)进行OTFS解调,包括以下步骤:
步骤4.2.1.应用魏格纳变换(Wigner Transform)将接收信号y(t)转换为时间-频率域信号Y[n,m];
魏格纳变换的公式为:
Figure BDA0003312208580000122
其中,(·)*表示共轭运算,grx(t)表示接收脉冲成形函数,使用矩形脉冲,该式为计算grx(t)与y(t)的互模糊函数;
以间隔t=nT和f=mΔf进行采样,得到时频域信号Y[n,m]:
Figure BDA0003312208580000136
步骤4.2.2.对时间-频率域信号Y[n,m]使用辛有限傅里叶变换(SymplecticFinite Fourier Transform,SFFT),将时间-频率域信号Y[n,m]映射为时延-多普勒域上符号y[k,l]:
Figure BDA0003312208580000131
步骤4.3.对时延-多普勒域符号y[k,l]进行检测,得到估计信号
Figure BDA0003312208580000132
并通过星座解映射获取对应星座点上的估计比特,进而恢复出原始发送的串行比特通信信号流;
步骤5.发射机处于雷达探测功能的工作时段,发送雷达信号进行目标探测,探测到目标后反射雷达回波至发射机,由发射机进行雷达信号处理;
步骤5.1.考虑到雷达探测可工作在Sub-6GHz、毫米波或太赫兹等高频段,高频段时发射机与探测目标之间通常视作视距LoS(line of sight)传输,则当雷达存在P个探测目标时,将时延-多普勒域雷达信道建模为:
Figure BDA0003312208580000133
其中,hp表示第p(p=1,...,P)个目标的复信道增益,τp表示发射机与第p个目标的双程往返时延,νp表示第p个目标的双程多普勒频移;
τp满足下式:
Figure BDA0003312208580000134
上式中,Rp(m)表示第p个目标与发射机的相对距离;
νp满足下式:
Figure BDA0003312208580000135
上式中,Vp(m/s)表示第p个目标与发射机的相对运动速度,速度Vp可为正数或负数以表示运动方向,c是光速,fc为工作的载波频率;
根据目标的时延τp和多普勒频移νp分别近似时延分辨率1/(MΔf)和多普勒分辨率1/(NT)的整数倍,即τp=lp/(MΔf),νp=kp/(NT),其中lp和kp为整数,则时延-多普勒雷达信道的离散形式表达为:
Figure BDA0003312208580000141
其中,h[k,l]表示多普勒抽头k和时延抽头l处的复信道增益,δ(·)为狄利克雷函数;
对应于多普勒频移(k)N/(NT)和时延l/(MΔf),当k=kp,l=lp时h[k,l]不为0,表示在多普勒抽头kp和时延抽头lp位置处探测到目标,当0≤k≠kp≤N-1且0≤l≠lp≤M-1时h[k,l]为0,表示在多普勒抽头k和时延抽头l位置处没有目标;
步骤5.2.发射的雷达信号经过时延-多普勒域雷达信道h(τ,ν),在探测到目标后反射雷达回波至发射机,生成接收信号y(t),发射机对y(t)进行OTFS解调,应用魏格纳变换与SFFT得到时延-多普勒域中的接收数据符号y[k,l],y[k,l]与发送端信息符号x[k,l]间的输入-输出关系为:
Figure BDA0003312208580000142
其中,<·>N表示模N运算,<·>M表示模M运算,w[k,l]为噪声矩阵,并且满足w[k,l]~CN(0,σ2),表示w[k,l]为服从均值为零、方差为σ2的复高斯分布的加性高斯白噪声,α[k,l]表示相移,其表达式如下:
Figure BDA0003312208580000143
/>
步骤5.3.对得到时延-多普勒域上的接收数据符号y[k,l]通过雷达匹配滤波接收处理,进而获取探测目标对应于相对距离的时延特性和对应于相对运动速度的多普勒特性;
步骤5.3.1.分别将N×M维接收数据信号符号矩阵y[k,l]、发送信号符号矩阵x[k,l]、噪声矩阵w[k,l]和雷达信道增益矩阵h[k,l](0≤k≤N-1,0≤l≤M-1)按列优先的顺序进行排列,得到MN×1维接收符号列矢量
Figure BDA0003312208580000144
发送符号列矢量/>
Figure BDA0003312208580000145
噪声列矢量/>
Figure BDA0003312208580000146
和雷达信道列矢量/>
Figure BDA0003312208580000147
其中,y[k,l]、x[k,l]、w[k,l]和h[k,l]分别为
Figure BDA0003312208580000148
的第k+Nl个元素(0≤k+Nl≤MN-1),将各矢量之间关系用矩阵与矢量形式表达为:
Figure BDA0003312208580000151
其中,
Figure BDA0003312208580000152
的第一列为/>
Figure BDA0003312208580000153
其余列为/>
Figure BDA0003312208580000154
的循环移位,MN×MN维矩阵/>
Figure BDA0003312208580000155
的第i行第j列元素表示为:
Figure BDA0003312208580000156
其中,i满足0≤i=k′+Nl′≤MN-1,j满足0≤j=k″+Nl″≤MN-1;
步骤5.3.2.对得到的接收数据符号y[k,l]进行匹配滤波处理,处理过程为下式:
Figure BDA0003312208580000157
其中,(·)H表示共轭转置,
Figure BDA0003312208580000158
为匹配滤波处理估计的雷达信道矢量,第k+Nl个元素为多普勒抽头k和时延抽头l处的复信道增益/>
Figure BDA0003312208580000159
非零元素所在位置代表探测到的目标;
若当
Figure BDA00033122085800001510
时/>
Figure BDA00033122085800001511
不为0,表示在多普勒抽头/>
Figure BDA00033122085800001512
和时延抽头/>
Figure BDA00033122085800001513
位置处探测到第p个目标,目标相应的多普勒频移为/>
Figure BDA00033122085800001514
对应的目标与发射机的相对运动速度/>
Figure BDA00033122085800001515
为:
Figure BDA00033122085800001516
目标相应的时延为
Figure BDA00033122085800001517
对应的目标与发射机的相对距离/>
Figure BDA00033122085800001518
为:
Figure BDA00033122085800001519
本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书界定。

Claims (10)

1.基于时分体制的OTFS通信雷达一体化波形设计方法,其特征在于,所述通信雷达一体化波形设计方法,具体包括以下步骤:
步骤1.对通信数据或雷达数据进行数字调制和OTFS调制,得到通信发射信号或雷达发射信号;
步骤1.1.通信信源根据待通信传输的内容产生串行比特通信信号流,雷达数据产生器按照伪随机PN序列生成雷达比特数据流;
步骤1.2.将串行比特通信信号流与雷达比特数据流通过星座映射进行数字调制,得到相对应的信息符号x[k,l],数字调制的方式包括正交相移键控和正交振幅键控;
步骤1.3.将信息符号x[k,l]放置于时延-多普勒域平面网格,进行OTFS调制;
步骤2.根据通信传输的覆盖能力和雷达探测的作用距离需求,发射机对雷达信号和通信信号进行功率控制,为雷达信号和通信信号分配发射功率;
步骤2.1.由雷达方程确定雷达探测目标的作用距离与各雷达参数间的关系;
步骤2.2.由弗里斯传输公式确定通信传输的覆盖范围与各通信参数间的关系;
步骤2.3.发射机进行发射信号的功率控制,为通信信号分配满足通信覆盖能力的发射功率,为雷达信号分配满足雷达探测距离需求的发射功率;
步骤3.通过选通转换开关控制通信与雷达的分时工作,为雷达探测功能和通信传输功能固定工作时隙,并插入防止雷达信号与通信信号互相干扰的保护带,实现雷达探测与通信传输以时分复用的方式共享天线、发射机和接收机平台;
步骤4.发射机处于通信传输的工作时段,则执行本步骤,发送OTFS通信信号与目标接收机进行信息传输,否则执行步骤5:
步骤4.1.发送OTFS连续时间通信信号,并经过信道传输后到达接收端,生成接收信号y(t);
步骤4.2.接收端对接收信号y(t)进行OTFS解调;
步骤4.3.对时延-多普勒域符号y[k,l]进行检测,得到估计信号
Figure FDA0003312208570000011
并通过星座解映射获取对应星座点上的估计比特,进而恢复出原始发送的串行比特通信信号流;
步骤5.发射机处于雷达探测功能的工作时段,发送雷达信号进行目标探测,探测到目标后反射雷达回波至发射机,由发射机进行雷达信号处理:
步骤5.1.建立时延-多普勒域雷达信道h(τ,ν);
步骤5.2.发射的雷达信号经过时延-多普勒域雷达信道h(τ,ν),在探测到目标后反射雷达回波至发射机,生成接收信号y(t),发射机对y(t)进行OTFS解调,应用魏格纳变换与辛有限傅里叶逆变换,得到时延-多普勒域上的接收数据符号y[k,l];
步骤5.3.对得到时延-多普勒域上的接收数据符号y[k,l]通过雷达匹配滤波接收处理,进而获取探测目标对应于相对距离的时延特性和对应于相对运动速度的多普勒特性。
2.根据权利要求1所述的通信雷达一体化波形设计方法,其特征在于,所述步骤1.3.中时延-多普勒域平面网格定义为Γ,并满足下式(1):
Figure FDA0003312208570000021
上式(1)中,k为多普勒域的索引,l为时延域的索引,1/(NT)表示多普勒轴的分辨率,1/(MΔf)表示时延轴的分辨率,(k)N表示对k取平移运算,并满足下式(2):
Figure FDA0003312208570000022
上式(2)中,目标运动方向将导致多普勒频移取值为正或者负,多普勒域索引k经过平移得到的(k)N的数值为正表示目标朝着靠近发射机的方向运动,数值为负则表示目标朝着远离发射机的方向运动。
3.根据权利要求2所述的通信雷达一体化波形设计方法,其特征在于,所述步骤1.3.包括以下步骤:
步骤1.3.1.使用辛有限傅里叶逆变换将时延-多普勒域中的信息符号x[k,l]映射为时间-频率域上的符号X[n,m],映射过程满足下式(3):
Figure FDA0003312208570000023
上式(3)中,X[n,m]表示时间-频率域网格上的信号,x[k,l]表示时延-多普勒域网格上的信号,N代表一帧中时域上OTFS符号的数量,M代表频率域上子载波的数量;
步骤1.3.2.对时间-频率域平面网格上的符号X[m,n]应用海森堡变换,使用发射脉冲成形函数获得时间域上准备发射的连续时间的通信信号或者雷达信号s(t),连续时间的通信信号或者雷达信号s(t)变换公式为下式(4):
Figure FDA0003312208570000024
上式(4)中,gtx(t)为发射矩阵脉冲成形滤波器的冲激响应,t代表当前时刻。
4.根据权利要求3所述的通信雷达一体化波形设计方法,其特征在于,所述步骤1.3.2.中时间-频率域平面网格定义为Λ,并满足下式(5):
Λ={(nT,mΔf),n=0,1,...,N-1,m=0,1,...,M-1} …(5),
上式(5)中,n为时间的索引,m为频率的索引,T是一帧中一个OTFS符号的持续时间,Δf是相邻子载波间的频率间隔,且Δf=1/T。
5.根据权利要求1所述的通信雷达一体化波形设计方法,其特征在于,所述步骤2.1.中,由雷达方程得到下式(6):
Figure FDA0003312208570000031
上式(6)中,Pr,rad表示发射机接收到的目标回波功率,Pt,rad表示雷达信号的发射功率,Gt表示发射天线的增益,λ表示工作的波长,且探测信号所用波长λ乘以工作的载波频率fc等于光速,σ表示目标的散射截面积,Rrad为雷达探测目标的作用距离;
由上式(6)得知,当雷达发射机接收到的回波功率等于发射机的最小可检测雷达信号功率时,即得到雷达检测目标的最大作用距离Rmax,rad如下式(7):
Figure FDA0003312208570000032
上式(7)中,Smin,rad表示发射机的最小可检测雷达信号功率。
6.根据权利要求1所述的通信雷达一体化波形设计方法,其特征在于,所述步骤2.2.中,由弗里斯传输公式得到下式(8):
Figure FDA0003312208570000033
上式(8)中,Pr,com表示目标接收机接收到的通信信号的功率,Pt,com表示通信信号的发射功率,Gr表示目标接收机的接收天线的增益,Ct表示发射天线的增益,Rcom为发射机与目标通信的通信覆盖半径;
由上式(8)得知,当目标接收机接收到的通信信号的功率等于目标接收机的通信信号接收功率阈值时,获得发射机与目标接收机正常通信的最大通信覆盖半径Rmaxcom,如下式(9):
Figure FDA0003312208570000034
上式(9)中,Smin,com表示目标接收机的通信信号接收功率阈值。
7.根据权利要求1所述的通信雷达一体化波形设计方法,其特征在于,所述步骤2.3.中,发射机进行发射信号的功率控制,具体包括以下三种情况:
当雷达探测距离和通信传输覆盖的需求为雷达探测的最大作用距离大于通信传输的最大通信覆盖半径时,则发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com满足下式(10):
Figure FDA0003312208570000041
当雷达探测距离和通信传输覆盖的需求为雷达最大作用距离等于最大通信覆盖半径时,发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com满足下式(11):
Figure FDA0003312208570000042
当雷达探测距离和通信传输覆盖的需求为雷达最大作用距离小于最大通信覆盖半径时,发射机为雷达信号和通信信号分配的发射功率Pt,rad和Pt,com满足下式(12):
Figure FDA0003312208570000043
8.根据权利要求1所述的通信雷达一体化波形设计方法,其特征在于,所述步骤3.中为通信信号和雷达信号分配固定的工作时隙,并插入防止雷达信号与通信信号互相干扰的保护带,具体方法如下:
设定发射机在T1时间段内发射OTFS雷达信号执行探测功能,设定发射机在T2时间段内发射OTFS通信信号执行通信功能,则T1和T2之间插入的保护带Tg满足下式(13):
Figure FDA0003312208570000044
上式(13)中,τmax表示雷达探测时发射机与待测目标的最大往返时延,Rmax,rad表示雷达探测目标的最大作用距离,c表示光速。
9.根据权利要求1所述的通信雷达一体化波形设计方法,其特征在于,所述步骤4.1.中,OTFS连续时间通信信号经过信道传输后到达接收端并生成接收信号y(t)的过程满足下式(14):
y(t)=∫∫h(τ,ν)s(t-τ)ej2πν(t-τ)dνdτ …(14),
上式(14)中,s(t)表示发送的通信信号,τ表示时延,ν表示多普勒频移,h(τ,ν)为时延-多普勒信道。
10.根据权利要求1所述的通信雷达一体化波形设计方法,其特征在于,所述步骤4.2.接收端对接收信号y(t)进行OTFS解调,包括以下步骤:
步骤4.2.1.应用魏格纳变换将接收信号y(t)转换为时间-频率域信号Y[n,m];
魏格纳变换的公式为下式(15):
Figure FDA0003312208570000051
上式(15)中,(·)*表示共轭运算,grx(t)表示接收脉冲成形函数,Y(t,f)为grx(t)与y(t)的互模糊函数;
以间隔t=nT和f=mΔf进行采样,得到时间-频率域信号Y[n,m]如下式(16):
Figure FDA0003312208570000052
步骤4.2.2.对时间-频率域信号Y[n,m]使用辛有限傅里叶变换,将时间-频率域信号Y[n,m]映射为时延-多普勒域上符号y[k,l],映射过程满足下式(17):
Figure FDA0003312208570000053
上式(17)中,N代表时域上OTFS符号的数量,M代表频率域上子载波的数量;
所述步骤5.1.中,时延-多普勒域雷达信道h(τ,ν)的建模表达式为下式(18):
Figure FDA0003312208570000054
上式(18)中,hp表示第p(p=1,...,P)个目标的复信道增益,τp表示发射机与第p个目标的双程往返时延,νp表示第p个目标的双程多普勒频移;
其中,τp满足下式(19):
Figure FDA0003312208570000055
上式(19)中,Rp(m)表示第p个目标与发射机的相对距离;
νp满足下式(20):
Figure FDA0003312208570000061
上式(20)中,Vp(m/s)表示第p个目标与发射机的相对运动速度,fc为工作的载波频率;
由式(18)~(20)得到离散形式的时延-多普勒域雷达信道h(τ,ν),如下式(21):
Figure FDA0003312208570000062
上式(21)中,h[k,l]表示多普勒抽头k和时延抽头l处的复信道增益,δ(·)为狄利克雷函数;
所述步骤5.2.中,时延-多普勒域上的接收数据符号y[k,l]与发送端信息符号x[k,l]间的输入输出关系满足下式(22):
Figure FDA0003312208570000063
上式(22)中,<·>N表示模N运算,<·>M表示模M运算,w[k,l]为噪声矩阵,并且满足w[k,l]~CN(0,σ2),表示w[k,l]为服从均值为零、方差为σ2的复高斯分布的加性高斯白噪声,α[k,l]表示相移;
α[k,l]的表达式如下式(23):
Figure FDA0003312208570000064
所述步骤5.3.中,对得到的接收数据符号y[k,l]进行雷达处理,具体包括如下步骤:
步骤5.3.1.分别将N×M维接收数据信号符号矩阵y[k,l]、发送信号符号矩阵x[k,l]、噪声矩阵w[k,l]和雷达信道增益矩阵h[k,l](0≤k≤N-1,0≤l≤M-1)按列优先的顺序进行排列,得到MN×1维接收符号列矢量
Figure FDA0003312208570000065
发送符号列矢量
Figure FDA0003312208570000066
噪声列矢量
Figure FDA0003312208570000067
和雷达信道列矢量
Figure FDA0003312208570000068
其中,y[k,l]、x[k,l]、w[k,l]和h[k,l]分别为y,x,w,
Figure FDA0003312208570000069
的第k+Nl个元素(0≤k+Nl≤MN-1),将各矢量之间关系用矩阵与矢量形式表达为下式(24):
Figure FDA0003312208570000071
上式(24)中,
Figure FDA0003312208570000072
的第一列为
Figure FDA0003312208570000073
其余列为
Figure FDA0003312208570000074
的循环移位;
将公式(24)中MN×MN维矩阵
Figure FDA0003312208570000075
的第i行第j列元素表示为下式(25):
Figure FDA0003312208570000076
上式(25)中,i满足0≤i=k′+Nl′≤MN-1,j满足0≤j=k″+Nl″≤MN-1;
步骤5.3.2.对得到的接收数据符号y[k,l]进行匹配滤波处理,处理过程为下式(26):
Figure FDA0003312208570000077
上式(26)中,(·)H表示共轭转置,
Figure FDA0003312208570000078
为匹配滤波处理估计的雷达信道矢量,第k+Nl个元素为多普勒抽头k和时延抽头l处的复信道增益
Figure FDA0003312208570000079
非零元素所在位置代表探测到的目标;
若当
Figure FDA00033122085700000710
时,
Figure FDA00033122085700000711
不为0,表示在多普勒抽头
Figure FDA00033122085700000712
和时延抽头
Figure FDA00033122085700000713
位置处探测到第p个目标,目标相应的多普勒频移为
Figure FDA00033122085700000714
对应的目标与发射机的相对运动速度
Figure FDA00033122085700000715
为下式(27):
Figure FDA00033122085700000716
目标相应的时延为
Figure FDA00033122085700000717
对应的目标与发射机的相对距离
Figure FDA00033122085700000718
为下式(28):
Figure FDA00033122085700000719
CN202111219969.4A 2021-10-20 2021-10-20 基于时分体制的otfs通信雷达一体化波形设计方法 Active CN114124238B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111219969.4A CN114124238B (zh) 2021-10-20 2021-10-20 基于时分体制的otfs通信雷达一体化波形设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111219969.4A CN114124238B (zh) 2021-10-20 2021-10-20 基于时分体制的otfs通信雷达一体化波形设计方法

Publications (2)

Publication Number Publication Date
CN114124238A CN114124238A (zh) 2022-03-01
CN114124238B true CN114124238B (zh) 2023-03-24

Family

ID=80376151

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111219969.4A Active CN114124238B (zh) 2021-10-20 2021-10-20 基于时分体制的otfs通信雷达一体化波形设计方法

Country Status (1)

Country Link
CN (1) CN114124238B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114594429B (zh) * 2022-03-24 2023-11-14 南京慧尔视智能科技有限公司 一种雷达抗干扰处理方法及装置
CN114944896B (zh) * 2022-04-08 2023-10-24 江苏科技大学 一种基于y型网络的索引调制otfs通信系统及方法
CN117336863A (zh) * 2022-06-22 2024-01-02 华为技术有限公司 一种通信方法和通信装置
WO2024007211A1 (en) * 2022-07-06 2024-01-11 Huawei Technologies Co., Ltd. Methods and apparatus for power domain multiplexing of communication and sensing signals
CN114935754A (zh) * 2022-07-20 2022-08-23 深圳芯盛思技术有限公司 采用数字调制技术的调频连续波雷达调制系统及使用方法
CN115396263B (zh) * 2022-07-29 2023-06-30 北京邮电大学 一种otfs通信感知一体化信号目标参数估计方法
CN115575916A (zh) * 2022-10-21 2023-01-06 平湖空间感知实验室科技有限公司 近地小行星地基雷达探测目标筛选方法、系统和电子设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107925446A (zh) * 2015-05-11 2018-04-17 凝聚技术股份有限公司 正交时间频率空间调制系统
CA3101010A1 (en) * 2018-07-27 2020-01-30 Raytheon Company Multidimensional shared spectrum access between radio and radar signals
CN111585644A (zh) * 2020-05-27 2020-08-25 北京邮电大学 雷达通信一体化系统、信号处理方法及设备、存储介质
CN112882005A (zh) * 2021-01-14 2021-06-01 东南大学 一种基于贝叶斯学习的otfs雷达目标参数估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9929783B2 (en) * 2012-06-25 2018-03-27 Cohere Technologies, Inc. Orthogonal time frequency space modulation system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107925446A (zh) * 2015-05-11 2018-04-17 凝聚技术股份有限公司 正交时间频率空间调制系统
CA3101010A1 (en) * 2018-07-27 2020-01-30 Raytheon Company Multidimensional shared spectrum access between radio and radar signals
CN111585644A (zh) * 2020-05-27 2020-08-25 北京邮电大学 雷达通信一体化系统、信号处理方法及设备、存储介质
CN112882005A (zh) * 2021-01-14 2021-06-01 东南大学 一种基于贝叶斯学习的otfs雷达目标参数估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
雷达通信一体化共享信号技术探讨;王洋;《信息通信》;20200731;203-204页 *

Also Published As

Publication number Publication date
CN114124238A (zh) 2022-03-01

Similar Documents

Publication Publication Date Title
CN114124238B (zh) 基于时分体制的otfs通信雷达一体化波形设计方法
Xiao et al. Waveform design and performance analysis for full-duplex integrated sensing and communication
Wei et al. Integrated sensing and communication signals toward 5G-A and 6G: A survey
Kumari et al. IEEE 802.11 ad-based radar: An approach to joint vehicular communication-radar system
Sturm et al. An OFDM system concept for joint radar and communications operations
Ozkaptan et al. OFDM pilot-based radar for joint vehicular communication and radar systems
WO2019166429A1 (en) Method for joint radar-communication
KR101814486B1 (ko) 레이더 및 무선 통신 기능을 가지는 차량용 레이더 장치 및 그것을 이용한 자원 조절 방법
CN111585644B (zh) 雷达通信一体化系统、信号处理方法及设备、存储介质
Dokhanchi et al. Joint automotive radar-communications waveform design
CN107728110B (zh) 一种基于lfm信号相位调制的共享信号设计方法
CN112363132A (zh) 一种基于fbmc的雷达通信一体化波形生成方法
US20040005016A1 (en) High bit rate ultra-wideband OFDM
CN113612716B (zh) 基于恒定包络正交频分复用的多用户通信雷达一体化系统
CN116056104A (zh) 感知通信一体化ofdm共享波形子载波的比特和功率分配方法
Wang et al. OFDM chirp waveform diversity for co-designed radar-communication system
Huang et al. FDA-OFDM for integrated navigation, sensing, and communication systems
Dokhanchi et al. Multicarrier phase modulated continuous waveform for automotive joint radar-communication system
Wu et al. A non-uniform multi-wideband OFDM system for terahertz joint communication and sensing
Tang et al. Self-interference-resistant IEEE 802.11 ad-based joint communication and automotive long range radar
Zacharia Performance analysis of OTFS signal with different pulse shapes for JCR systems
CN114866124A (zh) 基于时延对齐调制的通感一体化信号设计与波束赋形方法
Li et al. Adaptive filter bank multi-carrier waveform design for joint communication-radar system
Su et al. Improving Resource Efficiency of PMCW-Based JCRS Systems: Simultaneous Transmission of Pilot and Data via Orthogonal Codes
Al-Salehi et al. Dual-function radar–communications: information transmission during FDA radar listening mode

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