JP2022059422A - Arithmetic method and arithmetic device - Google Patents

Arithmetic method and arithmetic device Download PDF

Info

Publication number
JP2022059422A
JP2022059422A JP2020167162A JP2020167162A JP2022059422A JP 2022059422 A JP2022059422 A JP 2022059422A JP 2020167162 A JP2020167162 A JP 2020167162A JP 2020167162 A JP2020167162 A JP 2020167162A JP 2022059422 A JP2022059422 A JP 2022059422A
Authority
JP
Japan
Prior art keywords
value
exercise
input
physical exercise
time point
Prior art date
Legal status (The legal status 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 status listed.)
Pending
Application number
JP2020167162A
Other languages
Japanese (ja)
Inventor
隆広 和田
Takahiro Wada
勇起 佐藤
Yuki Sato
勇希 岡藤
Yuki Okafuji
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nara Institute of Science and Technology NUC
Original Assignee
Nara Institute of Science and Technology NUC
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 Nara Institute of Science and Technology NUC filed Critical Nara Institute of Science and Technology NUC
Priority to JP2020167162A priority Critical patent/JP2022059422A/en
Publication of JP2022059422A publication Critical patent/JP2022059422A/en
Pending legal-status Critical Current

Links

Images

Abstract

To provide an arithmetic device capable of more accurately calculating a level of kinesia.SOLUTION: An arithmetic device for calculating an estimated value of kinesia from an input value comprises: a sensor I/F 14 that functions as an input unit for inputting a first body motion value showing motion stimulation given to a person's body at a first point of time; and an arithmetic unit 110 for calculating an estimated value of kinesia of the person. The arithmetic unit uses a history of body motion values showing motion stimulation given to the person's body at points of time before the first point of time to predict a second body motion value, a prediction value for a body motion value showing motion stimulation to be given to the person's body at the first point of time; takes a value based on the first body motion value and the second body motion value as the input value and inputs the value to an arithmetic model for calculating the estimated value of kinesia from the input value to calculate the estimated value of kinesia of the person at the first point of time.SELECTED DRAWING: Figure 1

Description

本開示は、演算方法、及び、演算装置に関する。 The present disclosure relates to an arithmetic method and an arithmetic unit.

動揺病は平衡失調に起因する疾病であり、車酔い、船酔いなどが含まれる。そのため、車両や船舶などの乗り物の快適性(乗り心地)を向上させるために、動揺病の発生や程度を抑えることは効果的である。 Motion sickness is a disease caused by imbalance and includes motion sickness and sea sickness. Therefore, in order to improve the comfort (ride quality) of vehicles such as vehicles and ships, it is effective to suppress the occurrence and degree of motion sickness.

この点について、本願発明者による特開2020-18770号公報(以下、特許文献1)は、前庭感覚(半規管、耳石)と視覚、及び、それらの統合過程を模した数理モデルから構築された、動揺病の原因といわれる「感覚矛盾説」に基づく数理モデルを開示している。この数理モデルを用いることで、動揺病の程度を算出し、乗り物の揺れの制御などを行うことが可能になる。 In this regard, Japanese Patent Application Laid-Open No. 2020-18770 (hereinafter referred to as Patent Document 1) by the inventor of the present application was constructed from a mathematical model that imitates the vestibular sensation (semicircular canals, otoliths) and vision, and their integration process. , Discloses a mathematical model based on the "sensory contradiction theory" that is said to be the cause of agitation. By using this mathematical model, it is possible to calculate the degree of motion sickness and control the shaking of the vehicle.

特開2020-18770号公報Japanese Unexamined Patent Publication No. 2020-18770

乗り物の快適性をより向上させるには、動揺病の程度を高精度に算出することが望まれる。そこで、発明者は、自身が先に見出した数理モデルに基づいて、より高精度で動揺病の程度を算出するよう開発を重ねた。 In order to further improve the comfort of the vehicle, it is desirable to calculate the degree of motion sickness with high accuracy. Therefore, the inventor has repeatedly developed to calculate the degree of motion sickness with higher accuracy based on the mathematical model he found earlier.

ある実施の形態に従うと、演算方法は、入力値から動揺病の推定値を算出する演算モデルを用いた演算方法において、第1の時点において人の身体に与えられた運動刺激を表す第1の身体運動値を取得し、第1の時点において人の身体に与えられる運動刺激を表す身体運動値の予測値である第2の身体運動値を、第1の時点以前の時点に人の身体に与えられた運動刺激を表す身体運動値の履歴を用いて予測し、第1の身体運動値と第2の身体運動値とに基づく値を入力値として演算モデルに入力し、人の第1の時点における動揺病の推定値を算出する、ことを含む。 According to a certain embodiment, the calculation method is a calculation method using a calculation model for calculating an estimated value of motion sickness from an input value, and is a first method representing a motor stimulus given to a human body at a first time point. The second physical exercise value, which is a predicted value of the physical exercise value representing the exercise stimulus given to the human body at the first time point, is applied to the human body at the time point before the first time point. Prediction is made using the history of body exercise values representing a given exercise stimulus, and a value based on the first body exercise value and the second body exercise value is input to the calculation model as an input value, and the first person's first body exercise value. Includes calculating estimates of motion sickness at a time point.

他の実施の形態に従うと、演算装置は、入力値から動揺病の推定値を算出する演算装置であって、第1の時点において人の身体に与えられた運動刺激を表す第1の身体運動値を入力する入力部と、動揺病の推定値を算出する演算部と、を備え、演算部は、第1の時点において人の身体に与えられる運動刺激を表す身体運動値の予測値である第2の身体運動値を、第1の時点以前の時点に人の身体に与えられた運動刺激を表す身体運動値の履歴を用いて予測し、第1の身体運動値と第2の身体運動値とに基づく値を入力値として、入力値から動揺病の推定値を算出する演算モデルに入力し、人の第1の時点における動揺病の推定値を算出する、ことを含む。 According to another embodiment, the arithmetic device is an arithmetic device that calculates an estimated value of motion sickness from an input value, and is a first physical exercise representing a motor stimulus given to a human body at a first time point. It includes an input unit for inputting a value and an arithmetic unit for calculating an estimated value of motion sickness, and the arithmetic unit is a predicted value of a physical exercise value representing an exercise stimulus given to a human body at a first time point. The second physical exercise value is predicted using the history of the physical exercise value representing the exercise stimulus given to the human body before the first time point, and the first physical exercise value and the second physical exercise are predicted. It includes inputting a value based on a value into an arithmetic model for calculating an estimated value of motion sickness from an input value and calculating an estimated value of motion sickness at a person's first time point.

更なる詳細は、後述の実施形態として説明される。 Further details will be described in the embodiments described below.

図1は、実施の形態に係る演算装置の装置構成の概略を示したブロック図である。FIG. 1 is a block diagram showing an outline of an apparatus configuration of an arithmetic unit according to an embodiment. 図2は、演算装置での演算処理の概要を表した図である。FIG. 2 is a diagram showing an outline of arithmetic processing in the arithmetic unit. 図3は、演算装置の演算部の機能構成の一例を表したブロック図である。FIG. 3 is a block diagram showing an example of the functional configuration of the arithmetic unit of the arithmetic unit. 図4は、図3の各機能の詳細な処理の一例を表したブロック図である。FIG. 4 is a block diagram showing an example of detailed processing of each function of FIG. 図5は、人に与えられた運動刺激の履歴の一例である。FIG. 5 is an example of the history of exercise stimuli given to a person. 図6は、図5の運動刺激の履歴から得られた予測値の信頼度である。FIG. 6 shows the reliability of the predicted value obtained from the history of the motor stimulus of FIG. 図7は、人に与えられた運動刺激の履歴の一例である。FIG. 7 is an example of the history of exercise stimuli given to a person. 図8は、図7の運動刺激の履歴から得られた予測値の信頼度である。FIG. 8 shows the reliability of the predicted value obtained from the history of the motor stimulus of FIG. 7. 図9は、図5の運動刺激の履歴から得られた予測値、及び図6の信頼度を用いた、推定値である。9 is an estimated value using the predicted value obtained from the history of the motor stimulus of FIG. 5 and the reliability of FIG. 図10は、図7の運動刺激の履歴から得られた予測値、及び図8の信頼度を用いた、推定値である。FIG. 10 is an estimated value using the predicted value obtained from the history of the motor stimulus of FIG. 7 and the reliability of FIG.

<1.演算方法、及び、演算装置の概要> <1. Overview of arithmetic method and arithmetic unit>

(1)実施の形態に係る演算方法は、入力値から動揺病の推定値を算出する演算モデルを用いた演算方法において、第1の時点において人の身体に与えられた運動刺激を表す第1の身体運動値を取得し、第1の時点において人の身体に与えられる運動刺激を表す身体運動値の予測値である第2の身体運動値を、第1の時点以前の時点に人の身体に与えられた運動刺激を表す身体運動値の履歴を用いて予測し、第1の身体運動値と第2の身体運動値とに基づく値を入力値として演算モデルに入力し、人の第1の時点における動揺病の推定値を算出すること、を含む。 (1) The calculation method according to the embodiment is a calculation method using a calculation model for calculating an estimated value of motion sickness from an input value, and is a first calculation method representing a motor stimulus given to a human body at a first time point. The second physical exercise value, which is a predicted value of the physical exercise value representing the exercise stimulus given to the human body at the first time point, is obtained from the body exercise value of the human body before the first time point. Prediction is made using the history of the physical exercise value representing the exercise stimulus given to, and the value based on the first physical exercise value and the second physical exercise value is input to the calculation model as the input value, and the first of the person Includes calculating estimates of motion sickness at the time of.

