CN110519003A - A kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference - Google Patents

A kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference Download PDF

Info

Publication number
CN110519003A
CN110519003A CN201910808275.0A CN201910808275A CN110519003A CN 110519003 A CN110519003 A CN 110519003A CN 201910808275 A CN201910808275 A CN 201910808275A CN 110519003 A CN110519003 A CN 110519003A
Authority
CN
China
Prior art keywords
signal
frequency
unmanned plane
detection
time
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
CN201910808275.0A
Other languages
Chinese (zh)
Other versions
CN110519003B (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.)
Chinese people's liberation army 63892
Original Assignee
Chinese people's liberation army 63892
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 Chinese people's liberation army 63892 filed Critical Chinese people's liberation army 63892
Priority to CN201910808275.0A priority Critical patent/CN110519003B/en
Publication of CN110519003A publication Critical patent/CN110519003A/en
Application granted granted Critical
Publication of CN110519003B publication Critical patent/CN110519003B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04KSECRET COMMUNICATION; JAMMING OF COMMUNICATION
    • H04K3/00Jamming of communication; Counter-measures
    • H04K3/60Jamming involving special techniques
    • H04K3/62Jamming involving special techniques by exposing communication, processing or storing systems to electromagnetic wave radiation, e.g. causing disturbance, disruption or damage of electronic circuits, or causing external injection of faults in the information
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention discloses a kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference comprising following sequential steps: the antenna of broadband receiver receives electromagnetic wave signal, amplifies, filters, being mixed and IF process;Treated signal carries out carrier frequency measurement, carries out pulse detection, the pulse repetition period of measuring signal and pulse width parameter in conjunction with narrow-band filtering, with carry out unmanned plane signal whether there is or not anticipation;According to carrier frequency, pulse repetition period and width parameter that measurement obtains, unmanned plane signal can be judged whether there is;If judging, no unmanned plane signal is returned to previous step and continues to test, if tentatively determining whether man-machine signal, detects judgement with regard to the uplink and downlink signals carried out in next step;The uplink signal and downlink signal of unmanned plane are distinguished according to the bandwidth of unmanned plane signal and frequency hopping rate.Concealment of the present invention is strong, does not generate interference, being capable of round-the-clock, all weather operations.

Description

A kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference
Technical field
The present invention relates to electronic countermeasure technology fields, above and below a kind of unmanned plane based on signal characteristic difference Row communication link recognition methods.
Background technique
In recent years, the technological break-through along with unmanned plane on low latitude and miniaturization, all kinds of civilian unmanned planes are in the world Development upsurge is started in range.Unmanned plane is in logistics transportation, geologic prospect, movies-making, agrisilviculture, patrol monitoring, emergency The application demand rapid growth of rescue etc., attracts more and more science-and-technology enterprises to participate in emerging unmanned plane industry.With This simultaneously, consumer level civil micro unmanned plane, due to its is easy to operate, price is human-oriented, feeling of freshness and entertainment it is strong the features such as, it is general All over for it is big it is well-established with pursue, it is growing day by day to buy user.In future, inevitable more wide, the city of the application field of unmanned plane Field can not estimate.
However as the development of industry, accident relevant to unmanned plane is also frequently occurred among the visual field of people.2017 In April in year, Chengdu airport half a month, interior several unmanned planes of generation interfered flight events, caused society and discussed warmly;2015, U.S.'s law enforcement Department finds that a frame unmanned plane swarms into the White House;In addition, there is the phenomenon that carrying out criminal activity using unmanned plane in foreign countries, it is above each The behavior of the improper operation unmanned plane of kind belongs to " black to fly " in itself." black fly " unmanned plane caused various circles of society's concern and Worry, the demand effectively managed to it are very urgent.However, China relatively lags behind in management unmanned plane industry, mainly Be embodied in: operator lacks the training of system, and the development of specification unmanned plane, sale, the relevant law that uses are perfect not to the utmost, detection prison It is immature to control management technical means.
In face of " black fly " unmanned plane, not only specification is set, policies accomplish " not allowing winged ", it will also be on technological layer It studies anti-means and realizes " dare not fly ".Unmanned plane is broken through, high-power electromagnetic interference is relied primarily at present and destroys nothing Man-machine communication link.The communication link of unmanned plane includes uplink and downlink, wherein uplink signal is mainly used It is transmitted in information such as flight control instructions, down link signal is mainly used for the information such as drone status parameter, video image biography It is defeated.For the military unmanned air vehicle of high threat level, the working frequency of uplink and downlink link is inconsistent.It is breaking through During unmanned plane, the emphasis of high-power electromagnetic interference is the control instruction signal of uplink, threatens people to cut off Operation and control of the member to unmanned plane.If cannot the working frequency points of uplink be carried out with accurately frequency sweep compacting interference, arrive Up to unmanned plane electromagnetic interference signal equivalent radiation power with regard to inadequate, can not just play counter unmanned plane effect.Therefore, needle To the technological means of electromagnetism compacting interference, most primary problem is exactly the uplink and downlink link signal for perceiving unmanned plane, and area Divide uplink and downlink link.
Currently, carrying out the technical method of Detection location to unmanned plane mainly includes active radar detection localization method, photoelectricity Recognition and tracking method, passive acoustics detectiona localization method are detected, the radar target that above-mentioned each method passes through sensorcraft respectively returns The sound that wave, unmanned plane Infrared Image Features, unmanned plane motor issue is detected, and unmanned plane can not be detected and determine Uplink and downlink signal, therefore have no idea accurately to guide high-power electromagnetic interference counter unmanned plane, can only can according to unmanned plane Workable 5~6 frequency dot cycles transmitting, while electromagnetic interference is carried out, therefore interfering signal power is difficult to concentrate, interference is anti- The effect is unsatisfactory for system.
Summary of the invention
To solve the above problems, the object of the present invention is to provide a kind of unmanned plane uplink and downlinks based on signal characteristic difference to lead to Believe link recognition, small drone uplink and downlink communication chain under the conditions of non-cooperating is identified using radio detection means Road guides accurate interference, blocks the unmanned aerial vehicle (UAV) control signal of uplink, provides condition for next step unmanned plane adapter tube;Above-mentioned Small drone refers to not providing the noncooperative target of identity information, unknown transmission information under non-cooperating condition of acceptance.
For achieving the above object, the present invention adopts the following technical scheme:
A kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference comprising following steps:
S1, broadband receiver antenna receive electromagnetic wave signal, amplify, filter, being mixed and IF process;
S2, carrier frequency CF measurement is carried out to step S1 treated signal, carries out pulse detection in conjunction with narrow-band filtering, surveys Measure signal pulse repetition period PRI and pulse width PW parameter, with carry out unmanned plane signal whether there is or not anticipation;
S3, determine whether unmanned plane signal
Carrier frequency CF, the pulse repetition period PRI and pulse width PW parameter obtained according to step S2 measurement, can Judge whether there is unmanned plane signal;If judging, no unmanned plane signal is returned to step S2 and continues to test, if tentatively determining whether Man-machine signal detects judgement with regard to the uplink and downlink signals carried out in step S4;
S4, uplink and downlink signals detection judgement is carried out to unmanned plane signal
After identifying unmanned plane signal, bandwidth and frequency hopping rate to the unmanned plane signal carry out precise measurement and estimate Meter;
Firstly, accurately estimating the bandwidth of signal, perform the following operation:
1) power spectrum is sought
The data x (n) (n=0,1, L, N-1) for being N by length, is divided into L sections, and every section has M data, and wherein L is not to be overlapped When each data segment length, L=N/M, the i-th segment data xi(n)=x (n+iM-M) is then added to window function ω (n) each On data segment, each section of cyclic graph, i-th section of cyclic graph are found out
In formula, U is normalization factor,
Will be approximate irrelevant between each section of cyclic graph, last power Spectral Estimation is
2) wavelet decomposition obtains smooth power spectrum
Wavelet decomposition is carried out to the power spectrum that estimates, is partially separated by the detail section of signal and roughly, then extract slightly It omits the coefficient of part and power spectrum is reconstructed with this coefficient, eliminate the high frequency detail in power spectrum waveform, obtain smooth waveform;
3) the mobile covariance of smooth power spectrum is calculated
After obtaining smooth power spectrum, the starting point and ending point of passband is extracted to estimate bandwidth;Calculate mobile covariance:
R (k, k+1)=cov (d (k), d (k+1))
(k=1,2, L, N-1) (15)
In formula: d (k) is the reconstruction coefficients that rough part is extracted after wavelet decomposition, and the formula is utilized to traverse all data Point, what is found out is the mobile covariance between arbitrary neighborhood point;When every p, value can just change, so using Following formula calculates
In formula, k=1,1+p, L, N-p, the selection of p value, as p=1, are equivalent to by determining the characteristics of actual data point Formula;
4) estimate bandwidth
The position a and b where maximum 2 mobile covariance values are extracted, and is regarded as the starting point of bandwidth respectively And cut off, by difference | b-a-1 | as estimation bandwidth;
Cycle-index is set, is repeated the above steps 1)~4 to the L section for the data that length is N), L bandwidth is calculated, then ask It is average, find out the assembly average of estimation bandwidth, as the precise bandwidth estimated value of unmanned plane signal;
Then, accurately estimate frequency hopping rate, perform the following operation: analyzing signal using Short Time Fourier Transform, obtain at that time Frequency indicates, recycles wavelet transformation to extract the side information of the time-frequency representation, and then estimate frequency hopping rate using spectrum analysis.
Further, in above-mentioned step S2, the method for unmanned plane signal detection anticipation includes the following steps:
S2a, carrier frequency CF estimation
Rough estimate first is carried out to signal carrier frequency using spectrum peak detection calculations, determines frequency range;Then it recycles ZFFT further refines frequency spectrum;Finally " fence effect " bring frequency deviation is eliminated using quadratic interpolattion;
To step S1 treated signal, the presence or absence of carrier wave is judged by Threshold detection first, and rough estimate peak value Position;With xnUnmanned plane baseband signal is represented, if unmanned plane signal sequence to be estimated is
Wherein, A0WithThe respectively amplitude and first phase of carrier signal, if A0=1,f0And fsRespectively wait estimate The carrier frequency and data sampling rate of signal, N are number of sampling points, and r (n) is white Gaussian noise, if its variance is δr
Firstly, carrying out FFT transform to signal, it is as follows to obtain its frequency spectrum:
X (k)=FFT { x (n) }, k=0,1 ... N-1 (2)
Then, the ratio for calculating each frequency point signal strength and its neighborhood frequency point signal strength mean value is as follows:
Wherein
M is left and right adjoint point number, if above-mentioned ratio qjMore than setting threshold T, then be determined with carrier signal presence, then with This value of frequency point and its adjacent time big value of frequency point are reference carrier frequency;If this frequency point and its adjacent time big frequency point respective frequencies difference For f1 and f2, if f1 < f2, then have | f1-f2|=fs/N;[f1, f2] is signal(-) carrier frequency range to be estimated;
S2b, ultra-narrow low-pass filtering is carried out, specific treatment process is as follows:
Obtaining frequency range by FFT operation is [f1, f2], and the local oscillator complex signal frequency for calculating multiple modulation frequency displacement is
It is obtained shown in frequency shift signal such as formula (6) through multiple modulation, primary frequency feMoved zero-frequency, then at this time former frequency point f1, F2 is located at zero-frequency two sides and symmetrical about zero-frequency;
Low-pass filtering is carried out to signal using integral operation, in combination with data pick-up, sums to data sectional, takes each section Data mean value composition extracts signal sequence;If extraction yield is D, the mode of extraction is to take the mean value of D data point as extraction knot Fruit by way of being segmented summation and taking mean value while realizing low-pass filtering and data pick-up;Signal sampling rate becomes f ' after extractions =fs/ D, filtering extraction result are as follows:
Finally, envelope detection and pulse detection
After obtaining signal(-) carrier frequency, corresponding signal extraction is come out by envelope detection method, to realize pulse detection And parameter measurement;The step of envelope detection, is similar with front, i.e., signal is first downconverted to zero intermediate frequency, then passes through segmentation summation It takes the mode of mean value while realizing low-pass filtering and data pick-up, realize that narrow-band filtering and data volume reduce;I.e. by f in formula (5)e It changes into
Then repetitive (6)) and the step of formula (7), if obtained result is expressed as x2(n), envelope is taken to it, i.e., to x2 (n) modulus obtains x3(n)=| x2(n)|;
S2c, pulse detection is carried out to envelope waveform and measures related parameter values
Pulse detection is mainly that the detection of rising edge and failing edge is carried out to envelope waveform, using dynamic threshold detection method, Using the intersection point of envelope and thresholding as rising edge or failing edge;Envelope signal is first divided into M segment signal, Mei Yiduan by dynamic threshold method Signal is x4m(n), signal is divided into multiple uniform windows by m=1 ..., M, seek the part of signal in each window most Big value xmax(m)=max (x4m(n)), then generating corresponding thresholding is gm=α xmax(m);
It is respectively n that the pulse up and down obtained in one-time detection time width, which are located at, along positionupAnd ndown, rise in pulsewidth Edge and failing edge obtain up and down along position by the way of linear fit interpolation;Specific algorithm is as follows:
Wherein, tup、tdownRespectively rising edge, failing edge correspond to time, BdFor the bandwidth after filtering extraction;Need root The case where according to pulsewidth negative value, carries out displacement adjustment to rising edge or failing edge;Judge KupA rising edge and KdownA decline After, the time difference between adjacent rising edges or failing edge is the pulse repetition period, i.e.,
PRI=tup(k+1)-tup(k), k=1,2 ..., K (10)
The failing edge time of pairing is subtracted into rising time, pulsewidth can be obtained, i.e.,
PW=tdown(k)-tup(k), k=1,2 ..., K (11)
In addition, the rising edge of each pulse can be used as its arrival time;
It, can be according to measured signal(-) carrier frequency CF and pulse repetition period after estimating pulse relevant parameter PRI and pulsewidth PW parameter, are identified and are extracted to unmanned plane signal.
Further, in above-mentioned step S4, accurate the step of estimating frequency hopping rate, is as follows:
1) the Short Time Fourier Transform STFT for receiving signal x (t) is calculatedx
First assume hop rate estimationIt is to receive the time-frequency representation formula of signal Short Time Fourier Transform known to priori:
Wherein, h (τ-t) is window function, is 1 in (τ-t) ∈ [- Δ/2, Δ/2], other are 0, are set asWith Ensure at most there is a frequency hopping in the time zone of window function covering;
2) STFT is extractedxTime-frequency crestal line fx(t);
Secondary treatment is carried out for the time-frequency representation to FH signal, needs to extract its time-frequency crestal line, such as formula (18):
3) f is calculatedx(t) wavelet transformation W (a, t);
Using wavelet transformation
Wherein, a is scale parameter;It is Haar small echo;
Set the width Delta that scale a is equal to time frequency window, it is ensured that at most have one within the scope of the cover time of wavelet function Frequency hopping;
4) the amplitude sequence Abs [W (a, t)] of W (a, t) is calculated;
The amplitude sequence Abs [W (a, t)] of W (a, t) is calculated, which is a pseudo-random sequence, i.e. frequency hopping Occur in nTH(the n ∈ Z) moment;
5) the Fourier transformation FFT of Abs [W (a, t)] is calculated;
6) interval of detection spectrum spike, the interval are the estimation of corresponding hop rate: the wavelet transformation amplitude sequence of time-frequency crestal line There are discrete spectral lines at hop rate position, can be achieved with the accurate estimation of hop rate by detecting the discrete spectral line.
A kind of unmanned plane uplink and downlink communication link identification device based on signal characteristic difference comprising for receiving electromagnetism Wave signal simultaneously amplifies signal, filters, being mixed and the broadband receiver of IF process;For at through broadband receiver Signal after reason carries out carrier frequency CF estimation, pulse detection carried out in conjunction with narrow-band filtering, the pulse repetition period PRI of signal and The unmanned plane signal detection of pulse width PW parameter measurement prejudges module;For the pulse repetition period PRI to unmanned plane signal And the unmanned plane signal identification module that pulse width PW parameter is identified;For the bandwidth and frequency hopping rate progress to unmanned plane Precise measurement and the uplink and downlink signals of estimation detect judgment module.
Due to the adoption of the technical scheme as described above, the present invention has the following advantages:
The unmanned plane uplink and downlink communication link recognition based on signal characteristic difference, uses radio detection technology To identify unmanned plane uplink, downstream communications link under the conditions of non-cooperating, radio detection technology non-radiating electromagnetic signal, signal biography Defeated wave-path is the half of radar detection, therefore possesses farther detection range;Radio detection does not emit electromagnetic wave actively, hidden Covering property is strong, do not generate interference, can round-the-clock, all weather operations, while it is lower compared to active detection method cost;Firstly, According to unmanned plane signal, there are certain inherent features in time domain and frequency domain, have differences judgement with the other interference signals of same frequency range Unmanned plane signal whether there is, and then, further according to the feature difference of unmanned plane uplink and downlink link signal, identify uplink and downlink Link, and the frequency range of uplink and downlink signal is measured, and then advantage power can be concentrated targetedly on unmanned plane Downlink signal carries out compacting interference.Small drone refers to slipping into military area, test place, command post and security areas Wait military-civil unmanned vehicle, the wireless signal radiation source etc. of " black to fly " around critical facilities within the scope of 2km, and usual situation Lower empty weight is no more than 15kg, and take-off weight is no more than 25kg.
Detailed description of the invention
Fig. 1 is the flow chart of the unmanned plane uplink and downlink communication link recognition the present invention is based on signal characteristic difference;
Fig. 2 is unmanned plane uplink signal waveform and spectrogram;
Fig. 3 is unmanned plane downlink signal waveform and spectrogram;
Fig. 4 is unmanned plane uplink signal time frequency distribution map;
Fig. 5 is unmanned plane downlink signal time frequency distribution map.
Specific embodiment
Technical solution of the present invention is described in further detail with reference to the accompanying drawings and examples.
As shown in Figure 1, a kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference comprising with Lower specific steps:
S1, broadband receiver antenna receive electromagnetic wave signal, amplify, filter, being mixed and IF process;So as to Detection anticipation is carried out to signal in next step;
S2, carrier frequency CF measurement is carried out to step S1 treated signal, carries out pulse detection in conjunction with narrow-band filtering, surveys Measure signal pulse repetition period PRI and pulse width PW parameter, with carry out unmanned plane signal whether there is or not anticipation;Including having as follows Body step:
S2a, carrier frequency CF estimation
Rough estimate first is carried out to signal carrier frequency using spectrum peak detection calculations, frequency range is determined, then recycles (ZFFT is also referred to as Zoom-FFT to ZFFT, the Fast Fourier Transform referred to as refined, also known as choosing band Fast Fourier Transform, is one Kind complex electromagnetic parameters method) frequency spectrum is further refined, substantially reduce data volume while improving frequency resolution, most " fence effect " bring frequency deviation is eliminated using quadratic interpolattion afterwards;
To step S1 treated signal, the presence or absence of carrier wave is judged by Threshold detection first, and rough estimate peak value Position;With xnUnmanned plane baseband signal is represented, if unmanned plane signal sequence to be estimated is
Wherein, A0WithThe respectively amplitude and first phase of carrier signal, to simplify operation, if A0=1,f0And fs The carrier frequency and data sampling rate of signal respectively to be estimated, N are number of sampling points, and r (n) is white Gaussian noise, if its variance For δr
Firstly, carrying out FFT transform to signal, it is as follows to obtain its frequency spectrum:
X (k)=FFT { x (n) }, k=0,1 ... N-1 (2)
Then, the ratio for calculating each frequency point signal strength and its neighborhood frequency point signal strength mean value is as follows:
Wherein
M generally takes the frequency range of left and right 50kHz for left and right adjoint point number to be divided into reference between practical side frequency;On if State ratio qjMore than setting threshold T, T is 0.28~0.3, then is determined with carrier signal presence, then with this value of frequency point and its Adjacent time big value of frequency point is reference carrier frequency;If this frequency point and its adjacent time big frequency point respective frequencies are respectively f1 and f2, if F1 < f2, then have | f1-f2|=fs/N;Due to exist " fence effect ", the actual frequency of signal between f1 and f2, i.e., [f1, F2] it is signal(-) carrier frequency range to be estimated;
S2b, ultra-narrow low-pass filtering is carried out, specific treatment process is as follows:
Obtaining frequency range by FFT operation is [f1, f2], and the local oscillator complex signal frequency for calculating multiple modulation frequency displacement is
It is obtained shown in frequency shift signal such as formula (6) through multiple modulation, primary frequency feMoved zero-frequency, then at this time former frequency point f1, F2 is located at zero-frequency two sides and symmetrical about zero-frequency;
On the basis of not influencing filter effect, in order to simplify operation, improve efficiency, signal is carried out using integral operation Low-pass filtering sums to data sectional in combination with data pick-up, and each segment data mean value composition is taken to extract signal sequence;If taking out Taking rate is D, and the mode of extraction is to take the mean value of D data point as extraction as a result, same by way of being segmented summation and taking mean value Shi Shixian low-pass filtering and data pick-up improve processing speed;Signal sampling rate becomes f ' after extractions=fs/ D, filtering extraction knot Fruit is as follows:
Finally, envelope detection and pulse detection
After obtaining signal(-) carrier frequency, corresponding signal extraction is come out by envelope detection method, to realize pulse detection And parameter measurement;The step of envelope detection, is similar with front, i.e., signal is first downconverted to zero intermediate frequency, then passes through segmentation summation It takes the mode of mean value while realizing low-pass filtering and data pick-up, realize that narrow-band filtering and data volume reduce;I.e. by f in formula (5)e It changes into
Then repetitive (6)) and the step of formula (7), if obtained result is expressed as x2(n), envelope is taken to it, i.e., to x2 (n) modulus obtains x3(n)=| x2(n)|;
S2c, pulse detection is carried out to envelope waveform and measures related parameter values
Pulse detection is mainly that the detection of rising edge and failing edge is carried out to envelope waveform, using dynamic threshold detection method, Using the intersection point of envelope and thresholding as rising edge or failing edge;Envelope signal is first divided into M segment signal, Mei Yiduan by dynamic threshold method Signal is x4m(n), signal is divided into multiple uniform windows by m=1 ..., M, seek the part of signal in each window most Big value xmax(m)=max (x4m(n)), then generating corresponding thresholding is gm=α xmax(m), α=0.5 is generally taken;
It is respectively n that the pulse up and down obtained in one-time detection time width, which are located at, along positionupAnd ndown, in order to further The measurement accuracy for improving pulsewidth is more accurately obtained in pulsewidth rising edge and failing edge by the way of linear fit interpolation It rises, failing edge position;Specific algorithm is as follows:
Wherein, tup、tdownRespectively rising edge, failing edge correspond to time, BdFor the bandwidth after filtering extraction;Actually may be used Can rising edge and failing edge pairing it is wrong, leading to pulsewidth is negative value, it is therefore desirable to the case where according to pulsewidth negative value, to rising edge or Person's failing edge carries out displacement adjustment;Judge KupA rising edge and KdownAfter a failing edge, between adjacent rising edges or failing edge Time difference be the pulse repetition period, i.e.,
PRI=tup(k+1)-tup(k), k=1,2 ..., K (10)
The failing edge time of pairing is subtracted into rising time, pulsewidth can be obtained, i.e.,
PW=tdown(k)-tup(k), k=1,2 ..., K (11)
In addition, the rising edge of each pulse can be used as its arrival time;
It, can be according to measured signal(-) carrier frequency CF and pulse repetition period after estimating pulse relevant parameter PRI and pulsewidth PW parameter, are identified and are extracted to unmanned plane signal.
S3, determine whether unmanned plane signal
According to carrier frequency, pulse repetition period and width parameter that the measurement of step S1, S2 obtains, nobody can be judged whether there is Machine signal;If judging, no unmanned plane signal is returned to step S2 and continues to test, if tentatively determining whether man-machine signal, with regard to carrying out Uplink and downlink signals in step S4 detect judgement;
It is 14ms if there is pulse repetition period PRI, pulse width PW is 1ms, 2ms or 10ms;Spectral centroid CF is The signal of 2.4065GHz or 5.8GHz exists, it will be able to tentatively judge that unmanned plane signal exists;
S4, uplink and downlink signals detection judgement
After identifying unmanned plane signal according to step S3, need to carry out the bandwidth and frequency hopping rate of the unmanned plane signal Precise measurement and estimation;
Firstly, accurately estimating the bandwidth of signal, perform the following operation:
1) power spectrum is sought
The data x (n) (n=0,1, L, N-1) for being N by length, is divided into L sections, every section has M data, wherein L is not weigh The length of each data segment, L=N/M, the i-th segment data x when foldedi(n)=x (n+iM-M), is then added to window function ω (n) often On a data segment, each section of cyclic graph, i-th section of cyclic graph are found out
In formula, U is normalization factor,
Will be approximate irrelevant between each section of cyclic graph, last power Spectral Estimation is
2) wavelet decomposition obtains smooth power spectrum
Although there are many power spectrum smoothing that estimates above, but still include some high fdrequency components, wavelet transformation can be with The high frequency detail in signal is extracted, so carrying out wavelet decomposition for the detail section of signal and rough portion to the power spectrum estimated Separation, then extract the coefficient of rough part and power spectrum is reconstructed with this coefficient, so that it may eliminate the high frequency in power spectrum waveform Details obtains smooth waveform;
3) the mobile covariance of smooth power spectrum is calculated
After obtaining smooth power spectrum, the starting point and ending point of passband need to be extracted to estimate bandwidth, calculate mobile association side Difference:
R (k, k+1)=cov (d (k), d (k+1))
(k=1,2, L, N-1) (15)
In formula: d (k) is the reconstruction coefficients that rough part is extracted after wavelet decomposition, and the formula is utilized to traverse all data Point, what is found out is the mobile covariance between arbitrary neighborhood point;When every p, value can just change, so using Following formula calculates
In formula, k=1,1+p, L, N-p, the selection of p value, as p=1, are equivalent to by determining the characteristics of actual data point Formula;
4) estimate bandwidth
The position a and b where maximum 2 mobile covariance values are extracted, and is regarded as the starting point of bandwidth respectively And cut off, by difference | b-a-1 | as estimation bandwidth;Cycle-index is set, these steps is repeated and finds out estimation bandwidth Assembly average, as the precise bandwidth estimated value of unmanned plane signal;
Then, accurately estimate frequency hopping rate
The estimation of frequency hopping rate needs to analyze signal using Short Time Fourier Transform, obtains its time-frequency representation, recycles small echo The side information of the time-frequency representation is extracted in transformation, and then estimates hop rate using spectrum analysis;Detailed step is as follows:
1) the Short Time Fourier Transform STFT for receiving signal x (t) is calculatedx
First assume hop rate estimationIt is to receive the time-frequency representation formula of signal Short Time Fourier Transform known to priori:
Wherein, h (τ-t) is window function, is 1 in (τ-t) ∈ [- Δ/2, Δ/2], other are 0, are set asWith Ensure at most there is a frequency hopping in the time zone of window function covering;
2) STFT is extractedxTime-frequency crestal line fx(t);
Secondary treatment is carried out for the time-frequency representation to FH signal, needs to extract its time-frequency crestal line, such as formula (18):
3) f is calculatedx(t) wavelet transformation W (a, t);
In order to consider the identification of porch, using wavelet transformation,
Wherein, a is scale parameter;It is Haar small echo;
Set the width Delta that scale a is equal to time frequency window, it is ensured that at most have one within the scope of the cover time of wavelet function Frequency hopping;
4) the amplitude sequence Abs [W (a, t)] of W (a, t) is calculated;
The amplitude sequence Abs [W (a, t)] of W (a, t) is calculated, which is a pseudo-random sequence, i.e. frequency hopping Occur in nTH(the n ∈ Z) moment, and in nTH(the n ∈ Z) moment is different to establish a capital frequency hopping, the generation of jump whether be with Machine.
5) the Fourier transformation FFT of Abs [W (a, t)] is calculated;
6) interval of detection spectrum spike, the interval are the estimation of corresponding hop rate, the wavelet transformation amplitude sequence of time-frequency crestal line There are discrete spectral lines at hop rate position, and the accurate estimation of frequency hopping rate is achieved that by detecting the discrete spectral line.
The uplink signal and downlink letter of unmanned plane are distinguished according to signal bandwidth BW, the two characteristic quantities of frequency hopping rate HP Number: uplink signal is uplink signal, and general bandwidth is small, is 1.2MHz~2.8MHz, pulse width, and frequency hopping, frequency hopping rate Comparatively fast;Downlink signal band is roomy, is 9.8MHz~10.2MHz, and pulsewidth is wider, and frequency is stablized in bandwidth, and frequency hopping rate is slower.
Fig. 2 is remote signal (uplink signal) waveform and spectrogram of one of unmanned aerial vehicle example, is based on according to the present invention The unmanned plane uplink and downlink communication link recognition of signal characteristic difference can extract the recurrent pulse characteristic of remote signal, warp Estimation, the pulse repetition period is about 14ms, and pulse width respectively may be about 1ms and 2.17ms.
Fig. 3 is figure communication number (downlink signal) waveform and spectrogram of one of unmanned aerial vehicle example, is based on according to the present invention The unmanned plane uplink and downlink communication link recognition of signal characteristic difference, can extract the recurrent pulse characteristic of figure communication number, nothing Man-machine downlink signal also has apparent periodic pulse signal characteristic, and the pulse repetition period is about 14ms, and pulse width is about 10ms, spectral centroid are about 2.4065GHz.
Fig. 4, Fig. 5 are the time frequency distribution map of unmanned plane uplink signal and downlink signal respectively, are based on signal according to the present invention The unmanned plane uplink and downlink communication link recognition of feature difference, available uplink signal time-frequency figure, it is known that its center frequency point It is constantly jumping, and frequency hopping point is in rule variation;Unmanned plane uplink and downlink communication chain based on signal characteristic difference according to the present invention Road recognition methods, the frequency hopping rate of uplink signal are 2.14, and the frequency hopping rate of downlink signal is 0.71, the signal band of uplink Width is 1.2MHz, and the signal bandwidth of downlink is 10MHz.Therefore, according to the bandwidth and frequency hopping rate of signal to be estimated, flexibly Design a decision threshold, so that it may good differentiation and identification uplink signal and downlink signal.
Invention additionally discloses a kind of unmanned plane uplink and downlink communication link identification device based on signal characteristic difference comprising For receiving electromagnetic wave signal and amplifying, filter, be mixed and the broadband receiver of IF process to signal;For to warp Broadband receiver treated signal carries out carrier frequency CF estimation, carries out the pulse of pulse detection, signal in conjunction with narrow-band filtering Repetition period PRI and the unmanned plane signal detection of pulse width PW parameter measurement prejudge module;For the arteries and veins to unmanned plane signal Rush repetition period PRI and unmanned plane signal identification module that pulse width PW parameter is identified;For the bandwidth to unmanned plane The uplink and downlink signals for carrying out precise measurement and estimation with frequency hopping rate detect judgment module.
The foregoing is merely presently preferred embodiments of the present invention, rather than limitation of the present invention, is not departing from essence of the invention In the case where mind and range, equivalent changes and modifications made according to the patent scope of the present invention should all belong to of the invention special Within sharp protection scope.

