EP3123648B1 - Simultaneous call transmission detection - Google Patents
Simultaneous call transmission detection Download PDFInfo
- Publication number
- EP3123648B1 EP3123648B1 EP15721779.5A EP15721779A EP3123648B1 EP 3123648 B1 EP3123648 B1 EP 3123648B1 EP 15721779 A EP15721779 A EP 15721779A EP 3123648 B1 EP3123648 B1 EP 3123648B1
- Authority
- EP
- European Patent Office
- Prior art keywords
- frequency
- signal
- domain sum
- domain
- primary
- 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.)
- Active
Links
- 230000005540 biological transmission Effects 0.000 title description 35
- 238000001514 detection method Methods 0.000 title description 28
- 238000000034 method Methods 0.000 claims description 49
- 208000032370 Secondary transmission Diseases 0.000 claims description 29
- 238000004458 analytical method Methods 0.000 claims description 20
- 238000001228 spectrum Methods 0.000 claims description 16
- 208000032369 Primary transmission Diseases 0.000 claims description 15
- 230000001131 transforming effect Effects 0.000 claims description 10
- 230000002238 attenuated effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 description 17
- 238000006243 chemical reaction Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 16
- 230000000694 effects Effects 0.000 description 15
- 239000000872 buffer Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 13
- 238000004590 computer program Methods 0.000 description 9
- 238000005070 sampling Methods 0.000 description 8
- 230000001419 dependent effect Effects 0.000 description 7
- 238000001914 filtration Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000000654 additive Substances 0.000 description 3
- 230000000996 additive effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000012952 Resampling Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005236 sound signal Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000282461 Canis lupus Species 0.000 description 1
- 206010011469 Crying Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012505 colouration Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 230000002939 deleterious effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 231100001261 hazardous Toxicity 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/02—Amplitude-modulated carrier systems, e.g. using on-off keying; Single sideband or vestigial sideband modulation
- H04L27/06—Demodulator circuits; Receiver circuits
- H04L27/066—Carrier recovery circuits
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J11/00—Orthogonal multiplex systems, e.g. using WALSH codes
- H04J11/0023—Interference mitigation or co-ordination
- H04J11/0026—Interference mitigation or co-ordination of multi-user interference
- H04J11/0036—Interference mitigation or co-ordination of multi-user interference at the receiver
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0017—Arrangements for implementing traffic-related aircraft activities, e.g. arrangements for generating, displaying, acquiring or managing traffic information
- G08G5/0026—Arrangements for implementing traffic-related aircraft activities, e.g. arrangements for generating, displaying, acquiring or managing traffic information located on the ground
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0095—Aspects of air-traffic control not provided for in the other subgroups of this main group
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03D—DEMODULATION OR TRANSFERENCE OF MODULATION FROM ONE CARRIER TO ANOTHER
- H03D1/00—Demodulation of amplitude-modulated oscillations
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/02—Amplitude-modulated carrier systems, e.g. using on-off keying; Single sideband or vestigial sideband modulation
- H04L27/06—Demodulator circuits; Receiver circuits
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2649—Demodulators
- H04L27/265—Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
- H04L27/2651—Modification of fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators for performance improvement
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J11/00—Orthogonal multiplex systems, e.g. using WALSH codes
- H04J11/0023—Interference mitigation or co-ordination
- H04J11/0026—Interference mitigation or co-ordination of multi-user interference
Definitions
- This invention relates to an apparatus and method for the detection of simultaneous call transmissions, in particular to double-sideband amplitude modulation (DSB-AM) transmissions.
- DSB-AM double-sideband amplitude modulation
- SCT Simultaneous Call Transmission
- the first condition may not be met if one transmission originates much further away than the other (i.e. one plane is overhead and the other is several kilometres away); the second condition may not be met if the frequency difference between the two transmitters (defined by the precision of the quartz in the transmitting equipment) is minor.
- SCT can occur far outside of these values due to real world effects such as propagation loss, multipath error and frequency error.
- WO2015/008165A2 describes a method for the detection of more than one signals contained in a received signal. The method comprises: down-converting the received signal, thereby providing a down-converted signal in a complex IQ base band; and at least partially cancelling the strongest user in the down-converted signal, thereby allowing for the detection of a possible secondary user.
- US2010/0067570A1 and DE102007037105A1 describe an apparatus adapted to automatically detect when two transmissions occur simultaneously.
- This system phase demodulates the sum-signal of the dual transmission by converting the baseband In-Phase and Quadrature signal into unwrapped phase and amplitude.
- a periodic 'wobble' is present on the unwrapped phase if there are close frequency-separated transmissions, as the difference in frequency causes the phasors of the transmissions to rotate around one another.
- This phase time-series is then Fourier transformed using a bank of transformers with varying window lengths to determine whether any secondary transmission (i.e. a peak due to the unwrapped phase 'wobble') is present.
- a warning tone is added to the audio output if a secondary transmission is detected so as to alert the operator to the situation.
- the step of producing a phase time series is an intrinsically non-linear process (as it involves performing an arctan) and so errors propagate as the process progresses manifesting as intermodulation products in the output spectrum which are not physically present in the input.
- this solution can potentially generate large occurrences of 'false positives' where an alert is sounded when only one transmission is present. This is because this system has no way of suppressing common signal impairments such as sinusoidal mains hum, incidental FM (frequency modulation), frequency-selective multipath effects and 1/f 2 (reciprocal frequency-squared) phase noise. All of these effects can potentially introduce a certain amount of unwrapped phase 'wobble' and then be identified as secondary transmissions.
- a method of determining the presence of a secondary transmission in a time-domain sum-signal including a primary transmission and the secondary transmission comprising: transforming the time-domain sum-signal into a frequency domain sum-signal; wherein the frequency-domain sum-signal is a linear combination of the primary transmission and the secondary transmission; and wherein the transforming is based on a plurality of frequency bins; estimating a primary carrier frequency based on the frequency-domain sum-signal; shifting the frequency bins of the frequency-domain sum-signal based on the estimated primary carrier frequency; determining a conjugate of a sideband of the frequency-bin-shifted frequency-domain sum-signal; attenuating the frequency-bin-shifted frequency-domain sum-signal by using said conjugate of the sideband of the frequency-bin-shifted frequency-domain sum-signal; extracting at least one peak from the attenuated frequency-domain sum-signal; and determining the presence of the secondary transmission in the frequency-domain sum-signal
- an apparatus for determining the presence of a secondary transmission in a time-domain sum-signal including a primary transmission and the secondary transmission comprising: means for transforming the time-domain sum-signal into a frequency domain sum-signal, wherein the frequency-domain sum-signal is a linear combination of the primary transmission and the secondary transmission; and wherein the transforming is based on a plurality of frequency bins; means for estimating a primary carrier frequency based on the frequency-domain sum-signal; means for shifting the frequency bins of the frequency-domain sum-signal based on the estimated primary carrier frequency; means for determining a conjugate of a sideband of the frequency-bin-shifted frequency-domain sum-signal; means for attenuating the frequency-bin-shifted frequency-domain sum-signal by using said conjugate of the sideband of the frequency-bin-shifted frequency-domain sum-signal; means for extracting at least one peak from the attenuated frequency-domain sum-signal; and means for determining the presence of
- the invention also provides a computer program and a computer program product comprising software code adapted, when executed on a data processing apparatus, to perform any of the methods described herein, including any or all of their component steps.
- the invention also provides a computer program and a computer program product comprising software code which, when executed on a data processing apparatus, comprises any of the apparatus features described herein.
- the invention also provides a computer program and a computer program product having an operating system which supports a computer program for carrying out any of the methods described herein and/or for embodying any of the apparatus features described herein.
- the invention also provides a computer readable medium having stored thereon the computer program as aforesaid.
- the invention also provides a signal carrying the computer program as aforesaid, and a method of transmitting such a signal.
- Any apparatus feature as described herein may also be provided as a method feature, and vice versa.
- means plus function features may be expressed alternatively in terms of their corresponding structure, such as a suitably programmed processor and associated memory.
- the term 'primary transmission' or 'primary carrier' refers to the transmission with greatest power.
- the term 'secondary transmission' or 'secondary carrier' refers to any other (lower-powered) transmission occurring at the same time from another aircraft.
- Air traffic controller (ATC) to aircraft communication is generally very concise, each transmission is typically less than 10 seconds in duration, and can be as short as 2 or 3 seconds. For this reason, the latency in detecting a secondary call transmission (SCT) is preferably less than 2-3 seconds for this field of use.
- SCT secondary call transmission
- the term 'simultaneous' in this specification refers to the situation where two transmissions overlap in time, as in this scenario an ATC would either not hear the SCT or the audio would be filtered out by the radio.
- the primary carrier signal i.e. ⁇ c and c
- properties of the primary carrier signal i.e. ⁇ c and c
- the primary carrier Once the primary carrier has been identified, it can be isolated and removed, allowing subsequent analysis of Y( ⁇ ) to determine whether an SCT has occurred.
- Figure 3 illustrates a high-level dataflow of an exemplary method (referred to as 'Frequency Domain (FD) SCT detection') for detecting the presence of a secondary call transmission in an input signal.
- the input is a high sample rate IQ (in-phase / quadrature) baseband time-series (real or complex) and the outputs are SCT detection results (for example, as a tone inserted into the audio signal or a flag placed into a data stream).
- IQ in-phase / quadrature
- SCT detection results for example, as a tone inserted into the audio signal or a flag placed into a data stream.
- the term "sum-signal” preferably connotes a signal received by a receiver (for example an air traffic control system).
- a sum signal preferably comprises a primary carrier signal, noise, and potentially a secondary carrier signal (where present); the term “sum” merely implies, preferably, that such elements (primary carrier signal, noise, etc) are present in one and the same signal.
- a sum signal is received and converted to an IQ baseband sample.
- This signal is decimated 300 (downsampled) removing unnecessary frequency components so as to reduce the load on later processing steps.
- the decimated signal is sampled using an overlapping, sliding window buffer 302 which stores a given length of signal to be processed. Sampling the signal means that the entire signal does not need to be processed in one go, improving the latency of detection and reducing processor load.
- the length of the buffer is determined by the trade-off between processor load in analysing lengthy samples and the increased detection accuracy (i.e. high signal-to-noise ratio) longer samples afford.
- the sampling rate i.e. number of windows per second
- the sampling rate is determined by a trade-off between processor load in processing large numbers of windows within a short time-period and latency of SCT detection.
- Each window sample is then input into a Fast Fourier Transform (FFT) 304 - outputting X( ⁇ ) - for further processing in the frequency domain.
- FFT Fast Fourier Transform
- An FFT is preferable to a continuous Fourier Transform (FT) as it is much less processor-intensive.
- Other discrete transforms such as wavelet transforms or spectral line filters may be used.
- the decimation step may introduce 'decimator ripple' in the frequency domain output, this is corrected for at step 306 before any further processing.
- the frequency of the primary carrier transmission is estimated 308, and the FFT output is down-converted 310 based on this frequency.
- the in-phase elements to the primary carrier transmission are then cancelled by subtracting the conjugate of the negative sideband frequencies of the primary signal from their counterpart positive sideband frequencies in the (decimated, down-sampled, frequency domain) sum signal 312.
- the remaining signal - Y( ⁇ ) - is due to phase noise or other signals at different frequencies (which may be SCTs) to the main carrier transmission in the signal.
- This signal is analysed for peaks 316 (defined by a threshold) above an estimated noise floor 314.
- peaks 316 are indicative of an SCT being present as they represent significant magnitude parts of the signal which are not on the same carrier frequency as the primary transmission (e.g. peaks corresponding to a heterodyne tone).
- mains hum and phase noise may show up as peaks above a nominal noise floor.
- Noise effects such as these affecting the sidebands of the primary carrier are typically symmetric about the primary carrier frequency, thus an asymmetry analysis 318 is performed to determine whether a particular peak has a corresponding 'mirror image' peak. This analysis is performed using an 'asymmetry threshold'. The peak is also analysed for its magnitude (above the noise floor) as higher power peaks are more likely to be secondary transmissions rather than variations in noise. These two parameters (and/or others) are combined in a 'Feature Space Classification' 320 and an SCT can be signalled if a sample contains a peak exceeding the predetermined threshold(s).
- Figure 4 shows a schematic diagram of a radio receiver 104 adapted to perform the processes involved in the detection of a secondary call transmission as described above.
- a signal is received by an aerial and is input to an Analogue-to-Digital Converter (ADC) 402.
- ADC Analogue-to-Digital Converter
- the dotted section 400 represents a simplified digital radio receiver without any SCT detection capability.
- the digital signal is demodulated by demodulating unit 403, with assistance from a central processor 422 and memory 424. This is then passed to an audio output unit 420 and audio is outputted.
- An actual digital radio may include many additional components (such as tuning, filtering and amplifying circuitry), but such components are omitted for clarity in this figure.
- This audio extraction process occurs independently of the SCT detection, as this represents the primary purpose of the radio 104, to convert received signals into audio (or other useful information).
- the components used for SCT detection are shown outside of the dotted section 400.
- the digital signal is decimated (downsampled) by decimator 404 before being sampled in a sliding window buffer 406. Each window is then passed through a Fast Fourier Transform (FFT) 408.
- FFT Fast Fourier Transform
- the spectrum outputted from the FFT has filters / windows 410 applied to it so as to produce the signal defined by Y( ⁇ ) in equation 4.
- This output is passed to a comparator 412, which, with logic circuitry 412 and thresholds stored in memory 424, determines whether an SCT has occurred. If so, the operator is notified, for example by a tone being inserted into the audio output and/or a flag (such as an indicator on a user interface) is raised via tone / flag generator 418.
- Other information regarding the SCT such as indication of a confidence level, or timestamp of the event, may also be outputted.
- Figure 4 shows components separated for clarity whereas in reality many of these components may be combined as a single component (such as the comparator, logic combined with the processor) or further split into separate components.
- the intermediate digital signal after analogue to digital conversion is often at a higher sample rate than is required to support the DSB-AM sidebands of a primary signal.
- decimation is twofold: (1) to reduce the computational load and (2) reject signals outside the band of interest for SCT detection.
- the decimation step (and the subsequent ripple equalisation) would not be necessary if these issues are not of any relevance (e.g. if the analogue-to-digital converter has a low sample rate).
- the decimator design preferably has a narrow transition region (for example 10% of passband) with low passband ripple (0-3dB) and high stopband attenuation (for example more than 40dB) i.e. a typical specification for a high-quality decimator for audio applications.
- a typical low-pass mask specification given the third set of parameters in the above table would be ⁇ 1dB of passband ripple up to 5kHz, a transition region from 5kHz to 7kHz, and -60dB of gain in the decimation stopband.
- time-series is real-only, then a complex oscillator and mixer are required to down-convert the signal before the decimator.
- the decimation low-pass filter then requires sufficient stopband attenuation to adequately remove the frequency-shifted conjugate image.
- This stage presents the most recent T seconds block of sample data M times a second to the subsequent processing stages as illustrated in Figure 5 .
- Blocks preferably have a high degree of overlap to maximise the chance of detecting a secondary transmission as soon as it starts.
- not all SCT events may be detected in such a short time period as the signal may be too weak; in such circumstances, due to the overlap of windows, the next window would have 0.5 seconds of SCT to detect and so on.
- T*M>1 but ideally at around 8 windows overlap, so T*M ⁇ 8.
- the purpose is to allow strong signals to be detected quickly by the system at a coarse granularity in time, but also allow an adequate time-history to allow the coherent integration and detection of weak secondary signals.
- T can be in the range 1 to 4 seconds and M can be in the range 2 to 16.
- T governs the coherent integration period for detecting weak SCT signals and it is of advantage to be long and about the same length as typical primary transmission utterances.
- the value of (1/M) governs the maximum latency for detecting strong SCT signals and it is advantageous for M to be high for low latency.
- the processing load is proportional to the product T*M, thus there is a trade-off between performance and processor load when selecting values of T and M.
- Figure 6 illustrates how the current analysis window from the buffer is mapped into the FFT input with 'zero padding'.
- the mapping is unconventional in that the buffer is split into two halves and the first half is mapped into the final part of the FFT input and the second half is mapped into the start of the FFT input, with zero-value inputs occupying the intervening samples. This improves the operation of frequency domain down-conversion as described below.
- the FFT size N FFT is chosen to be around twice the size of the buffer window to provide sufficient oversampling for subsequent processing.
- the choice of oversampling ratio at approximately x 2 is a compromise between two conflicting factors: (1) critical sampling at around x 1 oversampling is not viable because down-conversion requires an unfeasibly long resampling filter for the required DSB-AM cancellation fidelity, (2) the system performance at say, > x 3 oversampling yields negligible performance benefit at the expense of increased computational complexity in the FFT.
- oversampling ratios of greater than x 3 may be used if computational complexity is not an issue - for example if the fidelity of the cancelation is paramount.
- oversampling ratio may be used. The larger the oversampling, the more processor-intensive the resulting analysis would be (due to the greater number of discrete frequency 'bins' in the frequency domain) but the system would be more accurate due to (at least) the spectrum having greater resolution.
- the low-pass filter discussed above with reference to decimation may have significant passband ripple in order to be implementable with realistic cost.
- Passband ripple is an artefact manifesting in the spectrum of a transformed signal having had imperfect (i.e. non-square) band-pass filters applied to it.
- the inverse transform 1/ H ( ⁇ ) is stored as a vector of N FFT complex weights which is applied to the output directly after the FFT has been computed.
- H ( ⁇ ) is symmetric about zero hertz, it is not symmetric about the primary carrier, so would not be cancelled out when calculating Y( ⁇ ) - which is described in more detail below.
- the highest magnitude FFT output bin (denoted as bin j ) is detected and its power and frequency are measured. This is asserted to be the primary carrier (i.e. strongest sinusoidal tone) and these measurements are passed on to the classification stage discussed in order to detect if any primary signal is present. Identifying the primary carrier frequency leads to identification of non-primary carrier signals (such as an SCT).
- a parabolic (quadratic) curve may be fitted to the points, for example using closed-form linear algebra.
- the fractional bin frequency f in the range of -0.5 to 0.5 of the maximum value of the fitted parabola is taken to be the best estimate of the true primary carrier frequency ⁇ c .
- the oversampling of the FFT e.g. twice oversampling
- the width of the primary mainlobe may be assessed by searching out from the peak in both negative and positive frequency until bins that are ⁇ 3dB (approximately ⁇ 0.5 in power) of the peak are identified (i.e. full-width, half maximum (FWHM) of the primary mainlobe).
- FWHM full-width, half maximum
- Frequency-domain down-conversion 310 is performed by generating a finite-impulse filter which shifts the frequency bins by -( j + f ) bins (i.e. by ⁇ c ) so that the underlying maximum of the primary carrier mainlobe is shifted exactly on to the zero frequency bin. This step effectively makes the primary carrier signal symmetric about zero hertz, making later computation and determination of SCT events simpler.
- Other windows may be used such as 'Kaiser' or 'Equiripple' windows, but cosine-family windows such as Hamming, Hann, Blackman family have the implementation benefit of combining good sidelobe performance with the precise and simple computation using cosines.
- x LIM sets the limits for the window (i.e. it is zero-valued for
- > x LIM ) , and hence defines the quality of the resampling (a typical value would be x LIM 5) .
- a small value is desirable in order to minimise the processing complexity of down-conversion.
- x LIM is discussed in more detail below with reference to Figures 7 and 8 .
- FIR Finite Impulse Response
- Frequency domain convolution of two signals is analogous to multiplication of their time-domain equivalents.
- Other values of f create envelopes intermediate between these extremes. Although this process would not be necessary if only integer values of f were used, doing so would introduce errors into the central frequency and thus mean that the later asymmetry analysis would carry through these errors.
- Equation 5 comprises the product of two terms; (1) a sin( x )/ x function with infinite support on x (which has too many terms to compute practically) and (2) a compactly supported window function (which makes g ( x ) economic to compute).
- this is the circular convolution of (1) an arbitrary frequency sinusoid and (2) a band-pass filter corresponding to the frequency shifted IFFT of the window function.
- the output of this filtering process is unit amplitude sinusoid except where a phase discontinuity passes through the filter where the two ends of the sinusoid are circularly "spliced" together. This creates the characteristic "dip" in the sinusoid envelope illustrated in Figure 8 which is worst-case when a 180 degree continuity passes through (as occurs with the half-bin case).
- Figure 8 also explains the utility of the unconventional zero-padding described above of mapping the "1 st half" and "2 nd half' of the input buffer to time-domain intervals where the envelope function is almost exactly unity.
- the mapping of the second half of the time window to the first part of the FFT input and vice versa means that the FFT input maintains its time-order as the end of the first half is effectively contiguous with the start of the second half (as the FFT can be visualised as wrapped around the surface of a cylinder).
- x LIM is a function of the oversampling ratio N FFT f s ⁇ T so as to be the smallest value to minimise the computational complexity of the window filter whilst not impinging on the 'flatness' of the envelope function. If x LIM is too small, the envelope function would begin to curve over the sections of the IFFT which contain the signal data, resulting in the signal being modified prior to DSB-AM cancellation.
- the value of x LIM which satisfies this trade-off has been empirically found to be approximately (12 / oversampling ratio).
- DSB-AM cancellation 312 as discussed above is effected by applying Equation 8 in order to generate an output vector y comprising N FFT /2+1 bins (the zero frequency bin and the right hand side of the spectrum).
- Y ( ⁇ ) is by mathematical definition conjugate symmetric about zero for an ideal primary carrier, only computation of the right hand side (i.e. positive frequency) is necessary. Only magnitude information is taken into y for the purpose of peak detection, hence the modulus is taken.
- y i x i ⁇ x N FFT ⁇ i mod N FFT ⁇ ; i ⁇ 0,1 , ... N FFT 2
- the quality of DSB-AM cancellation 312 is dependent on the temporal coherence of the primary signal.
- Phase noise on the primary carrier can lead to some feed through of tonal components in the sidebands which may appear as distinct tones in Y ( ⁇ ).
- ⁇ ⁇
- SCT tones in y are characterised by isolated narrowband peaks against a noise floor after DSB-AM cancellation of the primary carrier transmission.
- a noise floor estimate which is not biased by tonal peaks should be estimated.
- Noise levels may not be constant over the whole frequency range in question, so the noise level at every frequency bin is estimated in order to 1) capture secondary transmissions above the local noise level, but potentially below the noise level elsewhere, and 2) discount frequency bins with higher levels of noise than elsewhere.
- a single estimate of the noise level across the entire frequency spectrum would not be able to account for such circumstances, resulting in, in the case of 1) false negatives, and in the case of 2) false positives. Either of these scenarios is undesirable, false negatives particularly so in an ATC implementation as such events could result in a dangerous situation.
- An effective way to determine a frequency-dependent noise floor estimation is to calculate a moving-average of the magnitude across a range of bins centred around a particular frequency bin. If a large enough bin range is used and peaks are not frequent, this would be an accurate representation of the noise floor at that frequency bin.
- a short sliding-window rank-order statistic filter is applied which extracts the e.g. the median, power bin as the noise floor estimate.
- Analogous filters are used for removing impulsive noise from otherwise smooth functions in applications like image processing.
- n i median y i ⁇ N NFE ... i + N NFE ; i ⁇ 0,1 , ... N FFT 2 FFT
- N NFE D 1 D 2 .
- fewer, larger windows may be taken or alternatively, more, smaller windows.
- the choice of length of windows D 1 , D 2 is also dependent on the trade-off between too short being dominated by peaks and too long missing the trend of the noise, for example D 1 and D 2 could each vary between 4 and 64 as a general illustration in these circumstances.
- the median is the default rank-order statistic to draw out, but other measures of central tendency are possible, for example the 40 th centile, which will be less biased by peaks, but more susceptible to low power noise samples.
- Peaks are identified in y by identifying local maxima, where y i >y i-1 and y i >y i+1 . Performing just this analysis may pick up a lot of spurious fluctuations in the noise floor, for this reason only peaks (i.e. values of y i ) that satisfy a certain predefined threshold (peak_metric_thresh) are identified as SCT candidates.
- peak_metric_thresh Example values for peak_metric_thresh are provided below with reference to Figures 11 to 14 , but may vary from around 0.85 to 3 (or greater than 3) depending on the situation.
- this is where the distinct peaks in y are 10 peak_metric_thresh times higher than the (local) noise floor n.
- This is denoted the subset P of the set of all possible i values (bin indices) which satisfies Equation 10.
- Equation 11 log 10 y i n i ; i ⁇ P
- the threshold peak_metric_thresh is preferably a system-set parameter which may be calculated once upon calibration of the system; alternatively it may be dynamically calculated so as to result in a system with a specific false-positive rate. This may be useful if the variance in the noise floor (i.e. the accuracy of the noise floor estimate) changes over time so that the system becomes more prone to false negatives (if the variance decreases) or it becomes more prone to false-positives (the variance increases). In an average situation, a value for the peak metric threshold p ( i ) would be between 1 and 4, more preferably between 2 and 3 as a general illustration in these circumstances.
- Another metric that may be used to reduce the number of candidate peaks is to specify that two peaks must be separated by a minimum frequency otherwise they are treated as a single peak (i.e. the smaller peak is disregarded).
- the threshold min_freq_sep is defined. In one example this is between 5Hz and 50Hz, preferably between 7Hz and 15Hz, and preferably approximately 10Hz. Disregarding the smaller peak of a closely separated pair of peaks has negligible impact on the capability to detect genuine secondary tones when peak detections are sparsely separated. Such a feature allows strong peaks from e.g. 400Hz mains hum (which are highly conjugate-symmetric) to absorb their own sidelobe features which are much weaker in power but more asymmetric and thus can cause false positives.
- the method identifies the weaker peaks from the set P which are within +/- min_freq_sep of the current secondary tone candidate being analysed, and marks them for deletion from set P by placing them in the set Q as follows (with commentary accompanying each step):
- Equation 12 a non-negative real-valued asymmetry metric is computed using Equation 12. This is a measure of how asymmetric the power is between positive and negative frequencies (with respect to the down-converted primary carrier at zero frequency).
- SCT events have (by definition) a central frequency offset from that of the primary carrier and are thus asymmetric about primary carrier (and, after down-conversion, are asymmetric about zero hertz). There is a low probability of another tone precisely at the opposite frequency sign as this would correspond to a third SCT at a very specific frequency.
- the asymmetry metric a ( i ) provides a useful way to exploit values that are pre-computed elsewhere in the process (i.e. bins from the down-converted X ( ⁇ ) in vector x ) to reject false positives from poor quality primary transmitters.
- a threshold, asym_metric_thresh, for the value of a(i) is defined where peaks not meeting this threshold are discarded as being too symmetrical, and thus unlikely to be SCTs.
- the asymmetry threshold provides a means to discount peaks which have a high residual power following subtraction due to the fact that the symmetric peaks had a high power prior to subtraction - for example if the signal has a high level of noise (which is not perfectly symmetrical), or due to external effects such as mains hum.
- Figures 12 , 13(b) and 13(c) below show scenarios where the asymmetry threshold is utilised to reduce the false-positive rate by limiting the number of events above the power threshold which would otherwise be deemed to be SCT events.
- a number of checks may be performed.
- a threshold primary_pk_thresh is defined where SCT analysis is only undertaken if the primary peak is above this threshold. This threshold is corrected by the amount of gain applied to the signal (AGC_gain) so as to measure the absolute power of the primary signal.
- a threshold is also set for the maximum allowed width of the primary peak, primary_bw_thresh, where SCT analysis is only undertaken if the width of the primary carrier peak is greater than this threshold. This ensures that a certain lower bound is met on the mark-space ratio of the primary transmitter in the analysis window, for example it may be desirable for the primary transmission to occupy at least 50% of the time window. This can prevent some anomalies due to rising edges entering or trailing edges leaving the analysis window.
- the width of the primary carrier peak is an output which is simple to generate and which provides some clear information about the temporal activity of the primary transmitter.
- AGC Automatic Gain Control
- a more generalised analysis is to fit a suitable likelihood density function of the form prob(peak_metric, asymmetry_metric) given SCT present ("H1”) or SCT absent ("H0”) and then computing a likelihood ratio to make the decision.
- the exact form of the likelihood function would depend on the application, as well as other factors such as desired false positive rate.
- a more sophisticated algorithm than the decision logic described above, with some statistical modelling of the parameter density functions under different H1/H0 hypotheses e.g. Gaussian Mixtures Model, Fuzzy Clustering, Neural Network, or Support Vector Machines
- a confidence score for example between zero and one.
- Such a confidence level could be fed-back to the end user for information and/or calibration purposes.
- Figure 9 shows the spectrum X ( ⁇ ) of the input (a) and output (b) signals of Frequency Domain Down-Conversion.
- the primary carrier, voice sidebands, 400Hz mains tone sidebands and the secondary signal (creating an SCT-present scenario) are marked.
- the primary carrier is shifted to zero-frequency making the two voice and 400Hz mains sidebands and the carrier of the single secondary signal respectively symmetric and asymmetric about zero-frequency, as discussed above.
- Figure 10 illustrates the results of the DSB-AM cancellation spectrum Y ( ⁇ ) (a) and the noise floor estimated spectrum N ( ⁇ ) (b) in comparison to the superposed positive and negative frequency halves of X ( ⁇ ).
- DSB-AM cancellation has achieved around 25dB attenuation of the 400Hz tone with negligible attenuation of the secondary carrier. This is because the 400Hz mains hum modulates the primary carrier and is thus conjugate-symmetric with respect to the primary carrier. This means this feature is largely attenuated by the proposed frequency domain DSB-AM cancellation stage. However the secondary carrier is not conjugate symmetric with respect to the primary carrier and is not significantly attenuated.
- the noise floor estimate, N ( ⁇ ) - shown in Figure 10(b) follows the underlying spectral envelope of Y ( ⁇ ) without much bias from isolated peaks in Y ( ⁇ ). Note that imperfect DSB-AM cancellation of the (semi-coherent, poor quality) primary signal has led to some feed-through of primary voice spectrum which is followed by the noise floor estimate N ( ⁇ ).
- Figure 11 illustrates detected peaks (the peak_metric_thresh is set low to a value of 0.85 to allow false detections through for characterisation). Two peaks are correctly detected for respectively the 400Hz mains tone and secondary carrier. Though the peak metrics are of comparable magnitude (shown in Figure 11(a) ), the asymmetry metrics are different ( Figure 11(b) ).
- FIG. 14 illustrates the utility of such thresholds. For illustration, empirically setting the asym_metric_thresh ⁇ 0.4 and the peak_metric_thresh ⁇ 3.5 excludes most of the 400Hz false positives and still includes the majority of the genuine secondary signal cluster true positives as shown by Figure 12 .
- Figure 13 Three further scenarios are illustrated in Figure 13 , as described in the table below: Primary 400Hz Mains SCT absent SCT present Absent Figure 13(a) Figure 13(b) There are a few isolated points, well below the proposed thresholds. No false positives. Most of the SCT true positives lie in the top right quadrant described by the two proposed thresholds. Present Figure 13(c) Figure 12 There is a cluster of points, exceeding the peak threshold, but not the asymmetry threshold. This shows the value of the proposed dual threshold idea. No false positives would be generated. As discussed above.
- Such 'feature space classifications' may be provided to a user for system analysis, or the SCT determination may be performed directly on the data with no graphical output.
- Figure 14 shows a high-level flow diagram for the 'mixed domain' method; many of the steps having corresponding steps in the frequency domain SCT detection method. The detail relating to the corresponding steps described above applies to this alternative embodiment unless explicitly indicated otherwise.
- the first steps are as described previously, wherein the incoming signal is decimated 300 and 'chopped up' into overlapping windows 302.
- the method then branches, with one branch performing an FFT 500, estimating the frequency 502 and phase 504 of the primary transmission in the signal.
- the phase may be estimated by determining the phase of the samples used in determining the peak (e.g. the highest magnitude sample and the two either side). The highest magnitude samples would most likely be from the primary carrier so are most likely to have the primary phase.
- the primary carrier frequency and phase are used to down-convert 505 the time-domain windows by mixing each window with a complex sinusoid with the same frequency and phase-offset as the primary carrier transmission.
- the signal can be illustrated by Figure 15 where the in-phase (I) and quadrature (Q) components of a frequency down-converted signal (x'(t)) are plotted. If there were only perfect, phase noise free, primary carrier transmission, this vector would lie at constant ⁇ with its magnitude (i.e. length) changing with time. If there are any additive signals (such as SCTs or phase noise), the vector's angle would also change.
- the signal is phase-rotated by ⁇ and the part of the vector moving along the Q axis is measured.
- This step corresponds to the 'Quadrature split' 506 step in Figure 14 .
- This process is mathematically linear and so information is preserved and no artificial intermodulation effects are propagated through to the following processing steps.
- a real-only input FFT 508 is performed on the Q component of the phase-rotated signal. This provides a spectrum from which peaks are detected 510 corresponding to the out-of-phase components of the original signal.
- the specification above is primarily concerned with simultaneous voice transmissions received by an Air Traffic Controller, but it will be appreciated that the signal does not necessarily have to be voice transmissions.
- the signal may be digital information encoded into an AM radio transmission.
- the conjugate of the negative frequency sideband is subtracted from the related positive frequency sideband of the sum-signal so as to cancel out the primary carrier.
- the opposite operation is equally possible whereby the conjugate of the positive frequency sideband is subtracted from the related negative frequency sideband of the sum-signal.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Discrete Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Power Engineering (AREA)
- Algebra (AREA)
- Noise Elimination (AREA)
- Mobile Radio Communication Systems (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Interface Circuits In Exchanges (AREA)
- Rehabilitation Tools (AREA)
Description
- This invention relates to an apparatus and method for the detection of simultaneous call transmissions, in particular to double-sideband amplitude modulation (DSB-AM) transmissions.
- Simultaneous Call Transmission (SCT) is a situation when two or more transmissions occur simultaneously on the same frequency band. The end listener is then (usually) only able to understand the higher powered of the pair of transmissions. An example is illustrated in
Figure 1 . - This is a potentially hazardous situation as the sender of the weaker transmission may assume that they were actually heard by the end listener and take action accordingly. A situation where the consequent action would be incredibly dangerous would be where two planes transmit a signal to an air-traffic controller simultaneously, who then replies, and both planes believe the response is directed to them. Such a scenario may be noticed by an alert, trained human operator listening out for the characteristic phenomena of interfering voice and heterodyne tones arising from the frequency difference between the two transmitters. However this only reliably works when the weaker transmitter is (1) within the power range of the stronger transmitter e.g. 0 to -20dB and (2) the heterodyne lies in the filtered audio range (e.g. 300Hz to 3.5 kHz). The first condition may not be met if one transmission originates much further away than the other (i.e. one plane is overhead and the other is several kilometres away); the second condition may not be met if the frequency difference between the two transmitters (defined by the precision of the quartz in the transmitting equipment) is minor.
- Indeed, SCT can occur far outside of these values due to real world effects such as propagation loss, multipath error and frequency error.
- Hence automatic SCT detection is a desirable feature for radios.
-
WO2015/008165A2 describes a method for the detection of more than one signals contained in a received signal. The method comprises: down-converting the received signal, thereby providing a down-converted signal in a complex IQ base band; and at least partially cancelling the strongest user in the down-converted signal, thereby allowing for the detection of a possible secondary user. -
US2010/0067570A1 andDE102007037105A1 describe an apparatus adapted to automatically detect when two transmissions occur simultaneously. This system phase demodulates the sum-signal of the dual transmission by converting the baseband In-Phase and Quadrature signal into unwrapped phase and amplitude. A periodic 'wobble' is present on the unwrapped phase if there are close frequency-separated transmissions, as the difference in frequency causes the phasors of the transmissions to rotate around one another. This phase time-series is then Fourier transformed using a bank of transformers with varying window lengths to determine whether any secondary transmission (i.e. a peak due to the unwrapped phase 'wobble') is present. A warning tone is added to the audio output if a secondary transmission is detected so as to alert the operator to the situation. - This proposed solution has several significant drawbacks. The step of producing a phase time series is an intrinsically non-linear process (as it involves performing an arctan) and so errors propagate as the process progresses manifesting as intermodulation products in the output spectrum which are not physically present in the input. Furthermore, in real-world conditions this solution can potentially generate large occurrences of 'false positives' where an alert is sounded when only one transmission is present. This is because this system has no way of suppressing common signal impairments such as sinusoidal mains hum, incidental FM (frequency modulation), frequency-selective multipath effects and 1/f2 (reciprocal frequency-squared) phase noise. All of these effects can potentially introduce a certain amount of unwrapped phase 'wobble' and then be identified as secondary transmissions.
- False positives are very damaging because they cause an operator to lose faith in the equipment's reliability if it is 'crying wolf' too often. This may result in the operator turning off the automatic SCT feature completely or taking unnecessary mitigating action such as repetition of instructions. Because Air Traffic Control is a safety critical activity, an SCT detection system should be highly tolerant to real signal imperfections so that it is of the highest achievable reliability.
- On the other hand, false negatives are inevitable when the secondary transmission is very weak in power, when it becomes indiscernible from the noise floor, and also in the situation where the secondary transmission is superposed on the stronger transmission with negligible frequency difference.
- An improved solution is therefore needed.
- The invention is defined by the features of
independent claims 1 and 14. Further aspects of the invention are the subject of the dependent claims. - In one embodiment there is provided a method of determining the presence of a secondary transmission in a time-domain sum-signal including a primary transmission and the secondary transmission, the method comprising: transforming the time-domain sum-signal into a frequency domain sum-signal; wherein the frequency-domain sum-signal is a linear combination of the primary transmission and the secondary transmission; and wherein the transforming is based on a plurality of frequency bins; estimating a primary carrier frequency based on the frequency-domain sum-signal; shifting the frequency bins of the frequency-domain sum-signal based on the estimated primary carrier frequency; determining a conjugate of a sideband of the frequency-bin-shifted frequency-domain sum-signal; attenuating the frequency-bin-shifted frequency-domain sum-signal by using said conjugate of the sideband of the frequency-bin-shifted frequency-domain sum-signal; extracting at least one peak from the attenuated frequency-domain sum-signal; and determining the presence of the secondary transmission in the frequency-domain sum-signal based on said at least one peak.
- In another embodiment there is provided an apparatus for determining the presence of a secondary transmission in a time-domain sum-signal including a primary transmission and the secondary transmission, the apparatus comprising: means for transforming the time-domain sum-signal into a frequency domain sum-signal, wherein the frequency-domain sum-signal is a linear combination of the primary transmission and the secondary transmission; and wherein the transforming is based on a plurality of frequency bins; means for estimating a primary carrier frequency based on the frequency-domain sum-signal; means for shifting the frequency bins of the frequency-domain sum-signal based on the estimated primary carrier frequency; means for determining a conjugate of a sideband of the frequency-bin-shifted frequency-domain sum-signal; means for attenuating the frequency-bin-shifted frequency-domain sum-signal by using said conjugate of the sideband of the frequency-bin-shifted frequency-domain sum-signal; means for extracting at least one peak from the attenuated frequency-domain sum-signal; and means for determining the presence of the secondary transmission in the frequency-domain sum-signal based on said at least one peak.
- The invention extends to any novel aspects or features described and/or illustrated herein.
- Any feature in one aspect of the invention may be applied to other aspects of the invention, in any appropriate combination. In particular, method aspects may be applied to apparatus aspects, and vice versa.
- Furthermore, features implemented in hardware may be implemented in software, and vice versa. Any reference to software and hardware features herein should be construed accordingly.
- The invention also provides a computer program and a computer program product comprising software code adapted, when executed on a data processing apparatus, to perform any of the methods described herein, including any or all of their component steps.
- The invention also provides a computer program and a computer program product comprising software code which, when executed on a data processing apparatus, comprises any of the apparatus features described herein.
- The invention also provides a computer program and a computer program product having an operating system which supports a computer program for carrying out any of the methods described herein and/or for embodying any of the apparatus features described herein.
- The invention also provides a computer readable medium having stored thereon the computer program as aforesaid.
- The invention also provides a signal carrying the computer program as aforesaid, and a method of transmitting such a signal.
- Any apparatus feature as described herein may also be provided as a method feature, and vice versa. As used herein, means plus function features may be expressed alternatively in terms of their corresponding structure, such as a suitably programmed processor and associated memory.
- It should also be appreciated that particular combinations of the various features described and defined in any aspects of the invention can be implemented and/or supplied and/or used independently.
- In this specification the word 'or' can be interpreted in the exclusive or inclusive sense unless stated otherwise.
- The invention extends to methods and/or apparatus substantially as herein described with reference to the accompanying drawings.
- Purely by way of example, the present invention is illustrated by the accompanying drawings in which:
-
Figure 1 shows a Simultaneous Call Transmission (SCT) scenario; -
Figure 2 shows an example Double Side-Band Amplitude Modulated (DSB AM) signal; -
Figure 3 is an example flow diagram of an SCT detection method; -
Figure 4 is a schematic diagram of an apparatus operable to perform the method shown inFigure 3 ; -
Figure 5 shows overlapping windows of the sliding window buffer ofFigure 3 ; -
Figure 6 is an illustration of the 'zero padded FFT' ofFigure 3 ; -
Figure 7 shows example filters for use in an SCT detection method; -
Figure 8 shows the effect of the filters ofFigure 7 on an example time-domain equivalent (FFT input) envelope; -
Figure 9(a) shows an example frequency plot of a secondary signal, superposed on a primary signal carrying voice and 400Hz mains interference present; -
Figure 9(b) shows the signal ofFigure 9(a) following down-conversion; -
Figure 10(a) shows the left hand side of the signal ofFigure 9(b) reflected onto the right hand side; -
Figure 10(b) shows the signal ofFigure 10(a) following DSB-AM cancellation and a noise floor estimate; -
Figure 11(a) shows the signal ofFigure 10(b) following noise floor estimation; -
Figure 11(b) shows the asymmetry metric of the peaks detected from the signal ofFigure 10(a) ; -
Figure 12 shows a 'feature space' plot of the peak and asymmetry metrics of a large number of simulations; -
Figure 13(a) shows a scenario where an SCT and 400Hz mains noise are absent; -
Figure 13(b) shows a scenario where an SCT is present and 400Hz mains noise is absent; -
Figure 13(c) shows a scenario where an SCT and 400Hz mains noise are present; -
Figure 14 shows a flow diagram of an alternative method for SCT detection; and -
Figure 15 is an illustration of the in-phase and quadrature components of a signal where an SCT is present. - In this specification, the term 'primary transmission' or 'primary carrier' refers to the transmission with greatest power. The term 'secondary transmission' or 'secondary carrier' refers to any other (lower-powered) transmission occurring at the same time from another aircraft.
- Air traffic controller (ATC) to aircraft communication is generally very concise, each transmission is typically less than 10 seconds in duration, and can be as short as 2 or 3 seconds. For this reason, the latency in detecting a secondary call transmission (SCT) is preferably less than 2-3 seconds for this field of use. The term 'simultaneous' in this specification refers to the situation where two transmissions overlap in time, as in this scenario an ATC would either not hear the SCT or the audio would be filtered out by the radio.
-
- t is the time in seconds
- A is a gain constant (proportional to the transmitter Root Mean Square (RMS) power)
- v(t) is the real-valued audio signal, normalised to (-1,+1) peak to peak
- k is the modulation depth in the range (0,1) expressed as a percentage.
- ω c is the carrier frequency in radians/sec, typically approximately (2π)118MHz.
- θ is some notional phase offset (in radians) with respect to t=0
- j is √-1
- The spectrum (i.e. Fourier transform) of such a signal is shown in
Figure 2 - X(ω). The signal is centred on a theoretically infinitesimal carrier frequency ω c. In reality, this carrier band is broadened by system and transmission imperfections. There are 'conjugate symmetric' sidebands either side of ω c, such that, for within the bandwidth of v(t), X(ω c + ω) is the same as X(ω c - ω); equal magnitude and conjugate phase. -
-
-
- In order to perform such a calculation, properties of the primary carrier signal (i.e. ω c and c) must be first known. Once the primary carrier has been identified, it can be isolated and removed, allowing subsequent analysis of Y(ω) to determine whether an SCT has occurred.
-
Figure 3 illustrates a high-level dataflow of an exemplary method (referred to as 'Frequency Domain (FD) SCT detection') for detecting the presence of a secondary call transmission in an input signal. The input is a high sample rate IQ (in-phase / quadrature) baseband time-series (real or complex) and the outputs are SCT detection results (for example, as a tone inserted into the audio signal or a flag placed into a data stream). Each step is described in more detail below the following brief overview. - As used herein, the term "sum-signal" (and similar) preferably connotes a signal received by a receiver (for example an air traffic control system). Such a sum signal preferably comprises a primary carrier signal, noise, and potentially a secondary carrier signal (where present); the term "sum" merely implies, preferably, that such elements (primary carrier signal, noise, etc) are present in one and the same signal.
- A sum signal is received and converted to an IQ baseband sample. This signal is decimated 300 (downsampled) removing unnecessary frequency components so as to reduce the load on later processing steps.
- The decimated signal is sampled using an overlapping, sliding
window buffer 302 which stores a given length of signal to be processed. Sampling the signal means that the entire signal does not need to be processed in one go, improving the latency of detection and reducing processor load. - The length of the buffer is determined by the trade-off between processor load in analysing lengthy samples and the increased detection accuracy (i.e. high signal-to-noise ratio) longer samples afford. The sampling rate (i.e. number of windows per second) is determined by a trade-off between processor load in processing large numbers of windows within a short time-period and latency of SCT detection.
- Each window sample is then input into a Fast Fourier Transform (FFT) 304 - outputting X(ω) - for further processing in the frequency domain. An FFT is preferable to a continuous Fourier Transform (FT) as it is much less processor-intensive. Other discrete transforms such as wavelet transforms or spectral line filters may be used.
- The decimation step may introduce 'decimator ripple' in the frequency domain output, this is corrected for at
step 306 before any further processing. - The frequency of the primary carrier transmission is estimated 308, and the FFT output is down-converted 310 based on this frequency. The in-phase elements to the primary carrier transmission are then cancelled by subtracting the conjugate of the negative sideband frequencies of the primary signal from their counterpart positive sideband frequencies in the (decimated, down-sampled, frequency domain)
sum signal 312. - The remaining signal - Y(ω) - is due to phase noise or other signals at different frequencies (which may be SCTs) to the main carrier transmission in the signal. This signal is analysed for peaks 316 (defined by a threshold) above an estimated
noise floor 314. These peaks are indicative of an SCT being present as they represent significant magnitude parts of the signal which are not on the same carrier frequency as the primary transmission (e.g. peaks corresponding to a heterodyne tone). However, other effects, such as mains hum and phase noise may show up as peaks above a nominal noise floor. Noise effects such as these affecting the sidebands of the primary carrier are typically symmetric about the primary carrier frequency, thus anasymmetry analysis 318 is performed to determine whether a particular peak has a corresponding 'mirror image' peak. This analysis is performed using an 'asymmetry threshold'. The peak is also analysed for its magnitude (above the noise floor) as higher power peaks are more likely to be secondary transmissions rather than variations in noise. These two parameters (and/or others) are combined in a 'Feature Space Classification' 320 and an SCT can be signalled if a sample contains a peak exceeding the predetermined threshold(s). - This process is performed continuously for every sampled window of the incoming signal. It should be noted that the method has been shown as split into numerous discrete steps whereas in practice many of these steps may occur simultaneously or as part of a single step.
- The process described above (after sampling) is undertaken entirely on the spectrum of time windows of signal rather than any time-series. This is significant as the main (distinguishable) physical difference between primary and secondary transmissions is the slight difference in central frequency. Analysing the spectrum is thus addressing the fundamental problem. Furthermore, all of the above processes are mathematically linear, thus greatly reducing the potential of spurious artefacts being introduced, or for errors to propagate and become amplified as the process continues.
-
Figure 4 shows a schematic diagram of aradio receiver 104 adapted to perform the processes involved in the detection of a secondary call transmission as described above. - A signal is received by an aerial and is input to an Analogue-to-Digital Converter (ADC) 402. The dotted
section 400 represents a simplified digital radio receiver without any SCT detection capability. The digital signal is demodulated by demodulatingunit 403, with assistance from acentral processor 422 andmemory 424. This is then passed to anaudio output unit 420 and audio is outputted. An actual digital radio may include many additional components (such as tuning, filtering and amplifying circuitry), but such components are omitted for clarity in this figure. - This audio extraction process occurs independently of the SCT detection, as this represents the primary purpose of the
radio 104, to convert received signals into audio (or other useful information). The components used for SCT detection are shown outside of the dottedsection 400. - The digital signal is decimated (downsampled) by
decimator 404 before being sampled in a slidingwindow buffer 406. Each window is then passed through a Fast Fourier Transform (FFT) 408. - The spectrum outputted from the FFT has filters /
windows 410 applied to it so as to produce the signal defined by Y(ω) inequation 4. This output is passed to acomparator 412, which, withlogic circuitry 412 and thresholds stored inmemory 424, determines whether an SCT has occurred. If so, the operator is notified, for example by a tone being inserted into the audio output and/or a flag (such as an indicator on a user interface) is raised via tone /flag generator 418. Other information regarding the SCT, such as indication of a confidence level, or timestamp of the event, may also be outputted. -
Figure 4 shows components separated for clarity whereas in reality many of these components may be combined as a single component (such as the comparator, logic combined with the processor) or further split into separate components. - The following description further details the various steps briefly described above.
- In a DSB-AM radio receiver, the intermediate digital signal after analogue to digital conversion is often at a higher sample rate than is required to support the DSB-AM sidebands of a primary signal.
- The
decimation stage 300 represents the context-dependent low-pass filtering and down-sampling that may be required to reduce the signal bandwidth to a value of e.g. ±7kHz at a sample rate of fs =14kHz. For instance, such a bandwidth will support a primary DSB-AM signal with 4kHz audio bandwidth (A) and up to ±3kHz of frequency error (B). The relationship between A, B and fs are illustrated below, with additional example values:Audio bandwidth on Primary DSB signal (kHz) Maximum Primary Frequency Error (kHz) Complex Sample Rate after decimation (fs ) (kHz) A ±B fs=2*(A+B) 5 4 18 3 3 12 4 3 14 - The purpose of decimation is twofold: (1) to reduce the computational load and (2) reject signals outside the band of interest for SCT detection. The decimation step (and the subsequent ripple equalisation) would not be necessary if these issues are not of any relevance (e.g. if the analogue-to-digital converter has a low sample rate).
- The decimator design preferably has a narrow transition region (for example 10% of passband) with low passband ripple (0-3dB) and high stopband attenuation (for example more than 40dB) i.e. a typical specification for a high-quality decimator for audio applications. A typical low-pass mask specification given the third set of parameters in the above table would be ±1dB of passband ripple up to 5kHz, a transition region from 5kHz to 7kHz, and -60dB of gain in the decimation stopband.
- Note that if the time-series is real-only, then a complex oscillator and mixer are required to down-convert the signal before the decimator. The decimation low-pass filter then requires sufficient stopband attenuation to adequately remove the frequency-shifted conjugate image.
- This stage presents the most recent T seconds block of sample data M times a second to the subsequent processing stages as illustrated in
Figure 5 . Blocks preferably have a high degree of overlap to maximise the chance of detecting a secondary transmission as soon as it starts. Typical values for this stage (for the ATC example discussed above) are T=2 seconds, M=4Hz, in this case the longest time possible between an SCT event beginning and the end of one window - tmax is 1/M = 0.25 seconds. However, not all SCT events may be detected in such a short time period as the signal may be too weak; in such circumstances, due to the overlap of windows, the next window would have 0.5 seconds of SCT to detect and so on. Thus a detectable constant SCT would be detected within T + tmax (= T + 1/M) seconds. In order to ensure overlap the following identity must hold: T*M>1, but ideally at around 8 windows overlap, so T*M≥8. - The purpose is to allow strong signals to be detected quickly by the system at a coarse granularity in time, but also allow an adequate time-history to allow the coherent integration and detection of weak secondary signals.
- In use, M times a second, a buffer comprising the last T seconds of data is processed. This results in each data block of
length 1/M seconds being processed T*M times in total. For illustration T can be in therange 1 to 4 seconds and M can be in therange 2 to 16. T governs the coherent integration period for detecting weak SCT signals and it is of advantage to be long and about the same length as typical primary transmission utterances. The value of (1/M) governs the maximum latency for detecting strong SCT signals and it is advantageous for M to be high for low latency. The processing load is proportional to the product T*M, thus there is a trade-off between performance and processor load when selecting values of T and M. - Although in use a large number of window buffers will be processed, for clarity, the following description will solely focus on the processing of a single window.
-
Figure 6 illustrates how the current analysis window from the buffer is mapped into the FFT input with 'zero padding'. The mapping is unconventional in that the buffer is split into two halves and the first half is mapped into the final part of the FFT input and the second half is mapped into the start of the FFT input, with zero-value inputs occupying the intervening samples. This improves the operation of frequency domain down-conversion as described below. - The FFT size N FFT is chosen to be around twice the size of the buffer window to provide sufficient oversampling for subsequent processing. The choice of oversampling ratio at approximately x 2 is a compromise between two conflicting factors: (1) critical sampling at around x 1 oversampling is not viable because down-conversion requires an unfeasibly long resampling filter for the required DSB-AM cancellation fidelity, (2) the system performance at say, > x 3 oversampling yields negligible performance benefit at the expense of increased computational complexity in the FFT. Of course, oversampling ratios of greater than
x 3 may be used if computational complexity is not an issue - for example if the fidelity of the cancelation is paramount. - For example with a signal sampling rate, fs , of 14kHz and T=2 seconds, the buffer is 28,000 samples long. This indicates that an FFT size of N FFT=65,536 (the oversampling ratio being 2.34) is appropriate; presuming that a standard Digital Signal Processing (DSP) library function is used requiring a power-of-two size (i.e. N FFT = 2n, where
Sample Rate (kHz) T (seconds) Buffer Blocksize (kSamples) Radix-2 FFT Size (kSamples) Oversampling Ratio fs T C = fs*T N FFT N FFT / C 10 2 20 64 (≈216) 3.2 12 4 48 64 (≈216) 1.3 14 6 84 128 (≈217) 1.5 16 8 96 128 (≈217) 1.3 - Depending on operational requirements / constraints, a larger or smaller oversampling ratio may be used. The larger the oversampling, the more processor-intensive the resulting analysis would be (due to the greater number of discrete frequency 'bins' in the frequency domain) but the system would be more accurate due to (at least) the spectrum having greater resolution.
- For later convenience, the FFT output vector is denoted as the vector x with elements x i where i={0,1...,N FFT-1} counting up from the zero frequency bin.
- The low-pass filter discussed above with reference to decimation may have significant passband ripple in order to be implementable with realistic cost. Passband ripple is an artefact manifesting in the spectrum of a transformed signal having had imperfect (i.e. non-square) band-pass filters applied to it.
- The gain fluctuation across the band of interest can degrade the ability to perform primary carrier double side-band cancellation as it affects the conjugate symmetry property exploited in
equation 3. A low-cost and simple way to compensate this effect is to calculate the ripple across the decimated band H(ω) from the FFT of the impulse response caused by the decimation, and apply gain and phase compensation to the output of the FFT of 1/H(ω). - The
inverse transform 1/H(ω) is stored as a vector of N FFT complex weights which is applied to the output directly after the FFT has been computed. - Although H(ω) is symmetric about zero hertz, it is not symmetric about the primary carrier, so would not be cancelled out when calculating Y(ω) - which is described in more detail below.
- The highest magnitude FFT output bin (denoted as bin j) is detected and its power and frequency are measured. This is asserted to be the primary carrier (i.e. strongest sinusoidal tone) and these measurements are passed on to the classification stage discussed in order to detect if any primary signal is present. Identifying the primary carrier frequency leads to identification of non-primary carrier signals (such as an SCT).
- Taking the magnitude samples of the three FFT output bins {j-1, j, j+1} a parabolic (quadratic) curve may be fitted to the points, for example using closed-form linear algebra. The fractional bin frequency f in the range of -0.5 to 0.5 of the maximum value of the fitted parabola is taken to be the best estimate of the true primary carrier frequency ωc. The oversampling of the FFT (e.g. twice oversampling) provides an interpolated mainlobe of the primary carrier and thus facilitates an accurate peak position estimate. Accurate primary carrier peak estimation allows for a more accurate down-conversion, leading to improved subsequent DSB-AM cancellation as the centre point of the reflection is more accurate.
- At this stage the width of the primary mainlobe may be assessed by searching out from the peak in both negative and positive frequency until bins that are <3dB (approximately <0.5 in power) of the peak are identified (i.e. full-width, half maximum (FWHM) of the primary mainlobe). Leading and trailing edges of the primary transmitter in the analysis window cause wide mainlobes, and this measurement may be useful in the later 'feature classification' stage for assessing the time-domain activity of a primary transmitter.
- Frequency-domain down-
conversion 310 is performed by generating a finite-impulse filter which shifts the frequency bins by -(j+f) bins (i.e. by ωc) so that the underlying maximum of the primary carrier mainlobe is shifted exactly on to the zero frequency bin. This step effectively makes the primary carrier signal symmetric about zero hertz, making later computation and determination of SCT events simpler. - The formula for the filter is given in
Equation 5 where
Ncoeffs= 4, w=[0.3635819, -0.4891775, 0.1365995, -0.0106411]
w generates a Blackman-Nuttall window which has good sidelobe performance. Other windows may be used such as 'Kaiser' or 'Equiripple' windows, but cosine-family windows such as Hamming, Hann, Blackman family have the implementation benefit of combining good sidelobe performance with the precise and simple computation using cosines. - The value xLIM sets the limits for the window (i.e. it is zero-valued for |x| > xLIM ), and hence defines the quality of the resampling (a typical value would be xLIM=5). A small value is desirable in order to minimise the processing complexity of down-conversion. The choice of xLIM is discussed in more detail below with reference to
Figures 7 and 8 . - The components of this frequency domain down-conversion filter are shown in
Figure 7 . This figure also shows the 'half bin samples' which may be present due to the FFT oversampling described above with reference toFigure 6 . The effect of these fractional bins are described in more detail below with reference toFigure 8 . - Down-conversion is performed using Equation 6 as a circular convolution on the FFT output x but only including the non-zero terms from
Equation 5 to minimise computational cost. For example, with xLIM =5, only 2xLIM +1 =11 (-5 to +5) multiply / accumulates are needed per bin. This is analogous to implementing a short Finite Impulse Response (FIR) filter. - Frequency domain convolution of two signals is analogous to multiplication of their time-domain equivalents. In this case, the inverse Fourier transform of the g(x) term in Equation 6 is an arbitrary frequency sinusoid with a sampling-dependent envelope function: this is unity when samples are taken on a grid at integer values of x (i.e. f=0) and (worst-case) has ramping and a zero point when samples are taken on a grid at halfway between FFT bins (i.e. f=±0.5) as illustrated in
Figure 8 . Other values of f create envelopes intermediate between these extremes. Although this process would not be necessary if only integer values of f were used, doing so would introduce errors into the central frequency and thus mean that the later asymmetry analysis would carry through these errors. - The mathematical explanation for the envelope phenomenon shown in
Figure 8 is as follows.Equation 5 comprises the product of two terms; (1) a sin(x)/x function with infinite support on x (which has too many terms to compute practically) and (2) a compactly supported window function (which makes g(x) economic to compute). In time domain, by analogy, this is the circular convolution of (1) an arbitrary frequency sinusoid and (2) a band-pass filter corresponding to the frequency shifted IFFT of the window function. The output of this filtering process is unit amplitude sinusoid except where a phase discontinuity passes through the filter where the two ends of the sinusoid are circularly "spliced" together. This creates the characteristic "dip" in the sinusoid envelope illustrated inFigure 8 which is worst-case when a 180 degree continuity passes through (as occurs with the half-bin case). -
Figure 8 also explains the utility of the unconventional zero-padding described above of mapping the "1st half" and "2nd half' of the input buffer to time-domain intervals where the envelope function is almost exactly unity. The mapping of the second half of the time window to the first part of the FFT input and vice versa means that the FFT input maintains its time-order as the end of the first half is effectively contiguous with the start of the second half (as the FFT can be visualised as wrapped around the surface of a cylinder). Hence if the equivalent time-domain product is taken (by taking notional IFFTs of the frequency domain convolution), we have the desired effect of the signal multiplied by a unit-amplitude complex sinusoid in order effect high-quality, precise down-conversion of the primary signal. Slight deviation from unity over the non-zero-padded part of the envelope function is permitted as a perfect standard rectangular window is not necessary. A tolerance of deviation from the maximum of the envelope function of approximately 1% is preferable. - The choice of xLIM is a function of the oversampling ratio
-
- DSB-
AM cancellation 312 as discussed above is effected by applying Equation 8 in order to generate an output vector y comprising NFFT /2+1 bins (the zero frequency bin and the right hand side of the spectrum). As Y(ω) is by mathematical definition conjugate symmetric about zero for an ideal primary carrier, only computation of the right hand side (i.e. positive frequency) is necessary. Only magnitude information is taken into y for the purpose of peak detection, hence the modulus is taken. - The quality of DSB-
AM cancellation 312 is dependent on the temporal coherence of the primary signal. Phase noise on the primary carrier can lead to some feed through of tonal components in the sidebands which may appear as distinct tones in Y(ω). A simple technique for identifying such tones using the concept of power 'asymmetry' is described below. - Subtracting the conjugate of the negative frequencies from the positive frequencies of the sum signal (after down-conversion) effectively cancels out the part of the signal in-phase with the primary carrier (attenuating the frequency domain primary carrier within the sum-signal), leaving just signals which have introduced phase noise into the sum-signal. These signals include phase noise (which would generally be at a low-level across a wide range of frequencies) and specific tones, which would manifest as peaks in the frequency plot.
- SCT tones in y are characterised by isolated narrowband peaks against a noise floor after DSB-AM cancellation of the primary carrier transmission. Hence, in order to detect peaks, a noise floor estimate which is not biased by tonal peaks should be estimated. Noise levels may not be constant over the whole frequency range in question, so the noise level at every frequency bin is estimated in order to 1) capture secondary transmissions above the local noise level, but potentially below the noise level elsewhere, and 2) discount frequency bins with higher levels of noise than elsewhere. A single estimate of the noise level across the entire frequency spectrum would not be able to account for such circumstances, resulting in, in the case of 1) false negatives, and in the case of 2) false positives. Either of these scenarios is undesirable, false negatives particularly so in an ATC implementation as such events could result in a dangerous situation.
- An effective way to determine a frequency-dependent noise floor estimation is to calculate a moving-average of the magnitude across a range of bins centred around a particular frequency bin. If a large enough bin range is used and peaks are not frequent, this would be an accurate representation of the noise floor at that frequency bin. In one example a short sliding-window rank-order statistic filter is applied which extracts the e.g. the median, power bin as the noise floor estimate. Analogous filters are used for removing impulsive noise from otherwise smooth functions in applications like image processing.
- This concept is expressed simply in Equation 9 where the median window estimate is over ±NNFE bins (a typical value is NNFE =256 when NFFT =65536). If the window is too long, frequency-dependent changes in the noise floor are smoothed-out and the noise floor does not respond to local effects such as colouration from filters. On the other hand if NNFE is too short, legitimate SCT peaks may adversely bias the noise floor estimate leading to them being smoothed and subsequently discounted. A value for NNFE of approximately
-
- Median filtering is costly to compute in terms of procession time and power. A practical optimisation is to decimate y by summing contiguous blocks of D 1 samples and then using a much shorter median filter over ±D2 on the resulting decimated signal. The aggregate window size is NNFE = D 1 D 2 . For instance, D 1=16 and D 2=16 when NNFE =256. This has little performance loss when the time series y is dominated by the noise floor and has sparsely located peaks. In one
embodiment - The median is the default rank-order statistic to draw out, but other measures of central tendency are possible, for example the 40th centile, which will be less biased by peaks, but more susceptible to low power noise samples.
- Peaks are identified in y by identifying local maxima, where yi>yi-1 and yi>yi+1. Performing just this analysis may pick up a lot of spurious fluctuations in the noise floor, for this reason only peaks (i.e. values of yi) that satisfy a certain predefined threshold (peak_metric_thresh) are identified as SCT candidates. Example values for peak_metric_thresh are provided below with reference to
Figures 11 to 14 , but may vary from around 0.85 to 3 (or greater than 3) depending on the situation. -
-
- The threshold peak_metric_thresh is preferably a system-set parameter which may be calculated once upon calibration of the system; alternatively it may be dynamically calculated so as to result in a system with a specific false-positive rate. This may be useful if the variance in the noise floor (i.e. the accuracy of the noise floor estimate) changes over time so that the system becomes more prone to false negatives (if the variance decreases) or it becomes more prone to false-positives (the variance increases). In an average situation, a value for the peak metric threshold p(i) would be between 1 and 4, more preferably between 2 and 3 as a general illustration in these circumstances.
- Another metric that may be used to reduce the number of candidate peaks is to specify that two peaks must be separated by a minimum frequency otherwise they are treated as a single peak (i.e. the smaller peak is disregarded). The threshold min_freq_sep is defined. In one example this is between 5Hz and 50Hz, preferably between 7Hz and 15Hz, and preferably approximately 10Hz. Disregarding the smaller peak of a closely separated pair of peaks has negligible impact on the capability to detect genuine secondary tones when peak detections are sparsely separated. Such a feature allows strong peaks from e.g. 400Hz mains hum (which are highly conjugate-symmetric) to absorb their own sidelobe features which are much weaker in power but more asymmetric and thus can cause false positives. The method identifies the weaker peaks from the set P which are within +/- min_freq_sep of the current secondary tone candidate being analysed, and marks them for deletion from set P by placing them in the set Q as follows (with commentary accompanying each step):
-
- An asymmetry analysis favours 'legitimate' SCT events over other phase noise as SCT events have (by definition) a central frequency offset from that of the primary carrier and are thus asymmetric about primary carrier (and, after down-conversion, are asymmetric about zero hertz). There is a low probability of another tone precisely at the opposite frequency sign as this would correspond to a third SCT at a very specific frequency.
- In contrast, "worst-case" primary signals with the deleterious properties of (1) high phase noise and (2) voice sidebands contaminated with interference tones (from e.g. mains electricity) generate secondary carrier candidates which are very symmetric in power (from the core definition of DSB-AM conjugate symmetry).
- Hence the asymmetry metric a(i) provides a useful way to exploit values that are pre-computed elsewhere in the process (i.e. bins from the down-converted X(ω) in vector x) to reject false positives from poor quality primary transmitters.
- A threshold, asym_metric_thresh, for the value of a(i) is defined where peaks not meeting this threshold are discarded as being too symmetrical, and thus unlikely to be SCTs. The asymmetry threshold provides a means to discount peaks which have a high residual power following subtraction due to the fact that the symmetric peaks had a high power prior to subtraction - for example if the signal has a high level of noise (which is not perfectly symmetrical), or due to external effects such as mains hum.
Figures 12 ,13(b) and13(c) below show scenarios where the asymmetry threshold is utilised to reduce the false-positive rate by limiting the number of events above the power threshold which would otherwise be deemed to be SCT events. - Before a candidate peak from set P can be determined as an SCT event, a number of checks may be performed.
- In order for an SCT to be present, there must be first the presence of a primary peak. This eliminates false positives when there is no transmission being received. A threshold primary_pk_thresh is defined where SCT analysis is only undertaken if the primary peak is above this threshold. This threshold is corrected by the amount of gain applied to the signal (AGC_gain) so as to measure the absolute power of the primary signal.
- A threshold is also set for the maximum allowed width of the primary peak, primary_bw_thresh, where SCT analysis is only undertaken if the width of the primary carrier peak is greater than this threshold. This ensures that a certain lower bound is met on the mark-space ratio of the primary transmitter in the analysis window, for example it may be desirable for the primary transmission to occupy at least 50% of the time window. This can prevent some anomalies due to rising edges entering or trailing edges leaving the analysis window. The width of the primary carrier peak is an output which is simple to generate and which provides some clear information about the temporal activity of the primary transmitter.
- The following section describes logic which may implement the classification part of the method.
- The following additional inputs are used for detecting the presence of a primary signal (and have associated thresholds):
- primary_pk
Magnitude value of primary peak - primary_bw
Primary peak 3dB width in bins (FWHM) - AGC_gain
- The Automatic Gain Control magnitude gain applied elsewhere in the receiver.
- Automatic Gain Control (AGC) will modulate the dynamic range of signals; therefore the primary_pk value is scaled by the amount of applied AGC and thus needs to be re-scaled by the reciprocal of the AGC gain in order to have an absolute power in terms of dBm.
- The following decision logic is given as an example of how to generate a Boolean detection output.
if (primary_pk>(primary_pk_thresh/AGC_gain)) AND (primary_pk>primary_bw_thresh) AND there exists any a(i)>asym_metric_thresh; i∈P then SCT_detect=TRUE else SCT_detect=FALSE end (if)
- Primary DSB-AM signal carrying voice audio and additive loud 400Hz mains hum
- Primary carrier frequency error
- Significant phase noise on the primary carrier
- Secondary DSB-AM signal carrying voice
- Additive White Gaussian Noise (AWGN)
Primary 400Hz Mains | SCT absent | SCT present |
Absent | | |
There are a few isolated points, well below the proposed thresholds. No false positives. | Most of the SCT true positives lie in the top right quadrant described by the two proposed thresholds. | |
Present | | |
There is a cluster of points, exceeding the peak threshold, but not the asymmetry threshold. This shows the value of the proposed dual threshold idea. No false positives would be generated. | As discussed above. |
Claims (14)
- A method of determining the presence of a secondary transmission (200) in a time-domain sum-signal including a primary transmission (100) and the secondary transmission (200), the method comprising:transforming the time-domain sum-signal into a frequency-domain sum-signal (304),wherein the frequency-domain sum-signal is a linear combination of the primary transmission and the secondary transmission; andwherein the transforming is based on a plurality of frequency bins;estimating a primary carrier frequency (308) based on the frequency-domain sum-signal;shifting the frequency bins of the frequency-domain sum-signal based on the estimated primary carrier frequency;determining a conjugate of a sideband of the frequency-bin-shifted frequency-domain sum-signal (312);attenuating the frequency-bin-shifted frequency-domain sum-signal by using said conjugate of the sideband of the frequency-bin-shifted frequency-domain sum-signal (312);extracting at least one peak from the attenuated frequency-domain sum-signal (316); anddetermining the presence of the secondary transmission in the frequency-domain sum-signal based on said at least one peak.
- A method according to claim 1, comprising shifting the frequency bins so that the primary transmission (100) is effectively symmetric about zero Hertz.
- A method according to any preceding claim, wherein transforming the time-domain sum-signal into the frequency domain sum-signal comprises at least one of:oversampling;oversampling using an oversampling ratio of approximately 2; andoversampling using an oversampling ratio greater than 3.
- A method according to any preceding claim, further comprising estimating a phase of the primary carrier.
- A method according to any preceding claim, including phase-rotating the frequency-domain sum-signal.
- A method according to any preceding claim, wherein determining the presence of the secondary transmission comprises performing a symmetry analysis (318) on said at least one peak.
- A method according to any preceding claim, wherein the time-domain sum-signal is decimated (300) so as to reduce the bandwidth.
- A method according to claim 7, further comprising correcting for decimator ripple.
- A method according to claim 8, wherein correcting for decimator ripple comprises gain-transforming (306) the frequency-domain sum-signal.
- A method according to claim 9, wherein the gain transforming uses the reciprocal of a gain due to a magnitude spectrum of the decimator.
- A method according to any preceding claim wherein the time-domain sum-signal is sampled (302).
- A method according to claim 11, wherein the time-domain sum-signal is sampled in overlapping blocks.
- A method according to any preceding claim, further comprising at least one of:alerting an operator to the presence of the secondary transmission;inserting a tone into an audio output; andindicating the presence of the secondary transmission on a user interface.
- An apparatus for determining the presence of a secondary transmission (200) in a time-domain sum-signal including a primary transmission (100) and the secondary transmission (200), the apparatus comprising:means for transforming the time-domain sum-signal into a frequency-domain sum signal (408),wherein the frequency-domain sum-signal is a linear combination of the primary transmission and the secondary transmission; andwherein the transforming is based on a plurality of frequency bins;means for estimating a primary carrier frequency (308) based on the frequency-domain sum-signal;means for shifting the frequency bins of the frequency-domain sum-signal based on the estimated primary carrier frequency;means for determining a conjugate of a sideband of the frequency-bin-shifted frequency-domain sum-signal (312);means for attenuating the frequency-bin-shifted frequency-domain sum-signal by using said conjugate of the sideband of the frequency-bin-shifted frequency-domain sum-signal (312);means for extracting at least one peak from the attenuated frequency-domain sum-signal; andmeans for determining (414) the presence of the secondary transmission in the frequency-domain sum-signal based on said at least one peak.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PL15721779T PL3123648T3 (en) | 2014-03-24 | 2015-03-24 | Simultaneous call transmission detection |
HRP20211153TT HRP20211153T1 (en) | 2014-03-24 | 2021-07-19 | Simultaneous call transmission detection |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB1405259.1A GB2522083B (en) | 2014-03-24 | 2014-03-24 | Simultaneous call transmission detection |
PCT/GB2015/050878 WO2015145138A2 (en) | 2014-03-24 | 2015-03-24 | Simultaneous call transmission detection |
Publications (2)
Publication Number | Publication Date |
---|---|
EP3123648A2 EP3123648A2 (en) | 2017-02-01 |
EP3123648B1 true EP3123648B1 (en) | 2021-05-05 |
Family
ID=50686822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP15721779.5A Active EP3123648B1 (en) | 2014-03-24 | 2015-03-24 | Simultaneous call transmission detection |
Country Status (11)
Country | Link |
---|---|
US (2) | US10263721B2 (en) |
EP (1) | EP3123648B1 (en) |
CN (1) | CN106664148B (en) |
AU (1) | AU2015237965B2 (en) |
DK (1) | DK3123648T3 (en) |
ES (1) | ES2872332T3 (en) |
GB (1) | GB2522083B (en) |
HR (1) | HRP20211153T1 (en) |
PL (1) | PL3123648T3 (en) |
PT (1) | PT3123648T (en) |
WO (1) | WO2015145138A2 (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9825794B2 (en) | 2013-06-25 | 2017-11-21 | Thuy Duong NGUYEN | Weak signal detection in double transmission |
GB2522083B (en) | 2014-03-24 | 2016-02-10 | Park Air Systems Ltd | Simultaneous call transmission detection |
US11539632B2 (en) * | 2014-12-30 | 2022-12-27 | Research Electronics International, Llc | System and method for detecting constant-datagram-rate network traffic indicative of an unmanned aerial vehicle |
US10938682B2 (en) * | 2014-12-30 | 2021-03-02 | Research Electronics International, Llc | System and method for detecting constant-datagram-rate network traffic |
CN108886374B (en) * | 2016-01-18 | 2021-08-03 | 唯亚威通讯技术有限公司 | Method and apparatus for detecting distortion or deformation of cellular communication signals |
DE102016205609A1 (en) * | 2016-04-05 | 2017-10-05 | Rohde & Schwarz Gmbh & Co. Kg | Processing device and method |
JP6647424B2 (en) * | 2016-11-16 | 2020-02-14 | 株式会社日立国際電気 | transceiver |
CN108665889B (en) * | 2018-04-20 | 2021-09-28 | 百度在线网络技术(北京)有限公司 | Voice signal endpoint detection method, device, equipment and storage medium |
AT522205B1 (en) | 2019-03-08 | 2021-05-15 | Frequentis Ag | Method for recognizing and reproducing voice radio messages emitted by a plurality of transmitters via radio, as well as a device for this purpose |
US11855683B2 (en) * | 2021-06-29 | 2023-12-26 | Raytheon Company | Configurable acquisition engine for receiver of spread spectrum signals |
CN114124631B (en) * | 2021-11-15 | 2023-10-27 | 四川九洲空管科技有限责任公司 | Processing method suitable for audio synchronous control between embedded equipment of aircraft cabin |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4031462A (en) * | 1975-07-07 | 1977-06-21 | Motorola, Inc. | Frequency spectrum analyzer |
US5008939A (en) | 1989-07-28 | 1991-04-16 | Bose Corporation | AM noise reducing |
US7127008B2 (en) | 2003-02-24 | 2006-10-24 | Ibiquity Digital Corporation | Coherent AM demodulator using a weighted LSB/USB sum for interference mitigation |
US7184485B2 (en) * | 2004-07-01 | 2007-02-27 | Texas Instruments Incorporated | Time-domain windowing of multi-band OFDM system to enable spectral sculpting |
US7307700B1 (en) * | 2004-12-17 | 2007-12-11 | The Boeing Company | Ultra-linear signal processing for radar and laser radar |
US7587171B2 (en) * | 2005-03-09 | 2009-09-08 | Atc Technologies, Llc | Reducing interference in a wireless communications signal in the frequency domain |
CN101455008B (en) * | 2006-04-03 | 2012-10-24 | 伟俄内克斯研究公司 | Frequency offset correction for an ultrawideband communication system |
CN101277281B (en) * | 2007-03-29 | 2015-05-20 | 深圳赛意法微电子有限公司 | Method and equipment for estimating channel equalization for receiver |
DE102007037105A1 (en) * | 2007-05-09 | 2008-11-13 | Rohde & Schwarz Gmbh & Co. Kg | Method and device for detecting simultaneous double transmission of AM signals |
US8155641B2 (en) * | 2007-12-27 | 2012-04-10 | Michael Hirsch | System and method for preventing lapses of communication in radio voice communications |
CN102792655B (en) * | 2010-01-22 | 2016-10-12 | 索尼公司 | OFDM transmission in multi-carrier data transmission systems and reception equipment, method and system |
CA2791005C (en) * | 2010-02-25 | 2015-09-29 | Mitsubishi Electric Corporation | Interference wave suppressing apparatus, relay apparatus, relay system, and interference wave suppressing method |
DE102010056528B4 (en) * | 2010-12-29 | 2016-08-18 | Atlas Elektronik Gmbh | Method and device for increasing the DF accuracy of a receiver arrangement |
EP2533069A1 (en) * | 2011-06-10 | 2012-12-12 | Sony Corporation | Signal processing unit and method |
DE102011080999A1 (en) | 2011-08-16 | 2013-02-21 | Rohde & Schwarz Gmbh & Co. Kg | Method and device for detecting radio signals transmitted simultaneously in the same channel |
US9520144B2 (en) * | 2012-03-23 | 2016-12-13 | Dolby Laboratories Licensing Corporation | Determining a harmonicity measure for voice processing |
CN102608663A (en) * | 2012-03-23 | 2012-07-25 | 西安电子科技大学 | Interference canceller applied to detecting core quadrupole moment resonance signal |
US8891603B2 (en) * | 2012-06-25 | 2014-11-18 | Tektronix, Inc. | Re-sampling S-parameters for serial data link analysis |
US9304189B2 (en) * | 2013-03-08 | 2016-04-05 | Qualcomm, Incorporated | Systems and methods for detecting radar signals |
DE102013212067A1 (en) * | 2013-06-25 | 2015-01-08 | Rohde & Schwarz Gmbh & Co. Kg | Measuring device and measuring method for the detection of simultaneous double transmissions |
US9825794B2 (en) | 2013-06-25 | 2017-11-21 | Thuy Duong NGUYEN | Weak signal detection in double transmission |
US9319156B2 (en) * | 2013-12-04 | 2016-04-19 | Aruba Networks, Inc. | Analyzing a particular wireless signal based on characteristics of other wireless signals |
GB2522083B (en) * | 2014-03-24 | 2016-02-10 | Park Air Systems Ltd | Simultaneous call transmission detection |
US9628122B1 (en) * | 2016-07-25 | 2017-04-18 | The Aerospace Corporation | Circuits and methods for reducing interference that spectrally overlaps a desired signal based on dynamic gain control and/or equalization |
-
2014
- 2014-03-24 GB GB1405259.1A patent/GB2522083B/en active Active
-
2015
- 2015-03-24 US US15/128,974 patent/US10263721B2/en active Active
- 2015-03-24 AU AU2015237965A patent/AU2015237965B2/en active Active
- 2015-03-24 PT PT157217795T patent/PT3123648T/en unknown
- 2015-03-24 WO PCT/GB2015/050878 patent/WO2015145138A2/en active Application Filing
- 2015-03-24 CN CN201580027508.3A patent/CN106664148B/en active Active
- 2015-03-24 PL PL15721779T patent/PL3123648T3/en unknown
- 2015-03-24 DK DK15721779.5T patent/DK3123648T3/en active
- 2015-03-24 EP EP15721779.5A patent/EP3123648B1/en active Active
- 2015-03-24 ES ES15721779T patent/ES2872332T3/en active Active
-
2019
- 2019-02-27 US US16/286,647 patent/US20190199462A1/en not_active Abandoned
-
2021
- 2021-07-19 HR HRP20211153TT patent/HRP20211153T1/en unknown
Also Published As
Publication number | Publication date |
---|---|
HRP20211153T1 (en) | 2021-10-15 |
PT3123648T (en) | 2021-06-04 |
GB201405259D0 (en) | 2014-05-07 |
EP3123648A2 (en) | 2017-02-01 |
WO2015145138A2 (en) | 2015-10-01 |
PL3123648T3 (en) | 2021-11-22 |
GB2522083B (en) | 2016-02-10 |
US20170155460A1 (en) | 2017-06-01 |
WO2015145138A3 (en) | 2015-11-19 |
AU2015237965A1 (en) | 2016-11-10 |
CN106664148A (en) | 2017-05-10 |
US10263721B2 (en) | 2019-04-16 |
CN106664148B (en) | 2019-12-03 |
BR112016022034A2 (en) | 2017-08-15 |
DK3123648T3 (en) | 2021-06-28 |
US20190199462A1 (en) | 2019-06-27 |
AU2015237965B2 (en) | 2019-06-20 |
GB2522083A (en) | 2015-07-15 |
ES2872332T3 (en) | 2021-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3123648B1 (en) | Simultaneous call transmission detection | |
US8452245B2 (en) | Impedance measurement in an active radio frequency transmitter | |
CN110224721B (en) | Method and receiver for processing an analog signal from a transmission channel | |
US8385449B2 (en) | Method and device for detecting simultaneous double transmission of AM signals | |
US9374291B2 (en) | Downstream OFDM signal egress detection | |
JP6935425B2 (en) | Noise suppression device, noise suppression method, and receiving device and receiving method using these | |
EP3014832B1 (en) | Weak signal detection in double transmission | |
TW201127102A (en) | Applying spurious interference rejection to detect incumbent users of television channels | |
Huang et al. | Resolution doubled co-prime spectral analyzers for removing spurious peaks | |
US7460043B2 (en) | Analog-to-digital converter compensation system and method | |
JP2017129741A (en) | Noise reduction device and noise reduction method | |
US8249176B2 (en) | Method for determination of bridged taps on a transmission line | |
JP2004187135A (en) | Managing device and method | |
JP2006319815A (en) | Pulse noise canceler | |
EP1160969B1 (en) | Detection of distorsions in AM signals | |
Shishkin | Robust digital watermarks for audio signals | |
KR101948314B1 (en) | Apparatus for interference reduction used for polyphase filter bank-based repeater | |
JP4414804B2 (en) | Receiving apparatus and receiving method | |
Peesapaty et al. | Quantifying Uncertainty: Statistical Speech Signal Processing | |
RU2603493C2 (en) | Automated device for increasing of data transfer channel quality | |
KR20110080890A (en) | Dc offset voltage correction apparatus and method in rfid system | |
JP2019009564A (en) | Receiver, reception method and control program | |
JP2018196027A (en) | Reception signal processing apparatus, reception signal processing method, and program | |
JP2019033317A (en) | FSK demodulator | |
JP2012147094A (en) | Symbol estimation circuit and demodulator circuit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TPAC | Observations filed by third parties |
Free format text: ORIGINAL CODE: EPIDOSNTIPA |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20161024 |
|
AK | Designated contracting states |
Kind code of ref document: A2 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) | ||
TPAC | Observations filed by third parties |
Free format text: ORIGINAL CODE: EPIDOSNTIPA |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
17Q | First examination report despatched |
Effective date: 20180906 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R079 Ref document number: 602015068945 Country of ref document: DE Free format text: PREVIOUS MAIN CLASS: H04L0001000000 Ipc: H04J0011000000 |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
RIC1 | Information provided on ipc code assigned before grant |
Ipc: H04L 27/26 20060101ALI20200504BHEP Ipc: H04L 27/06 20060101ALI20200504BHEP Ipc: H04J 11/00 20060101AFI20200504BHEP |
|
INTG | Intention to grant announced |
Effective date: 20200526 |
|
GRAJ | Information related to disapproval of communication of intention to grant by the applicant or resumption of examination proceedings by the epo deleted |
Free format text: ORIGINAL CODE: EPIDOSDIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: EXAMINATION IS IN PROGRESS |
|
GRAP | Despatch of communication of intention to grant a patent |
Free format text: ORIGINAL CODE: EPIDOSNIGR1 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: GRANT OF PATENT IS INTENDED |
|
INTC | Intention to grant announced (deleted) | ||
INTG | Intention to grant announced |
Effective date: 20201019 |
|
GRAS | Grant fee paid |
Free format text: ORIGINAL CODE: EPIDOSNIGR3 |
|
GRAA | (expected) grant |
Free format text: ORIGINAL CODE: 0009210 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE PATENT HAS BEEN GRANTED |
|
AK | Designated contracting states |
Kind code of ref document: B1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
REG | Reference to a national code |
Ref country code: GB Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: EP |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: REF Ref document number: 1391156 Country of ref document: AT Kind code of ref document: T Effective date: 20210515 |
|
REG | Reference to a national code |
Ref country code: IE Ref legal event code: FG4D |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R096 Ref document number: 602015068945 Country of ref document: DE |
|
REG | Reference to a national code |
Ref country code: PT Ref legal event code: SC4A Ref document number: 3123648 Country of ref document: PT Date of ref document: 20210604 Kind code of ref document: T Free format text: AVAILABILITY OF NATIONAL TRANSLATION Effective date: 20210531 |
|
REG | Reference to a national code |
Ref country code: DK Ref legal event code: T3 Effective date: 20210621 |
|
REG | Reference to a national code |
Ref country code: HR Ref legal event code: TUEP Ref document number: P20211153T Country of ref document: HR |
|
REG | Reference to a national code |
Ref country code: NL Ref legal event code: FP |
|
REG | Reference to a national code |
Ref country code: LT Ref legal event code: MG9D |
|
REG | Reference to a national code |
Ref country code: NO Ref legal event code: T2 Effective date: 20210505 |
|
REG | Reference to a national code |
Ref country code: HR Ref legal event code: T1PR Ref document number: P20211153 Country of ref document: HR |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: FI Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: LT Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: BG Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210805 |
|
REG | Reference to a national code |
Ref country code: ES Ref legal event code: FG2A Ref document number: 2872332 Country of ref document: ES Kind code of ref document: T3 Effective date: 20211102 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LV Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: GR Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210806 Ref country code: SE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: RS Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: SK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: SM Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: EE Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: CZ Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: RO Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R097 Ref document number: 602015068945 Country of ref document: DE |
|
PLBE | No opposition filed within time limit |
Free format text: ORIGINAL CODE: 0009261 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT |
|
REG | Reference to a national code |
Ref country code: HR Ref legal event code: ODRP Ref document number: P20211153 Country of ref document: HR Payment date: 20220311 Year of fee payment: 8 |
|
26N | No opposition filed |
Effective date: 20220208 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: AL Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MC Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 |
|
REG | Reference to a national code |
Ref country code: CH Ref legal event code: PL |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: LU Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20220324 Ref country code: LI Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20220331 Ref country code: CH Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES Effective date: 20220331 |
|
REG | Reference to a national code |
Ref country code: HR Ref legal event code: ODRP Ref document number: P20211153 Country of ref document: HR Payment date: 20230310 Year of fee payment: 9 |
|
P01 | Opt-out of the competence of the unified patent court (upc) registered |
Effective date: 20230514 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: ES Payment date: 20230404 Year of fee payment: 9 |
|
REG | Reference to a national code |
Ref country code: AT Ref legal event code: UEP Ref document number: 1391156 Country of ref document: AT Kind code of ref document: T Effective date: 20210505 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: HU Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO Effective date: 20150324 |
|
REG | Reference to a national code |
Ref country code: HR Ref legal event code: ODRP Ref document number: P20211153 Country of ref document: HR Payment date: 20240312 Year of fee payment: 10 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: IS Payment date: 20240320 Year of fee payment: 10 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: NL Payment date: 20240319 Year of fee payment: 10 Ref country code: IE Payment date: 20240319 Year of fee payment: 10 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: AT Payment date: 20240325 Year of fee payment: 10 |
|
PG25 | Lapsed in a contracting state [announced via postgrant information from national office to epo] |
Ref country code: MK Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 Ref country code: CY Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT Effective date: 20210505 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: DE Payment date: 20240322 Year of fee payment: 10 Ref country code: PT Payment date: 20240308 Year of fee payment: 10 Ref country code: GB Payment date: 20240322 Year of fee payment: 10 |
|
PGFP | Annual fee paid to national office [announced via postgrant information from national office to epo] |
Ref country code: PL Payment date: 20240311 Year of fee payment: 10 Ref country code: NO Payment date: 20240223 Year of fee payment: 10 Ref country code: IT Payment date: 20240320 Year of fee payment: 10 Ref country code: HR Payment date: 20240312 Year of fee payment: 10 Ref country code: FR Payment date: 20240320 Year of fee payment: 10 Ref country code: DK Payment date: 20240320 Year of fee payment: 10 Ref country code: BE Payment date: 20240322 Year of fee payment: 10 |