WO2013118557A1 - マルチパス検出方法及び到来角度算出装置 - Google Patents
マルチパス検出方法及び到来角度算出装置 Download PDFInfo
- Publication number
- WO2013118557A1 WO2013118557A1 PCT/JP2013/050896 JP2013050896W WO2013118557A1 WO 2013118557 A1 WO2013118557 A1 WO 2013118557A1 JP 2013050896 W JP2013050896 W JP 2013050896W WO 2013118557 A1 WO2013118557 A1 WO 2013118557A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- unit
- multipath
- output
- arrival angle
- phase
- Prior art date
Links
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q3/00—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
- H01Q3/26—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
- G01S3/48—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems the waves arriving at the antennas being continuous or intermittent and the phase difference of signals derived therefrom being measured
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/318—Received signal strength
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B1/00—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor
- A61B1/04—Instruments for performing medical examinations of the interior of cavities or tubes of the body by visual or photographical inspection, e.g. endoscopes; Illuminating arrangements therefor combined with photographic or television appliances
- A61B1/041—Capsule endoscopes for imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/06—Devices, other than using radiation, for detecting or locating foreign bodies ; determining position of probes within or on the body of the patient
- A61B5/061—Determining position of a probe within the body employing means separate from the probe, e.g. sensing internal probe position employing impedance electrodes on the surface of the body
Definitions
- the present invention relates to an arrival angle calculation device that detects the phase of an incoming radio wave and uses it to calculate the radio wave arrival angle.
- Patent Document 1 proposes an arrival direction estimation device with a reduced calculation scale.
- coefficients of the direction of arrival are calculated by the complex conjugate circuit and the multiplication circuit for received signals received by two antennas, and arc tangent operation and inverse cosine By performing the calculation, the arrival direction of the received wave is estimated.
- Patent Document 2 determines that the point at which the signal power exhibits a local maximum value based on the power of the received signal is the desired wave, and then the signal power is large. A method is disclosed for determining a point as a first delayed wave.
- the delay time of the delayed wave with respect to the desired wave is short, and the delayed wave can not be determined in the received wave in which the desired wave and the delayed wave overlap in time. Therefore, if the direction of arrival is calculated using the received wave in which the desired wave and the delayed wave temporally overlap, the estimation accuracy of the direction of arrival is reduced due to the influence of the delayed wave.
- a method of determining the desired wave and the delayed wave a method of determining the desired wave when the signal power is larger than the threshold and determining the delayed wave when the signal power is smaller than the threshold may be considered. An erroneous determination occurs when there is a delayed wave of a signal level close to the desired wave.
- the present invention has been made in view of the above problems, and an object of the present invention is to provide a multipath detection method capable of detecting multipath accurately and an arrival angle calculation device capable of preventing deterioration in accuracy due to the influence of multipath.
- the multipath detection method is a multipath detection method for determining the presence or absence of multipath from signals received by a plurality of antennas, converting received signals of the respective antennas into frequency domains, and amplitudes of the frequency domains. It is characterized in that the presence or absence of multipath is determined by comparing between antennas for at least one of the characteristic and the frequency characteristic.
- the multipath can be detected with high accuracy by converting the received signal into the frequency domain and comparing the balance between the antennas for at least one of the amplitude characteristic and the phase characteristic.
- the multipath detection method of the present invention it is preferable to compare the balance of both the amplitude characteristic and the phase characteristic of the frequency domain to determine the presence or absence of multipath. According to this configuration, since the presence or absence of multipath can be determined by comparing the amplitude characteristic with the phase characteristic, multipath can be detected more accurately.
- the amplitude characteristic is indicated by a ratio of the signal power on the low band side to the signal power on the high band side at a predetermined width centering on the center frequency of the amplitude curve, and the phase characteristic Is preferably indicated by the inclination of the phase at a predetermined width around the center frequency of the amplitude curve.
- the ratio Pr of Pr1 to Pr2 is used to compare the balance of the amplitude characteristics between the two antennas. It is also possible to calculate and determine that there is a multipath if Pr is separated from 1 by a predetermined value or more.
- the phase tilt of one antenna is ⁇ 1
- the phase tilt of the other antenna is ⁇ 2
- the difference ⁇ between ⁇ 1 and ⁇ 2 is used to compare the balance of phase tilt between two antennas. It is also possible to calculate and determine that there is a multipath if ⁇ is separated from 0 by a predetermined value or more.
- the arrival angle calculation apparatus comprises a plurality of antennas, multipath detection means for determining the presence or absence of multipath from signals received by the plurality of antennas, and arrival of radio waves from the phase difference of the signals received by the plurality of antennas.
- Arrival angle calculation means for calculating an angle, wherein the multipath detection means converts the reception signal of each antenna into the frequency domain, and at least one of the amplitude characteristic and the frequency characteristic of the frequency domain is between antennas It is characterized in that the presence or absence of multipath is determined by comparing the balances of (1) and (2), and when there is no multipath, the arrival angle of radio waves is calculated using the phase difference of the signals received by the plurality of antennas.
- the reception signal of each antenna is converted into the frequency domain, and the presence or absence of multipath is determined by comparing the balance between the antennas for at least one of the amplitude characteristic and the frequency characteristic.
- Can. it is possible to appropriately detect the multipath that affects the calculation accuracy of the arrival angle, and calculate the arrival angle with high accuracy.
- the multipath detection means includes a peak detection unit that detects a peak of a received signal of each antenna, an FFT unit that converts a received signal of each antenna into a frequency domain, and the frequency domain And a multipath determination unit that determines the presence or absence of multipath by comparing the balance between the antennas for at least one of the amplitude characteristic and the frequency characteristic.
- the present invention it is possible to provide a multipath detection method capable of detecting multipath accurately and an arrival angle calculation device capable of preventing degradation of precision due to the influence of multipath.
- FIG. 5A is a diagram showing an example of an output waveform from the timing control unit.
- FIG. 5B is a diagram showing an example of an output waveform from the peak power detection unit.
- FIG. 5C is a diagram showing an example of an output waveform from the FFT unit.
- FIG. 6A is a diagram showing an example of the output waveform of the arctangent part.
- 6B is a diagram illustrating an example of an output waveform of the power calculation unit. It is a schematic diagram which shows the geometrical relationship of the electromagnetic wave which arrives at an antenna. It is a schematic diagram which shows the example of the position detection system containing an arrival angle calculation apparatus. It is a flowchart of arrival angle calculation in an arrival angle calculation device. It is a schematic diagram of the signal input into a peak detection part. When DSSS is used as a modulation system, it is a schematic diagram which shows the example of the signal input into a peak detection part. It is a figure which shows the example of the amplitude frequency curve and phase frequency curve which are output from a FFT part. It is a block diagram which shows another example of an arrival angle calculation part.
- FIG. 10 is a schematic view showing an outline of arrival angle calculation when the phase difference is around + 180 ° or ⁇ 180 °.
- FIG. 7 is a flow diagram of arrival angle calculation when the phase difference is around + 180 ° or ⁇ 180 °.
- FIG. 20A is a schematic view showing a configuration of symbols in OFDM.
- FIG. 20B is a schematic view showing the state of the correlation processing of the OFDM symbol string.
- FIGS. 21A and 21B are diagrams showing examples of output waveforms from the power calculation unit.
- FIG. 21C is a diagram illustrating an example of an output waveform from the adding unit.
- FIG. 21D is a diagram showing an example of an output waveform from each part of the inverse tangent part. It is a schematic diagram which shows the structural example of the capsule endoscope system which used the arrival angle calculation apparatus.
- FIG. 1 is a block diagram showing a configuration example of an arrival angle calculation device according to an embodiment of the present invention.
- the arrival angle calculation device 1 includes a reference signal generation unit 10 capable of oscillating a reference signal at a predetermined oscillation frequency, receiving antennas 11a and 11b arranged with a predetermined interval therebetween, and a receiving antenna 11a. , 11b are converted to reception signals using the reference signal output from the reference signal generation unit 10 and output from the reception units 12a and 12b and the reception signals output from the reception units 12a and 12b And an arithmetic unit 13 for performing various arithmetic processes for calculation.
- the arrival angle calculation device 1 receives radio waves having the same information at two points (or two or more points) separated by a predetermined distance, in order to calculate the arrival angle based on the phase delay caused by propagation delay of radio waves. There is a need. For this reason, it is necessary to provide two (or more) antennas and reception systems corresponding to received radio waves. However, as long as the same incoming radio wave (same information unit) can be received at two or more positions separated by a predetermined interval, the arrival angle calculation device 1 is not limited to the configuration including two or more receiving systems.
- the receiving units 12a and 12b include low noise amplifiers, mixers, band pass filters, and the like, and are configured to be able to receive radio waves of a predetermined frequency.
- the operation unit 13 includes correlation processing units 21a and 21b that perform correlation processing of the received signal, peak detection units 22a and 22b that detect peaks of the reception signal subjected to correlation processing, and peaks detected by the peak detection units 22a and 22b.
- Timing control units 23a and 23b that output signals from the correlation processing units 21a and 21b in accordance with the timing, a determination unit 24 that determines multipath interference based on the outputs of the timing control units 23a and 23b, and timing control And an arrival angle calculation unit 25 that calculates an arrival angle based on the signals from the units 23a and 23b and the determination unit 24.
- the configuration and functions of the arithmetic unit 13 may be realized by hardware or software.
- the correlation processing units 21a and 21b multiply the reception signals from the reception units 12a and 12b by signals having high correlation with the reception signals, and output the result. Since the signals multiplied by the correlation processing units 21a and 21b have high correlation with the received signal, the signals output from the correlation processing units 21a and 21b have peaks in the correlation section.
- the peak detectors 22a and 22b calculate the power of the output signal from the correlation processor 21a and 21b, and detect the power peak of the output signal.
- the timing control units 23a and 23b output the output signals from the correlation processing units 21a and 21b to the determination unit 24 and the arrival angle calculation unit 25 in accordance with the peak timings detected by the peak detection units 22a and 22b.
- the determination unit 24 converts the time domain signals output from the timing control units 23a and 23b into frequency domain signals, and evaluates the presence or absence of multipath in the received wave based on the frequency domain signals. Then, based on the evaluation, it is determined whether or not to use the received signal of a period corresponding to the target information unit (hereinafter referred to as an information unit period) for calculation of the arrival angle. The determination result is notified to the arrival angle calculation unit 25.
- FIG. 2 is a block diagram showing a specific configuration example of the arrival angle calculation apparatus in the case of using direct spread spectrum (DSSS) as a modulation method. Note that FIG. 2 shows only the configuration corresponding to the calculation unit 13 in FIG.
- DSSS direct spread spectrum
- the correlation processing unit 21a generates a spreading code, the spreading code generator 31, multipliers 32a and 32b for multiplying the reception signal by the spreading code, and outputs of the multipliers 32a and 32b for 1 bit ( It includes adders 33a and 33b that add together by the information unit period and output the result to the peak detection unit 22a and the timing control unit 23a.
- the peak detection unit 22a includes a power calculation unit 34a that calculates the power of the signals output from the adders 33a and 33b, and a peak power detection unit 35a that detects the power peak and outputs the detected power peak to the timing control unit 23a.
- the timing control unit 23a includes a buffer unit 36a that controls the timing of outputting the signals from the adders 33a and 33b to the determination unit 24 and the arrival angle calculation unit 25 based on the signal from the peak power detection unit 35a.
- the correlation processing unit 21b includes a spreading code generator 31, multipliers 32c and 32d, and adders 33c and 33d
- the peak detection unit 22b includes a power calculation unit 34b and a peak power detection unit 35b, and performs timing control.
- the unit 23 b includes a buffer unit 36 b.
- FIG. 3 is a block diagram showing a specific configuration example of the determination unit 24.
- the determination unit 24 cuts out the time domain signals output from the timing control units 23a and 23b at arbitrary intervals and converts the signals into frequency domain signals, and the FFT units 37a and 37b.
- the multi-path determination unit 38 determines whether multi-path exists in the target information unit period in the received wave based on the signal of the frequency domain output from 37a and 37b.
- the period in which the FFT units 37a and 37b convert the signal into the frequency domain can be appropriately changed in accordance with the period in which the multipath determination is desired.
- the entire information unit period may be converted to the frequency domain, and when it is desired to determine multipath interference in a period (peak period) near the peak point of the received signal. Only needs to convert the peak period into the frequency domain.
- the determination result (evaluation result) in the multipass determination unit 38 is notified to the averaging unit 45 of the arrival angle calculation unit 25.
- the arrival angle calculation unit 25 performs complex multiplication on the complex conjugate unit 41 that takes complex conjugate of the output of the buffer unit 36a, and complex multiplication of the output of the complex conjugate unit 41 and the output of the buffer unit 36b.
- 42 an inverse tangent unit 43 that performs an inverse tangent operation using the output of the complex multiplication unit 42
- an electric power calculation unit 44 that calculates electric power for each chip section from the output signal of the complex multiplication unit 42
- a determination unit 24 and electric power
- the averaging unit 45 averages the output of the arctangent unit 43 based on the information from the calculating unit 44
- the arrival angle conversion unit 46 converts the output into an arrival angle using the output of the averaging unit 45.
- the averaging unit 45 determines whether to use the output of the inverse tangent unit 43 corresponding to the target information unit period based on the determination result (evaluation result) notified from the multipath determination unit 38, and the determination is performed.
- the output of the arctangent portion 43 can be averaged according to.
- the spreading code generator 31 generates a spreading code for despreading the signal spread on the frequency axis by the DSSS.
- the spreading code corresponds to the spreading code used in code modulation (spreading) on the transmission side.
- the multipliers 32a and 32b multiply the received signal by the spreading code to perform despreading.
- the in-phase component I1 of the reception signal from the receiving unit 12a is input to the multiplier 32a.
- the quadrature component Q1 of the received signal from the receiving unit 12a is input to the multiplier 32b.
- the adders 33a and 33b add the periods (bit intervals) corresponding to the output of the chip intervals of the multipliers 32a and 32b to one bit and output the result.
- FIG. 4A shows an example of the output waveform from the adder 33a.
- FIG. 4B is a partially enlarged view of the output waveform shown in FIG. 4A.
- FIG. 4C shows an example of an output waveform from the adder 33b.
- FIG. 4D is a partially enlarged view of the output waveform shown in FIG. 4C.
- the output signal of the adder 33a and the output signal of the adder 33b are input to the power calculator 34a of the peak detector 22a and the buffer 36a of the timing controller 23a.
- the power calculation unit 34a calculates the power for each chip section from the output signals of the adders 33a and 33b. Specifically, the power calculation unit 34a adds the absolute value of the output signal of the adder 33a corresponding to the in-phase component and the absolute value of the output signal of the adder 33b corresponding to the quadrature component, and adds each chip interval. It is output to the peak power detection unit 35a as power information.
- the peak power detection unit 35a When receiving the power information for each chip section, the peak power detection unit 35a detects a power peak in the received signal and outputs it as power peak information to the buffer unit 36a of the timing control unit 23a.
- the square value of the output signal of the adder 33a and the square value of the output signal of the adder 33b may be added and output to the peak power detection unit 35a.
- the power peak information output from the peak detection unit 22a is information for determining the presence or absence of a peak of the received signal.
- the timing control unit 23a determines that the received signal has a peak at that timing and determines the signal Ia1 and the signal Qa1 for one bit. It is output to the unit 24 and the arrival angle calculation unit 25.
- the FFT unit 37a of the determination unit 24 cuts out the signal Ia1 and the signal Qa1 output from the timing control unit 23a (buffer unit 36a) in an arbitrary period, and converts the signal into a frequency domain signal by fast Fourier transform. As described above, the period during which the FFT unit 37a converts the signal into the frequency domain can be appropriately changed in accordance with the period in which the multipath determination is desired.
- FIG. 5A shows an example of an output waveform (one bit) from the timing control unit 23a (buffer unit 36a).
- FIG. 5B shows an example of an output waveform (one bit) from the peak power detection unit 35a.
- FIG. 5C shows an example of an output waveform from the FFT unit 37a. As shown in FIGS.
- the FFT unit 37a converts a signal in the time domain into a signal in the frequency domain (see FIG. 5C), and the multipath determination unit 38 uses the signal in this frequency domain. To determine the presence or absence of multipath in an arbitrary period.
- the balance between the antennas is compared for at least one of the amplitude characteristic and the phase characteristic of the signal in the frequency domain to determine the presence or absence of multipath.
- the balance of both the amplitude characteristic and the phase characteristic may be compared to determine the presence or absence of multipath.
- Correlation processing unit 21b (spread code generator 31, multipliers 32c and 32d, adders 33c and 33d), peak detection unit 22b (power calculation unit 34b, peak power detection unit 35b), timing control unit 23b (buffer unit 36b)
- the operation and function of the FFT unit 37b are the same as those of the correlation processing unit 21a (the spreading code generator 31, the multipliers 32a and 32b, and the adders 33a and 33b), the peak detection unit 22a (the power calculation unit 34a, and the peak power detection unit 35a).
- the timing control unit 23a (buffer unit 36a) and the FFT unit 37a.
- the FFT unit 37b of the determination unit 24 cuts out the signal Ia2 and the signal Qa2 output from the timing control unit 23b (buffer unit 36b) at an arbitrary time period and converts the signal Ia2 and Qa2 into signals in the frequency domain by fast Fourier transform. It is sent to the determination unit 38.
- the multipath determination unit 38 determines the presence or absence of multipath in the target period, the determination result is notified to the averaging unit 45 of the arrival angle calculation unit 25.
- the received signal input to the correlation processing unit 21 b and the received signal input to the correlation processing unit 21 a are signals obtained by receiving the same radio wave at two points separated by a predetermined interval, and the phases are slightly different. Therefore, the phases of the signal output from the timing control unit 23b and the signal output from the timing control unit 23a slightly differ. Assuming that the output O a1 of the timing control unit 23 a and the output O a2 of the timing control unit 23 b are represented by complex numbers with the signal corresponding to the in-phase component as the real part and the signal corresponding to the quadrature component as the imaginary part. It becomes like (2). Note that ⁇ 1 and ⁇ 2 represent the phase of each signal.
- the output O a1 of the timing control unit 23 a is input to the complex conjugate unit 41 of the arrival angle calculation unit 25.
- Complex conjugate unit 41 outputs the complex conjugate of the output O a1 of the timing controller 23a to the complex multiplier 42. That is, the complex conjugate unit 41 outputs the signal Ia1 and a signal obtained by inverting the sign of the signal Qa1.
- the output O a1 ′ of the complex conjugate unit 41 is expressed by a complex number, it becomes as in the following formula (3).
- Complex multiplier 42 the output O a1' complex conjugate unit 41, and an output O a2 of the timing controller 23b by complex multiplication, signal Ib and the signal Qb arctangent unit 43 and the power calculation unit is multiplication result Output to 44.
- the output O b of the complex multiplication unit 42, the in-phase component Ib and the quadrature component Qb of the output O b are represented by the following equations (4) to (6).
- the inverse tangent unit 43 performs an inverse tangent operation using the output of the complex multiplication unit 42. Specifically, the output signal Ib of the complex multiplication unit 42 is used as a denominator, and the arctangent calculation of the value with the output signal Qb as a numerator is performed.
- FIG. 6A shows an example of an output waveform from the arctangent portion 43. As shown in FIG. The output O arctan of the inverse tangent section 43 corresponds to the phase difference ⁇ 2 - ⁇ 1 and is expressed by the following equation (7).
- the power calculation unit 44 calculates the power for each chip section from the output signal of the complex multiplication unit 42. Specifically, the power calculation unit 44 adds the absolute value of Ib and the absolute value of Qb and outputs the result to the averaging unit 45 as power information for each chip section. Note that the squared value of Ib and the squared value of Qb may be added and output to the averaging unit 45. The example of the output waveform from the electric power calculation part 44 is shown to FIG. 6B.
- the averaging unit 45 averages the output O arctan of the inverse tangent unit 43 based on the determination result notified from the determination unit 24 and outputs the averaged output O arctan to the arrival angle conversion unit 46.
- the averaging unit 45 determines whether the target is the target.
- the output O arctan of the arc tangent unit 43 corresponding to the information unit period is not used for averaging.
- the averaging unit 45 determines the target information unit.
- the output O arctan of the arctangent portion 43 corresponding to the period is used for averaging. In this way, it is possible to calculate the arrival angle by excluding the received signal including the coherent multipath that overlaps with the peak of the desired wave in the time domain, so it is possible to improve the calculation accuracy of the arrival angle.
- the arrival angle conversion unit 46 converts the output into an arrival angle by inverse trigonometric function calculation using the output of the averaging unit 45. For example, an inverse sine operation can be applied as the inverse trigonometric function operation.
- the value obtained by the calculation, that is, the output of the arrival angle conversion unit 46 corresponds to the arrival angle ⁇ (rad).
- the output O arcsin of the incoming angle conversion unit 46 is expressed by the following equation (8). In the following equation, ⁇ (m) is the wavelength of the received wave, and d (m) is the distance between the receiving antennas.
- the arrival angle is obtained by the above-described processing because a geometrical relationship as shown in FIG. 7 is established.
- An angle between radio waves arriving at two receiving antennas 11a and 11b arranged at an interval d (m) with reference to a predetermined direction is defined as ⁇ (rad).
- the propagation distance of the radio wave arriving at the receiving antenna 11b is longer by ⁇ (m) than the propagation distance of the radio wave arriving at the receiving antenna 11a, and the phase delay (phase difference ⁇ 2 - ⁇ 1 (rad)) is It occurs.
- the relationship between the difference ⁇ in propagation distance and the phase difference ⁇ 2 - ⁇ 1 generated in this model can be expressed by the following equation (9) using the wavelength ⁇ (m) of the received wave. In the following formula, ⁇ ⁇ .
- Equation (10) is established from the geometrical relationship of the difference ⁇ of the propagation distance in the above model, the antenna interval d, and the arrival angle ⁇ .
- the arrival angle ⁇ is expressed as the following equation (11).
- Expression (11) corresponds to the processing in the arrival angle conversion unit 46.
- the arrival angle is calculated by the arrival angle calculation device of the present embodiment.
- the position detection system 101 shown in FIG. 8 includes an arrival angle calculation device 1a, another arrival angle calculation device 1b disposed apart from the arrival angle calculation device 1a by a predetermined distance D, and an access point 2 or a user terminal 3 It comprises.
- the access point 2 and the user terminal 3 respectively have a transmission system and a reception system (not shown), and are configured to be capable of bi-directional information transmission (communication). Further, the access point 2 and the user terminal 3 are configured to be able to transmit radio waves for arrival angle calculation to the arrival angle calculation device 1a and the arrival angle calculation device 1b by the transmission systems provided respectively.
- the target of position detection may be either the access point 2 or the user terminal 3.
- the arrival angle calculation device 1a receives radio waves transmitted from the transmission antenna of the access point 2 by the reception antennas 11aa and 11ab, and calculates an arrival angle based on the arrival angle calculation device 1a. Further, the arrival angle calculation device 1b receives radio waves transmitted from the transmission antenna of the access point 2 by the reception antennas 11ba and 11bb, and calculates an arrival angle based on the arrival angle calculation device 1b. If the positional relationship between the arrival angle calculation device 1a and the arrival angle calculation device 1b is known, the position of the access point 2 can be determined from the arrival angles relative to each other. In the case of position detection of the user terminal 3, the arrival angle calculation device 1 a and the arrival angle calculation device 1 b calculate the arrival angle of the radio wave transmitted from the user terminal 3.
- FIG. 9 is a flowchart of arrival angle calculation in the arrival angle calculation device 1 according to the present embodiment.
- the arrival angle calculation device 1 receives a radio wave of an arrival angle calculation target
- the reception units 12a and 12b output reception signals to the correlation processing units 21a and 21b.
- the correlation processing units 21a and 21b perform correlation processing and addition processing of the received signal (step S201).
- the peak detection units 22a and 22b detect the peak value P peak of the power from the output signals of the correlation processing units 21a and 21b. Then, the sum A of the power in the period (peak period t 1 ) near the peak point and the sum B of the power in the period t 2 excluding the peak period t 1 from the 1-bit period (period of information unit) are calculated.
- FIG. 10 schematically shows signals input to the peak detection units 22a and 22b.
- the desired wave to the peak period t 1 is the time period t 3 when adjacent multipath exists, respectively.
- the peak power P peak is the power of the peak point P in FIG. 10
- A is the sum of the power in the peak period t 1
- B is the power in the period t 2 excluding the peak period t 1 from one bit period Is the sum of
- FIG. 11 in the case of using DSSS as a modulation scheme, a peak having a time width about twice the period tc of the spreading code is formed. Therefore, it is possible to make the period of the 2 ⁇ tc and peak time t 1.
- the FFT units 37a and 37b of the determination unit 24 convert received signals represented in the time domain into signals in the frequency domain by fast Fourier transform, and obtain an amplitude frequency curve (amplitude curve) and a phase frequency curve (phase curve). It calculates (step S203). As shown in FIG. 10, in the case where the peak of the desired wave and the multipath overlap in the time domain, when the arrival angle is calculated using this received signal, the calculation accuracy of the arrival angle is reduced due to the influence of the multipath. . Therefore, as shown in the following steps, the arrival angle calculation device 1 according to the present embodiment uses the amplitude characteristic and the phase characteristic calculated from the amplitude frequency curve and the phase frequency curve to obtain multi-path in the target period. Determine the presence or absence. 12A and 12B show an example of an amplitude frequency curve and a phase frequency curve output from the FFT unit 37a, and an amplitude frequency curve and a phase frequency curve output from the FFT unit 37b, respectively.
- the values of Pr 1 and Pr 2 indicate the symmetry about the center frequency of the amplitude frequency curve output from the FFT unit 37a and the symmetry about the center frequency of the amplitude frequency curve output from the FFT unit 37b, respectively. . Since the influence of multipath appears in the direction in which these balances are broken, it can be determined whether the signal in the target period is affected by multipath depending on whether the value of Pr is within the predetermined range. . For example, when the effects of multipath are hardly affected, the values of Pr 1 and Pr 2 approach 1 and the value of Pr also approaches 1. On the other hand, when affected by multipath, the balance between Pr 1 and Pr 2 is lost, and the value of Pr deviates from 1. For example, when the value of Pr is 0.8 to 1.2, it can be determined that the system is not affected by multipath.
- the value of ⁇ represents the difference between the slope of the phase frequency curve output from the FFT unit 37a and the slope of the phase frequency curve output from the FFT unit 37b. Appear in the Therefore, it can be determined whether the signal of the target period is affected by multipath depending on whether the value of ⁇ is within the predetermined range. For example, when the value of ⁇ is from ⁇ 20 ° to + 20 °, it can be determined that it is not affected by multipath.
- step S205 If the value of ⁇ is within the predetermined range (step S205: YES), the multipath determination unit 38 determines that there is no multipath in the target period, and the averaging unit 45 of the arrival angle calculation unit 25 The determination result is notified (step S206).
- step S204 determines whether the value of Pr is out of the predetermined range (step S204: NO) or if the value of ⁇ is out of the predetermined range (step S205: NO).
- the multipath determination unit 38 performs multipath in the target period. Is determined, and the determination result is notified to the averaging unit 45 of the arrival angle calculation unit 25 (step S207).
- the averaging unit 45 that has received the notification of the determination result receives from the inverse tangent unit 43 corresponding to the information unit period including the target period.
- the arrival angle is calculated using the output of (step S208).
- the averaging unit 45 that has received the notification of the determination result starts from the arctangent portion 43 corresponding to the information unit period including the target period. The arrival angle is calculated without using the output of (step S208).
- the arrival angle calculation device 1 converts the received signal into the frequency domain and compares the balance among the antennas for at least one of the amplitude characteristic and the phase characteristic.
- multipath can be detected with high accuracy.
- FIG. 13 is a block diagram for explaining another aspect of the arrival angle calculation unit 25 in the arrival angle calculation device 1.
- Arrival angle calculator shown in FIG. 13. 25, a complex conjugate unit 51 which takes the complex conjugate of the output O a1 of the timing control unit 23a, and an output O a1' complex conjugate unit 51, the output O a2 of the timing controller 23b And an inverse tangent unit 53 that performs an arctangent operation using the output of the complex multiplication unit 52.
- the operations and functions of the complex conjugate unit 51, the complex multiplication unit 52, and the arc tangent unit 53 are the same as the operations and functions of the complex conjugate unit 41, the complex multiplier 42, and the arc tangent unit 43 described above.
- phase difference correction unit 54 that corrects the calculation result based on the calculation result (phase difference) of the inverse tangent unit 53, an averaging unit 55 that averages the output of the phase difference correction unit 54, and a phase difference correction unit 54
- the phase difference recorrection unit 56 corrects the calculation result (average value) of the averaging unit 55 when correction is performed in the above case
- the arrival angle conversion unit 57 converts the arrival angle using the output of the phase difference recorrection unit 56.
- the operation and function of the arrival angle conversion unit 57 are the same as the operation and function of the arrival angle conversion unit 46 described above.
- the phase difference correction unit 54 determines a predetermined value in the calculation result of the arctangent part.
- a process of adding an angle (phase difference) is performed.
- the arrival angle calculation unit 25 of this embodiment projects the phase difference on the coordinates of the phase difference range of ⁇ 180 ° to + 180 ° ( ⁇ to + ⁇ ). Therefore, for example, as shown in FIG. 15A, when the phase difference calculated by the arctangent portion 53 does not reach values around + 180 ° and ⁇ 180 °, the arrival angle is appropriately determined by averaging the values. Can be calculated.
- the phase difference correction unit 54 calculates the inverse tangent unit 53.
- a correction process is performed to add a predetermined angle (phase difference) to the result so that appropriate averaging can be performed. Whether or not the calculation result of the arc tangent unit 53 is a value near + 180 ° or ⁇ 180 ° can be determined based on the distribution of a plurality of phase differences obtained as the calculation result of the arc tangent unit 53.
- the calculation result of the inverse tangent unit 53 is a value near + 180 ° and ⁇ 180 °.
- the angle (phase difference) added by the phase difference correction unit 54 may be, for example, + 90 °, but is not limited to this as long as an appropriate averaging process can be performed. It may be either -90 °, + 180 ° or -180 °.
- the averaging unit 55 averages the output of the phase difference correction unit 54.
- the arrival angle calculation unit 25 performs correction to add a phase difference when a phase difference unsuitable for averaging is calculated, so that the averaging unit 55 can perform appropriate averaging processing.
- the point is that averaging is performed without using a target received signal when there is an interference multipath, as in the averaging unit 45.
- the phase difference re-correction unit 56 corrects the output of the averaging unit 55 when the phase difference correction unit 54 corrects the phase difference. Specifically, correction is performed to reduce the angle (phase difference) added as a correction value in the phase difference correction unit 54.
- FIG. 16 schematically shows the arrival angle calculation when the phase difference is around + 180 ° and ⁇ 180 °.
- the phase difference correction unit 54 adds a correction value (+ 90 °) to the phase difference to rotate the coordinate axis, Convert to the coordinate axis for calculating the average value.
- the averaging unit 55 calculates an average value ( ⁇ 92 °) based on the data.
- the phase difference re-correction unit 56 performs correction to reduce the correction value (+ 90 °) from the output data of the phase difference correction unit 54, and outputs the corrected data (+ 178 °) to the arrival angle conversion unit 57.
- FIG. 17 is a processing flow diagram of the arrival angle calculation unit 25.
- Complex conjugate unit 51 of the arrival angle calculation unit 25, in step 301, calculates the complex conjugate of the output O a1 of the timing controller 23a.
- complex multiplier 52 in step 302, multiplying the output O a1' output O a2 and complex conjugate unit 51 of the timing controller 23b.
- the arctangent unit 53 performs arctangent operation using the output of the complex multiplication unit 52 to calculate the phase difference between the received signals.
- step 304 the phase difference correction unit 54 determines whether the calculated phase difference is a value near + 180 ° and ⁇ 180 ° in the IQ plane. If the calculated phase difference is not a value near + 180 ° and ⁇ 180 °, the process proceeds to step 305, and the arrival angle calculation unit 25 calculates an arrival angle without correcting the phase difference. If the calculated phase difference is around + 180 ° or around ⁇ 180 °, the process proceeds to step 306. The determination is performed based on whether the number of retardations larger than + 90 ° or smaller than ⁇ 90 ° is larger than the number of retardations smaller than + 90 ° and larger than ⁇ 90 ° as described above Can.
- step 306 the phase difference correction unit 54 performs a process of adding 90 ° to the phase difference which is the calculation result of the arc tangent unit 53 (increases the phase difference by + 90 °).
- step 307 the averaging unit 55 averages the output of the phase difference correction unit 54.
- step 308 the phase difference recorrection unit 56 performs a process of subtracting 90 ° from the average value which is the calculation result of the averaging unit 55 (makes the phase difference ⁇ 90 °).
- the arrival angle conversion unit 57 calculates an arrival angle from the output of the phase difference recorrection unit 56.
- the arrival angle calculation unit 25 shown in FIG. 13 calculates an appropriate average value by performing a series of processes of adding and averaging a predetermined phase difference and reducing the predetermined phase difference. There is no need to reduce the calculation accuracy of. As a result, the calculation accuracy of the arrival angle can be sufficiently improved.
- phase difference correction unit 54 performs the process of adding a predetermined angle to the calculation result of the arc tangent unit 53
- the present invention is not limited to this as long as an appropriate averaging process can be realized.
- the arrival angle calculation unit 25 configured as shown in FIG. 18 can also be used.
- a complex multiplication unit 62 for performing complex multiplication.
- the operations and functions of the complex conjugate unit 61 and the complex multiplication unit 62 are the same as the operations and functions of the complex conjugate unit 41 and the complex multiplication unit 42 described above. Further, an IQ comparison unit 63 that compares the absolute value of the in-phase component (I component) of the output of the complex multiplication unit 62 with the absolute value of the quadrature component (Q component), and the IQ comparison unit using the output of the complex multiplication unit 62 An inverse tangent unit 64 which performs an inverse tangent operation by selecting and changing the operation method according to the output of 63.
- a correction unit 66 and an arrival angle conversion unit 67 that converts the output of the phase difference recorrection unit 66 into an arrival angle are provided.
- the operation and function of the arrival angle conversion unit 67 are the same as the operation and function of the arrival angle conversion unit 46 described above.
- the IQ comparison unit 63 determines whether or not the in-phase component (I component) of the output of the complex multiplication unit is negative, and the absolute value and quadrature component of the in-phase component (I component) of the output of the complex multiplication unit 62 The absolute value of the Q component is compared. Specifically, the IQ comparison unit 63 determines the sign of the in-phase component Ib, and whether the absolute value
- the in-phase component Ib becomes negative (Ib ⁇ 0), and the absolute value
- the inverse tangent unit 64 uses the output of the complex multiplication unit 62, selects an operation method according to the output of the IQ comparison unit 63, and performs an inverse tangent operation.
- of the in-phase component is approximately the same as or smaller than the absolute value
- An inverse tangent operation is performed on a value with the output Ib as a denominator and the output Qb as a numerator.
- the denominator of -Qb obtained by inverting the sign of the output Qb of the complex multiplication unit 62 Then, the arctangent operation of the value with the output Ib as the numerator is performed.
- of the quadrature component corresponds to the processing of rotating the coordinate axis by + 90 ° and performing the arctangent operation. That is, the phase difference obtained by the process is a value obtained by adding + 90 ° to the original phase difference.
- of the quadrature component is not limited to the above.
- an inverse tangent operation may be performed on a value with the output Qb of the complex multiplication unit 62 as a denominator and the sign of the output Ib inverted and -Ib as a numerator.
- the processing corresponds to processing of rotating the coordinate axis by -90 ° and performing arctangent calculation. That is, the phase difference obtained by the process is a value obtained by adding -90 ° to the original phase difference (a value obtained by subtracting + 90 °).
- the arc tangent operation may be performed by inverting the sign of the output Ib of the complex multiplication unit 62 and the sign of the output Qb.
- the process corresponds to the process of rotating the coordinate axis by + 180 ° (or -180 °) to perform the arctangent operation. That is, the phase difference obtained by the process is a value obtained by adding + 180 ° (or -180 °) to the original phase difference. An appropriate average value can be calculated also by such processing.
- the averaging unit 65 averages the output of the arctangent unit 64.
- the arrival angle calculation unit 25 performs correction to substantially add (or reduce) the phase difference when the phase difference unsuitable for averaging is calculated. Processing is possible.
- the phase difference recorrection unit 66 corrects the output of the averaging unit 65 when the arctangent unit 64 performs a process of rotating the coordinate axis by + 90 °. Specifically, correction is performed to reduce + 90 °.
- correction to reduce -90 ° that is, correction to add + 90 °
- correction is performed to reduce the + 180 ° (or -180 °).
- the arrival angle calculation unit 25 shown in FIG. 18 can also calculate an appropriate average value in the same manner as the arrival angle calculation unit 25 shown in FIG. As a result, the calculation accuracy of the arrival angle can be sufficiently improved.
- FIG. 19 is a block diagram showing a specific configuration example of the arrival angle calculation device in the case of using orthogonal frequency division multiplexing (OFDM) as a modulation scheme. Note that FIG. 19 shows only the configuration corresponding to the calculation unit 13 in FIG.
- OFDM orthogonal frequency division multiplexing
- the correlation processing unit 21a includes a complex conjugate unit 71a that takes a complex conjugate of the output of the receiving unit 12a, a delay unit 72a that delays the output of the receiving unit 12a by a predetermined period and outputs the delayed signal, and a complex conjugate unit 71a.
- a complex multiplication unit 73a that performs complex multiplication on the output of the delay unit 72a, and adders 74a and 74b that add the output of the complex multiplication unit 73a for a GI (guard interval) period and output the result.
- the peak detection unit 22a includes a power calculation unit 75a that calculates the power of the signals output from the adders 74a and 74b, and a peak power detection unit 76a that detects the power peak and outputs the detected power peak to the timing control unit 23a.
- the timing control unit 23a includes a delay unit 77a that controls the output timing of the signal from the reception unit 12a to the determination unit 24 and the arrival angle calculation unit 25 based on the signal from the peak power detection unit 76a.
- the correlation processing unit 21b includes a complex conjugate unit 71b, a delay unit 72b, a complex multiplication unit 73b, and adders 74c and 74d
- the peak detection unit 22b includes a power calculation unit 75b and a peak power detection unit 76b
- the timing control unit 23b includes a delay unit 77b.
- the determination unit 24 is the same as that shown in FIG. That is, the determination unit 24 outputs signals from the FFT units 37a and 37b, which cut out the time domain signals output from the timing control units 23a and 23b at arbitrary intervals and convert the signals into frequency domain signals, and the FFT units 37a and 37b. And a multipath determination unit 38 that determines whether or not there is a multipath in a target information unit period in the received wave based on the signal of the frequency domain to be detected. The period in which the FFT units 37a and 37b convert the signal into the frequency domain can be appropriately changed in accordance with the period in which the multipath determination is desired.
- the entire information unit period may be converted to the frequency domain, and when it is desired to determine multipath interference in a period (peak period) near the peak point of the received signal. Only needs to convert the peak period into the frequency domain.
- the determination result (evaluation result) in the multi-pass determining unit 38 is configured to be notified to the averaging unit 85 of the arrival angle calculating unit 25.
- the delay units 77a to Ia1 and Qa1 are similar to the case of FIG. The input is performed, and Ia2 and Qa2 are input from the delay unit 77b.
- the arrival angle calculation unit 25 includes a complex conjugate unit 81 that takes a complex conjugate of the output of the delay unit 77a, a complex multiplication unit 82 that complex-multiplies the output of the complex conjugate unit 81 and the output of the delay unit 77b, and a complex multiplication unit 42.
- averaging unit 85 averages the output of the arctangent unit 84 based on the above, and the arrival angle conversion unit 86 converts the output into an arrival angle using the output of the averaging unit 85.
- the averaging unit 85 determines whether to use the output of the inverse tangent unit 83 corresponding to the target information unit period based on the determination result (evaluation result) notified from the multipath determination unit 38, and the determination is performed.
- the output of the arc tangent unit 83 can be averaged according to the
- the delay units 72a and 72b delay the output of the reception unit 12a by a predetermined period and output the result, in order to obtain the autocorrelation of the OFDM symbol sequence.
- the delay units 72a and 72b input the complex multiplication unit 73a at the same timing that the end of the OFDM symbol output by the complex conjugate unit 71a and the GI (guard interval) output by the delay units 72a and 72b are the same.
- the output of the receiving unit 12a is delayed for a predetermined period and output.
- the complex multiplication unit 73a performs complex multiplication on the output of the complex conjugate unit 71a and the output of the delay unit 72a.
- the adders 74a and 74b add together the outputs of the complex multiplier 73a for each chip section for the GI period and output the result.
- FIG. 20A is a schematic view showing a configuration of an OFDM symbol string.
- the OFDM symbol sequence is configured of an OFDM symbol which is a data part and a GI placed at the beginning of the OFDM symbol.
- GI is data obtained by copying the end of an OFDM symbol, and is inserted to prevent interference between OFDM symbols.
- FIG. 20B is a schematic view showing the state of correlation processing (autocorrelation processing) of the OFDM symbol string in the correlation processing unit 21a.
- the output of the delay unit 72a lags the output of the complex conjugate unit 71a by the OFDM symbol length.
- autocorrelation can be obtained by multiplying the output of the complex conjugate unit 71a and the output of the delay unit 72a.
- the autocorrelation value shows a peak when the same data as GI appears in the output of the complex conjugate unit 71a and the output of the delay unit 72a. The beginning can be detected.
- the output signals of the adders 74a and 74b are input to the power calculator 75a of the peak detector 22a.
- the power calculator 75a calculates the power for each chip section from the output signals of the adders 74a and 74b.
- the power calculation unit 34 a adds the absolute value of the output signal corresponding to the in-phase component and the absolute value of the output signal corresponding to the quadrature component, and sets the peak power detection unit 76 a as power information for each chip section.
- Output to The square value of the output signal corresponding to the in-phase component and the square value of the output signal corresponding to the quadrature component may be added and output to the peak power detection unit 76a.
- FIG. 21A shows an example of an output waveform from the power calculation unit 75a.
- the peak power detection unit 76a detects a power peak in the received signal, and outputs the power peak information to the delay unit 77a of the timing control unit 23a.
- the power peak information output from the peak detection unit 22a is information for determining the presence or absence of a peak of the received signal.
- the peak period t 1 is equal to the GI period.
- one symbol period corresponds to a period obtained by adding the GI period and the data period (OFDM symbol period).
- the timing control unit 23a determines that the received signal from the receiving unit 12a is a signal that the received signal has a peak at that timing. And output to the arrival angle calculation unit 25.
- the FFT unit 37a of the determination unit 24 cuts out the signal Ia1 and the signal Qa1 output from the timing control unit 23a (the delay unit 77a) at an arbitrary time period, converts the signal into a frequency domain signal by fast Fourier transform, and performs multipath determination unit Send to 38
- the period during which the FFT unit 37a converts the signal into the frequency domain can be appropriately changed in accordance with the period in which the multipath determination is desired.
- the multipath determination unit 38 determines the presence or absence of multipath in the information unit period using the signal in the frequency domain. Specifically, the balance between the antennas is compared for at least one of the amplitude characteristic and the phase characteristic of the signal in the frequency domain to determine the presence or absence of multipath.
- the balance of both the amplitude characteristic and the phase characteristic may be compared to determine the presence or absence of multipath.
- Correlation processing unit 21b (complex conjugate unit 71b, delay unit 72b, complex multiplication unit 73b, adders 74c and 74d), peak detection unit 22b (power calculation unit 75b, peak power detection unit 76b), timing control unit 23b (delay unit) 77b), the operation and function of the FFT unit 37b are the correlation processing unit 21a (complex conjugate unit 71a, delay unit 72a, complex multiplication unit 73a, adders 74a and 74b), peak detection unit 22a (power calculation unit 75a, peak power
- the operations and functions of the detection unit 76a), the timing control unit 23a (delay unit 77a), and the FFT unit 37a are the same.
- the FFT unit 37b of the determination unit 24 cuts out the signal Ia2 and the signal Qa2 output from the timing control unit 23b (delay unit 77b) at an arbitrary time period and converts the signal Ia2 and Qa2 into signals in the frequency domain by fast Fourier transform. It is sent to the determination unit 38.
- the multipath determination unit 38 determines the presence or absence of multipath in the target period, the determination result is notified to the averaging unit 85 of the arrival angle calculation unit 25.
- the received signal input to the correlation processing unit 21b and the received signal input to the correlation processing unit 21a have slightly different phases because the same radio wave is received at two points separated by a predetermined interval. Therefore, the phases of the signal output from the timing control unit 23b and the signal output from the timing control unit 23a slightly differ.
- the output of the timing control unit 23a is input to the complex conjugate unit 81 of the arrival angle calculation unit 25.
- the complex conjugate unit 81 outputs the complex conjugate of the output of the timing control unit 23a to the complex multiplication unit 82.
- the complex multiplication unit 82 performs complex multiplication on the output of the complex conjugate unit 81 and the output of the timing control unit 23b, and outputs the operation result to the addition units 83a and 83b.
- the adders 83a and 83b add the outputs of the complex multipliers 82 for each chip interval by the GI period, and output the result to the arc tangent unit 84.
- FIG. 21C shows an example of output waveforms from the adding units 83a and 83b. In the drawing, the output waveform of the adding unit 83a is indicated by I, and the output waveform of the adding unit 83b is indicated by Q.
- the inverse tangent unit 84 performs an inverse tangent operation using the outputs of the adding units 83a and 83b to calculate the phase difference of the reception signal.
- FIG. 21D shows an example of an output waveform from the arc tangent unit 84.
- the averaging unit 85 averages the output of the arctangent unit 84 based on the determination result notified from the determination unit 24 and outputs the averaged output to the arrival angle conversion unit 86.
- the determination unit 24 notifies the averaging unit 85 of the determination result that the received signal in the target period includes a multipath
- the averaging unit 85 determines the target period.
- the output of the arc tangent unit 84 corresponding to the included information unit period is not used for averaging.
- the averaging unit 85 includes the information including the target period.
- the output of the arc tangent unit 84 corresponding to the unit period is used for averaging.
- the arrival angle can be calculated excluding the received signal including the coherent multipath, so that the calculation accuracy of the arrival angle can be improved.
- the arrival angle conversion unit 86 uses the output of the averaging unit 85 to convert into an arrival angle by inverse trigonometric function calculation. The value obtained by the calculation, that is, the output of the arrival angle conversion unit 86 corresponds to the arrival angle.
- the reception signal is converted into the frequency domain, and the balance between the antennas is compared for at least one of the amplitude characteristic and the phase characteristic.
- multipath can be detected with high accuracy.
- FIG. 22 is a schematic view showing a capsule endoscope system in which the arrival angle calculation device 1 is applied to position specification of a capsule endoscope.
- the capsule endoscope system shown in FIG. 22 includes a plurality of sensor arrays 401, and a data recorder 402 that records data from the sensor arrays 401.
- the sensor array 401 includes an antenna corresponding to the reception antenna of the arrival angle calculation device 1 and is configured to be able to receive radio waves from a capsule endoscope swallowed by a patient.
- the data recorder 402 specifies the position of the capsule endoscope swallowed by the patient from the phase information of the radio wave received by the sensor array 401.
- the capsule endoscope swallowed by the patient is moved by the peristaltic movement of the digestive tract.
- the position of the capsule endoscope is monitored, and it can be confirmed whether or not the examination site has been reached.
- the capsule endoscope images the situation of the examination site and transmits it to the data recorder 402, and the data recorder 402 records image information.
- the power of the camera etc. is turned on at the timing when the capsule endoscope reaches the examination site, and the power of the camera etc. can be turned off when the examination site is removed, so the battery capacity can be reduced. .
- the number of sensors (antennas) can be reduced.
- a large number of images can be transmitted as compared with the conventional capsule endoscope, and a clear image can be obtained.
- an excellent capsule endoscope system can be constructed.
- the received signal is converted into the frequency domain, and at least one of the amplitude characteristic and the phase characteristic is compared with each other to compare the balance among the antennas.
- the path can be detected accurately. As a result, it is possible to accurately calculate the arrival angle by appropriately removing the multipath.
- this invention is not limited to the description of the said embodiment, It can change suitably in the aspect which exhibits the effect. Further, in the above-described embodiment, the configuration and the like shown in the attached drawings are not limited to this, and can be appropriately modified within the range in which the effects of the present invention are exhibited.
- the arrival angle calculation device of the present invention can be used in a system for specifying the position of an object and various other applications.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Quality & Reliability (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Endoscopes (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Noise Elimination (AREA)
Abstract
Description
Claims (7)
- 複数のアンテナで受信した信号からマルチパスの有無を判定するマルチパス検出方法であって、
各アンテナの受信信号をそれぞれ周波数領域に変換し、前記周波数領域の振幅特性と位相特性との少なくとも一つについてアンテナ間で比較を行うことでマルチパスの有無を判定することを特徴とするマルチパス検出方法。 - 前記周波数領域の振幅特性と位相特性とをそれぞれ比較を行うことでマルチパスの有無を判定することを特徴とする請求項1に記載のマルチパス検出方法。
- 前記振幅特性は、振幅曲線の中心周波数を中心とした所定幅での低域側の信号電力と高域側の信号電力との比で示され、
前記位相特性は、振幅曲線の中心周波数を中心とした所定幅での位相の傾斜で示されることを特徴とする請求項1又は請求項2に記載のマルチパス検出方法。 - 一方のアンテナの振幅特性をPr1、他方のアンテナの振幅特性をPr2とし、2つのアンテナ間の振幅特性のバランスを比較するためにPr1とPr2との比Prを計算し、
Prが1から所定値以上離れていればマルチパスがあると判定することを特徴とする請求項3に記載のマルチパス検出方法。 - 一方のアンテナの位相傾斜をΔφ1、他方のアンテナの位相傾斜をΔφ2とし、2つのアンテナ間の位相傾斜のバランスを比較するためにΔφ1とΔφ2との差Δφを計算し、
Δφが0から所定値以上離れていればマルチパスがあると判定することを特徴とする請求項3又は請求項4に記載のマルチパス検出方法。 - 複数のアンテナと、複数のアンテナで受信した信号からマルチパスの有無を判定するマルチパス検出手段と、前記複数のアンテナで受信した信号の位相差から電波の到来角度を算出する到来角度算出手段と、を備え、
前記マルチパス検出手段は、各アンテナの受信信号をそれぞれ周波数領域に変換し、前記周波数領域の振幅特性と位相特性との少なくとも一つについてアンテナ間でのバランスを比較してマルチパスの有無を判定し、
マルチパスがない場合に、前記複数のアンテナで受信した信号の位相差を用いて電波の到来角度を算出することを特徴とする到来角度算出装置。 - 前記マルチパス検出手段は、
各アンテナの受信信号のピークを検出するピーク検出部と、
各アンテナの受信信号をそれぞれ周波数領域に変換するFFT部と、
前記周波数領域の振幅特性と位相特性との少なくとも一つについてアンテナ間でのバランスを比較してマルチパスの有無を判定するマルチパス判定部と、
を備えたことを特徴とする請求項6に記載の到来角度算出装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201380008644.9A CN104115021B (zh) | 2012-02-08 | 2013-01-18 | 多路径检测方法以及到达角度计算装置 |
JP2013557449A JP5860068B2 (ja) | 2012-02-08 | 2013-01-18 | マルチパス検出方法及び到来角度算出装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-025283 | 2012-02-08 | ||
JP2012025283 | 2012-02-08 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013118557A1 true WO2013118557A1 (ja) | 2013-08-15 |
Family
ID=48947326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2013/050896 WO2013118557A1 (ja) | 2012-02-08 | 2013-01-18 | マルチパス検出方法及び到来角度算出装置 |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP5860068B2 (ja) |
CN (1) | CN104115021B (ja) |
WO (1) | WO2013118557A1 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013217669A (ja) * | 2012-04-04 | 2013-10-24 | Toyota Motor Corp | 位相モノパルスレーダ装置 |
WO2019038088A1 (de) * | 2017-08-22 | 2019-02-28 | Continental Automotive Gmbh | Vorrichtung und verfahren zum bestimmen von qualitätseigenschaften für winkelmessung für ein kraftfahrzeug |
US10557914B2 (en) | 2015-03-31 | 2020-02-11 | Mitsubishi Heavy Industries Machinery Systems, Ltd. | Radio wave arrival angle detection device, vehicle detection system, radio wave arrival angle detection method, and vehicle erroneous detection prevention method |
WO2022034842A1 (ja) * | 2020-08-13 | 2022-02-17 | 株式会社ソシオネクスト | 推定方法および推定装置 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI609577B (zh) | 2016-01-29 | 2017-12-21 | 財團法人工業技術研究院 | 執行上行鏈路通道估測的方法及使用所述方法的基地台 |
CN106896378B (zh) * | 2017-03-14 | 2018-12-28 | 千寻位置网络有限公司 | 一种车载定位系统中多路径区域判别的方法 |
US11585916B2 (en) * | 2018-06-29 | 2023-02-21 | Apple Inc. | Ranging using beamforming antenna weight vectors and secure ranging packet with phase tracking |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10322303A (ja) * | 1997-05-15 | 1998-12-04 | Fujitsu Ten Ltd | 直交周波数分割多重方式用受信機 |
JP2006010345A (ja) * | 2004-06-22 | 2006-01-12 | Omron Corp | タグ通信装置、タグ通信装置の制御方法、タグ通信制御プログラム、および、タグ通信管理システム。 |
JP2007010639A (ja) * | 2004-08-16 | 2007-01-18 | Rcs:Kk | アクティブタグ装置 |
JP2008045954A (ja) * | 2006-08-11 | 2008-02-28 | Brother Ind Ltd | 方向検出装置 |
JP2010060459A (ja) * | 2008-09-04 | 2010-03-18 | Hitachi Ltd | モノパルスレーダ装置 |
JP2010276531A (ja) * | 2009-05-29 | 2010-12-09 | Hitachi Ltd | 到来方向推定装置及び到来方向推定方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL8502695A (nl) * | 1985-10-02 | 1987-05-04 | Philips Nv | Ontvanger voorzien van een multipaddetektor. |
US5694131A (en) * | 1996-03-01 | 1997-12-02 | Hughes Electronics | Method and apparatus for detecting multipath interference in a radar receiver |
US20030128660A1 (en) * | 2002-01-09 | 2003-07-10 | Atsushi Ito | OFDM communications apparatus, OFDM communications method, and OFDM communications program |
GB0206766D0 (en) * | 2002-03-22 | 2002-05-01 | Koninkl Philips Electronics Nv | Method of, and apparatus for, determining position |
-
2013
- 2013-01-18 JP JP2013557449A patent/JP5860068B2/ja active Active
- 2013-01-18 WO PCT/JP2013/050896 patent/WO2013118557A1/ja active Application Filing
- 2013-01-18 CN CN201380008644.9A patent/CN104115021B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10322303A (ja) * | 1997-05-15 | 1998-12-04 | Fujitsu Ten Ltd | 直交周波数分割多重方式用受信機 |
JP2006010345A (ja) * | 2004-06-22 | 2006-01-12 | Omron Corp | タグ通信装置、タグ通信装置の制御方法、タグ通信制御プログラム、および、タグ通信管理システム。 |
JP2007010639A (ja) * | 2004-08-16 | 2007-01-18 | Rcs:Kk | アクティブタグ装置 |
JP2008045954A (ja) * | 2006-08-11 | 2008-02-28 | Brother Ind Ltd | 方向検出装置 |
JP2010060459A (ja) * | 2008-09-04 | 2010-03-18 | Hitachi Ltd | モノパルスレーダ装置 |
JP2010276531A (ja) * | 2009-05-29 | 2010-12-09 | Hitachi Ltd | 到来方向推定装置及び到来方向推定方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013217669A (ja) * | 2012-04-04 | 2013-10-24 | Toyota Motor Corp | 位相モノパルスレーダ装置 |
US10557914B2 (en) | 2015-03-31 | 2020-02-11 | Mitsubishi Heavy Industries Machinery Systems, Ltd. | Radio wave arrival angle detection device, vehicle detection system, radio wave arrival angle detection method, and vehicle erroneous detection prevention method |
WO2019038088A1 (de) * | 2017-08-22 | 2019-02-28 | Continental Automotive Gmbh | Vorrichtung und verfahren zum bestimmen von qualitätseigenschaften für winkelmessung für ein kraftfahrzeug |
EP4246172A3 (de) * | 2017-08-22 | 2023-12-06 | Continental Automotive Technologies GmbH | Vorrichtung und verfahren zum bestimmen von qualitätseigenschaften für winkelmessung für ein kraftfahrzeug |
WO2022034842A1 (ja) * | 2020-08-13 | 2022-02-17 | 株式会社ソシオネクスト | 推定方法および推定装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104115021B (zh) | 2016-03-23 |
JP5860068B2 (ja) | 2016-02-16 |
CN104115021A (zh) | 2014-10-22 |
JPWO2013118557A1 (ja) | 2015-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2013118557A1 (ja) | マルチパス検出方法及び到来角度算出装置 | |
JP5677452B2 (ja) | 到来角度算出装置 | |
US8797213B2 (en) | Methods and systems for estimating angle of arrival | |
JP6035165B2 (ja) | レーダ装置 | |
JP5535024B2 (ja) | レーダ装置 | |
US20130088395A1 (en) | Method and apparatus for estimating direction of arrival | |
JP2019522220A (ja) | 速度の曖昧さの解消を含むmimoレーダにおける速度検出のための方法及び装置 | |
JP5250336B2 (ja) | タイミング同期方法およびその装置、並びにプリアンブルおよびその生成方法と装置 | |
JP3600459B2 (ja) | 電波到来方向推定方法及びその装置 | |
US20060007904A1 (en) | Apparatus and method for fading frequency estimation | |
US9124480B2 (en) | Method and apparatus for spectrum sensing | |
JP5185233B2 (ja) | 無線通信装置 | |
JP5729846B2 (ja) | マルチパス検出方法、マルチパス検出装置及び到来角度算出装置 | |
JP4486950B2 (ja) | Oqpsk復調器のタイミング推定器 | |
US20110032974A1 (en) | Carrier frequency synchronization detection circuit and correlation calculator | |
WO2004088867A1 (ja) | 受信装置 | |
KR100273714B1 (ko) | 희망신호 방향추정기 | |
JP5188204B2 (ja) | 測距通信装置 | |
JP3073919B2 (ja) | 同期装置 | |
KR20190087228A (ko) | 위상 배열 레이더의 부엽 차단 시스템 | |
JP4397732B2 (ja) | 測位装置 | |
KR20170128899A (ko) | 무선 통신 신호의 수신 시각 산출 장치 및 방법 | |
JP2009250704A (ja) | 距離測定装置用受信機および距離測定装置 | |
JP6090036B2 (ja) | スペクトル拡散受信機 | |
KR100540602B1 (ko) | 신호 에너지 누적을 이용한 신호 검출 방법 및 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201380008644.9 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13746377 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2013557449 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 13746377 Country of ref document: EP Kind code of ref document: A1 |