JP2008256501A - Signal analyzer - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a signal analyzer capable of reducing an influence when a deficiency of a periodical signal is generated. <P>SOLUTION: This analyzer is provided with a difference parameter setting part 1 for setting a difference parameter; an interval setting part 3 for setting a range of an arrival time interval or the like; a setting part 5 of the number of signals for setting the number of continuous periodical signals; an interval calculation part 2 for calculating a plurality of arrival time intervals by adding a duration of each signal and the difference parameter together; an interval subdivided section selection part 4 for selecting to which of a plurality of interval subdivided sections, each of the plurality of arrival time intervals corresponds; a plurality of phase accumulation calculation parts 6 for setting arrival time subdivided sections by using the number of the continuous periodical signals, selecting to which arrival time subdivided section, the arrival time and the arrival time interval are pertinent, calculating a complex number based on the arrival time and the arrival time interval in each arrival time subdivided section, and calculating the total value of complex numbers in each arrival time subdivided section; and periodicity determination part 7 for estimating a series number and the arrival time interval of the periodical signal and a time when the periodical signal arrives based on the total value of the complex numbers in each arrival time subdivided section. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

  The present invention estimates the number of series of signals (periodic signals) whose period of arrival (arrival time) is periodic from received signals in which radiated signals radiated from a plurality of radio wave sources are mixed. The present invention relates to a signal analysis apparatus for estimating the period of arrival time when there is an error.

  Conventionally, as a method for detecting such a periodic signal, an invention for obtaining the presence / absence of a periodic signal and the period of the periodic signal based on the frequency of the arrival time interval using the difference between the arrival times of the two signals as the arrival time interval. Yes (see, for example, Patent Document 1).

JP-A-56-80925

  However, since the conventional method for detecting a periodic signal as described above obtains two arrival time intervals based on the arrival time of one radiation signal, if one signal is lost, that is, it originally exists. If a signal is mistakenly recognized as not being present due to a low radiation signal level, the two arrival time intervals cannot be determined correctly. That is, there is a problem that it is impossible to correctly recognize arrival time intervals that are twice as many as the number of missing radiated signals, and the process is erroneous in a situation where many defects occur.

  In addition, the radio wave emitted by the radio wave source is not always received. For example, in the case where the radio wave source scans the antenna beam, the period is only the time when the main beam of the radio wave source is directed toward the antenna of the signal analyzer. The signal cannot be received. Furthermore, even when the radio wave source itself radiates the radio wave intermittently, the periodic signal can be received only during the time when the radio wave source radiates the radio wave. In such a case, the number of periodic signals that can be received is extremely reduced compared to the case where the periodic signals can always be received.

  On the other hand, a signal other than the periodic signal can be received even during a period when the periodic signal cannot be received. Therefore, when processing is performed by the conventional method including data of a time during which a periodic signal cannot be received, the absolute value of the desired complex number output from the accumulator circuit becomes small because the number of periodic signals is small, while the periodic signal The absolute value of complex numbers other than the desired value increases due to the influence of signals other than. As a result, the desired complex absolute value becomes relatively small, making it difficult to detect.

  The present invention has been made to solve the above-described problems, and an object of the present invention is to reduce the influence when a loss of a periodic signal occurs as compared with the conventional method. Even when a periodic signal cannot be received, the effect is reduced and processing is performed, so that the absolute value of the complex number output from the accumulator circuit is prevented from becoming relatively small, and the detection thereof is facilitated. A signal analyzing apparatus capable of performing the above is obtained.

  The signal analyzing apparatus according to the present invention has an arrival time interval in which the arrival time is a period of a periodic signal, and a duration obtained by subtracting the arrival time of the reception signal from the disappearance time, which is the time when the reception signal disappears. A difference parameter setting unit that sets a plurality of different parameters within a predetermined range, an interval setting unit that sets a range of the arrival time interval, and a width of an interval subdivision section obtained by dividing the range into a plurality of ranges A signal number setting unit that sets the number of continuous periodic signals on the assumption that the periodic signal can be continuously received, a plurality of durations of each received signal obtained as observation data, and the difference parameter setting unit An interval calculation unit that calculates a plurality of arrival time intervals by adding a plurality of set different parameters, respectively, and a plurality of arrival time intervals calculated by the interval calculation unit are included in the interval setting unit. An interval subdivision section selection unit that selects which one of the plurality of set interval subdivision sections corresponds to, and an arrival time subdivision section using the number of continuous periodic signals set by the signal number setting unit, Choose which arrival time subdivision interval the arrival time and arrival time interval correspond to, calculate a complex number based on the arrival time and arrival time interval for each arrival time subdivision interval, and calculate a complex number for each arrival time subdivision interval A plurality of phase accumulation calculation units for calculating the total value, and the number of periodic signal sequences and arrival time intervals, and a period for estimating the time at which the periodic signal arrives, based on the total value of complex numbers for each arrival time subdivision section And a sex determination unit.

  The signal analysis apparatus according to the present invention can determine the arrival time interval based on the duration of one signal, and can specify the presence / absence and period of the radiation signal included in the received signal. When one signal is lost, the number of arrival time intervals that cannot be obtained correctly can be set to 1, and the effect when a signal loss occurs can be reduced as compared with the conventional method.

  In addition, since the number of series of periodic signals is estimated by reducing the influence of the time when the radio wave source is not emitting periodic signals, the absolute value of complex numbers becomes relatively small as in the conventional method, making detection difficult. Thus, there is an effect that it is possible to improve the estimation accuracy of the number of sequences.