Claims (4)

1. a kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference, it is characterized in that: it includes following Step:
S1, broadband receiver antenna receive electromagnetic wave signal, amplify, filter, being mixed and IF process;
S2, carrier frequency CF measurement is carried out to step S1 treated signal, carries out pulse detection, measurement letter in conjunction with narrow-band filtering Number pulse repetition period PRI and pulse width PW parameter, with carry out unmanned plane signal whether there is or not anticipation;
S3, determine whether unmanned plane signal
Carrier frequency CF, the pulse repetition period PRI and pulse width PW parameter obtained according to step S2 measurement, can judge Whether unmanned plane signal is had;If judging, no unmanned plane signal is returned to step S2 and continues to test, if tentatively determining whether man-machine Signal detects judgement with regard to the uplink and downlink signals carried out in step S4;
S4, uplink and downlink signals detection judgement is carried out to unmanned plane signal
After identifying unmanned plane signal, bandwidth and frequency hopping rate to the unmanned plane signal carry out precise measurement and estimation;
Firstly, accurately estimating the bandwidth of signal, perform the following operation:
1) power spectrum is sought
The data x (n) (n=0,1, L, N-1) for being N by length, is divided into L sections, and every section has M data, and wherein L is every when not being overlapped The length of a data segment, L=N/M, the i-th segment data xi(n) window function ω (n), is then added to each data by=x (n+iM-M) Duan Shang finds out each section of cyclic graph, i-th section of cyclic graph
In formula, U is normalization factor,
Will be approximate irrelevant between each section of cyclic graph, last power Spectral Estimation is
2) wavelet decomposition obtains smooth power spectrum
Wavelet decomposition is carried out to the power spectrum that estimates, is partially separated by the detail section of signal and roughly, then extract rough portion Point coefficient and with this coefficient reconstruct power spectrum, eliminate power spectrum waveform in high frequency detail, obtain smooth waveform;
3) the mobile covariance of smooth power spectrum is calculated
After obtaining smooth power spectrum, the starting point and ending point of passband is extracted to estimate bandwidth;Calculate mobile covariance:
R (k, k+1)=cov (d (k), d (k+1))
(k=1,2, L, N-1) (15)
In formula: d (k) is the reconstruction coefficients that rough part is extracted after wavelet decomposition, traverses all data points using the formula, asks Out be mobile covariance between arbitrary neighborhood point;When every p, value can just change, so using following formula It calculates
In formula, k=1,1+p, L, N-p, the selection of p value, as p=1, are equivalent to formula by determining the characteristics of actual data point;
4) estimate bandwidth
The position a and b where maximum 2 mobile covariance values are extracted, and is regarded as the starting point of bandwidth respectively and cuts Stop, by difference | b-a-1 | as estimation bandwidth;
Cycle-index is set, is repeated the above steps 1)~4 to the L section for the data that length is N), L bandwidth is calculated, then be averaging, Find out the assembly average of estimation bandwidth, as the precise bandwidth estimated value of unmanned plane signal;
Then, accurately estimate frequency hopping rate, perform the following operation: analyzing signal using Short Time Fourier Transform, obtain its frequency schedule Show, recycles wavelet transformation to extract the side information of the time-frequency representation, and then estimate hop rate using spectrum analysis.
2. the unmanned plane uplink and downlink communication link recognition according to claim 1 based on signal characteristic difference, special Sign is: in its step S2, the method for unmanned plane signal detection anticipation includes the following steps:
S2a, carrier frequency CF estimation
Rough estimate first is carried out to signal carrier frequency using spectrum peak detection calculations, determines frequency range;Then ZFFT pairs is recycled Frequency spectrum further refines;Finally " fence effect " bring frequency deviation is eliminated using quadratic interpolattion;
To step S1 treated signal, the presence or absence of carrier wave, and the position of rough estimate peak value are judged by Threshold detection first; With xnUnmanned plane baseband signal is represented, if unmanned plane signal sequence to be estimated is
Wherein, A0WithThe respectively amplitude and first phase of carrier signal, if A0=1,f0And fsSignal respectively to be estimated Carrier frequency and data sampling rate, N are number of sampling points, and r (n) is white Gaussian noise, if its variance is δr
Firstly, carrying out FFT transform to signal, it is as follows to obtain its frequency spectrum:
X (k)=FFT { x (n) }, k=0,1 ... N-1 (2)
Then, the ratio for calculating each frequency point signal strength and its neighborhood frequency point signal strength mean value is as follows:
Wherein
M is left and right adjoint point number, if above-mentioned ratio qjMore than setting threshold T, then it is determined with carrier signal presence, then with this frequency Point value and its adjacent time big value of frequency point are reference carrier frequency;If this frequency point and its adjacent time big frequency point respective frequencies are respectively f1 And f2, if f1 < f2, then have | f1-f2|=fs/N;[f1, f2] is signal(-) carrier frequency range to be estimated;
S2b, ultra-narrow low-pass filtering is carried out, specific treatment process is as follows:
Obtaining frequency range by FFT operation is [f1, f2], and the local oscillator complex signal frequency for calculating multiple modulation frequency displacement is
It is obtained shown in frequency shift signal such as formula (6) through multiple modulation, primary frequency feZero-frequency is moved, then original frequency point f1, f2 is located at this time Zero-frequency two sides and symmetrical about zero-frequency;
Low-pass filtering is carried out to signal using integral operation, in combination with data pick-up, sums to data sectional, takes each segment data Mean value composition extracts signal sequence;If extraction yield is D, the mode of extraction is to take the mean value of D data point as extraction as a result, logical Segmentation summation is crossed to take the mode of mean value while realizing low-pass filtering and data pick-up;Signal sampling rate becomes f after extractions'=fs/ D, filtering extraction result are as follows:
Finally, envelope detection and pulse detection
After obtaining signal(-) carrier frequency, corresponding signal extraction is come out by envelope detection method, to realize pulse detection and ginseng Number measurement;The step of envelope detection, is similar with front, i.e., signal is first downconverted to zero intermediate frequency, is then taken by being segmented summation The mode of value realizes low-pass filtering and data pick-up simultaneously, realizes that narrow-band filtering and data volume reduce;I.e. by f in formula (5)eIt changes into
Then repetitive (6)) and the step of formula (7), if obtained result is expressed as x2(n), envelope is taken to it, i.e., to x2(n) it takes Mould obtains x3(n)=| x2(n)|;
S2c, pulse detection is carried out to envelope waveform and measures related parameter values
Pulse detection is mainly that the detection for carrying out rising edge and failing edge to envelope waveform will be wrapped using dynamic threshold detection method The intersection point of network and thresholding is as rising edge or failing edge;Envelope signal is first divided into M segment signal, each segment signal by dynamic threshold method For x4m(n), signal is divided into multiple uniform windows, seeks the local maximum of signal in each window by m=1 ..., M xmax(m)=max (x4m(n)), then generating corresponding thresholding is gm=α xmax(m);
It is respectively n that the pulse up and down obtained in one-time detection time width, which are located at, along positionupAnd ndown, in pulsewidth rising edge and Failing edge obtains up and down along position by the way of linear fit interpolation;Specific algorithm is as follows:
Wherein, tup、tdownRespectively rising edge, failing edge correspond to time, BdFor the bandwidth after filtering extraction;It needs according to arteries and veins The case where wide negative value, carries out displacement adjustment to rising edge or failing edge;Judge KupA rising edge and KdownAfter a failing edge, Time difference between adjacent rising edges or failing edge is the pulse repetition period, i.e.,
PRI=tup(k+1)-tup(k), k=1,2 ..., K (10)
The failing edge time of pairing is subtracted into rising time, pulsewidth can be obtained, i.e.,
PW=tdown(k)-tup(k), k=1,2 ..., K (11)
In addition, the rising edge of each pulse can be used as its arrival time;
After estimating pulse relevant parameter, can according to measured signal(-) carrier frequency CF and pulse repetition period PRI and Pulsewidth PW parameter is identified and is extracted to unmanned plane signal.
3. the unmanned plane uplink and downlink communication link recognition according to claim 1 based on signal characteristic difference, special Sign is: in its step S4, accurate the step of estimating frequency hopping rate, is as follows:
1) the Short Time Fourier Transform STFT for receiving signal x (t) is calculatedx
First assume hop rate estimationIt is to receive the time-frequency representation formula of signal Short Time Fourier Transform known to priori:
Wherein, h (τ-t) is window function, is 1 in (τ-t) ∈ [- Δ/2, Δ/2], other are 0, are set asTo ensure At most there is a frequency hopping in the time zone of window function covering;
2) STFT is extractedxTime-frequency crestal line fx(t);
Secondary treatment is carried out for the time-frequency representation to FH signal, needs to extract its time-frequency crestal line, such as formula (18):
3) f is calculatedx(t) wavelet transformation W (a, t);
Using wavelet transformation
Wherein, a is scale parameter;It is Haar small echo;
Set the width Delta that scale a is equal to time frequency window, it is ensured that at most there is a frequency within the scope of the cover time of wavelet function Jump;
4) the amplitude sequence Abs [W (a, t)] of W (a, t) is calculated;
The amplitude sequence Abs [W (a, t)] of W (a, t) is calculated, which is a pseudo-random sequence, i.e., frequency hopping occurs In nTH(the n ∈ Z) moment;
5) the Fourier transformation FFT of Abs [W (a, t)] is calculated;
6) interval of detection spectrum spike, the interval are the estimation of corresponding hop rate: the wavelet transformation amplitude sequence of time-frequency crestal line is being jumped There are discrete spectral lines at variable Rate position, can be achieved with the accurate estimation of hop rate by detecting the discrete spectral line.
4. a kind of implement in claims 1 to 33 on the unmanned plane based on signal characteristic difference of any claim the method Downstream communications link identification device, it is characterized in that: it include for receive electromagnetic wave signal and signal amplified, filter, The broadband receiver of mixing and IF process;For to through broadband receiver treated signal carry out carrier frequency CF estimation, The unmanned plane letter of pulse detection, the pulse repetition period PRI of signal and pulse width PW parameter measurement is carried out in conjunction with narrow-band filtering Number detection anticipation module;For the pulse repetition period PRI and the nothing that is identified of pulse width PW parameter to unmanned plane signal Man-machine signal identification module;Uplink and downlink signals inspection for bandwidth and frequency hopping rate progress precise measurement and estimation to unmanned plane Survey judgment module.
CN201910808275.0A 2019-08-29 2019-08-29 Unmanned aerial vehicle uplink and downlink communication link identification method and device based on signal characteristic difference Active CN110519003B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910808275.0A CN110519003B (en) 2019-08-29 2019-08-29 Unmanned aerial vehicle uplink and downlink communication link identification method and device based on signal characteristic difference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910808275.0A CN110519003B (en) 2019-08-29 2019-08-29 Unmanned aerial vehicle uplink and downlink communication link identification method and device based on signal characteristic difference