人の身体に与えられた運動刺激とは、例えば、頭部に与えられた運動刺激である。動揺病の推定値は、例えば、動揺病発症率(MSI:Motion Sickness Incidence)である。MSIは、物理刺激に対して嘔吐に至る人数の割合である。 The motor stimulus given to the human body is, for example, the motor stimulus given to the head. The estimated value of motion sickness is, for example, the incidence of motion sickness (MSI: Motion Sickness Incidence). MSI is the ratio of the number of people leading to vomiting to physical stimuli.

第1の時点は、例えば現時点である。この場合、第1の身体運動値は現時点で得られた身体運動値を指す。身体運動値の履歴は、第1の時点以前の時点に人の身体に与えられた運動刺激を表す、1又は複数の身体運動値である。第1の時点以前の時点の複数の身体運動値には第1の時点の身体運動値が含まれてもよい。第2の身体運動値は、過去の身体運動値から予測された現時点での身体運動値の予測値を指す。 The first time point is, for example, the present time. In this case, the first physical exercise value refers to the physical exercise value obtained at the present time. The history of physical exercise values is one or more physical exercise values representing the exercise stimuli given to the human body before the first time point. The plurality of physical exercise values at the time points before the first time point may include the physical exercise values at the first time point. The second physical exercise value refers to the predicted value of the current physical exercise value predicted from the past physical exercise value.

第1の身体運動値と第2の身体運動値とを用いて得られる入力値を演算モデルに入力することで、第1の時点での身体運動値(第1の身体運動値)に対して、第1の時点における身体運動値の予測値(第2の身体運動値)を考慮して動揺病の推定値が得られる。これにより、人の身体に与えられる運動が予測しやすい場合には酔いが少なく、予測しにくい場合には酔いやすいことを表した動揺病の推定値が得られる。 By inputting the input value obtained by using the first body exercise value and the second body exercise value into the calculation model, the body exercise value at the first time point (first body exercise value) can be obtained. , An estimated value of motion sickness can be obtained in consideration of the predicted value of the physical exercise value at the first time point (second physical exercise value). As a result, it is possible to obtain an estimated value of motion sickness, which indicates that the movement given to the human body is less sick when it is easy to predict, and it is easy to get sick when it is difficult to predict.

(2)好ましくは、第2の身体運動値を予測することは、身体運動値の履歴を解析することを含む。これにより、第2の身体運動値を算出することができる。 (2) Preferably, predicting the second physical exercise value includes analyzing the history of the physical exercise value. This makes it possible to calculate the second physical exercise value.

(3)好ましくは、第2の身体運動値を予測することにおいては、さらに、第1の時点に人の身体に与えられる運動刺激に影響を与える環境情報を用いる。これにより、第2の身体運動値を算出することができる。 (3) Preferably, in predicting the second physical exercise value, environmental information that affects the motor stimulus given to the human body at the first time point is used. This makes it possible to calculate the second physical exercise value.

(4)好ましくは、第1の身体運動値及び身体運動値の履歴は加速度を含み、第2の身体運動値は、加速度の予測値を含む。これにより、加速度の予測値を算出することができる。 (4) Preferably, the first body exercise value and the history of the body exercise value include the acceleration, and the second body exercise value includes the predicted value of the acceleration. This makes it possible to calculate the predicted value of acceleration.

(5)好ましくは、入力値は、第1の身体運動値に加速度の予測値を用いて得られる値を含む。第1の身体運動値に加速度の予測値を用いて得られる値は、一例として、第1の身体運動値と、信頼度を考慮した加速度の予測値とを含む。これにより、予測値を考慮して動揺病の推定値を算出することができる。 (5) Preferably, the input value includes a value obtained by using the predicted value of acceleration as the first body exercise value. The value obtained by using the predicted value of acceleration as the first body movement value includes, for example, the first body movement value and the predicted value of acceleration in consideration of reliability. This makes it possible to calculate an estimated value of motion sickness in consideration of the predicted value.

(6)好ましくは、第1の身体運動値及び身体運動値の履歴は角速度を含み、第2の身体運動値は、角速度の予測値を含む。これにより、角速度の予測値を算出することができる。 (6) Preferably, the first body exercise value and the history of the body exercise value include the angular velocity, and the second body exercise value includes the predicted value of the angular velocity. This makes it possible to calculate the predicted value of the angular velocity.

(7)好ましくは、入力値は、第1の身体運動値に角速度の予測値を用いて得られる値を含む。第1の身体運動値に加速度の予測値を用いて得られる値は、一例として、第1の身体運動値と、角速度の予測値とを含む。これにより、予測値を考慮して動揺病の推定値を算出することができる。 (7) Preferably, the input value includes a value obtained by using the predicted value of the angular velocity as the first body exercise value. The value obtained by using the predicted value of acceleration as the first body movement value includes, for example, the first body movement value and the predicted value of the angular velocity. This makes it possible to calculate an estimated value of motion sickness in consideration of the predicted value.

(8)好ましくは、第2の身体運動値を予測することは、身体運動値の履歴から身体運動値を予測する予測アルゴリズムを用いることを含む。予測アルゴリズムは、例えば、ガウス過程回帰である。これにより、第2の身体運動値を算出することができる。 (8) Preferably, predicting the second physical exercise value includes using a prediction algorithm for predicting the physical exercise value from the history of the physical exercise value. The prediction algorithm is, for example, Gaussian process regression. This makes it possible to calculate the second physical exercise value.

(9)好ましくは、演算モデルは、感覚器での運動感覚量に対応した値と、中枢神経系に構築されている感覚器の内部モデルで推定された運動感覚量である運動推定量に対応した値と、の差に基づいて動揺病の推定値を算出する。これにより、高精度で動揺病の推定値を算出できる。 (9) Preferably, the arithmetic model corresponds to a value corresponding to the amount of motion sensation in the sensory organ and a motion estimator which is the amount of motion sensation estimated by the internal model of the sensory organ constructed in the central nervous system. The estimated value of sway disease is calculated based on the difference between the value and the value. As a result, the estimated value of motion sickness can be calculated with high accuracy.

(10)好ましくは、入力値を演算モデルに入力することは、第2の身体運動値を、第2の身体運動値の予測の信頼度に応じて調整して入力することを含む。これにより、動揺病の推定値に適切に第2の身体運動値を反映させることができる。 (10) Preferably, inputting the input value into the calculation model includes adjusting and inputting the second body exercise value according to the reliability of the prediction of the second body exercise value. Thereby, the second physical exercise value can be appropriately reflected in the estimated value of motion sickness.

(11)好ましくは、第2の身体運動値を調整することは、信頼度が高いと演算モデルへの入力度合いを高くし、信頼度が低いと入力度合を低くすることを含む。これにより、動揺病の推定値に適切に第2の身体運動値を反映させることができる。 (11) Preferably, adjusting the second body exercise value includes increasing the degree of input to the calculation model when the reliability is high, and decreasing the degree of input when the reliability is low. Thereby, the second physical exercise value can be appropriately reflected in the estimated value of motion sickness.

(12)実施の形態に係る演算装置は、入力値から動揺病の推定値を算出する演算装置であって、第1の時点において人の身体に与えられた運動刺激を表す第1の身体運動値を入力する入力部と、人の動揺病の推定値を算出する演算部と、を備え、演算部は、第1の時点において人の身体に与えられる運動刺激を表す身体運動値の予測値である第2の身体運動値を、第1の時点以前の時点に人の身体に与えられた運動刺激を表す身体運動値の履歴を用いて予測し、第1の身体運動値と第2の身体運動値とに基づく値を入力値として、入力値から動揺病の推定値を算出する演算モデルに入力し、人の第1の時点における動揺病の推定値を算出する、ことを含む。これにより、演算装置では、(1)~(11)に記載の演算方法にて動揺病の推定値を算出することができる。 (12) The calculation device according to the embodiment is a calculation device that calculates an estimated value of motion sickness from an input value, and is a first physical exercise that represents an exercise stimulus given to a human body at a first time point. It includes an input unit for inputting values and a calculation unit for calculating an estimated value of human motion sickness, and the calculation unit is a predicted value of a physical exercise value representing an exercise stimulus given to the human body at a first time point. The second physical exercise value is predicted by using the history of the physical exercise value representing the exercise stimulus given to the human body before the first time point, and the first physical exercise value and the second physical exercise value are predicted. This includes inputting a value based on the physical exercise value into an arithmetic model for calculating an estimated value of motion sickness from the input value, and calculating an estimated value of motion sickness at a person's first time point. As a result, the arithmetic unit can calculate the estimated value of motion sickness by the arithmetic method described in (1) to (11).

<2.演算方法、及び、演算装置の例> <2. Calculation method and example of arithmetic unit>

[第1の実施の形態] [First Embodiment]

本実施の形態に係る演算装置は、動揺病の程度を算出する。動揺病の程度は、一例として動揺病の推定値などの指標値によって表される。本実施の形態に係る演算装置は動揺病の推定値を算出する。動揺病の推定値は、一例として、動揺病発症率(MSI:Motion Sickness Incidence)である。MSIは、物理刺激に対して嘔吐に至る人数の割合である。 The arithmetic unit according to this embodiment calculates the degree of motion sickness. The degree of motion sickness is represented by an index value such as an estimated value of motion sickness as an example. The arithmetic unit according to this embodiment calculates an estimated value of motion sickness. The estimated value of motion sickness is, for example, the incidence of motion sickness (MSI: Motion Sickness Incidence). MSI is the ratio of the number of people leading to vomiting to physical stimuli.

MSIは、例えば、車両などの乗り物の乗員の動揺病を防止、又は、程度を低減する制御に用いることができる。乗り物は、他の例として、船舶、電車、航空機、などである。動揺病が防止、又は、程度が軽減されることで、乗り物の快適性が向上する。そのため、より効果的に動揺病が防止、又は、程度が軽減されることが望ましい。そこで、一例として、乗り物の乗員を演算装置でMSIを算出するためのデータを測定する人とする。 MSI can be used, for example, for control to prevent or reduce motion sickness in occupants of vehicles such as vehicles. Vehicles are, for example, ships, trains, aircraft, and the like. The comfort of the vehicle is improved by preventing or reducing the degree of motion sickness. Therefore, it is desirable to prevent or reduce the degree of motion sickness more effectively. Therefore, as an example, let the occupant of the vehicle be a person who measures data for calculating MSI with an arithmetic unit.

発明者は、乗員の身体に与えられる運動が予測しやすい場合には酔いが少なく、予測しにくい場合には酔いやすいことに着目した。そこで、本実施の形態に係る演算装置では、従来の数理モデルに、このような運動の予想のしやすさを加味してMSIを算出する。 The inventor focused on less sickness when the movement given to the occupant's body is easy to predict, and easy to get sick when it is difficult to predict. Therefore, in the arithmetic unit according to the present embodiment, the MSI is calculated by adding the ease of predicting such a motion to the conventional mathematical model.

図1を参照して、本実施の形態に係る演算装置1は、プロセッサ11と、記憶部12と、を有する一般的なコンピュータなどで構成される。プロセッサ11は、例えば、CPU(Central Processing Unit)である。 With reference to FIG. 1, the arithmetic unit 1 according to the present embodiment is composed of a general computer or the like having a processor 11 and a storage unit 12. The processor 11 is, for example, a CPU (Central Processing Unit).

演算装置1は、センサインタフェース(I/F)14を備える。センサI/F14は、第1のセンサSE1及び第2のセンサSE2それぞれと有線、又は、無線で接続されて、センサ信号の入力を受け付ける。センサI/F14は、入力されたセンサ信号をプロセッサ11に入力する。 The arithmetic unit 1 includes a sensor interface (I / F) 14. The sensor I / F14 is connected to each of the first sensor SE1 and the second sensor SE2 by wire or wirelessly, and receives the input of the sensor signal. The sensor I / F 14 inputs the input sensor signal to the processor 11.

センサSE1は、前庭感覚器に対応したセンサである。センサSE1は、乗員の身体、例えば、頭部に与えられる運動刺激を検出し、運動刺激に対応した運動情報を出力する。センサSE1の出力信号I1は、運動情報を表す。以下、センサSE1からの出力を運動情報I1ともいう。 The sensor SE1 is a sensor corresponding to the vestibular sensory organ. The sensor SE1 detects the motor stimulus given to the occupant's body, for example, the head, and outputs the motor information corresponding to the motor stimulus. The output signal I1 of the sensor SE1 represents motion information. Hereinafter, the output from the sensor SE1 is also referred to as motion information I1.

運動情報I1は、頭部加速度及び頭部角度である。その他、頭部速度、又は、これらのうちの2以上の組み合わせ、などであってもよい。センサSE1は、例えば、乗員の頭部(例えば帽子)に取り付けられ、頭部加速度を計測する加速度センサ、及び、頭部角度を計測するジャイロセンサである。 The motion information I1 is the head acceleration and the head angle. In addition, the head speed or a combination of two or more of these may be used. The sensor SE1 is, for example, an acceleration sensor attached to the occupant's head (for example, a hat) and measuring the head acceleration, and a gyro sensor for measuring the head angle.

なお、運動情報I1は、乗員の頭部に与えられる運動刺激とみなされる運動に関する情報であってもよく、例えば、車両V自体の加速度及び角度であってもよい。車両Vの加速度及び角度が車両Vの乗員のそれと同じであると推定されるためである。この場合、第1のセンサSE1は車載の加速度センサ及びジャイロセンサである。 The motion information I1 may be information about a motion considered to be a motion stimulus given to the head of the occupant, and may be, for example, the acceleration and the angle of the vehicle V itself. This is because it is estimated that the acceleration and angle of the vehicle V are the same as those of the occupants of the vehicle V. In this case, the first sensor SE1 is an in-vehicle acceleration sensor and gyro sensor.

センサSE1は、予め規定された間隔(例えば10msec間隔)でセンシングを実行し、運動情報I1を演算装置1に対して出力する。 The sensor SE1 executes sensing at predetermined intervals (for example, 10 msec intervals) and outputs motion information I1 to the arithmetic unit 1.

記憶部12は、フラッシュメモリ、EEPROM(Electrically Erasable Programmable Read Only Memory)、ROM、RAM(Random Access Memory)などを含む。記憶部12は、1又は複数のプログラムからなる演算プログラム121を記憶する。 The storage unit 12 includes a flash memory, an EEPROM (Electrically Erasable Programmable Read Only Memory), a ROM, a RAM (Random Access Memory), and the like. The storage unit 12 stores an arithmetic program 121 composed of one or a plurality of programs.

プロセッサ11は、センサ信号に応じて記憶部12に記憶された演算プログラム121を読み出し、演算処理111を実行する演算部110として機能する。演算プログラム121は、CD-ROMや(Compact Disc Read only memory)やDVD-ROM(Digital Versatile Disk ROM)などの記録媒体に記録した状態で譲渡することもできるし、サーバコンピュータなどのコンピュータ装置からのダウンロードによって譲渡することもできる。演算プログラム121は、ウェブブラウザ上で動作するいわゆるウェブアプリケーションであってもよいし、プロセッサ11でのみ動作するいわゆる専用アプリケーションであってもよい。 The processor 11 functions as a calculation unit 110 that reads out the calculation program 121 stored in the storage unit 12 in response to the sensor signal and executes the calculation process 111. The arithmetic program 121 can be transferred in a state of being recorded on a recording medium such as a CD-ROM, (Compact Disc Read only memory) or DVD-ROM (Digital Versatile Disk ROM), or can be transferred from a computer device such as a server computer. It can also be transferred by downloading. The arithmetic program 121 may be a so-called web application that operates on a web browser, or may be a so-called dedicated application that operates only on the processor 11.

好ましくは、演算装置1は出力装置3に接続されている。出力装置3は、演算装置1での演算結果を出力可能な装置であって、例えば、ディスプレイである。プロセッサ11は、演算処理111の演算結果を出力装置3に渡し、出力させる。 Preferably, the arithmetic unit 1 is connected to the output unit 3. The output device 3 is a device capable of outputting the calculation result of the calculation device 1, and is, for example, a display. The processor 11 passes the calculation result of the calculation process 111 to the output device 3 and outputs the result.

演算処理111は、演算モデルである数理モデルMを用いてMSIを算出する処理である。数理モデルMは、人の身体に与えられた運動刺激を表す身体運動値から人の動揺病の推定値としてMSIを算出する演算モデルである。 The arithmetic processing 111 is a processing for calculating MSI using a mathematical model M which is an arithmetic model. The mathematical model M is an arithmetic model that calculates MSI as an estimated value of motion sickness in a person from a body movement value representing a movement stimulus given to the human body.

数理モデルMは、感覚器での運動感覚量に対応した値と、中枢神経系に構築されている感覚器の内部モデルで推定された運動感覚量である運動推定量に対応した値と、の差に基づいてMSIを算出する。感覚器は、前庭感覚を含む。 The mathematical model M has a value corresponding to the amount of motion sensation in the sensory organ and a value corresponding to the amount of motion estimation which is the amount of motion sensation estimated by the internal model of the sensory organ constructed in the central nervous system. The MSI is calculated based on the difference. Sensory organs include vestibular sensations.

図2を参照して、演算処理111を実行する演算部110には、センサSE1にて第1の時点で得られた、人の運動刺激を表す身体運動値(第1の身体運動値)が入力される。ここでの運動刺激は、物理的な運動刺激を含み、身体運動値は、運動刺激に対応する運動情報I1を含む。 With reference to FIG. 2, the calculation unit 110 that executes the calculation process 111 has a body movement value (first body movement value) that represents a human movement stimulus obtained at the first time point by the sensor SE1. Entered. The exercise stimulus here includes a physical exercise stimulus, and the body exercise value includes the exercise information I1 corresponding to the exercise stimulus.

演算処理111は感覚器処理(ステップS1)を含む。感覚器処理は、運動情報I1を用いて処理を行い、運動感覚量を示す運動感覚情報I3を出力することを含む。ここでの処理は、前庭感覚器の運動知覚の処理に対応した処理である。 The arithmetic processing 111 includes a sensory organ processing (step S1). The sensory organ processing includes processing using the motion information I1 and outputting the motion sensory information I3 indicating the amount of motion sensory. The processing here corresponds to the processing of the motion perception of the vestibular sensory organ.

演算処理111は内部モデル処理(ステップS2)を含む。内部モデル処理は、運動情報I1を用いて処理を行い、推定された運動感覚量である運動推定量を示す運動推定情報I4を出力することを含む。ここでの処理は、前庭感覚器の内部モデルでの運動推定の処理に対応した処理である。 The arithmetic processing 111 includes an internal model processing (step S2). The internal model processing includes processing using the motion information I1 and outputting the motion estimation information I4 indicating the motion estimation amount which is the estimated motion sensation amount. The processing here corresponds to the processing of motion estimation in the internal model of the vestibular sensory organ.

演算処理111は誤差算出処理(ステップS3)を含む。誤差算出処理は、運動感覚情報I3として出力された運動感覚量と、運動推定情報I4として出力された運動推定量との誤差を算出することを含む。 The calculation process 111 includes an error calculation process (step S3). The error calculation process includes calculating an error between the motion sensation amount output as the motion sensation information I3 and the motion estimation amount output as the motion estimation information I4.

