CN115842567B - 一种基于chirp通信的动态门限同步方法及装置 - Google Patents
一种基于chirp通信的动态门限同步方法及装置 Download PDFInfo
- Publication number
- CN115842567B CN115842567B CN202310147545.4A CN202310147545A CN115842567B CN 115842567 B CN115842567 B CN 115842567B CN 202310147545 A CN202310147545 A CN 202310147545A CN 115842567 B CN115842567 B CN 115842567B
- Authority
- CN
- China
- Prior art keywords
- condition
- value
- peak
- max
- synchronization
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000004891 communication Methods 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000005070 sampling Methods 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 4
- 238000012887 quadratic function Methods 0.000 claims description 4
- 230000007480 spreading Effects 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 238000012886 linear function Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 7
- 238000003775 Density Functional Theory Methods 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- 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
- Synchronisation In Digital Transmission Systems (AREA)
Abstract
本发明公开了一种基于CHIRP通信的动态门限同步方法及装置，针对现有技术中动态门限实现较为复杂，寻找一个无线通信***要兼容其所有应用场景的动态门限计算方法比较耗费时间。且在信噪比很低时，动态门限算法仍然存在一定的虚警率和漏警率。本发明的技术方案包括：首先在时域完成接收基带IQ信号与基频信号的除法运算，然后转到频域根据快速傅立叶变换的结果来完成解调门限计算、preamble搜索、帧同步、ID解调等的一些列后续处理。具有运算过程简洁、逻辑易于实现及资源占用率低等优点。
Description
技术领域
本发明属于物联网通信技术领域，具体涉及一种基于CHIRP通信的动态门限同步方法及装置。
背景技术
在无线通信中，想要正确解调接收到的数据包最重要也是最关键的一步就是同步，一种性能优异的同步算法不仅可以保证解调的正确性还可以提升性能，而同步算法中比较关键的便是门限的确定，门限的好坏直接关系同步的虚警率和漏警率，进而决定整个无线通信***灵敏度指标是否符合要求。因此同步算法的研究也是一直以来备受关注的方向。
在现有的同步技术中，门限主要分两类：1.固定门限，通过算法建模/仿真以及实际测试来确定最终应用的固定门限值。固定门限有实现简单、资源消耗低、易于获取的优点，但其对于大信号和小信号等不同的场景下虚警率和漏警率存在较大差异，相同条件下无法达到最优的性能。
2.动态门限，通过算法实时动态地计算门限，可在大信号和小信号等不同场景下的达到最优的解调性能。其存在的缺点是实现较为复杂，寻找一个无线通信***要兼容其所有应用场景的动态门限计算方法比较耗费时间。且在信噪比很低时，动态门限算法仍然存在一定的虚警率和漏警率。
发明内容
针对上述现有技术中存在的问题，本发明提出了一种基于CHIRP通信的动态门限同步方法及装置，其目的为：通过动态门限的计算方式，降低虚警率和漏警率，同时灵活调整因子可实现接近于期望值的灵敏度性能。
为实现上述目的本发明所采用的技术方案是：提供一种基于CHIRP通信的动态门限同步方法，包括：
S1：数模转化器对获取的初始数据进行处理，得到基带IQ信号，并将所述基带IQ信号输出到基带滤波器组；
S2：所述基带滤波器组根据输入的基带IQ信号的信号带宽，实施对应信号带宽的滤波和下采样，完成基带IQ信号采样率的转换，得到转换信号；
S3：将转换信号根据不同的扩频因子SF配置的点数分组，将不同的分组与基频信号进行除法运算；
S4：将除法运算得到的结果按照扩频因子SF配置进行可变点数的快速傅里叶变换；将快速傅里叶变换得到的结果用CORDIC函数求模，得到求模数据；
S5:找出所述求模数据的峰峰值位置及其数值，通过峰峰值数值位置及其数值计算动态门限；
S6：利用动态门限与峰峰值数值数值进行初步匹配同步，然后通过峰峰值位置进一步确认是否完成同步。
较优的，本发明S2中，基带滤波器组分为四级，分别为：一级2倍下采样滤波、两级4倍下采样滤波和一级通道滤波，基带IQ信号输入到基带滤波器组后，基带滤波器组根据不同信号带宽经过的基带滤波器组的路径，实施对应的滤波和下采样，完成基带IQ信号采样率的转换。
较优的，本发明S3具体为：
S3.1:在CSS中使用chirp符号的波形为：
其中，是一个关于时间的二次函数，，a和b是系数；
S3.2：定义chirp符号的瞬时频率为：
其中，是对时间t求导后的一次函数，a和b是系数；
因此，将chirp符号的瞬时频率特性模型化为为f0=at+b，当b=0时，为基频信号；
S3.3：其中chirp符号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M-1/2*m))]
其中，data_css 是产生的chirp符号数据，M为2^SF，m的取值为0～M-1；
S3.4：基频信号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M))]
其中，M为2^SF，m的取值为0～M-1；
S3.5：转换信号根据不同扩频因子SF以2^SF个采样点为单位分组，每组数据与本地产生的基频信号进行除法运算。
较优的，本发明S5中计算最终门限具体为：
S5.1：对整个分组的求模数据进行求和，结果为sum：
其中，SF为扩频因子，Dn为n点求模数据中各个采样点的数值；
S5.2：搜索获取峰峰值的数值，记为max_value；
S5.3：搜索获取峰峰值位置前面两个点的数值的平均值，记为max_p1；
S5.4：搜索获取峰峰值位置后面两个点的数值的平均值，记为max_p2；
S5.5：将步骤S5.1～S5.4的计算结果进行运算并求取平均，公式为：
Mean_thld=(sum - max_value - max_p1 - max_p2)/2^SF
S5.6：将求取平均得到的结果乘以系数F1得到动态门限，公式为：
Final_thld= Mean_thld * F1。
较优的，本发明S6中，进行初步匹配同步的具体步骤为：
S6.1：计算判决值1，记为judge1：
judge1= (（max_value + max_p1 + max_p2）/2^SF) * F2
S6.2：计算判决值2，记为judge2：
judge2= (max_value/2^SF) * F2
S6.3：计算判决值3，记为judge3：
judge3= (max_p1/2^SF )* F2
S6.4：计算判决值4，记为judge4：
judge4= (max_p2/2^SF) * F2
其中，F2为系数；
初步同步条件具体为：
S6.5：判断条件一是否匹配成功，若judge1>Final_thld则条件一满足，记为condition1=1，否则condition1=0；
S6.6：判断条件二是否匹配成功，若judge2>Final_thld则条件二满足，记为condition2=1，否则condition2=0；
S6.7：判断条件三是否匹配成功，若judge3>Final_thld则条件三满足，记为condition3=1，否则condition3=0；
S6.8：判断条件四是否匹配成功，若judge4>Final_thld则条件四满足，记为condition4=1，否则condition4=0；
S6.9：判断条件五是否匹配成功，若(max_value/2^SF)>K则条件五满足，记为condition5=1，否则condition5=0；
其中，K为系数；
S6.10：判断是否有同步产生，若condition1&&（condition2 || condition3）&&condition4&&condition5 结果为1，则完成初步匹配同步。
较优的，本发明S6中利用峰峰值位置进一步确定是否完成同步的具体步骤为：
S6.11：保存第一次得到求模数据后搜索的峰峰值位置，记为position(1)；
S6.12：保存第二次得到求模数据后搜索的峰峰值位置，记为position(2),若本次求模数据中2^SF点未满足初步同步条件则回到S6.11重新搜索；
S6.13：保存第n次得到求模数据后搜索的峰峰值位置，记为position(n),n为自然数，若本次求模数据中2^SF点未满足初步同步条件则回到S6.11重新搜索；
S6.14：步骤S6.11～S6.13以流水的形式持续搜索并保存峰峰值位置；若连续m次搜索的峰峰值位置均在某个位置的前后各一个点范围之内则判定为符合同步条件，其中m为系数；
S6.15：若步骤S6.10和步骤S6.14的条件同时成立即可判定同步完成，然后进行后续解调处理，若步骤S6.10和步骤S6.14的条件不能同时成立，则回到步骤S5，持续搜索峰峰值位置及其数值。
本发明还提出了一种基于CHIRP通信的动态门限同步装置，包括：
数模转换器：获取初始数据并进行预处理，输出基带IQ信号；
基带滤波器组：输入输出基带IQ信号，根据输入的基带IQ信号的信号带宽，实施对应信号带宽的滤波和下采样，完成基带IQ信号采样率的转换，得到转换信号；
除法运算单元：将转换信号根据不同的扩频因子SF配置的点数分组，将不同的分组与基频信号进行除法运算；
DFT求模单元：将除法运算得到的结果按照扩频因子SF配置进行可变点数的快速傅里叶变换；将快速傅里叶变换得到的结果用CORDIC函数求模，得到求模数据；
动态门限单元:找出所述求模数据的峰峰值位置及其数值，通过峰峰值数值位置及其数值计算动态门限；
同步单元：利用动态门限与峰峰值数值数值进行初步匹配同步，然后通过峰峰值位置进一步确认是否完成同步。
相比现有技术，本发明的技术方案具有如下优点/有益效果：
1.本发明通过在时域完成接收基带IQ信号与基频信号的除法运算，然后转到频域根据快速傅立叶变换的结果来完成解调门限计算、preamble搜索、帧同步、ID解调等的一些列后续处理，处理具有运算过程简洁、逻辑易于实现及资源占用率低等特点。
2.本发明通过灵活调整参数F1、参数F2和参数K，可实现接近于期望值的灵敏度性能。
附图说明
为了更清楚地说明本发明实施方式的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1是本发明实施例1的流程示意图。
图2是本发明实施例1的帧结构示意图。
图3是本发明实施例1的动态门限计算示意图。
具体实施方式
为使本发明目的、技术方案和优点更加清楚，下面对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式是本发明的一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。因此，以下提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。
实施例1：
如图1-图3所示，本实施例1提供一种基于CHIRP通信的动态门限同步方法，同步过程中所使用到的符号为preamble，图2所示为本发明***的帧结构。帧结构中preamble数量由上层协议配置，但至少为4个。记SF为10，以下为具体的实施步骤：
S1：数模转化器对获取的初始数据进行处理，得到基带IQ信号，并将所述基带IQ信号输出到基带滤波器组；
S2：所述基带滤波器组根据输入的基带IQ信号的信号带宽，实施对应信号带宽的滤波和下采样，完成基带IQ信号采样率的转换，得到转换信号；基带滤波器组分为四级，分别为：一级2倍下采样滤波、两级4倍下采样滤波和一级通道滤波，基带IQ信号输入到基带滤波器组后，基带滤波器组根据不同信号带宽经过的基带滤波器组的路径，实施对应的滤波和下采样，完成基带IQ信号采样率的转换。
S3：将转换信号以1024为单位分组，将不同的分组与基频信号进行除法运算；S3具体为：
S3.1:在CSS中使用chirp符号的波形为：
其中，是一个关于时间的二次函数，；
S3.2：定义chirp符号的瞬时频率为：
因此，将chirp符号的瞬时频率特性模型化为为f0=at+b，因此解调也就是找到每个chirp符号对应的b，方法就是将接收信号与本地产生的chirp基频信号相除，当b=0时，为基频信号；相除后得到结果即为频率为b的单音信号与信道畸变叠加的结果。然后通过DFT计算找到单音信号的频点，从而得到b的数值。
S3.3：其中chirp符号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M-1/2*m))]
其中，data_css 是产生的 chirp 符号数据，M为2^SF，m的取值范围为0～M-1；
S3.4：基频信号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M))]
其中，M为2^SF，m的取值为0～M-1；
S3.5：转换信号根据不同扩频因子SF以2^SF个采样点为单位分组，每组数据与本地产生的基频信号共轭相乘，
data_demod = rec_data*conj(basic_data)。
S4：将除法运算得到的结果按照扩频因子SF配置进行可变点数的快速傅里叶变换；将快速傅里叶变换得到的结果用CORDIC函数求模，即将除法运算得到的结果进行DFT运算和abs求模运算，得到求模数据；
S5:找出所述求模数据的峰峰值位置及其数值，通过峰峰值数值位置及其数值计算动态门限；计算最终门限具体为：
S5.1：对整个分组的求模数据进行求和，结果为sum：
其中，SF为扩频因子，Dn为n点求模数据中各个采样点的数值；
S5.2：搜索获取峰峰值的数值，记为max_value；
S5.3：搜索获取峰峰值位置前面两个点的数值的平均值，记为max_p1；
S5.4：搜索获取峰峰值位置后面两个点的数值的平均值，记为max_p2；
S5.5：将步骤S5.1～S5.4的计算结果进行运算并求取平均，公式为：
Mean_thld=(sum - max_value - max_p1 - max_p2)/2^SF
S5.6：将求取平均得到的结果乘以系数F1得到动态门限，公式为：
Final_thld= Mean_thld * F1。
S6：利用动态门限与峰峰值数值数值进行初步匹配同步，然后通过峰峰值位置进一步确认是否完成同步。进行初步匹配同步的具体步骤为：
S6.1：计算判决值1，记为judge1：
judge1= (（max_value + max_p1 + max_p2）/2^SF) * F2
S6.2：计算判决值2，记为judge2：
judge2= (max_value/2^SF) * F2
S6.3：计算判决值3，记为judge3：
judge3= (max_p1/2^SF )* F2
S6.4：计算判决值4，记为judge4：
judge4= (max_p2/2^SF) * F2
其中，F2为系数；
初步同步条件具体为：
S6.5：判断条件一是否匹配成功，若judge1>Final_thld则条件一满足，记为condition1=1，否则condition1=0；
S6.6：判断条件二是否匹配成功，若judge2>Final_thld则条件二满足，记为condition2=1，否则condition2=0；
S6.7：判断条件三是否匹配成功，若judge3>Final_thld则条件三满足，记为condition3=1，否则condition3=0；
S6.8：判断条件四是否匹配成功，若judge4>Final_thld则条件四满足，记为condition4=1，否则condition4=0；
S6.9：判断条件五是否匹配成功，若(max_value/2^SF)>K则条件五满足，记为condition5=1，否则condition5=0；
其中，K为系数；
S6.10：判断是否有同步产生，若condition1&&（condition2 || condition3）&&condition4&&condition5 结果为1，则完成初步匹配同步。
利用峰峰值位置进一步确定是否完成同步的具体步骤为：
S6.11：保存第一次得到求模数据后搜索的峰峰值位置，记为position(1)；
S6.12：保存第二次得到求模数据后搜索的峰峰值位置，记为position(2),若本次求模数据中2^SF点未满足初步同步条件则回到S6.11重新搜索；
S6.13：保存第n次得到求模数据后搜索的峰峰值位置，记为position(n),n为自然数，若本次求模数据中2^SF点未满足初步同步条件则回到S6.11重新搜索；
S6.14：步骤S6.11～S6.13以流水的形式持续搜索并保存峰峰值位置；若连续m次搜索的峰峰值位置均在某个位置的前后各一个点范围之内则判定为符合同步条件，其中m为系数；
S6.15：若步骤S6.10和步骤S6.14的条件同时成立即可判定同步完成，然后进行后续解调处理，若步骤S6.10和步骤S6.14的条件不能同时成立，则回到步骤S5，持续搜索峰峰值位置及其数值。
本实施例1中，以流水的形式持续搜索求模数据的最大值及位置，位置记为P，同时持续计算动态门限，动态门限计算的流程框图如图3所示,计算过程见S5.1-S5.6,同步过程见S6.1-S6.15。若连续3次最大值大于动态门限，且连续3次最大值的位置满足|P1-P2|<=1,并且|P2-P3|<=1，则同步完成，进行后续解调处理。
本实施例1还提出了一种基于CHIRP通信的动态门限同步装置，包括：
数模转换器：获取初始数据并进行预处理，输出基带IQ信号；
基带滤波器组：输入输出基带IQ信号，根据输入的基带IQ信号的信号带宽，实施对应信号带宽的滤波和下采样，完成基带IQ信号采样率的转换，得到转换信号；
除法运算单元：将转换信号根据不同的扩频因子SF配置的点数分组，将不同的分组与基频信号进行除法运算；
DFT求模单元：将除法运算得到的结果按照扩频因子SF配置进行可变点数的快速傅里叶变换；将快速傅里叶变换得到的结果用CORDIC函数求模，得到求模数据；
动态门限单元:找出所述求模数据的峰峰值位置及其数值，通过峰峰值数值位置及其数值计算动态门限；
同步单元：利用动态门限与峰峰值数值数值进行初步匹配同步，然后通过峰峰值位置进一步确认是否完成同步。
以上仅是本发明的优选实施方式，应当指出的是，上述优选实施方式不应视为对本发明的限制，本发明的保护范围应当以权利要求所限定的范围为准。对于本技术领域的普通技术人员来说，在不脱离本发明的精神和范围内，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (3)
1.一种基于CHIRP通信的动态门限同步方法，其特征在于，包括：
S1：数模转化器对获取的初始数据进行处理，得到基带IQ信号，并将所述基带IQ信号输出到基带滤波器组；
S2：所述基带滤波器组根据输入的基带IQ信号的信号带宽，实施对应信号带宽的滤波和下采样，完成基带IQ信号采样率的转换，得到转换信号；
S3：将转换信号根据不同的扩频因子SF配置的点数分组，将不同的分组与基频信号进行除法运算；具体为：
S3.1:在CSS中使用chirp符号的波形为：
其中，是一个关于时间的二次函数，，a和b是系数；
S3.2：定义chirp符号的瞬时频率为：
其中，是对时间t求导后的一次函数，a和b是系数；
因此，将chirp符号的瞬时频率特性模型化为f0=at+b，当b=0时，为基频信号；
S3.3：其中chirp符号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M-1/2*m))]
其中，data_css 是产生的chirp符号数据，M为2^SF，m的取值为0～M-1；
S3.4：基频信号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M))]
其中，M为2^SF，m的取值范围为0～M-1；
S3.5：转换信号根据不同扩频因子SF以2^SF个采样点为单位分组，每组数据与本地产生的基频信号进行除法运算；
S4：将除法运算得到的结果按照扩频因子SF配置进行可变点数的快速傅里叶变换；将快速傅里叶变换得到的结果用CORDIC函数求模，得到求模数据；
S5:找出所述求模数据的峰峰值位置及其数值，通过峰峰值数值位置及其数值计算动态门限；计算最终门限具体为：
S5.1：对整个分组的求模数据进行求和，结果为sum：
其中，SF为扩频因子，Dn为n点求模数据中各个采样点的数值；
S5.2：搜索获取峰峰值的数值，记为max_value；
S5.3：搜索获取峰峰值位置前面两个点的数值的平均值，记为max_p1；
S5.4：搜索获取峰峰值位置后面两个点的数值的平均值，记为max_p2；
S5.5：将步骤S5.1～S5.4的计算结果进行运算并求取平均，公式为：
Mean_thld=(sum - max_value - max_p1 - max_p2)/2^SF
S5.6：将求取平均得到的结果乘以系数F1得到动态门限，公式为：
Final_thld= Mean_thld * F1；
S6：利用动态门限与峰峰值数值进行初步匹配同步，然后通过峰峰值位置进一步确认是否完成同步；进行初步匹配同步的具体步骤为：
S6.1：计算判决值1，记为judge1：
judge1= (（max_value + max_p1 + max_p2）/2^SF) * F2
S6.2：计算判决值2，记为judge2：
judge2= (max_value/2^SF) * F2
S6.3：计算判决值3，记为judge3：
judge3= (max_p1/2^SF )* F2
S6.4：计算判决值4，记为judge4：
judge4= (max_p2/2^SF) * F2
其中，F2为系数；
初步同步条件具体为：
S6.5：判断条件一是否匹配成功，若judge1 > Final_thld则条件一满足，记为condition1=1，否则condition1=0；
S6.6：判断条件二是否匹配成功，若judge2 > Final_thld则条件二满足，记为condition2=1，否则condition2=0；
S6.7：判断条件三是否匹配成功，若judge3 > Final_thld则条件三满足，记为condition3=1，否则condition3=0；
S6.8：判断条件四是否匹配成功，若judge4 > Final_thld则条件四满足，记为condition4=1，否则condition4=0；
S6.9：判断条件五是否匹配成功，若(max_value/2^SF) > K则条件五满足，记为condition5=1，否则condition5=0；
其中，K为系数；
S6.10：判断是否有同步产生，若condition1 && （condition2 || condition3）&&condition4 && condition5 结果为1，则完成初步匹配同步；
利用峰峰值位置进一步确定是否完成同步的具体步骤为：
S6.11：保存第一次得到求模数据后搜索的峰峰值位置，记为position(1)；
S6.12：保存第二次得到求模数据后搜索的峰峰值位置，记为position(2),若本次求模数据中2^SF点未满足初步同步条件则回到S6.11重新搜索；
S6.13：保存第n次得到求模数据后搜索的峰峰值位置，记为position(n),n为自然数，若本次求模数据中2^SF点未满足初步同步条件则回到S6.11重新搜索；
S6.14：步骤S6.11～S6.13以流水的形式持续搜索并保存峰峰值位置；若连续m次搜索的峰峰值位置均在某个位置的前后各一个点范围之内则判定为符合同步条件，其中m为系数；
S6.15：若步骤S6.10和步骤S6.14的条件同时成立即可判定同步完成，然后进行后续解调处理，若步骤S6.10和步骤S6.14的条件不能同时成立，则回到步骤S5，持续搜索峰峰值位置及其数值。
2.根据权利要求1所述的一种基于CHIRP通信的动态门限同步方法，其特征在于，S2中，基带滤波器组分为四级，分别为：一级2倍下采样滤波、两级4倍下采样滤波和一级通道滤波，基带IQ信号输入到基带滤波器组后，基带滤波器组根据不同信号带宽经过的基带滤波器组的路径，实施对应的滤波和下采样，完成基带IQ信号采样率的转换。
3.一种基于CHIRP通信的动态门限同步装置，其特征在于，包括：
数模转换器：获取初始数据并进行预处理，输出基带IQ信号；
基带滤波器组：输入输出基带IQ信号，根据输入的基带IQ信号的信号带宽，实施对应信号带宽的滤波和下采样，完成基带IQ信号采样率的转换，得到转换信号；
除法运算单元：将转换信号根据不同的扩频因子SF配置的点数分组，将不同的分组与基频信号进行除法运算；具体为：
步骤3.1:在CSS中使用chirp符号的波形为：
其中，是一个关于时间的二次函数，，a和b是系数；
步骤3.2：定义chirp符号的瞬时频率为：
其中，是对时间t求导后的一次函数，a和b是系数；
因此，将chirp符号的瞬时频率特性模型化为f0=at+b，当b=0时，为基频信号；
步骤3.3：其中chirp符号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M-1/2*m))]
其中，data_css 是产生的chirp符号数据，M为2^SF，m的取值为0～M-1；
步骤3.4：基频信号产生的matlab简化模型为：
data_css(m+1) = [exp(j*2*pi*(m^2/2/M))]
其中，M为2^SF，m的取值范围为0～M-1；
步骤3.5：转换信号根据不同扩频因子SF以2^SF个采样点为单位分组，每组数据与本地产生的基频信号进行除法运算；
DFT求模单元：将除法运算得到的结果按照扩频因子SF配置进行可变点数的快速傅里叶变换；将快速傅里叶变换得到的结果用CORDIC函数求模，得到求模数据；
动态门限单元:找出所述求模数据的峰峰值位置及其数值，通过峰峰值数值位置及其数值计算动态门限；计算最终门限具体为：
步骤5.1：对整个分组的求模数据进行求和，结果为sum：
其中，SF为扩频因子，Dn为n点求模数据中各个采样点的数值；
步骤5.2：搜索获取峰峰值的数值，记为max_value；
步骤5.3：搜索获取峰峰值位置前面两个点的数值的平均值，记为max_p1；
步骤5.4：搜索获取峰峰值位置后面两个点的数值的平均值，记为max_p2；
步骤5.5：将步骤5.1～步骤5.4的计算结果进行运算并求取平均，公式为：
Mean_thld=(sum - max_value - max_p1 - max_p2)/2^SF
步骤5.6：将求取平均得到的结果乘以系数F1得到动态门限，公式为：
Final_thld= Mean_thld * F1；
同步单元：利用动态门限与峰峰值数值进行初步匹配同步，然后通过峰峰值位置进一步确认是否完成同步；进行初步匹配同步的具体步骤为：
步骤6.1：计算判决值1，记为judge1：
judge1= (（max_value + max_p1 + max_p2）/2^SF) * F2
步骤6.2：计算判决值2，记为judge2：
judge2= (max_value/2^SF) * F2
步骤6.3：计算判决值3，记为judge3：
judge3= (max_p1/2^SF )* F2
步骤6.4：计算判决值4，记为judge4：
judge4= (max_p2/2^SF) * F2
其中，F2为系数；
初步同步条件具体为：
步骤6.5：判断条件一是否匹配成功，若judge1 > Final_thld则条件一满足，记为condition1=1，否则condition1=0；
步骤6.6：判断条件二是否匹配成功，若judge2 > Final_thld则条件二满足，记为condition2=1，否则condition2=0；
步骤6.7：判断条件三是否匹配成功，若judge3 > Final_thld则条件三满足，记为condition3=1，否则condition3=0；
步骤6.8：判断条件四是否匹配成功，若judge4 > Final_thld则条件四满足，记为condition4=1，否则condition4=0；
步骤6.9：判断条件五是否匹配成功，若(max_value/2^SF) > K则条件五满足，记为condition5=1，否则condition5=0；
其中，K为系数；
步骤6.10：判断是否有同步产生，若condition1 && （condition2 || condition3）&&condition4 && condition5 结果为1，则完成初步匹配同步；
利用峰峰值位置进一步确定是否完成同步的具体步骤为：
步骤6.11：保存第一次得到求模数据后搜索的峰峰值位置，记为position(1)；
步骤6.12：保存第二次得到求模数据后搜索的峰峰值位置，记为position(2),若本次求模数据中2^SF点未满足初步同步条件则回到步骤6.11重新搜索；
步骤6.13：保存第n次得到求模数据后搜索的峰峰值位置，记为position(n),n为自然数，若本次求模数据中2^SF点未满足初步同步条件则回到步骤6.11重新搜索；
步骤6.14：步骤6.11～步骤6.13以流水的形式持续搜索并保存峰峰值位置；若连续m次搜索的峰峰值位置均在某个位置的前后各一个点范围之内则判定为符合同步条件，其中m为系数；
步骤6.15：若步骤6.10和步骤6.14的条件同时成立即可判定同步完成，然后进行后续解调处理，若步骤6.10和步骤6.14的条件不能同时成立，则回到动态门限单元，持续搜索峰峰值位置及其数值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310147545.4A CN115842567B (zh) | 2023-02-22 | 2023-02-22 | 一种基于chirp通信的动态门限同步方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310147545.4A CN115842567B (zh) | 2023-02-22 | 2023-02-22 | 一种基于chirp通信的动态门限同步方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115842567A CN115842567A (zh) | 2023-03-24 |
| CN115842567B true CN115842567B (zh) | 2023-05-02 |
Family
ID=85580047
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310147545.4A Active CN115842567B (zh) | 2023-02-22 | 2023-02-22 | 一种基于chirp通信的动态门限同步方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115842567B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114257270A (zh) * | 2021-12-24 | 2022-03-29 | 浙江大学 | 高灵敏星间扩频通信***与快速捕获方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001268066A (ja) * | 2000-03-23 | 2001-09-28 | Matsushita Electric Ind Co Ltd | 同期検出装置及び同期検出方法 |
| US6961364B1 (en) * | 2000-04-18 | 2005-11-01 | Flarion Technologies, Inc. | Base station identification in orthogonal frequency division multiplexing based spread spectrum multiple access systems |
| US7746939B2 (en) * | 2005-07-29 | 2010-06-29 | Itt Manufacturing Enterprises, Inc. | Methods and apparatus for encoding information in a signal by spectral notch modulation |
| CN101702704A (zh) * | 2009-11-16 | 2010-05-05 | 清华大学 | 时域同步的多载波扩频收发方法、装置及*** |
| CN103364802A (zh) * | 2013-07-12 | 2013-10-23 | 北京大学 | 一种基于双通道联合利用的tddm扩频导航信号捕获方法 |
| CN107769816B (zh) * | 2017-11-01 | 2020-10-09 | 中山大学花都产业科技研究院 | 一种Chirp扩频通信***接收机时间同步***及方法 |
| CN114006649A (zh) * | 2021-10-29 | 2022-02-01 | 重庆两江卫星移动通信有限公司 | 一种卫星通信双chirp信号的判决方法 |
-
2023
- 2023-02-22 CN CN202310147545.4A patent/CN115842567B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114257270A (zh) * | 2021-12-24 | 2022-03-29 | 浙江大学 | 高灵敏星间扩频通信***与快速捕获方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115842567A (zh) | 2023-03-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104125190B (zh) | 适于低信噪比信道环境的ofdm***符号定时同步实现方法 | |
| CN103475621B (zh) | 一种多载波同步***和同步方法 | |
| CN110191071B (zh) | 一种窄带物联网***中基于信道估计的测量方法及装置 | |
| CN107911329B (zh) | 一种信号分析仪ofdm信号解调方法 | |
| CN111132272A (zh) | 一种5g nr***小区搜索定时同步方法 | |
| CN101997807B (zh) | 一种信道估计方法及装置 | |
| CN101232472A (zh) | 一种导频与数据混合叠加的ofdm信道检测方法 | |
| CN103036820B (zh) | 一种基于参考信号的多小区信道估计方法和装置 | |
| CN104363194A (zh) | 基于波形变换的psk调制识别方法 | |
| CN115952434A (zh) | 一种面向信号识别深度学习模型训练的数据增强方法 | |
| CN112910816A (zh) | 基于相位加权共轭对称的滤波器组多载波定时同步方法 | |
| CN107247276A (zh) | 基于重叠多块补零算法的弱信号捕获方法 | |
| CN112738000B (zh) | 一种pss分块互相关检测方法 | |
| CN108881084A (zh) | 一种基于gp分布的bpsk/qpsk信号识别方法 | |
| CN107181558A (zh) | 一种基于lte***的邻区检测方法及装置 | |
| CN115842567B (zh) | 一种基于chirp通信的动态门限同步方法及装置 | |
| CN109104215A (zh) | 一种基于小波变换的跳频信号码元速率盲估计方法 | |
| CN112714448A (zh) | 动态频谱共享中的主同步信号检测方法 | |
| CN111611686A (zh) | 一种通信信号时频域的检测方法 | |
| CN106330792A (zh) | 一种基于dft变换域的带内滤噪信道估计算法及*** | |
| CN114172771A (zh) | 一种基于突发通信的快速频偏估计和补偿的方法及*** | |
| CN113938368A (zh) | 基于fpga的lte-v收端同步的方法及*** | |
| CN106453165B (zh) | 一种dft－base信道估计方法和装置 | |
| CN101984696B (zh) | 无线局域网中不同格式帧的检测方法及检测装置 | |
| CN114071705B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |