CN110535543A - A kind of radio-frequency spectrum signal detecting threshold calculation method - Google Patents
A kind of radio-frequency spectrum signal detecting threshold calculation method Download PDFInfo
- Publication number
- CN110535543A CN110535543A CN201910392957.8A CN201910392957A CN110535543A CN 110535543 A CN110535543 A CN 110535543A CN 201910392957 A CN201910392957 A CN 201910392957A CN 110535543 A CN110535543 A CN 110535543A
- Authority
- CN
- China
- Prior art keywords
- frequency spectrum
- frequency
- subsegment
- level
- data block
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The present invention provides a kind of radio-frequency spectrum signal detecting threshold calculation methods.The method of the present invention removes valley and exponential smoothing operation to the frequency spectrum data execution that radio-frequency transmitter exports；Frequency spectrum data is divided into wavelength coverage by frequency separation, then wavelength coverage is divided into subsegment, maximum spectral line level in each subsegment is calculated, takes the minimum value of subsegment maximum spectral line level as wavelength coverage with reference to ambient level；Finally, being minimized between frequency spectrum segment data and section are with reference to ambient level as signal detecting threshold.The calculated signal detecting threshold of the method for the present invention can adapt to various bandwidth signals and Testing of Feeble Signals, and algorithm committed memory is few, and computing cost is low, adapt to various fixations, the requirement that airborne and onboard wireless pyroelectric monitor system is to signal detecting threshold calculated performance.
Description
Technical field
The present invention relates to monitoring radio-frequency spectrum technical field, specially a kind of radio-frequency spectrum signal detecting threshold is calculated
Method.
Background technique
The groundwork process of monitoring radio-frequency spectrum system is to execute spectrum scan in monitoring frequency range, detects in frequency spectrum
Existing signal, and statistics and analysis is carried out to detecting signal.
The groundwork process of monitoring radio-frequency spectrum system is to execute spectrum scan in monitoring frequency range, detects in frequency spectrum
Existing signal, and statistics and analysis is carried out to detecting signal.Detection threshold is the foundation of signal detection, when some frequency separation frequency
When spectrum level is higher than detection threshold, it is meant that signal occurs in this frequency separation, and frequency spectrum monitoring system enrolls its frequency, level
And the parameters such as bandwidth, and further analysis and classification are executed to it.
Requirement of the frequency spectrum monitoring system to signal detecting threshold algorithm is can to adapt to narrowband and Broadband Detection simultaneously,
It can adapt to Testing of Feeble Signals, and have dynamic updating ability to adapt to the variation of electromagnetic environment.
Level measurement precision, dynamic range and the polydispersity that mutually reconciles of radio-frequency transmitter hardware decision frequency spectrum monitoring system
Energy.Executing the Fast Fourier Transform (FFT) of time domain/frequency-domain transform is mature technology.Signal detecting threshold superiority and inferiority determines spectrum monitoring system
System intercepts and captures the ability of various bandwidth signals, intercepts and captures the ability of weak signal, the ability of electromagnetic environment variation is adapted to, to finally determine
Frequency spectrum monitoring system performance.
In radio monitoring technical field, signal detecting threshold can be classified as level threshold, ambient threshold and automatic threshold:
1, level threshold: specifying a fixed level value as signal detecting threshold by operator, is only applicable in a small number of special
Surely task is monitored.
2, ambient threshold: the maximum of period is used to keep frequency spectrum as signal detecting threshold.It is mainly used for this period
Emerging signal detection afterwards.
3, automatic threshold: the signal detecting threshold calculated automatically according to frequency spectrum data detects foundation as signal.Automatically-controlled door
Limit is a kind of widest signal detecting threshold of the scope of application.
Mainly pass through " rolling average " or " clipping section mean value " two methods in existing frequency spectrum monitoring system to calculate automatically
Signal detecting threshold, " rolling average " algorithm calculate the mean value in the adjacent section of each frequency spectrum point, and as signal detection door
Limit；Thresholding calculated in this way is only used for narrow band signal detection, and computing cost is larger, it is difficult to which dynamic updates." clipping Duan Jun
Value " algorithm calculates section mean value first, using section mean value to section mean value is calculated again after r level clipping, by the Duan Jun after clipping
Value is used as signal detecting threshold；When signal is more or when signal is wider in frequency range, threshold level calculated in this way and practical frequency
It is larger to compose background difference, is easy to omit weak signal in signal detection.
Summary of the invention
It is an object of the invention to provide a kind of radio-frequency spectrum signal detecting threshold calculation methods, to solve existing skill
The problem of wide band energy detection, weak energy measuring and detection threshold dynamic update cannot be had both in art simultaneously.
To achieve the above object, the invention provides the following technical scheme: a kind of radio-frequency spectrum signal detecting threshold calculates
Method, steps are as follows:
S1, parameter setting: setting median filter window width constant group WMED, section bandwidth constant group WSEG, subsegment bandwidth
Constant group WBKG, the value and frequency dependence of member node in this three array；Indexing smoothing filter factor beta.
S2, data block reception and data processing；
S3, data block segmentation: according to section bandwidth array WSEGWith the frequency separation f of frequency spectrum data blockB1~fBK, by spectrum number
N number of frequency spectrum data section SEGn, 1≤n≤N are divided into according to block；There is P data point in frequency spectrum data section SEGn, each data point
Level is LS1~LSP；
S4, subsegment segmentation: in each frequency spectrum data section in SEGn, according to subsegment bandwidth array WBKGIn the wavelength coverage of sum
Frequency of heart fSEG, frequency spectrum data section is divided into M subsegment；
S5, subsegment maximum level are sought: the maximum spectral line level sought in M subsegment obtains subsegment maximum level Lmax1~
LmaxM；
S6, wavelength coverage are sought with reference to ambient level: seeking maximum spectral line level L in M subsegmentmax1~LmaxMMinimum value,
It obtains wavelength coverage and refers to ambient level LREF；
S7, wavelength coverage signal detecting threshold are sought: in the data point level value L of frequency spectrum data section SEGnS1~LSPAnd frequency spectrum
Section refers to ambient level LREFIn be minimized, obtain wavelength coverage signal detecting threshold LT1~LTPAs output；
It the step of between S8, repetition S4 to S7, is finished until all frequency spectrum data sections are all processed.
Preferably, data block reception and data processing include: in the step S2
S21, data block reception: the wide-band frequency spectrum data block B of radio-frequency transmitter output is receivedR, frequency spectrum data points
For K, the frequency of data point is fB1~fBK, the level of data point is LB1~LBK；
S22, valley point filtering: median filtering is executed to the valley point in frequency spectrum data block BR, acquisition eliminates valley point
Frequency spectrum data block BU, median filter length of window depends on median filtering constant group WMED and the frequency spectrum at spectrum number strong point divides
Resolution；
S23, exponential smoothing: to frequency spectrum data block BUExponential smoothing is executed, filter coefficient β obtains smoothed out frequency
Modal data block BE。
Compared with prior art, the beneficial effects of the present invention are: the embodiment of the invention provides a kind of radio-frequency spectrum letters
Number detection threshold calculation method, valley point filtering is for eliminating in frequency spectrum data since spectral line caused by frequency content missing collapses
It falls into；Exponential smoothing is for eliminating the fluctuating of spectral line caused by noise；Exponential smoothing operation only stores a data, committed memory
Few, computing cost is low；Data block segmentation, which plays, to be avoided calculating ambient level in frequency spectrum section and rises and falls excessive effect；Subsegment point
Cut, subsegment maximum level is sought and wavelength coverage with reference to ambient level finding process be the operation unrelated with signal bandwidth, frequency
There are can obtain reasonably to examine ambient level when narrowband or broadband signal in modal data；Wavelength coverage signal detecting threshold was sought
Be minimized between Cheng Kao ambient level and section spectrum number strong point, guarantee the signal detecting threshold of output close to ambient level,
To improve Testing of Feeble Signals ability；Inventive algorithm adapts to various fixations, airborne and onboard wireless pyroelectric monitor system to letter
The requirement of number detection threshold calculated performance.
Detailed description of the invention
Fig. 1 is the radio-frequency spectrum signal detecting threshold calculation method flow diagram in the embodiment of the present invention.
Fig. 2 is the method for the present invention schematic diagram.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The present invention provides a kind of technical solution referring to FIG. 1-2: a kind of radio-frequency spectrum signal detecting threshold calculating side
Method, steps are as follows:
S1, parameter setting: setting median filter window width constant group WMEDSection bandwidth constant group WSEG, subsegment bandwidth it is normal
Array WBKG, the value and frequency dependence of member node in this three array；Indexing smoothing filter factor beta.
Wherein；
Radio-frequency transmitter working frequency section is 30MHz~6GHZ, intermediate-frequency bandwidth BWIF=40MHZ in this example；Radio frequency
Receiver API exports not averaged Fast Fourier Transform frequency spectrum data block, and the spectral resolution that monitoring task uses is
2.4KHz；According to radio-frequency transmitter operating condition, it is 0.0625 that exponential smoothing filter factor beta, which is arranged, and median filter is arranged
Window width constant group WMEDFor table 1, section bandwidth constant group W is setSEGFor table 2, subsegment bandwidth constant group W is setBKGFor table 3:
Table (1)
| Frequency values section (unit: MHz) | Section bandwidth constant group WSEG (unit: KHz) |
| 30~80 | 500.0 |
| 80~110 | 1000.0 |
| 110~500 | 1500 |
| 300~500 | 5000 |
| 500~6000 | 40000 |
Table (2)
Table (3)
S2, data block reception and data processing；
Wherein, S21, data block reception: the wide-band frequency spectrum data block B of radio-frequency transmitter output is receivedR, frequency spectrum data
Points are K, and the frequency of data point is fB1~fBK, the level of data point is LB1~LBK；
Wherein；Data block centre frequency is FC, the frequency of data point is fB(k), wherein 1≤k≤K, low-limit frequency LB1
Highest frequency is LBK, the level of data point is LB(k)。
S22, valley point filtering: to frequency spectrum data block BRIn valley point execute median filtering, acquisition eliminate valley point
Frequency spectrum data block BU, median filter length of window is depending on median filtering constant group WMEDWith the frequency spectrum point at spectrum number strong point
Resolution；
Realize the pseudo-program of median filtering are as follows:
For k=2:K-1
if LB(k)<LB(k-1)&&LB(k)<LB(k+1)
To with WMEDFor window, with LB(k)Centered on the r level put be ranked up；
LU(k)Level intermediate value after sequence；
else
LU(k)=LB(k)；
end
end
S23, exponential smoothing: to frequency spectrum data block BUExponential smoothing is executed, filter coefficient β obtains smoothed out frequency
Modal data block BE。
Wherein: the calculation formula of exponential smoothing are as follows:
Wherein:It is the last index sharpening result of record；
LE(k)For frequency spectrum data level；
LU(k)For the frequency spectrum data level after median filtering；
β is filter coefficient；
K is frequency spectrum data points.
S3, data block segmentation: according to the frequency separation f of section bandwidth array WSEG and frequency spectrum data blockB1~fBK, by frequency spectrum
Data block is divided into N number of frequency spectrum data section SEGn, 1≤n≤N；There are P data point, each data point in frequency spectrum data section SEGn
Level be Ls(P)；
S4, subsegment segmentation: in each frequency spectrum data section in SEGn, according to subsegment bandwidth array WBKGIn the wavelength coverage of sum
Frequency of heart fSEG, frequency spectrum data section is divided into M subsegment；
Wherein: M=BWSEG/WEKG[j]；Wherein j is wavelength coverage data block centre frequency fSEGCorresponding WBKGArray index.
S5, subsegment maximum level are sought: the maximum spectral line level sought in M subsegment obtains subsegment maximum level Lmax1~
LmaxM；
S6, wavelength coverage are sought with reference to ambient level: seeking maximum spectral line level L in M subsegmentmax1~LmaxMMinimum value,
It obtains wavelength coverage and refers to ambient level LREF；
S7, wavelength coverage signal detecting threshold are sought: in the data point level value L of frequency spectrum data section SEGnS1~LSPAnd frequency spectrum
Section obtains wavelength coverage signal detecting threshold L with reference to being minimized in ambient level LREFT1~LTPAs output；
It the step of between S8, repetition S4 to S7, is finished until all frequency spectrum data sections are all processed.
Although the present invention is described in detail referring to the foregoing embodiments, for those skilled in the art,
It is still possible to modify the technical solutions described in the foregoing embodiments, or part of technical characteristic is carried out etc.
With replacement, all within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in this
Within the protection scope of invention.
Claims (2)
1. a kind of radio-frequency spectrum signal detecting threshold calculation method, it is characterised in that: steps are as follows:
S1, parameter setting: setting median filter window width constant group WMED, section bandwidth constant group WSEG, subsegment bandwidth constant
Group WBKG, the value and frequency dependence of member node in this three array；Indexing smoothing filter factor beta.
S2, data block reception and data processing；
S3, data block segmentation: according to section bandwidth array WSEGWith the frequency separation f of frequency spectrum data blockB1~fBK, by frequency spectrum data block
It is divided into N number of frequency spectrum data section SEGn, 1≤n≤N；There are P data point, the level of each data point in frequency spectrum data section SEGn
For LS1~LSP；
S4, subsegment segmentation: in each frequency spectrum data section in SEGn, according to subsegment bandwidth array WBKGThe wavelength coverage center frequency of sum
Rate fSEG, frequency spectrum data section is divided into M subsegment；
S5, subsegment maximum level are sought: the maximum spectral line level sought in M subsegment obtains subsegment maximum level Lmax1~LmaxM；
S6, wavelength coverage are sought with reference to ambient level: seeking maximum spectral line level L in M subsegmentmax1~LmaxMMinimum value, obtain
Wavelength coverage refers to ambient level LREF；
S7, wavelength coverage signal detecting threshold are sought: in the data point level value L of frequency spectrum data section SEGnS1~LSPJoin with wavelength coverage
Examine ambient level LREFIn be minimized, obtain wavelength coverage signal detecting threshold LT1~LTPAs output；
It the step of between S8, repetition S4 to S7, is finished until all frequency spectrum data sections are all processed.
2. radio-frequency spectrum signal detecting threshold calculation method according to claim 1, it is characterised in that: in the step S2
Data block reception and data processing include:
S21, data block reception: receiving the wide-band frequency spectrum data block BR of radio-frequency transmitter output, and frequency spectrum data points are K,
The frequency of data point is fB1~fBK, the level of data point is LB1~LBK；
S22, valley point filtering: median filtering is executed to the valley point in frequency spectrum data block BR, acquisition eliminates the frequency of valley point
Modal data block BU, median filter length of window is depending on median filtering constant group WMEDWith the spectral resolution at spectrum number strong point；
S23, exponential smoothing: to frequency spectrum data block BUExponential smoothing is executed, filter coefficient β obtains smoothed out frequency spectrum data
Block BE。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910392957.8A CN110535543B (en) | 2019-05-13 | 2019-05-13 | Radio frequency spectrum signal detection threshold calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910392957.8A CN110535543B (en) | 2019-05-13 | 2019-05-13 | Radio frequency spectrum signal detection threshold calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110535543A true CN110535543A (en) | 2019-12-03 |
| CN110535543B CN110535543B (en) | 2022-01-28 |
Family
ID=68659798
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910392957.8A Active CN110535543B (en) | 2019-05-13 | 2019-05-13 | Radio frequency spectrum signal detection threshold calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110535543B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102759658A (en) * | 2011-04-29 | 2012-10-31 | 特克特朗尼克公司 | Method for automatically setting frequency span in a spectrum analyzer |
| CN103698699A (en) * | 2013-12-06 | 2014-04-02 | 西安交通大学 | Asynchronous motor fault monitoring and diagnosing method based on model |
| US20150257036A1 (en) * | 2013-03-15 | 2015-09-10 | Isco International, Llc | Method and apparatus for avoiding interference |
| CN105492913A (en) * | 2013-05-23 | 2016-04-13 | 滨特尔热能管理有限责任公司 | Sub-harmonic arc fault detection system and method |
| CN109474355A (en) * | 2018-01-17 | 2019-03-15 | 国家无线电频谱管理研究所有限公司 | Adaptive noise THRESHOLD ESTIMATION and method for extracting signal based on spectrum monitoring data |
-
2019
- 2019-05-13 CN CN201910392957.8A patent/CN110535543B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102759658A (en) * | 2011-04-29 | 2012-10-31 | 特克特朗尼克公司 | Method for automatically setting frequency span in a spectrum analyzer |
| US20150257036A1 (en) * | 2013-03-15 | 2015-09-10 | Isco International, Llc | Method and apparatus for avoiding interference |
| CN105492913A (en) * | 2013-05-23 | 2016-04-13 | 滨特尔热能管理有限责任公司 | Sub-harmonic arc fault detection system and method |
| CN103698699A (en) * | 2013-12-06 | 2014-04-02 | 西安交通大学 | Asynchronous motor fault monitoring and diagnosing method based on model |
| CN109474355A (en) * | 2018-01-17 | 2019-03-15 | 国家无线电频谱管理研究所有限公司 | Adaptive noise THRESHOLD ESTIMATION and method for extracting signal based on spectrum monitoring data |
Non-Patent Citations (1)
| Title |
|---|
| 齐佩汉: "基于功率谱分段对消频谱感知算法研究及性能分析", 《电子与信息学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110535543B (en) | 2022-01-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100518012C (en) | Authorization user signal detecting method for cognitive radio system | |
| CN105182070B (en) | A kind of signal detecting method | |
| CN111157953B (en) | Two-stage threshold constant false alarm detection algorithm under strong ground clutter | |
| CN111510255B (en) | Frequency hopping signal blind detection and parameter estimation method based on broadband frequency spectrum data | |
| KR102221839B1 (en) | Real-time cardiac rate detection apparatus in noisy environment and method thereof | |
| CN109274438B (en) | Self-adaptive double-threshold cooperative interference detection method and device | |
| CN112684251A (en) | Target signal frequency domain detection method based on power spectrum template | |
| WO2016060620A1 (en) | Frame based spike detection module | |
| CN103778921A (en) | Method for eliminating nonuniform noise in speech collected by radar | |
| CN108847910B (en) | Spectrum sensing method and device and spectrum sensing equipment | |
| CN106656372B (en) | Frequency band interference detection method of frequency hopping system | |
| CN113447893B (en) | Radar pulse signal frequency spectrum automatic detection method, system and medium | |
| CN110535543A (en) | A kind of radio-frequency spectrum signal detecting threshold calculation method | |
| CN109596354A (en) | Band-pass filtering method based on the identification of adaptive resonance frequency band | |
| US8175829B2 (en) | Analyzer for signal anomalies | |
| CN106073822B (en) | A kind of superonic spectrum adaptive noise estimation envelope extraction method | |
| CN109004996B (en) | Signal detection method based on multi-sine-window power spectrum peak value | |
| Aziz et al. | Spectrum sensing for cognitive radio using multicoset sampling | |
| Al-Badrawi et al. | An EMD-based double threshold detector for spectrum sensing in cognitive radio networks | |
| CN113824518A (en) | Non-cooperative signal detection method based on frequency spectrum | |
| CN113054945B (en) | Effective excitation detection method of surface acoustic wave resonator based on entropy analysis | |
| CN111147168B (en) | Signal detection method with power spectrum and statistics fused | |
| Hezarkhani et al. | Performance analysis of a CA-CFAR detector in the interfering target and homogeneous background | |
| CN112504429A (en) | High-precision demodulation algorithm for strong interference DVS | |
| CN110830130A (en) | Statistical frequency spectrum detection method in cognitive radio |
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 |