演算処理111は動揺病程度算出処理(ステップS4)を含む。動揺病程度算出処理は、誤差算出処理で算出された誤差からMSIを算出することを含む。動揺病程度算出処理で得られた演算結果は、出力装置3に渡される。これにより、出力装置3で演算結果としてMSIが出力される。 The calculation process 111 includes a motion sickness degree calculation process (step S4). The motion sickness degree calculation process includes calculating MSI from the error calculated by the error calculation process. The calculation result obtained in the motion sickness degree calculation process is passed to the output device 3. As a result, the output device 3 outputs MSI as a calculation result.

演算処理111は、さらに、運動予測処理(ステップS5)を含む。運動予測処理は、センサSE1が第1の身体運動値を得た第1の時点における身体運動値の予測値(第2の身体運動値)を算出し、運動予測値I5を出力することを含む。予測値は、第1の時点以前の複数の時点それぞれに人の身体に与えられた運動刺激(運動刺激の履歴)を表す複数の身体運動値を用いて算出される。内部モデル処理においては、運動予測値I5をさらに用いて処理を行う。 The arithmetic process 111 further includes a motion prediction process (step S5). The exercise prediction process includes calculating a predicted value of the body movement value (second body movement value) at the first time point when the sensor SE1 obtains the first body movement value, and outputting the movement prediction value I5. .. The predicted value is calculated using a plurality of physical exercise values representing the exercise stimuli (history of exercise stimuli) given to the human body at each of the plurality of time points prior to the first time point. In the internal model processing, the processing is further performed using the motion prediction value I5.

好ましくは、運動予測処理は、さらに、予測値の信頼度を算出し、信頼度I6を出力することを含む。内部モデル処理においては、運動予測値I5に信頼度I6を考慮して用いてもよい。 Preferably, the motion prediction process further includes calculating the reliability of the predicted value and outputting the reliability I6. In the internal model processing, the motion prediction value I5 may be used in consideration of the reliability I6.

ここで、数理モデルMにて採用した、動揺病の程度の推定原理について説明する。感覚器の1つである前庭器において頭部運動が検知され、頭部運動に対応する信号(運動感覚信号)が生じる。また、中枢神経系に同様の過程を模擬した内部モデルが構築されていると考えられており、内部モデルでは、頭部運動の運動感覚量を示す信号(運動推定信号)が生成される。運動感覚信号と運動推定信号との間に誤差がある場合、つまり、感覚器で得られた運動感覚量と内部モデルで推定された運動感覚量との間に誤差がある場合、その誤差は感覚矛盾となる。感覚矛盾が動揺病を引き起こすとされている。 Here, the principle of estimating the degree of motion sickness adopted in the mathematical model M will be described. Head movement is detected in the vestibular organ, which is one of the sensory organs, and a signal corresponding to the head movement (motion sensory signal) is generated. Further, it is considered that an internal model simulating a similar process is constructed in the central nervous system, and the internal model generates a signal (motion estimation signal) indicating the amount of motor sensation of head movement. If there is an error between the kinesthetic signal and the kinesthetic estimation signal, that is, if there is an error between the kinesthetic amount obtained by the sensory organ and the kinesthetic amount estimated by the internal model, the error is sensory. It becomes a contradiction. Sensory contradictions are said to cause motion sickness.

さらに、演算処理111では、運動の予想のしやすさを加味するため、数理モデルMに入力する入力値として、第1の時点に得られた身体運動値(第1の身体運動値)と、第1の時点の身体運動値の予測値(第2の身体運動値)とを用いて得られる値を用いる。予測値は所定のゲインを乗じて用いられてもよい。所定のゲインは、一例として予測値の信頼度に基づくものであってもよい。以下の説明では、入力値は、身体運動値と、予測値を信頼度に応じて調整した値と、を含む。 Further, in the arithmetic processing 111, in order to take into account the ease of predicting the movement, the body movement value (first body movement value) obtained at the first time point is used as the input value to be input to the mathematical model M. The value obtained by using the predicted value of the physical exercise value at the first time point (second physical exercise value) is used. The predicted value may be used by multiplying it by a predetermined gain. The predetermined gain may be based on the reliability of the predicted value as an example. In the following description, the input value includes a physical exercise value and a value obtained by adjusting the predicted value according to the reliability.

図3を用いて、数理モデルMに従って演算処理111を実行する演算装置1の演算部110の機能構成を説明する。図3に示された各部は、プロセッサ11によって実現される。 FIG. 3 will explain the functional configuration of the arithmetic unit 110 of the arithmetic unit 1 that executes the arithmetic processing 111 according to the mathematical model M. Each part shown in FIG. 3 is realized by the processor 11.

図3を参照して、演算部110は処理部101を含む。処理部101は、感覚器としての前庭感覚器での運動知覚の処理に対応した処理を実行する。処理部101には、第1の時点において得られた運動情報I1が入力される。運動情報I1は、頭部座標系から見た頭部の重力と慣性力との合力である頭部加速度と、頭部角速度と、を含む。 With reference to FIG. 3, the arithmetic unit 110 includes a processing unit 101. The processing unit 101 executes a process corresponding to the process of motion perception in the vestibular sensory organ as a sensory organ. The motion information I1 obtained at the first time point is input to the processing unit 101. The motion information I1 includes a head acceleration, which is a resultant force of gravity and inertial force of the head as seen from the head coordinate system, and a head angular velocity.

処理部101での処理は、耳石器での運動知覚に対応した処理、及び、半規管での運動知覚に対応した処理を含む。処理部101は、運動情報I1が入力されると、前庭感覚器から出力される、運動感覚量を示す運動感覚信号に相当する運動感覚情報を出力する。 The process in the processing unit 101 includes a process corresponding to the motion perception in the otolith device and a process corresponding to the motion perception in the semicircular canal. When the motor information I1 is input, the processing unit 101 outputs the motor sensory information corresponding to the motor sensory signal indicating the motor sensory amount output from the vestibular sensory organ.

演算部110は運動予測部120を含む。運動予測部120には、運動情報I1が入力される。運動予測部120は、運動情報I1が入力されると、運動予測処理を実行して、第1の時点における身体運動の予測値(第2の身体運動値)を出力する。好ましくは、さらに、予測値の信頼度を出力してもよい。身体運動の予測値及び予測値の信頼度は、図2の運動予測値I5及び信頼度I6に相当する。 The calculation unit 110 includes a motion prediction unit 120. The motion information I1 is input to the motion prediction unit 120. When the exercise information I1 is input, the exercise prediction unit 120 executes the exercise prediction process and outputs the predicted value of the physical exercise (second physical exercise value) at the first time point. Preferably, the reliability of the predicted value may be further output. The predicted values of physical exercise and the reliability of the predicted values correspond to the predicted exercise values I5 and the reliability I6 in FIG.

演算部110は身体運動推定部103を含む。身体運動推定部103は、神経系において身体に与えられる刺激の遠心性コピーや体性感覚などに対応する情報(以下、コピー情報と略する)を生成する処理に対応した処理を実行する。詳しくは、身体運動推定部103には運動情報I1が入力される。また、身体運動推定部103には運動予測値I5及び信頼度I6も入力される。 The calculation unit 110 includes a body movement estimation unit 103. The body movement estimation unit 103 executes a process corresponding to a process of generating information (hereinafter, abbreviated as copy information) corresponding to an efferent copy of a stimulus given to the body in the nervous system, somatosensory, and the like. Specifically, the exercise information I1 is input to the body exercise estimation unit 103. Further, the exercise prediction value I5 and the reliability I6 are also input to the body exercise estimation unit 103.

身体運動推定部103は、運動予測値I5に対してゲインKを乗じる。ゲインKは一例として、予め記憶されているものである。これにより、運動予測値I5が定数倍されて、後述する内部モデルに入力される。 The body movement estimation unit 103 multiplies the motion prediction value I5 by the gain K. The gain K is stored in advance as an example. As a result, the motion prediction value I5 is multiplied by a constant and input to the internal model described later.

好ましくは、身体運動推定部103は、信頼度I6とゲインKとの関係を予め記憶しており、身体運動推定部103は、入力された信頼度I6に応じたゲインKを決定する。そして、運動予測値I5に対して、決定されたゲインKを乗じる。信頼度I6とゲインKとの関係は、信頼度をパラメータとした演算式であってもよい。ここでは、信頼度が高いほどゲインが大きくなる関係が記憶されている。 Preferably, the body movement estimation unit 103 stores the relationship between the reliability I6 and the gain K in advance, and the body movement estimation unit 103 determines the gain K according to the input reliability I6. Then, the determined gain K is multiplied by the motion prediction value I5. The relationship between the reliability I6 and the gain K may be an arithmetic expression with the reliability as a parameter. Here, the relationship that the gain increases as the reliability increases is stored.

身体運動推定部103は、運動情報I1に対して予め記憶しているゲインを乗じる。そして、身体運動推定部103は、運動情報I1から得られた値と運動予測値I5から得られた値とに基づく値をコピー情報とする。すなわち、身体運動推定部103は、運動予測値I5を信頼度I6で調整した情報を、上記のコピー情報とする。 The body movement estimation unit 103 multiplies the movement information I1 by a gain stored in advance. Then, the body exercise estimation unit 103 uses the value based on the value obtained from the exercise information I1 and the value obtained from the exercise prediction value I5 as copy information. That is, the body movement estimation unit 103 uses the information obtained by adjusting the movement prediction value I5 with the reliability I6 as the above-mentioned copy information.

