CN116660856A - 一种基于5g时隙同步的外辐射源雷达信号处理方法 - Google Patents
一种基于5g时隙同步的外辐射源雷达信号处理方法 Download PDFInfo
- Publication number
- CN116660856A CN116660856A CN202310962814.2A CN202310962814A CN116660856A CN 116660856 A CN116660856 A CN 116660856A CN 202310962814 A CN202310962814 A CN 202310962814A CN 116660856 A CN116660856 A CN 116660856A
- Authority
- CN
- China
- Prior art keywords
- signal
- representing
- radiation source
- signals
- ssb
- 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
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 37
- 238000003672 processing method Methods 0.000 title claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 22
- 230000005540 biological transmission Effects 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 27
- 125000004122 cyclic group Chemical group 0.000 claims description 20
- 230000001360 synchronised effect Effects 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 11
- 238000000605 extraction Methods 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 20
- 238000012544 monitoring process Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 5
- 230000006978 adaptation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 102220541155 Tumor necrosis factor receptor superfamily member 5_S35G_mutation Human genes 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details 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/411—Identification of targets based on measurements of radar reflectivity
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明公开了一种基于5G时隙同步的外辐射源雷达信号处理方法，针对5G信号采用TDD传输的问题，通过搜索并解码SSB信号实现5G信号帧的时频同步，在移除上行时隙后，进行距离‑多普勒处理，可以避免上行信号对目标探测的影响，有效提高目标分辨率和回波信噪比，分别得到目标的时延和多普勒频移信息，实现对目标的有效探测。本发明提供的一种基于5G时隙同步的外辐射源雷达信号处理方法，利用5G信号作为外辐射源雷达信号进行目标探测，低空覆盖能力好，无电磁污染，可全天候连续监控；相对于常用的FM、4G等辐射源信号，距离分辨率和速度分辨率更好。
Description
技术领域
本发明涉及一种基于5G时隙同步的外辐射源雷达信号处理方法，属于外辐射源雷达目标探测技术领域。
背景技术
目前，低空目标防控手段主要包括无线电侦测、光电探测、声波探测和雷达探测等。其中，外辐射源雷达是一种利用非配合辐射源信号进行目标探测的设备，又称被动雷达，属于雷达设备中的一种。但与常规雷达不同，其本身并不发射信号，而是通过接收目标反射的非合作照射源（广播、电视、通信基站等）回波信号进行探测，能以无线电静默的方式对多种类型和多批次的低空目标进行探测，具有成本低无辐射、组网探测能力强等优势，特别适合于对电磁辐射有限制同时又对探测性能要求较高的城市环境。
随着第五代移动通信（5G）的迅速发展，相对于其他辐射源而言，将5G信号作为非合作照射源进行目标检测具有如下优点：
1）、支持更灵活的子载波间隔，将时隙作为调度单位，能够根据不同需求配置帧结构参数，适配多样化业务需求和资源部署。
2）、通信基站布设密集，采用MIMO天线设计和波束赋形等技术，通信覆盖率更高。
3）、频段更宽，通信带宽更大，最大频谱宽度可达100MHz，使得距离分辨率和速度分辨率更高。
但与其他辐射源不同，5G基站采用时分双工（Time Division Duplexing，TDD）传输模式，共有三种“类型”的时隙：1）、下行链路时隙（DL），即5G基站的发射信号；2）、上行链路时隙（UL），即通信终端的发射信号；以及3）、用于配置DL/UL模式的特殊（SP）时隙。由于上行时隙和下行时隙的发射终端不同，在外辐射源雷达系统距离-多普勒处理的过程中，上行时隙会对目标探测结果产生影响。
因此，基于5G的外辐射源雷达系统信号处理方法不同，需要加以改进，以实现对低空目标的稳定探测。
发明内容
目的：为了克服现有技术中存在的不足，本发明提供一种基于5G时隙同步的外辐射源雷达信号处理方法。
技术方案：为解决上述技术问题，本发明采用的技术方案为：
一种基于5G时隙同步的外辐射源雷达信号处理方法，包括如下步骤：
步骤S1、获取中心频率为、采样率为/>的5G基站直达波信号/>，目标回波信号/>。
步骤S2、根据5G协议生成3种本地主同步信号，并将3种本地主同步信号调制为对应的时域信号/>，考虑整数倍的最小子载波间隔/>的频偏，将时域信号/>分别与5G基站直达波信号/>进行互相关运算得到PSS互相关结果，根据PSS互相关结果/>的最大值确定SSB同步信号块时域信号的时间偏移和小区ID号/>。
步骤S3、根据时间偏移从5G基站直达波信号/>中提取4个CP-OFDM符号长度的SSB同步信号块时域信号/>，并获取SSB同步信号块时域信号/>的循环前缀部分/>和重复的数据部分/>，将循环前缀部分/>与重复的数据部分进行互相关运算得到频偏互相关结果/>，根据频偏互相关结果/>估计5G信号的频率偏移/>。
步骤S4、解码SSB同步信号块时域信号，提取SSS序列/>，根据5G协议生成336种本地辅同步信号/>，分别与SSS序列/>进行互相关运算得到SSS互相关结果/>，根据SSS互相关结果/>的最大值确定小区ID号/>。
步骤S5、根据小区ID号和小区ID号/>求取联合小区ID号/>，解码SSB同步信号块时域信号/>，根据联合小区ID号/>提取DMRS序列/>，根据5G协议生成8种本地解调参考信号/>，分别与/>进行互相关运算得到DMRS互相关结果，根据DMRS互相关结果/>的最大值确定SSB信号索引/>。
步骤S6、根据SSB传输模式确定索引为的SSB信号开始符号位置，并根据SSB信号开始符号位置与SSB信号的时间偏移/>联合确定信号帧的时间偏移/>，并通过/>和频率偏移/>调整5G基站直达波信号/>和目标回波信号/>，分别得到同步后的5G基站直达波信号/>和目标回波信号/>。
步骤S7、根据不同5G基站TDD传输模式将同步后的5G基站直达波信号中的灵活时隙和上行时隙部分置0，得到下行信号/>。
步骤S8、对下行信号和同步后的目标回波信号/>计算互模糊函数，根据互模糊函数/>最大值得到目标的时延/>和多普勒频移/>信息。
作为优选方案，所述时间偏移和小区ID号/>的计算公式如下：
；
式中，
；
其中，表示3种本地主同步信号的序号，/>表示取共轭，/>表示时延，/>表示整数倍的最小子载波间隔/>，/>表示取最大值，/>表示/>取最大值时对应的/>，表示/>取最大值时对应的/>，/>表示e的指数函数，/>表示/>取最大值时对应的整数倍，/>表示虚数单位，/>表示圆周率，/>表示时间。
作为优选方案，所述5G信号的频率偏移计算公式如下：
；
式中，；
其中，表示取共轭，/>表示求和，/>表示时延，/>为角度运算符，/>为信号子载波间隔。
作为优选方案，所述SSB同步信号块时域信号计算公式如下：
；
所述循环前缀部分和重复的数据部分/>计算公式如下：
；
；
其中，为采样率，/>为OFDM符号循环前缀CP短长度，/>为OFDM符号长度，为CP-OFDM符号短长度，/>为CP-OFDM符号序数，对应/>。
作为优选方案，所述小区ID号计算公式如下：
；
式中，；
其中，表示离散序列的自变量，/>表示离散序列的滑动点数，/>表示SSS序列的长度，/>表示取共轭，/>表示336种本地辅同步信号的序号，/>表示取最大值，/>表示取最大值时对应的/>。
作为优选方案，所述SSB信号索引计算公式如下：
；
式中，；
其中，表示离散序列的自变量，/>表示离散序列的滑动点数，/>表示DMRS序列的长度，/>表示取共轭，/>表8种本地解调参考信号的序号，/>表示取最大值，/>表示取最大值时对应的/>。
作为优选方案，所述同步后的5G基站直达波信号和目标回波信号/>，计算公式如下：
；
；
式中，当SSB传输模式为Case C时，计算公式如下：
；
其中，为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为采样率。
作为优选方案，所述下行信号计算公式如下：
当TDD传输模式为2.5ms单周期时，下行信号计算公式如下：
；
式中：
；
；
其中，为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为采样率，为周期数。
当TDD传输模式为2.5ms双周期时，下行信号计算公式如下：
；
式中：
；
当TDD传输模式为5ms单周期时，下行信号计算公式如下：
；
式中：
；
。
作为优选方案，，/>；
其中，为正交频分复用OFDM符号长度，/>为正交频分复用OFDM循环前缀短长度，/>为正交频分复用OFDM循环前缀长长度。
作为优选方案，所述目标的时延和多普勒频移/>计算公式如下：
；
式中：；
其中，表示取共轭，/>表示时延，/>表示多普勒频率，/>表示取最大值，/>表示/>取最大值时对应的/>，/>表示/>取最大值时对应的/>。
有益效果：本发明提供的一种基于5G时隙同步的外辐射源雷达信号处理方法，利用5G信号作为外辐射源雷达信号进行目标探测，低空覆盖能力好，无电磁污染，可全天候连续监控；相对于常用的FM、4G等辐射源信号，距离分辨率和速度分辨率更好。针对5G信号采用TDD传输的问题，通过搜索并解码SSB信号实现5G信号帧的时频同步，在移除上行时隙后，进行距离-多普勒处理，可以避免上行信号对目标探测的影响，有效提高目标分辨率和回波信噪比，分别得到目标的时延和多普勒频移信息，实现对目标的有效探测。
附图说明
图1是本发明使用的外辐射源雷达系统结构示意图。
图2是本发明使用的5G信号SSB传输模式（Case C）示意图，对应于中国移动5G信号（2.515~2.615GHz）。
图3是本发明使用的5G信号TDD传输模式（5ms单周期）示意图，对应于中国移动5G信号（2.515~2.615GHz）。
图4是本发明使用的一种基于5G时隙同步的外辐射源雷达信号处理方法的流程图。
图5是本发明实施例中参考通道信号时频图。
图6是本发明实施例中同步后的参考通道信号时频图。
图7是本发明实施例中上行时隙移除后的参考通道信号时频图。
图8是本发明实施例中距离-多普勒处理结果图。
具体实施方式
下面结合本发明实例中的附图，对本发明实例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域技术人员在没有做创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。
下面结合具体实施例对本发明作更进一步的说明。
本发明的技术方案针对的对象是5G下行信号外辐射源雷达，其一般结构如图1所示。其中，外辐射源雷达包括：指向5G基站方向的参考通道天线用于获取5G基站直达波信号，指向探测空域的监测通道天线用于获取目标回波信号/>，和完成后续信号处理的信号处理单元。
在本实施例中，外辐射源雷达接收到的参考信号和目标回波信号/>为中国移动5G信号（2.515GHz~2.615GHz），信号子载波间隔为/>，最小子载波间隔/>，信号帧长为10ms，信号子帧时长为1ms，信号时隙时长为0.5ms，每个时隙中包含14个CP-OFDM（Cyclic Prefix Orthogonal Frequency Division Multiplexing，循环前缀正交频分复用）符号。根据5G技术协议，中国移动5G信号的SSB（SynchronizationSignal/Physical Broadcast Channel block，同步信号/广播信道块）传输模式（Case C）示意图如图2所示，在20ms周期内传输8个SSB信号，对应索引/>为0~7，每个SSB开始符号位置为2，8，16，22，30，36，44，50（第一个符号位置记为0），并分别占据4个CP-OFDM符号长度；此外，中国移动5G信号TDD传输模式（5ms单周期）示意图如图3所示，在5ms周期内共有10个0.5ms时隙，其中前7个时隙为下行时隙，第8个时隙为灵活时隙，最后2个时隙为上行时隙，在灵活时隙中前6个CP-OFDM符号为下行符号，中间4个为灵活符号，最后4个为上行符号。
本发明一种基于5G时隙同步的外辐射源雷达信号处理方法实施过程和信号处理流程如图4所示，具体描述为以下过程：
步骤S1、外辐射源雷达接收机通过参考通道获取中心频率为、采样率为/>的5G基站直达波信号/>，并通过监测通道获取目标回波信号/>，表示为：
;
其中，表示5G基站的发射信号，/>表示外辐射源雷达接收机参考通道的噪声，/>表示目标回波的幅度，/>表示目标的时延，/>表示目标的多普勒频率，/>表示外辐射源雷达接收机监测通道的噪声，/>表示虚数单位，/>表示圆周率，/>表示时间。其中，/>为所要求得的时延，/>为所要求得的多普勒频率。
在本实施例中，中心频率为、采样率为/>。图5是本发明实施例中参考通道信号时频图。
步骤S2、参数获取：从参考通道的5G基站直达波信号中读取后续处理所需参数值，包括：OFDM（Orthogonal Frequency Division Multiplexing）正交频分复用符号长度/>、CP（Cyclic Prefix）循环前缀短长度/>和循环前缀长长度/>、CP-OFDM（Cyclic Prefix Orthogonal Frequency Division Multiplexing）循环前缀正交频分复用符号短长度/>，/>和循环前缀正交频分复用符号长长度，/>。
在本实施例中，OFDM符号长度、循环前缀CP短长度/>和循环前缀CP长长度/>、CP-OFDM符号短长度/>和CP-OFDM符号长长度/>。
步骤S3、PSS（Primary Synchronization Signal）主同步信号检测：根据5G协议生成3种本地主同步信号，并将/>调制为对应的时域信号/>，考虑整数倍的最小子载波间隔/>的频偏，将时域信号/>分别与S1中的参考通道的5G基站直达波信号/>进行滑动互相关运算得到PSS互相关结果/>，根据/>的最大值确定SSB（Synchronization Signal/Physical Broadcast Channel block）同步信号块时域信号的时间偏移/>和小区ID号/>，计算公式表示为：
；
其中，表示离散序列的自变量，/>表示3种本地主同步信号的序号，/>表示取共轭，/>表示时延，/>表示整数倍的最小子载波间隔/>，/>表示取最大值，/>表示取最大值时对应的/>，/>表示/>取最大值时对应的/>，/>表示e的指数函数，/>表示/>取最大值时对应的整数倍。
步骤S4、频率偏移估计：根据S3得到的时间偏移从S1的参考通道的5G基站直达波信号/>中提取4个CP-OFDM符号长度的SSB同步信号块时域信号/>，计算公式表示为：
；
的循环前缀部分/>和重复的数据部分/>表示为：
；
；
其中，为采样率，/>为OFDM符号循环前缀CP短长度，/>为OFDM符号长度，为CP-OFDM符号短长度，/>为CP-OFDM符号序数，对应/>。
将的循环前缀部分/>和重复的数据部分/>进行互相关运算得到频偏互相关结果/>，根据/>估计5G信号的频率偏移/>，计算公式表示为：
；
；
其中，表示取共轭，/>表示求和，/>表示时延，/>为角度运算符，/>为信号子载波间隔。
步骤S5、SSS（Secondary Synchronization Signal）辅同步信号检测：解码S4中的SSB同步信号块时域信号，提取SSS序列/>，根据5G协议生成336种本地辅同步信号/>，分别与SSS序列/>进行互相关运算得到SSS互相关结果/>，根据的最大值确定小区ID号/>，计算公式表示为：
；
其中，表示离散序列的自变量，/>表示离散序列的滑动点数，/>表示SSS序列的长度，/>表示取共轭，/>表示336种本地辅同步信号的序号，/>表示取最大值，/>表示取最大值时对应的/>。
步骤S6、DMRS（Demodulation Reference Signal）解调参考信号估计：求取联合小区ID号，/>，解码S4中的SSB时域信号/>，根据联合小区ID号/>提取DMRS序列/>，根据5G协议生成8种本地解调参考信号/>，分别与进行互相关运算得到DMRS互相关结果/>，根据/>的最大值确定SSB信号索引/>，计算公式表示为：
；
其中，表示离散序列的自变量，/>表示离散序列的滑动点数，/>表示DMRS序列的长度，/>表示取共轭，/>表8种本地解调参考信号的序号，/>表示取最大值，/>表示取最大值时对应的/>。
步骤S7、信号时频同步：根据SSB传输模式确定S6中索引为的SSB信号开始符号位置，与S3中SSB信号的时间偏移/>联合确定信号帧的时间偏移/>，并通过/>和S4得到的频率偏移/>调整S1中5G基站直达波信号/>和目标回波信号/>，分别得到同步后的5G基站直达波信号/>和目标回波信号/>，计算公式表示为：
；
；
式中，当SSB传输模式为Case C时，信号帧的时间偏移，表示为：
；
其中，为SSB信号的时间偏移，/>为SSB信号索引，/>为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为采样率。
在本实施例中，SSB传输模式为Case C，SSB信号索引。图6所示为基于本实施例的信号处理结果，具体为同步后的参考通道信号时频图。
步骤S8、上行时隙移除：根据不同运营商5G基站使用的2.5ms单周期、2.5ms双周期和5ms单周期三种TDD传输模式将S7得到的中的灵活时隙和上行时隙部分置0，得到下行信号/>。
当TDD传输模式为2.5ms单周期时，将S7得到的中灵活时隙和上行时隙置0，得到下行信号/>，计算公式表示为：
；
；
；
其中，为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为一个周期的信号长度，/>为第一个周期的下行时隙长度，/>为同步后的5G基站直达波信号，/>为采样率，/>为周期数。
当TDD传输模式为2.5ms双周期时，将S7得到的中灵活时隙和上行时隙置0，得到下行信号/>，计算公式表示为：/>
；
；
其中，为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为一个周期的信号长度，/>为第一个周期的下行时隙长度，/>为第二个周期的下行时隙长度，/>为同步后的5G基站直达波信号，/>为采样率，/>为周期数。
当TDD传输模式为5ms单周期时，将S7得到的中灵活时隙和上行时隙置0，得到下行信号/>，计算公式表示为：
；
；
；
其中，为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为一个周期的信号长度，/>为第一个周期的下行时隙长度，/>为同步后的5G基站直达波信号，/>为采样率，/>为周期数。
在本实施例中，中国移动5G基站使用的TDD传输模式为5ms单周期。图7所示为基于本实施例的信号处理结果，具体为上行时隙移除后的参考通道信号时频图。
步骤S9、距离-多普勒处理：对S8得到的下行信号和S7得到的同步后的目标回波信号/>计算互模糊函数/>，根据互模糊函数/>最大值得到目标的时延/>和多普勒频移/>信息，实现目标探测，计算公式表示为：
；/>
；
其中，表示取共轭，/>表示时延，/>表示多普勒频率，/>表示取最大值，/>表示/>取最大值时对应的/>，/>表示/>取最大值时对应的/>。
图8所示为基于本实施例的信号处理结果，具体为距离-多普勒处理得到的互模糊函数图。可见结合上述步骤，上行信号导致的虚假目标副峰消失，探测目标峰值被保留。通过图8所示的互模糊函数最大值即可得到目标的时延和多普勒频移信息，验证了本发明方法的可行性和有效性。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：包括如下步骤：
步骤S1、获取中心频率为、采样率为/>的5G基站直达波信号/>，目标回波信号；
步骤S2、根据5G协议生成3种本地主同步信号，并将3种本地主同步信号调制为对应的时域信号/>，考虑整数倍的最小子载波间隔/>的频偏，将时域信号/>分别与5G基站直达波信号/>进行互相关运算得到PSS互相关结果，根据PSS互相关结果/>的最大值确定SSB同步信号块时域信号的时间偏移和小区ID号/>；
步骤S3、根据时间偏移从5G基站直达波信号/>中提取4个CP-OFDM符号长度的SSB同步信号块时域信号/>，并获取SSB同步信号块时域信号/>的循环前缀部分和重复的数据部分/>，将循环前缀部分/>与重复的数据部分/>进行互相关运算得到频偏互相关结果/>，根据频偏互相关结果/>估计5G信号的频率偏移/>；
步骤S4、解码SSB同步信号块时域信号，提取SSS序列/>，根据5G协议生成336种本地辅同步信号/>，分别与SSS序列/>进行互相关运算得到SSS互相关结果/>，根据SSS互相关结果/>的最大值确定小区ID号/>；
步骤S5、根据小区ID号和小区ID号/>求取联合小区ID号/>，解码SSB同步信号块时域信号/>，根据联合小区ID号/>提取DMRS序列/>，根据5G协议生成8种本地解调参考信号/>，分别与/>进行互相关运算得到DMRS互相关结果，根据DMRS互相关结果/>的最大值确定SSB信号索引/>；
步骤S6、根据SSB传输模式确定索引为的SSB信号开始符号位置，并根据SSB信号开始符号位置与SSB信号的时间偏移/>联合确定信号帧的时间偏移/>，并通过/>和频率偏移/>调整5G基站直达波信号/>和目标回波信号/>，分别得到同步后的5G基站直达波信号/>和目标回波信号/>；
步骤S7、根据不同5G基站TDD传输模式将同步后的5G基站直达波信号中的灵活时隙和上行时隙部分置0，得到下行信号/>；
步骤S8、对下行信号和同步后的目标回波信号/>计算互模糊函数/>，根据互模糊函数/>最大值得到目标的时延/>和多普勒频移/>信息。
2.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述时间偏移和小区ID号/>的计算公式如下：
；
式中，；
其中，表示3种本地主同步信号的序号，/>表示取共轭，/>表示时延，/>表示整数倍的最小子载波间隔/>，/>表示取最大值，/>表示/>取最大值时对应的/>，表示/>取最大值时对应的/>，/>表示e的指数函数，/>表示/>取最大值时对应的整数倍，/>表示虚数单位，/>表示圆周率，/>表示时间。
3.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述5G信号的频率偏移计算公式如下：
；
式中，；
其中，表示取共轭，/>表示求和，/>表示时延，/>为角度运算符，/>为信号子载波间隔。
4.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述SSB同步信号块时域信号计算公式如下：
；
所述循环前缀部分和重复的数据部分/>计算公式如下：
；
；
其中，为采样率，/>为OFDM符号循环前缀CP短长度，/>为OFDM符号长度，为CP-OFDM符号短长度，/>为CP-OFDM符号序数，对应/>。
5.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述小区ID号计算公式如下：
；
式中，；
其中，表示离散序列的自变量，/>表示离散序列的滑动点数，/>表示SSS序列的长度，/>表示取共轭，/>表示336种本地辅同步信号的序号，/>表示取最大值，/>表示取最大值时对应的/>。
6.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述SSB信号索引计算公式如下：
；
式中，；
其中，表示离散序列的自变量，/>表示离散序列的滑动点数，/>表示DMRS序列的长度，/>表示取共轭，/>表8种本地解调参考信号的序号，/>表示取最大值，/>表示取最大值时对应的/>。
7.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述同步后的5G基站直达波信号和目标回波信号/>，计算公式如下：
；
；
式中，当SSB传输模式为Case C时，计算公式如下：
；
其中，为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为采样率。
8.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述下行信号计算公式如下：
当TDD传输模式为2.5ms单周期时，下行信号计算公式如下：
；
式中：
；
；
其中，为CP-OFDM符号长长度，/>为CP-OFDM符号短长度，/>为采样率，/>为周期数；
当TDD传输模式为2.5ms双周期时，下行信号计算公式如下：
；
式中：
；
当TDD传输模式为5ms单周期时，下行信号计算公式如下：
；
式中：
；
。
9.根据权利要求4、7、8任一项所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：
；
；
其中，为正交频分复用OFDM符号长度，/>为正交频分复用OFDM循环前缀短长度，/>为正交频分复用OFDM循环前缀长长度。
10.根据权利要求1所述的一种基于5G时隙同步的外辐射源雷达信号处理方法，其特征在于：所述目标的时延和多普勒频移/>计算公式如下：
；
式中：；
其中，表示取共轭，/>表示时延，/>表示多普勒频率，/>表示取最大值，/>表示取最大值时对应的/>，/>表示/>取最大值时对应的/>。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310962814.2A CN116660856B (zh) | 2023-08-02 | 2023-08-02 | 一种基于5g时隙同步的外辐射源雷达信号处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310962814.2A CN116660856B (zh) | 2023-08-02 | 2023-08-02 | 一种基于5g时隙同步的外辐射源雷达信号处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116660856A true CN116660856A (zh) | 2023-08-29 |
| CN116660856B CN116660856B (zh) | 2023-11-21 |
Family
ID=87721100
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310962814.2A Active CN116660856B (zh) | 2023-08-02 | 2023-08-02 | 一种基于5g时隙同步的外辐射源雷达信号处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116660856B (zh) |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3341849A (en) * | 1966-01-26 | 1967-09-12 | Bendix Corp | Self-calibrating, self-testing radio altimeter |
| JPH0375581A (ja) * | 1989-08-18 | 1991-03-29 | Mitsubishi Electric Corp | レーザレーダ装置 |
| AU4330696A (en) * | 1995-02-01 | 1996-08-08 | Nec Corporation | Multi-directional multiple access communication system |
| EP0727674A1 (en) * | 1995-02-14 | 1996-08-21 | Telefonaktiebolaget Lm Ericsson | Method and device for signal-processing in a radar system |
| WO2000049426A1 (en) * | 1999-02-17 | 2000-08-24 | Raytheon Company | Mprf interpulse phase modulation for maximizing doppler clear space |
| CA2901610A1 (en) * | 2001-11-12 | 2003-05-12 | Telstra Corporation Limited | Surface wave radar |
| WO2011046436A1 (en) * | 2009-10-14 | 2011-04-21 | Nederlandse Organisatie Voor Toegepast-Natuurwetenschappeliijk Onderzoek Tno | Fmcw distance measuring method and devices |
| WO2012027880A1 (en) * | 2010-08-30 | 2012-03-08 | Motorola Mobility, Inc. | Methods and apparatus for carrier frequency offset estimation and carrier frequency offset correction |
| CN106772309A (zh) * | 2017-03-31 | 2017-05-31 | 西安电子科技大学 | 基于通信卫星辐射源的雷达目标检测方法 |
| WO2017147439A1 (en) * | 2016-02-25 | 2017-08-31 | Cohere Technologies | Reference signal packing for wireless communications |
| CN108398676A (zh) * | 2018-05-04 | 2018-08-14 | 电子科技大学 | 一种外辐射源雷达微弱运动目标探测方法 |
| CN108680910A (zh) * | 2018-05-15 | 2018-10-19 | 河海大学 | 基于波形认知的调频广播外辐射源雷达目标检测方法 |
| CN109799484A (zh) * | 2019-01-31 | 2019-05-24 | 河海大学 | 一种外辐射源雷达系统多径抑制方法、系统及存储介质 |
| CN111650563A (zh) * | 2020-06-15 | 2020-09-11 | 桂林电子科技大学 | 外辐射源雷达同频干扰时延和能量快速估计系统及方法 |
| CN111988840A (zh) * | 2020-08-31 | 2020-11-24 | 西安电子科技大学 | 基于同步信号块联合统计的物理小区标识提取方法 |
| CN113608227A (zh) * | 2021-08-02 | 2021-11-05 | 中国科学院空天信息创新研究院 | 光子辅助雷达混频与直达波自干扰对消一体化装置及方法 |
| CN113820677A (zh) * | 2021-08-02 | 2021-12-21 | 西安电子科技大学 | 一种Link16信号外辐射源雷达的副峰抑制方法 |
| CN115412413A (zh) * | 2022-07-14 | 2022-11-29 | 南京信息工程大学 | 一种基于5g ofdm信号的外辐射源雷达杂波抑制方法 |
| CN116184346A (zh) * | 2022-12-27 | 2023-05-30 | 南京信息工程大学 | 一种5g下行信号外辐射源雷达相参积累处理方法 |
-
2023
- 2023-08-02 CN CN202310962814.2A patent/CN116660856B/zh active Active
Patent Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3341849A (en) * | 1966-01-26 | 1967-09-12 | Bendix Corp | Self-calibrating, self-testing radio altimeter |
| JPH0375581A (ja) * | 1989-08-18 | 1991-03-29 | Mitsubishi Electric Corp | レーザレーダ装置 |
| AU4330696A (en) * | 1995-02-01 | 1996-08-08 | Nec Corporation | Multi-directional multiple access communication system |
| EP0727674A1 (en) * | 1995-02-14 | 1996-08-21 | Telefonaktiebolaget Lm Ericsson | Method and device for signal-processing in a radar system |
| WO2000049426A1 (en) * | 1999-02-17 | 2000-08-24 | Raytheon Company | Mprf interpulse phase modulation for maximizing doppler clear space |
| CA2901610A1 (en) * | 2001-11-12 | 2003-05-12 | Telstra Corporation Limited | Surface wave radar |
| WO2011046436A1 (en) * | 2009-10-14 | 2011-04-21 | Nederlandse Organisatie Voor Toegepast-Natuurwetenschappeliijk Onderzoek Tno | Fmcw distance measuring method and devices |
| WO2012027880A1 (en) * | 2010-08-30 | 2012-03-08 | Motorola Mobility, Inc. | Methods and apparatus for carrier frequency offset estimation and carrier frequency offset correction |
| WO2017147439A1 (en) * | 2016-02-25 | 2017-08-31 | Cohere Technologies | Reference signal packing for wireless communications |
| CN106772309A (zh) * | 2017-03-31 | 2017-05-31 | 西安电子科技大学 | 基于通信卫星辐射源的雷达目标检测方法 |
| CN108398676A (zh) * | 2018-05-04 | 2018-08-14 | 电子科技大学 | 一种外辐射源雷达微弱运动目标探测方法 |
| CN108680910A (zh) * | 2018-05-15 | 2018-10-19 | 河海大学 | 基于波形认知的调频广播外辐射源雷达目标检测方法 |
| CN109799484A (zh) * | 2019-01-31 | 2019-05-24 | 河海大学 | 一种外辐射源雷达系统多径抑制方法、系统及存储介质 |
| CN111650563A (zh) * | 2020-06-15 | 2020-09-11 | 桂林电子科技大学 | 外辐射源雷达同频干扰时延和能量快速估计系统及方法 |
| CN111988840A (zh) * | 2020-08-31 | 2020-11-24 | 西安电子科技大学 | 基于同步信号块联合统计的物理小区标识提取方法 |
| CN113608227A (zh) * | 2021-08-02 | 2021-11-05 | 中国科学院空天信息创新研究院 | 光子辅助雷达混频与直达波自干扰对消一体化装置及方法 |
| CN113820677A (zh) * | 2021-08-02 | 2021-12-21 | 西安电子科技大学 | 一种Link16信号外辐射源雷达的副峰抑制方法 |
| CN115412413A (zh) * | 2022-07-14 | 2022-11-29 | 南京信息工程大学 | 一种基于5g ofdm信号的外辐射源雷达杂波抑制方法 |
| CN116184346A (zh) * | 2022-12-27 | 2023-05-30 | 南京信息工程大学 | 一种5g下行信号外辐射源雷达相参积累处理方法 |
Non-Patent Citations (3)
| Title |
|---|
| CHALOUPKA, Z.等: "Phase synchronization for 5G using mass market GNSS receivers", 2018 EUROPEAN FREQUENCY AND TIME FORUM (EFTF), pages 192 - 196 * |
| 王锐, 等: "基于5G基站信号的被动雷达直达波抑制技术研究", 舰船电子工程, no. 8, pages 58 - 60 * |
| 罗铱镅, 等: "基于改进的多因素SVM算法的太阳射电暴强度检测（英文）", JOURNAL OF SOUTHEAST UNIVERSITY(ENGLISH EDITION), vol. 38, no. 1, pages 20 - 26 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116660856B (zh) | 2023-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6208172B2 (ja) | 可変ガードバンドを使用する柔軟な帯域幅を達成する方法および装置 | |
| US7613104B2 (en) | Method, apparatus and computer program product providing synchronization for OFDMA downlink signal | |
| RU2365055C2 (ru) | Точное вхождение в синхронизм | |
| CN108924943B (zh) | 基于窄带物联网随机接入信道的最大相关估计检测方法 | |
| CN110531323B (zh) | 一种适用于mimo/ofdm外辐射源雷达的参考信号重构方法 | |
| WO2017191522A1 (en) | Preamble detection and time-of-arrival estimation for a single- tone frequency hopping random access preamble | |
| US20060098749A1 (en) | Apparatus and method for detecting ranging signal in an orthogonal frequency division multiple access mobile communication system | |
| US10177941B2 (en) | Method and apparatus for estimating and correcting phase error in wireless communication system | |
| KR20060022642A (ko) | 멀티 캐리어 전송을 행하는 무선 통신 시스템, 수신 장치및 수신 방법, 송신 장치 및 송신 방법, 및 지연 시간 산출장치 및 지연 시간 산출 방법 | |
| US20060215779A1 (en) | Receiving method and receiving apparatus | |
| US20130010903A1 (en) | Method and System for Diversity and Mask Matching in Channel Estimation in OFDM Communication Networks Using Circular Convolution | |
| CN101883412A (zh) | Lte在高速移动条件下初始小区搜索方法 | |
| KR20130029687A (ko) | 근거리 무선 통신 시스템에서의 수신 동기 획득 방법 및 그 장치 | |
| CN100369394C (zh) | 一种在WiMAX系统基站接收端进行帧同步定时的方法 | |
| CN115412413B (zh) | 一种基于5g ofdm信号的外辐射源雷达杂波抑制方法 | |
| KR20090084382A (ko) | 광대역 무선접속 통신 시스템에서 주파수 동기를 위한 상관장치 및 방법 | |
| CN111817814A (zh) | 一种基于ssb的5g nr信号屏蔽方法及系统 | |
| CN111147207B (zh) | 信号的发送方法及终端 | |
| CN116184346B (zh) | 一种5g下行信号外辐射源雷达相参积累处理方法 | |
| CN101895370B (zh) | 一种ofdm通信系统的干扰检测方法 | |
| Klöck et al. | LTE-signal processing for passive radar air traffic surveillance | |
| CN116660856B (zh) | 一种基于5g时隙同步的外辐射源雷达信号处理方法 | |
| CN111147405A (zh) | 一种信号的发送方法及终端 | |
| CN109738868B (zh) | 一种基于信道识别的外辐射源雷达非平稳杂波抑制方法 | |
| CN113810323B (zh) | 融合5g nr的leo卫星多普勒频偏变化率估计方法 |
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 |