CN106656887A - Method of estimating signal center frequency point and bandwidth - Google Patents

Method of estimating signal center frequency point and bandwidth Download PDF

Info

Publication number
CN106656887A
CN106656887A CN201610913626.0A CN201610913626A CN106656887A CN 106656887 A CN106656887 A CN 106656887A CN 201610913626 A CN201610913626 A CN 201610913626A CN 106656887 A CN106656887 A CN 106656887A
Authority
CN
China
Prior art keywords
spike
frequency
signal
peak
frequency spectrum
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
Application number
CN201610913626.0A
Other languages
Chinese (zh)
Other versions
CN106656887B (en
Inventor
王剑
王超
张旭东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201610913626.0A priority Critical patent/CN106656887B/en
Publication of CN106656887A publication Critical patent/CN106656887A/en
Application granted granted Critical
Publication of CN106656887B publication Critical patent/CN106656887B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/01Equalisers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0014Carrier regulation
    • H04L2027/0024Carrier regulation at the receiver end
    • H04L2027/0026Correction of carrier offset

Abstract

The invention discloses a method of estimating a signal center frequency point and a bandwidth. The method comprises steps: a receiver antenna is measured, and a gain curve for the antenna in a receiving frequency band is obtained; the received signals at each time are subjected to Fourier transform to obtain a frequency spectrum, and the frequency spectrum is compensated according to the gain curve; the average value for frequency spectrum absolute values after multiple times of received signal compensation is calculated, a first frequency spectrum is obtained through smoothing processing, a second frequency spectrum is obtained through left-right turning and then smoothing filtering, and by such analogy, an Nth frequency spectrum is obtained, wherein N is an even number; positions and peak values of all peaks in the Nth frequency spectrum are acquired; as for any peak, a first position and a second position when the Nth frequency spectrum falls to a peak setting threshold multiple for the first time to the left and to the right are found out, a peak extending interval length is thereafter calculated, and peaks whose extending intervals are zero are removed; and the bandwidth and the center frequency point of any left peak are calculated. The method of the invention has the following advantage that the signal center frequency point and bandwidth estimation difficulty in a low-signal-to-noise ratio environment can be solved.

Description

A kind of evaluation method of signal center's frequency and bandwidth
Technical field
The invention belongs to wireless communication technology field, more particularly to the evaluation method of a kind of signal center's frequency and bandwidth.
Background technology
In wireless communication technology field, receiver is possible to receiving in the range of signal bandwidth while receiving several not Same signal.These signals may adopt different modulation systems, with different carrier frequency and signal bandwidth.Therefore, to this A little signals are carried out before parameter Estimation, Modulation Identification and demodulation it may first have to estimate the center frequency point and bandwidth of signal.
Existing signal bandwidth method of estimation is all directly to be estimated on the frequency spectrum or power spectrum of signal, due to letter Number spectrum itself rises and falls violent, and burr is a lot, is difficult to distinguish signal with noise in the case where signal to noise ratio is relatively low.
The content of the invention
Present invention seek to address that above-mentioned technical problem.
For this purpose, it is an object of the invention to propose the evaluation method of a kind of signal center's frequency and bandwidth, by signal Spectral smoothing filtering is solving the problems, such as that spectral ripple is violent, burr is a lot;By peak value searching accurately determining in signal Heart frequency and bandwidth, so as to solve signal center's frequency and a bandwidth estimation difficult problem under low signal-to-noise ratio environment.
To achieve these goals, embodiment of the invention discloses that the estimation side of a kind of signal center's frequency and bandwidth Method, comprises the following steps:S1:Receiver antenna is measured, obtains receiving the gain curve of frequency range internal antenna;S2:To every The secondary signal Fourier transformation for receiving obtains receiving the frequency spectrum of signal;According to the gain curve of the antenna to receiving every time Signal compensate, to eliminate impact of the antenna gain to signal spectrum;S3:Step S2 is performed a plurality of times, after taking repeatedly compensation The mean value of the absolute value of signal;S4:The mean value of the absolute value of signal after the multiple compensation is carried out into smooth filter for the first time Ripple obtains the first frequency spectrum, and left and right overturns, and again smothing filtering obtains the second frequency spectrum, and N frequency spectrums are obtained by that analogy, and wherein N is Even number;S5:The position that all spike peak values in the N frequency spectrums and all spikes are located is obtained, for arbitrary spike, according to Given threshold starts respectively to scan for N frequency spectrums to from left to right from spike position, finds N frequency spectrums and drops to first The position of given threshold multiple described in spike peak value, is designated as first position and the second place, according to the first position and described The second place calculates the extension siding-to-siding block length of spike, removes the interval null spike of extension;Calculating arbitrarily remains spike Center frequency point and bandwidth.
Signal center's frequency according to embodiments of the present invention and the evaluation method of bandwidth, by signal spectrum smothing filtering To solve the problems, such as that spectral ripple is violent, burr is a lot;By peak value searching accurately to determine the center frequency point and bandwidth of signal, So as to solve signal center's frequency and a bandwidth estimation difficult problem under low signal-to-noise ratio environment.
In addition, the evaluation method of signal center's frequency according to the above embodiment of the present invention and bandwidth, can also have such as Lower additional technical characteristic:
Further, the receiver antenna is measured using vector network analyzer.
Further, the frequency according to the gain curve of the antenna and the frequency spectrum of the reception signal to the reception signal Spectrum is compensated.
Further, step S4 is further included:S401:F (k) is obtained into Y (k) by smoothing filter, wherein, F (k) For the mean value of the absolute value of signal after compensation, Y (k) is the frequency spectrum obtained through smoothing filter process;S402:Will be flat Frequency spectrum Y (k) after cunning carries out left and right upset and obtains Y ' (k), and Y ' (k) is assigned to Y (k) again;S403:Repeat S401 N frequency spectrums are obtained with S402.
Further, step S5 is further included:Find all spikes peak value Peak (j) and the spike institute of smooth rear frequency spectrum Position Posi (j), j=1,2 ..., K, K for spike number, for any one spike Peak (j):It is located from spike Position Posi (j) start, respectively to from left to right search first less than C × Peak (j) point, its position is designated as respectively BandL (j) and bandR (j), wherein C is the threshold value of setting, 0<C<1;If spike left side is not searched less than C × Peak (j) Point, or spike left side search procedure in occur in that frequency spectrum value more than spike peak (j) situation, make bandL (j) be equal to posi(j);If spike right side does not search in the point less than C × Peak (j), or spike right side search procedure and occurs in that frequency Spectrum value then makes bandR (j) be equal to posi (j) more than the situation of spike peak (j).
Further, if bandL (j) is equal to posi (j) and bandR (j) is not equal to posi (j), bandR is made J () is equal to posi (j);If bandR (j) is equal to posi (j) and bandL (j) is not equal to posi (j), bandL (j) is made Equal to posi (j).
Further, step S6 is further included:Siding-to-siding block length bandW (j) for calculating spike peak (j) extension is equal to BandR (j) deducts bandL (j);Remove the null spikes of bandW (j);For the spike for arbitrarily remaining, with reference to frequency sweep Frequency band range, center frequency point and bandwidth of the interval computation spike to induction signal is extended according to spike.
The additional aspect and advantage of the present invention will be set forth in part in the description, and partly will become from the following description Obtain substantially, or recognized by the practice of the present invention.
Description of the drawings
The above-mentioned and/or additional aspect and advantage of the present invention will become from the description with reference to accompanying drawings below to embodiment It is substantially and easy to understand, wherein:
Fig. 1 is the flow chart of the evaluation method of signal center's frequency and bandwidth of one embodiment of the invention.
Fig. 2 is 200 frequency sweeps are obtained in 860MHz~880MHz frequency bands signal through Fourier transformation and after compensating Absolute spectra mean value.
Fig. 3 is 12nd frequency spectrum of the frequency spectrum of Fig. 2 after multiple smothing filtering and upset.
Fig. 4 is that the rear all of spike of frequency spectrum and the corresponding extension interval of all spikes are smoothed shown in Fig. 3.
Fig. 5 is to remove in Fig. 4 spike and extend siding-to-siding block length to be after zero spike, the final center frequency point for obtaining and bandwidth The result of estimation.
Fig. 6 is 200 frequency sweeps are obtained in 950MHz~970MHz frequency bands signal through Fourier transformation and after compensating Absolute spectra mean value.
Fig. 7 is 12nd frequency spectrum of the frequency spectrum of Fig. 6 after multiple smothing filtering and upset.
Fig. 8 is that the rear all of spike of frequency spectrum and the corresponding extension interval of all spikes are smoothed shown in Fig. 7.
Fig. 9 is to remove in Fig. 8 spike and extend siding-to-siding block length to be after zero spike, the final center frequency point for obtaining and bandwidth The result of estimation.
Specific embodiment
Embodiments of the invention are described below in detail, the example of the embodiment is shown in the drawings, wherein from start to finish Same or similar label represents same or similar element or the element with same or like function.Below with reference to attached The embodiment of figure description is exemplary, is only used for explaining the present invention, and is not considered as limiting the invention.
In describing the invention, it is to be understood that term " first ", " second " are only used for describing purpose, and can not It is interpreted as indicating or implying relative importance.
With reference to explained below and accompanying drawing, it will be clear that these and other aspects of embodiments of the invention.In these descriptions In accompanying drawing, specifically disclose some particular implementations in embodiments of the invention to represent the enforcement for implementing the present invention Some modes of the principle of example, but it is to be understood that the scope of embodiments of the invention is not limited.Conversely, the present invention Embodiment includes all changes, modification and the equivalent fallen in the range of the spirit and intension of attached claims.
Below in conjunction with the Description of Drawings present invention.
Fig. 1 is the flow chart of the evaluation method of signal center's frequency and bandwidth of one embodiment of the invention.
S1:Receiver antenna is measured, obtains receiving the gain curve of frequency range internal antenna.
In one embodiment of the invention, receiver antenna is measured using vector network analyzer.Vector net Network analyzer can be while measuring amplitude information and phase information, can effectively suppress interference and spuious, dynamic range big.
S2:Signal to receiving every time does Fourier transformation and obtains receiving the frequency spectrum of signal;It is bent according to the gain of antenna Line is compensated to the frequency spectrum for obtaining, to eliminate impact of the antenna gain to signal spectrum.
In one embodiment of the invention, signal compensation is carried out using below equation:
F′m(k)=Fm(k)/s21(k), k=1,2 ..., K
Wherein, FmK () is described the m time receives the frequency spectrum of signal, F 'mK () is the frequency spectrum after compensation, K is reception signal Frequency spectrum is counted, s21K () is the gain curve of the antenna.
S3:Step S2 is performed a plurality of times, the mean value of the absolute value of signal after repeatedly compensation is taken.
In one embodiment of the invention, the mean value of the absolute value of signal after compensation is calculated by below equation:
Wherein, M is number of repetition, and F (k) is the mean value of the absolute value of signal after compensation.Frequency is centered on Fig. 2 870Hz, the signal obtained for 200 times with frequency sweep in a width of 20MHz (i.e. 860MHz~880MHz) frequency band passes through frequency spectrum after compensation The mean value of absolute value.As can be seen that 860MHz~880MHz bandwidth memories are in several narrow band signals and a broadband signal, And this broadband signal is likely to be and is made up of multiple subbands;Due to the impact of noise so that the frequency spectrum of signal exists substantial amounts of Burr.
S4:The mean value of the absolute value of signal after repeatedly compensation is carried out smoothing for the first time and obtains the first frequency spectrum, left and right is turned over Turn, again smothing filtering obtains the second frequency spectrum, obtain N frequency spectrums by that analogy, wherein N is even number and for adjustable parameter.
In one embodiment of the invention, step S4 is further included:
S401:Smoothing filter by F (k) by a length for M obtains Y (k), and a kind of computational methods of Y (k) are:
S402:To prevent smothing filtering from causing center frequency point to shift, frequency spectrum Y (k) after smoothing carries out left and right and turns over Turn to obtain Y ' (k), and Y ' (k) is assigned to Y (k) again;
S403:Repeat S401 and S402 obtains N frequency spectrums.
Fig. 3 gives the 12nd frequency spectrum obtained through 12 smothing filterings and upset of frequency spectrum shown in Fig. 2.As can be seen that Frequency spectrum after smothing filtering is affected to significantly reduce by burr.
S5:The position that all spike peak values in N frequency spectrums and all spikes are located is obtained, for arbitrary spike, according to Given threshold starts respectively to scan for N frequency spectrums to from left to right from spike position, finds N frequency spectrums and drops to first The position of given threshold multiple described in spike peak value, is designated as first position and the second place, according to first position and the second place The extension siding-to-siding block length of spike is calculated, removes the interval null spike of extension.Fig. 4 gives all points of frequency spectrum of Fig. 3 the 12nd The extension siding-to-siding block length at peak and spike, wherein thicker rectangle frame represents that the interval spike more than zero of extension is corresponding prolonging Interval is stretched, other vertical lines represent the corresponding extension interval for extending interval null spike (in order to represent extension in figure Interval be zero the corresponding extension of spike it is interval, mapping season extends the corresponding extension siding-to-siding block length of interval null spike For 1).
In one embodiment of the invention, step S5 is further included:Find all spike peak values of smooth rear frequency spectrum Position Posi (j) that Peak (j) and spike are located, j=1,2 ..., K, K is the number of spike, for any one spike Peak (j):
From the beginning of position Posi (j) being located from spike, the point of C × Peak (j) is less than first to from left to right search respectively, its Position is designated as respectively bandL (j) and bandR (j), and wherein C is the threshold value of setting, and 0<C<1.If spike left side is not searched for Feelings of the frequency spectrum value more than spike peak (j) are occurred in that in the point less than C × Peak (j), or spike left side search procedure Condition, makes bandL (j) be equal to posi (j);Point less than C × Peak (j), or spike right side are not searched by spike right side Situation of the frequency spectrum value more than spike peak (j) is occurred in that in search procedure, then makes bandR (j) be equal to posi (j).
To exclude many inspection situations that the frequency spectrum left and right sides is likely to occur, if bandL (j) is equal to posi (j) and bandR J () is not equal to posi (j), then make bandR (j) be equal to posi (j);If bandR (j) is equal to posi (j) and bandL (j) It is not equal to posi (j), then makes bandL (j) be equal to posi (j).
S6:For the spike for arbitrarily remaining, with reference to the frequency band range of frequency sweep, according to the spike extension interval computation point Center frequency point and bandwidth of the peak to induction signal.
In one embodiment of the invention, step S6 is further included:Calculate the siding-to-siding block length that spike peak (j) extends BandW (j)=bandR (j)-bandL (j), removes null those spikes of bandW (j).Assume the frequency band that receiver is received Scope is [WL,WR], then for the spike for remaining, bandwidth bandwidth (j) is calculated according to below equation:
Center frequency point freqcenter(j)It is calculated according to below equation:
Fig. 5 gives center frequency point and the bandwidth that all estimations of frequency spectrum shown in Fig. 2 are obtained.Rectangle frame therein is with regard to table Show the center frequency point and bandwidth that each signal is detected in 860MHz~880MHz bandwidth.Broadband signal therein is divided into The narrow band signal of 4 relative narrowers, this correspondence two kinds of situation in practice, the first situation is:860MHz~880MHz is true Real storage has the broadband signal that multiple subbands are constituted at one;Second situation be:It is implicitly present in the range of 860MHz~880MHz One broadband signal, but because signal to noise ratio is too low, the impact of noise is excessive, causes it to have been partitioned into multiple relative narrowers Narrow band signal.
In order to further illustrate the high efficiency of algorithm involved in the present invention, Fig. 6~Fig. 9 gives signal according to the present invention Center and another group of measured result of Survey of Bandwidth Estimation Scheme.The frequency sweep center frequency point for arranging receiver is 960MHz, and Fig. 6 gives The mean value of absolute value of 200 receipt signal frequency spectrums after compensation, it can be seen that exist two in the frequency band range Individual narrow band signal.Fig. 7 gives 12nd frequency spectrum of the frequency spectrum described in Fig. 6 after 12 smothing filterings and upset, it can be seen that Frequency spectrum obvious reduction affected by noise after smoothing.Fig. 8 gives all spikes of the 12nd frequency spectrum described in Fig. 7 and right The spike answered extends interval, it can be seen that have two to extend the interval spike more than zero and several interval null points of extension Peak.Fig. 9 gives to remove and extend described in Fig. 8 after interval null spike, signal center's frequency for eventually detecting and Correspondence bandwidth.
In addition, the evaluation method of signal center's frequency of the embodiment of the present invention and bandwidth other constitute and act on for All it is known for those skilled in the art, in order to reduce redundancy, does not repeat.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show The description of example " or " some examples " etc. means to combine specific features, structure, material or spy that the embodiment or example are described Point is contained at least one embodiment of the present invention or example.In this manual, to the schematic representation of above-mentioned term not Necessarily refer to identical embodiment or example.And, the specific features of description, structure, material or feature can be any One or more embodiments or example in combine in an appropriate manner.
Although an embodiment of the present invention has been shown and described, it will be understood by those skilled in the art that:Not These embodiments can be carried out with various changes, modification, replacement and modification in the case of the principle and objective that depart from the present invention, this The scope of invention is limited by claim and its equivalent.

Claims (7)

1. the evaluation method of a kind of signal center's frequency and bandwidth, it is characterised in that comprise the following steps:
S1:Receiver antenna is measured, obtains receiving the gain curve of frequency range internal antenna;
S2:Signal to receiving every time does Fourier transformation and obtains receiving the frequency spectrum of signal;It is bent according to the gain of the antenna Line is compensated to the frequency spectrum for obtaining, to eliminate impact of the antenna gain to signal spectrum;
S3:Step S2 is performed a plurality of times, the mean value of the absolute value of signal spectrum after repeatedly compensation is taken;
S4:The mean value of the absolute value of signal after the multiple compensation is carried out into first time smothing filtering and obtains the first frequency spectrum, it is left Right upset, again smothing filtering obtain the second frequency spectrum, obtain N frequency spectrums by that analogy, wherein N is even number;
S5:The position that all spike peak values and all spikes are located in the N frequency spectrums is obtained, for arbitrary spike, according to setting Determine threshold value to start respectively to scan for N frequency spectrums to from left to right from spike position, find N frequency spectrums and drop to point first The position of given threshold multiple described in peak-to-peak value, is designated as first position and the second place, according to the first position and described The extension siding-to-siding block length of two position calculation spikes, removes the interval null spike of extension;
S6:Calculating arbitrarily remains the center frequency point and bandwidth of spike.
2. the evaluation method of signal center's frequency according to claim 1 and bandwidth, it is characterised in that in step sl, The receiver antenna is measured using vector network analyzer.
3. the evaluation method of signal center's frequency according to claim 1 and 2 and bandwidth, it is characterised in that according to described The frequency spectrum of the gain curve of antenna and the reception signal is compensated to the frequency spectrum of the reception signal.
4. the evaluation method of signal center's frequency according to claim 1 and bandwidth, it is characterised in that step S4 is further Including:
S401:F (k) is obtained into Y (k) by a smoothing filter, wherein, F (k) is the average of the absolute value of signal after compensation Value, Y (k) is the frequency spectrum obtained through smoothing filter process;
S402:Y (k) is carried out into left and right upset and obtains Y ' (k), and Y ' (k) is assigned to Y (k) again;
S403:Repeat S401 and S402 obtains N frequency spectrums.
5. the evaluation method of signal center's frequency according to claim 1 and bandwidth, it is characterised in that step S5 is further Including:Position Posi (j) that all spikes peak value Peak (j) and spike for finding smooth rear frequency spectrum are located, j=1,2 ..., K, K For the number of spike, for any one spike Peak (j):
From the beginning of position Posi (j) being located from spike, the point of C × Peak (j) is less than first to from left to right search respectively, its place Position is designated as respectively bandL (j) and bandR (j), and wherein C is the threshold value of setting, and 0<C<1;If spike left side does not search low Situation of the frequency spectrum value more than spike peak (j), order are occurred in that in the point of C × Peak (j), or spike left side search procedure BandL (j) is equal to posi (j);If spike right side does not search the point less than C × Peak (j), or spike right side was searched for Situation of the frequency spectrum value more than spike peak (j) is occurred in that in journey, then makes bandR (j) be equal to posi (j).
6. the evaluation method of signal center's frequency according to claim 5 and bandwidth, it is characterised in that
If bandL (j) is equal to posi (j) and bandR (j) is not equal to posi (j), bandR (j) is made to be equal to posi (j); If bandR (j) is equal to posi (j) and bandL (j) is not equal to posi (j), bandL (j) is made to be equal to posi (j).
7. the evaluation method of signal center's frequency according to claim 6 and bandwidth, it is characterised in that step S6 is further Including:
Siding-to-siding block length bandW (j) for calculating spike peak (j) extension deducts bandL (j) equal to bandR (j);Remove bandW (j) null spike;
For the spike for arbitrarily remaining, with reference to the frequency band range of frequency sweep, according to spike extension interval computation spike correspondence The center frequency point and bandwidth of signal.
CN201610913626.0A 2016-10-19 2016-10-19 A kind of evaluation method of signal center's frequency point and bandwidth Active CN106656887B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610913626.0A CN106656887B (en) 2016-10-19 2016-10-19 A kind of evaluation method of signal center's frequency point and bandwidth

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610913626.0A CN106656887B (en) 2016-10-19 2016-10-19 A kind of evaluation method of signal center's frequency point and bandwidth

Publications (2)

Publication Number Publication Date
CN106656887A true CN106656887A (en) 2017-05-10
CN106656887B CN106656887B (en) 2019-08-09

Family

ID=58856380

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610913626.0A Active CN106656887B (en) 2016-10-19 2016-10-19 A kind of evaluation method of signal center's frequency point and bandwidth

Country Status (1)

Country Link
CN (1) CN106656887B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113726348A (en) * 2021-07-21 2021-11-30 湖南艾科诺维科技有限公司 Smoothing filtering method and system for radio signal frequency spectrum
CN115276693A (en) * 2022-09-26 2022-11-01 中孚信息股份有限公司 Multi-frequency-point multi-bandwidth identification method based on welch

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104777355A (en) * 2015-04-21 2015-07-15 上海昕石电子科技有限公司 Method and device for quickly detecting frequency
WO2016045388A1 (en) * 2014-09-23 2016-03-31 华为技术有限公司 Self-adaptive narrow band interference cancellation method and apparatus

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016045388A1 (en) * 2014-09-23 2016-03-31 华为技术有限公司 Self-adaptive narrow band interference cancellation method and apparatus
CN104777355A (en) * 2015-04-21 2015-07-15 上海昕石电子科技有限公司 Method and device for quickly detecting frequency

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈源等: "利用IIR自适应频率估计器实现多频点估计", 《电讯技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113726348A (en) * 2021-07-21 2021-11-30 湖南艾科诺维科技有限公司 Smoothing filtering method and system for radio signal frequency spectrum
CN113726348B (en) * 2021-07-21 2022-06-21 湖南艾科诺维科技有限公司 Smoothing filtering method and system for radio signal frequency spectrum
CN115276693A (en) * 2022-09-26 2022-11-01 中孚信息股份有限公司 Multi-frequency-point multi-bandwidth identification method based on welch

Also Published As

Publication number Publication date
CN106656887B (en) 2019-08-09

Similar Documents

Publication Publication Date Title
EP3259613B1 (en) System and method for estimating time of arrival (toa)
US10270626B2 (en) Adaptive narrowband interference cancellation method, and apparatus
CN102843328B (en) Method for detecting PRACH (Physical Random Access Channel) sequence in TDD LTE (Time Division Duplex Long Term Evolution) system
RU2593146C1 (en) Method for adaptive detection of signals of moving targets on background of multicomponent passive interference
CN114025379B (en) Broadband multi-signal detection method, device and equipment
CN105897643B (en) Based on joint narrowband Interference Detection, elimination and the channel estimation methods of continuous cyclic shift training sequence in ofdm system
CN110837079B (en) Target detection method and device based on radar
WO2016035439A1 (en) Interference identification device, radio communication device, and interference identification method
CN113259029B (en) Real-time automatic detection and identification method suitable for unmanned aerial vehicle signals
CN106656887A (en) Method of estimating signal center frequency point and bandwidth
CN112684251A (en) Target signal frequency domain detection method based on power spectrum template
CN106093896B (en) A kind of quick detection and the method for estimating pulse Doppler radar pulse jitter
US8804691B1 (en) Phase based preamble detection
CN113447893B (en) Radar pulse signal frequency spectrum automatic detection method, system and medium
US20160192217A1 (en) Fast signal surveyor
CN110020677B (en) Continuous passerby number detection method based on WiFi Doppler frequency shift
CN104363065B (en) The wireless communication system frequency spectrum sensing method estimated based on non-Gaussian system
CN108718223B (en) Blind spectrum sensing method for non-cooperative signals
JP5463534B2 (en) Ionogram ionospheric echo automatic polarization separation and automatic reading method
KR101652193B1 (en) Target detection method and apparatus using RADAR
CN110233687A (en) A kind of multichannel data combined signal judgement detection method
KR101694521B1 (en) Apparatus and method for generating radio fingerprint map
CN111641934B (en) Method for estimating electromagnetic signal situation distribution based on vehicle-mounted equipment monitoring data
CN116359851A (en) Radar active interference detection and identification method and device based on converged network
DE102021109930A1 (en) Communication device and position estimation method

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