Embodiment 1 FIG.
A signal analyzing apparatus according to Embodiment 1 of the present invention will be described with reference to FIGS. FIG. 1 is a block diagram showing a configuration of a signal analyzing apparatus according to Embodiment 1 of the present invention. FIG. 2 is a block diagram showing the configuration of the phase accumulation calculation unit of the signal analyzer according to Embodiment 1 of the present invention. In the following, in each figure, the same reference numerals indicate the same or corresponding parts.

  The signal analyzing apparatus according to Embodiment 1 of the present invention extracts, as observation data, the arrival time and duration of each signal included in the received signal from a received signal in which radiation signals from a plurality of radio wave sources are mixed. Based on the extraction result, the number of periodic signals having a periodic arrival time is estimated, and the arrival time interval of the periodic signal having a periodic arrival time and the time at which the periodic signal is received are estimated.

  First, the principle of the present invention will be described.

Consider a case where a radiation signal is received from only one radio wave source. The radiation signal from this single radio wave source is a periodic signal in which the arrival time of each signal is a periodic signal, and the interval between the arrival times of the periodic signals is r. When the total number of received signals is N (N is a natural number of 1 or more) and the arrival time of each periodic signal is _{tTOA, n} (1 ≦ n ≦ N), the arrival time _{tTOA, n} is _expressed by the following formula ( 1).

t _{TOA, n} = t _{TOA, 1} + (n−1) · r (1)

Further, a time obtained by subtracting the arrival time of each signal from the time when each signal disappeared (hereinafter referred to as the disappearance time) is referred to as a duration, and is assumed to be t _{DOC, n} (1 ≦ n ≦ N). First, the arrival time interval Δt _{n, m} shown in the following equation (2) is calculated.

Δt _{n, m} = t _{DOC, n} + tm (2)

  Here, tm is a parameter (hereinafter, referred to as a difference parameter) having a predetermined M value (M is a natural number of 1 or more, m is 1 ≦ m ≦ M) shown in the following equation (3). To do.

tm = tm _min , (tm _min + Δtm), (tm _min + 2 · Δtm),..., (t _max −Δtm), tm _max (3)

Since the difference parameter tm takes M values, the arrival time interval Δt _{n, m} is calculated as N × M for N periodic signals. Next, consider complex numbers _{pn, m} as shown in equation (4).

p _{n, m} = exp (2πjt _{TOA, n} / Δt _{n, m} ) (4)

Note that since the arrival time interval Δt _{n, m} takes N × M values, the complex numbers _{pn, m} are similarly calculated as N × M for N signals. Here, consider the case of a difference parameter tm ₁ (m ₁ is 1 ≦ m ₁ ≦ M) satisfying the following equation (5) among the M difference parameters.