Publications (2)

Publication Number Publication Date
CN110519003A true CN110519003A (en) 2019-11-29
CN110519003B CN110519003B (en) 2023-04-28

Family

ID=68627912

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910808275.0A Active CN110519003B (en) 2019-08-29 2019-08-29 Unmanned aerial vehicle uplink and downlink communication link identification method and device based on signal characteristic difference

Country Status (1)

Country Link
CN (1) CN110519003B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111224912A (en) * 2020-01-16 2020-06-02 西安理工大学 Method for determining time difference of arrival of air-ground link signal, interception station and storage medium
CN111600824A (en) * 2020-04-22 2020-08-28 中国人民解放军战略支援部队信息工程大学 Unmanned aerial vehicle graph signaling signal identification method and device
US10948568B1 (en) 2019-09-12 2021-03-16 Cypress Semiconductor Corporation Radar detection in a wireless LAN
CN113242114A (en) * 2021-05-08 2021-08-10 成都空间矩阵科技有限公司 System and method for detecting burst type narrow pulse signal
CN114189310A (en) * 2021-12-07 2022-03-15 中国人民解放军32802部队 Unmanned aerial vehicle measurement and control signal accurate interference method based on signal reconnaissance and prediction
CN114244450A (en) * 2021-12-03 2022-03-25 北京蓝玛星际科技有限公司 Signal identification method and device, electronic equipment and storage medium
CN114696922A (en) * 2022-02-22 2022-07-01 电子科技大学 Frequency hopping signal detection method suitable for unmanned aerial vehicle communication
CN115065580A (en) * 2022-07-28 2022-09-16 成都华日通讯技术股份有限公司 Link16 data chain identification and parameter estimation method under broadband
CN117786398A (en) * 2024-02-28 2024-03-29 杰创智能科技股份有限公司 Frequency hopping signal characteristic identification method, system, electronic equipment and storage medium
CN117786398B (en) * 2024-02-28 2024-05-31 杰创智能科技股份有限公司 Frequency hopping signal characteristic identification method, system, electronic equipment and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106549893A (en) * 2015-09-16 2017-03-29 中国航空工业第六八研究所 A kind of Unmanned Aerial Vehicle Data link frequency deviation detection method
CN108700668A (en) * 2017-06-29 2018-10-23 深圳市大疆创新科技有限公司 Detect method, the unmanned plane of the positioning device of unmanned plane
US20190051046A1 (en) * 2017-08-11 2019-02-14 Jing Jin Incident Site Investigation and Management Support System Based on Unmanned Aerial Vehicles
CN109490969A (en) * 2018-09-21 2019-03-19 浙江大学 One kind being based on the periodic unmanned plane detection method of unmanned plane downlink signal
WO2019099384A1 (en) * 2017-11-16 2019-05-23 Kyocera Corporation Uplink interference-based monitoring of downlink signals by unmanned aerial vehicle

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106549893A (en) * 2015-09-16 2017-03-29 中国航空工业第六八研究所 A kind of Unmanned Aerial Vehicle Data link frequency deviation detection method
CN108700668A (en) * 2017-06-29 2018-10-23 深圳市大疆创新科技有限公司 Detect method, the unmanned plane of the positioning device of unmanned plane
US20190051046A1 (en) * 2017-08-11 2019-02-14 Jing Jin Incident Site Investigation and Management Support System Based on Unmanned Aerial Vehicles
WO2019099384A1 (en) * 2017-11-16 2019-05-23 Kyocera Corporation Uplink interference-based monitoring of downlink signals by unmanned aerial vehicle
CN109490969A (en) * 2018-09-21 2019-03-19 浙江大学 One kind being based on the periodic unmanned plane detection method of unmanned plane downlink signal

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10948568B1 (en) 2019-09-12 2021-03-16 Cypress Semiconductor Corporation Radar detection in a wireless LAN
CN111224912B (en) * 2020-01-16 2020-09-29 西安理工大学 Method for determining time difference of arrival of air-ground link signal, interception station and storage medium
CN111224912A (en) * 2020-01-16 2020-06-02 西安理工大学 Method for determining time difference of arrival of air-ground link signal, interception station and storage medium
CN111600824A (en) * 2020-04-22 2020-08-28 中国人民解放军战略支援部队信息工程大学 Unmanned aerial vehicle graph signaling signal identification method and device
CN111600824B (en) * 2020-04-22 2023-04-07 中国人民解放军战略支援部队信息工程大学 Unmanned aerial vehicle graph signaling signal identification method and device
CN113242114B (en) * 2021-05-08 2023-01-31 成都空间矩阵科技有限公司 Burst type narrow pulse signal detection system and method
CN113242114A (en) * 2021-05-08 2021-08-10 成都空间矩阵科技有限公司 System and method for detecting burst type narrow pulse signal
CN114244450A (en) * 2021-12-03 2022-03-25 北京蓝玛星际科技有限公司 Signal identification method and device, electronic equipment and storage medium
CN114189310B (en) * 2021-12-07 2022-06-21 中国人民解放军32802部队 Unmanned aerial vehicle measurement and control signal accurate interference method based on signal reconnaissance and prediction
CN114189310A (en) * 2021-12-07 2022-03-15 中国人民解放军32802部队 Unmanned aerial vehicle measurement and control signal accurate interference method based on signal reconnaissance and prediction
CN114696922A (en) * 2022-02-22 2022-07-01 电子科技大学 Frequency hopping signal detection method suitable for unmanned aerial vehicle communication
CN115065580A (en) * 2022-07-28 2022-09-16 成都华日通讯技术股份有限公司 Link16 data chain identification and parameter estimation method under broadband
CN115065580B (en) * 2022-07-28 2024-04-02 成都华日通讯技术股份有限公司 Link16 data Link identification and parameter estimation method under broadband
CN117786398A (en) * 2024-02-28 2024-03-29 杰创智能科技股份有限公司 Frequency hopping signal characteristic identification method, system, electronic equipment and storage medium
CN117786398B (en) * 2024-02-28 2024-05-31 杰创智能科技股份有限公司 Frequency hopping signal characteristic identification method, system, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN110519003B (en) 2023-04-28

Similar Documents

Publication Publication Date Title
CN110519003A (en) A kind of unmanned plane uplink and downlink communication link recognition based on signal characteristic difference
CN109981192B (en) Frequency spectrum monitoring method of airspace blackout flying unmanned aerial vehicle
US9991908B2 (en) Blind source separation of signals having low signal-to-noise ratio
US20090075590A1 (en) Method and System for Estimating Time of Arrival of Signals Using Multiple Different Time Scales
CN112684251B (en) Target signal frequency domain detection method based on power spectrum template
CN102636775B (en) Wind profile radar echo spectrum reconfiguration method based on fuzzy logic recognition
CN111896926B (en) Low-altitude target detection method and system based on strong clutter suppression
CN108872962A (en) Laser radar weak signal extraction and decomposition method based on Fourier Transform of Fractional Order
CN109709512A (en) A kind of single station of unmanned plane detecting and unmanned plane detecting system
CN109614930A (en) A kind of unmanned plane spectrum detection method based on deep learning
CN110289918A (en) A kind of wireless signal detection method and system
CN111948618A (en) Forward scattering target detection method and system based on satellite external radiation source
Schneckenburger et al. Reflector localization for geometrical modeling the air–ground channel
CN112818876A (en) Electromagnetic signal extraction and processing method based on deep convolutional neural network
CN113176542B (en) Time domain self-adaptive main lobe self-defense interference resisting method and system
CN111046025A (en) Unmanned aerial vehicle signal detection method and device
CN105336340A (en) Wind noise rejection method and device for low altitude target acoustic detection system
CN112929141A (en) Unmanned aerial vehicle detection and identification method and system based on graph-borne signal matching
Djurovic et al. High-resolution time-frequency representations based on the local polynomial Fourier transform for over-the-horizon radars
Müller et al. Acoustic flow perception in cf-bats: extraction of parameters
CN102685049B (en) Fractional order channelized separation method for reaching two linear frequency modulation (LFM) signals at the same time
CN105911546B (en) A kind of recognition methods of sea clutter and device
Fang et al. DME pulse interference suppression based on NLS for GPS
CN103885044A (en) Method for suppressing clutter and noise of narrow-band radar echoes based on CLEAN algorithm
CN109347580B (en) Self-adaptive threshold signal detection method with known duty ratio

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