Communication signal automatic detection and monitoring method and device based on frequency domain threshold progressive segmentation
阅读说明：本技术 基于频域阈值递进分割的通信信号自动检测与监视方法及装置 (Communication signal automatic detection and monitoring method and device based on frequency domain threshold progressive segmentation ) 是由 杨蛟龙 屈会鹏 孙亮亮 杨文洁 任颖 教富龙 吕波 张玉峰 于 2019-10-25 设计创作，主要内容包括：本申请实施例公开了一种基于频域阈值递进分割的通信信号自动检测与监视方法及装置,方法包括：定义初始化功率阈值对带宽内的频点进行分割,在低于功率阈值的频点集合中利用直方图法统计获取底噪分布区域,通过计算底噪分布区域中的各频点对应功率的一阶中心矩,估计底噪功率；定义功率阈值变化范围,在范围内从大至小,逐步递减阈值,对带宽内的频点进行分割；比较相邻两次分割得到的频率集合和功率值,更新频点集；根据阈值逐步递减得到的最终频点分割结果,提取各个载波并分别进行载波信号的关键特征参数估计。具有自动化、低成本、易实现、实时性好的特点,适用于多监测点布署的场景,扩展了通信信号载波监视系统的实现方式。(The embodiment of the application discloses a communication signal automatic detection and monitoring method and device based on frequency domain threshold progressive segmentation, wherein the method comprises the following steps: defining an initialization power threshold value to divide frequency points in a bandwidth, counting and acquiring a background noise distribution area in a frequency point set lower than the power threshold value by using a histogram method, and estimating background noise power by calculating a first-order central moment of power corresponding to each frequency point in the background noise distribution area; defining a power threshold variation range, gradually decreasing the threshold from large to small in the range, and segmenting frequency points in the bandwidth; comparing the frequency set and the power value obtained by two adjacent partitions, and updating the frequency point set; and extracting each carrier according to a final frequency point segmentation result obtained by gradually decreasing the threshold value and respectively estimating the key characteristic parameters of the carrier signals. The method has the characteristics of automation, low cost, easy realization and good real-time performance, is suitable for a scene of deployment of multiple monitoring points, and expands the realization mode of a communication signal carrier monitoring system.)
1. A method for automatic detection and monitoring of a communication signal based on progressive division by frequency domain thresholds, the method comprising:
acquiring a data frame to be analyzed, and dividing the frequency domain of the data frame into frequency points (f)n，Pn) Expressed by way of set PSD, wherein PSD { (f)n，Pn)}，fnFrequency of n +1 point, PnFor the measured power within the resolution bandwidth at that frequency, N is 0, 1, 2, …, N;
searching for a power value P within a set PSDnDetermining the maximum value PmaxAnd a minimum value PminWhen (P)max-Pmin) And when the data frame to be analyzed contains the communication signal, determining that the data frame to be analyzed contains the communication signal when the data frame is more than or equal to 2 dB.
2. The method of claim 1, wherein after determining that the data frame to be analyzed contains a communication signal, the method further comprises:
initialization power threshold Pthn0＝(Pmax-Pmin)/4+PminBy the use of Pthn0Dividing frequency points in the PSD set to obtain a set NSI { (f)n，Pn) In which P isn≤Pth0，(fn，Pn) E is the PSD; performing histogram statistics on power values of frequency points in the NSI set according to a preset gear, determining the first gear with the most distributed points, and obtaining the first gear starting and stopping power with the most distributed points as PnstAnd Pnsp；
According to PnstAnd PnspSearching frequency points in the PSD set to obtain a set NSF { (f)n，Pn)}，Pnst≤Pn≤Pnsp，(fn，Pn) E, PSD, calculating a first central moment E (NSF) of the power value in the NSF set, namely considering the background noise power level P of the frame frequency domain datanoise＝E(NSF)；
Reinitializing the initial power threshold Pthf＝Pmax-2, with PthfFor frequency point (f) in the setx，Px) Dividing the image into sets LP { (fp)x，Ppx) And RP { (fp)x，Ppx) In which Ppx≤Pthf≤Ppx+1，(fpx，Ppx)∈PSD；Ppx≥Pthf≥Ppx+1，(fpx，Ppx) E is PSD, x is 0, 1, 2, …, N-1; reserving frequency points which are paired in LP and RP, wherein the paired points are corresponding to a point in one RP, and the two points are respectively corresponding to the intersection of the left roll-off edge and the right roll-off edge of the same carrier signal and the threshold powerPoint, searching M to form point pairs;
taking M pairs of point pairs as elements in LP and RP respectively, corresponding M carriers to each pair of point pair, and using bandwidth boundary frequency fpl on the left side of the carrierxRight side bandwidth boundary frequency fprxThe combination defines one carrier, and the sets LP and RP will constitute the set of carrier boundaries { (fpl)x，fprx)}，Pplx≤Pthf≤Pplx+1，Pprx≥Pthf≥Pprx+1，(fplx，Pplx)∈LP，(fprx，Pprx) E.g., RP, the carriers formed by all elements in the carrier boundary set are confirmed one by one, element (fpl)k，fprk) Determined carrier region PCarrier { (f)n，Pn) Wherein, fplk≤fn≤fprk，(fn，Pn) Belongs to PSD, and searches power value P in set PCarriernFinding the maximum value PPCmaxAnd minimum value PPCminWhen (PPC)max-PPCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier valid carrier boundaries form a new carrier boundary set carrier border.
3. The method of claim 2, further comprising:
will power threshold PthfDecreasing by steps of-1 dB, i.e. Pthf_New＝Pthf-1, using a decreasing threshold value PthfA New pair of frequency points (f) in the set PSDn，Pn) The division is carried out again to obtain a carrier left side bandwidth boundary point set LN and a left side bandwidth boundary point set RN, and the sets LN and RN form a new carrier set NCarrirBorder { (fnl)x，fnrx)}；Pnlx≤Pthf_New≤Pnlx+1，Pnrx≥Pthf_New≥Pnrx+1，(fnlx，Pnlx)∈LN，(fnrx，Pnrx) E.g. RN, confirming the carriers formed by all elements in the set one by one, element (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) Belongs to PSD, and searches for power value P in set NCarriernTo obtain the maximum value PNCmaxAnd minimum value PNCminWhen (PNC)max-PNCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier-valid carrier boundaries form a new set of carrier boundaries NCarrierBorder.
4. The method of claim 3, further comprising:
for the sets Carrier Border and NCarrier Border, the carrier boundary in NCarrier Border (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) E PSD, if NCarrier does not contain a carrier region determined by any carrier boundary in Carrier Border, then this carrier boundary (fnl)k，fnrk) Adding into a New carrier boundary carrier _ New set; if the NCarrier includes a carrier region defined by a carrier boundary in the carrier border, it is judged that the NCarrier includes several carrier regions defined by carrier boundaries in the carrier border, and if it includes only one, the carrier border is determined (fnl)k，fnrk) Adding into a Carrier Border _ New set, if two or more, adding the contained multiple carrier boundaries into the Carrier Border _ New set;
the carrier border set CarrierBorder is updated so that CarrierBorder is CarrierBorder _ New.
5. The method of claim 4, further comprising:
when power threshold value Pthf＝PnoiseAt +2, the power threshold P is stoppedthfPerforming step-1 dB decrement operation;
finally, a carrier boundary set is obtainedWith carrierBorder, each element in the set constitutes a carrier, if any (fpl)k，fprk) Determining a carrier region PCarrier { (f)n，Pn)}，fplk≤fn≤fprk，(fn，Pn) Determining the carrier bandwidth BW as fpr by the element PSDk-fplkCenter frequency ofSignal powerWherein (f)i，Pi) E is the number of elements in the carrier region PCarrier set, i is 1, 2.
6. An apparatus for automatic detection and monitoring of communication signals based on progressive division of frequency domain thresholds, the apparatus comprising:
an acquisition unit for acquiring a data frame to be analyzed and dividing the frequency domain of the data frame into frequency points (f)n，Pn) Expressed by way of set PSD, wherein PSD { (f)n，Pn)}，fnFrequency of n +1 point, PnFor the measured power within the resolution bandwidth at that frequency, N is 0, 1, 2, …, N;
a processing unit for searching the power value P in the set PSDnDetermining the maximum value PmaxAnd a minimum value PminWhen (P)max-Pmin) And when the data frame to be analyzed contains the communication signal, determining that the data frame to be analyzed contains the communication signal when the data frame is more than or equal to 2 dB.
7. The apparatus of claim 6, wherein the processing unit, after determining that the data frame to be analyzed contains the communication signal, is further configured to:
initialization power threshold Pthn0＝(Pmax-Pmin)/4+PminBy the use of Pthn0Dividing frequency points in the PSD set to obtain a set NSI { (f)n，Pn) In which P isn≤Pth0，(fn，Pn) E is the PSD; performing histogram statistics on power values of frequency points in the NSI set according to a preset gear, determining the first gear with the most distributed points, and obtaining the first gear starting and stopping power with the most distributed points as PnstAnd Pnsp；
According to PnstAnd PnspSearching frequency points in the PSD set to obtain a set NSF { (f)n，Pn)}，Pnst≤Pn≤Pnsp，(fn，Pn) E, PSD, calculating a first central moment E (NSF) of the power value in the NSF set, namely considering the background noise power level P of the frame frequency domain datanoise＝E(NSF)；
Reinitializing the initial power threshold Pthf＝Pmax-2, with PthfFor frequency point (f) in the setx，Px) Dividing the image into sets LP { (fp)x，Ppx) And RP { (fp)x，Ppx) In which Ppx≤Pthf≤Ppx+1，(fpx，Ppx)∈PSD；Ppx≥Pthf≥Ppx+1，(fpx，Ppx) E is PSD, x is 0, 1, 2, …, N-1; reserving frequency points which are paired in the LPs and the RPs, wherein the frequency points are paired, namely, each point in the LPs corresponds to a point in one of the RPs, the two points respectively correspond to intersection points of the left roll-off edge and the right roll-off edge of the same carrier signal and threshold power, and M is searched to form a point pair;
taking M pairs of point pairs as elements in LP and RP respectively, corresponding M carriers to each pair of point pair, and using bandwidth boundary frequency fpl on the left side of the carrierxRight side bandwidth boundary frequency fprxThe combination defines one carrier, and the sets LP and RP will constitute the set of carrier boundaries { (fpl)x，fprx)}，Pplx≤Pthf≤Pplx+1，Pprx≥Pthf≥Pprx+1，(fplx，Pplx)∈LP，(fprx，Pprx) E.g., RP, the carriers formed by all elements in the carrier boundary set are confirmed one by one, element (fpl)k，fprk) Determined carrier region PCarrier { (f)n，Pn) Wherein, fplk≤fn≤fprk，(fn，Pn) Belongs to PSD, and searches power value P in set PCarriernFinding the maximum value PPCmaxAnd minimum value PPCminWhen (PPC)max-PPCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier valid carrier boundaries form a new carrier boundary set carrier border.
8. The apparatus of claim 7, wherein the processing unit is further configured to:
will power threshold PthfDecreasing by steps of-1 dB, i.e. Pthf_New＝Pthf-1, using a decreasing threshold value PthfA New pair of frequency points (f) in the set PSDn，Pn) The division is carried out again to obtain a carrier left side bandwidth boundary point set LN and a left side bandwidth boundary point set RN, and the sets LN and RN form a new carrier set NCarrirBorder { (fnl)x，fnrx)}；Pnlx≤Pthf_New≤Pnlx+1，Pnrx≥Pthf_New≥Pnrx+1，(fnlx，Pnlx)∈LN，(fnrx，Pnrx) E.g. RN, confirming the carriers formed by all elements in the set one by one, element (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) Belongs to PSD, and searches for power value P in set NCarriernTo obtain the maximum value PNCmaxAnd minimum value PNCminWhen (PNC)max-PNCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier-valid carrier boundaries form a new set of carrier boundaries NCarrierBorder.
9. The apparatus of claim 8, wherein the processing unit is further configured to:
for the sets Carrier Border and NCarrier Border, the carrier boundary in NCarrier Border (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) E PSD, if NCarrier does not contain a carrier region determined by any carrier boundary in Carrier Border, then this carrier boundary (fnl)k，fnrk) Adding into a New carrier boundary carrier _ New set; if the NCarrier includes a carrier region defined by a carrier boundary in the carrier border, it is judged that the NCarrier includes several carrier regions defined by carrier boundaries in the carrier border, and if it includes only one, the carrier border is determined (fnl)k，fnrk) Adding into a Carrier Border _ New set, if two or more, adding the contained multiple carrier boundaries into the Carrier Border _ New set;
the carrier border set CarrierBorder is updated so that CarrierBorder is CarrierBorder _ New.
10. The apparatus of claim 9, wherein the processing unit is further configured to:
when power threshold value Pthf＝PnoiseAt +2, the power threshold P is stoppedthfPerforming step-1 dB decrement operation;
finally, a carrier boundary set Carrier Border is obtained, each element in the set constitutes a carrier, and if the element (fpl)k，fprk) Determining a carrier region PCarrier { (f)n，Pn)}，fplk≤fn≤fprk，(fn，Pn) Determining the carrier bandwidth BW as fpr by the element PSDk-fplkCenter frequency of
Technical Field
The embodiment of the application relates to a communication signal automatic detection and monitoring method and device based on frequency domain threshold progressive segmentation.
Background
In a communication system, an operator or a regulatory body needs to monitor and analyze communication signals so as to find signal interference or other abnormalities in time, which is called carrier monitoring. The carrier monitoring mainly comprises two key processes of signal detection and signal analysis, wherein the signal detection completes blind analysis on signals in a monitoring bandwidth, and extracts all carrier signals; the signal analysis is a parameter estimation for each carrier signal based on the result of the signal detection.
Currently, there are two main types of carrier monitoring methods: time domain based analysis and frequency domain based analysis. In the former, a communication signal needs to be changed into an intermediate frequency, the intermediate frequency signal is sampled by a high-speed data acquisition technology, the data obtained by sampling is subjected to deep analysis by a software radio method, signals in a frequency band range are detected, and key parameters of the signals are estimated, such as a Kratos's Monics system, an Atos's SkyMon system and the like. The method has the advantages that the system is simple in structure, low in cost and high in analysis speed, and the defects that the method is based on frequency domain analysis, lacks time domain data, and is difficult to carry out deep analysis of signals, such as modulation mode identification, Carrier Under Carrier (CUC) analysis and the like.
With the gradual development of communication technologies, especially the application of technologies such as satellite communication spot beams and High Throughput (HTS), the carrier monitoring system based on time domain analysis is expensive in cost, and is difficult to meet the requirement for full monitoring of multi-spot beams.
The signal blind detection technology in the carrier monitoring system is almost developed based on time domain analysis, and signal detection based on frequency domain mostly needs to rely on artificial judgment to find signal abnormality in engineering.
Disclosure of Invention
In order to solve the above technical problems, embodiments of the present application desirably provide an automatic detection and monitoring method and apparatus for a communication signal based on frequency domain threshold progressive segmentation.
The technical scheme of the invention is realized as follows:
the embodiment of the application provides a communication signal automatic detection and monitoring method based on frequency domain threshold progressive segmentation, which comprises the following steps:
acquiring a data frame to be analyzed, and dividing the frequency domain of the data frame into frequency points (f)n，Pn) Expressed by way of set PSD, wherein PSD { (f)n，Pn)}，fnFrequency of n +1 point, PnFor the measured power within the resolution bandwidth at that frequency, N is 0, 1, 2, …, N;
searching for a power value P within a set PSDnDetermining the maximum value PmaxAnd a minimum value PminWhen (P)max-Pmin) And when the data frame to be analyzed contains the communication signal, determining that the data frame to be analyzed contains the communication signal when the data frame is more than or equal to 2 dB.
In the above technical solution, after determining that the data frame to be analyzed includes the communication signal, the method further includes:
initialization power threshold Pthn0＝(Pmax-Pmin)/4+PminBy the use of Pthn0Dividing frequency points in the PSD set to obtain a set NSI { (f)n，Pn) In which P isn≤Pth0，(fn，Pn) E is the PSD; carrying out histogram statistics on the power values of the frequency points in the NSI set according to a preset gear to determine the most distributed pointsThe starting power and the stopping power of the first gear with the maximum number of distributed points are respectively PnstAnd Pnsp；
According to PnstAnd PnspSearching frequency points in the PSD set to obtain a set NSF { (f)n，Pn)}，Pnst≤Pn≤Pnsp，(fn，Pn) E, PSD, calculating a first central moment E (NSF) of the power value in the NSF set, namely considering the background noise power level P of the frame frequency domain datanoise＝E(NSF)；
Reinitializing the initial power threshold Pthf＝Pmax-2, with PthfFor frequency point (f) in the setx，Px) Dividing the image into sets LP { (fp)x，Ppx) And Rp { (fp)x，Ppx) In which Ppx≤Pthf≤Ppx+1，(fpx，Ppx)∈PSD；Ppx≥Pthf≥Ppx+1，(fpx，Ppx) E is PSD, x is 0, 1, 2, …, N-1; reserving frequency points which are paired in the LPs and the RPs, wherein the frequency points are paired, namely, each point in the LPs corresponds to a point in one of the RPs, the two points respectively correspond to intersection points of the left roll-off edge and the right roll-off edge of the same carrier signal and threshold power, and M is searched to form a point pair;
taking M pairs of point pairs as elements in LP and RP respectively, corresponding M carriers to each pair of point pair, and using bandwidth boundary frequency fpl on the left side of the carrierxRight side bandwidth boundary frequency fprxThe combination defines one carrier, and the sets LP and RP will constitute the set of carrier boundaries { (fpl)x，fprx)}，Pplx≤Pthf≤Pplx+1，Pprx≥Pthf≥Pprx+1，(fplx，Pplx)∈LP，(fprx，Pprx) E.g., RP, the carriers formed by all elements in the carrier boundary set are confirmed one by one, element (fpl)k，fprk) Determined carrier region PCarrier { (f)n，Pn) Wherein, fplk≤fn≤fprk，(fn，Pn) E.g. PSD, searchSet of power values P within PCarriernFinding the maximum value PPCmaxAnd minimum value PPCminWhen (PPC)max-PPCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier valid carrier boundaries form a new carrier boundary set carrier border.
In the above technical solution, the method further includes:
will power threshold PthfDecreasing by steps of-1 dB, i.e. Pthf_New＝Pthf-1, using a decreasing threshold value PthfA New pair of frequency points (f) in the set PSDn，Pn) The division is carried out again to obtain a carrier left side bandwidth boundary point set LN and a left side bandwidth boundary point set RN, and the sets LN and RN form a new carrier set NCarrirBorder { (fnl)x，fnrx)}；Pnlx≤Pthf_New≤Pnlx+1，Pnrx≥Pthf_New≥Pnrx+1，(fnlx，Pnlx)∈LN，(fnrx，Pnrx) E.g. RN, confirming the carriers formed by all elements in the set one by one, element (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) Belongs to PSD, and searches for power value P in set NCarriernTo obtain the maximum value PNCmaxAnd minimum value PNCminWhen (PNC)max-PNCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier-valid carrier boundaries form a new set of carrier boundaries NCarrierBorder.
In the above technical solution, the method further includes:
for the sets Carrier Border and NCarrier Border, the carrier boundary in NCarrier Border (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) E.g., PSD, if NCarrier does not contain the code of Carrier BordThe carrier region determined by any carrier boundary in the er, the carrier boundary (fnl)k，fnrk) Adding into a New carrier boundary carrier _ New set; if the NCarrier includes a carrier region defined by a carrier boundary in the carrier border, it is judged that the NCarrier includes several carrier regions defined by carrier boundaries in the carrier border, and if it includes only one, the carrier border is determined (fnl)k，fnrk) Adding into a Carrier Border _ New set, if two or more, adding the contained multiple carrier boundaries into the Carrier Border _ New set;
the carrier border set CarrierBorder is updated so that CarrierBorder is CarrierBorder _ New.
In the above technical solution, the method further includes:
when power threshold value Pthf＝PnoiseAt +2, the power threshold P is stoppedthfPerforming step-1 dB decrement operation;
finally, a carrier boundary set Carrier Border is obtained, each element in the set constitutes a carrier, and if the element (fpl)k，fprk) Determining a carrier region PCarrier { (f)n，Pn)}，fplk≤fn≤fprk，(fn，Pn) Determining the carrier bandwidth BW as fpr by the element PSDk-fplkCenter frequency of
The invention also provides a communication signal automatic detection and monitoring device based on frequency domain threshold progressive segmentation, which comprises:
an acquisition unit for acquiring a data frame to be analyzed and dividing the frequency domain of the data frame into frequency points (f)n，Pn) Expressed by the way of integrating PSD, wherein PSD is equal to{(fn，Pn)}，fnFrequency of n +1 point, PnFor the measured power within the resolution bandwidth at that frequency, N is 0, 1, 2, …, N;
a processing unit for searching the power value P in the set PSDnDetermining the maximum value PmaxAnd a minimum value PminWhen (P)max-Pmin) And when the data frame to be analyzed contains the communication signal, determining that the data frame to be analyzed contains the communication signal when the data frame is more than or equal to 2 dB.
In the above technical solution, after the processing unit determines that the data frame to be analyzed includes the communication signal, the processing unit is further configured to:
initialization power threshold Pthn0＝(Pmax-Pmin)/4+PminBy the use of Pthn0Dividing frequency points in the PSD set to obtain a set NSI { (f)n，Pn) In which P isn≤Pth0，(fn，Pn) E is the PSD; performing histogram statistics on power values of frequency points in the NSI set according to a preset gear, determining the first gear with the most distributed points, and obtaining the first gear starting and stopping power with the most distributed points as PnstAnd Pnsp；
According to PnstAnd PnspSearching frequency points in the PSD set to obtain a set NSF { (f)n，Pn)}，Pnst≤Pn≤Pnsp，(fn，Pn) E, PSD, calculating a first central moment E (NSF) of the power value in the NSF set, namely considering the background noise power level P of the frame frequency domain datanoise＝E(NSF)；
Reinitializing the initial power threshold Pthf＝Pmax-2, with PthfFor frequency point (f) in the setx，Px) Dividing the image into sets LP { (fp)x，Ppx) And RP { (fp)x，Ppx) In which Ppx≤Pthf≤Ppx+1，(fpx，Ppx)∈PSD；Ppx≥Pthf≥Ppx+1，(fpx，Ppx) E is PSD, x is 0, 1, 2, …, N-1; reserving the frequency points in LP and RP which appear in pairsThe point in each LP corresponds to a point in one RP, the two points respectively correspond to the intersection points of the left roll-off edge and the right roll-off edge of the same carrier signal and the threshold power, and M is searched to form a point;
taking M pairs of point pairs as elements in LP and RP respectively, corresponding M carriers to each pair of point pair, and using bandwidth boundary frequency fpl on the left side of the carrierxRight side bandwidth boundary frequency fprxThe combination defines one carrier, and the sets LP and RP will constitute the set of carrier boundaries { (fpl)x，fprx)}，Pplx≤Pthf≤Pplx+1，Pprx≥Pthf≥Pprx+1，(fplx，Pplx)∈LP，(fprx，Pprx) E.g., RP, the carriers formed by all elements in the carrier boundary set are confirmed one by one, element (fpl)k，fprk) Determined carrier region PCarrier { (f)n，Pn) Wherein, fplk≤fn≤fprk，(fn，Pn) Belongs to PSD, and searches power value P in set PCarriernFinding the maximum value PPCmaxAnd minimum value PPCminWhen (PPC)max-PPCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier valid carrier boundaries form a new carrier boundary set carrier border.
In the foregoing technical solution, the processing unit is further configured to:
will power threshold PthfDecreasing by steps of-1 dB, i.e. Pthf_New＝Pthf-1, using a decreasing threshold value PthfA New pair of frequency points (f) in the set PSDn，Pn) The division is carried out again to obtain a carrier left side bandwidth boundary point set LN and a left side bandwidth boundary point set RN, and the sets LN and RN form a new carrier set NCarrirBorder { (fnl)x，fnrx}；Pnlx≤Pthf_New≤Pnlx+1，Pnrx≥Pthf_New≥Pnrx+1，(fnlx，Pnlx)∈LN，(fnrx，Pnrx) E.g. RN, for all in the setThe carriers comprising the elements are identified one by one, the elements (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) Belongs to PSD, and searches for power value P in set NCarriernTo obtain the maximum value PNCmaxAnd minimum value PNCminWhen (PNC)max-PNCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier-valid carrier boundaries form a new set of carrier boundaries NCarrierBorder.
In the foregoing technical solution, the processing unit is further configured to:
for the sets Carrier Border and NCarrier Border, the carrier boundary in NCarrier Border (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) E PSD, if NCarrier does not contain a carrier region determined by any carrier boundary in Carrier Border, then this carrier boundary (fnl)k，fnrk) Adding into a New carrier boundary carrier _ New set; if the NCarrier includes a carrier region defined by a carrier boundary in the carrier border, it is judged that the NCarrier includes several carrier regions defined by carrier boundaries in the carrier border, and if it includes only one, the carrier border is determined (fnl)k，fnrk) Adding into a Carrier Border _ New set, if two or more than two are contained, adding a plurality of contained carrier boundaries into the Carrier Border _ New set;
the carrier border set CarrierBorder is updated so that CarrierBorder is CarrierBorder _ New.
In the foregoing technical solution, the processing unit is further configured to:
when power threshold value Pthf＝PnoiseAt +2, the power threshold P is stoppedthfPerforming step-1 dB decrement operation;
finally, a carrier boundary set Carrier Border is obtained, and each element in the set formsOne carrier wave, several elements (fpl)k，fprk) Determining a carrier region PCarrier { (f)n，Pn)}，fplk≤fn≤fprk，(fn，Pn) Determining the carrier bandwidth BW as fpr by the element PSDk-fplkCenter frequency of
According to the technical scheme of the embodiment of the application, the communication signal automatic detection and analysis based on the frequency domain data are realized, the method and the device have the characteristics of automation, low cost, easiness in implementation and good real-time performance, are suitable for a scene of deployment of multiple monitoring points, and expand the implementation mode of a communication signal carrier monitoring system.
Drawings
Fig. 1 is a schematic flowchart of an automatic detection and monitoring method for a communication signal based on frequency domain threshold progressive segmentation according to an embodiment of the present application;
fig. 2 is a schematic flow chart illustrating an implementation of a communication signal automatic detection and monitoring method based on frequency domain threshold progressive segmentation according to an embodiment of the present application;
fig. 3 is a schematic structural diagram of an automatic communication signal detection and monitoring apparatus based on frequency domain threshold progressive division according to an embodiment of the present application.
Detailed Description
The embodiments described in the present invention can be combined without conflict.
Fig. 1 is a schematic flowchart of an automatic detection and monitoring method for a communication signal based on frequency domain threshold progressive segmentation according to an embodiment of the present application, and as shown in fig. 1, the automatic detection and monitoring method for a communication signal based on frequency domain threshold progressive segmentation according to the present example includes the following steps:
In the above technical solution, after determining that the data frame to be analyzed includes the communication signal, the method further includes:
initialization power threshold Pthn0＝(Pmax-Pmin)/4+PminBy the use of Pthn0Dividing frequency points in the PSD set to obtain a set NSI { (f)n，Pn) In which P isn≤Pth0，(fn，Pn) E is the PSD; performing histogram statistics on power values of frequency points in the NSI set according to a preset gear, determining the first gear with the most distributed points, and obtaining the first gear starting and stopping power with the most distributed points as PnstAnd Pnsp(ii) a Here, the preset gear may be 10 th gear, 15 th gear, 20 th gear, 5 th gear, etc., and may be set according to process requirements.
According to PnstAnd PnspSearching frequency points in the PSD set to obtain a set NSF { (f)n，Pn)}，Pnst≤Pn≤Pnsp，(fn，Pn) E, PSD, calculating a first central moment E (NSF) of the power value in the NSF set, namely considering the background noise power level P of the frame frequency domain datanoise＝E(NSF)；
Reinitializing the initial power threshold Pthf＝Pmax-2, with PthfFor frequency point (f) in the setx，Px) Dividing the image into sets LP { (fp)x，Ppx) And RP { (fp)x，Ppx) In which Ppx≤Pthf≤Ppx+1，(fpx，Ppx)∈PSD；Ppx≥Pthf≥Ppx+1，(fpx，Ppx) E is PSD, x is 0, 1, 2, …, N-1; reserving frequency points which are paired in the LPs and the RPs, wherein the frequency points are paired, namely, each point in the LPs corresponds to a point in one of the RPs, the two points respectively correspond to intersection points of the left roll-off edge and the right roll-off edge of the same carrier signal and threshold power, and M is searched to form a point pair;
taking M pairs of point pairs as elements in LP and RP respectively, corresponding M carriers to each pair of point pair, and using bandwidth boundary frequency fpl on the left side of the carrierxRight side bandwidth boundary frequency fprxThe combination defines one carrier, and the sets LP and RP will constitute the set of carrier boundaries { (fpl)x，fprx)}，Pplx≤Pthf≤Pplx+1，Pprx≥Pthf≥Pprx+1，(fplx，Pplx)∈LP，(fprx，Pprx) E.g., RP, the carriers formed by all elements in the carrier boundary set are confirmed one by one, element (fpl)k，fprk) Determined carrier region PCarrier { (f)n，Pn) Wherein, fplk≤fn≤fprk，(fn，Pn) Belongs to PSD, and searches power value P in set PCarriernFinding the maximum value PPCmaxAnd minimum value PPCminWhen (PPC)max-PPCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier valid carrier boundaries form a new carrier boundary set carrier border.
In specific implementation, the method can further comprise the following processing steps:
will power threshold PthfDecreasing by steps of-1 dB, i.e. Pthf_New＝Pthf-1, using a decreasing threshold value PthfA New pair of frequency points (f) in the set PSDn，Pn) The division is carried out again to obtain a carrier left side bandwidth boundary point set LN and a left side bandwidth boundary point set RN, and the sets LN and RN form a new carrier set NCarrirBorder { (fnl)x，fnrx)}；Pnlx≤Pthf_New≤Pnlx+1，Pnrx≥Pthf_New≥Pnrx+1，(fnlx，Pnlx)∈LN，(fnrx，Pnrx) E.g. RN, confirming the carriers formed by all elements in the set one by one, element (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) Belongs to PSD, and searches for power value P in set NCarriernTo obtain the maximum value PNCmaxAnd minimum value PNCminWhen (PNC)max-PNCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier-valid carrier boundaries form a new set of carrier boundaries NCarrierBorder.
In specific implementation, the method can further comprise the following processing steps:
for the sets Carrier Border and NCarrier Border, the carrier boundary in NCarrier Border (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) E PSD, if NCarrier does not contain a carrier region determined by any carrier boundary in Carrier Border, then this carrier boundary (fnl)k，fnrk) Adding into a New carrier boundary carrier _ New set; if the NCarrier includes a carrier region defined by a carrier boundary in the carrier border, it is judged that the NCarrier includes several carrier regions defined by carrier boundaries in the carrier border, and if it includes only one, the carrier border is determined (fnl)k，fnrk) Adding into a Carrier Border _ New set, if two or more, adding the contained multiple carrier boundaries into the Carrier Border _ New set;
the carrier border set CarrierBorder is updated so that CarrierBorder is CarrierBorder _ New.
In specific implementation, the method can further comprise the following processing steps:
when power threshold value Pthf＝PnoiseAt +2, stopStopping on power threshold PthfPerforming step-1 dB decrement operation;
finally, a carrier boundary set Carrier Border is obtained, each element in the set constitutes a carrier, and if the element (fpl)k，fprk) Determining a carrier region PCarrier { (f)n，Pn)}，fplk≤fn≤fprk，(fn，Pn) Determining the carrier bandwidth BW as fpr by the element PSDk-fplkCenter frequency ofSignal power
According to the technical scheme of the embodiment of the application, the communication signal automatic detection and analysis based on the frequency domain data are realized, the method and the device have the characteristics of automation, low cost, easiness in implementation and good real-time performance, are suitable for a scene of deployment of multiple monitoring points, and expand the implementation mode of a communication signal carrier monitoring system.
The embodiment of the present application takes the russian communication satellite project as an example, and further illustrates the essence of the technical solution of the embodiment of the present invention.
The white Russian communication satellite comprises white Russian regional beams and African regional beams, the ground carrier monitoring system consists of two parts, one part is deployed in a white Russian ground station, and the system is built by a Zodiac product based on time domain analysis; the other part is deployed in an African Nigeria ground station, and the system is composed of LPT-3000R based on frequency domain analysis. LPT-3000R is a low cost simple spectral analysis instrument offered by LPT corporation, and each frame of power spectral data consists of 401 points. The monitoring client is connected with LPT-3000R through VPN network in the white Russian ground station.
Fig. 2 is a schematic flow chart of an implementation process of the communication signal automatic detection and monitoring method based on frequency domain threshold progressive segmentation according to the embodiment of the present application, and as shown in fig. 2, the implementation of the communication signal automatic detection and monitoring method based on frequency domain threshold progressive segmentation according to the embodiment of the present application includes the following processing steps:
1. setting the RBW of LPT-3000R to be 100kHz, the VBW to be 1kHz, and setting the SPAN bandwidth to be 40MHz according to the frequency range of the satellite transponder;
2. acquiring a frame of power spectrum data, and searching for a maximum value P of powermaxAnd a minimum value Pmin，(Pmax-Pmin) When the data is more than or equal to 2dB, the frame data is considered to contain communication signals;
3. defining an initialization power threshold as Pthn0＝(Pmax-Pmin)/4+PminDividing 401 (frequency, power) points by using the threshold, performing histogram statistics of 10-grade or 15-grade for the area below the threshold, searching the first-grade histogram with the most distributed points, performing statistics on the corresponding power distributed in the first-grade histogram, calculating the first-order central moment of the first-grade histogram, and using the first-order central moment as the background noise power level P of the frame power spectrum datanoise；
4. Redefining the Power initial threshold Pthf＝Pmax-2, using the threshold to segment 401 (frequency, power) points, eliminating the segment points that do not appear in pairs, using the segment points that appear in pairs as the left and right boundary component carriers, searching for the maximum power value PPC of the points in the carrier regionmaxAnd minimum value PPCmin，(PPCmax-PPCmin) And when the number of the carrier elements is more than or equal to 2dB, considering the carrier formed by the carrier boundary to be effective, otherwise, discarding the carrier boundary elements. The effective paired segmentation points form a carrier boundary set Carrier Border;
5. for power threshold PthfDecreasing by steps of-1 dB, i.e. Pthf_New＝Pthf1, dividing 401 (frequency, power) points by using a new threshold, and screening according to the method in the step 4 to obtain a new carrier boundary set NCarrierBorder;
6. the sets carrierBorder and NCarrierBorder in step 4 and step 5, previously described, are compared for carrier boundaries in NCarrierBorder (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)|fnlk≤fn≤fnrk，(fn，Pn) E.g., PSD, N is 0, 1, 2,.., N }, which carrier boundary if NCarrier does not contain a carrier region determined by any carrier boundary in carrier border (fnl)k，fnrk) Adding into a New carrier boundary carrier _ New set; if the NCarrier includes a carrier region defined by a carrier boundary in the carrier border, it is determined that the NCarrier includes several carrier regions defined by carrier boundaries in the carrier border, and if only one carrier region is included, the carrier boundary is determined (fnl)k，fnrk) Adding into a Carrier Border _ New set, if a plurality (two or more) are contained, adding the contained plurality of carrier boundaries into the Carrier Border _ New set;
7. updating a carrier boundary set carrier border in the step 4, wherein carrier border is CarrierBorder _ New;
8. repeating the processes in the step 5, the step 6 and the step 7, and when the power threshold value P is in the step 5thf＝PnoiseAt +2, the power threshold P is stoppedthfDecreasing progressively;
9. the result is a Carrier Border set, where each element in the set constitutes a carrier, e.g., element (fpl)k，fprk) Determining a carrier region PCarrier { (f)n，Pn)|fplk≤fn≤fprk，(fn，Pn) E.g., PSD, N is 0, 1, 2,.., N }, and the carrier bandwidth BW is fprk-fplkCenter frequency of
Fig. 3 is a schematic structural diagram of a communication signal automatic detection and monitoring device based on frequency domain threshold progressive division according to an embodiment of the present application, and as shown in fig. 3, the communication signal automatic detection and monitoring device based on frequency domain threshold progressive division according to the embodiment of the present application includes:
an obtaining unit 30 for obtaining a data frame to be analyzed and calculating a frequency domain of the data frame by a frequency point (f)n，Pn) Expressed by way of set PSD, wherein PSD { (f)n，Pn)}，fnFrequency of n +1 point, PnFor the measured power within the resolution bandwidth at that frequency, N is 0, 1, 2, …, N;
processing unit 31, searching for power value P in set PSDnDetermining the maximum value PmaxAnd a minimum value PminWhen (P)max-Pmin) And when the data frame to be analyzed contains the communication signal, determining that the data frame to be analyzed contains the communication signal when the data frame is more than or equal to 2 dB.
In the above technical solution, after the processing unit 31 determines that the data frame to be analyzed includes the communication signal, it is further configured to:
initialization power threshold Pthn0＝(Pmax-Pmin)/4+PminBy the use of Pthn0Dividing frequency points in the PSD set to obtain a set NSI { (f)n，Pn) In which P isn≤Pth0，(fn，Pn) E is the PSD; performing histogram statistics on power values of frequency points in the NSI set according to a preset gear, determining the first gear with the most distributed points, and obtaining the first gear starting and stopping power with the most distributed points as PnstAnd Pnsp；
According to PnstAnd PnspSearching frequency points in the PSD set to obtain a set NSF { (f)n，Pn)}，Pnst≤Pn≤Pnsp，(fn，Pn) E, PSD, calculating a first central moment E (NSF) of the power value in the NSF set, namely considering the background noise power level P of the frame frequency domain datanoise＝E(NSF)；
Reinitializing the initial power threshold Pthf＝Pmax-2, with PthfFor frequency point (f) in the setx，Px) Dividing the image into sets LP { (fp)x，Ppx) And RP { (fp)x，Ppx) In which Ppx≤Pthf≤Ppx+1，(fpx，Ppx)∈PSD；Ppx≥Pthf≥Ppx+1，(fpx，Ppx) E is PSD, x is 0, 1, 2, …, N-1; reserving frequency points which are paired in the LPs and the RPs, wherein the frequency points are paired, namely, each point in the LPs corresponds to a point in one of the RPs, the two points respectively correspond to intersection points of the left roll-off edge and the right roll-off edge of the same carrier signal and threshold power, and M is searched to form a point pair;
taking M pairs of point pairs as elements in LP and RP respectively, corresponding M carriers to each pair of point pair, and using bandwidth boundary frequency fpl on the left side of the carrierxRight side bandwidth boundary frequency fprxThe combination defines one carrier, and the sets LP and RP will constitute the set of carrier boundaries { (fpl)x，fprx)}，Pplx≤Pthf≤Pplx+1，Pprx≥Pthf≥Pprx+1，(fplx，Pplx)∈LP，(fprx，Pprx) E.g., RP, the carriers formed by all elements in the carrier boundary set are confirmed one by one, element (fpl)k，fprk) Determined carrier region PCarrier { (f)n，Pn) Wherein, fplk≤fn≤fprk，(fn，Pn) Belongs to PSD, and searches power value P in set PCarriernFinding the maximum value PPCmaxAnd minimum value PPCminWhen (PPC)max-PPCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier valid carrier boundaries form a new carrier boundary set carrier border.
In the above technical solution, the processing unit 31 is further configured to:
will power threshold PthfDecreasing by steps of-1 dB, i.e. Pthf_New＝Pthf-1, using a decreasing threshold value PthfA New pair of frequency points (f) in the set PSDn，Pn) Re-dividing to obtain a carrier left side bandwidth boundary point set LN and a left side bandwidthSet of boundary points RN, sets LN and RN will constitute a new set of carriers NCarrierBorder { (fnl)x，fnrx)}；Pnlx≤Pthf_New≤Pnlx+1，Pnrx≥Pthf_New≥Pnrx+1，(fnlx，Pnlx)∈LN，(fnrx，Pnrx) E.g. RN, confirming the carriers formed by all elements in the set one by one, element (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) Belongs to PSD, and searches for power value P in set NCarriernTo obtain the maximum value PNCmaxAnd minimum value PNCminWhen (PNC)max-PNCmin) When the data rate is more than or equal to 2dB, determining that the carrier formed by the carrier boundary is effective;
the carrier-valid carrier boundaries form a new set of carrier boundaries NCarrierBorder.
In the above technical solution, the processing unit 31 is further configured to:
for the sets Carrier Border and NCarrier Border, the carrier boundary in NCarrier Border (fnl)k，fnrk) Determined carrier region NCarrier { (f)n，Pn)}，fnlk≤fn≤fnrk，(fn，Pn) E PSD, if NCarrier does not contain a carrier region determined by any carrier boundary in Carrier Border, then this carrier boundary (fnl)k，fnrk) Adding into a New carrier boundary carrier _ New set; if the NCarrier includes a carrier region defined by a carrier boundary in the carrier border, it is judged that the NCarrier includes several carrier regions defined by carrier boundaries in the carrier border, and if it includes only one, the carrier border is determined (fnl)k，fnrk) Adding into a Carrier Border _ New set, if two or more, adding the contained multiple carrier boundaries into the Carrier Border _ New set;
the carrier border set CarrierBorder is updated so that CarrierBorder is CarrierBorder _ New.
In the above technical solution, the processing unit 31 is further configured to:
when power threshold value Pthf＝PnoiseAt +2, the power threshold P is stoppedthfPerforming step-1 dB decrement operation;
finally, a carrier boundary set Carrier Border is obtained, each element in the set constitutes a carrier, and if the element (fpl)k，fprk) Determining a carrier region PCarrier { (f)n，Pn)}，fplk≤fn≤fprk，(fn，Pn) Determining the carrier bandwidth BW as fpr by the element PSDk-fplkCenter frequency of
According to the technical scheme of the embodiment of the application, the communication signal automatic detection and analysis based on the frequency domain data are realized, the method and the device have the characteristics of automation, low cost, easiness in implementation and good real-time performance, are suitable for a scene of deployment of multiple monitoring points, and expand the implementation mode of a communication signal carrier monitoring system.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of a hardware embodiment, a software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only a preferred embodiment of the present application, and is not intended to limit the scope of the present application.
- 上一篇：一种医用注射器针头装配设备
- 下一篇：一种调制方式识别方法及装置