tm ₁ = rt _{DOC, n} (5)

Substituting Equation (5) into Equation (2), the arrival time interval Δt _{n, m1} becomes as shown in Equation (6) below, and Equations (1) and (6) are substituted into Equation (4). Then, the complex numbers _{pn and m1} are as shown in the following equation (7), and have the same value regardless of the number n of the periodic signal. In the equation (7), θ is expressed by the following equation (8).

Δt _{n, m1} = t _{DOC, n} + (rt− _{DOC, n} ) = r (6)
p _{n, m1} = exp (2πj (t _{TOA, 1} + (n−1) · r) / Δt _{n, m1} )
= Exp (2πjt _{TOA, 1} / r)
= Exp (jθ) (7)
θ = 2πt _{TOA, 1} / r (8)

Therefore, the sum of the complex numbers _{pn, m} for each size of the arrival time interval Δt _{n, m} is calculated. Then, since the desired complex numbers _{pn, m} (1 ≦ n ≦ N) in which the arrival time interval Δt _{n, m} is equal to r are all the values of the above equation (7), the complex numbers _{pn, m} The absolute value of the sum is N.

On the other hand, when the arrival time interval Δt _{n, m} is not equal to r, the complex numbers _{pn, m} (1 ≦ n ≦ N) of the signals 1 to N are not all equal to each other. Even if they are calculated, they cancel each other, so they do not pile up, and their absolute values do not increase. Also, signals other than periodic signals do not have the same complex number _{pn, m} as in the above equation (7), and therefore the absolute value of the sum of complex numbers does not increase.

  Therefore, from the magnitude of the absolute value of the complex number calculated for each magnitude of the arrival time interval, it can be determined whether there is a periodic signal for that arrival time interval, and the arrival when the absolute value of the complex number has increased. From the time interval, the cycle of the periodic signal (arrival time interval) can be specified.

However, the above is a case where a periodic signal can always be received while a signal is being received. That is, _assuming that the signal reception time is t _rec , there is a relationship of the following equation (9) between the interval r of the arrival time of the periodic signal and the number N of periodic signals.

t _rec ≈ N · r (9)

  In such a case, the accumulation of complex numbers in the arrival time intervals including r is larger than the increase in complex numbers in other arrival time intervals.

However, for example, if the radio wave source itself does not emit radio waves and can receive a periodic signal for only a fraction of the time compared with the signal reception time t _rec , the number of periodic signals is a fraction of the number. Therefore, the desired complex absolute value also becomes a fraction. On the other hand, since the radiation signal is received even while the periodic signal cannot be received, the non-desired complex absolute value gradually increases, and the difference between the desired complex absolute value and the non-desired complex absolute value is reduced. End up. As a result, there is a case where the number of periodic signal sequences is erroneously estimated and the value of the arrival item interval is erroneous.

  By the way, the periodic signal may be received when the radio source emits radio waves or when the radio wave is stopped (the same applies when the radio source scans the beam). It will be divided into times that cannot be done. Therefore, when calculating the complex number, consider that the total is calculated by dividing based on the arrival time. Then, when the periodic signal can be continuously received, the ratio of the number of periodic signals to the number of radiated signals in the section increases, and the detection becomes easy. Therefore, as a specific method, first, the number of continuous periodic signals s (integer) that can continuously receive periodic signals is set. When the complex numbers are summed for the arrival time interval r, the complex numbers are divided into I shown in the following formulas (10) and (11) according to the arrival time, and the absolute value of the total value is calculated for each.

I = (t _rec − (s + 1/2) · r) / (r / 2) (10)
(I-1) .1 / _{2.r.ltoreq.t TOA, n} <(i-1) .1 / 2.r + (1 + 1/2) r
(1 ≦ i ≦ I) (11)

In the case of this method, the ratio of the number of periodic signals to the number of radiated signals can be increased to calculate the total of complex numbers, so that a periodic signal can be received for only a fraction of the time compared to _{trec as} described above. However, in the section corresponding to the time when the periodic signal can be received, the ratio of the periodic signal increases, and the detection becomes possible. This is the principle of the present invention. Based on this principle, the signal analyzer according to the first embodiment of the present invention will be described in detail.

