Method of grouping pulses intercepted by a listening system implemented by computing means of a computer, by analysing a file of recordings comprising measurements of magnitudes associated with each of the intercepted pulses, the magnitudes associated with the intercepted pulse of index k comprising at least one level n k and an arrival time t k , said method grouping the intercepted pulses into C groups, C being an integer equal to or greater than 1, the pulses grouped in a group being considered as having been transmitted by the same source provided with a radar antenna operating in scanning mode, characterized in that the grouping method comprises the steps consisting of: - selecting, from the record file, analysis windows of limited duration to reduce the number of pulses to be analysed simultaneously, a window comprising K pulses, K being an integer equal to or greater than 3; and, for each analysis window, - implementing, from a modelling based on a quadratic regression of the levels of the pulses as a function of the arrival time of the intercepted pulses belonging to the same group, coupled with a mixing model of the groups, said modelling being able to be put in the form: n k = { T k ′ θ 1 + w k si z k , 1 = 1 ⋯ T k ′ θ C + w k si z k , C = 1 , k = 1 ... K where T k is the vector of size t k 2 t k 1 where θ ( c ) = ( α ( c ) ,β ( c ) ,γ ( c ) ), c = 1 ...C, where w k is a sequence of noise samples, and where the indicator variables z k, c are 1 if the k-th pulse belongs to group c and 0 otherwise, an iterative algorithm for estimating the modelling parameters according to a likelihood maximisation principle, the modelling parameters consisting of the number of groups C , the coefficients of the quadratic regression θ ( c ) for each group c, a variance of the measurement noise σ 2 , and the indicator variables of the mixture z k,c ; and, - group the pulses according to the estimated values of the mixture indicator variables. - selecting, from the record file, analysis windows of limited duration to reduce the number of pulses to be analysed simultaneously, a window comprising K pulses, K being an integer equal to or greater than 3; and, for each analysis window, - implementing, from a modelling based on a quadratic regression of the levels of the pulses as a function of the arrival time of the intercepted pulses belonging to the same group, coupled with a mixing model of the groups, said modelling being able to be put in the form: n k = { T k ′ θ 1 + w k si z k , 1 = 1 ⋯ T k ′ θ C + w k si z k , C = 1 , k = 1 ... K where T k is the vector of size t k 2 t k 1 where θ ( c ) = ( α ( c ) ,β ( c ) ,γ ( c ) ), c = 1 ...C, where w k is a sequence of noise samples, and where the indicator variables z k, c are 1 if the k-th pulse belongs to group c and 0 otherwise, an iterative algorithm for estimating the modelling parameters according to a likelihood maximisation principle, the modelling parameters consisting of the number of groups C , the coefficients of the quadratic regression θ ( c ) for each group c, a variance of the measurement noise σ 2 , and the indicator variables of the mixture z k,c ; and, - group the pulses according to the estimated values of the mixture indicator variables. <EOS>