WO2020110393A1 - 温度計測装置、計測装置、音波受信装置、及び、プログラム - Google Patents

温度計測装置、計測装置、音波受信装置、及び、プログラム Download PDF

Info

Publication number
WO2020110393A1
WO2020110393A1 PCT/JP2019/033965 JP2019033965W WO2020110393A1 WO 2020110393 A1 WO2020110393 A1 WO 2020110393A1 JP 2019033965 W JP2019033965 W JP 2019033965W WO 2020110393 A1 WO2020110393 A1 WO 2020110393A1
Authority
WO
WIPO (PCT)
Prior art keywords
sound wave
waveform data
measuring device
temperature
temperature measuring
Prior art date
Application number
PCT/JP2019/033965
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 JP2020558099A priority Critical patent/JP6898022B2/ja
Priority to CN201980078514.XA priority patent/CN113167660A/zh
Publication of WO2020110393A1 publication Critical patent/WO2020110393A1/ja
Priority to US17/229,120 priority patent/US20210231507A1/en

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/22—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects
    • G01K11/24—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects of the velocity of propagation of sound
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P13/00—Indicating or recording presence, absence, or direction, of movement
    • G01P13/02—Indicating direction only, e.g. by weather vane
    • G01P13/04—Indicating positive or negative direction of a linear movement or clockwise or anti-clockwise direction of a rotational movement
    • G01P13/045—Indicating positive or negative direction of a linear movement or clockwise or anti-clockwise direction of a rotational movement with speed indication
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/245—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by measuring transit time of acoustical waves

Definitions

  • the present invention relates to a temperature measuring device, a measuring device, a sound wave receiving device, and a program.
  • Patent Document 1 discloses a technique of arranging a plurality of sensor units capable of transmitting and receiving ultrasonic waves in a space to measure the temperature of the space from the propagation time of the ultrasonic waves.
  • the sensor unit since the ultrasonic wave has no directivity, the sensor unit receives not only the desired ultrasonic wave component but also the ultrasonic wave component reflected by the reflecting member (for example, a wall) existing in the space. .. Therefore, ultrasonic waves may be erroneously detected.
  • the S/N ratio when measuring the temperature of the space is low. As a result, the accuracy of the measurement result of the temperature of the space decreases.
  • the object of the present invention is to improve the S/N ratio of the measurement result of the air characteristic (for example, temperature) of the space.
  • One aspect of the present invention is A sound wave transmitting device arranged in the target space, and a temperature measuring device connectable to the sound wave receiving device arranged in the target space, A means for controlling the sound wave transmitting device to transmit a sound wave, From a sound wave receiving device that receives the sound wave, a means for acquiring received waveform data relating to the waveform of the sound wave, A means for calculating a temperature distribution of the target space based on the received waveform data, It is a temperature measuring device.
  • FIG. 9 is a schematic diagram showing a configuration of a sound wave transmitting device and a sound wave receiving device according to a first modification.
  • FIG. 9 is a schematic diagram showing configurations of a sound wave transmitting device and a sound wave receiving device according to a modification 4; It is explanatory drawing of the principle of the example 4 of a change. It is an explanatory view of the outline of the modification 5.
  • 13 is a flowchart of a temperature measurement process of modification 5; FIG.
  • FIG. 11 is a schematic diagram of a sound wave transmitting device and a sound wave receiving device according to a modification 6; It is a figure which shows an example of the sensor arrangement of the modification 8. It is a figure which shows the structure of the sound wave transmission device of the modification 9. It is a figure which shows the structure of the sound wave receiving device of the modification 9. It is a figure which shows the outline
  • FIG. 1 is a block diagram showing the configuration of the measurement system of the first embodiment.
  • FIG. 2 is a block diagram showing a detailed configuration of the measurement system of the first embodiment.
  • the measurement system 1 includes a temperature measurement device 10, a sound wave transmission device 20, a sound wave reception device 30, an air conditioner 40, and a thermometer 50.
  • the temperature measuring device 10 is connected to the sound wave transmitting device 20, the sound wave receiving device 30, the air conditioner 40, and the thermometer 50.
  • the temperature measuring device 10, the sound wave transmitting device 20, the sound wave receiving device 30, the air conditioner 40, and the thermometer 50 are arranged in the target space SP.
  • the temperature measuring device 10 has the following functions. -Function of controlling the sound wave transmitting device 20-Function of acquiring received waveform data from the sound wave receiving device 30-Function of measuring temperature distribution of the target space SP-Function of controlling the air conditioner 40-Function of the thermometer 50 of the target space SP Function of Acquiring Reference Temperature Information Regarding Temperature Measurement Results
  • the temperature measuring device 10 is, for example, a smartphone, a tablet terminal, or a personal computer.
  • the sound wave transmission device 20 is configured to transmit an ultrasonic beam having directivity (an example of “sound wave”) under the control of the temperature measurement device 10. Further, the sound wave transmission device 20 is configured to change the transmission direction of ultrasonic waves.
  • the sound wave receiving device 30 is configured to receive the ultrasonic beam transmitted from the sound wave transmitting device 20 and generate received waveform data according to the received ultrasonic beam.
  • the sound wave receiving device 30 is, for example, an omnidirectional microphone or a directional microphone.
  • the air conditioner 40 is configured to adjust the temperature of the target space SP under the control of the temperature measuring device 10.
  • the thermometer 50 is configured to measure the temperature of the target space SP (hereinafter referred to as “reference temperature”).
  • the temperature measuring device 10 includes a storage device 11, a processor 12, an input/output interface 13, and a communication interface 14.
  • the storage device 11 is configured to store programs and data.
  • the storage device 11 is, for example, a combination of a ROM (Read Only Memory), a RAM (Random Access Memory), and a storage (for example, a flash memory or a hard disk).
  • the programs include, for example, the following programs.
  • -OS Operating System
  • program-Application program that executes information processing (for example, information processing for measuring the temperature distribution of the target space SP)-Data relating to sound wave velocity characteristics relating to the speed of sound waves with respect to space temperature
  • the data includes, for example, the following data.
  • -Database referred to in information processing-Data obtained by executing information processing that is, execution result of information processing
  • the processor 12 is configured to realize the function of the temperature measuring device 10 by activating a program stored in the storage device 11.
  • the processor 12 is an example of a computer.
  • the input/output interface 13 is configured to acquire a user instruction from an input device connected to the temperature measuring device 10 and output information to an output device connected to the temperature measuring device 10.
  • the input device is, for example, a keyboard, a pointing device, a touch panel, or a combination thereof.
  • the input device also includes a thermometer 50.
  • the output device is, for example, a display.
  • the output device also includes the air conditioner 40.
  • the communication interface 14 is configured to control communication with an external device (for example, a server).
  • an external device for example, a server.
  • FIG. 3 is a schematic diagram showing the configuration of the sound wave transmitting device of the first embodiment.
  • the sound wave transmission device 20 includes a plurality of ultrasonic vibrators (an example of “vibration element”) 21 and a control circuit 22.
  • the control circuit 22 vibrates the plurality of ultrasonic transducers 21 under the control of the temperature measuring device 10.
  • the plurality of ultrasonic transducers 21 vibrate, ultrasonic beams are transmitted in the transmission direction (Z-axis direction) orthogonal to the transmission surface (XY plane).
  • FIG. 4 is a schematic diagram showing the configuration of the sound wave receiving apparatus of the first embodiment.
  • the sound wave receiving device 30 includes an ultrasonic vibrator 31 and a control circuit 32.
  • the ultrasonic transducer 31 vibrates when receiving the ultrasonic beam transmitted from the sound wave transmitting device 20.
  • the control circuit 32 is configured to generate received waveform data according to the vibration of the ultrasonic transducer 31.
  • FIG. 5 is an explanatory diagram of the outline of the first embodiment.
  • a temperature measurement device 10 (not shown), sound wave transmission devices 20a to 20b, and sound wave reception devices 30a to 30b are provided in a space (hereinafter referred to as “target space”) SP that is a target of temperature measurement. And are placed.
  • the temperature measuring device 10 can be connected to the sound wave transmitting device 20 and the sound wave receiving device 30.
  • the temperature measuring device 10 controls the sound wave transmitting devices 20a to 20b so as to transmit a sound wave.
  • the temperature measuring device 10 acquires received waveform data regarding the waveform of the received sound wave from the sound wave receiving devices 30a and 30b.
  • the temperature measuring device 10 calculates the temperature distribution of the target space SP based on the received waveform data.
  • the temperature distribution of the target space SP is calculated from the propagation time of sound waves. Therefore, the S/N ratio of the measurement result of the temperature of the space can be improved.
  • FIG. 6 is a diagram showing the data structure of the spatial data table of the first embodiment.
  • the spatial data table of FIG. 6 stores spatial information regarding the space in which the sound wave transmitting device 20 and the sound wave receiving device 30 are arranged (hereinafter referred to as “target space”).
  • the spatial data table includes a "coordinate” field and a "reflection characteristic” field. Each field is associated with each other.
  • the "coordinates” field stores the coordinates of the reflecting member existing in the target space (hereinafter referred to as “reflecting member coordinates”).
  • the reflection member coordinates are represented by a coordinate system having an arbitrary reference point in the target space as an origin (hereinafter referred to as “spatial coordinate system”).
  • “Reflection characteristic” field stores reflection characteristic information regarding the reflection characteristic of the reflection member.
  • the "reflection characteristic” field includes a “reflection type” field, a “reflectance” field, and a “normal angle” field.
  • the reflection type is one of the following. ⁇ Diffuse reflection ⁇ Specular reflection
  • the value of the reflectance of the reflecting member is stored in the "Reflectance" field.
  • the value of the normal angle of the reflection surface of the reflection member is stored in the "normal angle" field.
  • FIG. 7 is a diagram showing a data structure of the sensor data table of the first embodiment.
  • the sensor data table stores information about the sound wave transmitting device 20 and the sound wave receiving device 30 (hereinafter referred to as “sensor information”).
  • the sensor data table includes a “sensor ID” field, a “coordinate” field, and a “sensor type” field. Each field is associated with each other.
  • “Sensor ID” field stores sensor identification information for identifying the sound wave transmitting device 20 or the sound wave receiving device 30.
  • the “coordinates” field stores coordinates indicating the position of the sound wave transmitting device 20 or the sound wave receiving device 30 (hereinafter referred to as “sensor coordinates”).
  • the sensor coordinates are represented by a spatial coordinate system.
  • a tag “transmission” indicating the sound wave transmitting device 20 or a tag “reception” indicating the sound wave receiving device 30 is stored.
  • FIG. 8 is a diagram showing a data structure of the route data table of the first embodiment.
  • the route data table includes a "route ID” field, a "transmission sensor” field, and a "reception sensor” field.
  • the route identification information for identifying the route is stored in the "route ID" field.
  • the "transmission sensor” field stores sensor identification information of the sound wave transmitting device 20 that constitutes the route.
  • the "reception sensor” field stores the sensor identification information of the sound wave receiving device 30 that constitutes the route.
  • FIG. 9 is a diagram showing the data structure of the mesh data table of the first embodiment.
  • FIG. 10 is an explanatory diagram of the filter according to the first embodiment.
  • the mesh data table includes a “mesh ID” field, a “coordinate” field, a “route ID” field, and a “filter” field.
  • the mesh identification information for identifying the virtual mesh is stored in the “mesh ID” field.
  • the "coordinates" field stores mesh coordinates indicating the position of the virtual mesh.
  • the mesh coordinates are represented by a spatial coordinate system.
  • the route identification information of the route is stored in the "route ID" field.
  • the target space SP is divided into a plurality of virtual meshes Mi (i is an argument).
  • Each virtual mesh Mi has a three-dimensional shape.
  • the virtual mesh M1 includes a plurality of paths P101 and P200.
  • the path P101 is a path from the sound wave transmitting device 20a to the sound wave receiving device 30a.
  • the path P200 is a path from the sound wave transmitting device 20b to the sound wave receiving device 30b.
  • the “filter” field stores filter information regarding a filter for extracting a specific waveform from the waveform of the ultrasonic beam reproduced by the received waveform data received by the sound wave receiving device 30.
  • the filter information is associated with the route identification information stored in the “route ID” field.
  • the “filter” field includes a "time filter” field and an "amplitude filter” field.
  • the “time filter” field stores information about a time filter for extracting a specific waveform along the time axis.
  • the temporal filter is, for example, at least one of the following (FIG. 10).
  • ⁇ Lower limit time threshold THtb ⁇ Upper limit time threshold THtt
  • the “amplitude filter” field stores information about the amplitude filter for extracting the specific waveform along the amplitude axis.
  • the amplitude filter is, for example, at least one of the following (FIG. 10).
  • ⁇ Lower limit amplitude threshold THab ⁇ Upper limit amplitude threshold THat
  • FIG. 11 is a diagram showing an example of the sensor arrangement of the first embodiment.
  • FIG. 12 is a flowchart of the temperature measurement process of the first embodiment.
  • FIG. 13 is a detailed flowchart of calculating the path temperature of FIG.
  • FIG. 14 is a diagram showing an example of a screen displayed in the process of FIG.
  • a plurality of sound wave transmitting devices 20a to 20e and a plurality of sound wave receiving devices 30a to 30e are arranged in the target space SP.
  • the plurality of sound wave transmitters 20a to 20e face the plurality of sound wave receivers 30a to 30e, respectively.
  • the sound wave transmitting device 20a faces the sound wave receiving device 30a.
  • FIG. 11 shows an example in which five sensor pairs are formed.
  • FIG. 11 shows an example in which the mesh temperatures of the virtual meshes M1 to M4 can be measured.
  • step S110 the temperature measuring device 10 calculates the route temperature (S111) according to a predetermined route temperature calculation model. Details of step S111 will be described with reference to FIG.
  • the temperature measuring device 10 executes the determination of the target route (S1110). Specifically, the processor 12 refers to the mesh data table (FIG. 9) and stores the information of the “route ID” field associated with the mesh identification information determined in step S110 (that is, the route passing through the target mesh Mt ( Hereinafter, referred to as a “target route”) Pi (i is a route argument) is specified.
  • the processor 12 refers to the mesh data table (FIG. 9) and stores the information of the “route ID” field associated with the mesh identification information determined in step S110 (that is, the route passing through the target mesh Mt ( Hereinafter, referred to as a “target route”) Pi (i is a route argument) is specified.
  • the temperature measuring device 10 outputs an ultrasonic beam (S1111).
  • the processor 12 refers to the route data table (FIG. 8 ), and the information in the “transmission sensor” field associated with the route identification information identified in step S1110 (that is, the sound wave transmission device to be controlled). (Hereinafter, referred to as “target sound wave transmitting device”) 20) and information in the “reception sensor” field (that is, the sound wave receiving device to be controlled (hereinafter referred to as “target sound wave receiving device”) 30) are specified.
  • the processor 12 transmits an ultrasonic wave control signal to the target sound wave transmission device 20.
  • the target sound wave transmission device 20 transmits an ultrasonic beam in accordance with the ultrasonic wave control signal transmitted from the temperature measurement device 10. Specifically, the plurality of ultrasonic transducers 21 simultaneously vibrate according to the ultrasonic control signal. Thereby, the ultrasonic beam that travels in the transmission direction (Z-axis direction) is transmitted from the target sound wave transmitting device 20 toward the target sound wave receiving device 30.
  • the temperature measuring device 10 acquires the reception waveform data (S1112). Specifically, the ultrasonic transducer 31 of the target sound wave receiving device 30 vibrates by receiving the ultrasonic beam transmitted from the target sound wave transmitting device 20 in step S1111.
  • the control circuit 32 generates reception waveform data (FIG. 10) according to the vibration of the ultrasonic transducer 31.
  • the control circuit 32 transmits the generated reception waveform data to the temperature measuring device 10.
  • the processor 12 of the temperature measuring device 10 acquires the received waveform data transmitted from the sound wave receiving device 30.
  • the temperature measuring device 10 executes filtering (S1113).
  • the processor 12 refers to the mesh data table (FIG. 9) and specifies the “filter” field associated with the route identification information of the target route Pi determined in step S110. For example, when the mesh identification information of the target mesh is “M001”, the following filter information is specified.
  • the processor 12 extracts the component of the ultrasonic beam that has traveled along the target path Pi from the components included in the received waveform data based on the specified filter information.
  • the temperature measurement device 10 executes the calculation of the path temperature (S1114).
  • the processor 12 refers to the “coordinates” field of the sensor data table (FIG. 7) and specifies the coordinates of the sound wave transmitting device 20 and the coordinates of the sound wave receiving device 30 that configure the sensor pair for each sensor pair. ..
  • the processor 12 determines the distance between the sound wave transmitting device 20 and the sound wave receiving device 30 (hereinafter referred to as “inter-sensor distance”) based on the identified combination of the sound wave transmitting device 20 and the sound wave receiving device 30. Calculate Ds.
  • the processor 12 identifies the time (hereinafter referred to as “propagation time”) t corresponding to the peak value of the component extracted in step S1113.
  • the propagation time t is the time required for the ultrasonic beam that has traveled along the target path Pi to reach the sound wave receiving apparatus 30 after the ultrasonic wave transmitting device 20 transmits the ultrasonic beam (that is, from the starting point of the target path). It means the time during which the ultrasonic beam propagates to the end point).
  • the processor 12 calculates the path temperature TEMPpathi of the target path Pi using the sound velocity C of the ultrasonic waves, the inter-sensor distance Ds, the propagation time t, and the reference temperature T0.
  • step S1114 If step S1114 has not been completed for all target routes Pi (S1115-NO), the temperature measuring device 10 executes step S1110.
  • step S1114 When step S1114 is completed for all target paths Pi (S1115-YES), the temperature measuring device 10 executes the mesh temperature calculation (S112) of FIG. Specifically, the processor 12 calculates the mesh temperature TEMPmesht of the target mesh Mt using the path temperatures TEMPpathi of all the target paths Pi calculated in step S1114 (FIG. 11) (Equation 1).
  • TEMPmesht AVE(TEMPpathi)... (Formula 1)
  • AVE(x) a function for obtaining the average value of x
  • step S112 If step S112 has not been completed for all target meshes Mt (S113-NO), the temperature measuring device 10 executes step S110.
  • step S112 When step S112 is completed for all target meshes Mt (S113-YES), the temperature measuring device 10 presents the measurement result (S114). Specifically, the processor 12 displays the screen P10 (FIG. 14) on the display.
  • the screen P10 includes a display object A10.
  • the image IMG10 is displayed on the display object A10.
  • the image IMG10 shows the mesh temperature TEMPmesht calculated in step S112 for each of the plurality of virtual meshes forming the target space SP.
  • the temperature distribution of the target space SP is calculated from the propagation time of the ultrasonic beam. Therefore, the S/N ratio of the measurement result of the temperature of the space can be improved.
  • the second embodiment is an example in which the transmission direction of the ultrasonic beam of the acoustic wave transmission device 20 is variable.
  • FIG. 15 is a schematic diagram showing the configuration of the sound wave transmitting device of the second embodiment.
  • the sound wave transmission device 20 includes a plurality of ultrasonic vibrators 21, a control circuit 22, and an actuator 23.
  • the plurality of ultrasonic transducers 21 are two-dimensionally arranged on the transmission surface (XY plane). That is, the ultrasonic transducers 21 form a transducer array TA.
  • the actuator 23 is configured to change the orientation of the transmission surface (XY plane) with respect to the transmission axis (Z axis).
  • the actuator 23 directs the transmission surface toward the transmission axis (Z axis)
  • the ultrasonic beam USW0 is transmitted.
  • the actuator 23 tilts the transmission surface with respect to the transmission axis (Z axis)
  • the ultrasonic beam USW1 is transmitted.
  • FIG. 16 is a diagram showing a data structure of the route data table of the second embodiment.
  • FIG. 17 is an explanatory diagram of the transmission angle according to the second embodiment.
  • the route data table includes a “transmission angle” field in addition to the fields of FIG. 8 (“route ID” field, “transmission sensor” field, and “reception sensor” field).
  • the value of the transmission angle of the ultrasonic beam with respect to the transmission axis (Z axis) of the sound wave transmission device 20 is stored in the “transmission angle” field.
  • the paths P1 to P3 are specified by the combination of the sound wave transmitting device 20, the sound wave receiving device 30, and the transmission angle ⁇ .
  • FIG. 18 is a diagram showing an example of the sensor arrangement of the second embodiment.
  • FIG. 19 is a detailed flowchart of calculating the path temperature according to the second embodiment.
  • a plurality of sound wave transmitting devices 20a to 20b and a plurality of sound wave receiving devices 30a to 30b are arranged in the target space SP.
  • Each of the plurality of sound wave transmitting devices 20a to 20b is capable of transmitting an ultrasonic beam along a path reaching the plurality of sound wave receiving devices 30a to 30b.
  • the sound wave transmitting device 20a can transmit an ultrasonic beam along a path P20a reaching the sound wave receiving device 30a and an ultrasonic beam along a path P21a reaching the sound wave receiving device 30b.
  • the mesh temperature of the virtual mesh including the path of the ultrasonic beam. Since the paths P20b and P21b in which the ultrasonic beam is reflected by the reflection member (for example, at least one of the wall and the ceiling) of the target space SP are included, it is possible to measure the mesh temperature of more virtual meshes.
  • the temperature measurement device 10 of the second embodiment executes the determination of the transmission direction (S2110) after step S1110.
  • the processor 12 refers to the route data table (FIG. 8 ), and the information in the “transmission sensor” field associated with the route identification information identified in step S1110 (that is, the sound wave transmission device to be controlled). (Hereinafter referred to as "target sound wave transmitting device") 20), information in the "reception sensor” field (that is, the sound wave receiving device to be controlled (hereinafter referred to as “target sound wave receiving device”) 30), and "transmission angle” Field information (that is, the value of the transmission angle of the target sound wave transmission device 20) is specified.
  • the processor 12 transmits an ultrasonic wave control signal to the target sound wave transmission device 20.
  • the ultrasonic control signal includes the value of the transmission angle.
  • the target sound wave transmission device 20 transmits an ultrasonic beam in the direction indicated by the value of the transmission angle included in the ultrasonic control signal transmitted from the temperature measurement device 10.
  • the actuator 23 changes the orientation of the transmission surface (XY plane) with respect to the transmission axis (Z axis) with reference to the value of the transmission angle included in the ultrasonic control signal.
  • the control circuit 22 vibrates the plurality of ultrasonic transducers 21 at the same time. Thereby, the ultrasonic beam traveling in the direction indicated by the value of the transmission angle included in the ultrasonic control signal is transmitted.
  • step S2110 the temperature measuring device 10 executes steps S1111 to S1116 as in FIG.
  • the transmission angle of the sound wave transmitting device 20 is variable.
  • the paths of ultrasonic beams transmitted from one sound wave transmission device 20 increase.
  • Modification 1 is a modification regarding the sound wave transmitting device 20 and the sound wave receiving device 30.
  • FIG. 20 is a schematic diagram showing the configurations of the sound wave transmitting device and the sound wave receiving device of the first modification.
  • the sensor unit SU is arranged in the target space SP of the first modification.
  • the sensor unit SU includes a sound wave transmitting device 20 and a sound wave receiving device 30.
  • the sound wave transmission device 20 includes a plurality of ultrasonic transducers 21.
  • the plurality of ultrasonic transducers 21 are arranged in a fixed direction (X direction).
  • the sound wave receiving device 30 includes a plurality of ultrasonic transducers 31.
  • the plurality of ultrasonic transducers 31 are arranged in a fixed direction (Y direction).
  • the arrangement direction of the ultrasonic transducers 21 of the sound wave transmitting device 20 and the arrangement direction of the ultrasonic transducers 31 of the sound wave receiving device 30 are different from each other.
  • the array direction (X direction) of the ultrasonic transducers 21 is orthogonal to the array direction (Y direction) of the ultrasonic transducers 31.
  • the sound pressure distribution of the ultrasonic beam transmitted from the sound wave transmitting device 20 spreads in the direction (Y direction) orthogonal to the array direction (X direction) of the ultrasonic transducers 21. This ultrasonic beam is reflected by the reflecting member in the target space SP.
  • the sound wave receiving device 30 receives the ultrasonic beam reflected by the reflecting member.
  • the sensitivity distribution of the ultrasonic transducers 31 to the ultrasonic waves of the ultrasonic transducers 31 spreads in the direction (X direction) orthogonal to the array direction (Y direction) of the ultrasonic transducers 31.
  • waveform data corresponding to the vibration of the ultrasonic transducer 31 having the sensitivity distribution spreading in the X direction can be obtained for the ultrasonic beam having the sound pressure distribution spreading in the Y direction. That is, the sound wave receiving device 30 most strongly acquires the component at the point where the sound pressure distribution and the sensitivity distribution intersect. This can increase the resolution of space temperature measurement.
  • Modification 2 is an example in which the transmission direction is changed without using the actuator 23.
  • the control circuit 22 of the second modification calculates the phase difference required to realize the transmission angle with respect to the Z axis, according to the value of the transmission angle included in the ultrasonic control signal.
  • the control circuit 22 vibrates the plurality of ultrasonic transducers 21 at different timings so that the calculated phase difference of the ultrasonic beams occurs.
  • the difference in the timing of vibration of each ultrasonic transducer 21 forms the phase difference of the ultrasonic waves transmitted from each ultrasonic transducer 21. Thereby, the ultrasonic beam traveling in the direction indicated by the value of the transmission angle included in the ultrasonic control signal is transmitted.
  • the transmission direction can be changed without using the actuator 23 (that is, the mechanical mechanism for changing the transmission direction).
  • Modification 3 is an example in which the filter is corrected according to the temperature of the target space SP.
  • the processor 12 of the modification 3 acquires the reference temperature from the thermometer 50 in step S1113.
  • the processor 12 corrects the identified filter information based on the reference temperature.
  • the processor 12 extracts the component of the ultrasonic beam that has traveled along the target path Pi from the components included in the received waveform data based on the corrected filter information.
  • filtering is executed based on the temperature of the target space SP. This makes it possible to more reliably extract the component of the ultrasonic beam traveling along the target path Pi from the received waveform data. As a result, the accuracy of mesh temperature calculation can be improved.
  • Modification 4 is an example in which the influence of the air flow in the target space SP on the measurement result is reduced.
  • FIG. 21 is a schematic diagram showing the configurations of the sound wave transmitting device and the sound wave receiving device of the fourth modification.
  • FIG. 22 is an explanatory diagram of the principle of the modified example 4.
  • At least two sensor units SUa and SUb are arranged in the target space SP of Modification 4 (FIG. 21A).
  • the sensor units SUa and SUb each include a sound wave transmitting device 20 and a sound wave receiving device 30 (FIG. 21B).
  • the ultrasonic beam transmitted from the sound wave transmitting device 20a of the sensor unit SUa is received by the sound wave receiving device 30b of the sensor unit SUb.
  • the sound wave receiving device 30b generates reception waveform data (an example of “second reception waveform data”) corresponding to the received ultrasonic beam.
  • the ultrasonic beam transmitted from the sound wave transmitting device 20b of the sensor unit SUb is received by the sound wave receiving device 30a of the sensor unit SUa.
  • the sound wave receiving device 30a generates reception waveform data (an example of “first reception waveform data”) corresponding to the received ultrasonic beam.
  • the temperature measuring device 10 calculates the average speed
  • (
  • the temperature measuring device 10 calculates the temperature corresponding to the average velocity
  • the temperature corresponding to the average speed not including the wind speed component between the pair of sensor units SUa and SUb is calculated.
  • the S/N ratio of the measurement result of the temperature of the space can be further improved.
  • Modification 4 an example of using two sensor units SUa and SUb has been described, but the present embodiment is not limited to this.
  • the present embodiment for example, as long as the sound waves pass through substantially the same traveling path, and the condition that the vector sum of the wind speed components included in those sound waves is satisfied, there are three or more kinds of sound speeds or signals corresponding thereto. It is also applicable when the average of is used.
  • Modification 4 an example of using two sensor units SUa and SUb has been described, but the present embodiment is not limited to this. Modification 4 is also applicable to the case where the reflection/scattering path is used once or more.
  • Modification 5 is an example of a temperature measurement algorithm using a time series filter.
  • FIG. 23 is an explanatory diagram of the outline of the modified example 5.
  • the processor 12 of Modification 5 is configured to execute the path temperature calculation model Mpt(t) and the time series filter FIL.
  • the route temperature calculation model Mpt(t) outputs the route temperature PD(t
  • the time series filter FIL outputs the path temperature calculation model Mpt(t) (path temperature PD(t
  • the time series filter FIL includes at least one of the following, for example. ⁇ Kalman filter ⁇ Extended Kalman filter ⁇ Unscented Kalman filter ⁇ Particle filter
  • FIG. 24 is a flowchart of the temperature measurement process of the modified example 5.
  • the temperature measuring device 10 of the modified example 5 executes steps S110 to S113 as in FIG.
  • step S112 If step S112 has not been completed for all target meshes Mt (S113-NO), the temperature measuring device 10 executes step S110.
  • step S112 When step S112 is completed for all target meshes Mt (S113-YES), the temperature measuring device 10 executes time series filtering (S310). Specifically, the processor 12 acquires the reference temperature Tref(t) at time t from the thermometer 50. The processor 12 sets the path temperature Tp(t
  • the S/N ratio of the measurement result of the temperature of the space can be further improved by executing the time series filtering.
  • the time-series filter FIL of Modification 5 may further calculate the temperature distribution D(t) at time t by referring to the external environment information at time t-1.
  • the external environment information at time t-1 includes, for example, the following information. -Information on the heat quantity of the air conditioner 40-Information on the outside air temperature around the target space SP-Information on the three-dimensional shape of the target space SP-Information on the heat insulation performance of the target space SP-The number of people existing in the target space SP Information-Information on the movement of a person existing in the target space SP-Information on the wind of the air conditioner 40-Information on the wind in the target space SP
  • Modification 6 is an example of a combination of N (N is an integer of 2 or more) sound wave transmitting devices 20 and M (M is an integer of 2 or more) sound wave receiving devices 30.
  • FIG. 25 is a schematic diagram of a sound wave transmitting device and a sound wave receiving device of the sixth modification.
  • N sound wave transmitters 20(1) to 20(N) transmit ultrasonic beams.
  • the M acoustic wave receiving devices 30(1) to 30(M) receive the ultrasonic beams transmitted from the N acoustic wave transmitting devices 20(1) to 20(N), respectively, and Received waveform data corresponding to the beam is generated.
  • the processor 12 of the first example of the modified example 6 executes the filtering (S1113) of FIG. 13 on the reception waveform data generated by the sound wave receiving devices 30(1) to 30(M) to obtain the reception waveform.
  • the path of the ultrasonic beam corresponding to the data that is, the combination of the acoustic wave transmitting device 20(n) (n is an integer of 2 to N) and the acoustic wave receiving device 30(m) that is the source of the ultrasonic beam) Identify.
  • the processor 12 of the second example of the modified example 6 converts the received waveform data into an envelope and specifies the rising time of the envelope.
  • the processor 12 determines the path of the ultrasonic beam corresponding to the received waveform data (that is, the combination of the sound wave transmitting device 20(n) and the sound wave receiving device 30(m), which are the sources of the ultrasonic beam, depending on the time. Identify.
  • the sound wave transmitter 20(n) of the third example of the modification 6 transmits ultrasonic wave beams having different frequencies.
  • the processor 12 refers to the frequency of the received waveform data and refers to the path of the ultrasonic beam corresponding to the received waveform data (that is, the sound wave transmitting device 20(n) and the sound wave receiving device 30 (which are the source of the ultrasonic beam). m) combination) is specified.
  • one sound wave receiving device 30 receives the ultrasonic beam transmitted from the plurality of sound wave transmitting devices 20(n).
  • the temperature measuring device 10 identifies the sound wave transmitting device 20(n) that oscillated the ultrasonic beam by filtering the received waveform data. Thereby, the S/N ratio of the measurement result of the temperature of the space can be improved.
  • Modification 7 Modification 7 will be described. Modification 7 is an example of measuring the distribution of the wind vector in addition to the temperature.
  • the temperature measuring device 10 (not shown) of FIG. 5, at least two sound wave transmitting devices 20a to 20b, and at least two sound wave receiving devices 30a to 30b are arranged. ing.
  • the temperature measurement device 10 controls the sound wave transmission device 20a so as to transmit an ultrasonic beam having a predetermined transmission frequency Fs.
  • the sound wave transmission device 20a transmits an ultrasonic beam having a transmission frequency Fs under the control of the temperature measurement device 10.
  • the sound wave receiving device 30a When receiving the ultrasonic beam, the sound wave receiving device 30a generates received waveform data.
  • the Doppler effect caused by the wind between the sound wave transmitting device 20a and the sound wave receiving device 30a is generated in the ultrasonic beam received by the sound wave receiving device 30a. Therefore, the reception frequency Fra of the ultrasonic beam received by the sound wave receiving device 30a is different from the transmission frequency Fs.
  • the temperature measuring device 10 acquires the reception waveform data from the sound wave receiving device 30a, and refers to the reception waveform data to specify the reception frequency Fra.
  • the relationship of Expression 3.1 holds between the transmission frequency Fs and the reception frequency Fra.
  • Fs (C+Vwa)/C ⁇ Fra (Equation 3.1)
  • Vwa theoretical wind velocity on the path between the sound wave transmitting device 20a and the sound wave receiving device 30a
  • the theoretical wind speed Vwa can be expressed as the formula 3.2.
  • the processor 12 calculates the theoretical wind speed Vwa using equation 3.1.
  • Vwa C ⁇ Fs/Fra ⁇ C (Equation 3.2)
  • the processor 12 uses Equation 3.3 to calculate the theoretical wind speed Vwb on the path between the sound wave transmitting device 20b and the sound wave receiving device 30b.
  • Vwb C ⁇ Fs/Frb ⁇ C (Equation 3.3)
  • the storage device 11 stores a correlation function G(x) of temperature and wind speed in advance.
  • the processor 12 uses the mesh temperature TEMPmesh and the correlation function obtained in step S112 to calculate the corrected wind velocities Vrwa and Vrwb as in equations 3.4 to 3.5.
  • Vrwa G(Tmesh) ⁇ Vwa (Equation 3.4)
  • Vrwb G(Tmesh) ⁇ Vwb (Equation 3.5)
  • the temperature measuring device 10 obtains from the formulas 3.4 to 3.5 by referring to the angles formed by the traveling path of the sound wave transmitting device 20a to the sound wave receiving device 30a and the traveling path of the sound wave transmitting device 20b to the sound wave receiving device 30b.
  • the wind vector is calculated by synthesizing the corrected wind velocities Vrwa and Vrwb.
  • the wind vector of the mesh including the intersection where the traveling paths of the ultrasonic waves intersect is obtained. Therefore, when a plurality of intersections are formed, the distribution of wind vectors in the target space SP (that is, the distribution of wind speed and wind direction) is obtained.
  • the wind vector obtained has two dimensions.
  • the modified example 7 is also applicable to a case where three traveling routes form an intersection. In this case, the obtained wind vector has three dimensions.
  • Modification 8 is an example of measuring the distribution of wind vectors using ultrasonic waves. 26: is a figure which shows an example of the sensor arrangement of the modification 8. FIG.
  • At least four sensor units SUa to SUd have ultrasonic wave traveling paths of the sensor units SUa to SUb and ultrasonic wave traveling of the sensor units SUc to SUd.
  • the paths are arranged so as to intersect with each other (preferably orthogonal to each other).
  • the configuration of each sensor unit SUa to SUd is the same as that of the modification 4 (FIG. 21B).
  • the temperature measuring device 10 uses the equation 4.1 to remove the temperature factor from the absolute value
  • Vwab
  • the temperature measuring device 10 uses the equation 4.2 to remove the temperature factor from the absolute value
  • the wind speed component Vwcd during the period is calculated.
  • Vwcd
  • the temperature measuring device 10 refers to the angles formed by the traveling paths of the sensor units SUa to SUb and the traveling paths of the sensor units SUc to SUd, and obtains the wind speed component Vwab obtained from the equation 4.1 and the equation obtained from the equation 4.2.
  • the wind vector is calculated by synthesizing the obtained wind speed component Vwcd.
  • the wind vector of the mesh in which the traveling path of the ultrasonic waves of the sensor unit arranged in the target space SP includes the intersection is obtained. Therefore, when a plurality of intersections are formed, the distribution of the wind vector (that is, the wind speed and the wind direction) is obtained instead of the temperature distribution of the target space SP.
  • the wind vector is calculated by using the absolute values
  • the pair of sensor units forms an intersection.
  • the wind vector obtained has two dimensions.
  • the modified example 8 is also applicable to a case where three sensor units form an intersection. In this case, the obtained wind vector has three dimensions.
  • Modification 9 Modification 9 will be described. Modification 9 is an example of measuring temperature using a transducer array.
  • the sound wave transmitting device 20 includes a plurality of ultrasonic vibrators 21 and a control circuit 22.
  • each ultrasonic transducer 21 transmits an ultrasonic beam USW traveling in the Z direction.
  • the sound wave receiving device 30 includes a plurality of ultrasonic vibrators 31 and a control circuit 32.
  • the plurality of ultrasonic transducers 31 are two-dimensionally arranged on the transmission surface (XY plane). That is, the plurality of ultrasonic transducers 31 form a transducer array TA.
  • each ultrasonic transducer 31 vibrates when receiving the ultrasonic beam USW transmitted from the acoustic wave transmission device 20.
  • the first example of Modification 9 is an example in which the temperature is measured using a pair of the sound wave receiving device 30 and the sound wave transmitting device 20.
  • FIG. 29 is a diagram showing an outline of the first example of the ninth modification.
  • the transducers 21a to 21c of the sound wave transmitting device 20 of the first example of the modification 9 respectively transmit the ultrasonic beams USW0 to USW2 in the Z+ direction.
  • the ultrasonic beams USW0 to USW2 shift in the X+ direction under the influence of the airflow AF. As a result, the ultrasonic beam USW2 travels toward the outside of the acoustic wave receiving device 30. On the other hand, the ultrasonic beams USW0 to USW1 are received by the ultrasonic transducers 31b to 31c.
  • the sound wave transmitting device 20 and the sound wave receiving device 30 both form the transducer array TA.
  • the ultrasonic beams USW0 to USW2 emitted from the transducer array TA of the sound wave transmitting device 20 can easily reach the transducer array TA of the sound wave receiving device 30 even if affected by the airflow AF.
  • the same effect as this embodiment can be obtained without being affected by the airflow AF.
  • the second example of Modification 9 is an example of measuring the temperature using a pair of sensor units SU (a combination of the sound wave receiving device 30 and the sound wave transmitting device 20).
  • 30: is a figure which shows the structure of the sensor unit of the modification 9.
  • FIG. 31: is a figure which shows the outline
  • the sensor unit SU of the second example of modification 9 includes a sound wave transmitting device 20 and a sound wave receiving device 30.
  • the sound wave transmission device 20 includes a plurality of ultrasonic vibrators 21 and a control circuit 22.
  • a plurality of ultrasonic transducers 21 are two-dimensionally arranged on the transmission surface (XY plane). That is, the ultrasonic transducers 21 form a transducer array TA.
  • the sound wave receiving device 30 includes a plurality of ultrasonic vibrators 31 and a control circuit 32.
  • the plurality of ultrasonic transducers 31 are two-dimensionally arranged on the transmission surface (XY plane). That is, the plurality of ultrasonic transducers 31 form a transducer array TA.
  • a pair of sensor units SUa and SUb are arranged in the target space SP.
  • the sound wave transmission device 20 of the sensor unit SUa emits ultrasonic waves toward the sound wave reception device 30 of the sensor unit SUb.
  • the sound wave transmitting device 20 of the sensor unit SUb emits ultrasonic waves toward the sound wave receiving device 30 of the sensor unit SUa.
  • the transducers 21a to 21c of the sound wave transmitting device 20 of the sensor unit SUa respectively transmit the ultrasonic beams Uswa0 to Uswa2 in the Z+ direction.
  • the transducers 21a to 21c of the sound wave transmitting device 20 of the sensor unit SUb respectively transmit the ultrasonic beams USWb0 to USWb2 in the Z-direction.
  • the ultrasonic beams Uswa0 to Uswa2 and uswb0 to uswb2 shift in the X+ direction under the influence of the airflow AF.
  • the ultrasonic beam Uswa2 travels toward the outside of the acoustic wave receiving device 30 of the sensor unit SUb.
  • the ultrasonic beam USWb2 travels toward the outside of the acoustic wave receiving device 30 of the sensor unit SUa.
  • the ultrasonic beams Uswa0 to Uswa1 are received by the ultrasonic transducers 31b to 31c of the sound wave receiving device 30 of the sensor unit SUb.
  • the ultrasonic beams USWb0 to USWb1 are received by the ultrasonic transducers 31b to 31c of the acoustic wave receiving device 30 of the sensor unit SUa.
  • the sound wave transmitting device 20 and the sound wave receiving device 30 both form the transducer array TA.
  • the ultrasonic beams USW0 to USW2 emitted from the transducer array TA of the sound wave transmitting device 20 can easily reach the transducer array TA of the sound wave receiving device 30 even if affected by the airflow AF.
  • the temperature of the space between the sensor units SUa and SUb is measured using both the ultrasonic beam traveling in the Z- direction and the ultrasonic beam traveling in the Z+ direction. As a result, the same effect as this embodiment can be obtained without being affected by the airflow AF.
  • Modification 10 is an example in which a waveform to be referred to is specified from among the waveforms of the sound waves received by one sound wave receiving device 30.
  • FIG. 32 is an explanatory diagram of a filter according to the modification 10.
  • the sound wave receiving device 30 of the first example of modification 1 has a time filter for each path.
  • the temporal filter is, for example, at least one of the following (FIG. 32).
  • ⁇ Lower limit time threshold THtb ⁇ Upper limit time threshold THtt
  • the lower limit time threshold THtb and the upper limit time threshold THtt are determined by at least one of the following. -Distance between the sound wave transmitting device 20 and the sound wave receiving device 30-Distance that the ultrasonic beam propagates before the ultrasonic wave beam transmitted from the sound wave transmitting device 20 is received by the sound wave receiving device 30
  • step S1112 the control circuit 32 applies the time filter (for example, the time window Wt) to the received waveform data after generating the received waveform data (FIG. 32) according to the vibration of the ultrasonic transducer 31.
  • the waveform WF2 to be referred to is extracted from the waveforms WF1 to WF3 included in the received waveform data.
  • the control circuit 32 transmits the received waveform data including the extracted waveform WF2 to the temperature measuring device 10.
  • the temperature measurement device 10 of the second example of the modification 10 has the same time filter as that of the first example.
  • step S1113 the processor 12 applies a time filter (for example, the time window Wt) to the received waveform data, thereby selecting the waveform WF2 to be referenced from the waveforms WF1 to WF3 included in the received waveform data. Extract.
  • a time filter for example, the time window Wt
  • the storage device 11 may be connected to the temperature measuring device 10 via the network NW.
  • the sound wave receiving device 30 including the ultrasonic transducer 31 is shown.
  • the sound wave receiving device 30 may include a plurality of ultrasonic transducers 31 as with the sound wave transmitting device 20.
  • one acoustic wave transmission device 20 transmits ultrasonic beams along a plurality of paths
  • one acoustic wave reception device 30 receives ultrasonic beams along a plurality of paths.
  • n is an integer of 2 or more
  • sound wave transmission devices 20 has an ultrasonic beam along one path (that is, n sound wave transmission devices 20 have ultrasonic beams along n paths).
  • each of the n sound wave receiving devices 30 may receive an ultrasonic beam along each path (that is, the n sound wave receiving devices 30 transmit the ultrasonic waves along the n paths). Beam may be received).
  • the sound wave transmitting device 20 may send an ultrasonic beam including an autocorrelation signal (for example, an M-sequence signal) having a relatively strong autocorrelation. Thereby, the S/N ratio of the measurement result of the temperature of the space can be further improved.
  • an autocorrelation signal for example, an M-sequence signal
  • the sound wave transmitting device 20 may individually transmit the ultrasonic beam including the different autocorrelation signals, so that the sound wave receiving device 30 may identify the sound wave transmitting device 20 that is the source of the ultrasonic beam. Alternatively, the sound wave receiving device 30 may identify the sound wave transmitting device 20 that is the transmission source of the ultrasonic beam by transmitting an ultrasonic beam having an oscillation frequency different for each sound wave transmitting device 20.
  • the temperature measuring device 10 can measure the following air characteristic distribution. ⁇ Distribution of concentration of chemical substances (eg CO 2 ) in air ⁇ Distribution of humidity ⁇ Distribution of odor ⁇ Distribution of toxic gas
  • the sound wave transmitting device 20 and the sound wave receiving device 30 are defined separately, but the scope of the present embodiment is not limited to this.
  • one ultrasonic transducer may have a function of transmitting ultrasonic waves and a function of receiving ultrasonic waves.
  • At least one of the equation used to calculate the route temperature TEMPpathi in step S1114 (FIG. 13) and the equation used to calculate the mesh temperature TEMPmesht in step S112 (FIG. 12) is the external environment information (for example, At least one of the outside air temperature, the outside air humidity, and the outside air pressure) may be included as a parameter.
  • the S/N ratio of the measurement result of the air characteristic of the space can be improved regardless of the external environment information.
  • the sound wave transmitting device 20 shows an example of transmitting an ultrasonic beam having directivity, but the present embodiment is not limited to this.
  • the present embodiment is also applicable to the case where the sound wave transmission device 20 transmits an audible sound beam (that is, a sound wave having a frequency different from that of the ultrasonic beam).
  • the temperature distribution is not limited to the mesh temperature TEMPmesh.
  • the temperature distribution also includes at least one of the following: ⁇ Temperature at multiple points on the route ⁇ Average temperature on the route

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

対象空間に配置された音波送信装置、及び、対象空間に配置された音波受信装置と接続可能な温度計測装置は、音波を送信させるように、音波送信装置を制御する手段を備え、超音波ビームを受信する音波受信装置から、音波の波形に関する受信波形データを取得する手段を備え、受信波形データに基づいて、対象空間の温度分布を計算する手段を備える。

Description

温度計測装置、計測装置、音波受信装置、及び、プログラム
　本発明は、温度計測装置、計測装置、音波受信装置、及び、プログラムに関する。
　空気中を伝搬する音波の速度が温度に応じて変化する原理を用いて、音波の伝搬時間から空間の温度を計測することが可能である。
　例えば、特許文献１は、超音波を送受信可能な複数のセンサユニットを空間に配置することにより、超音波の伝搬時間から空間の温度を計測する技術を開示している。
特開２０１４－０９５６００号公報
　特許文献１では、超音波に指向性がないため、センサユニットは、所望の超音波成分だけでなく、空間に存在する反射部材（例えば、壁）で反射した超音波の成分を受信してしまう。そのため、超音波を誤検出する可能性がある。
　つまり、従来、空間の温度を計測する場合のＳ／Ｎ比が低い。その結果、空間の温度の計測結果の精度が低下する。
　本発明の目的は、空間の空気特性（例えば、温度）の計測結果のＳ／Ｎ比を向上させることである。
　本発明の一態様は、
　対象空間に配置された音波送信装置、及び、前記対象空間に配置された音波受信装置と接続可能な温度計測装置であって、
　音波を送信させるように、前記音波送信装置を制御する手段を備え、
　前記音波を受信する音波受信装置から、前記音波の波形に関する受信波形データを取得する手段を備え、
　前記受信波形データに基づいて、前記対象空間の温度分布を計算する手段を備える、
温度計測装置である。
　本発明によれば、空間の空気特性の計測結果のＳ／Ｎ比を向上させることができる。
第１実施形態の計測システムの構成を示すブロック図である。 第１実施形態の計測システムの詳細な構成を示すブロック図である。 第１実施形態の音波送信装置の構成を示す概略図である。 第１実施形態の音波受信装置の構成を示す概略図である。 第１実施形態の概要の説明図である。 第１実施形態の空間データテーブルのデータ構造を示す図である。 第１実施形態のセンサデータテーブルのデータ構造を示す図である。 第１実施形態の経路データテーブルのデータ構造を示す図である。 第１実施形態のメッシュデータテーブルのデータ構造を示す図である。 第１実施形態のフィルタの説明図である。 第１実施形態のセンサ配置の一例を示す図である。 第１実施形態の温度計測の処理のフローチャートである。 図１２の経路温度の計算の詳細なフローチャートである。 図１２の処理において表示される画面例を示す図である。 第２実施形態の音波送信装置の構成を示す概略図である。 第２実施形態の経路データテーブルのデータ構造を示す図である。 第２実施形態の送信角の説明図である。 第２実施形態のセンサ配置の一例を示す図である。 第２実施形態の経路温度の計算の詳細なフローチャートである。 変形例１の音波送信装置及び音波受信装置の構成を示す概略図である。 変形例４の音波送信装置及び音波受信装置の構成を示す概略図である。 変化例４の原理の説明図である。 変形例５の概要の説明図である。 変形例５の温度計測の処理のフローチャートである。 変形例６の音波送信装置及び音波受信装置の概略図である。 変形例８のセンサ配置の一例を示す図である。 変形例９の音波送信装置の構成を示す図である。 変形例９の音波受信装置の構成を示す図である。 変形例９の第１例の概要を示す図である。 変形例９のセンサユニットの構成を示す図である。 変形例９の第２例の概要を示す図である。 変形例１０のフィルタの説明図である。
　以下、本発明の一実施形態について、図面に基づいて詳細に説明する。なお、実施形態を説明するための図面において、同一の構成要素には原則として同一の符号を付し、その繰り返しの説明は省略する。
（１）第１実施形態
　第１実施形態を説明する。
（１－１）計測システムの構成
　第１実施形態の計測システムの構成について説明する。図１は、第１実施形態の計測システムの構成を示すブロック図である。図２は、第１実施形態の計測システムの詳細な構成を示すブロック図である。
　図１及び図２に示すように、計測システム１は、温度計測装置１０と、音波送信装置２０と、音波受信装置３０と、空調装置４０と、温度計５０と、を備える。
　温度計測装置１０は、音波送信装置２０、音波受信装置３０、空調装置４０、及び、温度計５０に接続されている。
　温度計測装置１０、音波送信装置２０、音波受信装置３０、空調装置４０、及び、温度計５０は対象空間ＳＰに配置されている。
　温度計測装置１０は、以下の機能を備える。
　・音波送信装置２０を制御する機能
　・音波受信装置３０から受信波形データを取得する機能
　・対象空間ＳＰの温度分布を計測する機能
　・空調装置４０を制御する機能
　・温度計５０から対象空間ＳＰの温度の測定結果に関する基準温度情報を取得する機能
　温度計測装置１０は、例えば、スマートフォン、タブレット端末、又は、パーソナルコンピュータである。
　音波送信装置２０は、温度計測装置１０の制御に従い、指向性を有する超音波ビーム（「音波」の一例）を送信するように構成される。また、音波送信装置２０は、超音波の送信方向を変更するように構成される。
　音波受信装置３０は、音波送信装置２０から送信された超音波ビームを受信し、且つ、受信した超音波ビームに応じた受信波形データを生成するように構成される。音波受信装置３０は、例えば、無指向性マイクロフォン又は指向性マイクロフォンである。
　空調装置４０は、温度計測装置１０の制御に従い、対象空間ＳＰの温度を調整するように構成される。
　温度計５０は、対象空間ＳＰの温度（以下「基準温度」という）を測定するように構成される。
（１－１－１）温度計測装置の構成
　第１実施形態の温度計測装置１０の構成について説明する。
　図２に示すように、温度計測装置１０は、記憶装置１１と、プロセッサ１２と、入出力インタフェース１３と、通信インタフェース１４と、を備える。
　記憶装置１１は、プログラム及びデータを記憶するように構成される。記憶装置１１は、例えば、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、及び、ストレージ（例えば、フラッシュメモリ又はハードディスク）の組合せである。
　プログラムは、例えば、以下のプログラムを含む。
　・ＯＳ（Operating System）のプログラム
　・情報処理（例えば、対象空間ＳＰの温度分布を計測するための情報処理）を実行するアプリケーションのプログラム
　・空間の温度に対する音波の速度に関する音波速度特性に関するデータ
　データは、例えば、以下のデータを含む。
　・情報処理において参照されるデータベース
　・情報処理を実行することによって得られるデータ（つまり、情報処理の実行結果）
　プロセッサ１２は、記憶装置１１に記憶されたプログラムを起動することによって、温度計測装置１０の機能を実現するように構成される。プロセッサ１２は、コンピュータの一例である。
　入出力インタフェース１３は、温度計測装置１０に接続される入力デバイスからユーザの指示を取得し、かつ、温度計測装置１０に接続される出力デバイスに情報を出力するように構成される。
　入力デバイスは、例えば、キーボード、ポインティングデバイス、タッチパネル、又は、それらの組合せである。また、入力デバイスは、温度計５０を含む。
　出力デバイスは、例えば、ディスプレイである。また、出力デバイスは、空調装置４０を含む。
　通信インタフェース１４は、外部装置（例えば、サーバ）との間の通信を制御するように構成される。
（１－１－２）音波送信装置の構成
　第１実施形態の音波送信装置２０の構成を説明する。図３は、第１実施形態の音波送信装置の構成を示す概略図である。
　図３Ａに示すように、音波送信装置２０は、複数の超音波振動子（「振動素子」の一例）２１と、制御回路２２と、を備える。
　図３Ｂに示すように、制御回路２２は、温度計測装置１０の制御に従って、複数の超音波振動子２１を振動させる。複数の超音波振動子２１が振動すると、送信面（ＸＹ平面）に対して直交する送信方向（Ｚ軸方向）に向かって、超音波ビームが送信される。
（１－１－３）音波受信装置の構成
　第１実施形態の音波受信装置３０の構成を説明する。図４は、第１実施形態の音波受信装置の構成を示す概略図である。
　図４に示すように、音波受信装置３０は、超音波振動子３１と、制御回路３２と、を備える。
　超音波振動子３１は、音波送信装置２０から送信された超音波ビームを受信すると振動する。
　制御回路３２は、超音波振動子３１の振動に応じた受信波形データを生成するように構成される。
（１－２）実施形態の概要
　第１実施形態の概要について説明する。図５は、第１実施形態の概要の説明図である。
　図５に示すように、温度計測の対象となる空間（以下「対象空間」という）ＳＰには、温度計測装置１０（不図示）と、音波送信装置２０ａ～２０ｂと、音波受信装置３０ａ～３０ｂと、が配置されている。温度計測装置１０は、音波送信装置２０及び音波受信装置３０と接続可能である。
　温度計測装置１０は、音波を送信させるように、音波送信装置２０ａ～２０ｂを制御する。
　温度計測装置１０は、音波受信装置３０ａ～３０ｂから、受信された音波の波形に関する受信波形データを取得する。
　温度計測装置１０は、受信波形データに基づいて、対象空間ＳＰの温度分布を計算する。
　本実施形態によれば、音波の伝搬時間から対象空間ＳＰの温度分布を計算する。これにより、空間の温度の計測結果のＳ／Ｎ比を向上させることができる。
（１－３）データテーブル
　第１実施形態のデータテーブルについて説明する。
（１－３－１）空間データテーブル
　第１実施形態の空間データテーブルについて説明する。図６は、第１実施形態の空間データテーブルのデータ構造を示す図である。
　図６の空間データテーブルには、音波送信装置２０及び音波受信装置３０が配置された空間（以下「対象空間」という）に関する空間情報が格納される。
　空間データテーブルは、「座標」フィールドと、「反射特性」フィールドと、を含む。各フィールドは、互いに関連付けられている。
　「座標」フィールドには、対象空間に存在する反射部材の座標（以下「反射部材座標」という）が格納される。反射部材座標は、対象空間の任意の基準点を原点とする座標系（以下「空間座標系」という）で表される。
　「反射特性」フィールドには、反射部材の反射特性に関する反射特性情報が格納される。「反射特性」フィールドは、「反射種別」フィールドと、「反射率」フィールドと、「法線角」フィールドと、を含む。
　「反射種別」フィールドには、反射種別に関する情報が格納される。反射種別は、以下の何れかである。
　・拡散反射
　・鏡面反射
　「反射率」フィールドには、反射部材の反射率の値が格納される。
　「法線角」フィールドには、反射部材の反射面の法線角度の値が格納される。
（１－３－２）センサデータテーブル
　第１実施形態のセンサデータテーブルについて説明する。図７は、第１実施形態のセンサデータテーブルのデータ構造を示す図である。
　図７に示すように、センサデータテーブルには、音波送信装置２０及び音波受信装置３０に関する情報（以下「センサ情報」という）が格納される。
　センサデータテーブルは、「センサＩＤ」フィールドと、「座標」フィールドと、「センサタイプ」フィールドと、を含む。
　各フィールドは、互いに関連付けられている。
　「センサＩＤ」フィールドには、音波送信装置２０又は音波受信装置３０を識別するセンサ識別情報が格納される。
　「座標」フィールドには、音波送信装置２０又は音波受信装置３０の位置を示す座標（以下「センサ座標」という）が格納される。センサ座標は、空間座標系で表される。
　「センサタイプ」フィールドには、音波送信装置２０であることを示すタグ「送信」、又は、音波受信装置３０であることを示すタグ「受信」が格納される。
（１―３－３）経路データテーブル
　第１実施形態の経路データテーブルについて説明する。図８は、第１実施形態の経路データテーブルのデータ構造を示す図である。
　図８に示すように、経路データテーブルには、経路に関する経路情報が格納される。
　経路データテーブルは、「経路ＩＤ」フィールドと、「送信センサ」フィールドと、「受信センサ」フィールドと、を含む。
　「経路ＩＤ」フィールドには、経路を識別する経路識別情報が格納される。
　「送信センサ」フィールドには、経路を構成する音波送信装置２０のセンサ識別情報が格納される。
　「受信センサ」フィールドには、経路を構成する音波受信装置３０のセンサ識別情報が格納される。
（１－３－４）メッシュデータテーブル
　第１実施形態のメッシュデータテーブルについて説明する。図９は、第１実施形態のメッシュデータテーブルのデータ構造を示す図である。図１０は、第１実施形態のフィルタの説明図である。
　図９に示すように、メッシュデータテーブルには、仮想メッシュに関するメッシュ情報が格納される。
　メッシュデータテーブルは、「メッシュＩＤ」フィールドと、「座標」フィールドと、「経路ＩＤ」フィールドと、「フィルタ」フィールドと、を含む。
　「メッシュＩＤ」フィールドには、仮想メッシュを識別するメッシュ識別情報が格納される。
　「座標」フィールドには、仮想メッシュの位置を示すメッシュ座標が格納される。メッシュ座標は、空間座標系で表される。
　「経路ＩＤ」フィールドには、経路の経路識別情報が格納される。
　図１０に示すように、対象空間ＳＰは、複数の仮想メッシュＭｉ（ｉは引数）に分割される。各仮想メッシュＭｉは、３次元形状を有する。
　例えば、仮想メッシュＭ１は、複数の経路Ｐ１０１及びＰ２００を含む。経路Ｐ１０１は、音波送信装置２０ａから音波受信装置３０ａに至る経路である。経路Ｐ２００は、音波送信装置２０ｂから音波受信装置３０ｂに至る経路である。
　「フィルタ」フィールドには、音波受信装置３０によって受信された受信波形データによって再現される超音波ビームの波形から特定波形を抽出するためのフィルタに関するフィルタ情報が格納される。フィルタ情報は、「経路ＩＤ」フィールドに格納された経路識別情報に関連付けられる。「フィルタ」フィールドは、「時間フィルタ」フィールドと、「振幅フィルタ」フィールドと、を含む。
　「時間フィルタ」フィールドには、時間軸に沿って特定波形を抽出するための時間フィルタに関する情報が格納される。時間フィルタは、例えば、以下の少なくとも１つである（図１０）。
　・下限時間閾値ＴＨｔｂ
　・上限時間閾値ＴＨｔｔ
　・下限時間閾値ＴＨｔｂと上限時間閾値ＴＨｔｔとによって規定される時間ウインドウＷｔ
　「振幅フィルタ」フィールドには、振幅軸に沿って特定波形を抽出するための振幅フィルタに関する情報が格納される。振幅フィルタは、例えば、以下の少なくとも１つである（図１０）。
　・下限振幅閾値ＴＨａｂ
　・上限振幅閾値ＴＨａｔ
　・下限振幅閾値ＴＨａｂと上限振幅閾値ＴＨａｔとによって規定される振幅ウインドウＷａ
（１－４）温度計測の処理
　第１実施形態の温度計測の処理について説明する。図１１は、第１実施形態のセンサ配置の一例を示す図である。図１２は、第１実施形態の温度計測の処理のフローチャートである。図１３は、図１２の経路温度の計算の詳細なフローチャートである。図１４は、図１２の処理において表示される画面例を示す図である。
　図１１に示すように、対象空間ＳＰには、複数の音波送信装置２０ａ～２０ｅと、複数の音波受信装置３０ａ～３０ｅと、が配置される。
　複数の音波送信装置２０ａ～２０ｅは、それぞれ、複数の音波受信装置３０ａ～３０ｅに対向している。例えば、音波送信装置２０ａは、音波受信装置３０ａに対向している。これは、音波送信装置２０ａ及び音波受信装置３０ａが、センサペアを形成することを意味している。
　図１１は、５つのセンサペアが形成される例を示している。
　複数の経路を含む仮想メッシュのメッシュ温度の計測が可能である。
　図１１は、仮想メッシュＭ１～Ｍ４のメッシュ温度の計測が可能である例を示している。
　図１２に示すように、温度計測装置１０は、対象メッシュの決定（Ｓ１１０）を実行する。
　具体的には、プロセッサ１２は、図５に示すように、対象空間ＳＰを構成する複数の仮想メッシュの中から対象メッシュＭｔ（ｔ＝１～４）のメッシュ識別情報を決定する。
　ステップＳ１１０の後、温度計測装置１０は、所定の経路温度計算モデルに従って、経路温度の計算（Ｓ１１１）を実行する。
　図１３を参照して、ステップＳ１１１の詳細を説明する。
　温度計測装置１０は、対象経路の決定（Ｓ１１１０）を実行する。
　具体的には、プロセッサ１２は、メッシュデータテーブル（図９）を参照して、ステップＳ１１０で決定したメッシュ識別情報に関連付けられた「経路ＩＤ」フィールドの情報（つまり、対象メッシュＭｔを通る経路（以下「対象経路」という）Ｐｉ（ｉは、経路の引数）の経路識別情報）を特定する。
　ステップＳ１１１０の後、温度計測装置１０は、超音波ビームの出力（Ｓ１１１１）を実行する。
　具体的には、プロセッサ１２は、経路データテーブル（図８）を参照して、ステップＳ１１１０で特定した経路識別情報に関連付けられた「送信センサ」フィールドの情報（つまり、制御対象となる音波送信装置（以下、「対象音波送信装置」という）２０）と、「受信センサ」フィールドの情報（つまり、制御対象となる音波受信装置（以下「対象音波受信装置」という）３０）と、を特定する。
　プロセッサ１２は、対象音波送信装置２０に超音波制御信号を送信する。
　対象音波送信装置２０は、温度計測装置１０から送信された超音波制御信号に応じて超音波ビームを送信する。
　具体的には、複数の超音波振動子２１は、超音波制御信号に応じて同時に振動する。
　これにより、対象音波送信装置２０から対象音波受信装置３０に向かって、送信方向（Ｚ軸方向）に進行する超音波ビームが送信される。
　ステップＳ１１１１の後、温度計測装置１０は、受信波形データの取得（Ｓ１１１２）を実行する。
　具体的には、対象音波受信装置３０の超音波振動子３１は、ステップＳ１１１１で対象音波送信装置２０から送信された超音波ビームを受信することにより振動する。
　制御回路３２は、超音波振動子３１の振動に応じた受信波形データ（図１０）を生成する。
　制御回路３２は、生成した受信波形データを温度計測装置１０に送信する。
　温度計測装置１０のプロセッサ１２は、音波受信装置３０から送信された受信波形データを取得する。
　ステップＳ１１１２の後、温度計測装置１０は、フィルタリング（Ｓ１１１３）を実行する。
　具体的には、プロセッサ１２は、メッシュデータテーブル（図９）を参照して、ステップＳ１１０で決定した対象経路Ｐｉの経路識別情報に関連付けられた「フィルタ」フィールドを特定する。
　例えば、対象メッシュのメッシュ識別情報が「Ｍ００１」である場合、以下のフィルタ情報が特定される。
　・経路識別情報「Ｐ００１」：時間閾値ＴＨｔ１以上及び振幅閾値ＴＨａ１以上
　・経路識別情報「Ｐ００２」：時間閾値ＴＨｔ２以上及び振幅ウインドウＷａ２内
　・経路識別情報「Ｐ００３」：時間ウインドウＷｔ３内及び振幅閾値ＴＨａ３以上
　・経路識別情報「Ｐ００４」：時間ウインドウＷｔ４内及び振幅ウインドウＷａ４内
　プロセッサ１２は、特定したフィルタ情報に基づいて、受信波形データに含まれる成分のうち、対象経路Ｐｉに沿って進行した超音波ビームの成分を抽出する。
　ステップＳ１１１３の後、温度計測装置１０は、経路温度の計算（Ｓ１１１４）を実行する。
　具体的には、プロセッサ１２は、センサデータテーブル（図７）の「座標」フィールドを参照して、センサペア毎に、センサペアを構成する音波送信装置２０の座標及び音波受信装置３０の座標を特定する。
　プロセッサ１２は、特定した音波送信装置２０の座標及び音波受信装置３０の座標の組合せに基づいて、当該音波送信装置２０と当該音波受信装置３０との間の距離（以下「センサ間距離」という）Ｄｓを計算する。
　プロセッサ１２は、ステップＳ１１１３で抽出された成分のピーク値に対応する時間（以下「伝搬時間」という）ｔを特定する。伝搬時間ｔは、音波送信装置２０が超音波ビームを送信してから、対象経路Ｐｉに沿って進行した超音波ビームが音波受信装置３０に到達するまでの所要時間（つまり、対象経路の始点から終点までを超音波ビームが伝搬する時間）を意味する。
　プロセッサ１２は、超音波の音速Ｃ、センサ間距離Ｄｓと、伝搬時間ｔと、基準温度Ｔ０と、を用いて、対象経路Ｐｉの経路温度ＴＥＭＰｐａｔｈｉを計算する。
　全ての対象経路ＰｉについてステップＳ１１１４が終了していない場合（Ｓ１１１５－ＮＯ）、温度計測装置１０は、ステップＳ１１１０を実行する。
　全ての対象経路ＰｉについてステップＳ１１１４が終了すると（Ｓ１１１５－ＹＥＳ）、温度計測装置１０は、図１２のメッシュ温度の計算（Ｓ１１２）を実行する。
　具体的には、プロセッサ１２は、ステップＳ１１１４（図１１）において計算された全ての対象経路Ｐｉの経路温度ＴＥＭＰｐａｔｈｉを用いて、対象メッシュＭｔのメッシュ温度ＴＥＭＰｍｅｓｈｔを計算する（式１）。
　ＴＥＭＰｍｅｓｈｔ　＝　ＡＶＥ（ＴＥＭＰｐａｔｈｉ）…（式１）
　・ＡＶＥ（ｘ）：ｘの平均値を求める関数
　全ての対象メッシュＭｔについてステップＳ１１２が終了していない場合（Ｓ１１３－ＮＯ）、温度計測装置１０は、ステップＳ１１０を実行する。
　全ての対象メッシュＭｔについてステップＳ１１２が終了すると（Ｓ１１３－ＹＥＳ）、温度計測装置１０は、計測結果の提示（Ｓ１１４）を実行する。
　具体的には、プロセッサ１２は、画面Ｐ１０（図１４）をディスプレイに表示する。
　画面Ｐ１０は、表示オブジェクトＡ１０を含む。
　表示オブジェクトＡ１０には、画像ＩＭＧ１０が表示される。
　画像ＩＭＧ１０は、対象空間ＳＰを構成する複数の仮想メッシュのそれぞれについて、ステップＳ１１２で計算されたメッシュ温度ＴＥＭＰｍｅｓｈｔを示している。
　第１実施形態によれば、超音波ビームの伝搬時間から対象空間ＳＰの温度分布を計算する。これにより、空間の温度の計測結果のＳ／Ｎ比を向上させることができる。
（２）第２実施形態
　第２実施形態を説明する。第２実施形態は、音波送信装置２０の超音波ビームの送信方向が可変である例である。
（２－１）音波送信装置の構成
　音波送信装置２０の構成を説明する。図１５は、第２実施形態の音波送信装置の構成を示す概略図である。
　図１５に示すように、音波送信装置２０は、複数の超音波振動子２１と、制御回路２２と、アクチュエータ２３と、を備える。
　図１５Ａに示すように、複数の超音波振動子２１は、送信面（ＸＹ平面）上に二次元配列される。つまり、複数の超音波振動子２１は、振動子アレイＴＡを形成する。
　アクチュエータ２３は、送信軸（Ｚ軸）に対する送信面（ＸＹ平面）の向きを変更するように構成される。
　アクチュエータ２３が送信面を送信軸（Ｚ軸）方向に向けると、超音波ビームＵＳＷ０が送信される。
　アクチュエータ２３が送信面を送信軸（Ｚ軸）に対して傾斜させると、超音波ビームＵＳＷ１が送信される。
（２－２）経路データテーブル
　第２実施形態の経路データテーブルについて説明する。図１６は、第２実施形態の経路データテーブルのデータ構造を示す図である。図１７は、第２実施形態の送信角の説明図である。
　図１６に示すように、経路データテーブルは、図８のフィールド（「経路ＩＤ」フィールド、「送信センサ」フィールド、及び、「受信センサ」フィールド）に加えて、「送信角度」フィールドを含む。
　「送信角度」フィールドには、音波送信装置２０の送信軸（Ｚ軸）に対する超音波ビームの送信角の値が格納される。
　図１７に示すように、経路Ｐ１～Ｐ３は、音波送信装置２０、音波受信装置３０、及び、送信角度θの組合せによって特定される。
（２－３）情報処理
　第２実施形態の情報処理について説明する。図１８は、第２実施形態のセンサ配置の一例を示す図である。図１９は、第２実施形態の経路温度の計算の詳細なフローチャートである。
　図１８に示すように、対象空間ＳＰには、複数の音波送信装置２０ａ～２０ｂと、複数の音波受信装置３０ａ～３０ｂと、が配置される。
　複数の音波送信装置２０ａ～２０ｂは、何れも、複数の音波受信装置３０ａ～３０ｂに到達する経路に沿った超音波ビームを送信可能である。例えば、音波送信装置２０ａは、音波受信装置３０ａに到達する経路Ｐ２０ａに沿った超音波ビームと、音波受信装置３０ｂに到達する経路Ｐ２１ａに沿った超音波ビームと、を送信可能である。
　超音波ビームの経路を含む仮想メッシュのメッシュ温度の計測が可能である。
　対象空間ＳＰの反射部材（例えば、壁及び天井の少なくとも１つ）で超音波ビームが反射する経路Ｐ２０ｂ及びＰ２１ｂを含むので、より多くの仮想メッシュのメッシュ温度の計測が可能である。
　図１９に示すように、第２実施形態の温度計測装置１０は、ステップＳ１１１０の後、送信方向の決定（Ｓ２１１０）を実行する。
　具体的には、プロセッサ１２は、経路データテーブル（図８）を参照して、ステップＳ１１１０で特定した経路識別情報に関連付けられた「送信センサ」フィールドの情報（つまり、制御対象となる音波送信装置（以下、「対象音波送信装置」という）２０）と、「受信センサ」フィールドの情報（つまり、制御対象となる音波受信装置（以下「対象音波受信装置」という）３０）と、「送信角度」フィールドの情報（つまり、対象音波送信装置２０の送信角の値）と、を特定する。
　プロセッサ１２は、対象音波送信装置２０に超音波制御信号を送信する。超音波制御信号は、送信角の値を含む。
　対象音波送信装置２０は、温度計測装置１０から送信された超音波制御信号に含まれる送信角度の値が示す方向に超音波ビームを送信する。
　具体的には、アクチュエータ２３は、超音波制御信号に含まれる送信角の値を参照して、送信軸（Ｚ軸）に対する送信面（ＸＹ平面）の向きを変更する。
　制御回路２２は、複数の超音波振動子２１を同時に振動させる。
　これにより、超音波制御信号に含まれる送信角の値が示す方向に進行する超音波ビームが送信される。
　ステップＳ２１１０の後、温度計測装置１０は、図１２と同様に、ステップＳ１１１１～Ｓ１１１６を実行する。
　第２実施形態によれば、音波送信装置２０の送信角が可変である。これにより、１つの音波送信装置２０から送信される超音波ビームの経路が増える。その結果、対象空間ＳＰの温度分布の計測に必要な音波送信装置２０の数を低減することができ、且つ、音波送信装置２０及び音波受信装置３０の配置の自由度を向上させることができる。
（３）変形例
　本実施形態の変形例について説明する。
（３－１）変形例１
　変形例１について説明する。変形例１は、音波送信装置２０及び音波受信装置３０に関する変形例である。図２０は、変形例１の音波送信装置及び音波受信装置の構成を示す概略図である。
　図２０に示すように、変形例１の対象空間ＳＰには、センサユニットＳＵが配置される。
　センサユニットＳＵは、音波送信装置２０と、音波受信装置３０と、を備える。
　音波送信装置２０は、複数の超音波振動子２１を備える。
　複数の超音波振動子２１は、一定方向（Ｘ方向）に配列される。
　音波受信装置３０は、複数の超音波振動子３１を備える。
　複数の超音波振動子３１は、一定方向（Ｙ方向）に配列される。
　音波送信装置２０の超音波振動子２１の配列方向と、音波受信装置３０の超音波振動子３１の配列方向と、は互いに異なる。好ましくは、超音波振動子２１の配列方向（Ｘ方向）は、超音波振動子３１の配列方向（Ｙ方向）と直交する。
　音波送信装置２０から送信された超音波ビームの音圧分布は、超音波振動子２１の配列方向（Ｘ方向）に対する直交方向（Ｙ方向）に広がっている。
　この超音波ビームは、対象空間ＳＰの反射部材で反射する。
　音波受信装置３０は、反射部材で反射した超音波ビームを受信する。超音波振動子３１の超音波振動子３１の超音波に対する感度分布は、超音波振動子３１の配列方向（Ｙ方向）に対する直交方向（Ｘ方向）に広がっている。
　変形例１によれば、Ｙ方向に広がる音圧分布を有する超音波ビームを、Ｘ方向に広がる感度分布を有する超音波振動子３１の振動に応じた波形データが得られる。つまり、音波受信装置３０は、音圧分布と感度分布が交わる点の成分を最も強く取得する。これにより、空間温度の計測の解像度を増加させることができる。
（３－２）変形例２
　変形例２について説明する。変形例２は、アクチュエータ２３を用いずに、送信方向を変更する例である。
　変形例２の制御回路２２は、超音波制御信号に含まれる送信角の値に応じて、Ｚ軸に対する送信角を実現するために必要な位相差を計算する。
　制御回路２２は、超音波ビームが計算した位相差が発生するように、複数の超音波振動子２１をそれぞれ異なるタイミングで振動させる。各超音波振動子２１の振動のタイミングの差は、各超音波振動子２１から送信される超音波の位相差を形成する。
　これにより、超音波制御信号に含まれる送信角の値が示す方向に進行する超音波ビームが送信される。
　変形例２によれば、アクチュエータ２３（つまり、送信方向を変えるための機械的機構）を用いることなく、送信方向を変えることができる。
（３－３）変形例３
　変形例３について説明する。変形例３は、対象空間ＳＰの温度に応じて、フィルタを補正する例である。
　変形例３のプロセッサ１２は、ステップＳ１１１３において、温度計５０から基準温度を取得する。
　プロセッサ１２は、基準温度に基づいて、特定したフィルタ情報を補正する。
　プロセッサ１２は、補正したフィルタ情報に基づいて、受信波形データに含まれる成分のうち、対象経路Ｐｉに沿って進行した超音波ビームの成分を抽出する。
　変形例３によれば、対象空間ＳＰの温度に基づいてフィルタリングを実行する。これにより、受信波形データから、対象経路Ｐｉに沿って進行した超音波ビームの成分をより確実に抽出することができる。その結果、メッシュ温度の計算の精度を向上させることができる。
（３－４）変形例４
　変形例４について説明する。変形例４は、対象空間ＳＰ内の気流が計測結果に与える影響を軽減する例である。図２１は、変形例４の音波送信装置及び音波受信装置の構成を示す概略図である。図２２は、変化例４の原理の説明図である。
　図２１に示すように、変形例４の対象空間ＳＰには、少なくとも２つのセンサユニットＳＵａ及びＳＵｂが配置される（図２１Ａ）。
　センサユニットＳＵａ及びＳＵｂは、それぞれ、音波送信装置２０と、音波受信装置３０と、を備える（図２１Ｂ）。
　センサユニットＳＵａの音波送信装置２０ａから送信された超音波ビームは、センサユニットＳＵｂの音波受信装置３０ｂによって受信される。
　音波受信装置３０ｂは、受信した超音波ビームに応じた受信波形データ（「第２受信波形データ」の一例）を生成する。
　センサユニットＳＵｂの音波送信装置２０ｂから送信された超音波ビームは、センサユニットＳＵａの音波受信装置３０ａによって受信される。
　音波受信装置３０ａは、受信した超音波ビームに応じた受信波形データ（「第１受信波形データ」の一例）を生成する。
　図２２に示すように、音波送信装置２０ａから音波受信装置３０ｂに向かう経路（以下「往路」という）を進行する超音波ビームの実速度の絶対値｜Ｖａｂ｜については、式２ａの関係が成立する。
　｜Ｖａｂ｜　＝　Ｄａｂ／ｔａｂ　＝　Ｃ＋Ｖｗａｂ　…（式２ａ）
　・Ｄａｂ：センサユニットＳＵａ及びＳＵｂの間のセンサ間距離
　・ｔａｂ：往路の伝搬時間
　・Ｖｗａｂ：センサユニットＳＵａ及びＳＵｂの間の風速成分
　音波送信装置２０ｂから音波受信装置３０ａに向かう経路（以下「復路」という）は、往路の逆方向である。したがって、復路を進行する超音波ビームの実速度の絶対値｜Ｖｂａ｜については、式２ｂの関係が成立する。
　｜Ｖｂａ｜　＝　Ｄａｂ／ｔｂａ　＝　Ｃ－Ｖｗａｂ　…（式２ｂ）
　・ｔｂａ：復路の伝搬時間
　温度計測装置１０は、式２ｃを用いて、センサユニットＳＵａ及びＳＵｂの平均速度｜Ｖａ：ｂ｜を計算する。
　｜Ｖａ：ｂ｜　＝　（｜Ｖａｂ｜＋｜Ｖｂａ｜）／２　…（式２ｃ）
　式２ｃにより、往路の風速成分Ｖａｂと復路の風速成分Ｖｂａが互いに相殺される。したがって、平均速度｜Ｖａ：ｂ｜は、風速成分Ｖｗａｂを含まない。
　温度計測装置１０は、記憶装置１１に記憶された音波速度特性に関するデータを参照して、平均速度｜Ｖａ：ｂ｜に対応する温度を計算する。
　変形例４によれば、１対のセンサユニットＳＵａ及びＳＵｂの間の風速成分を含まない平均速度に対応する温度を計算する。これにより、空間の温度の計測結果のＳ／Ｎ比をさらに向上させることができる。
　変形例４では、２つのセンサユニットＳＵａ及びＳＵｂを使用する例を説明したが、本実施形態はこれに限られない。本実施形態は、例えば、音波がほぼ同一の進行経路を通っており、且つ、それらの音波に含まれる風速成分のベクトル和がキャンセルされる条件を満たす限り、３種類以上の音速又はそれに準ずる信号の平均を用いる場合にも適用可能である。
　変形例４では、２つのセンサユニットＳＵａ及びＳＵｂを使用する例を説明したが、本実施形態はこれに限られない。変形例４は、1回以上の反射散乱経路を使用する場合にも適用可能である。
（３－５）変形例５
　変形例５について説明する。変形例５は、時系列フィルタを用いた温度計測アルゴリズムの例である。
（３－５－１）変形例５の概要
　変形例５の概要について説明する。図２３は、変形例５の概要の説明図である。
　図２３に示すように、変形例５のプロセッサ１２は、経路温度計算モデルＭｐｔ（ｔ）と、時系列フィルタＦＩＬと、を実行するように構成される。
　経路温度計算モデルＭｐｔ（ｔ）は、時刻ｔの受信波形データＲＷ（ｔ｜ｘ，ｙ，ｚ）に応じて、時刻ｔの経路温度ＰＤ（ｔ｜ｘ，ｙ，ｚ）を出力するように構成される。
　時系列フィルタＦＩＬは、経路温度計算モデルＭｐｔ（ｔ）の出力（経路温度ＰＤ（ｔ｜ｘ，ｙ，ｚ））、温度計５０によって計測された時刻ｔの基準温度Ｔｒｅｆ（ｔ）、及び、時刻ｔ－１の温度分布Ｄ（ｔ－１）の組合せに応じて、時刻ｔの温度分布Ｄｔ（ｔ）を出力するように構成される。
　時系列フィルタＦＩＬは、例えば、以下の少なくとも１つを含む。
　・カルマンフィルタ
　・拡張カルマンフィルタ
　・無香カルマンフィルタ
　・パーティクルフィルタ
（３－５－２）温度計測の処理
　変形例５の温度計測の処理について説明する。図２４は、変形例５の温度計測の処理のフローチャートである。
　図２４に示すように、変形例５の温度計測装置１０は、図１２と同様に、ステップＳ１１０～Ｓ１１３を実行する。
　全ての対象メッシュＭｔについてステップＳ１１２が終了していない場合（Ｓ１１３－ＮＯ）、温度計測装置１０は、ステップＳ１１０を実行する。
　全ての対象メッシュＭｔについてステップＳ１１２が終了すると（Ｓ１１３－ＹＥＳ）、温度計測装置１０は、時系列フィルタリング（Ｓ３１０）を実行する。
　具体的には、プロセッサ１２は、温度計５０から時刻ｔの基準温度Ｔｒｅｆ（ｔ）を取得する。
　プロセッサ１２は、ステップＳ１１１で得られた時刻ｔの経路温度Ｔｐ（ｔ｜ｘ，ｙ，ｚ）と、基準温度Ｔｒｅｆ（ｔ）と、を時刻ｔ－１の温度分布（ｔ－１）と、を時系列フィルタＦＩＬに入力することにより、時刻ｔの温度分布Ｄ（ｔ）を計算する。
　温度分布Ｄ（ｔ）は、時刻ｔ＋１の温度分布Ｄ（ｔ＋１）の計算において参照される。
　変形例５によれば、時系列フィルタリングを実行することにより、空間の温度の計測結果のＳ／Ｎ比を更に向上させることができる。
　なお、変形例５の時系列フィルタＦＩＬは、更に、時刻ｔ－１の外部環境情報を参照して、時刻ｔの温度分布Ｄ（ｔ）を計算しても良い。時刻ｔ－１の外部環境情報は、例えば、以下の情報を含む。
　・空調装置４０の熱量に関する情報
　・対象空間ＳＰの周辺の外気温に関する情報
　・対象空間ＳＰの３次元形状に関する情報
　・対象空間ＳＰの断熱性能に関する情報
　・対象空間ＳＰ内に存在する人の数に関する情報
　・対象空間ＳＰ内に存在する人の動きに関する情報
　・空調装置４０の風に関する情報
　・対象空間ＳＰ内の風に関する情報
（３－６）変形例６
　変形例６について説明する。変形例６は、Ｎ（Ｎは２以上の整数）個の音波送信装置２０、及び、Ｍ（Ｍは２以上の整数）個の音波受信装置３０の組合せの例である。図２５は、変形例６の音波送信装置及び音波受信装置の概略図である。
　図２５に示すように、Ｎ個の音波送信装置２０（１）～２０（Ｎ）が超音波ビームを送信する。
　Ｍ個の音波受信装置３０（１）～３０（Ｍ）は、それぞれ、Ｎ個の音波送信装置２０（１）～２０（Ｎ）から送信された超音波ビームを受信し、且つ、当該超音波ビームに対応する受信波形データを生成する。
　変形例６の第１例のプロセッサ１２は、音波受信装置３０（１）～３０（Ｍ）によって生成された受信波形データに対して、図１３のフィルタリング（Ｓ１１１３）を実行することにより、受信波形データに対応する超音波ビームの経路（つまり、当該超音波ビームの発信源である音波送信装置２０（ｎ）（ｎは、２～Ｎの整数）及び音波受信装置３０（ｍ）の組合せ）を特定する。
　変形例６の第２例のプロセッサ１２は、受信波形データを包絡線に変換し、且つ、包絡線の立ち上がりの時刻を特定する。プロセッサ１２は、当該時刻によって、受信波形データに対応する超音波ビームの経路（つまり、当該超音波ビームの発信源である音波送信装置２０（ｎ）及び音波受信装置３０（ｍ）の組合せ）を特定する。
　変形例６の第３例の音波送信装置２０（ｎ）は、互いに、周波数が異なる超音波ビームを送信する。プロセッサ１２は、受信波形データの周波数を参照して、受信波形データに対応する超音波ビームの経路（つまり、当該超音波ビームの発信源である音波送信装置２０（ｎ）及び音波受信装置３０（ｍ）の組合せ）を特定する。
　変形例６によれば、１つの音波受信装置３０が複数の音波送信装置２０（ｎ）から送信された超音波ビームを受信する。温度計測装置１０は、受信波形データをフィルタリングすることにより、超音波ビームを発振した音波送信装置２０（ｎ）を特定する。これにより、空間の温度の計測結果のＳ／Ｎ比を向上させることができる。
（３－７）変形例７
　変形例７について説明する。変形例７は、温度に加えて、風ベクトルの分布を計測する例である。
　変形例７の対象空間ＳＰには、例えば、図５の温度計測装置１０（不図示）と、少なくとも２つの音波送信装置２０ａ～２０ｂと、少なくとも２つの音波受信装置３０ａ～３０ｂと、が配置されている。
　温度計測装置１０は、所定の送信周波数Ｆｓを有する超音波ビームを送信するように、音波送信装置２０ａを制御する。
　音波送信装置２０ａは、温度計測装置１０の制御に従い、送信周波数Ｆｓを有する超音波ビームを送信する。
　音波受信装置３０ａは、超音波ビームを受信すると、受信波形データを生成する。音波受信装置３０ａが受信する超音波ビームには、音波送信装置２０ａと音波受信装置３０ａとの間の風に起因するドップラー効果が発生する。したがって、音波受信装置３０ａが受信する超音波ビームの受信周波数Ｆｒａは、送信周波数Ｆｓとは異なる。
　温度計測装置１０は、音波受信装置３０ａから受信波形データを取得し、且つ、受信波形データを参照して受信周波数Ｆｒａを特定する。
　ドップラー効果を考慮すると、送信周波数Ｆｓと受信周波数Ｆｒａとの間には、式３．１の関係が成立する。
　Ｆｓ　＝　（Ｃ＋Ｖｗａ）／Ｃ×Ｆｒａ　…（式３．１）
　・Ｃ：超音波の音速
　・Ｖｗａ：音波送信装置２０ａと音波受信装置３０ａとの間の経路上の理論上の風速
　式３．１を展開すると、理論上の風速Ｖｗａは、式３．２のように表すことができる。プロセッサ１２は、式３．１を用いて、理論上の風速Ｖｗａを計算する。
　Ｖｗａ　＝　Ｃ×Ｆｓ／Ｆｒａ－Ｃ　…（式３．２）
　プロセッサ１２は、式３．３を用いて、音波送信装置２０ｂと音波受信装置３０ｂとの間の経路上の理論上の風速Ｖｗｂを計算する。
　Ｖｗｂ　＝　Ｃ×Ｆｓ／Ｆｒｂ－Ｃ　…（式３．３）
　・Ｆｒｂ：音波受信装置３０ｂが受信した超音波ビームの受信周波数
　記憶装置１１には、温度と風速の相関関数Ｇ（ｘ）が予め格納されている。
　プロセッサ１２は、ステップＳ１１２で得られたメッシュ温度ＴＥＭＰｍｅｓｈ及び相関関数を用いて、式３．４～式３．５のように、補正風速Ｖｒｗａ及びＶｒｗｂを計算する。
　Ｖｒｗａ　＝　Ｇ（Ｔｍｅｓｈ）×Ｖｗａ　…（式３．４）
　Ｖｒｗｂ　＝　Ｇ（Ｔｍｅｓｈ）×Ｖｗｂ　…（式３．５）
　温度計測装置１０は、音波送信装置２０ａ～音波受信装置３０ａの進行経路及び音波送信装置２０ｂ～音波受信装置３０ｂの進行経路の成す角度を参照して、式３．４～式３．５から得られた補正風速Ｖｒｗａ及びＶｒｗｂを合成することにより、風ベクトルを計算する。
　変形例７によれば、超音波の進行経路が交差する交差点を含むメッシュの風ベクトルが得られる。したがって、複数の交差点が形成される場合、対象空間ＳＰの風ベクトルの分布（つまり、風速及び風向きの分布）が得られる。
　変形例７では、２本の進行経路が交差点を形成する例を示した。この場合、得られる風ベクトルの次元は２次元である。
　なお、変形例７は、３本の進行経路が交差点を形成する場合にも適用可能である。この場合、得られる風ベクトルの次元は３次元である。
（３－８）変形例８
　変形例８について説明する。変形例８は、超音波を用いて風ベクトルの分布を計測する例である。図２６は、変形例８のセンサ配置の一例を示す図である。
　図２６に示すように、変形例８の対象空間ＳＰには、少なくとも４つのセンサユニットＳＵａ～ＳＵｄが、センサユニットＳＵａ～ＳＵｂの超音波の進行経路と、センサユニットＳＵｃ～ＳＵｄの超音波の進行経路とが互いに交わる（好ましくは、直交する）ように配置される。
　各センサユニットＳＵａ～ＳＵｄの構成は、変形例４（図２１Ｂ）と同様である。
　変形例８の温度計測装置１０は、式４．１を用いて、センサユニットＳＵａ～ＳＵｂ間の往路を進行する超音波ビームの速度の絶対値｜Ｖａｂ｜から温度因子を除去することにより、センサユニットＳＵａ及びＳＵｂの間の風速成分Ｖｗａｂを計算する。
　Ｖｗａｂ　＝　｜Ｖａｂ｜－｜Ｖａ：ｂ｜　…（式４．１）
　温度計測装置１０は、式４．２を用いて、センサユニットＳＵｃ～ＳＵｄ間の往路を進行する超音波ビームの速度の絶対値｜Ｖｃｄ｜から温度因子を除去することにより、センサユニットＳＵｃ及びＳＵｄの間の風速成分Ｖｗｃｄを計算する。
　Ｖｗｃｄ　＝　｜Ｖｃｄ｜－｜Ｖｃ：ｄ｜　…（式４．２）
　・｜Ｖｃ：ｄ｜：センサユニットＳＵｃ及びＳＵｄの平均速度
　温度計測装置１０は、センサユニットＳＵａ～ＳＵｂの進行経路及びセンサユニットＳＵｃ～ＳＵｄの進行経路の成す角度を参照して、式４．１から得られた風速成分Ｖｗａｂと、式４．２から得られた風速成分Ｖｗｃｄと、を合成することにより、風ベクトルを計算する。
　変形例８によれば、対象空間ＳＰに配置されたセンサユニットの超音波の進行経路が交差点を含むメッシュの風ベクトルが得られる。したがって、複数の交差点が形成される場合、対象空間ＳＰの温度分布に代えて、風ベクトル（つまり、風速及び風向き）の分布が得られる。
　変形例８では、往路を進行する超音波ビームの絶対値｜Ｖａｂ｜及び｜Ｖｃｄ｜に代えて、復路を進行する超音波ビームの絶対値｜Ｖｂａ｜及び｜Ｖｄｃ｜を用いて、風ベクトルを推定しても良い。
　変形例８では、一対のセンサユニットが交差点を形成する例を示した。この場合、得られる風ベクトルの次元は２次元である。
　なお、変形例８は、３つのセンサユニットが交差点を形成する場合にも適用可能である。この場合、得られる風ベクトルの次元は３次元である。
（３－９）変形例９
　変形例９を説明する。変形例９は、振動子アレイを用いて温度を計測する例である。
（３－９－１）変形例９の計測システムの構成
　変形例９の計測システム１の構成を説明する。
（３－９－１－１）変形例９の音波送信装置の構成
　変形例９の音波送信装置２０の構成を説明する。図２７は、変形例９の音波送信装置の構成を示す図である。
　図２７に示すように、音波送信装置２０は、複数の超音波振動子２１と、制御回路２２と、を備える。
　図２７Ａに示すように、複数の超音波振動子２１は、送信面（ＸＹ平面）上に二次元配列される。つまり、複数の超音波振動子２１は、振動子アレイＴＡを形成する。
　図２７Ｂに示すように、各超音波振動子２１は、Ｚ方向に沿って進行する超音波ビームＵＳＷを送信する。
（３－９－１－２）変形例９の音波受信装置の構成
　変形例９の音波受信装置３０の構成を説明する。図２８は、変形例９の音波受信装置の構成を示す図である。
　図２８に示すように、音波受信装置３０は、複数の超音波振動子３１と、制御回路３２と、を備える。
　図２８Ａに示すように、複数の超音波振動子３１は、送信面（ＸＹ平面）上に二次元配列される。つまり、複数の超音波振動子３１は、振動子アレイＴＡを形成する。
　図２８Ｂに示すように、各超音波振動子３１は、音波送信装置２０から送信された超音波ビームＵＳＷを受信すると振動する。
（３－９－２）変形例９の具体例
　変形例９の具体例を説明する。
（３－９－２－１）変形例９の第１例
　変形例９の第１例を説明する。変形例９の第１例は、一対の音波受信装置３０及び音波送信装置２０を用いて温度を計測する例である。図２９は、変形例９の第１例の概要を示す図である。
　以下、対象空間ＳＰには、図２９のＸ＋方向に向かう気流ＡＦが存在する。
　変形例９の第１例の音波送信装置２０の振動子２１ａ～２１ｃは、それぞれ、Ｚ＋方向に向かって超音波ビームＵＳＷ０～ＵＳＷ２を送信する。
　超音波ビームＵＳＷ０～ＵＳＷ２は、気流ＡＦの影響により、Ｘ＋方向にシフトする。
　その結果、超音波ビームＵＳＷ２は音波受信装置３０の外側に向かって進行する。
　一方、超音波ビームＵＳＷ０～ＵＳＷ１は、超音波振動子３１ｂ～３１ｃに受信される。
　このように、音波送信装置２０及び音波受信装置３０が、何れも、振動子アレイＴＡを形成する。これにより、音波送信装置２０の振動子アレイＴＡから放射された超音波ビームＵＳＷ０～ＵＳＷ２は、気流ＡＦの影響を受けたとしても、音波受信装置３０の振動子アレイＴＡに到達し易くなる。これにより、気流ＡＦの影響を受けることなく、本実施形態と同様の効果が得られる。
（３－９－２－２）変形例９の第２例
　変形例９の第２例を説明する。変形例９の第２例は、一対のセンサユニットＳＵ（音波受信装置３０及び音波送信装置２０の組合せ）を用いて温度を計測する例である。図３０は、変形例９のセンサユニットの構成を示す図である。図３１は、変形例９の第２例の概要を示す図である。
　図３０Ａに示すように、変形例９の第２例のセンサユニットＳＵは、音波送信装置２０と、音波受信装置３０と、を備える。
　音波送信装置２０は、複数の超音波振動子２１と、制御回路２２と、を備える。
　複数の超音波振動子２１は、送信面（ＸＹ平面）上に二次元配列される。つまり、複数の超音波振動子２１は、振動子アレイＴＡを形成する。
　音波受信装置３０は、複数の超音波振動子３１と、制御回路３２と、を備える。
　複数の超音波振動子３１は、送信面（ＸＹ平面）上に二次元配列される。つまり、複数の超音波振動子３１は、振動子アレイＴＡを形成する。
　図３０Ｂに示すように、対象空間ＳＰには、一対のセンサユニットＳＵａ及びＳＵｂが配置される。
　センサユニットＳＵａの音波送信装置２０は、センサユニットＳＵｂの音波受信装置３０に向かって超音波を放射する。
　センサユニットＳＵｂの音波送信装置２０は、センサユニットＳＵａの音波受信装置３０に向かって超音波を放射する。
　以下、対象空間ＳＰには、図３１のＸ＋方向に向かう気流ＡＦが存在する。
　センサユニットＳＵａの音波送信装置２０の振動子２１ａ～２１ｃは、それぞれ、Ｚ＋方向に向かって超音波ビームＵＳＷａ０～ＵＳＷａ２を送信する。
　センサユニットＳＵｂの音波送信装置２０の振動子２１ａ～２１ｃは、それぞれ、Ｚ－方向に向かって超音波ビームＵＳＷｂ０～ＵＳＷｂ２を送信する。
　超音波ビームＵＳＷａ０～ＵＳＷａ２及びＵＳＷｂ０～ＵＳＷｂ２は、気流ＡＦの影響により、Ｘ＋方向にシフトする。
　その結果、超音波ビームＵＳＷａ２は、センサユニットＳＵｂの音波受信装置３０の外側に向かって進行する。超音波ビームＵＳＷｂ２は、センサユニットＳＵａの音波受信装置３０の外側に向かって進行する。
　一方、超音波ビームＵＳＷａ０～ＵＳＷａ１は、センサユニットＳＵｂの音波受信装置３０の超音波振動子３１ｂ～３１ｃに受信される。超音波ビームＵＳＷｂ０～ＵＳＷｂ１は、センサユニットＳＵａの音波受信装置３０の超音波振動子３１ｂ～３１ｃに受信される。
　このように、音波送信装置２０及び音波受信装置３０が、何れも、振動子アレイＴＡを形成する。これにより、音波送信装置２０の振動子アレイＴＡから放射された超音波ビームＵＳＷ０～ＵＳＷ２は、気流ＡＦの影響を受けたとしても、音波受信装置３０の振動子アレイＴＡに到達し易くなる。更に、センサユニットＳＵａ～ＳＵｂの間の空間の温度を、Ｚ－方向に進行する超音波ビーム及びＺ＋方向に進行する超音波ビームの両方を用いて計測する。これにより、気流ＡＦの影響を受けることなく、本実施形態と同様の効果が得られる。
（３－１０）変形例１０
　変形例１０を説明する。変形例１０は、１つの音波受信装置３０が受信した音波の波形の中から、参照すべき波形を特定する例である。図３２は、変形例１０のフィルタの説明図である。
（３－１０－１）変形例１０の第１例
　変形例１０の第１例を説明する。
　変形例１の第１例の音波受信装置３０は、経路毎の時間フィルタを有する。時間フィルタは、例えば、以下の少なくとも１つである（図３２）。
　・下限時間閾値ＴＨｔｂ
　・上限時間閾値ＴＨｔｔ
　・下限時間閾値ＴＨｔｂと上限時間閾値ＴＨｔｔとによって規定される時間ウインドウＷｔ
　下限時間閾値ＴＨｔｂ及び上限時間閾値ＴＨｔｔは、以下の少なくとも１つによって決まる。
　・音波送信装置２０と音波受信装置３０との間の距離
　・音波送信装置２０から送信された超音波ビームが音波受信装置３０によって受信されるまでの間に超音波ビームが伝搬する距離
　ステップＳ１１１２（図１３）において、制御回路３２は、超音波振動子３１の振動に応じた受信波形データ（図３２）を生成した後、時間フィルタ（例えば、時間ウインドウＷｔ）を受信波形データに適用することにより、受信波形データに含まれる波形ＷＦ１～ＷＦ３の中から参照すべき波形ＷＦ２を抽出する。
　制御回路３２は、抽出された波形ＷＦ２を含む受信波形データを温度計測装置１０に送信する。
（３－１０－１）変形例１０の第２例
　変形例１０の第２例を説明する。
　変形例１０の第２例の温度計測装置１０は、第１例と同様の時間フィルタを有する。
　ステップＳ１１１３（図１３）において、プロセッサ１２は、時間フィルタ（例えば、時間ウインドウＷｔ）を受信波形データに適用することにより、受信波形データに含まれる波形ＷＦ１～ＷＦ３の中から参照すべき波形ＷＦ２を抽出する。
（４）その他の変形例
　その他の変形例を説明する。
　記憶装置１１は、ネットワークＮＷを介して、温度計測装置１０と接続されてもよい。
　図４の例では、超音波振動子３１を備える音波受信装置３０の例を示した。しかし、音波受信装置３０は、音波送信装置２０と同様に、複数の超音波振動子３１を備えても良い。
　図５の例では、１個の音波送信装置２０が複数の経路に沿った超音波ビームを送信し、且つ、１個の音波受信装置３０が複数の経路に沿った超音波ビームを受信する例を示した。しかし、本実施形態はこれに限られない。ｎ（ｎは２以上の整数）個の音波送信装置２０のそれぞれが１本の経路に沿った超音波ビーム（つまり、ｎ個の音波送信装置２０がｎ本の経路に沿った超音波ビーム）を送信し、且つ、ｎ個の音波受信装置３０のそれぞれが各経路に沿った超音波ビームを受信しても良い（つまり、ｎ個の音波受信装置３０がｎ本の経路に沿った超音波ビームを受信しても良い）。
　上記の実施形態では、メッシュ温度ＴＥＭＰｍｅｓｈｔの計算に平均値を求める関数を用いる例を示したが、本実施形態のメッシュ温度ＴＥＭＰｍｅｓｈｔの計算方法はこれに限られるものではない。
　音波送信装置２０は、自己相関が比較的強い自己相関信号（例えば、Ｍ系列信号）を含む超音波ビームを送信しても良い。これにより、空間の温度の計測結果のＳ／Ｎ比を更に向上させることができる。
　音波送信装置２０が個別に異なる自己相関信号を含む超音波ビームを送信することにより、音波受信装置３０が、超音波ビームの発信源となる音波送信装置２０を識別しても良い。
　また、音波送信装置２０毎に異なる発振周波数を有する超音波ビームを送信することにより、音波受信装置３０が、超音波ビームの発信源となる音波送信装置２０を識別しても良い。
　温度計測装置１０は、温度分布及び風ベクトル分布以外に、以下の空気特性の分布を計測することも可能である。
　・空気中の化学物質（例えば、ＣＯ２）の濃度の分布
　・湿度の分布
　・臭気の分布
　・有毒ガスの分布
　本実施形態では、音波送信装置２０及び音波受信装置３０を区別して規定したが、本実施形態の範囲は、これに限られない。本実施形態は、１つの超音波振動子が超音波を送信する機能及び超音波を受信する機能を備えても良い。
　本実施形態では、ステップＳ１１１４（図１３）において経路温度ＴＥＭＰｐａｔｈｉの計算に用いる式、及び、ステップＳ１１２（図１２）においてメッシュ温度ＴＥＭＰｍｅｓｈｔの計算に用いる式の少なくとも１つは、外部環境情報（例えば、外気温、外気の湿度、及び、外気圧の少なくとも１つ）をパラメータとして含んでも良い。この場合、外部環境情報に関わらず、空間の空気特性の計測結果のＳ／Ｎ比を向上させることができる。
　本実施形態では、音波送信装置２０は、指向性を有する超音波ビームを送信する例を示したが、本実施形態は、これに限られない。本実施形態は、音波送信装置２０が可聴音ビーム（つまり、超音波ビームとは異なる周波数を有する音波）を送信する場合にも適用可能である。
　本実施形態において、温度分布とは、メッシュ温度ＴＥＭＰｍｅｓｈに限られない。温度分布は、以下の少なくとも１つも含む。
　・経路上の複数点の温度
　・経路上の平均温度
　以上、本発明の実施形態について詳細に説明したが、本発明の範囲は上記の実施形態に限定されない。また、上記の実施形態は、本発明の主旨を逸脱しない範囲において、種々の改良や変更が可能である。また、上記の実施形態及び変形例は、組合せ可能である。
１　　　　　：計測システム
１０　　　　：温度計測装置
１１　　　　：記憶装置
１２　　　　：プロセッサ
１３　　　　：入出力インタフェース
１４　　　　：通信インタフェース
２０　　　　：音波送信装置
２１　　　　：超音波振動子
２２　　　　：制御回路
２３　　　　：アクチュエータ
３０　　　　：音波受信装置
３１　　　　：超音波振動子
３２　　　　：制御回路
４０　　　　：空調装置
５０　　　　：温度計

Claims (20)

  1. 　対象空間に配置された音波送信装置、及び、前記対象空間に配置された音波受信装置と接続可能な温度計測装置であって、
    　音波を送信させるように、前記音波送信装置を制御する手段を備え、
    　前記音波を受信する音波受信装置から、前記音波の波形に関する受信波形データを取得する手段を備え、
    　前記受信波形データに基づいて、前記対象空間の温度分布を計算する手段を備える、
    温度計測装置。
  2. 　対象空間に配置された少なくとも２つのセンサユニットと接続可能な温度計測装置であって、
    　　第１センサユニットは、第１音波送信装置、及び、第１音波受信装置を備え、
    　　第２センサユニットは、第２音波送信装置、及び、第２音波受信装置を備え、
    　音波を送信させるように、各音波送信装置を制御する手段を備え、
    　前記第１音波受信装置から、前記第２音波送信装置から送信された音波の波形に関する第１受信波形データを取得する手段を備え、
    　前記第２音波受信装置から、前記第１音波送信装置から送信された音波の波形に関する第２受信波形データを取得する手段を備え、
    　前記第１受信波形データ及び前記第２受信波形データに基づいて、前記第１センサユニット及び前記第２センサユニットの間の対象空間の温度分布を計算する手段を備える、
    温度計測装置。
  3. 　前記計算する手段は、
    　　前記第１受信波形データに基づいて、前記第１音波送信装置によって送信された音波の第１実速度を計算し、
    　　前記第２受信波形データに基づいて、前記第２音波送信装置によって送信された音波の第２実速度を計算し、
    　　前記第１実速度及び前記第２実速度の平均速度を計算し、
    　　前記平均速度に基づいて、前記温度分布を計算する、
    請求項２に記載の温度計測装置。
  4. 　対象空間に配置されたＮ（Ｎは２以上の整数）個の音波送信装置及びＭ（Ｍは２以上の整数）個の音波受信装置と接続可能な温度計測装置であって、
    　音波を送信させるように、各音波送信装置を制御する手段を備え、
    　各音波受信装置から、前記音波の波形に関する受信波形データを取得する手段を備え、
    　Ｍ個の受信波形データに基づいて、前記受信波形データに対応する音波の発信源である音波送信装置及び前記音波受信装置の組合せを特定する手段を備え、
    　前記特定された組合せ毎の受信波形データに基づいて、前記対象空間の温度分布を計算する手段を備える、
    温度計測装置。
  5. 　前記対象空間を複数の仮想メッシュに分割する手段を備え、
    　各仮想メッシュについて、各仮想メッシュに含まれる複数の経路を決定する手段を備え、
    　前記決定された複数の経路に沿って進行する音波を送信させるように、前記音波送信装置を制御する手段を備え、
    　前記決定された複数の経路のそれぞれについて、前記受信波形データに基づいて、前記音波の伝搬時間を計算する手段を備え、
    　各仮想メッシュを通る複数の経路の伝搬時間に基づいて、各仮想メッシュのメッシュ温度を計算する手段を備え、
    請求項１～４の何れかに記載の温度計測装置。
  6. 　前記経路を決定する手段は、前記対象空間に存在する反射部材の反射特性に関する情報に基づいて、前記経路を決定する、
    請求項５に記載の温度計測装置。
  7. 　前記反射特性は、前記反射部材の反射種別、反射率、及び、反射面の法線角度の少なくとも１つを含む、
    請求項６に記載の温度計測装置。
  8. 　前記制御する手段は、前記音波送信装置の送信面の向きを変更することにより、前記音波の送信角を変更する、
    請求項１～７の何れかに記載の温度計測装置。
  9. 　前記制御する手段は、前記音波送信装置に含まれる複数の振動素子のそれぞれを異なるタイミングで振動させることにより、前記音波の送信角を変更する、
    請求項１～８の何れかに記載の温度計測装置。
  10. 　前記受信波形データに対してフィルタを適用することにより、前記経路に沿って進行する音波の成分を抽出する手段を備え、
    　前記計算する手段は、前記抽出された成分に基づいて、前記温度分布を計算する、
    請求項１～９の何れかに記載の温度計測装置。
  11. 　前記対象空間の温度に基づいて前記フィルタを補正する手段を備える、
    請求項１０に記載の温度計測装置。
  12. 　前記計算する手段は、時系列フィルタを用いて、前記温度分布を計算する、
    請求項１～１１の何れかに記載の温度計測装置。
  13. 　前記音波送信装置から送信された音波の送信周波数と、前記音波受信装置により受信された音波の受信周波数と、前記温度分布と、を参照して、前記対象空間の風ベクトルの分布を計算する手段を備える、
    請求項１～１２のいずれかに記載の温度計測装置。
  14. 　複数の受信波形データに時間フィルタを適用することにより、参照すべき波形を抽出する手段を備え、
    　前記計算する手段は、前記抽出された波形に基づいて、前記分布を計算する、
    請求項１～１３の何れかに記載の温度計測装置。
  15. 　対象空間に配置された少なくとも２つのセンサユニットと接続可能な温度計測装置であって、
    　　第１センサユニットは、第１音波送信装置、及び、第１音波受信装置を備え、
    　　第２センサユニットは、第２音波送信装置、及び、第２音波受信装置を備え、
    　音波を送信させるように、各音波送信装置を制御する手段を備え、
    　前記第１音波受信装置から、前記第２音波送信装置から送信された音波の波形に関する第１受信波形データを取得する手段を備え、
    　前記第２音波受信装置から、前記第１音波送信装置から送信された音波の波形に関する第２受信波形データを取得する手段を備え、
    　前記第１受信波形データ及び前記第２受信波形データに基づいて、前記対象空間の空気特性の分布を計算する手段を備える、
    計測装置。
  16. 　前記計算する手段は、
    　　前記第１受信波形データに基づいて、前記第１音波送信装置によって送信された音波の第１実速度を計算し、
    　　前記第２受信波形データに基づいて、前記第２音波送信装置によって送信された音波の第２実速度を計算し、
    　　前記第１実速度及び前記第２実速度の平均速度を計算し、
    　　前記平均速度に基づいて、前記空気特性の分布を計算する、
    請求項１５に記載の計測装置。
  17. 　前記空気特性は、前記対象空間内を流れる風の風ベクトルである、
    請求項１５又は１６に記載の計測装置。
  18. 　前記空気特性は、前記対象空間の温度、並びに、前記対象空間内の空気中の化学物質の濃度、湿度、臭気、及び、有毒ガスの少なくとも１つである、
    請求項１５～１７の何れかに記載の計測装置。
  19. 　対象空間に配置される音波受信装置であって、
    　前記対象空間に配置された音波送信装置から送信された音波ビームを受信する手段を備え、
    　前記音波ビームの受信波形データを生成する手段を備え、
    　前記受信波形データに時間フィルタを適用することにより、前記受信波形データに含まれる複数の波形の中から特定の波形を抽出する手段を備え、
    　前記抽出された波形を含む受信波形データを温度計測装置に送信する手段を備える、
    音波受信装置。
  20. 　コンピュータに、請求項１～１８の何れかに記載の各手段を実現させるためのプログラム。
PCT/JP2019/033965 2018-11-30 2019-08-29 温度計測装置、計測装置、音波受信装置、及び、プログラム WO2020110393A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2020558099A JP6898022B2 (ja) 2018-11-30 2019-08-29 計測装置、音波受信装置、プログラム、及び、計測方法
CN201980078514.XA CN113167660A (zh) 2018-11-30 2019-08-29 温度测量装置、测量装置、声波接收装置以及程序
US17/229,120 US20210231507A1 (en) 2018-11-30 2021-04-13 Measuring apparatus, and measuring method

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2018-225630 2018-11-30
JP2018225630 2018-11-30
JP2019-022911 2019-02-12
JP2019022911 2019-02-12

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/229,120 Continuation US20210231507A1 (en) 2018-11-30 2021-04-13 Measuring apparatus, and measuring method

Publications (1)

Publication Number Publication Date
WO2020110393A1 true WO2020110393A1 (ja) 2020-06-04

Family

ID=70852767

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2019/033965 WO2020110393A1 (ja) 2018-11-30 2019-08-29 温度計測装置、計測装置、音波受信装置、及び、プログラム

Country Status (4)

Country Link
US (1) US20210231507A1 (ja)
JP (3) JP6898022B2 (ja)
CN (1) CN113167660A (ja)
WO (1) WO2020110393A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021199877A1 (ja) * 2020-04-01 2021-10-07 ピクシーダストテクノロジーズ株式会社 制御システム、制御方法、およびプログラム
WO2022091887A1 (ja) * 2020-10-29 2022-05-05 ダイキン工業株式会社 環境状態測定装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0280922A (ja) * 1988-09-16 1990-03-22 Babcock Hitachi Kk ガス流体の流速および温度を測定する方法
JPH05240719A (ja) * 1992-02-28 1993-09-17 Japan Radio Co Ltd 超音波遠隔水温測定装置
JPH07280666A (ja) * 1994-04-07 1995-10-27 Ishikawajima Harima Heavy Ind Co Ltd 炉内ガス温度測定装置
JP2006118921A (ja) * 2004-10-20 2006-05-11 Arata Kaneko 実時間型水中音響トモグラフィ装置及び水中音響トモグラフィシステム
JP2007222438A (ja) * 2006-02-24 2007-09-06 Aichi Tokei Denki Co Ltd ショーケース、温度維持ケース及び温度維持ケースの温度検出方法
JP2014102091A (ja) * 2012-11-16 2014-06-05 Panasonic Corp 超音波の伝搬時間計測のための超音波受信装置、及び超音波の伝搬時間計測のための超音波受信方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH628983A5 (fr) * 1978-09-19 1982-03-31 Pierre Emile Charles Ravussin Appareil pour la mesure automatique a distance du profil vertical de la temperature de l'atmosphere.
JP2758680B2 (ja) * 1989-12-26 1998-05-28 オーバル機器工業株式会社 超音波流量計
US5343744A (en) * 1992-03-06 1994-09-06 Tsi Incorporated Ultrasonic anemometer
JP3224286B2 (ja) * 1992-11-02 2001-10-29 株式会社日本自動車部品総合研究所 超音波を用いた温度測定装置
JP3424024B2 (ja) * 1994-10-17 2003-07-07 川崎重工業株式会社 ガス温度推定方法および装置
JP2000205971A (ja) * 1999-01-07 2000-07-28 Babcock Hitachi Kk 音響式ガス温度計測装置及びその方法
JP2002181598A (ja) 2000-12-18 2002-06-26 Yokogawa Electric Corp 超音波渦流量計
JP2003090770A (ja) * 2001-09-20 2003-03-28 Babcock Hitachi Kk 音波式ガス温度計測装置と方法
JP2003194604A (ja) 2001-12-26 2003-07-09 Babcock Hitachi Kk 音波式ガス温度・流量計とガス中音波の伝播時間の検出方法
JP2012026654A (ja) 2010-07-23 2012-02-09 Daikin Industries Ltd 空調機
CN102288324B (zh) * 2011-07-26 2013-12-18 沈阳工业大学 仓储粮食温度分布声学监测法
JP5875492B2 (ja) * 2012-09-13 2016-03-02 三菱電機株式会社 風計測装置
CN106771347B (zh) * 2016-12-30 2019-08-09 内蒙古工业大学 一种扫频式超声波风速测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0280922A (ja) * 1988-09-16 1990-03-22 Babcock Hitachi Kk ガス流体の流速および温度を測定する方法
JPH05240719A (ja) * 1992-02-28 1993-09-17 Japan Radio Co Ltd 超音波遠隔水温測定装置
JPH07280666A (ja) * 1994-04-07 1995-10-27 Ishikawajima Harima Heavy Ind Co Ltd 炉内ガス温度測定装置
JP2006118921A (ja) * 2004-10-20 2006-05-11 Arata Kaneko 実時間型水中音響トモグラフィ装置及び水中音響トモグラフィシステム
JP2007222438A (ja) * 2006-02-24 2007-09-06 Aichi Tokei Denki Co Ltd ショーケース、温度維持ケース及び温度維持ケースの温度検出方法
JP2014102091A (ja) * 2012-11-16 2014-06-05 Panasonic Corp 超音波の伝搬時間計測のための超音波受信装置、及び超音波の伝搬時間計測のための超音波受信方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021199877A1 (ja) * 2020-04-01 2021-10-07 ピクシーダストテクノロジーズ株式会社 制御システム、制御方法、およびプログラム
WO2022091887A1 (ja) * 2020-10-29 2022-05-05 ダイキン工業株式会社 環境状態測定装置
JP2022074026A (ja) * 2020-10-29 2022-05-17 ダイキン工業株式会社 環境状態測定装置
JP7152692B2 (ja) 2020-10-29 2022-10-13 ダイキン工業株式会社 環境状態測定装置

Also Published As

Publication number Publication date
JP7030368B2 (ja) 2022-03-07
JP6898022B2 (ja) 2021-07-07
US20210231507A1 (en) 2021-07-29
JP2021144049A (ja) 2021-09-24
CN113167660A (zh) 2021-07-23
JPWO2020110393A1 (ja) 2021-02-15
JP2022078086A (ja) 2022-05-24

Similar Documents

Publication Publication Date Title
Mao et al. CAT: High-precision acoustic motion tracking
JP7030368B2 (ja) 情報処理装置、情報処理方法、及び、プログラム
JP2019530049A (ja) 触覚振動子デバイスによる反響定位
JP2011179940A (ja) 超音波流量計
JP5321106B2 (ja) 超音波計測器
JP6697982B2 (ja) ロボットシステム
JP5683397B2 (ja) ポインティングシステム
JP7001289B2 (ja) 情報処理装置、情報処理方法、及び、プログラム
JP6043026B1 (ja) 流場測定装置および流場測定方法
WO2021145440A1 (ja) 計測装置、計測方法、及び、プログラム
JP2022124347A (ja) 計測装置、プログラム、および計測方法
JP2021076459A (ja) 音波受信装置、音波送信装置、および計測システム
WO2021199877A1 (ja) 制御システム、制御方法、およびプログラム
EP3654041B1 (en) Volume acceleration sensor calibration
JP2021113808A (ja) 計測装置、計測方法、及び、プログラム
JP2021165626A (ja) 計測装置、コントローラ、プログラム、および計測方法
JP2002054999A (ja) 音波発生検出装置、環境状態測定装置及び方法
JP2004257877A (ja) 音源検出方法、音源検出装置及びロボット
WO2022249882A1 (ja) 情報処理装置、情報処理方法、およびプログラム
JP2022110654A (ja) 計測装置、プログラム、および計測方法
JP2002267732A (ja) パッシブ位置標定方法およびパッシブ位置標定装置
JP2022125505A (ja) 計測装置、計測方法、およびプログラム
JP5647070B2 (ja) ポインティングシステム
JP5665656B2 (ja) ポインティングシステム
JP5462656B2 (ja) ポインティングシステム

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: 19891417

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2020558099

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19891417

Country of ref document: EP

Kind code of ref document: A1