In FIG. 1, the signal analyzer according to the first embodiment includes a difference parameter setting unit 1, an interval calculation unit 2, an interval setting unit 3, an interval subdivision section selection unit 4, a signal number setting unit 5, K phase accumulation calculation units 6 (6 ₁ , 6 ₂ ,..., 6 _K ) and a periodicity determination unit 7 are provided.

In FIG. 2, the phase accumulation calculation unit 6 _k includes an arrival time subdivision section selection unit 61, I arrival phase calculation units 62 (62 ₁ , 62 ₂ ,..., 62 _I ), I pieces Cumulative calculation unit 63 (63 ₁ , 63 ₂ ,..., 63 _I ) is provided.

  The different parameter setting unit 1 sets M different parameters tm, which is the difference between the arrival time interval and the duration, within a predetermined range. Further, the interval calculation unit 2 obtains the arrival time interval by adding the difference parameter tm to the duration of each signal. Further, the interval setting unit 3 sets the range of the arrival time interval to be analyzed and the width of the interval subdivision section obtained by dividing the range into a plurality. Further, the interval subdivision section selection unit 4 selects an interval subdivision section corresponding to the size of the arrival time interval output from the interval calculation unit 2.

  The signal number setting unit 5 sets the number of continuous periodic signals that can continuously receive periodic signals, and the time width of an arrival time subdivision section described later is determined according to the setting result. The phase accumulation calculation unit 6 selects a corresponding arrival time subdivision section based on the arrival time interval value and the arrival time output from the interval subdivision section selection unit 4, and based on the arrival time interval value and the arrival time. The complex number is obtained, and the total of the complex numbers is obtained for each arrival time subdivision interval. Further, the periodicity determination unit 7 estimates the number of sequences of the periodic signal, its period, and the time when the periodic signal starts to arrive based on the complex number output from the phase accumulation calculation unit 6.

  Next, the operation of the signal analyzing apparatus according to the first embodiment will be described with reference to the drawings. FIG. 3 is a diagram for explaining a method for selecting an interval subdivision section of the signal analyzing apparatus according to the first embodiment of the present invention. Moreover, FIG. 4 is a figure for demonstrating the selection method of the arrival time subdivision section of the signal analyzer which concerns on Embodiment 1 of this invention.

When the total number of signals included in the received signal is N, the arrival time of each signal is observed as _{tTOA, n} (1 ≦ n ≦ N), and the duration is observed as _{tDOC, n} (1 ≦ n ≦ N) Will be described in detail. Here, the relationship between the arrival time, disappearance time, and duration of the received signal is as shown in FIG.

The difference parameter setting unit 1 sets a minimum value tm _min , a maximum value tm _max, and a step interval Δtm of the difference parameter tm, which is a difference between the arrival time interval and the duration of the signal. The minimum value tm _min and the maximum value tm _max are a value obtained by multiplying the duration of each signal by a preset constant (hereinafter referred to as a multiplied value), or a value obtained by dividing the duration by a preset value (hereinafter referred to as division). And so on).

Further, when the minimum value tm _min and the maximum value tm _max are assumed in advance, the values (hereinafter referred to as assumed values) may be used, and among the multiplied value, the divided value, and the assumed value, A maximum value or a minimum value may be used. Further, the step interval Δtm is set assuming a signal to be analyzed. In the following description, it is assumed that the different parameter setting unit 1 sets M (M is a natural number of 1 or more) different parameters as shown in the following equation (12).

tm = tm _min , tm _min + Δtm, tm _min + 2Δtm,..., tm _min + (M−1) Δtm (12)

