CN106680795A - Time-domain modulation domain parameter combination measure method of frequency agile signal - Google Patents

Time-domain modulation domain parameter combination measure method of frequency agile signal Download PDF

Info

Publication number
CN106680795A
CN106680795A CN201611246663.7A CN201611246663A CN106680795A CN 106680795 A CN106680795 A CN 106680795A CN 201611246663 A CN201611246663 A CN 201611246663A CN 106680795 A CN106680795 A CN 106680795A
Authority
CN
China
Prior art keywords
frequency
signal
pulse
value
agile
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.)
Pending
Application number
CN201611246663.7A
Other languages
Chinese (zh)
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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN201611246663.7A priority Critical patent/CN106680795A/en
Publication of CN106680795A publication Critical patent/CN106680795A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418—Theoretical aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a time-domain modulation domain parameter combination measurement method of frequency agile signal which mainly solves the problem low precision of signal parameter and large amount of calculation of frequency agile signal measured by the present technology measures. The implementation steps comprise the following: 1. the real part of the frequency signal is obtained by the signal analyzer; 2. a signal envelope of the frequency agile signal is calculated according to the discrete numerical value; 3. the signal envelope estimates the approximate amplitude of the frequency agile signal; 4. by the approximate amplitude to set the threshold of judgment, the number of pulses to detect the frequency of agile signal of the frequency is measured; 5. the pulse width of each pulse is calculated; 6. the time domain parameters are calculated by the number of pulses and the pulse width;7.a instantaneous frequency of the frequency agile signal is calculated based on the time domain parameters; 8. the instantaneous frequency calculated the modulation domain parameters of the frequency agile signal. The system has the advantages of small calculating amount and higher accuracy and can be used radar signal processing.

Description

The temporal modulation field parameter union measuring method of frequency agility signal
Technical field
The present invention relates to Radar Technology field, the temporal modulation field parameter joint survey of more particularly to a kind of frequency agility signal Amount method, can be used for Radar Signal Processing.
Background technology
Frequency agility signal, refers to one kind of carrier frequency random quick change in certain frequency band of radar signal adjacent pulse Signal type.The parameter of measurement frequency agile signal is one of emphasis of Radar Signal Processing engineering field.Traditional measurement side Only for time domain or modulation domain, time domain and modulation domain are not combined together carries out parameter measurement to method, and measurement frequency is prompt Varying signal modulation field parameter needs to carry out time frequency analysis, conventional Time-Frequency Analysis Method, such as Short Time Fourier Transform STFT and Wigner-Vile is distributed, and they have many defects when to frequency agility signal time frequency analysis.STFT exist window length and Window function selects difficult problem, and due to having used Fourier transformation, operand is larger, and can not solve temporal resolution and The contradiction of frequency resolution.Although Wigner-Vile distributions do not use window function, not there is a problem of STFT, to many points Amount signal carries out there is " cross term " interference during time-frequency distributions, and time-frequency distributions will thicken, inaccurate to Frequency Estimation, and It employs integral operation, and operand is larger.
The content of the invention
Present invention aims to the deficiency of above-mentioned prior art, proposes a kind of temporal modulation of frequency agility signal Field parameter union measuring method, to reduce operand and error, obtains accurate agile frequency modulation domain parameter, improves measurement essence Degree.
The technical thought for realizing the object of the invention is joint time domain, modulation domain measurement, using statistics pulse during time-domain analyses The method of signal envelope numerical value estimates decision threshold, the measurement of the time domain parameters such as pulse number, pulse width is realized, in modulation domain The instantaneous frequency profile of signal is obtained during analysis using Phase difference, the modulation field parameter such as accurate agile frequency is obtained, Its technical step includes as follows：
(1) discrete values of measurement frequency agile signal；
1a) obtain frequency agility signal real part I (t) to be measured with signal analyzer reading big with the numerical value of imaginary part Q (t) It is little；
1b) set sample rate fs, the frequency agility signal real part to be measured and imaginary part of sampling 1 second time span of selection, Obtain this discrete value I (n) of frequency agility signal real part and discrete value Q (n) of imaginary part；
(2) time domain parameter of frequency agility signal is obtained by discrete value：
(2a) according to the real of frequency agility signal, discrete value I (n) of imaginary part and Q (n), frequency agility signal envelope is calculated：
(2b) approximate top value A' of frequency agility signal is estimated according to frequency agility signal envelope y (n)topWith approximate bottom Value A'ba, obtaining approximate range measured value is：A'=A'top-A'ba；
(2c) with 10% approximate range measured value 10%A' as rising edge and the approximate detection threshold of trailing edge, detection The pulse number n of frequency agility signal is obtained, with the points and trailing edge at the approximate detection threshold place of each rising edge of a pulse Approximate detection threshold be located points be respectively divided by sample rate, obtain the approximate judgement of the rising edge of each corresponding pulse Moment t 'RiseWith the approximate judgement moment t ' of trailing edgeDrop, the pulse width for calculating pulse signal each pulse is：τ=t 'Rise-t′Drop；
(2d) in the pulse width of each pulse, with Density Distribution averaging method, the accurate top value of each pulse is calculated AtopWith floors Aba, and try to achieve each pulse Precise amplitudes value A and be：A=Atop-Aba；
(2e) according to the maximum V of each pulsemaxWith minima Vmin, calculate overshoot S of each pulseoverAnd undershoot Sunder：
Sover=Vmax-Atop, Sunder=Aba-Vmin,
(2f) according to accurate top value A of each pulse of frequency agility signaltop, floors AbaWith Precise amplitudes value A, calculate Reference line M in the pulse amplitude of frequency agility signalOnWith reference line M under pulse amplitudeUnder：
MOn=Aba+ 90% × A, MUnder=Aba+ 10% × A；
(2g) according to reference line M in the pulse amplitude of frequency agility signalOn, reference line M under pulse amplitudeUnderAnd pulse width τ, searches respectively in rising edge corresponding moment value t of reference line in pulse amplitude in pulse widthOn r, refer under pulse amplitude Corresponding moment value t of lineUnder rWith corresponding moment value t of reference line in pulse amplitude in trailing edgeOn f, reference line pair under pulse amplitude Moment value t answeredUnder f, and the moment t' in adjacent rising edge next time in pulse amplitude corresponding to reference lineUnder r, calculate frequency The rise time t of agile signal pulser, fall time tf, pulse period T, shut-in time toffWith dutycycle dt：
tr=tOn r-tUnder r, tf=tUnder f-tOn f, T=t'Under r-tUnder r, toff=T- τ, dt=τ/T；
(3) the modulation field parameter of frequency agility signal is obtained by discrete value：
(3a) according to the real of frequency agility signal, discrete value I (n) of imaginary part and Q (n), the wink of frequency agility signal is calculated Shi Xiangwei：
(3b) according to the instantaneous phase of frequency agility signalWhen being carried out to frequency agility signal using Phase difference Frequency analysis, obtain the instantaneous frequency of frequency agility signal：Wherein fsFor the sampling of signal Rate；
(3c) according to instantaneous frequency f of frequency agility signalcN (), calculates the meansigma methodss of frequency agility signal transient frequency Favg；
(3d) according to meansigma methodss F of frequency agility signal transient frequencyavgWith pulse number n, statistic frequency agile signal wink When number N of the frequency near meansigma methodss, if meeting 90%n≤N≤n, the agile mode for judging the frequency agility signal is Agile in arteries and veins, is otherwise judged as agile between arteries and veins；
(3e) according to frequency agility signal reality, discrete value I (n) of imaginary part and Q (n), respectively to the point in each pulse Real part, imaginary part does fast Fourier transform and obtains corresponding real part I (ω) and imaginary part Q (ω), and calculate its instantaneous power：
(3f) according to instantaneous power spec (ω) of frequency agility signal, agile mode between agile mode or arteries and veins in calculating arteries and veins Agile number C and agile frequency freq (i), i=1,2 ..., C.
The present invention has the advantage that compared with prior art：
1st, the present invention estimates frequency agility during the approximate range of pulse is estimated using Density Distribution averaging method Approximate top value A' of signaltopWith approximate floors A'ba, precision is higher compared with existing Density Distribution mode method, works as frequency agility When the sampling number of signal is more, pulse envelope maximum ymaxWith pulse envelope minima yminBetween divide equally interval number get over Approximate floors A' that are many, being tried to achievetopWith approximate top value A'baIt is more accurate, more preferable required precision can be met.
2nd, the present invention carries out adopting Phase difference during time frequency analysis, and relative to STFT methods, operand is greatly reduced, and Without window function, also do not limited by indeterminacy criterion；Relative to Wigner-Vile distribution methods, there is no the interference of " cross term " Problem, it is little with measurement parameter operand relative to additive method, the advantages of error is little, it is suitable for engineering practice.
Below in conjunction with the accompanying drawings, the present invention is described in further detail.
Description of the drawings
Fig. 1 is the flowchart of the present invention；
Fig. 2 is the time-frequency curve chart for carrying out time frequency analysis to frequency agility signal with existing STFT methods；
Fig. 3 is the time-frequency curve chart for carrying out time frequency analysis to frequency agility signal with existing Wigner-Vile methods；
Fig. 4 is the instantaneous phase curve and time-frequency curve chart that time frequency analysis are carried out with the inventive method.
Specific embodiment
In radar is received, the time domain of radar return, modulation field parameter are equipped with many useful informations, and radar return Parameter is unknown, can not pass through radar receiving antenna direct access.For different types of radar signal, its parameter measurement Method is also not quite similar, and frequency agility signal is a kind of conventional transmission signal of radar, and measurement frequency agile signal parameter is thunder One of emphasis up to signal processing engineering field.Frequency agility signal makes parameter measurement particularly due to the complicated characteristic of intra-pulse modulation Difficulty, operand is huge, and Jing often occurs measuring inaccurate problem.In order to accurately quick measurement frequency agile signal Time domain, modulation field parameter, it is necessary to take the measuring method that a kind of operand is little and precision is higher.The present invention is exactly by signal The numerical values recited of analyser reading frequency agile signal, and then using joint time domain, the thought of modulation domain, it is average with Density Distribution The design parameter of the more accurate measurement frequency agile signal in detail of method and Phase difference.
With reference to Fig. 1, the implementation steps of the present invention are as follows：
Step 1：The discrete values of measurement frequency agile signal.
1a) the radar return that radar passes through reception antenna receives frequency agile signal form, and the frequency for receiving is prompt Real part I (t) and imaginary part Q (t) of varying signal is saved in .dat files with data mode；
1b) by .dat files passes to signal analyzer, reading obtains real part I (t) and imaginary part Q of frequency agility signal The numerical values recited of (t)；
1c) set sample rate fs, the frequency agility signal to be measured of 1 second time span is extracted, discrete adopting is carried out to it Sample, obtains discrete value I (n) of real part and discrete value Q (n) of imaginary part of the frequency agility signal of this 1 second length.
Step 2：Calculate the time domain parameter of frequency agility signal.
2a) according to the real of frequency agility signal, discrete value I (n) of imaginary part and Q (n), frequency agility signal envelope is calculated：
2b) estimate approximate top value A' of frequency agility signaltopWith approximate floors A'ba, and approximate range A'；
2b1) search the maximum y for obtaining signal envelope y (n)maxWith minima ymin, wherein n=1,2 ..., N, N are frequency Total points obtained by rate agile signal sampling；
Signal envelope y (n) 2b2) is divided into from small to large 100 intervals, wherein k-th interval E (k) span is most Little value is EminK (), maximum is Emax(k)：
2b3) statistical signal envelope y (n) falls the number in each interval E (k), is designated as count value C (k) respectively, k=1, 2…100；
The position that the maximum of front 50 count values is located 2b4) is searched in count value C (k), k is designated as1, 50 after lookup The position that the maximum of individual count value is located, is designated as k2；
2b5) computation interval E (k1) meansigma methodss, as approximate top value A' of frequency agility signaltop, computation interval E (k2) Meansigma methodss, the as approximate floors A' of frequency agility signalba；
2b6) according to the approximate floors A' of frequency agility signalbaWith approximate top value A'top, obtain the near of frequency agility signal Like amplitude A '=A'top-A'ba；
2c) the pulse number m of frequency agility signal is detected, and record the approximate rising edge judgement moment t' of each pulseRise、 The judgement moment t' of approximate trailing edgeDrop；
2c1) pulse number of frequency agility signal is designated as into m, and it is 0 to initialize m, from the start bit of frequency agility signal Put and begin look for first approximate rising edge, lookup result is designated as into flag bit flag, if searching successfully, by flag bit flag 1 is put, and records this approximate rising edge and adjudicate moment t'Rise, otherwise, flag bit flag is set to 0；
After 2c2) finding approximate rising edge, approximate trailing edge is searched in continuation backward, if searching successfully, by flag bit Flag puts 1, and records this approximate trailing edge judgement moment t'Drop, otherwise, flag bit flag is set to 0；
After 2c3) finding approximate trailing edge, approximate rising edge is searched in continuation backward, if searching successfully, by flag bit Flag puts 1, records this approximate rising edge and adjudicates moment t'Rise, and make the pulse number m of frequency agility signal plus 1, and otherwise, will Flag bit flag sets to 0；
The currency of flag bit flag 2c4) is detected, if flag bit flag=1, return to step 2c2), if flag bit Flag=0, then terminate to search；
2d) calculate the width of each pulse of frequency agility signal：τ=t'Drop-t'Rise；
2e) in the pulse width of each pulse, with Density Distribution averaging method, it is calculated each pulse and accurately pushes up Value AtopWith floors Aba, and try to achieve each pulse Precise amplitudes value A and be：A=Atop-Aba；
2f) detect the maximum V of each pulse of frequency agility signalmaxWith minima Vmin, calculate the overshoot of each pulse Sover=Vmax-AtopWith undershoot Sunder=Aba-Vmin；
2g) calculate the rise time t of pulser, fall time tf, pulse period T, shut-in time toffAnd dutycycle dt：
2g1) reference line M in the pulse amplitude of setpoint frequency agile signalOnWith reference line M under pulse amplitudeUnder：
MOn=Aba+ 90% × A,
MUnder=Aba+ 10% × A；
2g2) search reference line M in pulse amplitude in pulse width τOnWith reference line M under pulse amplitudeUnderThe corresponding moment Value, the corresponding moment value includes：Moment t in rising edge in pulse amplitude corresponding to reference lineOn r, pulse width in rising edge Moment t corresponding to the lower reference line of degreeUnder r, moment t in trailing edge in pulse amplitude corresponding to reference lineOn f, pulse in trailing edge Moment t under amplitude corresponding to reference lineUnder f, and the moment in adjacent rising edge next time in pulse amplitude corresponding to reference line t'Under r；
2g3) according to 2g2) search the moment value for going out, calculate the following parameter of pulse：
Rise time tr：tr=tOn r-tUnder r,
Fall time tf：tf=tUnder f-tOn f,
Pulse period T：T=t'Under r-tUnder r,
Shut-in time toff：toff=T- τ,
Dutycycle dt：dt=τ/T.
Step 3：Calculate the modulation field parameter of frequency agility signal.
3a) according to the real of frequency agility signal, discrete value I (n) of imaginary part and Q (n), frequency agility signal each point is calculated Instantaneous phase：N=1,2 ..., N, wherein N are the total points obtained by frequency agility signal sampling；
3b) according to instantaneous phaseTime frequency analysis are carried out to frequency agility signal using Phase difference, frequency is extracted The instantaneous frequency of agile signal：N=1,2 ..., N, wherein fsFor the sample rate of signal, N is the total points obtained by frequency agility signal sampling；
3c) according to instantaneous frequency f of frequency agility signal each pointc(n), calculated rate agile signal transient frequency it is average Value；
3d) number M of the instantaneous frequency of statistic frequency agile signal near meansigma methodss, if meeting 90%m≤M≤m, The agile mode of the frequency agility signal being judged for frequency agility in arteries and veins, being otherwise judged as frequency agility between arteries and veins, m is frequency agility signal Pulse number；
3e) according to frequency agility signal reality, discrete value I (n) of imaginary part and Q (n), respectively to the point in each pulse Real part, imaginary part do fast Fourier transform and obtain corresponding real part I (ω) and imaginary part Q (ω), and calculate its instantaneous power：
3f) according to instantaneous power spec (ω) of frequency agility signal, agile number C of agile mode and victory in arteries and veins are calculated Frequency freq (i), i=1,2 ..., C；
3f1) statistics obtains codomain spec (j) of instantaneous power spec (ω), and j=1,2 ..., N, wherein N are frequency agility Total points obtained by signal sampling；
3f2) it is calculated momentary average value avg：Avg=(spec (1)+spec (2)+...+spec (j)+...+spec (N))/ N；
3f3) value of each codomain spec (j) judged successively, if meeting spec (j)>Avg and spec (j-2)< Spec (j-1), spec (j-1)<Spec (j), spec (j)>Spec (j+1), spec (j+1)>Spec (j+2), then by the point note For a spectrum peak pSpec (ω)；
The all of spectrum peak pSpec (ω) for finding out 3f4) is counted, each spectrum peak pSpec (ω) is compared successively Size, determine maximum spectrum peak value pSpecmax(ω) and maximum spectrum peak value number, maximum spectrum peak number is Agile number C；
3f5) pulse width of each pulse is divided into into agile number section, its fitting frequency is calculated each section respectively and is Agile frequency freq (i), i=1,2 ..., C；
3g) according to instantaneous power spec (ω) of frequency agility signal, agile number C of agile mode and victory between arteries and veins are calculated Frequency freq (i), i=1,2 ..., C；
3g1) calculate each pulse of frequency agility signal its fitting frequency freq (k) respectively, k=1,2 ..., m, its Middle m is the pulse number of frequency agility signal；
3g2) search the maximum freq drawn in fitting frequency freq (k)maxAnd its position of place pulse；
3g3) according to the maximum freq of frequency agility signal fitting frequencymaxThe pulse position at place, calculates two neighboring The difference of maximum place pulse position, the difference is agile number C；
3g4) take the fitting frequency of front agile number C, as agile frequency freq (i), i=1,2 ..., C.
The present invention can enter one to the time domain of frequency agility signal, the effect of modulation domain parameter measurement by following emulation experiment Step is demonstrate,proved.
Experiment 1, with time domain parameter measurement of the inventive method to frequency agility signal Matlab emulation, parameter setting are carried out For：Agile mode is agile between arteries and veins, and pulse width is 10us, and the pulse repetition period is 60us, and complete pulse number is 6, sampling Rate is 20MHz, bandwidth 20MHz.The time domain parameter for measuring is as shown in table 1：
The frequency agility signal time domain parameter table of table 1
As it can be seen from table 1 the pulse number of measurement is 6, pulse width is 10us, and the repetition period is 60us, during rising Between be 0.05us, fall time is 0.05us, and the measurement parameter is consistent with the parameter of actual set.
Experiment 2, with modulation domain parameter measurement of the inventive method to frequency agility signal Matlab emulation is carried out, and parameter sets It is set to：Bandwidth 20Mhz, peak frequency, minimum frequency are 10.2e+008 and 9.8e+008 in arteries and veins.The each intra-pulse modulation of measurement gained Parameter is as shown in table 2:
The frequency agility signal intra-pulse modulation parameter list of table 2
From table 2 it can be seen that agile mode is agile between arteries and veins, agile number is 3, and agile frequency is respectively 0MHz, 1.5MHz, 3MHz, accurately, effect is fine for measurement result.
Experiment 3, carries out time frequency analysis and obtains with existing STFT methods and Wigner-Vile methods to frequency agility signal respectively To corresponding time-frequency curve, measurement result is as shown in Figure 2 and Figure 3.Wherein Fig. 2 be the frequency agility signal that measures of STFT methods when Frequency division Butut, Fig. 3 is the time frequency distribution map of the frequency agility signal that Wigner-Vile methods are measured.
From figure 2 it can be seen that the resolution of the time frequency distribution map obtained with existing STFT methods is poor.
From figure 3, it can be seen that using existing Wigner-Vile location modes, the time-frequency distributions of frequency agility signal go out Show " falseness ", useful time varying spectrum figure thickens, therefore STFT and Wigner-Vile distributions are unsuitable for analyzing frequency victory Varying signal.
Experiment 4, with the inventive method, time frequency analysis is carried out to frequency agility signal and obtains instantaneous phase curve and time-frequency song Line, measurement result as shown in figure 4, wherein Fig. 4 (a) for frequency agility signal instantaneous phase curve, Fig. 4 (b) be frequency agility The time-frequency curve of signal, Fig. 4 (c) for frequency agility signal instantaneous phase curve frequency be 1.5MHz when enlarged drawing, Fig. 4 D () is enlarged drawing of the instantaneous phase curve of frequency agility signal when frequency is 3MHz.
Can be seen that from Fig. 4 simulation results, and the measuring method of the frequency agility signal parameter of the present invention, with algorithm letter Singly, there is no " cross term " interference in the features such as operand is little, be particularly suitable for the measurement of big data quantity in engineering practice.
In sum, the measuring method of frequency agility signal parameter of the invention more precisely can be believed frequency agility Number time frequency analysis are carried out, be more quickly obtained each time domain parameter and modulation field parameter.

Claims (5)

1. the temporal modulation field parameter union measuring method of a kind of frequency agility signal, including：
(1) discrete values of measurement frequency agile signal；
(1a) numerical values recited for obtaining frequency agility signal real part I (t) to be measured and imaginary part Q (t) is read with signal analyzer；
(1b) sample rate f is sets, the frequency agility signal real part to be measured and imaginary part of sampling 1 second time span of selection are obtained Discrete value I (n) of real part and discrete value Q (n) of imaginary part of this frequency agility signal；
(2) time domain parameter of frequency agility signal is obtained by discrete value：
(2a) according to discrete value I (n) of real part and discrete value Q (n) of imaginary part of frequency agility signal, frequency agility signal is calculated Envelope：
(2b) approximate top value A' of frequency agility signal is estimated according to frequency agility signal envelope y (n)topWith approximate floors A 'ba, obtaining approximate range measured value is：A'=A'top-A'ba；
(2c) with 10% approximate range measured value 10%A ' as rising edge and the approximate detection threshold of trailing edge, detection is obtained The pulse number m of frequency agility signal, the points being located with the approximate detection threshold of the rising edge of each pulse and trailing edge The points that approximate detection threshold is located are respectively divided by sample rate, during the approximate judgement of the rising edge for obtaining each corresponding pulse Carve t'RiseWith the approximate judgement moment t' of trailing edgeDrop, the pulse width for calculating pulse signal each pulse is：τ=t'Drop-t'Rise；
(2d) in the pulse width of each pulse, with Density Distribution averaging method, each pulse of frequency agility signal is calculated Accurate top value AtopWith floors Aba, and try to achieve Precise amplitudes value A of each frequency agility signal and be：A=Atop-Aba；
(2e) according to each pulse maximum V of frequency agility signalmaxWith minima Vmin, calculate the overshoot value of each pulse SoverWith undershoot value Sunder：
Sover=Vmax-Atop, Sunder=Aba-Vmin；
(2f) according to accurate top value A of each pulse of frequency agility signaltop, floors AbaWith Precise amplitudes value A, frequency is calculated Reference line M in the pulse amplitude of agile signalOnWith reference line M under pulse amplitudeUnder：
MOn=Aba+ 90% × A, MUnder=Aba+ 10% × A；
(2g) according to reference line M in the pulse amplitude of frequency agility signalOn, reference line M under pulse amplitudeUnderWith pulse width τ, Corresponding moment value t of reference line in pulse amplitude in the interior rising edge of lookup respectively of pulse widthOn r, reference line pair under pulse amplitude Moment value t answeredUnder rWith corresponding moment value t of reference line in pulse amplitude in trailing edgeOn f, reference line is corresponding under pulse amplitude Moment value tUnder f, and the moment t' in adjacent rising edge next time in pulse amplitude corresponding to reference lineUnder r, calculate frequency agility The rise time t of signal pulser, fall time tf, pulse period T, shut-in time toffWith dutycycle dt：
tr=tOn r-tUnder r, tf=tUnder f-tOn f, T=t'Under r-tUnder r, toff=T- τ, dt=τ/T；
(3) the modulation field parameter of frequency agility signal is obtained by discrete value：
(3a) according to the real of frequency agility signal, discrete value I (n) of imaginary part and Q (n), the instantaneous phase of frequency agility signal is calculated Position：
(3b) according to the instantaneous phase of frequency agility signalFrequency division when being carried out to frequency agility signal using Phase difference Analysis, obtains the instantaneous frequency of frequency agility signal：Wherein fsFor the sample rate of signal；
(3c) according to instantaneous frequency f of frequency agility signalcN (), calculates meansigma methodss F of frequency agility signal transient frequencyavg；
(3d) according to meansigma methodss F of frequency agility signal transient frequencyavgWith pulse number m, statistic frequency agile signal transient frequency Number M of the rate near meansigma methodss, if meeting 90%m≤M≤m, judges the agile mode of the frequency agility signal in arteries and veins Agile, is otherwise judged as agile between arteries and veins；
(3e) according to frequency agility signal reality, discrete value I (n) of imaginary part and Q (n), respectively to the reality of the point in each pulse Portion, imaginary part do fast Fourier transform and obtain corresponding real part I (ω) and imaginary part Q (ω), and calculate its instantaneous power：
(3f) according to instantaneous power spec (ω) of frequency agility signal, the victory of agile mode between agile mode or arteries and veins in arteries and veins is calculated Become number C and agile frequency freq (i), i=1,2 ..., C.
2. method according to claim 1, wherein estimating pulse according to signal envelope y (n) in the step (2b) Approximate top value A'topWith approximate floors A'ba, carry out as follows：
2b1) search and draw maximum y in signal envelope y (n)maxWith minima ymin, wherein n=1,2 ..., N, N are that frequency is prompt Total points obtained by varying signal sampling；
Signal envelope y (n) 2b2) is divided into from small to large 100 intervals, the minima of k-th interval E (k) span is calculated Emin(k) and maximum Emax(k)：
E m i n ( k ) = y m i n + ( k - 1 ) y m a x - y min 100 , E m a x ( k ) = y min + k y m a x - y min 100 , k = 1 , 2 , ... , 100 ;
2b3) statistical signal envelope y (n) falls the number in each interval E (k), is designated as count value C (k) respectively, k=1, and 2 ... 100；
The position that the maximum of front 50 count values is located 2b4) is searched in count value C (k), k is designated as1, 50 countings after lookup The position that the maximum of value is located, is designated as k2；
2b5) computation interval E (k1) meansigma methodss, the as approximate floors A' of pulseba, computation interval E (k2) meansigma methodss, as Approximate top value A' of pulsetop。
3. method according to claim 1, wherein being detected as rising edge, trailing edge using 10%A' in the step (2c) Thresholding, the pulse number of measurement frequency agile signal, is carried out as follows：
2c1) pulse number of frequency agility signal is designated as into m, and it is 0 to initialize m, is opened from the original position of frequency agility signal Begin to search first approximate rising edge, lookup result be designated as into flag bit flag, if searching successfully, flag bit flag is put into 1, Otherwise, flag bit flag is set to 0；
After 2c2) finding approximate rising edge, approximate trailing edge is searched in continuation backward, if searching successfully, flag bit flag is put 1, otherwise, flag bit flag is set to 0；
After 2c3) finding approximate trailing edge, approximate rising edge is searched in continuation backward, if searching successfully, flag bit flag is put 1, and make the pulse number m of frequency agility signal plus 1, otherwise, flag bit flag is set to 0；
The currency of flag bit flag 2c4) is detected, if flag bit flag=1, return to step 2c2), if flag bit flag= 0, then terminate to search.
4. method according to claim 1, wherein calculate in the step (3f) in arteries and veins agile number C of agile mode and Agile frequency freq (i), is carried out as follows：
3f1) statistics obtains codomain spec (j) of instantaneous power spec (ω) of frequency agility signal, wherein j=1,2 ..., N, N Total points obtained by frequency agility signal sampling；
3f2) calculate momentary average value：Avg=(spec (1)+spec (2)+...+spec (j)+...+spec (N))/N；
3f3) value of each codomain spec (j) judged successively, if meeting spec (j)>Avg and spec (j-2)<spec (j-1), spec (j-1)<Spec (j), spec (j)>Spec (j+1), spec (j+1)>Spec (j+2), then be designated as one by the point Individual spectrum peak pSpec (ω)；
The all of spectrum peak pSpec (ω) for finding out 3f4) is counted, the big of each spectrum peak pSpec (ω) is compared successively It is little, determine maximum spectrum peak value pSpecmax(ω) and maximum spectrum peak value number, maximum spectrum peak number is agile Number C；
It is 3f5) agile number section by pulse width τ of each pulse point, its fitting frequency is calculated each section respectively and is victory Frequency freq (i), i=1,2 ..., C.
5. method according to claim 1, wherein calculate in the step (3f) between arteries and veins agile number C of agile mode and Agile frequency freq (i), is carried out as follows：
3f6) respectively its is calculated to each pulse of frequency agility signal and be fitted frequency freq (k), k=1,2 ..., m, wherein m are The pulse number of frequency agility signal；
3f7) search the maximum freq drawn in fitting frequency freq (k)maxAnd its position of place pulse；
3f8) according to the maximum freq of frequency agility signal fitting frequencymaxThe pulse position at place, calculates two neighboring maximum The difference of value place pulse position, the difference is agile number C；
3f9) take the fitting frequency of front agile number C, as agile frequency freq (i), i=1,2 ..., C.
CN201611246663.7A 2016-12-29 2016-12-29 Time-domain modulation domain parameter combination measure method of frequency agile signal Pending CN106680795A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611246663.7A CN106680795A (en) 2016-12-29 2016-12-29 Time-domain modulation domain parameter combination measure method of frequency agile signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611246663.7A CN106680795A (en) 2016-12-29 2016-12-29 Time-domain modulation domain parameter combination measure method of frequency agile signal