身体運動推定部103には、後述する適応処理部109から出力された情報が入力される。身体運動推定部103は、コピー情報、及び、適応処理部109から出力された情報を統合する統合処理を実行する。統合処理は、自身の身体運動を把握するために人が行っている様々な感覚情報を統合する感覚統合に対応した処理である。身体運動推定部103は、統合処理によって出力される情報を内部モデルに出力する。そのため、ゲインが大きいほど、運動予測値I5が内部モデルに入力される度合が大きくなり、ゲインが小さくほど、運動予測値I5が内部モデルに入力される度合が小さくなる。コピー情報は、誤差を含む運動情報I1に相当する。 Information output from the adaptive processing unit 109, which will be described later, is input to the body movement estimation unit 103. The body movement estimation unit 103 executes an integrated process of integrating the copy information and the information output from the adaptive processing unit 109. The integrated process is a process corresponding to sensory integration that integrates various sensory information performed by a person in order to grasp one's own physical movement. The body movement estimation unit 103 outputs the information output by the integrated processing to the internal model. Therefore, the larger the gain, the larger the degree to which the motion prediction value I5 is input to the internal model, and the smaller the gain, the smaller the degree to which the motion prediction value I5 is input to the internal model. The copy information corresponds to the motion information I1 including the error.

なお、信頼度を用いた調整方法は、信頼度が高いほど入力される度合が大きくなり、信頼度が低いほど入力される度合が小さくなるような調整方法であれば他のどのような方法であってもよい。 As for the adjustment method using reliability, any other method can be used as long as the adjustment method is such that the higher the reliability, the larger the degree of input, and the lower the reliability, the smaller the degree of input. There may be.

演算部110は内部モデルでの運動推定の処理に対応した処理を実行する前庭感覚器モデル107を含む。前庭感覚器モデル107には、身体運動推定部103から出力された情報が入力される。身体運動推定部103から情報が入力されることで、前庭感覚器モデル107は、図2の運動推定情報I4に相当する情報を出力する。 The arithmetic unit 110 includes a vestibular sensory organ model 107 that executes a process corresponding to the process of motion estimation in the internal model. Information output from the body movement estimation unit 103 is input to the vestibular sensory organ model 107. When the information is input from the physical motion estimation unit 103, the vestibular sensory organ model 107 outputs the information corresponding to the motion estimation information I4 in FIG.

詳しくは、前庭感覚器モデル107は、身体運動推定部103から入力された情報を用いて内部モデルの前庭感覚器モデルでの運動推定の処理に対応した処理を実行する。前庭感覚器モデル107は、この処理によって運動感覚量を推定し、推定された運動感覚量である運動推定量を示す運動推定情報を出力する。 Specifically, the vestibular sensory organ model 107 executes processing corresponding to the processing of motion estimation in the vestibular sensory organ model of the internal model using the information input from the body movement estimation unit 103. The vestibular sensory organ model 107 estimates the kinesthetic amount by this processing, and outputs the motion estimation information indicating the motion estimation amount which is the estimated motion sensation amount.

演算部110は比較部108を含む。比較部108は加算器である。比較部108には、処理部101から出力された運動感覚情報、及び、前庭感覚器モデル107から出力された運動推定情報が入力される。 The arithmetic unit 110 includes a comparison unit 108. The comparison unit 108 is an adder. The kinesthetic information output from the processing unit 101 and the kinesthetic estimation information output from the vestibular sensory organ model 107 are input to the comparison unit 108.

比較部108は、運動感覚情報と運動推定情報とを用いて比較処理を実行する。これにより、比較部108は、処理部101で得られた運動感覚量と前庭感覚器モデル107で推定された運動感覚量(運動推定量)とを比較し、それらの差分(誤差)を算出する。算出された誤差は、感覚矛盾情報として出力される。 The comparison unit 108 executes a comparison process using the motion sensory information and the motion estimation information. As a result, the comparison unit 108 compares the kinesthetic amount obtained by the processing unit 101 with the kinesthetic sensory amount (motor estimation amount) estimated by the vestibular sensory organ model 107, and calculates the difference (error) between them. .. The calculated error is output as sensory contradiction information.

演算部110は適応処理部109を含む。適応処理部109には、比較部108から出力された感覚矛盾情報が入力される。適応処理部109は、入力された感覚矛盾情報を蓄積し、内部モデルでの処理に適応させるための適応処理を施した情報を出力する。 The arithmetic unit 110 includes an adaptive processing unit 109. The sensory contradiction information output from the comparison unit 108 is input to the adaptation processing unit 109. The adaptation processing unit 109 accumulates the input sensory contradiction information and outputs the information subjected to the adaptation processing for adapting to the processing in the internal model.

演算部110は感覚矛盾処理部125を含む。感覚矛盾処理部125には、比較部108から出力された感覚矛盾情報が入力される。感覚矛盾処理部125は、感覚矛盾情報を用いて動揺病の程度を示す指標値を算出するための演算式を予め記憶している。感覚矛盾処理部125は、記憶している演算式に比較部108から入力された感覚矛盾情報を代入することによってMSIを算出する。 The arithmetic unit 110 includes a sensory contradiction processing unit 125. The sensory contradiction information output from the comparison unit 108 is input to the sensory contradiction processing unit 125. The sensory contradiction processing unit 125 stores in advance an arithmetic expression for calculating an index value indicating the degree of agitation using the sensory contradiction information. The sensory contradiction processing unit 125 calculates the MSI by substituting the sensory contradiction information input from the comparison unit 108 into the stored arithmetic expression.

図4を参照して、数理モデルMに従って演算処理111を実行する演算部110には、センサSE1からの運動情報I1として、頭部にかかる重力加速度gと慣性加速度aとの和である頭部加速度f(f=g+a)と、頭部角速度ωと、慣性加速度aと、の第1の時点におけるセンシング結果が入力される。 With reference to FIG. 4, the arithmetic unit 110 that executes the arithmetic processing 111 according to the mathematical model M has the head, which is the sum of the gravitational acceleration g applied to the head and the inertial acceleration a, as the motion information I1 from the sensor SE1. The sensing results of the acceleration f (f = g + a), the head angular velocity ω, and the inertial acceleration a at the first time point are input.

処理部101は、処理部OTO、及び、処理部SCCを含む。処理部OTOには頭部加速度fが入力される。処理部OTOは、耳石器での運動知覚の処理に対応した処理を実行する。処理部SCCには頭部角速度ωが入力される。処理部SCCは、半規管での運動知覚の処理に対応した処理を実行する。 The processing unit 101 includes a processing unit OTO and a processing unit SCC. The head acceleration f is input to the processing unit OTO. The processing unit OTO executes a process corresponding to the process of motion perception in the otolith device. The head angular velocity ω is input to the processing unit SCC. The processing unit SCC executes processing corresponding to the processing of motion perception in the semicircular canals.

処理部OTO及び処理部SCCの伝達特性は、いずれも頭部固定座標系で記述された式で表される。処理部OTOの伝達特性は単位行列で表される。処理部SCCの伝達特性は図中の式(1)で表される。なお、式(1)中のτa、τdは時定数である。 The transmission characteristics of the processing unit OTO and the processing unit SCC are both expressed by the equations described in the head fixed coordinate system. The transfer characteristics of the processing unit OTO are represented by an identity matrix. The transmission characteristics of the processing unit SCC are represented by the equation (1) in the figure. Note that τa and τd in Eq. (1) are time constants.

処理部OTOからは信号f’が出力される。処理部SCCからは頭部角速度ωの感覚量ω’が出力される。信号f’及び信号ω’に対してローパスフィルタLPを適用することで、重力方向速度の感覚量vsが算出される。感覚量vsを用いて慣性加速度aの感覚量asが算出される。ローパスフィルタLPの特性は図中の式(2)で表される。なお、式(2)中のτは時定数である。感覚量as、感覚量ωs、及び、感覚量vsは、運動感覚情報に相当する。 The signal f'is output from the processing unit OTO. The processing unit SCC outputs the sensory amount ω'of the head angular velocity ω. By applying the low-pass filter LP to the signal f'and the signal ω', the perceived amount vs. the velocity in the direction of gravity is calculated. The sensory quantity as of the inertial acceleration a is calculated using the sensory quantity vs. The characteristics of the low-pass filter LP are represented by the equation (2) in the figure. Note that τ in Eq. (2) is a time constant. The sensory amount as, the sensory amount ωs, and the sensory amount vs. correspond to the kinesthetic information.

前庭感覚器モデル107は、耳石器の内部モデルに相当する処理部<OTO>と、半規管の内部モデルに相当する処理部<SCC>と、を含む。表記「<>」で挟まれた符号は、内部モデルに相当する処理を行う処理部を指している。以降の説明でも同様である。 The vestibular sensory organ model 107 includes a processing unit <OTO> corresponding to an internal model of an otolith and a processing unit <SCC> corresponding to an internal model of a semicircular canal. The code sandwiched between the notation "<>" refers to the processing unit that performs the processing corresponding to the internal model. The same applies to the following description.

処理部<OTO>の伝達特性は単位行列で表される。処理部<SCC>の伝達特性は図中の式(3)で表される。処理部<OTO>及び処理部<SCC>それぞれから出力される信号に対してローパスフィルタ<LP>を適用することで、慣性加速度aの推定量as^、頭部角速度ωの推定量ωs^、及び、重力方向速度の推定量vs^が得られる。ローパスフィルタ<LP>の特性もまた、図中の式(2)で表される。処理部<OTO>及び処理部<SCC>それぞれから出力される信号は、運動推定情報に相当する。 The transfer characteristics of the processing unit <OTO> are represented by an identity matrix. The transmission characteristics of the processing unit <SCC> are represented by the equation (3) in the figure. By applying the low-pass filter <LP> to the signals output from the processing unit <OTO> and the processing unit <SCC>, the estimated amount of inertial acceleration a, the estimated amount of head angular velocity ω, ωs ^, And the estimator vs ^ of the velocity in the direction of gravity is obtained. The characteristics of the low-pass filter <LP> are also represented by the equation (2) in the figure. The signals output from each of the processing unit <OTO> and the processing unit <SCC> correspond to motion estimation information.