The interval setting unit 3 sets a range of arrival time intervals to be analyzed below, a width obtained by dividing the range into a plurality of ranges, and the like. Note that the divided intervals of arrival time intervals are referred to as interval subdivision intervals. The range of the arrival time interval to be analyzed is set assuming a periodic signal to be analyzed. In addition, as the width of the interval subdivision section, a value equal to the resolution of the arrival time or duration, or a value obtained by multiplying or dividing the resolution by a certain value is used. In the following, it is assumed that the range of arrival time intervals to be analyzed is set to τ ₀ to τ _K , and that there are K interval subdivision sections. The width of the interval subdivision section represents τ _k −τ _(k−1) (1 ≦ k ≦ K).

  The signal number setting unit 5 sets the number of continuous periodic signals that can be received continuously. The number of continuous signal settings is set assuming a periodic signal to be analyzed. In the following description, it is assumed that the signal number setting unit 5 sets the number of continuous periodic signals to s.

As shown in FIG. 3, the interval calculation unit 2 includes N duration values t _{DOC, n of} each signal obtained as observation data _, and M difference parameters tm obtained from the difference parameter setting unit 1. N × M arrival time intervals Δt _{n and m} are calculated and output to the interval subdivision section selection unit 4. Specifically, the interval calculation unit 2 calculates and outputs the arrival time interval Δt _{n, m} by the following equation (13).

Δt _{n, m} = t _{DOC, n} + tm (13)

The interval subdivision section selection unit 4 includes K interval subdivision sections [τ ₀ , τ ₁ , in which N × M arrival time intervals Δt _{n, m} output from the interval calculation unit 2 are set according to the interval setting unit 3. ), [Τ ₁ , τ ₂ ),..., [Τ _(K−1) , τ _K ). However, the k-th (1 ≦ k ≦ K) section [τ _(k−1) , τ _k ) means a range of interval subdivision sections, and means a section of the following equation (14).

τ _(k−1) ≦ Δt _{n, m} <τ _k (14)

Further, k is the number of the interval subdivision section, and the interval subdivision section [τ _(k−1) , τ _k ) and the interval subdivision section k are synonymous. The interval subdivision section selection unit 4 selects the interval subdivision sections [τ _(k−1) , τ _k ) corresponding to the arrival time intervals Δt _{n, m} according to the above equation (14).

Since the arrival time intervals Δt _{n, m} output from the interval calculation unit 2 are N × M, the interval subdivision section selection unit 4 selects an interval subdivision section corresponding to each of N × M. In FIG. 3B, M different parameters are added to the arrival time of one signal, and an interval subdivision section k ₁ to an interval subdivision section k _(M−1) corresponding to the obtained arrival time interval are selected. It shows how to do.

The phase accumulation calculation unit 6 obtains a complex number using the arrival time interval Δt _{n, m} output from the interval subdivision selection unit 4 and the arrival time t _{TOA, n} (1 ≦ n ≦ N) of the signal, Is calculated for each arrival time subdivision interval corresponding to the complex number obtained based on the arrival time _{tTOA, n} .

Here, it is assumed that the phase accumulation calculation unit 6 _k (1 ≦ k ≦ K) corresponds to the interval subdivision section k, and the operation will be described by taking the _kth phase accumulation calculation unit 6 _k as an example.

First, as shown in FIG. 4, the arrival time subdivision section selection unit 61 constituting the phase accumulation calculation unit 6 _k receives the input arrival time t _{TOA, n} and the arrival time interval Δt _n calculated by the interval calculation unit 2. _{, M} is selected to which arrival time subdivision. For example, FIG. 4A shows a state in which the first received signal corresponds to the arrival time subdivided sections i _{4 to} i ₆ . It is assumed that the arrival time subdivision section is set in advance using the number of continuous periodic signals s set by the signal number setting unit 5. Now, since the interval subdivision section is k, the center value T _k of the arrival time interval is expressed by the following equation (15).

T _k = (τ _(k−1) + τ _k ) / 2 (15)

In this case, the step interval of the arrival time subdivision section is T _k / 2. Further, the range (width) of the arrival time subdivision section is set to (s + 1/2) T _k using the number of continuous periodic signals s as shown in FIG. 4B. Therefore, if the number of arrival time subdivisions is I, the arrival time subdivision i (1 ≦ i ≦ I) is expressed by the following equation (16).

