WO2023276592A1 - 推定方法、および、推定装置 - Google Patents

推定方法、および、推定装置 Download PDF

Info

Publication number
WO2023276592A1
WO2023276592A1 PCT/JP2022/023087 JP2022023087W WO2023276592A1 WO 2023276592 A1 WO2023276592 A1 WO 2023276592A1 JP 2022023087 W JP2022023087 W JP 2022023087W WO 2023276592 A1 WO2023276592 A1 WO 2023276592A1
Authority
WO
WIPO (PCT)
Prior art keywords
complex transfer
transfer function
combined complex
antenna elements
signal
Prior art date
Application number
PCT/JP2022/023087
Other languages
English (en)
French (fr)
Inventor
武司 中山
翔一 飯塚
尚樹 本間
アブドゥサイミ アブドゥアイニ
信之 白木
Original Assignee
パナソニックＩｐマネジメント株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by パナソニックＩｐマネジメント株式会社 filed Critical パナソニックＩｐマネジメント株式会社
Priority to CN202280044197.1A priority Critical patent/CN117546046A/zh
Priority to EP22832749.0A priority patent/EP4365621A1/en
Priority to JP2023531747A priority patent/JPWO2023276592A1/ja
Publication of WO2023276592A1 publication Critical patent/WO2023276592A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/003—Bistatic radar systems; Multistatic radar systems
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/536—Discriminating between fixed and moving objects or between objects moving at different speeds using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50—Systems of measurement based on relative movement of target
    • G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/56—Discriminating between fixed and moving objects or between objects moving at different speeds for presence detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415—Identification of targets based on measurements of movement associated with the target

Definitions

  • the present disclosure relates to an estimation method and an estimation device for estimating the direction or position of a living body by using radio signals.
  • Patent Document 1 or Patent Document 2 discloses a method of using radio signals as a method of knowing the position of a person.
  • Patent Literature 1 discloses a method of locating a person to be detected by using the frequency response of a complex transfer function of a wireless device.
  • Patent Document 2 discloses a method of detecting the position or state of a person to be detected by using difference information between two complex transfer functions of a radio device at two points in time at predetermined intervals. disclosed.
  • Patent Document 3 discloses a method of performing timing correction between a plurality of sensors using a timer correction command.
  • Patent Literature 4 discloses a method of measuring the degree of fatigue based on the statistical results of two sensors.
  • US Pat. No. 6,201,203 discloses a method for distinguishing between human-induced Doppler shift and clock frequency error between a transmitter and a receiver.
  • the present disclosure has been made in view of the above circumstances, and is compatible with radio equipment that outputs at least part of CSI, specifically SIMO (Single Input Multiple Output) or MISO (Multiple Input Single Output). It is an object of the present invention to provide an estimation method and the like capable of accurately estimating the direction or position of a living body with respect to an estimating device using a wireless device.
  • SIMO Single Input Multiple Output
  • MISO Multiple Input Single Output
  • an estimation method includes a first radio device including a first antenna unit having M (M is a natural number of 2 or more) first antenna elements; and a second radio device comprising a second antenna unit having second antenna elements (N is a natural number of 2 or more).
  • a plurality of first reception signals received by each of the N second antenna elements wherein the first transmission signal is transmitted a plurality of times using one of the first antenna elements of observing a plurality of first received signals including reflected signals of the first transmitted signals transmitted a plurality of times and reflected by the living body during a first period corresponding to the cycle of the movement of the living body, , transmitting a second transmission signal multiple times using one of the N second antenna elements, and a plurality of second receptions received by each of the M first antenna elements; a signal, wherein a plurality of second received signals including reflected signals of the second transmitted signals transmitted a plurality of times and reflected by the living body are observed during the first period of time, and observed during the first period of time a first combination representing propagation characteristics between each of the M first antenna elements and each of the N second antenna elements using the plurality of first received signals and the plurality of second received signals; A plurality of complex transfer functions are calculated, and the effects of the clock frequency error and the transmission power error occurring in the first radio and the second
  • a method for accurately estimating the direction or position of a living body with respect to an estimating device using a wireless device that outputs at least part of CSI can be realized.
  • FIG. 1 is a block diagram showing a configuration of an estimation device and an example of a detection target according to Embodiment 1.
  • FIG. FIG. 2 is a diagram conceptually showing how signal waves are transmitted in the antenna section shown in FIG.
  • FIG. 3 is a block diagram showing an example of a configuration of an estimation device and a detection target according to Embodiment 1.
  • FIG. 4 is a conceptual diagram showing an example of frequency components of received I and Q signals in Embodiment 1.
  • FIG. 5 is a flowchart showing estimation processing of the estimation device according to Embodiment 1.
  • FIG. FIG. 6 is a block diagram showing a configuration of an estimation device and an example of a detection target according to Embodiment 2.
  • FIG. 7 is a block diagram showing a configuration of an estimation device and an example of a detection target according to Embodiment 2.
  • FIG. FIG. 8 is a diagram showing the concept of experiments using the estimation method according to the second embodiment.
  • FIG. 9 is a diagram showing experimental conditions using the estimation method according to the second embodiment.
  • FIG. 10 is a diagram showing experimental results using the estimation method according to the second embodiment.
  • FIG. 11 is a diagram showing another experimental result using the estimation method according to the second embodiment.
  • Patent Document 1 discloses a method of knowing the position or state of a person to be detected by analyzing components including Doppler shift using Fourier transform. More specifically, time changes of the elements of the complex transfer function are recorded, and the time waveform is Fourier transformed. A living body such as a person gives a slight Doppler effect to the reflected wave due to living body activity such as breathing or heartbeat. Therefore, the component containing the Doppler shift contains the human influence. On the other hand, the non-Doppler-shifted components correspond to radio signals unaffected by people, ie reflected waves from fixed objects or direct waves between transmitting and receiving antennas. In view of the above, Patent Literature 1 discloses that the position or state of a person to be detected can be known by analyzing the component including the Doppler shift.
  • Patent Document 2 discloses a method for detecting the position or state of a person to be detected by using difference information between two complex transfer functions of a radio device at two points in time at predetermined intervals. is disclosed.
  • Patent Document 3 discloses a method for correcting the timing by sending a timer correction command between a plurality of sensors, but it is difficult to adjust the phase of the RF signal of the wireless device by command.
  • Patent Document 4 is a technology intended for an acceleration sensor, and it is difficult to apply it to a wireless device.
  • Patent Document 3 describes frequency fluctuations derived from the transmitter and the receiver, such as the clock frequency error between the transmitter and the receiver, and the sampling clock frequency error used in the AD converter and the DA converter.
  • a method for estimating a person's position by suppressing the components caused by the Doppler shift and clock frequency error caused by the person is disclosed.
  • wireless devices that output only a part of CSI, such as M ⁇ 1 elements or 1 ⁇ N elements.
  • a wireless device that outputs only a part of these CSIs cannot acquire CSI having M ⁇ N elements, so there is a problem that it is difficult to estimate the direction or position of a living body.
  • the inventors acquired CSI of the forward and backward paths, and combined the acquired CSI in consideration of phase rotation, as shown in FIGS. 1 and 3. , the correlation function (more specifically, the correlation matrix) of the estimator can be calculated.
  • the inventors have found that, in addition to the clock frequency error described in Patent Document 4, the difference in CSI acquisition time, phase shift, reversal of phase rotation due to forward and backward paths, etc. should be taken into consideration.
  • the inventors found that it is possible to generate a channel by, for example, focusing on the first element of CSI in the forward and backward paths, matching the phase of one and then combining with the other, and , and found a method for estimating the direction or position of a living body using a radio that outputs only a part of CSI.
  • An estimation method includes a first radio device including a first antenna unit having M (M is a natural number of 2 or more) first antenna elements, and N (N is a natural number of 2 or more) ), the estimation method performed by an estimation device comprising a second radio device having a second antenna unit having a second antenna element of ), wherein the measurement target area includes A first transmission signal is transmitted a plurality of times using one first antenna element, and a plurality of first reception signals received by each of the N second antenna elements, wherein the first transmission signal is transmitted a plurality of times.
  • a plurality of first received signals including reflected signals of one transmitted signal reflected by a living body are observed during a first period corresponding to the period of movement of the living body, and the N second antennas are arranged in the area to be measured. transmitting a second transmitted signal multiple times using a second antenna element of one of the elements, and a plurality of second received signals received by each of said M first antenna elements, said multiple transmitted observing, during the first period, a plurality of second received signals including a reflected signal of the second transmitted signal reflected by the living body, and the plurality of first received signals observed during the first period; and using the plurality of second received signals to calculate a plurality of first combined complex transfer functions representing propagation characteristics between each of the M first antenna elements and the N second antenna elements, By subjecting a plurality of the first combined complex transfer functions to a predetermined method for suppressing the effects of clock frequency errors and transmission power errors occurring in the first radio and the second radio, calculating a plurality of second combined complex transfer functions, calculating biological information corresponding to a pre
  • the CSI on the outbound path and the CSI on the inbound path are determined by the clock frequency error between the first radio and the second radio and
  • the direction or position of the living body can be estimated by combining while performing calculations to suppress the influence of the transmission power error.
  • the forward and backward CSI may be a portion of the CSI output by a radio that outputs only a portion of the CSI. Therefore, according to the above estimation method, the direction or position of the living body can be estimated using a wireless device that outputs at least part of the CSI.
  • the first element of the first transmission signal and the first element of the second transmission signal are used to calculate from the plurality of first reception signals and the phase of one of the second complex transfer function calculated from the plurality of second received signals is adjusted to the phase of the other, and then the first complex transfer function and the second complex transfer
  • the first combined complex transfer function may be calculated by combining functions.
  • the first radio when the CSI on the outbound path and the CSI on the return path are combined, by matching the phases of the first complex transfer function and the second complex transfer function before combining, the first radio This contributes to suppressing the effects of clock frequency error and transmission power error between the radio and the second radio. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • a second combined complex transfer function may be calculated in which a component corresponding to at least one of timing fluctuations of analog-to-digital conversion of is suppressed.
  • the clock fluctuation between the transmission section and the reception section of the transmission signal and the timing fluctuation of digital-analog conversion or analog-digital conversion At least one is used to more easily calculate the second combined complex transfer function. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • phase correction by the first element of the first transmission signal and the first element of the second transmission signal is performed on the first transmission signal and the second transmission signal.
  • the second combined complex transfer function may be calculated by applying a phase rotation of one of the first transmitted signal and the second transmitted signal to one of the second transmitted signals to eliminate the phase rotation of the other.
  • an estimation method includes a first radio device including a first antenna unit having M (M is a natural number of 2 or more) first antenna elements, and N (N is 2 or more) natural number) of the M first antenna elements in the area to be measured.
  • a plurality of third received signals including reflected signals in which one transmitted signal is reflected by a living body during a first period corresponding to the cycle of the movement of the living body, transmitting the second transmit signal from each of the elements multiple times, and a plurality of fourth received signals received by one first antenna element of the M first antenna elements, wherein A plurality of fourth received signals, in which the second transmitted signal includes a reflected signal reflected by the living body, are observed during the first period, the plurality of third received signals observed during the first period, and the Using a plurality of fourth received signals, a plurality of first combined complex transfer functions representing propagation characteristics between each of the M first antenna elements and the N second antenna elements are calculated, and a predetermined method, by subjecting the plurality of first combined complex transfer functions to a calculation that suppresses the effects of the clock frequency error and the transmission power error occurring in the first radio device and the second radio device, calculating a second combined complex transfer function, calculating biometric information corresponding to a predetermined frequency range corresponding to a component
  • the CSI of the forward path and the CSI of the return path are the clock frequency error between the first radio and the second radio and
  • the direction or position of the living body can be estimated by combining while performing calculations to suppress the influence of the transmission power error.
  • the forward and backward CSI may be a portion of the CSI output by a radio that outputs only a portion of the CSI. Therefore, according to the above estimation method, the direction or position of the living body can be estimated using a wireless device that outputs at least part of the CSI.
  • the first element of the third received signal and the first element of the fourth received signal are used to calculate from the plurality of third received signals and the phase of one of the first complex transfer function calculated from the plurality of fourth received signals is adjusted to the phase of the other, and then the first complex transfer function and the second complex transfer
  • the first combined complex transfer function may be calculated by combining functions.
  • the first radio when the CSI on the outbound path and the CSI on the return path are combined, by matching the phases of the first complex transfer function and the second complex transfer function before combining, the first radio This contributes to suppressing the effects of clock frequency error and transmission power error between the radio and the second radio. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • a second combined complex transfer function may be calculated in which a component corresponding to at least one of timing fluctuations of analog-to-digital conversion of is suppressed.
  • the clock fluctuation between the transmission section and the reception section of the transmission signal and the timing fluctuation of digital-analog conversion or analog-digital conversion At least one is used to more easily calculate the second combined complex transfer function. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • phase correction by the first element of the third received signal and the first element of the fourth received signal is performed on the third received signal and the fourth received signal.
  • the second combined complex transfer function may be calculated by applying a phase rotation to one of the third received signal and the fourth received signal to eliminate the phase rotation of the other.
  • the second combined complex transfer function extracting one element from the first combined complex transfer function and dividing all the elements of the first combined complex transfer function by the extracted one element. , the second combined complex transfer function may be calculated.
  • an average value of all elements of the first combined complex transfer function is calculated, and all elements of the first combined complex transfer function are divided by the calculated average value. to calculate the second combined complex transfer function.
  • a correlation matrix is calculated by a predetermined method for the first combined complex transfer function, and eigenvalue decomposition is performed on the correlation matrix to obtain one or more eigenvalues and the one or more calculating an eigenvector corresponding to each eigenvalue of and multiplying the first combined complex transfer function by an eigenvector corresponding to a maximum eigenvalue among the one or more calculated eigenvectors, among the one or more calculated eigenvectors
  • the third complex transfer function which is the channel component of the direct wave
  • the eigenvalues and eigenvectors calculated from the first combined complex transfer function can be used to more easily calculate the second combined complex transfer function from the first combined complex transfer function. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • a left singular vector and a right singular vector are calculated by singular value decomposition of the first combined complex transfer function by a predetermined method, and the left singular vector and the right singular vector multiplying the first combined complex transfer function by a vector to calculate a fourth complex transfer function that is a channel component of the direct wave, and multiplying all elements of the first combined complex transfer function by the fourth complex transfer function;
  • the second combined complex transfer function may be calculated by dividing.
  • the singular values and singular vectors calculated from the first combined complex transfer function can be used to more easily calculate the second combined complex transfer function from the first combined complex transfer function. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • the predetermined direction-of-arrival estimation method may be an estimation method that uses at least one of the MUSIC (MUltiple Signal Classification) method, the beamformer method, and the Capon method.
  • MUSIC MUltiple Signal Classification
  • the direction or position of the living body can be more easily estimated by any one of the MUSIC method, the beamformer method, and the Capon method. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • the first received signal or the second received signal may be a signal modulated by an OFDM (Orthogonal Frequency Division Multiplexing) signal.
  • OFDM Orthogonal Frequency Division Multiplexing
  • the direction or position of the living body can be more easily estimated using the first received signal or the second received signal modulated by the OFDM signal. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • the third received signal or the fourth received signal may be a signal modulated by an OFDM (Orthogonal Frequency Division Multiplexing) signal.
  • OFDM Orthogonal Frequency Division Multiplexing
  • the direction or position of the living body can be more easily estimated using the first received signal or the second received signal modulated by the OFDM signal. Therefore, according to the above estimation method, it is possible to more easily estimate the direction or position of a living body using a wireless device that outputs at least part of CSI.
  • an estimating device includes a first radio device including a first antenna unit having M (M is a natural number of 2 or more) first antenna elements, and N (N is 2 or more) and a second radio device having a second antenna unit having second antenna elements of (a natural number of ), wherein the first radio device is arranged in a measurement target area to include the M first antennas transmitting a first transmit signal multiple times using a first antenna element of one of the elements, the second radio receiving a plurality of first received signals received by each of the N second antenna elements; wherein a plurality of first received signals including reflected signals of the first transmitted signals transmitted a plurality of times and reflected by a living body are observed during a first period corresponding to a period of movement of the living body; A second radio device transmits a second transmission signal multiple times using one of the N second antenna elements to the area to be measured, and the first radio device transmits the M a plurality of second received signals received by each of the first antenna elements, wherein the second transmitted
  • a second combined complex transfer function calculating unit for calculating; a biological information calculating unit for extracting biological information corresponding to a predetermined frequency range corresponding to the component related to the living body in the plurality of second combined complex transfer functions; and a position estimation processing unit for estimating the direction or position of the living body using the biological information by a direction estimation method.
  • an estimating device includes a first radio device including a first antenna unit having M (M is a natural number of 2 or more) first antenna elements, and N (N is 2 or more) and a second radio device having a second antenna unit having second antenna elements of (a natural number of ), wherein the first radio device is arranged in a measurement target area to include the M first antennas transmitting a first transmit signal multiple times from each of the N second antenna elements; observing a plurality of third received signals including reflected signals of the first transmitted signals transmitted a plurality of times and reflected by a living body during a first period corresponding to a period of movement of the living body; transmits a second transmission signal multiple times from each of the N second antenna elements to the area to be measured, and the first radio transmits one of the M first antenna elements a plurality of fourth received signals received by the first antenna element, the plurality of fourth received signals including reflected signals reflected by the living body from the second transmitted signals transmitted a plurality of times; period, and the
  • the present disclosure can be implemented not only as a device, but also as an integrated circuit including processing means included in such a device, as a method using processing means that constitute the device as steps, or as steps can be realized as a program to be executed by a computer, or as information, data or a signal indicating the program.
  • These programs, information, data and signals may be distributed via recording media such as CD-ROMs and communication media such as the Internet.
  • the estimation device 10 and the estimation method combine forward CSI and return CSI, extract a direct wave that does not pass through a living body, detect phase rotation caused by frequency error, and remove phase rotation caused by frequency error from the entire channel. By doing so, it becomes possible to estimate the direction in which the living body exists.
  • the estimation apparatus 10 estimates the direction of a living body, which is a detection target, using outbound and inbound CSI obtained from two predetermined wireless devices. .
  • FIG. 1 is a block diagram showing an example of the configuration of estimation device 10 according to the present embodiment.
  • FIG. 1 also shows a living body 30 that is a detection target of the estimation device 10 shown in FIG.
  • the estimation device 10 shown in FIG. A transfer function calculator 24 , a second combined complex transfer function calculator 25 , a biological information calculator 26 , and a position estimation processor 27 are provided.
  • the estimating device 10 estimates the direction in which the living body 30 exists relative to the estimating device 10 .
  • the estimating device 10 is not limited to having all of the above components in one housing.
  • the estimating device 10 may include a plurality of devices communicably connected to each other, and the above components may be separately arranged in each of the plurality of devices. In that case, the estimating device 10 may be called an estimating system.
  • the estimating device 10 uses the first radio 12 having M antenna elements and the second radio 11 having N antenna elements to transmit a known signal from the first radio 12 to the second radio 11
  • FIG. 1 shows a case of acquiring CSI (corresponding to forward CSI) when transmitting .
  • FIG. 3 shows a case where the estimating device 10 acquires CSI (corresponding to return CSI) when a known signal is transmitted from the second radio 11 to the first radio 12 .
  • the estimating apparatus 10 has a full MIMO configuration, but as an example, a method of estimating the position of the living body 30 when CSI with a SIMO configuration is obtained for the outward and return paths will be described.
  • the first radio device 12 generates a high-frequency signal used for estimating the direction of the living body 30 using the internal clock fTX of the first radio device 12 . For example, as shown in FIG. 2 , the first radio device 12 transmits the generated signal (transmission wave) from the antenna element provided in the antenna section 20 . Also, the first radio 12 processes the received signal received by each of the antenna elements during reception.
  • FIG. 4 shows a conceptual diagram showing an example of the frequency components of the I and Q signals received by the estimation device 10 in this embodiment. More specifically, the dashed line 100 shown in FIG. 4 indicates the DC component of the signal received by the estimator 10, the solid line 101 indicates the AC component of the I signal received by the estimator 10, and the dashed line 102 indicates 4 shows the AC component of the Q signal received by the estimator 10;
  • the antenna section 20 has M (M is a natural number of 2 or more) antenna elements (specifically, antenna elements #1 to #M).
  • the first transmission/reception unit 22 transmits or receives a predetermined signal using M antenna elements.
  • the first transmitting/receiving unit 22 transmits a known signal for biometric estimation during transmission.
  • the first transmitting/receiving unit 22 receives a received signal received by each of the antenna elements, which is transmitted from the antenna element on the transmitting side and includes a reflected signal reflected by the living body 30, to the living body 30. Observations are made for a first period corresponding to a cycle derived from 30 activities.
  • the antenna element on the transmitting side may be an antenna element of the antenna section 21 of the second radio device 11 .
  • the living body 30 is the living body 30 as shown in FIG. 2, but is not limited thereto.
  • a cycle derived from the activity of the living body 30 is a living body-derived cycle (biological variation cycle) including at least one of breathing, heartbeat, and body movement of the living body 30 .
  • the first transceiver 22 sends the observed received signal to the first combined complex transfer function calculator 24 as a second received signal.
  • the second radio 11 generates a high-frequency signal used for estimating the direction of the living body 30 using the internal clock fRX like the first radio 12 .
  • the second radio device 11 transmits the generated signal (transmission wave) from the antenna element included in the antenna section 21 .
  • the second wireless device 11 processes the received signal received by each of the antenna elements during reception.
  • the antenna section 21 has N (N is a natural number of 2 or more) antenna elements (specifically, antenna elements #1 to #N).
  • the second transmitting/receiving unit 23 transmits received signals received by each of the N antenna elements, which are transmitted from the antenna elements on the transmitting side and include reflected signals reflected by the living body 30, to the living body 30. Observation is made for the first period corresponding to the cycle derived from the activity.
  • the first period may be the same as the first period during which the first transmitting/receiving section 22 observes the received signal.
  • the second transceiver 23 sends the observed received signal to the first combined complex transfer function calculator 24 as the first received signal.
  • the first transmitting/receiving unit 22 and the second transmitting/receiving unit 23 each have a plurality of antenna elements, and the high-frequency signal received by the corresponding antenna element is transmitted using the internal clock fTX or fRX, Converts to a low-frequency signal that can be processed.
  • the first transmitting/receiving section 22 transmits a known signal for acquiring CSI from one predetermined antenna element of the antenna section 20 .
  • the transmitted signal reflected by the living body 30 is received by the N antenna elements of the antenna section 21 of the second wireless device 11 .
  • the second transmitting/receiving section 23 sends the received signal to the first combined complex transfer function calculating section 24 as a first received signal.
  • the second transmitting/receiving section 23 transmits a known signal for acquiring CSI from one predetermined antenna element of the antenna section 21 .
  • the transmitted signal reflected by the living body 30 is received by the M antenna elements of the antenna section 20 of the first wireless device 12 .
  • the first transmitting/receiving section 22 sends the received signal to the first combined complex transfer function calculating section 24 as a second received signal.
  • the first combined complex transfer function calculator 24 calculates a plurality of complex transfer functions representing propagation characteristics between the first radio device 12 and the second radio device 11 from the first received signal and the second received signal. do.
  • the first combined complex transfer function calculation unit 24 calculates the first radio device 12 and the second A complex transfer function representing the propagation characteristics with the wireless device 11 is calculated.
  • FIG. 2 is a diagram conceptually showing how signal waves are transmitted in the antenna section 21 shown in FIG.
  • transmission waves transmitted from the antenna elements of the antenna section 20 are reflected by the living body 30 and reach the array antenna of the antenna section 21 .
  • the array antenna is a linear array consisting of M antenna elements with an element spacing of d.
  • ⁇ be the direction of the living body 30 viewed from the front of the array antenna. It is assumed that the distance between the living body 30 and the array antenna is sufficiently large, and the reflected wave originating from the living body 30 and arriving at the array antenna can be regarded as a plane wave.
  • the first combined complex transfer function calculator 24 calculates the first received signal vector observed using the array antenna , a first complex transfer function vector representing the propagation characteristics between the antenna element on the transmitting side and the array antenna on the receiving side can be calculated.
  • the first complex transfer function vector h F is, for example, It can be calculated by where s is the complex transmitted signal and is assumed to be known. Note that the first complex transfer function vector h F that is the forward CSI is an example of the first complex transfer function.
  • the first combined complex transfer function calculator 24 calculates the second received signal vector observed using the array antenna as in the case of forward CSI.
  • a second complex transfer function vector representing the propagation characteristics between the antenna element on the transmitting side and the array antenna on the receiving side can be calculated.
  • the second complex transfer function vector h B is, for example, It can be calculated by note that, represents the transpose.
  • h F11 and h B11 are essentially the same value (in other words, by design), but the phases differ due to clock frequency error, power error, and acquisition time difference. (ie, combine the first complex transfer function vector and the second complex transfer function vector).
  • the first combined complex transfer function calculator 24 uses the first element of the first transmission signal and the first element of the second transmission signal to calculate the first complex transfer function calculated from the plurality of first received signals. By matching the phase of one of the function and the second complex transfer function calculated from the plurality of second received signals to the phase of the other, and then combining the first complex transfer function and the second complex transfer function, Compute the first combined complex transfer function.
  • the first combined complex transfer function vector h(t) is expressed as follows.
  • the first combined complex transfer function vector h(t) includes frequency variation components from the transmitter and receiver, and Doppler shift from the living body 30 .
  • Transmitter- and receiver-derived frequency variation components include, for example, (i) attenuation or phase rotation due to spatial propagation of the transmitted signal, (ii) clock frequency error between the transmitter and receiver (fRX-fTX), (iii) ) Sampling clock frequency error and the like used in the radio such as DA conversion are included.
  • the second combined complex transfer function calculator 25 calculates h(t) Extract any one element h l in .
  • the second combined complex transfer function calculator 25 divides all elements of the first combined complex transfer function vector h(t) as direct wave components by one element h l extracted as described above.
  • a second combined complex transfer function vector h' is calculated by removing the phase rotation.
  • any element of the first combined complex transfer function vector h(t), such as the element h1 may be used as the element of the direct wave component.
  • the second combined complex transfer function vector h' is an example of the second combined complex transfer function.
  • the second combined complex transfer function calculation unit 25 performs a predetermined operation using one or more elements of the first combined complex transfer function vector h(t) to obtain the first combined complex transfer function vector h From (t), (1) clock fluctuation between the first radio 12 and the second radio 11, and (2) timing fluctuation of digital-analog conversion of the transmission signal or analog-digital conversion of the reception signal. Compute a second combined complex transfer function vector h' in which at least one corresponding component is suppressed.
  • the biological information calculation unit 26 uses the calculated second combined complex transfer function vector h′ to calculate two second combined complex transfer functions at two points in time at predetermined intervals in the same manner as the difference information calculation unit described in Patent Document 1.
  • Biological information is calculated by calculating the difference information of the vector h'. That is, the biological information calculation unit 26 extracts the biological information corresponding to the components related to the living body 30 by extracting the biological information corresponding to the predetermined frequency range in the calculated plurality of second combined complex transfer function vectors h'. .
  • the biometric information calculator 26 extracts biometric information corresponding to components affected by vital activity including at least one of respiration, heartbeat, and body movement of the living body 30 .
  • this biological information can be obtained by calculating the biological component from the second combined complex transfer function vector h', and it goes without saying that the same effect can be obtained by using the frequency response.
  • the number of differential values (difference information) of the second complex transfer function vector h' corresponding to the antenna section 21 is also plural. These are collectively defined as a complex difference channel vector. If the number of antenna elements is N, the complex differential channel vector is and is. Also, l or m is a positive integer representing a measurement number and a sample time.
  • the position estimation processing unit 27 uses the biological information extracted from the second combined complex transfer function vector h′ to estimate the direction in which the living body 30 exists, using the estimation device 10 as a direction reference, by a predetermined direction-of-arrival estimation method. do. More specifically, the predetermined direction-of-arrival estimation method may be an estimation method based on the MUSIC (MUltipleSignal Classification) algorithm, such as the direction estimation processing unit described in Patent Document 1, and this case will be described as an example. An estimation method based on the beamformer method or the Capon method may be used.
  • MUSIC MUltipleSignal Classification
  • aF is the forward steering vector
  • aB is the return steering vector
  • the signal transmitted from the first radio 12 may be a continuous signal (CW signal) or an encoded signal such as an OFDM (Orthogonal Frequency Division Multiplexing) signal.
  • CW signal continuous signal
  • OFDM Orthogonal Frequency Division Multiplexing
  • the second combined complex transfer function calculator 25 of the present embodiment uses the first combined complex transfer function vector h(t) as a direct wave component for any given value in the first combined complex transfer function vector h(t).
  • the second combined complex transfer function vector h' is calculated by dividing by one element hl , the division is not limited to one element hl .
  • the second combined complex transfer function calculator 25 adds any one element h l , the element h lmin with the smallest time variation of
  • the average value of all elements of the first complex transfer function vector h0 may be used.
  • singular value decomposition of the complex transfer function at an arbitrary time t is performed to obtain a singular vector, and thereby the second combined complex transfer function vector h' is calculated. good.
  • the complex transfer function may be observed for a certain period of time, and the correlation matrix for the entire observation time may be subjected to eigenvalue decomposition to obtain eigenvectors, thereby calculating the second combined complex transfer function vector h'.
  • the direct wave component may be obtained by eigenvalue decomposition of the correlation matrix for the entire observation time to obtain an eigenvector, thereby calculating the second combined complex transfer function vector h'.
  • the first combined complex transfer function calculation unit 24 uses the first complex transfer function vector h F and the second complex transfer function vector h B , and utilizes the fact that the phase rotation appears to be opposite between the forward and return paths. , which gives A second complex transfer function vector may be calculated as This eliminates the phase rotation and is equivalent to observing at an intermediate frequency between the two stations.
  • the first radio device 12, the second radio device 11, the first combined complex transfer function calculation unit 24, the second combined complex transfer function calculation unit 25, the biological information calculation unit 26, and the position estimation processing unit 27 are It may be implemented by one or more processors executing a program stored in a memory provided by, or may be implemented by one or more dedicated circuits. That is, the above components may be realized by software or by hardware.
  • FIG. 5 is a flow chart showing estimation processing of the estimation device 10 according to the present embodiment.
  • the estimation device 10 transmits transmission signals to the area to be measured using the first radio 12 and the second radio 11, respectively, and observes the received signals during the first period (S10). More specifically, the estimating device 10 observes a received signal including a reflected signal transmitted from one antenna element and reflected by the living body 30 for a first period corresponding to a period derived from the activity of the living body 30. do. Observation is made for each of the outward route and the return route.
  • step S10 the estimating apparatus 10 obtains the propagation characteristics between one transmitting antenna element and each of the receiving antenna elements from a plurality of received signals observed in the first period.
  • a plurality of combined complex transfer functions are calculated (S20). Since the details are as described above, the description here is omitted. The same applies to the following.
  • the estimation device 10 divides the first combined complex transfer function calculated in step S20 by the direct wave component to calculate the second combined complex transfer function (S30).
  • the estimation device 10 calculates biological information from the second combined complex transfer function calculated in step S30 (S40).
  • the estimating device 10 uses the two or more pieces of biological information calculated in step S40 to estimate the position where the living body 30 exists (S50).
  • the estimation apparatus 10 and the estimation method of the present embodiment the direct wave component that does not pass through the living body 30 is extracted from the received signal, the phase rotation caused by the frequency error is detected, and the Since the clock fluctuation and the phase rotation caused by the frequency error are removed, it is possible to accurately estimate the direction in which the living body 30 exists. Therefore, the estimating device 10 can accurately estimate the direction in which the living body 30 exists relative to the estimating device 10 .
  • an estimation device 10A and an estimation method for estimating the direction in which a living body exists by combining forward CSI and return CSI will be described, even if the CSI is a MISO configuration.
  • FIG. 6 and 7 are block diagrams each showing an example of the configuration of estimation device 10A according to the present embodiment. Elements similar to those in FIG. 1 are denoted by the same reference numerals, and detailed description thereof is omitted. In the present embodiment, the same elements as in the first embodiment have the same operations and variations, unless otherwise mentioned, and redundant descriptions are omitted.
  • the estimating device 10A shown in FIG. A function calculator 24A, a second combined complex transfer function calculator 25, a biological information calculator 26, and a position estimation processor 27 are provided.
  • the estimating device 10A estimates the direction in which the living body 30 exists relative to the estimating device 10A.
  • the estimating device 10A uses the first radio 12 having M antenna elements and the second radio 11 having N antenna elements to transmit known signals from the first radio 12 to the second radio 11.
  • FIG. 6 shows a case of acquiring CSI (corresponding to forward CSI) when transmitting .
  • FIG. 7 shows a case where the estimating device 10A acquires CSI (corresponding to return CSI) when a known signal is transmitted from the second radio device 11 to the first radio device 12.
  • the estimating device 10A has a full MIMO configuration. As an example, a method of estimating the position of the living body 30 by the estimating device 10A when CSI with a MISO configuration is obtained for the forward and return paths will be described. do.
  • the first transmission/reception unit 22A transmits or receives a predetermined signal using M antenna elements.
  • the first transmission/reception unit 22A transmits a known signal for biometric estimation during transmission.
  • the first transmission/reception unit 22A receives a received signal received by each of the antenna elements, which is transmitted from the antenna element on the transmitting side and includes a reflected signal reflected by the living body 30, to the living body 30. A first period corresponding to a cycle derived from 30 activities is observed.
  • the first transmitting/receiving section 22A sends the observed received signal as a fourth received signal to the first combined complex transfer function calculating section 24A.
  • the second transmitting/receiving unit 23A transmits received signals received by each of the N antenna elements, which are transmitted from the antenna elements on the transmitting side and include reflected signals reflected by the living body 30, to the living body 30. Observation is made for the first period corresponding to the cycle derived from the activity.
  • the second transmitting/receiving section 23A sends the observed received signal as a third received signal to the first combined complex transfer function calculating section 24A.
  • the first transmitting/receiving unit 22A transmits known signals for acquiring CSI from the M antenna elements of the antenna unit 20, and the signals reflected by the living body 30 are is received by one predetermined antenna element of the antenna section 21 of the second radio 11 .
  • the second transmitting/receiving section 23A sends the received signal as a third received signal to the first combined complex transfer function calculating section 24A.
  • the second transmitting/receiving section 23A transmits known signals for acquiring CSI from the N antenna elements of the antenna section 21 .
  • a predetermined single antenna element of the antenna section 20 of the first wireless device 12 receives the transmitted signal reflected by the living body 30 .
  • the first transmission/reception section 22A sends the received signal as a fourth reception signal to the first combined complex transfer function calculation section 24A.
  • the first combined complex transfer function calculator 24A calculates a plurality of complex transfer functions representing propagation characteristics between the first radio device 12 and the second radio device 11 from the third received signal and the fourth received signal. do.
  • the first combined complex transfer function calculation unit 24A calculates the first wireless device 12 and the second A complex transfer function representing the propagation characteristics with the wireless device 11 is calculated.
  • the first combined complex transfer function calculator 24A calculates the third received signal vector observed using the array antenna , it is possible to calculate the third complex transfer function vector representing the propagation characteristics between the antenna element on the transmitting side and the array antenna on the receiving side.
  • the third complex transfer function vector h F is, for example, It can be calculated by
  • the first combined complex transfer function calculator 24A calculates the fourth received signal vector , it is possible to calculate a fourth complex transfer function vector representing the propagation characteristics between the antenna element on the transmitting side and the array antenna on the receiving side.
  • the fourth complex transfer function vector h B is, for example, It can be calculated by
  • h F11 and h B11 are essentially the same value (in other words, by design), but the phases differ due to clock frequency error, power error, and acquisition time difference. (ie, combine the first complex transfer function vector and the second complex transfer function vector).
  • the first combined complex transfer function calculator 24A uses the first element of the third received signal and the first element of the fourth received signal to calculate the first complex transfer calculated from the plurality of third received signals. By matching the phase of one of the function and the second complex transfer function calculated from the plurality of fourth received signals to the phase of the other, and then combining the first complex transfer function and the second complex transfer function, Compute the first combined complex transfer function.
  • the first combined complex transfer function vector h(t) of the SIMO configuration is expressed as follows.
  • the first combined complex transfer function vector h(t) includes frequency fluctuation components derived from the transmitter and receiver and Doppler shift derived from the living body 30 .
  • Frequency fluctuation components originating from the transmitter and receiver may be removed by the same method as described in the first embodiment. Also, after the above removal, the method for estimating the position of the living body 30 is the same as the method described in the first embodiment.
  • the direction in which the living body 30 exists can be estimated by combining the forward CSI and the backward CSI even with MISO-configured CSI. becomes possible.
  • FIG. 8 is a diagram showing the concept of a simulation experiment using the estimation method according to this embodiment.
  • FIG. 9 is a diagram showing experimental conditions using the estimation method according to the present embodiment.
  • Both the transmitting array antenna (Tx array) and the receiving array antenna (Rx array) shown in FIG. 8 have a 4 ⁇ 4 MIMO (Multiple Input Multiple Output) configuration using a four-element patch array antenna. It also has an SP4T (Single-Pole-4-Throw) switch on the transmission side and four receivers on the reception side.
  • Tx array transmitting array antenna
  • Rx array receiving array antenna
  • the distance between the array elements of the transmitting and receiving antennas is 0.5 wavelength
  • the distance D between the transmitting and receiving antennas is 4.0 m
  • the antenna height h is the chest height of a human (Living-Body) standing upright. It was set to a certain 1.0m.
  • An unmodulated continuous wave (CW: Continuous Wave) of 2.47125 GHz was transmitted from the transmitter, and the channel measurement time was 33 seconds. During the channel measurement, there was no one other than the subject, and the subject faced the wall on the antenna side.
  • CW Continuous Wave
  • FIG. 10 is a diagram showing simulation results using the estimation method according to the second embodiment.
  • (a) of FIG. 10 shows the results of a conventional simulation using full MIMO CSI.
  • FIG. 10(b) shows the result of a conventional simulation using CSI with a SIMO configuration.
  • Fig. 10 shows the results of body position estimation when there are two subjects in the space.
  • FIG. 11 is a diagram showing another experimental result using the estimation method according to the second embodiment.
  • FIG. 11 shows a cumulative probability distribution (CDF: Cumulative Distribution Function) of position estimation errors when there are two subjects.
  • CDF Cumulative Distribution Function
  • the horizontal axis indicates the position estimation error (unit: m)
  • the vertical axis indicates the CDF for the position estimation error.
  • FIG. 11(a) shows the CDF when full MIMO CSI is used
  • FIG. 11(b) shows the CDF when SIMO configured CSI is used.
  • the direction in which the living body exists can be determined by acquiring the CSI of the forward path and the return path, combining them, and processing the frequency error in the same manner. It is possible to estimate
  • Embodiments 1 and 2 direction estimation or position estimation of the living body 30 has been described as an example, but the subject of estimation processing is not limited to the living body 30 .
  • the target of the estimation process is applicable to various living organisms (machines, etc.) that, when irradiated with high-frequency signals, give a Doppler effect to the reflected waves due to their activity or movement.
  • the present disclosure can be implemented not only as an estimating device having such characteristic components, but also as an estimating method or the like in which the characteristic components included in the estimating device are used as steps. can. It can also be implemented as a computer program that causes a computer to execute each characteristic step included in such a method. It goes without saying that such a computer program can be distributed via a non-temporary computer-readable recording medium such as a CD-ROM or a communication network such as the Internet.
  • the present disclosure can be used for an estimation device that estimates the direction and position of a living body using radio signals. It can be used for an estimation device installed in a home appliance that performs detection, a monitoring device that detects the intrusion of a living body, and the like.

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

推定方法は、第一無線機のＭ個の第一アンテナ素子のうちの一つを用いて第一送信信号を送信し、第一送信信号が生体によって反射された反射信号を含む第一受信信号を第二無線機のＮ個の第二アンテナ素子で観測し、Ｎ個の第二アンテナ素子のうちの一つを用いて第二送信信号を送信し、第二送信信号が生体によって反射された反射信号を含む第二受信信号を観測し（Ｓ１０）、Ｍ個の第一アンテナ素子とＮ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を算出し（Ｓ２０）、第一結合複素伝達関数に、第一無線機と第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで第二結合複素伝達関数を算出し（Ｓ３０）、複数の第二結合複素伝達関数における生体に関する成分に対応する所定周波数範囲に対応する生体情報を抽出し（Ｓ４０）、生体情報を用いて生体の存在する方向または位置を推定する（Ｓ５０）。

Description

推定方法、および、推定装置
　本開示は、無線信号を利用することで生体の方向または位置を推定する推定方法、および、推定装置に関する。
　特許文献１または特許文献２では、人物の位置などを知る方法として、無線信号を利用する方法が開示されている。特許文献１には、無線機の複素伝達関数の周波数応答を利用して検出対象となる人物の位置を知る方法が開示されている。
　また、特許文献２には、無線機の複素伝達関数のうち、所定の間隔の２つの時点における２つの複素伝達関数の差分情報を利用して検出対象となる人物の位置または状態を知る方法が開示されている。
　また、特許文献３には、複数のセンサ間でタイマー補正命令によりタイミング補正を行う方法が開示されている。特許文献４には、２つのセンサの統計結果に基づいて疲労度を測定する方法が開示されている。特許文献５には、人により発生したドップラシフトと、送信機と受信機との間のクロック周波数誤差とを区別する方法が開示されている。
特開２０１５－７２１７３号公報 特開２０１７－１２９５５８号公報 特開２０１４－２１６７８６号公報 特開２０１１－０１９８４５号公報 特開２０２０－１０９３８９号公報
　しかしながら、特許文献１～５に示される技術では、例えばＭ本のアンテナを持つ送信機とＮ本のアンテナを持つ受信機とによって構成される推定装置に対する生体の方向または位置を推定するには、Ｍ×Ｎの要素をもつチャネル状態情報（ＣＳＩ：Ｃｈａｎｎｅｌ　Ｓｔａｔｅ　Ｉｎｆｏｒｍａｔｉｏｎ、以下ＣＳＩと呼称する）の取得が必要である。
　昨今の無線機では、例えば、Ｍ×１、または、１×ＮといったＣＳＩの一部だけを出力する無線機が存在する。これらのＣＳＩの一部だけを出力する無線機では、生体の方向または位置を推定することが難しいという問題がある。
　本開示は、上述の事情を鑑みてなされたもので、ＣＳＩの一部を少なくとも出力する無線機、具体的には、ＳＩＭＯ（Ｓｉｎｇｌｅ　Ｉｎｐｕｔ　Ｍｕｌｔｉｐｌｅ　Ｏｕｔｐｕｔ）またはＭＩＳＯ（Ｍｕｌｔｉｐｌｅ　Ｉｎｐｕｔ　Ｓｉｎｇｌｅ　Ｏｕｔｐｕｔ）などに対応した無線機を用いて、推定装置に対する生体の方向または位置を精度よく推定することができる推定方法などを提供することを目的とする。
　上記目的を達成するために、本開示の一形態に係る推定方法は、Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置が実行する推定方法であって、測定対象の領域に、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子を用いて第一送信信号を複数回送信し、前記Ｎ個の第二アンテナ素子のそれぞれにより受信された複数の第一受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第一受信信号を、前記生体の動きの周期に相当する第一期間に観測し、測定対象の領域に、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子を用いて第二送信信号を複数回送信し、前記Ｍ個の第一アンテナ素子のそれぞれにより受信された複数の第二受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第二受信信号を、前記第一期間に観測し、前記第一期間に観測された前記複数の第一受信信号、および前記複数の第二受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出し、所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出し、複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を算出し、所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する推定方法である。
　本開示によれば、ＣＳＩの一部を少なくとも出力する無線機を用いて推定装置に対する生体の方向または位置を精度よく推定する方法などを実現できる。
図１は、実施の形態１における推定装置の構成、および、検出対象の一例を示すブロック図である。 図２は、図１に示すアンテナ部における信号波の伝達の様子を概念的に示す図である。 図３は、実施の形態１における推定装置の構成、および、検出対象の一例を示すブロック図である。 図４は、実施の形態１における受信したＩ信号およびＱ信号の周波数成分の一例を示す概念図である。 図５は、実施の形態１における推定装置の推定処理を示すフローチャートである。 図６は、実施の形態２における推定装置の構成、および、検出対象の一例を示すブロック図である。 図７は、実施の形態２における推定装置の構成、および、検出対象の一例を示すブロック図である。 図８は、実施の形態２に係る推定方法を用いた実験の概念を示す図である。 図９は、実施の形態２に係る推定方法を用いた実験の条件を示す図である。 図１０は、実施の形態２に係る推定方法を用いた実験結果を示す図である。 図１１は、実施の形態２に係る推定方法を用いた別の実験結果を示す図である。
　（本開示の基礎となった知見）
　人物の位置などを知る方法として、無線信号を利用する方法が検討されている。
　例えば特許文献１には、フーリエ変換を用いてドップラシフトを含む成分を解析することで検出対象となる人物の位置または状態を知る方法が開示されている。より具体的には、複素伝達関数の要素の時間変化を記録し、その時間波形をフーリエ変換する。人物などの生体は、呼吸または心拍などの生体活動により、反射波に僅かなドップラ効果を与える。したがって、ドップラシフトを含む成分は、人物の影響を含んでいる。一方、ドップラシフトのない成分は、人物の影響を受けていない無線信号、つまり固定物からの反射波または送受信アンテナ間の直接波に対応する。以上のことから、特許文献１では、ドップラシフトを含む成分を解析することで検出対象となる人物の位置または状態を知ることができることが開示されている。
　また、特許文献２には、無線機の複素伝達関数のうち、所定の間隔の２つの時点における２つの複素伝達関数の差分情報を利用して、検出対象となる人物の位置または状態を知る方法が開示されている。
　また、特許文献３には、複数のセンサ間にタイマー補正命令を送りタイミング補正する方法が開示されているが、無線機のＲＦ信号の位相調整をコマンドで行うのは難しい。また、特許文献４は、加速度センサを対象とした技術であって、無線機に適用することは難しい。
　また、特許文献３には、送信機と受信機との間のクロック周波数誤差、ならびに、ＡＤ変換機およびＤＡ変換機内で使われるサンプリングクロック周波数誤差など、送信機および受信機由来の周波数変動に関して、人により発生したドップラシフトとクロック周波数誤差により発生した成分を抑制して人の位置を推定する方法が開示されている。
　しかしながら、特許文献１～５に示される技術では、例えばＭ本のアンテナを持つ送信機とＮ本のアンテナを持つ受信機とによって構成される推定装置に対する生体の方向または位置を推定するには、Ｍ×Ｎの要素をもつＣＳＩの取得が必要である。
　昨今の無線機では、例えばＭ×１の要素、または、１×Ｎの要素というように、ＣＳＩの一部だけを出力する無線機が存在する。これらのＣＳＩの一部だけを出力する無線機では、Ｍ×Ｎの要素をもつＣＳＩを取得することができないので、生体の方向または位置を推定することが難しいという問題がある。
　発明者らは、上記問題に対して研究を重ねた結果、図１および図３に示すように、往路と復路とのＣＳＩを取得し、取得したＣＳＩを位相回転に配慮して結合することで、推定装置の相関関数（より具体的には相関行列）を算出できることを見出した。ここでは、発明者らは、特許文献４記載のクロック周波数誤差に加え、ＣＳＩ取得時刻の違いまたは位相ズレ、往路および復路による位相回転の反転などを考慮すればよいことを見出した。これにより、発明者らは、例えば、往路と復路とのＣＳＩの第一エレメントに着目し、一方の位相に合わせてから他方と結合することで、チャネルを生成することが可能であること、そして、ＣＳＩの一部だけを出力する無線機を用いて生体の方向または位置を推定する方法等を見出すことに至った。
　本開示の一態様に係る推定方法は、Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置が実行する推定方法であって、測定対象の領域に、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子を用いて第一送信信号を複数回送信し、前記Ｎ個の第二アンテナ素子のそれぞれにより受信された複数の第一受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第一受信信号を、前記生体の動きの周期に相当する第一期間に観測し、測定対象の領域に、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子を用いて第二送信信号を複数回送信し、前記Ｍ個の第一アンテナ素子のそれぞれにより受信された複数の第二受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第二受信信号を、前記第一期間に観測し、前記第一期間に観測された前記複数の第一受信信号、および前記複数の第二受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出し、所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出し、複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を算出し、所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する推定方法である。
　上記態様によれば、１つのアンテナ素子により送信した信号を複数のアンテナ素子により受信する場合における、往路のＣＳＩと復路のＣＳＩとを、第一無線機と第二無線機とのクロック周波数誤差および送信電力誤差の影響を抑制する演算を施しながら結合することで、生体の存在する方向または位置を推定することができる。往路と復路とのＣＳＩは、ＣＳＩの一部だけを出力する無線機によって出力された、ＣＳＩの一部であってよい。よって、上記推定方法によれば、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第一結合複素伝達関数を算出する際には、前記第一送信信号の第一エレメントと前記第二送信信号の第一エレメントとを用いて、前記複数の第一受信信号から算出される第一複素伝達関数と、前記複数の第二受信信号から算出される第二複素伝達関数との一方の位相を他方の位相に合わせてから、前記第一複素伝達関数と前記第二複素伝達関数とを結合することで、前記第一結合複素伝達関数を算出してもよい。
　上記態様によれば、往路のＣＳＩと復路のＣＳＩとの結合のときに、第一複素伝達関数と第二複素伝達関数との位相を合わせてから結合することによって、より容易に、第一無線機と第二無線機とのクロック周波数誤差および送信電力誤差の影響を抑制することに寄与する。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数を算出する際には、前記第一結合複素伝達関数の１以上の要素を用いて所定の演算を行うことにより、前記第一結合複素伝達関数から、（１）アンテナ部から送信する送信信号を生成する送信部と受信部との間のクロック揺らぎ、および、（２）前記送信信号のデジタル－アナログ変換、または、前記第一受信信号もしくは前記第二受信信号のアナログ－デジタル変換のタイミング揺らぎ、の少なくとも一方に対応する成分が抑制された第二結合複素伝達関数を算出してもよい。
　上記態様によれば、往路のＣＳＩと復路のＣＳＩとの結合のときに、送信信号の送信部と受信部との間のクロック揺らぎ、および、デジタル－アナログ変換またはアナログ－デジタル変換のタイミング揺らぎの少なくとも一方を用いて、より容易に、第二結合複素伝達関数を算出する。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数を算出する際には、前記第一送信信号の第一エレメントと前記第二送信信号の第一エレメントとによる位相補正を前記第一送信信号と前記第二送信信号とに施すことによって、前記第一送信信号と前記第二送信信号との一方に対する他方の位相回転を消失させることで、前記第二結合複素伝達関数を算出してもよい。
　上記態様によれば、第一送信信号と第二送信信号との一方に対する他方の位相回転を消失させることで、より容易に、第二結合複素伝達関数を算出することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　また、本開示の一態様に係る推定方法は、Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置が実行する推定方法であって、測定対象の領域に、前記Ｍ個の第一アンテナ素子のそれぞれより第一送信信号を複数回送信し、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子により受信された複数の第三受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第三受信信号を、前記生体の動きの周期に相当する第一期間に観測し、測定対象の領域に、前記Ｎ個の第二アンテナ素子のそれぞれより第二送信信号を複数回送信し、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子により受信された複数の第四受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第四受信信号を、前記第一期間に観測し、前記第一期間に観測された前記複数の第三受信信号、および前記複数の第四受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出し、所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出し、複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を算出し、所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する推定方法である。
　上記態様によれば、複数のアンテナ素子により送信した信号を１つのアンテナ素子により受信する場合における、往路のＣＳＩと復路のＣＳＩとを、第一無線機と第二無線機とのクロック周波数誤差および送信電力誤差の影響を抑制する演算を施しながら結合することで、生体の存在する方向または位置を推定することができる。往路と復路とのＣＳＩは、ＣＳＩの一部だけを出力する無線機によって出力された、ＣＳＩの一部であってよい。よって、上記推定方法によれば、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第一結合複素伝達関数を算出する際には、前記第三受信信号の第一エレメントと前記第四受信信号の第一エレメントとを用いて、前記複数の第三受信信号から算出される第一複素伝達関数と、前記複数の第四受信信号から算出される第二複素伝達関数との一方の位相を他方の位相に合わせてから、前記第一複素伝達関数と前記第二複素伝達関数とを結合することで、前記第一結合複素伝達関数を算出してもよい。
　上記態様によれば、往路のＣＳＩと復路のＣＳＩとの結合のときに、第一複素伝達関数と第二複素伝達関数との位相を合わせてから結合することによって、より容易に、第一無線機と第二無線機とのクロック周波数誤差および送信電力誤差の影響を抑制することに寄与する。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数を算出する際には、前記第一結合複素伝達関数の１以上の要素を用いて所定の演算を行うことにより、前記第一結合複素伝達関数から、（１）アンテナ部から送信する送信信号を生成する送信部と受信部との間のクロック揺らぎ、および、（２）前記送信信号のデジタル－アナログ変換、または、前記第三受信信号もしくは前記第四受信信号のアナログ－デジタル変換のタイミング揺らぎ、の少なくとも一方に対応する成分が抑制された第二結合複素伝達関数を算出してもよい。
　上記態様によれば、往路のＣＳＩと復路のＣＳＩとの結合のときに、送信信号の送信部と受信部との間のクロック揺らぎ、および、デジタル－アナログ変換またはアナログ－デジタル変換のタイミング揺らぎの少なくとも一方を用いて、より容易に、第二結合複素伝達関数を算出する。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数を算出する際には、前記第三受信信号の第一エレメントと前記第四受信信号の第一エレメントとによる位相補正を前記第三受信信号と前記第四受信信号とに施すことによって、前記第三受信信号と前記第四受信信号との一方に対する他方の位相回転を消失させることで、前記第二結合複素伝達関数を算出してもよい。
　上記態様によれば、第一送信信号と第二送信信号との一方に対する他方の位相回転を消失させることで、より容易に、第二結合複素伝達関数を算出することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数の算出では、前記第一結合複素伝達関数から一つの要素を抽出し、抽出した前記一つの要素で、前記第一結合複素伝達関数の全要素を除算することで、前記第二結合複素伝達関数を算出してもよい。
　上記態様によれば、第一結合複素伝達関数から抽出した一つの要素による除算を用いて、より容易に、第一結合複素伝達関数から第二結合複素伝達関数を算出することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数の算出では、前記第一結合複素伝達関数の全要素の平均値を算出し、算出した前記平均値で、前記第一結合複素伝達関数の全要素を除算することで、前記第二結合複素伝達関数を算出してもよい。
　上記態様によれば、第一結合複素伝達関数の全要素の平均値による除算を用いて、より容易に、第一結合複素伝達関数から第二結合複素伝達関数を算出することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数の算出では、前記第一結合複素伝達関数について所定の方法で相関行列を算出し、前記相関行列を固有値分解することで、１以上の固有値と、前記１以上の固有値それぞれに対応する固有ベクトルを算出し、算出した１以上の前記固有ベクトルのうち、算出した前記１以上の固有値のうちの最大の固有値に対応する固有ベクトルを前記第一結合複素伝達関数に乗算することで、直接波のチャネル成分である第三複素伝達関数を算出し、前記第一結合複素伝達関数の全要素を前記第三複素伝達関数にて除算することで、前記第二結合複素伝達関数を算出してもよい。
　上記態様によれば、第一結合複素伝達関数から算出した固有値および固有ベクトルを用いて、より容易に、第一結合複素伝達関数から第二結合複素伝達関数を算出することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第二結合複素伝達関数の算出では、前記第一結合複素伝達関数について所定の方法で特異値分解することで左特異ベクトルおよび右特異ベクトルを算出し、前記左特異ベクトルおよび前記右特異ベクトルを前記第一結合複素伝達関数に乗算することで、直接波のチャネル成分である第四複素伝達関数を算出し、前記第一結合複素伝達関数の全要素を前記第四複素伝達関数にて除算することで、前記第二結合複素伝達関数を算出してもよい。
　上記態様によれば、第一結合複素伝達関数から算出した特異値および特異ベクトルを用いて、より容易に、第一結合複素伝達関数から第二結合複素伝達関数を算出することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記所定の到来方向推定手法は、ＭＵＳＩＣ（ＭＵｌｔｉｐｌｅＳＩｇｎａｌ　Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法、ビームフォーマ法、およびＣａｐｏｎ法のいずれかを少なくとも用いる推定手法であってもよい。
　上記態様によれば、ＭＵＳＩＣ法、ビームフォーマ法、およびＣａｐｏｎ法のいずれかにより、生体の方向または位置を、より容易に推定することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第一受信信号または前記第二受信信号は、ＯＦＤＭ（Ｏｒｔｈｏｇｏｎａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ）信号に変調されている信号であってもよい。
　上記態様によれば、ＯＦＤＭ信号に変調された第一受信信号または第二受信信号を用いて、より容易に、生体の方向または位置を推定することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　例えば、前記第三受信信号または前記第四受信信号は、ＯＦＤＭ（Ｏｒｔｈｏｇｏｎａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ）信号に変調されている信号であってもよい。
　上記態様によれば、ＯＦＤＭ信号に変調された第一受信信号または第二受信信号を用いて、より容易に、生体の方向または位置を推定することができる。よって、上記推定方法によれば、より容易に、ＣＳＩの一部を少なくとも出力する無線機を用いて生体の方向または位置を推定することができる。
　また、本開示の一態様に係る推定装置は、Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置であって、前記第一無線機が、測定対象の領域に、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子を用いて第一送信信号を複数回送信し、前記第二無線機が、前記Ｎ個の第二アンテナ素子のそれぞれにより受信された複数の第一受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第一受信信号を、前記生体の動きの周期に相当する第一期間に観測し、前記第二無線機が、測定対象の領域に、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子を用いて第二送信信号を複数回送信し、前記第一無線機が、前記Ｍ個の第一アンテナ素子のそれぞれにより受信された複数の第二受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第二受信信号を、前記第一期間に観測し、前記推定装置は、さらに、前記第一期間に観測された前記複数の第一受信信号、および前記複数の第二受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出する第一結合複素伝達関数算出部と、所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出する第二結合複素伝達関数算出部と、複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を抽出する生体情報算出部と、所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する位置推定処理部とを備える推定装置である。
　上記態様によれば、上記推定方法と同様の効果を奏する。
　また、本開示の一態様に係る推定装置は、Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置であって、前記第一無線機が、測定対象の領域に、前記Ｍ個の第一アンテナ素子のそれぞれより第一送信信号を複数回送信し、前記第二無線機が、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子により受信された複数の第三受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第三受信信号を、前記生体の動きの周期に相当する第一期間に観測し、前記第二無線機が、測定対象の領域に、前記Ｎ個の第二アンテナ素子のそれぞれより第二送信信号を複数回送信し、前記第一無線機が、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子により受信された複数の第四受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第四受信信号を、前記第一期間に観測し、前記推定装置は、さらに、前記第一期間に観測された前記複数の第三受信信号、および前記複数の第四受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出する第一結合複素伝達関数算出部と、所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出する第二結合複素伝達関数算出部と、複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を抽出する生体情報算出部と、所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する位置推定処理部とを備える推定装置である。
　上記態様によれば、上記推定方法と同様の効果を奏する。
　なお、本開示は、装置として実現するだけでなく、このような装置が備える処理手段を備える集積回路として実現したり、その装置を構成する処理手段をステップとする方法として実現したり、それらステップをコンピュータに実行させるプログラムとして実現したり、そのプログラムを示す情報、データまたは信号として実現したりすることもできる。そして、それらプログラム、情報、データおよび信号は、ＣＤ－ＲＯＭ等の記録媒体やインターネット等の通信媒体を介して配信してもよい。
　以下、本開示の実施の形態について、図面を用いて詳細に説明する。なお、以下で説明する実施の形態は、いずれも本開示の好ましい一具体例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、本開示を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、本開示の最上位概念を示す独立請求項に記載されていない構成要素については、より好ましい形態を構成する任意の構成要素として説明される。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
　（実施の形態１）
　本実施の形態では、ＣＳＩの一部を少なくとも出力する無線機を用いて推定装置に対する生体の方向または位置を精度よく推定する推定装置および推定方法について説明する。
　具体的には、本実施の形態では、ＳＩＭＯ構成のＣＳＩであっても、往路のＣＳＩと復路のＣＳＩとを結合することを用いて、生体の存在する方向を推定する推定装置１０および推定方法を説明する。推定装置１０および推定方法では、往路のＣＳＩと復路のＣＳＩを結合し、生体を経由しない直接波を抽出し、周波数誤差によって生じる位相回転を検出し、チャネル全体から周波数誤差由来の位相回転を除去することで、生体の存在する方向を推定することが可能となる。
　以下では、図面を参照しながら、本実施の形態における推定装置１０が、所定の２つ無線機により得る往路および復路のＣＳＩを用いて、検出対象である生体の方向を推定することについて説明する。
　［推定装置１０の構成］
　図１は、本実施の形態における推定装置１０の構成の一例を示すブロック図である。図１には、図１に示す推定装置１０の検出対象である生体３０が合わせて示されている。
　図１に示す推定装置１０は、第二無線機１１と、第一無線機１２と、アンテナ部２０とアンテナ部２１と、第一送受信部２２と、第二送受信部２３と、第一結合複素伝達関数算出部２４と、第二結合複素伝達関数算出部２５と、生体情報算出部２６と、位置推定処理部２７とを備える。推定装置１０は、推定装置１０に対する相対的な、生体３０の存在する方向を推定する。
　なお、推定装置１０は、一の筐体内に上記の構成要素のすべてを備えるものに限定されない。推定装置１０は、互いに通信可能に接続された複数の装置を備え、複数の装置のそれぞれに、上記の構成要素が分かれて配置されていてもよい。その場合、推定装置１０は、推定システムと呼ばれてもよい。
　推定装置１０が、Ｍ個のアンテナ素子を持つ第一無線機１２、および、Ｎ個のアンテナ素子を持つ第二無線機１１を用いて、第一無線機１２から第二無線機１１に既知信号を送信した時のＣＳＩ（往路ＣＳＩに相当）を取得する場合を図１に示す。また、推定装置１０が、第二無線機１１から第一無線機１２に既知信号を送信した時のＣＳＩ（復路ＣＳＩに相当）を取得する場合を図３に示す。推定装置１０は、フルＭＩＭＯ構成を有しているが、一例として、ＳＩＭＯ構成のＣＳＩが往路と復路との分だけ得られた場合に生体３０の位置を推定する方法を説明する。
　［第一無線機１２］
　第一無線機１２は、生体３０の方向を推定するために用いる高周波の信号を、第一無線機１２の内部クロックｆＴＸを用いて生成する。例えば、図２に示すように、第一無線機１２は、生成した信号（送信波）を、アンテナ部２０が備えるアンテナ素子から送信する。また、第一無線機１２は、受信時は、上記アンテナ素子のそれぞれにより受信された受信信号を処理する。
　図４に本実施の形態における、推定装置１０が受信したＩ信号およびＱ信号の周波数成分の一例を示す概念図を示す。より具体的には、図４に示される点線１００は、推定装置１０が受信した信号のＤＣ成分を示し、実線１０１は、推定装置１０が受信したＩ信号のＡＣ成分を示し、破線１０２は、推定装置１０が受信したＱ信号のＡＣ成分を示す。
　［アンテナ部２０］
　アンテナ部２０は、Ｍ個（Ｍは２以上の自然数）のアンテナ素子（具体的には、アンテナ素子＃１～＃Ｍ）を有する。
　［第一送受信部２２］
　第一送受信部２２は、Ｍ個のアンテナ素子により、所定の信号を送信、または、受信する。
　第一送受信部２２は、送信の際は、生体推定用の既知の信号を送信する。
　第一送受信部２２は、受信の際は、アンテナ素子のそれぞれにより受信された受信信号であって、送信側のアンテナ素子から送信され、生体３０によって反射された反射信号を含む受信信号を、生体３０の活動に由来する周期に相当する第一期間に観測する。
　ここで、送信側のアンテナ素子は、第二無線機１１が備えるアンテナ部２１が有するアンテナ素子であり得る。生体３０は、図２に示すような生体３０であるがこれに限られない。また、生体３０の活動に由来する周期は、生体３０の呼吸、心拍および体動の少なくとも一つを含む生体由来の周期（生体変動周期）である。第一送受信部２２は、観測した受信信号を第二受信信号として第一結合複素伝達関数算出部２４に送る。
　［第二無線機１１］
　第二無線機１１は、生体３０の方向を推定するために用いる高周波の信号を、第一無線機１２と同様に内部クロックｆＲＸを用いて生成する。例えば、図２に示すように、第二無線機１１は、生成した信号（送信波）を、アンテナ部２１が備えるアンテナ素子から送信する。
　また、第二無線機１１は、受信時は、上記アンテナ素子のそれぞれにより受信された受信信号を処理する。
　［アンテナ部２１］
　アンテナ部２１は、Ｎ個（Ｎは２以上の自然数）のアンテナ素子（具体的には、アンテナ素子＃１～＃Ｎ）を有する。
　［第二送受信部２３］
　第二送受信部２３は、Ｎ個のアンテナ素子のそれぞれにより受信された受信信号であって、送信側のアンテナ素子から送信され、生体３０によって反射された反射信号を含む受信信号を、生体３０の活動に由来する周期に相当する第一期間について観測する。ここで、上記第一期間は、第一送受信部２２が受信信号を観測する第一期間と同じであってよい。
　第二送受信部２３は、観測した受信信号を第一受信信号として、第一結合複素伝達関数算出部２４に送る。
　本実施の形態では、第一送受信部２２と第二送受信部２３とは、それぞれ複数のアンテナ素子を持ち、対応するアンテナ素子で受信された高周波の信号を、内部クロックｆＴＸまたはｆＲＸを用いて、信号処理が可能な低周波の信号に変換する。
　より具体的には、ＳＩＭＯ構成のＣＳＩを利用する場合、第一送受信部２２は、アンテナ部２０が有する所定の１つのアンテナ素子よりＣＳＩ取得用の既知信号を送信する。送信した信号が生体３０によって反射された信号を第二無線機１１のアンテナ部２１が有するＮ個のアンテナ素子が受信する。受信した信号を、第二送受信部２３が第一受信信号として、第一結合複素伝達関数算出部２４に送る。
　さらに、第二送受信部２３は、アンテナ部２１が有する所定の１つのアンテナ素子よりＣＳＩ取得用の既知信号を送信する。送信した信号が生体３０によって反射された信号を第一無線機１２のアンテナ部２０が有するＭ個のアンテナ素子が受信する。受信した信号を、第一送受信部２２が第二受信信号として、第一結合複素伝達関数算出部２４に送る。
　［第一結合複素伝達関数算出部２４］
　第一結合複素伝達関数算出部２４は、第一受信信号と第二受信信号とから、第一無線機１２と第二無線機１１のそれぞれとの間の伝搬特性を表す複素伝達関数を複数算出する。
　第一結合複素伝達関数算出部２４は、第一送受信部２２および第二送受信部２３から伝達された低周波の第一受信信号、または、第二受信信号から、第一無線機１２と第二無線機１１との間の伝搬特性を表す複素伝達関数を算出する。
　以降において、往路ＣＳＩの取得について図２を用いて具体的に説明する。
　図２は、図１に示すアンテナ部２１における信号波の伝達の様子を概念的に示す図である。図２に示すように、アンテナ部２０のアンテナ素子から送信される送信波は、生体３０によって反射され、アンテナ部２１のアレーアンテナに到達する。ここで、アレーアンテナは、Ｍ個のアンテナ素子からなり、素子間隔ｄのリニアアレーである。また、アレーアンテナの正面から見た生体３０の方向をθとする。生体３０とアレーアンテナとの距離は十分に大きく、アレーアンテナに到来する生体３０由来の反射波は平面波と見なせるものとする。
　この場合、第一結合複素伝達関数算出部２４は、アレーアンテナを使って観測された第一受信信号ベクトル
Figure JPOXMLDOC01-appb-M000001
 から、送信側のアンテナ素子と受信側のアレーアンテナとの間の伝搬特性を表す第一複素伝達関数ベクトルを算出することができる。第一複素伝達関数ベクトルｈFは、例えば、
Figure JPOXMLDOC01-appb-M000002
Figure JPOXMLDOC01-appb-M000003
 により算出できる。ここで、sは複素送信信号であり、既知であるものとする。なお、往路ＣＳＩである第一複素伝達関数ベクトルｈFは、第一複素伝達関数の一例である。
　次に、復路ＣＳＩの取得について説明する。
　第一結合複素伝達関数算出部２４は、往路ＣＳＩの場合と同様に、アレーアンテナを使って観測された第二受信信号ベクトル
Figure JPOXMLDOC01-appb-M000004
 から、送信側のアンテナ素子と受信側のアレーアンテナとの間の伝搬特性を表す第二複素伝達関数ベクトルを算出することができる。第二複素伝達関数ベクトルｈBは、例えば、
Figure JPOXMLDOC01-appb-M000005
Figure JPOXMLDOC01-appb-M000006
 により算出できる。なお、
Figure JPOXMLDOC01-appb-M000007
 は転置を表す。ここで、本来（言い換えれば設計上）、ｈＦ１１とｈＢ１１とは同じ値であるが、クロック周波数誤差、電力誤差および取得時間差などにより位相が異なるので、下記により一方の位相に合わせてからチャネルを結合（つまり、第一複素伝達関数ベクトルと第二複素伝達関数ベクトルとを結合）する。
Figure JPOXMLDOC01-appb-M000008
　言い換えれば、第一結合複素伝達関数算出部２４は、第一送信信号の第一エレメントと第二送信信号の第一エレメントとを用いて、複数の第一受信信号から算出される第一複素伝達関数と、複数の第二受信信号から算出される第二複素伝達関数との一方の位相を他方の位相に合わせてから、第一複素伝達関数と第二複素伝達関数とを結合することで、第一結合複素伝達関数を算出する。
　これにより、第一結合複素伝達関数ベクトルｈ（ｔ）が下記のように表現される。
Figure JPOXMLDOC01-appb-M000009
　［第二結合複素伝達関数算出部２５］
　更に、第一結合複素伝達関数ベクトルｈ（ｔ）には、送信機および受信機由来の周波数変動成分と、生体３０由来のドップラシフトとが含まれる。送信機および受信機由来の周波数変動成分には、例えば、（ｉ）送信信号の空間伝播による減衰または位相回転、（ｉｉ）送信機および受信機間のクロック周波数誤差（ｆＲＸ－ｆＴＸ）、（ｉｉｉ）ＤＡ変換など無線機内にて用いられるサンプリングクロック周波数誤差などが含まれる。第一結合複素伝達関数ベクトルｈ（ｔ）から送信機および受信機由来の周波数変動成分の位相回転を除去するために、第二結合複素伝達関数算出部２５は、直接波成分としてｈ（ｔ）内の任意の一つの要素ｈｌを抽出する。
Figure JPOXMLDOC01-appb-M000010
　そして、第二結合複素伝達関数算出部２５は、第一結合複素伝達関数ベクトルｈ（ｔ）の全要素を直接波成分として、上記のように抽出した一つの要素ｈｌにて除算することで上記位相回転を除去することによって、第二結合複素伝達関数ベクトルｈ′を算出する。ここで、直接波成分の要素は、要素ｈ１など、第一結合複素伝達関数ベクトルｈ（ｔ）のうちの一つの要素であればどの要素を用いてもよい。なお、第二結合複素伝達関数ベクトルｈ′は、第二結合複素伝達関数の一例である。
　このように、第二結合複素伝達関数算出部２５は、第一結合複素伝達関数ベクトルｈ（ｔ）の１以上の要素を用いて所定の演算を行うことにより、第一結合複素伝達関数ベクトルｈ（ｔ）から、（１）第一無線機１２と第二無線機１１間とのクロック揺らぎ、および、（２）送信信号のデジタル－アナログ変換または受信信号のアナログ－デジタル変換のタイミング揺らぎ、の少なくとも一方に対応する成分が抑制された第二結合複素伝達関数ベクトルｈ′を算出する。
　［生体情報算出部２６］
　生体情報算出部２６は、算出された第二結合複素伝達関数ベクトルｈ′を用い、特許文献１に記載の差分情報算出部と同様に所定間隔の２つの時点における２つの第二結合複素伝達関数ベクトルｈ′の差分情報を算出することで、生体情報を算出する。つまり、生体情報算出部２６は、算出された複数の第二結合複素伝達関数ベクトルｈ′における所定周波数範囲に対応する生体情報を抽出することで、生体３０に関する成分に対応する生体情報を抽出する。例えば、生体情報算出部２６は、生体情報として、生体３０の呼吸、心拍および体動の少なくともいずれかを含むバイタル活動の影響を受けた成分に対応する生体情報を抽出する。
　なお、この生体情報は第二結合複素伝達関数ベクトルｈ′より生体由来の成分を算出すればよく、周波数応答を用いても同様の効果が得られることは言うまでもない。
　本実施の形態では、アンテナ素子が複数（Ｎ個）存在するので、アンテナ部２１に対応する第二複素伝達関数ベクトルｈ′の差分値（差分情報）の個数も複数である。これらをまとめて複素差分チャネルベクトルと定義する。アンテナ素子の個数をＮとすると複素差分チャネルベクトルは、
Figure JPOXMLDOC01-appb-M000011
 と表せ、
Figure JPOXMLDOC01-appb-M000012
 である。また、ｌまたはｍは、測定番号を表す正の整数であり、サンプル時間である。
　［位置推定処理部２７］
　位置推定処理部２７は、第二結合複素伝達関数ベクトルｈ′より抽出した当該生体情報を用いて、所定の到来方向推定手法により、推定装置１０を方向の基準として生体３０の存在する方向を推定する。より具体的には、所定の到来方向推定手法は、特許文献１記載の方向推定処理部のようなＭＵＳＩＣ（ＭＵｌｔｉｐｌｅＳＩｇｎａｌ　Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）アルゴリズムに基づく推定手法であってよく、この場合を例として説明するが、ビームフォーマ法またはＣａｐｏｎ法に基づく推定手法であってもよい。
　位置推定処理部２７は、例えば（式１）に示した複素差分チャネルベクトルを用いて、ＭＵＳＩＣアルゴリズムを適用する場合は、
Figure JPOXMLDOC01-appb-M000013
 により相関行列化し固有値分解する。ここでＵは固有ベクトル行列を表す。Ｅは、ｌおよびｍそれぞれを変化させた場合の平均を得る演算を意味する。更に、ＭＵＳＩＣスペクトラムＰについて下記（式２）が適用される。
Figure JPOXMLDOC01-appb-M000014
　ここでＵＮは雑音空間固有ベクトル、ａはステアリングベクトルである。
Figure JPOXMLDOC01-appb-M000015
　ａＦは往路のステアリングベクトル、ａＢは復路のステアリングベクトルである。
Figure JPOXMLDOC01-appb-M000016
　なお、第一無線機１２より送信される信号は、連続信号（ＣＷ信号）であってもよいし、ＯＦＤＭ（Ｏｒｔｈｏｇｏｎａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ）信号のような符号化された信号でもよい。
　［直接波成分の他の例］
　なお、本実施の形態の第二結合複素伝達関数算出部２５は、第一結合複素伝達関数ベクトルｈ（ｔ）を、直接波成分として第一結合複素伝達関数ベクトルｈ（ｔ）内の任意の一つの要素ｈｌにて除算することで、第二結合複素伝達関数ベクトルｈ′を算出したが、一つの要素ｈｌにて除算することに限らない。
　例えば、第二結合複素伝達関数算出部２５は、第一結合複素伝達関数ベクトルｈ（ｔ）の除算に用いる直接波成分には、従来の非同期吸収法にあるように任意の一つの要素ｈｌの代わりに、｜ｈｌ｜の時変動が最も小さい要素ｈｌｍｉｎを用いてもよい。また、第一複素伝達関数ベクトルｈ０の全要素の平均値を用いてもよい。また、複素伝達関数を一定期間観測することで、任意の時間ｔの複素伝達関数を特異値分解することで特異ベクトルを求め、これにより、第二結合複素伝達関数ベクトルｈ′を算出してもよい。また、複素伝達関数を一定期間観測し、当該観測時間全体の相関行列を固有値分解することで固有ベクトルを求め、これにより、第二結合複素伝達関数ベクトルｈ′を算出してもよい。また、当該直接波成分は、上記観測時間全体の相関行列を固有値分解することで固有ベクトルを求め、これにより第二結合複素伝達関数ベクトルｈ′を算出してもよい。
　なお、第一結合複素伝達関数算出部２４は、第一複素伝達関数ベクトルｈFと第二複素伝達関数ベクトルｈBを用いて、往路と復路で位相回転が逆に見えることを利用して
Figure JPOXMLDOC01-appb-M000017
Figure JPOXMLDOC01-appb-M000018
 を算出し、これにより
Figure JPOXMLDOC01-appb-M000019
 として第二複素伝達関数ベクトルを算出してもよい。これにより位相回転が消失し、２局の中間の周波数で観測したことと等価になる。
　なお、第一無線機１２、第二無線機１１、第一結合複素伝達関数算出部２４、第二結合複素伝達関数算出部２５、生体情報算出部２６および位置推定処理部２７は、推定装置１０が備えるメモリに格納されているプログラムを、１以上のプロセッサが実行することにより実現されてもよいし、１以上の専用回路により実現されてもよい。つまり、上記の構成要素は、ソフトウェアで実現されてもよいし、ハードウェアにより実現されてもよい。
　［推定装置１０の動作］
　以上のように構成された推定装置１０の推定処理の動作について説明する。図５は、本実施の形態における推定装置１０の推定処理を示すフローチャートである。
　まず、推定装置１０は、測定対象の領域に第一無線機１２、および、第二無線機１１を用いてそれぞれ送信信号を送信し、第一期間に受信信号を観測する（Ｓ１０）。より具体的には、推定装置１０は、１個のアンテナ素子から送信され、生体３０によって反射された反射信号を含む受信信号を、生体３０の活動に由来する周期に相当する第一期間について観測する。観測は、往路分と復路分とのそれぞれについてなされる。
　次に、推定装置１０は、ステップＳ１０において第一期間に観測された複数の受信信号から、１個の送信側アンテナ素子と複数の受信側アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出する（Ｓ２０）。詳細は上述した通りであるため、ここでの説明は省略する。以下も同様である。
　次に、推定装置１０は、ステップＳ２０で算出した第一結合複素伝達関数を直接波成分にて除算することで、第二結合複素伝達関数を算出する（Ｓ３０）。
　次に、推定装置１０は、ステップＳ３０で算出した第二結合複素伝達関数より生体情報を算出する（Ｓ４０）。
　そして、推定装置１０は、ステップＳ４０で算出した２以上の生体情報を用いて、生体３０の存在する位置を推定する（Ｓ５０）。
　［効果等］
　本実施の形態の推定装置１０および推定方法によれば、受信信号から生体３０を経由しない直接波成分を抽出し、周波数誤差によって生じる位相回転を検出し、チャネル全体から送信部および受信部間のクロック揺らぎ、および、周波数誤差由来の位相回転を除去するため、生体３０の存在する方向を精度よく推定することが可能となる。このため、推定装置１０は、推定装置１０に対する相対的な、生体３０が存在する方向を精度よく推定することができる。
　（実施の形態２）
　本実施の形態では、ＣＳＩの一部を少なくとも出力する無線機を用いて推定装置に対する生体の方向または位置を精度よく推定する推定装置および推定方法について、実施の形態１とは異なる形態を説明する。
　本実施の形態では、ＭＩＳＯ構成のＣＳＩであっても、往路のＣＳＩと復路のＣＳＩとを結合することを用いて、生体の存在する方向を推定する推定装置１０Ａおよび推定方法を説明する。
　［推定装置１０Ａの構成］
　図６および図７は、それぞれ、本実施の形態における推定装置１０Ａの構成の一例を示すブロック図である。図１と同様の要素には同一の符号を付しており、詳細な説明を省略する。本実施の形態において、実施の形態１と同じ要素については特に言及がない限り、同様の動作および同様のバリエーションがあるとし、重複した説明は省略する。
　図６に示す推定装置１０Ａは、第二無線機１１と、第一無線機１２と、アンテナ部２０とアンテナ部２１と、第一送受信部２２Ａと第二送受信部２３Ａと、第一結合複素伝達関数算出部２４Ａと、第二結合複素伝達関数算出部２５と、生体情報算出部２６と、位置推定処理部２７とを備える。推定装置１０Ａは、推定装置１０Ａに対する相対的な、生体３０の存在する方向を推定する。
　推定装置１０Ａが、Ｍ個のアンテナ素子を持つ第一無線機１２、および、Ｎ個のアンテナ素子を持つ第二無線機１１を用いて、第一無線機１２から第二無線機１１に既知信号を送信した時のＣＳＩ（往路ＣＳＩに相当）を取得する場合を図６に示す。また、推定装置１０Ａが、第二無線機１１から第一無線機１２に既知信号を送信した時のＣＳＩ（復路ＣＳＩに相当）を取得する場合を図７に示す。推定装置１０Ａは、フルＭＩＭＯ構成を有しているが、一例として、ＭＩＳＯ構成のＣＳＩが往路と復路との分だけ得られた場合に、推定装置１０Ａが生体３０の位置を推定する方法を説明する。
　［第一送受信部２２Ａ］
　第一送受信部２２Ａは、Ｍ個のアンテナ素子により、所定の信号を送信、または、受信する。
　第一送受信部２２Ａは、送信の際は、生体推定用の既知の信号を送信する。
　第一送受信部２２Ａは、受信の際は、アンテナ素子のそれぞれにより受信された受信信号であって、送信側のアンテナ素子から送信され、生体３０によって反射された反射信号を含む受信信号を、生体３０の活動に由来する周期に相当する第一期間について観測する。
　第一送受信部２２Ａは、観測した受信信号を第四受信信号として第一結合複素伝達関数算出部２４Ａに送る。
　［第二送受信部２３Ａ］
　第二送受信部２３Ａは、Ｎ個のアンテナ素子のそれぞれにより受信された受信信号であって、送信側のアンテナ素子から送信され、生体３０によって反射された反射信号を含む受信信号を、生体３０の活動に由来する周期に相当する第一期間について観測する。
　第二送受信部２３Ａは、観測した受信信号を第三受信信号として第一結合複素伝達関数算出部２４Ａに送る。
　より具体的には、ＭＩＳＯ構成のＣＳＩを利用する場合、第一送受信部２２Ａは、アンテナ部２０が有するＭ個のアンテナ素子がＣＳＩ取得用の既知信号を送信し、生体３０によって反射された信号を第二無線機１１のアンテナ部２１が有する所定の１個のアンテナ素子が受信する。受信した信号を、第二送受信部２３Ａが第三受信信号として、第一結合複素伝達関数算出部２４Ａに送る。
　さらに、第二送受信部２３Ａは、アンテナ部２１が有するＮ個のアンテナ素子よりＣＳＩ取得用の既知信号を送信する。送信した信号が生体３０によって反射された信号を第一無線機１２のアンテナ部２０が有する所定の１個のアンテナ素子が受信する。受信した信号を、第一送受信部２２Ａが第四受信信号として、第一結合複素伝達関数算出部２４Ａに送る。
　［第一結合複素伝達関数算出部２４Ａ］
　第一結合複素伝達関数算出部２４Ａは、第三受信信号と第四受信信号とから、第一無線機１２と第二無線機１１のそれぞれとの間の伝搬特性を表す複素伝達関数を複数算出する。
　第一結合複素伝達関数算出部２４Ａは、第一送受信部２２Ａおよび第二送受信部２３Ａから伝達された低周波の第三受信信号、または、第四受信信号から、第一無線機１２と第二無線機１１との間の伝搬特性を表す複素伝達関数を算出する。
　以降において、往路ＣＳＩの取得について説明する。
　この場合、第一結合複素伝達関数算出部２４Ａは、アレーアンテナを使って観測された第三受信信号ベクトル
Figure JPOXMLDOC01-appb-M000020
 から、送信側のアンテナ素子と受信側のアレーアンテナとの間の伝搬特性を表す第三複素伝達関数ベクトルを算出することができる。第三複素伝達関数ベクトルｈFは、例えば、
Figure JPOXMLDOC01-appb-M000021
Figure JPOXMLDOC01-appb-M000022
 により算出できる。
　次に、復路ＣＳＩの取得について説明する。
　第一結合複素伝達関数算出部２４Ａは、復路ＣＳＩの場合と同様に、アレーアンテナを使って観測された第四受信信号ベクトル
Figure JPOXMLDOC01-appb-M000023
 から、送信側のアンテナ素子と受信側のアレーアンテナとの間の伝搬特性を表す第四複素伝達関数ベクトルを算出することができる。第四複素伝達関数ベクトルｈBは、例えば、
Figure JPOXMLDOC01-appb-M000024
Figure JPOXMLDOC01-appb-M000025
 により算出できる。
　ここで、本来（言い換えれば設計上）、ｈＦ１１とｈＢ１１とは同じ値であるが、クロック周波数誤差、電力誤差および取得時間差などにより位相が異なるので、下記により一方の位相に合わせてからチャネルを結合（つまり、第一複素伝達関数ベクトルと第二複素伝達関数ベクトルとを結合）する。
　言い換えれば、第一結合複素伝達関数算出部２４Ａは、第三受信信号の第一エレメントと第四受信信号の第一エレメントとを用いて、複数の第三受信信号から算出される第一複素伝達関数と、複数の第四受信信号から算出される第二複素伝達関数との一方の位相を他方の位相に合わせてから、第一複素伝達関数と第二複素伝達関数とを結合することで、第一結合複素伝達関数を算出する。
Figure JPOXMLDOC01-appb-M000026
　これにより、ＳＩＭＯ構成の第一結合複素伝達関数ベクトルｈ（ｔ）が下記のように表現される。
Figure JPOXMLDOC01-appb-M000027
　更に第一結合複素伝達関数ベクトルｈ（ｔ）は、送信機および受信機由来の周波数変動成分と、生体３０由来のドップラシフトとが含まれる。送信機および受信機由来の周波数変動成分は、実施の形態１で説明した方法と同様の方法で除去すればよい。また、上記の除去以降、生体３０の位置を推定する方法は、実施の形態１で説明した方法と同様である。
　［効果等］
　本実施の形態の推定装置１０Ａおよび推定方法によれば、ＭＩＳＯ構成のＣＳＩであっても、往路のＣＳＩと復路のＣＳＩとを結合することを用いて、生体３０の存在する方向を推定することが可能となる。
　ここで、本実施の形態に係る効果を確かめるために実験による評価を行った。以下に、実験について説明する。
　（実験）
　図８は、本実施の形態に係る推定方法を用いたシミュレーション実験の概念を示す図である。図９は、本実施の形態に係る推定方法を用いた実験の条件を示す図である。
　図８に示す送信アレーアンテナ（Ｔｘ　ａｒｒａｙ）と受信アレーアンテナ（Ｒｘ　ａｒｒａｙ）との双方は、４素子パッチアレーアンテナを用いた４×４ＭＩＭＯ（Ｍｕｌｔｉｐｌｅ　Ｉｎｐｕｔ　Ｍｕｌｔｉｐｌｅ　Ｏｕｔｐｕｔ）構成である。また、送信側にＳＰ４Ｔ（Ｓｉｎｇｌｅ－Ｐｏｌｅ－４－Ｔｈｒｏｗ）スイッチを有し、受信側に４系統の受信機を有する。
　本実験では、これらの機器を用いてＭＩＭＯチャネルの測定を行った。
　ここで、図９に示すように、送受信アンテナのアレー素子間隔を０．５波長、送受信間距離Ｄを４．０ｍ、アンテナ高ｈを人間（Ｌｉｖｉｎｇ－Ｂｏｄｙ）の直立時の胸の高さである１．０ｍに設定した。送信機からは２．４７１２５ＧＨｚの無変調連続波（ＣＷ：Ｃｏｎｔｉｎｕｏｕｓ　Ｗａｖｅ）が送信され、チャネル測定時間は３３秒とした。チャネル測定時、被験者以外無人とし、被験者はアンテナ側の壁に対して正面を向いた状態とした。
　図１０は、実施の形態２に係る推定方法を用いたシミュレーション結果を示す図である。図１０の（ａ）は、フルＭＩＭＯのＣＳＩを用いて従来法によりシミュレーションした結果を示している。図１０の（ｂ）は、ＳＩＭＯ構成のＣＳＩを用いて従来法によりシミュレーションした結果を示している。
　図１０では、空間に被験者が２人存在している場合の生体位置推定の結果が示されている。実験時の被験者の立ち位置は、被験者１が（Ｘ＝０．４ｍ，Ｙ＝１．８ｍ）、被験者２が（Ｘ＝３．２ｍ，Ｙ＝３．０ｍ）であった。
　図１０の（ａ）および（ｂ）において、実際の被験者の位置が×印で示されている。また、空間内の各位置についてのＭＵＳＩＣスペクトラムＰ（（式２）参照）がグレースケールで示されており、黒に近い色を有する位置に被験者が存在すると推定されることを示している。
　図１０の（ｂ）に示されるＳＩＭＯ構成のＣＳＩの場合において、空間内の各位置において被験者が存在する確率は、図１０の（ａ）に示される場合とほぼ同じである。よって、図１０の（ｂ）よりＳＩＭＯ構成のＣＳＩの場合でも、実施の形態２に係る推定方法により人の生体位置推定が可能であることが分かる。
　図１１は、実施の形態２に係る推定方法を用いた別の実験結果を示す図である。図１１には、被験者が２人の場合の位置推定誤差の累積確率分布（ＣＤＦ：Ｃｕｍｕｌａｔｉｖｅ　Ｄｉｓｔｒｉｂｕｔｉｏｎ　Ｆｕｎｃｔｉｏｎ）が示されている。図１１において、横軸が位置推定誤差（単位：ｍ）を示し、位置推定誤差に対するＣＤＦが縦軸に示されている。図１１の（ａ）がフルＭＩＭＯのＣＳＩを用いた時のＣＤＦを示し、図１１の（ｂ）がＳＩＭＯ構成のＣＳＩを用いた時のＣＤＦを示している。
　図１１より、位置推定誤差０．１ｍのＣＤＦ値が０．９１と得られており、全体の９１％が誤差０．１ｍ以内に収まっていることを示している。したがって、実施の形態２に係る推定方法は、従来のフルＭＩＭＯのＣＳＩを用いた時の手法と同等に精度よく推定できていることが分かる。これにより、本実施の形態によってＳＩＭＯ構成のＣＳＩを用いた生体の位置推定であっても高い精度で生体位置を推定できることが示された。
　以上のように、本開示によれば、ＭＩＳＯ構成のＣＳＩであっても往路と復路のＣＳＩを取得することで、それらを結合し周波数誤差の処理を同様に行うことで生体の存在する方向を推定することが可能となる。
　以上、本開示の一態様に係る推定装置および推定方法について、実施の形態に基づいて説明したが、本開示は、これらの実施形態に限定されるものではない。本開示の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施形態に施したもの、あるいは異なる実施形態における構成要素を組み合わせて構築される形態も、本開示の範囲内に含まれる。
　例えば、実施の形態１および２では、生体３０の方向推定または位置推定を例として説明したが、推定処理の対象は生体３０に限らない。推定処理の対象は、高周波の信号が照射された場合に、その活動または動きによって反射波にドップラ効果を与える種々の生体（機械等）に適用可能である。
　また、本開示は、このような特徴的な構成要素を備える、推定装置として実現することができるだけでなく、推定装置に含まれる特徴的な構成要素をステップとする推定方法などとして実現することもできる。また、そのような方法に含まれる特徴的な各ステップをコンピュータに実行させるコンピュータプログラムとして実現することもできる。そして、そのようなコンピュータプログラムを、ＣＤ－ＲＯＭ等のコンピュータで読取可能な非一時的な記録媒体あるいはインターネット等の通信ネットワークを介して流通させることができるのは、言うまでもない。
　本開示は、無線信号を利用した生体の方向や位置を推定する推定装置に利用でき、特に、生体と機械を含む生体の方向や位置を測定する測定器、生体の方向や位置に応じた制御を行う家電機器、生体の侵入を検知する監視装置などに搭載される推定装置等に利用できる。
　１０、１０Ａ　推定装置
　１１　第二無線機
　１２　第一無線機
　２０、２１　アンテナ部
　２２、２２Ａ　第一送受信部
　２３、２３Ａ　第二送受信部
　２４、２４Ａ　第一結合複素伝達関数算出部
　２５　第二結合複素伝達関数算出部
　２６　生体情報算出部
　２７　位置推定処理部
　３０　生体

Claims (17)

  1. 　Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置が実行する推定方法であって、
    　測定対象の領域に、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子を用いて第一送信信号を複数回送信し、
    　前記Ｎ個の第二アンテナ素子のそれぞれにより受信された複数の第一受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第一受信信号を、前記生体の動きの周期に相当する第一期間に観測し、
    　測定対象の領域に、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子を用いて第二送信信号を複数回送信し、
    　前記Ｍ個の第一アンテナ素子のそれぞれにより受信された複数の第二受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第二受信信号を、前記第一期間に観測し、
    　前記第一期間に観測された前記複数の第一受信信号、および前記複数の第二受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出し、
    　所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出し、
    　複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を算出し、
    　所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する
    　推定方法。
  2. 　前記第一結合複素伝達関数を算出する際には、
    　前記第一送信信号の第一エレメントと前記第二送信信号の第一エレメントとを用いて、前記複数の第一受信信号から算出される第一複素伝達関数と、前記複数の第二受信信号から算出される第二複素伝達関数との一方の位相を他方の位相に合わせてから、前記第一複素伝達関数と前記第二複素伝達関数とを結合することで、前記第一結合複素伝達関数を算出する
    　請求項１記載の推定方法。
  3. 　前記第二結合複素伝達関数を算出する際には、
    　前記第一結合複素伝達関数の１以上の要素を用いて所定の演算を行うことにより、前記第一結合複素伝達関数から、
    　（１）アンテナ部から送信する送信信号を生成する送信部と受信部との間のクロック揺らぎ、および、
    　（２）前記送信信号のデジタル－アナログ変換、または、前記第一受信信号もしくは前記第二受信信号のアナログ－デジタル変換のタイミング揺らぎ、
    　の少なくとも一方に対応する成分が抑制された第二結合複素伝達関数を算出する
    　請求項２記載の推定方法。
  4. 　前記第二結合複素伝達関数を算出する際には、
    　前記第一送信信号の第一エレメントと前記第二送信信号の第一エレメントとによる位相補正を前記第一送信信号と前記第二送信信号とに施すことによって、前記第一送信信号と前記第二送信信号との一方に対する他方の位相回転を消失させることで、前記第二結合複素伝達関数を算出する
    　請求項１記載の推定方法。
  5. 　Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置が実行する推定方法であって、
    　測定対象の領域に、前記Ｍ個の第一アンテナ素子のそれぞれより第一送信信号を複数回送信し、
    　前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子により受信された複数の第三受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第三受信信号を、前記生体の動きの周期に相当する第一期間に観測し、
    　測定対象の領域に、前記Ｎ個の第二アンテナ素子のそれぞれより第二送信信号を複数回送信し、
    　前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子により受信された複数の第四受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第四受信信号を、前記第一期間に観測し、
    　前記第一期間に観測された前記複数の第三受信信号、および前記複数の第四受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出し、
    　所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出し、
    　複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を算出し、
    　所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する
    　推定方法。
  6. 　前記第一結合複素伝達関数を算出する際には、
    　前記第三受信信号の第一エレメントと前記第四受信信号の第一エレメントとを用いて、前記複数の第三受信信号から算出される第一複素伝達関数と、前記複数の第四受信信号から算出される第二複素伝達関数との一方の位相を他方の位相に合わせてから、前記第一複素伝達関数と前記第二複素伝達関数とを結合することで、前記第一結合複素伝達関数を算出する
    　請求項５記載の推定方法。
  7. 　前記第二結合複素伝達関数を算出する際には、
    　前記第一結合複素伝達関数の１以上の要素を用いて所定の演算を行うことにより、前記第一結合複素伝達関数から、
    　（１）アンテナ部から送信する送信信号を生成する送信部と受信部との間のクロック揺らぎ、および、
    　（２）前記送信信号のデジタル－アナログ変換、または、前記第三受信信号もしくは前記第四受信信号のアナログ－デジタル変換のタイミング揺らぎ、
    　の少なくとも一方に対応する成分が抑制された第二結合複素伝達関数を算出する
    　請求項６記載の推定方法。
  8. 　前記第二結合複素伝達関数を算出する際には、
    　前記第三受信信号の第一エレメントと前記第四受信信号の第一エレメントとによる位相補正を前記第三受信信号と前記第四受信信号とに施すことによって、前記第三受信信号と前記第四受信信号との一方に対する他方の位相回転を消失させることで、前記第二結合複素伝達関数を算出する
    　請求項５記載の推定方法。
  9. 　前記第二結合複素伝達関数の算出では、
    　前記第一結合複素伝達関数から一つの要素を抽出し、抽出した前記一つの要素で、前記第一結合複素伝達関数の全要素を除算することで、前記第二結合複素伝達関数を算出する
    　請求項３または７に記載の推定方法。
  10. 　前記第二結合複素伝達関数の算出では、
    　前記第一結合複素伝達関数の全要素の平均値を算出し、算出した前記平均値で、前記第一結合複素伝達関数の全要素を除算することで、前記第二結合複素伝達関数を算出する
    　請求項３または７に記載の推定方法。
  11. 　前記第二結合複素伝達関数の算出では、
    　前記第一結合複素伝達関数について所定の方法で相関行列を算出し、前記相関行列を固有値分解することで、１以上の固有値と、前記１以上の固有値それぞれに対応する固有ベクトルを算出し、
    　算出した１以上の前記固有ベクトルのうち、算出した前記１以上の固有値のうちの最大の固有値に対応する固有ベクトルを前記第一結合複素伝達関数に乗算することで、直接波のチャネル成分である第三複素伝達関数を算出し、
    前記第一結合複素伝達関数の全要素を前記第三複素伝達関数にて除算することで、前記第二結合複素伝達関数を算出する
    　請求項３または７に記載の推定方法。
  12. 　前記第二結合複素伝達関数の算出では、
    　前記第一結合複素伝達関数について所定の方法で特異値分解することで左特異ベクトルおよび右特異ベクトルを算出し、前記左特異ベクトルおよび前記右特異ベクトルを前記第一結合複素伝達関数に乗算することで、直接波のチャネル成分である第四複素伝達関数を算出し、
    　前記第一結合複素伝達関数の全要素を前記第四複素伝達関数にて除算することで、前記第二結合複素伝達関数を算出する
    　請求項３または７に記載の推定方法。
  13. 　前記所定の到来方向推定手法は、ＭＵＳＩＣ（ＭＵｌｔｉｐｌｅＳＩｇｎａｌ　Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法、ビームフォーマ法、およびＣａｐｏｎ法のいずれかを少なくとも用いる推定手法である
    　請求項１～８のいずれか１項に記載の推定方法。
  14. 　前記第一受信信号または前記第二受信信号は、ＯＦＤＭ（Ｏｒｔｈｏｇｏｎａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ）信号に変調されている信号である、
    　請求項１～４のいずれか１項に記載の推定方法。
  15. 　前記第三受信信号または前記第四受信信号は、ＯＦＤＭ（Ｏｒｔｈｏｇｏｎａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ）信号に変調されている信号である、
    　請求項５～８のいずれか１項に記載の推定方法。
  16. 　Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置であって、
    　前記第一無線機が、測定対象の領域に、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子を用いて第一送信信号を複数回送信し、
    　前記第二無線機が、前記Ｎ個の第二アンテナ素子のそれぞれにより受信された複数の第一受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第一受信信号を、前記生体の動きの周期に相当する第一期間に観測し、
    　前記第二無線機が、測定対象の領域に、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子を用いて第二送信信号を複数回送信し、
    　前記第一無線機が、前記Ｍ個の第一アンテナ素子のそれぞれにより受信された複数の第二受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第二受信信号を、前記第一期間に観測し、
    　前記推定装置は、さらに、
    　前記第一期間に観測された前記複数の第一受信信号、および前記複数の第二受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出する第一結合複素伝達関数算出部と、
    　所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出する第二結合複素伝達関数算出部と、
    　複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を抽出する生体情報算出部と、
    　所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する位置推定処理部とを備える
    　推定装置。
  17. 　Ｍ個（Ｍは２以上の自然数）の第一アンテナ素子を有する第一アンテナ部を備えた第一無線機と、Ｎ個（Ｎは２以上の自然数）の第二アンテナ素子を有する第二アンテナ部を備えた第二無線機とを備える推定装置であって、
    　前記第一無線機が、測定対象の領域に、前記Ｍ個の第一アンテナ素子のそれぞれより第一送信信号を複数回送信し、
    　前記第二無線機が、前記Ｎ個の第二アンテナ素子のうちの一つの第二アンテナ素子により受信された複数の第三受信信号であって、複数回送信された前記第一送信信号が生体によって反射された反射信号を含む複数の第三受信信号を、前記生体の動きの周期に相当する第一期間に観測し、
    　前記第二無線機が、測定対象の領域に、前記Ｎ個の第二アンテナ素子のそれぞれより第二送信信号を複数回送信し、
    　前記第一無線機が、前記Ｍ個の第一アンテナ素子のうちの一つの第一アンテナ素子により受信された複数の第四受信信号であって、複数回送信された前記第二送信信号が前記生体によって反射された反射信号を含む複数の第四受信信号を、前記第一期間に観測し、
    　前記推定装置は、さらに、
    　前記第一期間に観測された前記複数の第三受信信号、および前記複数の第四受信信号を用いて、前記Ｍ個の第一アンテナ素子と前記Ｎ個の第二アンテナ素子のそれぞれとの間の伝搬特性を表す第一結合複素伝達関数を複数算出する第一結合複素伝達関数算出部と、
    　所定の方法で、複数の前記第一結合複素伝達関数に、前記第一無線機と前記第二無線機とに生じたクロック周波数誤差と送信電力誤差との影響を抑制する演算を施すことで、複数の第二結合複素伝達関数を算出する第二結合複素伝達関数算出部と、
    　複数の前記第二結合複素伝達関数における、前記生体に関する成分に対応する所定周波数範囲に対応する生体情報を抽出する生体情報算出部と、
    　所定の到来方向推定手法により、前記生体情報を用いて、前記生体の存在する方向または位置を推定する位置推定処理部とを備える
    　推定装置。
PCT/JP2022/023087 2021-06-29 2022-06-08 推定方法、および、推定装置 WO2023276592A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202280044197.1A CN117546046A (zh) 2021-06-29 2022-06-08 推定方法及推定装置
EP22832749.0A EP4365621A1 (en) 2021-06-29 2022-06-08 Estimation method and estimation device
JP2023531747A JPWO2023276592A1 (ja) 2021-06-29 2022-06-08

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2021-108202 2021-06-29
JP2021108202 2021-06-29

Publications (1)

Publication Number Publication Date
WO2023276592A1 true WO2023276592A1 (ja) 2023-01-05

Family

ID=84690278

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2022/023087 WO2023276592A1 (ja) 2021-06-29 2022-06-08 推定方法、および、推定装置

Country Status (4)

Country Link
EP (1) EP4365621A1 (ja)
JP (1) JPWO2023276592A1 (ja)
CN (1) CN117546046A (ja)
WO (1) WO2023276592A1 (ja)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012112712A (ja) * 2010-11-22 2012-06-14 Mitsubishi Electric Corp 送受信装置
JP2014216786A (ja) 2013-04-24 2014-11-17 株式会社チノー センシングシステム
JP2015072173A (ja) 2013-10-02 2015-04-16 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． 位置推定装置及び位置推定方法
JP2017129558A (ja) 2016-01-15 2017-07-27 パナソニックＩｐマネジメント株式会社 推定装置および推定方法
JP2020109389A (ja) 2018-12-28 2020-07-16 パナソニックＩｐマネジメント株式会社 推定方法、推定装置およびプログラム
KR102243445B1 (ko) * 2020-09-24 2021-04-22 한화시스템 주식회사 안테나 장치, 이를 구비하는 능동 위상배열 레이더, 및 이의 오차 교정방법

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012112712A (ja) * 2010-11-22 2012-06-14 Mitsubishi Electric Corp 送受信装置
JP2014216786A (ja) 2013-04-24 2014-11-17 株式会社チノー センシングシステム
JP2015072173A (ja) 2013-10-02 2015-04-16 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． 位置推定装置及び位置推定方法
JP2017129558A (ja) 2016-01-15 2017-07-27 パナソニックＩｐマネジメント株式会社 推定装置および推定方法
JP2019132850A (ja) * 2016-01-15 2019-08-08 パナソニックＩｐマネジメント株式会社 プログラム
JP2020190566A (ja) * 2016-01-15 2020-11-26 パナソニックＩｐマネジメント株式会社 プログラム
JP2020109389A (ja) 2018-12-28 2020-07-16 パナソニックＩｐマネジメント株式会社 推定方法、推定装置およびプログラム
KR102243445B1 (ko) * 2020-09-24 2021-04-22 한화시스템 주식회사 안테나 장치, 이를 구비하는 능동 위상배열 레이더, 및 이의 오차 교정방법

Also Published As

Publication number Publication date
CN117546046A (zh) 2024-02-09
JPWO2023276592A1 (ja) 2023-01-05
EP4365621A1 (en) 2024-05-08

Similar Documents

Publication Publication Date Title
JP6865394B2 (ja) プログラム
US11150336B2 (en) Positioning sensor and direction estimation method
US10895631B2 (en) Sensor and estimating method
US10928496B2 (en) Sensor and method for estimating position of living body
US11163057B2 (en) Estimation device, living body count estimation device, estimation method, and recording medium
JP7349661B2 (ja) 推定方法、推定装置およびプログラム
US20220214421A1 (en) Estimation device, estimation method, and recording medium
JP2014228291A (ja) 無線検出装置及び無線検出方法
JP6832534B2 (ja) 推定装置および推定方法
JP7162192B2 (ja) 生体数推定装置、生体数推定方法、及び、プログラム
JP2018112539A (ja) センサおよび位置推定方法
WO2023276592A1 (ja) 推定方法、および、推定装置
US11047968B2 (en) Estimating method and estimating device
CN111381229B (zh) 推测方法、推测装置以及记录介质
US20220221549A1 (en) Sensor and position estimation method
EP4365622A1 (en) Estimation device and estimation method
US20220350014A1 (en) Sensor
EP4088656A1 (en) Estimation device and estimation method
Valls et al. Vital Signs Estimation Using a 26 GHz Multi-Beam Communication Testbed
Hu et al. Respiration estimation for indoor human via commodity WiFi based on time-frequency analysis

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22832749

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2023531747

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 18572093

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 202280044197.1

Country of ref document: CN

WWE Wipo information: entry into national phase

Ref document number: 2022832749

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2022832749

Country of ref document: EP

Effective date: 20240129