比較部108は、処理部OTO及び処理部SCCそれぞれから出力された感覚量as、感覚量vs、及び、感覚量ωsと、処理部<OTO>及び処理部<SCC>それぞれから出力された推定量as^、推定量ωs^、及び、推定量vs^との誤差Δa、Δv、Δωを算出し、それぞれを示す信号を出力する。比較部108からの出力信号は、感覚矛盾情報に相当する。 The comparison unit 108 includes the sensory amount as, the sensory amount vs, and the sensory amount ωs output from each of the processing unit OTO and the processing unit SCC, and the estimated amount output from each of the processing unit <OTO> and the processing unit <SCC>. The errors Δa, Δv, and Δω with as ^, the estimator ωs ^, and the estimator vs ^ are calculated, and signals indicating each are output. The output signal from the comparison unit 108 corresponds to sensory contradiction information.

適応処理部109は、比較部108からの出力信号に示される誤差Δa、Δv、Δωをそれぞれ積分し、それぞれにゲインKac、Kωc、Kvcを乗じて得られた処理後の誤差Δa’、Δv’、Δω’を身体運動推定部103に入力する。 The adaptive processing unit 109 integrates the errors Δa, Δv, and Δω shown in the output signal from the comparison unit 108, respectively, and multiplies the gains Kac, Kωc, and Kvc, respectively, to obtain the error after processing Δa', Δv'. , ΔΩ'is input to the body movement estimation unit 103.

運動予測部120には、頭部加速度fが入力される。運動予測部120は、第1の時点における運動刺激を予測する処理を実行する。運動刺激を予測する処理は後述する。一例として、運動予測部120は、運動刺激のうちの慣性加速度aを予測する。この場合、運動予測部120は、その処理によって、予測された慣性加速度である予測慣性加速度a’を出力する。さらに、運動予測部120は、予測慣性加速度a’の信頼度を算出し、信頼度を示す信号dを出力する。 The head acceleration f is input to the motion prediction unit 120. The motion prediction unit 120 executes a process of predicting the motion stimulus at the first time point. The process of predicting motor stimuli will be described later. As an example, the motion prediction unit 120 predicts the inertial acceleration a in the motion stimulus. In this case, the motion prediction unit 120 outputs the predicted inertial acceleration a', which is the predicted inertial acceleration by the processing. Further, the motion prediction unit 120 calculates the reliability of the predicted inertial acceleration a'and outputs a signal d indicating the reliability.

身体運動推定部103には、予測慣性加速度a’及び信号dが入力される。また、身体運動推定部103には、頭部角速度ωが入力される。身体運動推定部103は、信号dに基づいてゲインKaを決定する。身体運動推定部103は、頭部角速度ωに、予め記憶しているゲインKωを乗じて頭部角速度ω~を得る。身体運動推定部103は、予測慣性加速度a’に決定されたゲインKaを乗じて慣性加速度a~を得る。頭部角速度ω~及び慣性加速度a~は、コピー情報に相当する。 The predicted inertial acceleration a'and the signal d are input to the body motion estimation unit 103. Further, the head angular velocity ω is input to the body movement estimation unit 103. The body movement estimation unit 103 determines the gain Ka based on the signal d. The body movement estimation unit 103 obtains the head angular velocity ω ~ by multiplying the head angular velocity ω by the gain Kω stored in advance. The body motion estimation unit 103 obtains the inertial acceleration a ~ by multiplying the predicted inertial acceleration a'by the determined gain Ka. The head angular velocity ω ~ and the inertial acceleration a ~ correspond to copy information.

身体運動推定部103の実行する統合処理は、一例として、加算処理である。そのため、身体運動推定部103は、加算器64,66,67をさらに含む。加算器64は、慣性加速度a~と誤差Δa’とを加算し、加算器66に渡す。加算器66は、慣性加速度a~と誤差Δa’との加算結果、及び、誤差Δv’を加算して得られた結果を示す信号を処理部<OTO>に入力する。加算器67は、頭部角速度ω~と誤差Δω’とを加算した結果を処理部<SCC>に入力する。 The integrated process executed by the body movement estimation unit 103 is, for example, an addition process. Therefore, the body movement estimation unit 103 further includes adders 64, 66, 67. The adder 64 adds the inertial acceleration a ~ and the error Δa'and passes it to the adder 66. The adder 66 inputs to the processing unit <OTO> a signal indicating the addition result of the inertial acceleration a ~ and the error Δa'and the result obtained by adding the error Δv'. The adder 67 inputs the result of adding the head angular velocity ω ~ and the error Δω'to the processing unit <SCC>.

感覚矛盾処理部125には、比較部108から出力された誤差Δvが入力される。誤差Δvは、主観的重力方向誤差(Subjective Vertical Conflict)とも呼ばれる。感覚矛盾処理部125は、誤差Δvに2次のHill関数及び2次遅れ伝達関数を適用する処理を実行する。この処理によって、MSIが出力される。 The error Δv output from the comparison unit 108 is input to the sensory contradiction processing unit 125. The error Δv is also called the Subjective Vertical Conflict. The sensory contradiction processing unit 125 executes a process of applying a quadratic Hill function and a quadratic lag transfer function to the error Δv. By this process, MSI is output.

なお、この例では、数理モデルMが、前庭感覚器で表された、1つの感覚器を用いるものとしている。他の例として、数理モデルMは、さらに他の感覚器も用いてもよい。すなわち、前庭感覚器に加えて他の感覚器についても運動感覚量を算出させ、その誤差に対応した値も用いてMSIを算出させてもよい。人には前庭感覚器のみならず複数種類の感覚器が存在し、人は、神経系において様々な感覚情報を統合して自身の身体運動を把握していると考えられるためである。他の感覚器は、例えば、視覚器である。視覚器では、視覚刺激に対する視覚器での運動感覚量に対応する信号(視覚感覚信号)が生じる。この場合、内部モデルには視覚器などの他の感覚器に対応した内部モデルも含まれると考えられる。このようにすることで、より高精度にMSIが算出されると考えられる。 In this example, it is assumed that the mathematical model M uses one sensory organ represented by the vestibular sensory organ. As another example, the mathematical model M may also use other sensory organs. That is, in addition to the vestibular sensory organ, the kinesthetic amount may be calculated for other sensory organs, and the MSI may be calculated using the value corresponding to the error. This is because humans have not only vestibular sensory organs but also multiple types of sensory organs, and it is thought that humans integrate various sensory information in the nervous system to grasp their own physical movements. Other sensory organs are, for example, visual organs. In the visual organ, a signal (visual sensory signal) corresponding to the amount of kinesthetic sensory in the visual organ to the visual stimulus is generated. In this case, the internal model is considered to include an internal model corresponding to other sensory organs such as a visual organ. By doing so, it is considered that MSI is calculated with higher accuracy.

運動予測部120での運動刺激を予測し、予測値の信頼度を算出する処理について説明する。一例として、運動予測部120での処理は、人の運動刺激の履歴の解析を含む。以前の時点において人に与えられる運動刺激を指し、1又は複数の身体運動値である。複数の身体運動値は、第1の時点の身体運動値を含んでもよい。 A process of predicting the motion stimulus in the motion prediction unit 120 and calculating the reliability of the predicted value will be described. As an example, the processing in the exercise prediction unit 120 includes analysis of the history of human exercise stimulation. It refers to the motor stimulus given to a person at a previous time point, and is one or more physical exercise values. The plurality of physical exercise values may include the physical exercise values at the first time point.

運動予測部120は、一例として、身体運動値の履歴から運動刺激の予測値が得られる予測アルゴリズムを予め記憶している。好ましくは、予測値の信頼度が得られる予測アルゴリズムもさらに記憶している。予測値及びその信頼度が得られる予測アルゴリズムは、一例としてガウス過程回帰である。 As an example, the exercise prediction unit 120 stores in advance a prediction algorithm for obtaining a predicted value of exercise stimulus from a history of physical exercise values. Preferably, the prediction algorithm for obtaining the reliability of the predicted value is also stored. The prediction algorithm for obtaining the predicted value and its reliability is, for example, Gaussian process regression.

例えば、図5に示されたような運動刺激の履歴が入力されたとする。図5は、人の身体に与えられる慣性加速度の時間変化を表している。図5の履歴では、慣性加速度の時間変化がパターン化され、一定周期で同じ変化が繰り返されている。 For example, it is assumed that the history of exercise stimulation as shown in FIG. 5 is input. FIG. 5 shows the time change of the inertial acceleration given to the human body. In the history of FIG. 5, the time change of the inertial acceleration is patterned, and the same change is repeated at a fixed cycle.

この場合、ガウス過程回帰モデルによって図6の信頼度が得られる。図6を参照して、慣性加速度の履歴の1周期目のうちは信頼度が低いものの、2周期目以降は周期的な慣性加速度の時間変化であることから、そのパターンに基づいた予測値の信頼度が高くなっていることがわかる。この場合、身体運動推定部103の決定する、予測値に乗じるためのゲインは大きくなる。 In this case, the Gaussian process regression model gives the reliability of FIG. With reference to FIG. 6, the reliability is low in the first cycle of the history of inertial acceleration, but since the second and subsequent cycles are periodic changes in inertial acceleration, the predicted values based on the pattern are obtained. It can be seen that the reliability is high. In this case, the gain for multiplying the predicted value determined by the body movement estimation unit 103 becomes large.

他の例として、図7に示された運動刺激の履歴が入力されたとする。図7も、人の身体に与えられる慣性加速度の時間変化を表している。図7の履歴では、慣性加速度の時間変化は周期的に加速度の変化が生じてはいるものの、パターン化されていない。 As another example, it is assumed that the history of the motor stimulus shown in FIG. 7 is input. FIG. 7 also shows the time change of the inertial acceleration given to the human body. In the history of FIG. 7, the time change of the inertial acceleration is not patterned although the change of the acceleration occurs periodically.