(I−1) T _k / 2 ≦ t _{TOA, n} <(i−1) T _k / 2 + (s + 1/2) T _k
(16)

The arrival time subdivision section selection unit 61 selects which arrival time subdivision section the input arrival time _{tTOA, n} and the arrival time interval Δtn _{, m} correspond to according to the equation (16). It outputs to the arrival phase calculation part 62 _i (1 <= i <= I) provided for every. Note that there may be a plurality of arrival time subdivisions satisfying the above equation (16) for one arrival time. In this case, the arrival time and the arrival time interval are set in the plurality of arrival time subdivisions. Output.

In the following description, the arrival time subdivision section selection unit 61 will be described assuming that the corresponding arrival time subdivision section is the arrival time subdivision section i ₁ that satisfies the following equation (17).

(I ₁ −1) T _k / 2 ≦ t _{TOA, n} <(i ₁ −1) T _k / 2 + (s + _½ ) T _k
(17)

The arrival phase calculation unit 62i ₁ corresponding to the arrival time subdivision section i ₁ uses the following equation (2) based on the arrival time t _{TOA, n} and the arrival time interval Δt _{n, m} output from the arrival time subsection selection unit 61: complex _{p n} according _18), to calculate the _m, and outputs the cumulative calculation unit 63i ₁ corresponding to the arrival time subdivision intervals _{i 1.}

_{pn, m} = exp (2πjt _{TOA, n} / Δt _{n, m} ) (18)

The cumulative calculation unit 63i _1, the complex _{p n} output from the incoming phase calculation section 62i _1, calculates the sum of _m, and outputs the periodicity determining unit 7.

The periodicity determination unit 7 calculates the absolute value of the total value of I × K (I: arrival time subdivision section number, K: interval subdivision section number) complex numbers output from the phase accumulation calculation unit 6 _k. At the same time, it is determined whether or not the absolute value is larger than a preset threshold value. If it is larger, it is determined that the periodic signal is included in the subdivision section.

Now, it is assumed that the absolute value exceeds the threshold value in the interval subdivision section k ₂ and the arrival time subdivision section i ₂ . In this case, it can be seen from the interval subdivision section k ₂ that the arrival time interval is greater than or _equal to τ _(k2-1) and less than τ _k2 , and the time when the periodic signal arrives from the arrival time subdivision section i ₂ is ( It can also be specified that i ₂ −1) T _k2 / 2 or more and less than (i ₂ −1) T _k / 2 + (s + _½ ) T _k2 . The threshold value used in the periodicity determination unit 7 is a value obtained by multiplying or dividing a value (set value) in the continuous periodic signal number s set in the signal number setting unit 5, or a value in s is added or subtracted. Things are used.

In the above description, the description has been made by paying attention to the arrival time of the radiation signal. However, in the arrival time subdivision section selection unit 61, instead of the arrival time, the disappearance time of the radiation signal (arrival time (t _{TOA, n} ) And the duration (t _{DOC, n} ) added value: t _{TOA, n} + t _{DOC, n} ) may be used to select the arrival time subdivision interval.

The arrival time subdivision section selection unit 61 adds the center time (arrival time (t _{TOA, n} ) and disappearance time (t _{TOA, n} + t _{DOC, n} ) of the radiation signal instead of the arrival time to 2 An arrival time subdivision may be selected using a value obtained by division: (2t _{TOA, n} + t _{DOC, n} ) / 2).

  In the arrival phase calculation unit 62, the same effect as described above can be obtained even if the complex number is calculated by the following equation (19) using the disappearance time instead of the arrival time.

p _{n, m} = exp (2πj (t _{TOA, n} + t _{DOC, n} ) / Δt _{n, m} ) (19)

  Further, even when the arrival phase calculation unit 62 calculates the complex number by the following equation (20) using the center time instead of the arrival time, the same effect as described above can be obtained.

_{pn, m} = exp (πj (2t _{TOA, n} + t _{DOC, n} ) / Δt _{n, m} ) (20)

  Furthermore, this method can be applied not only to radio waves but also to signals that receive sound waves, signals that receive light waves, and the like.