Publications (1)

Publication Number Publication Date
CN106680795A true CN106680795A (en) 2017-05-17

Family

ID=58872459

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611246663.7A Pending CN106680795A (en) 2016-12-29 2016-12-29 Time-domain modulation domain parameter combination measure method of frequency agile signal

Country Status (1)

Country Link
CN (1) CN106680795A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110579631A (en) * 2019-10-08 2019-12-17 中国航空工业集团公司北京长城计量测试技术研究所 Accurate measurement method for time established by agile signal based on sine fitting
CN110749878A (en) * 2019-09-24 2020-02-04 北京理工大学 Intra-pulse frequency coding signal parameter measurement method based on digital channelization
CN114089290A (en) * 2021-11-12 2022-02-25 扬州宇安电子科技有限公司 Method for rapidly and efficiently measuring frequency and guiding noise interference
CN115993484A (en) * 2023-02-27 2023-04-21 南京天朗防务科技有限公司 Band-limited signal instantaneous frequency measurement method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1418441A1 (en) * 2002-09-23 2004-05-12 Thales Nederland B.V. Agile PRT Deconvolution method and systems, and its uses
CN103064073A (en) * 2012-12-04 2013-04-24 上海无线电设备研究所 Method based on frequency agility for changing radar target properties
CN104330782A (en) * 2014-11-04 2015-02-04 西安电子科技大学 Time domain and modulation domain parameter combined measuring method of triangular frequency-modulation pulse signals
US9069066B2 (en) * 2013-05-20 2015-06-30 Faran Awais Butt Radar deception jamming prevention using bi-static and mono-static radars
CN105785338A (en) * 2016-03-09 2016-07-20 西安电子科技大学 Method for optimizing carrier frequency of frequency-agile radar

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1418441A1 (en) * 2002-09-23 2004-05-12 Thales Nederland B.V. Agile PRT Deconvolution method and systems, and its uses
CN103064073A (en) * 2012-12-04 2013-04-24 上海无线电设备研究所 Method based on frequency agility for changing radar target properties
US9069066B2 (en) * 2013-05-20 2015-06-30 Faran Awais Butt Radar deception jamming prevention using bi-static and mono-static radars
CN104330782A (en) * 2014-11-04 2015-02-04 西安电子科技大学 Time domain and modulation domain parameter combined measuring method of triangular frequency-modulation pulse signals
CN105785338A (en) * 2016-03-09 2016-07-20 西安电子科技大学 Method for optimizing carrier frequency of frequency-agile radar

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
苗雪平: ""基于LabWindows/CVI的雷达脉冲信号分析软件设计及实现"", 《万方数据》 *
谭信: ""捷变频雷达信号处理仿真"", 《舰船电子对抗》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110749878A (en) * 2019-09-24 2020-02-04 北京理工大学 Intra-pulse frequency coding signal parameter measurement method based on digital channelization
CN110749878B (en) * 2019-09-24 2021-08-13 北京理工大学 Intra-pulse frequency coding signal parameter measurement method based on digital channelization
CN110579631A (en) * 2019-10-08 2019-12-17 中国航空工业集团公司北京长城计量测试技术研究所 Accurate measurement method for time established by agile signal based on sine fitting
CN110579631B (en) * 2019-10-08 2021-09-21 中国航空工业集团公司北京长城计量测试技术研究所 Accurate measurement method for time established by agile signal based on sine fitting
CN114089290A (en) * 2021-11-12 2022-02-25 扬州宇安电子科技有限公司 Method for rapidly and efficiently measuring frequency and guiding noise interference
CN115993484A (en) * 2023-02-27 2023-04-21 南京天朗防务科技有限公司 Band-limited signal instantaneous frequency measurement method