この場合、ガウス過程回帰モデルによって図8の信頼度が得られる。図8を参照して、この場合、全体に信頼度が一定せず、図6の信頼度と比較して低いことがわかる。この場合、身体運動推定部103の決定する、予測値に乗じるためのゲインは、図5の場合よりも小さくなる。 In this case, the Gaussian process regression model gives the reliability of FIG. With reference to FIG. 8, in this case, it can be seen that the reliability is not constant as a whole and is lower than the reliability of FIG. In this case, the gain for multiplying the predicted value determined by the body movement estimation unit 103 is smaller than that in the case of FIG.

図5及び図7の運動刺激の履歴を用いて、同じ運動刺激に対するMSIを算出すると、それぞれ、図9、図10の結果が得られた。すなわち、図5の運動刺激の履歴はパターン化されているため図6のように信頼度が高い。その結果、決定されるゲインが大きくなる。つまり、推定値の内部モデルに入力される度合が高くなる。この場合、図9に示されたように、MSIの大きさは全体的に0.6以下の低さであり、時間経過に沿った傾きも緩やかとなる。 When the MSI for the same exercise stimulus was calculated using the history of the exercise stimuli of FIGS. 5 and 7, the results of FIGS. 9 and 10 were obtained, respectively. That is, since the history of the motor stimulus in FIG. 5 is patterned, the reliability is high as shown in FIG. As a result, the determined gain is large. That is, the degree of input to the internal model of the estimated value is high. In this case, as shown in FIG. 9, the size of MSI is as low as 0.6 or less as a whole, and the inclination with the passage of time becomes gentle.

一方、図7の運動刺激の履歴の場合、パターン化されたものでないために図8のように信頼度が低い。その結果、決定されるゲインが小さくなる。つまり、推定値の内部モデルに入力される度合が低くなる。この場合、図10に示されたように、MSIの大きさは図0のMSIと比較して全体的に高く、時間経過に沿った傾きも大きくなる。 On the other hand, in the case of the history of the motor stimulus of FIG. 7, the reliability is low as shown in FIG. 8 because it is not patterned. As a result, the determined gain becomes smaller. That is, the degree of input to the internal model of the estimated value is low. In this case, as shown in FIG. 10, the size of the MSI is generally higher than that of the MSI of FIG. 0, and the inclination with the passage of time is also large.

図9及び図10の計算結果より、運動予測部120で予測値とその信頼度とを算出し、身体運動推定部103が信頼度に応じて上記のように予測値を用いることで、演算処理111の演算処理によって、身体に与えられる運動が予測しやすい場合には酔いが少なく、予測しにくい場合には酔いやすいことを表すMSIが得られることが分かった。 From the calculation results of FIGS. 9 and 10, the exercise prediction unit 120 calculates the predicted value and its reliability, and the body exercise estimation unit 103 uses the predicted value as described above according to the reliability, thereby performing arithmetic processing. It was found that the arithmetic processing of 111 can obtain MSI indicating that the movement given to the body is less sick when it is easy to predict, and it is easy to get sick when it is difficult to predict.

[第2の実施の形態] [Second Embodiment]

運動予測部120における予測値を算出する処理では、さらに、第1の時点において運動刺激に影響を与える環境情報を用いてもよい。また、予測値の信頼度を算出する処理も同様に、第1の時点において運動刺激に影響を与える環境情報を用いてもよい。 In the process of calculating the predicted value in the exercise prediction unit 120, environmental information that affects the exercise stimulus at the first time point may be further used. Similarly, in the process of calculating the reliability of the predicted value, environmental information that affects the motor stimulus at the first time point may be used.

第1の時点において運動刺激に影響を与える環境情報は、例えば、第1の時点に得られた交通状況画像や、交通情報や、地図情報や、走行計画などである。これら情報は、第1の時点において乗員や乗り物に取り付けられたセンサから得られたり、図示しない通信装置を用いてサーバなどの他の装置から取得したりされる。 The environmental information that affects the motor stimulus at the first time point is, for example, a traffic situation image, a traffic information, a map information, a traveling plan, etc. obtained at the first time point. This information is obtained from a sensor attached to a occupant or a vehicle at the first time point, or is obtained from another device such as a server using a communication device (not shown).

具体的に、交通状況画像の場合、演算装置1は、図1に示されたように、交通状況画像を得るためのセンサSE2にさらに接続されている。センサSE2は、一例としてカメラであって、例えば、乗員の頭部(例えば帽子)や車両などの乗り物の前方に取り付けられ、乗員の顔の向く方向を撮影する。センサSE2は、所定のタイミングで撮影を行い、画像情報I2を演算装置に対して出力する。 Specifically, in the case of a traffic situation image, the arithmetic unit 1 is further connected to the sensor SE2 for obtaining the traffic situation image, as shown in FIG. The sensor SE2 is, for example, a camera, which is attached to the front of a vehicle such as a occupant's head (for example, a hat) or a vehicle, and photographs the direction in which the occupant's face faces. The sensor SE2 takes an image at a predetermined timing and outputs the image information I2 to the arithmetic unit.

画像情報I2は、図2に示されたように、感覚器処理(ステップS1)及び運動予測処理(ステップS5)で用いられる。この場合、画像情報I2は、運動情報I1とともに処理部101、身体運動推定部103、及び、運動予測部120に入力される。 The image information I2 is used in the sensory organ processing (step S1) and the motion prediction processing (step S5) as shown in FIG. In this case, the image information I2 is input to the processing unit 101, the body exercise estimation unit 103, and the exercise prediction unit 120 together with the exercise information I1.

この場合、運動予測部120は、画像情報I2で表される画像に対して画像解析処理を実行する。ここでは、一例として、運動予測部120は、進行方向を抽出してもよい。例えば、運動予測部120は、センターラインを抽出することで、カーブの有無や大きさを識別する。これにより、乗員の身体に与えられる運動刺激を予測可能になる。 In this case, the motion prediction unit 120 executes an image analysis process on the image represented by the image information I2. Here, as an example, the motion prediction unit 120 may extract the traveling direction. For example, the motion prediction unit 120 identifies the presence / absence and size of a curve by extracting the center line. This makes it possible to predict the motor stimulus given to the occupant's body.

例えば、右にカーブしていることが識別された場合、運動予測部120は、乗員の身体に与えられる慣性加速度の向きを左向きと決定する。また、運動予測部120は、識別されたカーブの曲率が小さいほど乗員の身体に与えられる慣性加速度の大きさを大きくするよう、曲率に応じて慣性加速度の大きさを決定する。このようにすることでも、運動予測部120は運動刺激を予測することができる。 For example, when it is identified that the vehicle is curved to the right, the motion prediction unit 120 determines that the direction of the inertial acceleration given to the occupant's body is to the left. Further, the motion prediction unit 120 determines the magnitude of the inertial acceleration according to the curvature so that the smaller the curvature of the identified curve, the larger the magnitude of the inertial acceleration given to the occupant's body. By doing so, the motion prediction unit 120 can predict the motion stimulus.

なお、以上の説明では、運動予測部120で予測する運動刺激を表す身体運動値を慣性加速度としているが、慣性加速度に替えて、又は加えて、他の指標であってもよい。他の指標は、例えば、頭部加速度や角速度である。この場合も、上記と同様にして運動刺激を予測することができる。また、さらに、信頼度を算出することができる。 In the above description, the body movement value representing the movement stimulus predicted by the movement prediction unit 120 is defined as the inertial acceleration, but it may be replaced with or in addition to the inertial acceleration as another index. Other indicators are, for example, head acceleration and angular velocity. In this case as well, the motor stimulus can be predicted in the same manner as described above. Further, the reliability can be calculated.

<3.付記>
本発明は、上記実施形態に限定されるものではなく、様々な変形が可能である。
<3. Addendum>
The present invention is not limited to the above embodiment, and various modifications are possible.

1 :演算装置
3 :出力装置
11 :プロセッサ
12 :記憶部
14 :センサI/F
64 :加算器
66 :加算器
67 :加算器
101 :処理部
103 :身体運動推定部
107 :前庭感覚器モデル
108 :比較部
109 :適応処理部
110 :演算部
111 :演算処理
120 :運動予測部
121 :演算プログラム
125 :感覚矛盾処理部
I1 :運動情報
I2 :画像情報
I3 :運動感覚情報
I4 :運動推定情報
I5 :運動予測値
I6 :信頼度
K :ゲイン
Ka :ゲイン
Kac :ゲイン
Kvc :ゲイン
Kω :ゲイン
Kωc :ゲイン
LP :ローパスフィルタ
M :数理モデル
OTO :処理部
SCC :処理部
SE1 :センサ
SE2 :センサ
a :慣性加速度
a’ :予測慣性加速度
as :感覚量
as^ :推定量
d :信号
f :頭部加速度
f’ :信号
g :重力加速度
vs :感覚量
vs^ :推定量
Δa :誤差
Δa’ :誤差
Δv :誤差
Δv’ :誤差
Δω’ :誤差
ω :頭部角速度
ω’ :感覚量
ωs :感覚量
ωs^ :推定量
1: Arithmetic logic unit 3: Output device 11: Processor 12: Storage unit 14: Sensor I / F
64: Adder 66: Adder 67: Adder 101: Processing unit 103: Body movement estimation unit 107: Front yard sensory device model 108: Comparison unit 109: Adaptive processing unit 110: Calculation unit 111: Calculation processing 120: Motion prediction unit 121: Arithmetic program 125: Sensory inconsistency processing unit I1: Motion information I2: Image information I3: Motion sensory information I4: Motion estimation information I5: Motion prediction value I6: Reliability K: Gain Ka: Gain Kac: Gain Kvc: Gain Kω : Gain Kωc: Gain LP: Low pass filter M: Mathematical model OTO: Processing unit SCC: Processing unit SE1: Sensor SE2: Sensor a: Inertial acceleration a': Predicted inertial acceleration as: Sensory amount as ^: Estimated amount d: Signal f : Head acceleration f': Signal g: Gravity acceleration vs: Sensory amount vs ^: Estimated amount Δa: Error Δa': Error Δv: Error Δv': Error Δω': Error ω: Head angular velocity ω': Sensory amount ωs : Sensory amount ωs ^: Estimated amount