It is a block diagram which shows the structure of the signal analyzer based on Embodiment 1 of this invention. It is a block diagram which shows the structure of the phase accumulation calculation part of the signal analyzer which concerns on Embodiment 1 of this invention. It is a figure for demonstrating the selection method of the interval subdivision section of the signal analyzer which concerns on Embodiment 1 of this invention. It is a figure for demonstrating the selection method of the arrival time subdivision section of the signal analyzer which concerns on Embodiment 1 of this invention.

Explanation of symbols

  1 difference parameter setting unit, 2 interval calculation unit, 3 interval setting unit, 4 interval subdivision section selection unit, 5 signal number setting unit, 6 phase accumulation calculation unit, 7 periodicity determination unit, 61 arrival time subdivision section selection unit, 62 Arrival phase calculation unit, 63 cumulative calculation unit.

Claims (8)

The difference parameter that is the difference between the arrival time interval where the arrival time is the period of the periodic signal and the continuation time, which is the time obtained by subtracting the arrival time of the reception signal from the disappearance time, which is the time when the reception signal disappeared, is predetermined. A difference parameter setting unit for setting a plurality of parameters within a range;
An interval setting unit for setting a range of the arrival time interval and an interval subdivision section obtained by dividing the range into a plurality of sections;
A signal number setting unit that sets the number of continuous periodic signals on the assumption that the periodic signal can be continuously received;
An interval calculation unit for calculating a plurality of arrival time intervals by adding a plurality of durations of each received signal obtained as observation data and a plurality of difference parameters set by the difference parameter setting unit; ,
An interval subdivision selection unit that selects which of the plurality of interval subdivisions set according to the interval setting unit corresponds to the plurality of arrival time intervals calculated by the interval calculation unit;
Set the arrival time subdivision using the number of continuous periodic signals set in the signal number setting unit, select which arrival time subdivision corresponds to the arrival time and the arrival time interval, and for each arrival time subdivision Calculating a complex number based on the arrival time and the arrival time interval, and calculating a total number of complex numbers for each arrival time subdivision section;
A signal analysis comprising: a periodicity determination unit configured to estimate the number of sequences of periodic signals, the arrival time interval, and the time at which the periodic signal arrives based on a total value of complex numbers for each arrival time subdivision section apparatus. Each of the plurality of phase accumulation calculation units is
The arrival time subdivision section is set using the arrival time and the number of continuous periodic signals set by the signal number setting unit, and the arrival time subdivision section corresponding to the arrival time and the arrival time interval is selected. Subdivision section selection section,
A plurality of arrival phase calculation units that are provided for each arrival time subdivision section and calculate a complex number based on the arrival time and the arrival time interval output from the arrival time subdivision section selection unit;
The signal analysis device according to claim 1, further comprising: a plurality of cumulative calculation units that calculate a total of complex numbers output from the corresponding arrival phase calculation units. The periodicity determining unit calculates the absolute value of the total of the complex numbers, compares the absolute value with a threshold value calculated based on the number of continuous periodic signals set by the signal number setting unit, and calculates the absolute value. The signal analysis device according to claim 1, wherein when the value is larger than the threshold, it is determined that a periodic signal is included in the arrival time subdivision section. The signal analysis apparatus according to claim 2, wherein the arrival time subdivision section selection unit selects an arrival time subdivision section using an annihilation time obtained by adding the arrival time and the duration time instead of the arrival time. The said arrival time subdivision section selection part selects an arrival time subdivision section using the center time which added the arrival time and the annihilation time, and divided by 2 instead of the said arrival time. Signal analysis equipment. The signal analysis device according to claim 2, wherein the arrival phase calculation unit calculates a complex number based on an disappearance time obtained by adding an arrival time and a duration and an arrival time interval instead of the arrival time. The said arrival phase calculation part calculates a complex number based on the center time and arrival time interval which added the arrival time and the annihilation time, and divided by 2 instead of the said arrival time. Signal analyzer. The signal analysis device according to claim 2, wherein ranges of arrival time subdivision sections corresponding to the arrival phase calculation unit and the accumulation calculation unit are different from each other by a half of the corresponding arrival time intervals.

Images