Similar Documents

Publication Publication Date Title
CN104330782B (en) Time domain and modulation domain parameter combined measuring method of triangular frequency-modulation pulse signals
CN106680795A (en) Time-domain modulation domain parameter combination measure method of frequency agile signal
CN109459732A (en) A kind of method for identifying and classifying of wideband radar interference signal
CN115047417B (en) Method and system for simulating radar echo and interference
CN111580091B (en) Sea surface weak target detection method based on AR spectrum singular intensity function
CN110988834B (en) Pulse arrival time measuring method based on self-adaptive threshold value double thresholds
WO2022170665A1 (en) Radar signal waveform uncertainty test system
EP2375663A1 (en) Method for classifying signals according to the modulation type used
CN109521412B (en) Radar networking airspace target detection method based on local statistic fusion
CN104122538A (en) Method for determining noise power of wind profile radar
CN104483668A (en) High-accuracy radar signal detecting and tracking system and method
JP5752481B2 (en) Passive radar device, guidance device, and radio wave detection method
CN111239704A (en) Atmospheric detection radar target echo identification processing method, device, equipment and medium
CN102638290A (en) Method and device for extracting multipath signal component based on channel measurement
CN112764000B (en) Method for estimating parameters of sea clutter amplitude model of lognormal texture based on quantiles
CN106772271A (en) A kind of radar signal TOA methods of estimation based on matching auto-correlation algorithm
Allen et al. Robust statistics for deterministic and stochastic gravitational waves in non-Gaussian noise. II. Bayesian analyses
CN110333506B (en) Method for extracting inhaul cable position parameters of cable force measurement radar
CN101545969B (en) Method for evaluating point target resolution of side-glance SAR
CN106597410A (en) Joint measurement method of time domain and modulation domain parameters of barker code signal
CN101404546A (en) Frequency spectrum occupancy rate measuring method based on signal moment characteristics
CN105959035A (en) Direct sequence spread spectrum signal interception detection method
CN109031260A (en) A kind of LFM signal time delay measurement method based on the analysis of fractional Fourier modulation rate
US20230168338A1 (en) Pulse compressed radar signal detection method
CN116449304A (en) SAR emission pulse arrival time measurement method based on frequency measurement

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170517

WD01 Invention patent application deemed withdrawn after publication