Claims (12)

入力値から動揺病の推定値を算出する演算モデルを用いた演算方法において、
第1の時点において人の身体に与えられた運動刺激を表す第1の身体運動値を取得し、
前記第1の時点において前記人の身体に与えられる運動刺激を表す身体運動値の予測値である第2の身体運動値を、前記第1の時点以前の時点に前記人の身体に与えられた運動刺激を表す身体運動値の履歴を用いて予測し、
前記第1の身体運動値と前記第2の身体運動値とに基づく値を前記入力値として前記演算モデルに入力し、前記人の前記第1の時点における前記動揺病の推定値を算出する、ことを含む
演算方法。
In the calculation method using the calculation model that calculates the estimated value of motion sickness from the input value
Obtain the first physical exercise value representing the exercise stimulus given to the human body at the first time point,
A second body exercise value, which is a predicted value of a body exercise value representing an exercise stimulus given to the person's body at the first time point, was given to the person's body at a time point before the first time point. Predict using the history of physical exercise values that represent exercise stimuli,
A value based on the first physical exercise value and the second physical exercise value is input to the calculation model as the input value, and an estimated value of the motion sickness at the first time point of the person is calculated. Calculation method including that.
前記第2の身体運動値を予測することは、前記身体運動値の履歴を解析することを含む
請求項1に記載の演算方法。
The calculation method according to claim 1, wherein predicting the second physical exercise value includes analyzing the history of the physical exercise value.
前記第2の身体運動値を予測することにおいては、さらに、前記第1の時点に前記人の身体に与えられる前記運動刺激に影響を与える環境情報を用いる
請求項1に記載の演算方法。
The calculation method according to claim 1, further comprising environmental information affecting the exercise stimulus given to the human body at the first time point in predicting the second physical exercise value.
前記第1の身体運動値及び前記身体運動値の履歴は加速度を含み、
前記第2の身体運動値は、前記加速度の予測値を含む
請求項1~3のいずれか一項に記載の演算方法。
The first physical exercise value and the history of the physical exercise value include acceleration, and the history includes acceleration.
The calculation method according to any one of claims 1 to 3, wherein the second body exercise value includes a predicted value of the acceleration.
前記入力値は、前記第1の身体運動値に前記加速度の予測値を用いて得られる値を含む
請求項4に記載の演算方法。
The calculation method according to claim 4, wherein the input value includes a value obtained by using the predicted value of the acceleration in the first body exercise value.
前記第1の身体運動値及び前記身体運動値の履歴は角速度を含み、
前記第2の身体運動値は、前記角速度の予測値を含む
請求項1~3のいずれか一項に記載の演算方法。
The first physical exercise value and the history of the physical exercise value include the angular velocity.
The calculation method according to any one of claims 1 to 3, wherein the second physical exercise value includes a predicted value of the angular velocity.
前記入力値は、前記第1の身体運動値に前記角速度の予測値を用いて得られる値を含む
請求項6に記載の演算方法。
The calculation method according to claim 6, wherein the input value includes a value obtained by using the predicted value of the angular velocity in the first body exercise value.
前記第2の身体運動値を予測することは、身体運動値の履歴から身体運動値を予測する予測アルゴリズムを用いることを含む
請求項1~7のいずれか一項に記載の演算方法。
The calculation method according to any one of claims 1 to 7, wherein predicting the second physical exercise value includes using a prediction algorithm for predicting the physical exercise value from the history of the physical exercise value.
前記演算モデルは、感覚器での運動感覚量に対応した値と、中枢神経系に構築されている前記感覚器の内部モデルで推定された運動感覚量である運動推定量に対応した値と、の差に基づいて前記動揺病の推定値を算出する
請求項1~8のいずれか一項に記載の演算方法。
The arithmetic model includes a value corresponding to the amount of motor sensation in the sensory organ, and a value corresponding to the amount of motor sensation estimated by the internal model of the sensory organ constructed in the central nervous system. The calculation method according to any one of claims 1 to 8, wherein the estimated value of the sway disease is calculated based on the difference between the two.
前記入力値を前記演算モデルに入力することは、前記第2の身体運動値を、前記第2の身体運動値の予測の信頼度に応じて調整して入力することを含む
請求項1~9のいずれか一項に記載の演算方法。
Inputting the input value into the calculation model includes adjusting and inputting the second body exercise value according to the reliability of prediction of the second body exercise value. Claims 1 to 9. The calculation method described in any one of the above items.
前記第2の身体運動値を調整することは、前記信頼度が高いと前記演算モデルへの入力度合いを高くし、前記信頼度が低いと前記入力度合を低くすることを含む
請求項10に記載の演算方法。
The tenth aspect of claim 10, wherein adjusting the second body exercise value includes increasing the input degree to the calculation model when the reliability is high and decreasing the input degree when the reliability is low. Calculation method.
入力値から動揺病の推定値を算出する演算装置であって、
第1の時点において人の身体に与えられた運動刺激を表す第1の身体運動値を入力する入力部と、
前記人の前記動揺病の推定値を算出する演算部と、を備え、
前記演算部は、
前記第1の時点において前記人の身体に与えられる運動刺激を表す身体運動値の予測値である第2の身体運動値を、前記第1の時点以前の時点に前記人の身体に与えられた運動刺激を表す身体運動値の履歴を用いて予測し、
前記第1の身体運動値と前記第2の身体運動値とに基づく値を前記入力値として、前記入力値から前記動揺病の推定値を算出する演算モデルに入力し、前記人の前記第1の時点における前記動揺病の推定値を算出する、ことを含む
演算装置。
It is an arithmetic unit that calculates the estimated value of motion sickness from the input value.
An input unit for inputting a first physical exercise value representing an exercise stimulus given to the human body at the first time point,
A calculation unit for calculating the estimated value of the motion sickness of the person is provided.
The arithmetic unit
A second body exercise value, which is a predicted value of a body exercise value representing an exercise stimulus given to the person's body at the first time point, was given to the person's body at a time point before the first time point. Predict using the history of physical exercise values that represent exercise stimuli,
The value based on the first physical exercise value and the second physical exercise value is used as the input value, and is input to the arithmetic model for calculating the estimated value of the motion sickness from the input value, and the person's first body exercise value is calculated. An arithmetic unit that includes calculating the estimated value of the motion sickness at the time point in time.
JP2020167162A 2020-10-01 2020-10-01 Arithmetic method and arithmetic device Pending JP2022059422A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020167162A JP2022059422A (en) 2020-10-01 2020-10-01 Arithmetic method and arithmetic device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020167162A JP2022059422A (en) 2020-10-01 2020-10-01 Arithmetic method and arithmetic device

Publications (1)

Publication Number Publication Date
JP2022059422A true JP2022059422A (en) 2022-04-13

Family

ID=81124217

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020167162A Pending JP2022059422A (en) 2020-10-01 2020-10-01 Arithmetic method and arithmetic device

Country Status (1)

Country Link
JP (1) JP2022059422A (en)

Similar Documents

Publication Publication Date Title
Asadi et al. A review on otolith models in human perception
CN112673378B (en) Device for generating estimator, monitoring device, method for generating estimator, and program for generating estimator
Irmak et al. Objective and subjective responses to motion sickness: the group and the individual
JP7046259B2 (en) Motion sickness estimation device, motion sickness suppression device and motion sickness estimation method
US20210154430A1 (en) Systems and methods for predicting and preventing motion sickness
US10766483B2 (en) Active vehicle virtual reality prevention of predictive motion sickness
US20220001894A1 (en) System for controlling autonomous vehicle for reducing motion sickness
JP7177465B2 (en) Evaluation device, control device, motion sickness reduction system, evaluation method, and computer program
CN110770107B (en) Method and device for evaluating the fatigue of an occupant of a vehicle
JP2022059422A (en) Arithmetic method and arithmetic device
JPH10281945A (en) Riding comfort evaluating device
JP7217874B2 (en) EVALUATION DEVICE, CONTROL DEVICE, SYSTEM FOR IMPROVING OPERATING SECONDITY, EVALUATION METHOD, AND COMPUTER PROGRAM
JP7183782B2 (en) Emotion estimation device, environment providing system, vehicle, emotion estimation method, and information processing program
JP6832005B2 (en) Subject judgment device, method, and program
CN116244954A (en) Vehicle comfort evaluation method, device, computer equipment and storage medium
CN107531249B (en) Method and device for detecting a change in the behaviour of a driver of a motor vehicle
US10299715B2 (en) Stress monitoring for individuals in moving structures
CN111209645B (en) Dizziness somatosensory modeling data acquisition device and method, terminal and storage medium
JP4681571B2 (en) Driving simulator and driving simulation method
Atsumi et al. Evaluation of vehicle motion sickness due to vehicle vibration
JP2024018713A (en) Information processing device, information processing method, and information processing program
CN117113856B (en) Dummy model determination method and system suitable for pre-collision scene
Kim et al. Modeling motion sickness using a four-wheel vehicle model augmented with a passenger model
CN117064386A (en) Method, apparatus, device, medium and program product for determining perceived reaction time
WO2022230629A1 (en) Electronic device, control method for electronic device, and program

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20210716

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20210716

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230523

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231117

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231212

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240122

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240416