JPWO2020220564A5 - - Google Patents

Download PDF

Info

Publication number
JPWO2020220564A5
JPWO2020220564A5 JP2021562859A JP2021562859A JPWO2020220564A5 JP WO2020220564 A5 JPWO2020220564 A5 JP WO2020220564A5 JP 2021562859 A JP2021562859 A JP 2021562859A JP 2021562859 A JP2021562859 A JP 2021562859A JP WO2020220564 A5 JPWO2020220564 A5 JP WO2020220564A5
Authority
JP
Japan
Prior art keywords
abr
detected
signal
current
time lag
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
JP2021562859A
Other languages
Japanese (ja)
Other versions
JP2022534851A (en
Publication date
Priority claimed from CN201910362550.0A external-priority patent/CN109998539B/en
Application filed filed Critical
Publication of JP2022534851A publication Critical patent/JP2022534851A/en
Publication of JPWO2020220564A5 publication Critical patent/JPWO2020220564A5/ja
Priority to JP2023178139A priority Critical patent/JP2023179702A/en
Pending legal-status Critical Current

Links

Claims (11)

開始時の音響強度で対応するABR記録を取得するステップA1と、
適応平均法の演算に基づき、ABR信号の検出の条件を満たすまで、反復によって記録の平均回数を段階的に増加させ、
前記ABR信号の検出の条件が、現在の反復時に現在の音響強度に対応するABR記録をグループ化した後、現在の平均回数に基づいて各グループの平均曲線をそれぞれ計算し、グループの間の相互相関関数の最大値に対応するタイムラグを取得し、タイムラグのずれが所定範囲内にあるか否かに基づき、タイムロック特性を有するABR信号が存在するか否かを判断するステップA2と、
ステップA2でABR信号が検出されているか否かを判断し、ABR信号が検出されている場合、ステップA4を実行し、ABR信号が検出されていない場合、ステップA5を実行するステップA3と、
現在の音響強度が設定された最低測定音響強度に達しているか否かを判断し、最低測定音響強度に達していない場合、新たな音響強度で対応するABR記録を取得して元データと組み合わせ、ステップA2を再実行し、前記新たな音響強度は、現在の音響強度から設定された間隔を差し引いた後に得られるものであり、最低測定音響強度に達している場合、A6を実行するステップA4と、
連続するP個の音響強度でABR信号が検出されていないか否かを判断し、連続するP個の音響強度でABR信号が検出されていない場合、ステップA6を実行し、そうでなければ、ステップA4の実行に移行するステップA5と、
ABR信号を検出するために必要とする最低音響強度を可聴閾値、または、適応平均法の演算において、開始時の音響強度から測定した音響強度の最小値までの各音響強度で使用される反復回数に対して関数フィッティングおよび補間法によって可聴閾値、に対応する音響強度を取得するステップA6と、
を含むことを特徴とする聴性脳幹反応の自動測定方法。
step A1 of acquiring a corresponding ABR recording at the starting sound intensity;
stepwise increasing the number of recordings averaged by iteration until the conditions for detection of the ABR signal are met, based on adaptive averaging calculations;
After the ABR signal detection conditions group the ABR recordings corresponding to the current sound intensity at the current iteration, respectively calculate the average curve of each group based on the current average number of times, and the interaction between the groups a step A2 of acquiring a time lag corresponding to the maximum value of the correlation function and determining whether or not there is an ABR signal having a time-lock characteristic based on whether or not the deviation of the time lag is within a predetermined range;
Step A3 of determining whether or not the ABR signal is detected in step A2, executing step A4 if the ABR signal is detected, and executing step A5 if the ABR signal is not detected;
determining whether the current sound intensity has reached the set minimum measured sound intensity, and if not, obtaining the corresponding ABR recording at the new sound intensity and combining it with the original data; re-executing step A2, wherein said new sound intensity is obtained after subtracting a set interval from the current sound intensity, and if the lowest measured sound intensity has been reached, step A4 of performing A6; ,
determining whether the ABR signal is not detected in the P consecutive sound intensities, if the ABR signal is not detected in the P consecutive sound intensities, perform step A6; otherwise, a step A5 that transitions to the execution of step A4;
The minimum sound intensity required to detect an ABR signal is the audible threshold, or the number of iterations used at each sound intensity from the starting sound intensity to the minimum measured sound intensity in the adaptive averaging calculation. obtaining a sound intensity corresponding to the audible threshold, A6, by function fitting and interpolation for
A method for automatically measuring an auditory brainstem response, comprising:
前記ステップA2の適応平均法の演算が、
前記ABR記録が、入力タイプが繰り返される単一記録である時間曲線であり、それをランダムに二つのグループに分けかつ平均曲線をそれぞれ計算するステップS1’と、
二つのグループのそれぞれが平均された後の相互相関関数を計算するステップS2’と、
相互相関関数の最大値に対応するタイムラグを取得するステップS3’と、
タイムラグの絶対値がk値よりも小さい場合、タイムラグずれがk個のデータポイント内にあることを意味し、安定したタイムロック信号が検出されていることを表し、この場合、データの取得と収集を停止し、ABR信号が検出されている場合に対応し、
タイムラグの絶対値がk値より大きい場合、タイムラグずれがk個のデータポイント内にないことを意味し、安定したタイムロック信号が検出されていないことを表し、この場合、さらにステップS5’を実行することを特徴とする、
タイムラグの絶対値とk値とを比較し、タイムラグずれがk個のデータポイント内にあるか否かを判断するステップS4’と、
最大反復回数に達している場合、データの取得と収集を停止し、ABR信号が検出されていない場合に対応し、最大反復回数に達していない場合、さらに現在収集されたABR記録に新たに収集されたABR記録を加え、現在の音響強度でS1’~S5’の演算プロセスを繰り返し実行することを特徴とする、
現在の音響強度での最大反復回数に達しているか否かを判断するステップS5’と、を含むことと、
kが予め設定された固定値であり、または、全てのデータポイントを予め設定された比率で計算した後に得られる値であることと、を特徴とする請求項1に記載の聴性脳幹反応の自動測定方法。
The calculation of the adaptive averaging method in step A2 is
a step S1' of dividing the ABR record into two groups randomly and calculating an average curve respectively, wherein the ABR record is a time curve whose input type is a single record repeated;
a step S2' of calculating the cross-correlation function after each of the two groups has been averaged;
Step S3' of obtaining a time lag corresponding to the maximum value of the cross-correlation function;
If the absolute value of the time lag is smaller than the k value, it means that the time lag deviation is within k data points, indicating that a stable time-locked signal has been detected, in which case data acquisition and collection and respond if an ABR signal is detected,
If the absolute value of the time lag is greater than the k value, it means that the time lag deviation is not within k data points, indicating that a stable time lock signal has not been detected, in which case further execute step S5′. characterized by
a step S4' of comparing the absolute value of the time lag with the k value and determining whether the time lag deviation is within k data points;
If the maximum number of iterations has been reached, stop data acquisition and collection, corresponding to the case where no ABR signal is detected, and if the maximum number of iterations has not been reached, then a new acquisition over the currently acquired ABR record. adding the recorded ABR recording and repeatedly performing the calculation process of S1' to S5' at the current acoustic intensity,
determining whether the maximum number of iterations at the current acoustic intensity has been reached;
2. The automatic auditory brainstem response of claim 1, wherein k is a preset fixed value or a value obtained after calculating all data points at a preset ratio. Measuring method.
前記ステップA2の適応平均法の演算が、
前記ABR記録が、入力タイプが繰り返される単一記録である時間曲線であり、それをQ回並行して判断し、各回に以下のステップS1~S4を実行することと、
現在収集されたABR記録をランダムに二つのグループに分けかつ平均曲線をそれぞれ計算するステップS1と、
二つのグループのそれぞれが平均された後の相互相関関数を計算するステップS2と、
相互相関関数の最大値に対応するタイムラグを取得するステップS3と、
タイムラグの絶対値とk値とを比較し、タイムラグずれがk個のデータポイント内にあるか否かを判断し、
タイムラグの絶対値がk値よりも小さい場合、タイムラグずれがk個のデータポイント内にあることを意味し、安定したタイムロック信号が検出されていることを表し、タイムラグの絶対値がk値以上の場合、タイムラグずれがk個のデータポイント内にないことを意味し、安定したタイムロック信号が検出されていないことを表すステップS4と、
Q回の並行判断のうちの各回にいずれも安定したタイムロック信号が検出されているか否かを判断し、
Q回のうちの各回にいずれも安定したタイムロック信号が検出されている場合、データの取得と収集を停止し、ABR信号が検出されている場合に対応し、
Q回のうちの各回に必ずしも安定したタイムロック信号が検出されていない場合、さらにステップS6を実行するステップS5と、
現在の音響強度での最大反復回数に達しているか否かを判断し、
最大反復回数に達している場合、データの取得と収集を停止し、ABR信号が検出されていない場合に対応し、
最大反復回数に達していない場合、さらに現在収集されたABR記録に新たに収集されたABR記録を加え、現在の音響強度でS1~S6の演算プロセスを繰り返し実行するステップS6と
を含むことと、
kが、予め設定された固定値であり、または、全てのデータポイントを予め設定された比率で計算した後に得られる値であることと、
を特徴とする請求項1に記載の聴性脳幹反応の自動測定方法。
The calculation of the adaptive averaging method in step A2 is
wherein said ABR record is a time curve whose input type is a single record repeated, determining it Q times in parallel, each time performing steps S1-S4 below;
step S1 of randomly dividing the currently collected ABR recordings into two groups and calculating mean curves respectively;
a step S2 of calculating the cross-correlation function after each of the two groups has been averaged;
Step S3 of obtaining a time lag corresponding to the maximum value of the cross-correlation function;
comparing the absolute value of the time lag and the k value to determine whether the time lag deviation is within k data points;
If the absolute value of the time lag is smaller than the k value, it means that the time lag deviation is within k data points, indicating that a stable timelock signal is detected, and the absolute value of the time lag is greater than or equal to the k value. if , step S4, meaning that the time lag deviation is not within k data points, indicating that a stable timelock signal has not been detected;
determining whether or not a stable time lock signal is detected in each of the Q parallel determinations;
stopping data acquisition and collection when a stable timelock signal is detected each of the Q times, corresponding to when an ABR signal is detected;
step S5 of further executing step S6 if a stable timelock signal is not necessarily detected each of the Q times;
Determine whether the maximum number of iterations at the current acoustic intensity has been reached,
if the maximum number of iterations has been reached, stop data acquisition and collection, and respond if no ABR signal is detected;
if the maximum number of iterations has not been reached, further adding newly acquired ABR records to the currently acquired ABR records and repeatedly performing the computational processes of S1-S6 at the current acoustic intensity;
k is a preset fixed value or a value obtained after calculating all data points with a preset ratio;
The method for automatically measuring an auditory brainstem response according to claim 1, characterized by:
前記ステップA2の適応平均法の演算が、
前記ABR記録が、入力タイプが複数回の平均記録である時間曲線であり、今回の反復で新たに加えられる平均曲線avgA加を、前回の反復の平均曲線avgA旧と組み合わせ、重みに従って現在の平均曲線avgA新を計算し、
Figure 2020220564000001
ここでは、avgA新、avgA旧、avgA加がそれぞれQグループがあり、Mが、現在の反復回数であるステップS1”と、
Qグループの現在の平均曲線avgA新に対して、二グループ間毎の相互相関関数を計算するステップS2”と、
相互相関関数の最大値に対応するQグループのタイムラグを取得するステップS3”と、
Q回の並行判断において、各グループのタイムラグの絶対値をそれぞれk値と比較し、タイムラグずれがk個のデータポイント内にあるか否かを判断し、タイムラグの絶対値がk値よりも小さい場合、タイムラグずれがk個のデータポイント内にあることを意味し、安定したタイムロック信号が検出されていることを表し、タイムラグの絶対値がk値以上の場合、タイムラグずれがk個のデータポイント内にないことを意味し、安定したタイムロック信号が検出されていないことを示すステップS4”と、
Q回の並行判断のうちの各回にいずれも安定したタイムロック信号が検出されているか否かを判断し、
Q回のうちの各回にいずれも安定したタイムロック信号が検出されている場合、データの取得と収集を停止し、ABR信号が検出されている場合に対応し、
Q回のうちの各回に必ずしも安定したタイムロック信号が検出されていない場合、さらにステップS6”を実行するステップS5”と、
現在の音響強度での最大反復回数に達しているか否かを判断し、
最大反復回数に達している場合、データの取得と収集を停止し、ABR信号が検出されていない場合に対応し、最大反復回数に達していない場合、さらに現在収集されたABR測定平均データに新たに収集されたABR測定平均データを加え、現在の音響強度でS1”~S6”の演算プロセスを繰り返し実行するステップS6”とを含むことと、
kが、予め設定された固定値であり、または、全てのデータポイントを予め設定された比率で計算した後に得られる値であることと、を特徴とする請求項1に記載の聴性脳幹反応の自動測定方法。
The calculation of the adaptive averaging method in step A2 is
The ABR record is a time curve whose input type is multiple average records, and the new average curve avgA added in this iteration is combined with the average curve avgA of the previous iteration, and the current average Calculate the curve avgA new,
Figure 2020220564000001
Here, there are Q groups of avgA new, avgA old, and avgA addition, and M is the current number of iterations in step S1'';
Step S2'' of calculating the cross-correlation function between each two groups for the current average curve avgA new of the Q group;
Step S3″ of obtaining the time lag of the Q group corresponding to the maximum value of the cross-correlation function;
In Q parallel judgments, the absolute value of the time lag of each group is respectively compared with the k value, and it is judged whether the time lag deviation is within k data points, and the absolute value of the time lag is less than the k value. If the time lag is within k data points, it means that a stable time lock signal has been detected. step S4'', meaning not within a point, indicating that a stable timelock signal has not been detected;
determining whether or not a stable time lock signal is detected in each of the Q parallel determinations;
stopping data acquisition and collection when a stable timelock signal is detected each of the Q times, corresponding to when an ABR signal is detected;
step S5'' of further executing step S6'' if a stable time-lock signal is not necessarily detected each of the Q times;
Determine whether the maximum number of iterations at the current acoustic intensity has been reached,
If the maximum number of iterations has been reached, data acquisition and collection is stopped, corresponding to the case where no ABR signal is detected, and if the maximum number of iterations has not been reached, the currently collected ABR measurement average data is updated. adding the collected averaged ABR measurement data to the current acoustic intensity and repeatedly performing the computational process of S1″-S6″ at the current acoustic intensity;
2. The auditory brainstem response of claim 1, wherein k is a preset fixed value or a value obtained after calculating all data points at a preset ratio. Automatic measurement method.
各音響強度で前記ABR記録平均曲線を計算する時に使用される平均回数が、反復回数とN値との積であり、N値が、各回の反復時に新たに加えられるABR記録のグループ数に対応することを特徴とする請求項1~4のいずれか一項に記載の聴性脳幹反応の自動測定方法。 The number of averages used in calculating the ABR recording average curve at each sound intensity is the product of the number of iterations and the N value, where N corresponds to the number of groups of ABR recordings newly added at each iteration. The method for automatically measuring an auditory brainstem response according to any one of claims 1 to 4, characterized in that: 前記ステップA6の各音響強度で使用される反復回数を正規化した後、シグモイド(Sigmoid)関数フィッティングおよび補間法によって可聴閾値に対応する音響強度を得ることを特徴とする請求項1~4のいずれか一項に記載の聴性脳幹反応の自動測定方法。 The sound intensity corresponding to the audible threshold is obtained by sigmoid function fitting and interpolation after normalizing the number of iterations used in each sound intensity in step A6. 1. The method for automatically measuring auditory brainstem response according to claim 1. 各音響強度でのABR記録は、動物実験データまたは臨床データであり、リアルタイム測定によって取得されるかまたはオフラインで記憶されるデータであることと、
前記ABR記録が、
信号増幅と、バンドパスフィルタリングと、
前記ABR波形を収集する時間間隔を調整し、音声刺激が開始した後に対応する時間間隔のABR時間曲線を選択して分析対象とすることと、
バックグラウンドノイズに対応するABR時間曲線を排除することと、
平滑化スプラインフィッティング関数によって、低周波数バックグラウンドノイズを除去することと、からなるグループの少なくとも一つの前処理を経ることと、
を特徴とする請求項1~4のいずれか一項に記載の聴性脳幹反応の自動測定方法。
ABR recordings at each acoustic intensity are animal experimental data or clinical data, data acquired by real-time measurements or stored offline;
The ABR recording is
signal amplification and bandpass filtering;
adjusting the time interval for collecting the ABR waveforms and selecting an ABR time curve for the corresponding time interval after the onset of the audio stimulus for analysis;
eliminating ABR time curves corresponding to background noise;
removing low frequency background noise by a smoothing spline fitting function; and undergoing at least one preprocessing of the group consisting of:
The method for automatically measuring an auditory brainstem response according to any one of claims 1 to 4, characterized by:
請求項1~7のいずれか一項に記載の聴性脳幹反応の自動測定方法に用いられる聴性脳幹反応の自動測定装置であって、
各音響強度の刺激で収集されたABR記録を取得する入力モジュールと、
前記入力モジュールが開始音響強度から測定した音響強度の最小値までの間に、各音響強度の刺激でバッチによって収集されたABR記録を取得するように駆動する制御モジュールと、
適応平均法の演算を実行することによって、ABR信号の検出の条件を満たすまで、反復によって記録の平均回数を段階的に増加させ、ABR信号の検出の条件は、現在の反復時に現在の音響強度に対応するABR記録をグループ化した後、現在の平均回数に基づいて各グループの平均曲線をそれぞれ計算し、グループの間の相互相関関数の最大値に対応するタイムラグを取得し、タイムラグのずれが所定のデータポイント範囲内にあるか否かに基づき、タイムロック特性を有するABR信号が存在するか否かを判断することを含む適応平均法演算モジュールと、
適応平均法の演算時に各音響強度で使用される反復回数を保存する記憶モジュールと、
適応平均法演算モジュールによって出力された結果に基づき、適応平均法の演算時にABR信号が検出されているか否かを判断し、かつ、ABR信号が検出されていると判断した場合、さらに現在の音響強度が設定された測定した音響強度の最小値に達しているか否かを判断し、またはABR信号が検出されていないと判断した場合、さらに連続するP個の音響強度でABR信号が検出されていないか否かを判断するメイン判断モジュールとを含むことと、
測定した音響強度の最小値に達していないと判断した場合、メイン判断モジュールが、制御モジュールに指令を出し、制御モジュールによってさらに入力モジュールが新たな音響強度で対応するABR記録を取得するように駆動し、新たな音響強度は、現在の音響強度から間隔を差し引いた後に得られるものであり、または、測定した音響強度の最小値に達していると判断した場合、メイン判断モジュールは、制御モジュールに指令を出し、制御モジュールによって記憶モジュールから適応平均法の演算時に使用される反復回数を呼び出し、ABR信号を検出するために必要とする音響強度の最小値を可聴閾値とすることと、を特徴とする聴性脳幹反応の自動測定装置。
An automatic auditory brainstem response measuring device for use in the method for automatically measuring an auditory brainstem response according to any one of claims 1 to 7,
an input module for obtaining ABR recordings collected at each acoustic intensity stimulus;
a control module that drives the input module to acquire batch-acquired ABR recordings at each acoustic intensity stimulus from a starting acoustic intensity to the minimum measured acoustic intensity;
By performing an adaptive averaging operation, the number of recording averaging steps is increased by iterations until the condition for ABR signal detection is met, the condition for ABR signal detection being the current sound intensity at the current iteration. After grouping the ABR records corresponding to , calculate the average curve of each group respectively based on the current average times, obtain the time lag corresponding to the maximum value of the cross-correlation function between the groups, and determine if the time lag deviation is an adaptive averaging computation module including determining whether an ABR signal having time-locked characteristics is present based on whether it is within a predetermined range of data points;
a storage module for storing the number of iterations used at each acoustic intensity when computing the adaptive averaging;
Based on the result output by the adaptive averaging calculation module, it is determined whether or not the ABR signal is detected during the adaptive averaging calculation, and if it is determined that the ABR signal is detected, the current sound It is determined whether the intensity has reached a set minimum value of measured sound intensities, or if it is determined that no ABR signal has been detected, further ABR signals have been detected at P consecutive sound intensities. a main determination module for determining whether there is
If it determines that the minimum measured sound intensity has not been reached, the main decision module commands the control module, which further drives the input module to acquire a corresponding ABR record at the new sound intensity. and the new sound intensity is obtained after subtracting the interval from the current sound intensity, or if it determines that the minimum measured sound intensity has been reached, the main decision module tells the control module issuing a command and recalling by the control module from the storage module the number of iterations used in computing the adaptive averaging method, and setting the audible threshold to be the minimum sound intensity required to detect the ABR signal. An automatic measurement device for auditory brainstem responses.
適応平均法の演算時に使用される反復回数に基づいて関数フィッティングおよび補間法によって可聴閾値に対応する音響強度を取得する関数フィッティングモジュールをさらに含むことを特徴とする請求項8に記載の聴性脳幹反応の自動測定装置。 9. The auditory brainstem response of claim 8, further comprising a function fitting module for obtaining the acoustic intensity corresponding to the audible threshold by function fitting and interpolation based on the number of iterations used when computing the adaptive averaging method. automatic measuring device. 前記適応平均法演算モジュールが、制御モジュールにデータを収集する指令を出し、入力モジュールによって入力されるデータタイプが繰り返される単一記録である場合、前記適応平均法演算モジュールが、各音響強度の各回の反復に対応して1回の判断を実行するかまたは複数回の並行判断を実行し、1回の判断または複数回の並行判断のうちの各回の判断時に、現在収集されたABR記録をランダムに二つのグループに分けて平均曲線をそれぞれ計算した後に相互相関関数の最大値に対応するタイムラグを計算し、タイムラグの絶対値とk値とを比較し、タイムラグずれがk個のデータポイント内にあるか否かを判断し、タイムラグの絶対値がk値よりも小さい場合、安定したタイムロック信号が検出されていることを表し、タイムラグの絶対値がk値以上の場合、安定したタイムロック信号が検出されていないことを表すことと、
1回の判断または複数回の並行判断のうちの各回の判断がいずれも、現在の音響強度で安定したタイムロック信号が検出されていることを示す場合、適応平均法演算モジュールは、制御モジュールにデータの収集を停止する指令を出し、現在の音響強度でABR信号が検出されていることを示す戻り情報を出し、かつ記憶モジュールによって現在の反復回数を保存することと、
1回の判断が現在の音響強度で安定したタイムロック信号が検出されていないことを示すか、または複数回の並行判断のうちの各回の判断が必ずしも、現在の音響強度で安定したタイムロック信号が検出されていることを示していない場合、適応平均法演算モジュールは、さらに現在の音響強度での最大反復回数に達しているか否かを判断することと、
最大反復回数に達している場合、適応平均法演算モジュールは、制御モジュールにデータの収集を停止する指令を出し、現在の音響強度でABR信号が検出できないことを示す戻り情報を出し、かつ記憶モジュールによって現在の反復回数を保存することと、
最大反復回数に達していない場合、適応平均法演算モジュールは、さらに制御モジュールに新たなデータを加える指令を出し、かつ元データに新たに収集されたデータを加えた後に現在の音響強度で適応平均法の演算プロセスを繰り返し実行することと、
kが、予め設定された固定値であり、または、全てのデータポイントを予め設定された比率で計算した後に得られる値であることと、を特徴とする請求項8または9に記載の聴性脳幹反応の自動測定装置。
If the adaptive averaging module commands the control module to collect data and the data type input by the input module is a repeated single record, the adaptive averaging module will performing a single decision or performing multiple parallel decisions corresponding to the iteration of , and randomizing the currently collected ABR records at each decision of the single decision or multiple parallel decisions. After calculating the average curve for each of the two groups, calculate the time lag corresponding to the maximum value of the cross-correlation function, compare the absolute value of the time lag and the k value, and find that the time lag deviation is within k data points If the absolute value of the time lag is smaller than the k value, it indicates that a stable time lock signal has been detected; if the absolute value of the time lag is greater than the k value, the stable time lock signal is not detected; and
If the single determination or each of the multiple parallel determinations both indicate that a stable time-locked signal is being detected at the current acoustic intensity, the adaptive averaging computation module instructs the control module to issuing a command to stop collecting data, issuing return information indicating that an ABR signal is being detected at the current acoustic intensity, and storing the current number of iterations by the storage module;
Either a single determination indicates that a stable time-locked signal has not been detected at the current sound intensity, or each of multiple parallel determinations necessarily indicates a stable time-locked signal at the current sound intensity. has not been detected, the adaptive averaging computation module further determines whether the maximum number of iterations at the current acoustic intensity has been reached;
If the maximum number of iterations has been reached, the adaptive averaging computation module commands the control module to stop collecting data, issues return information indicating that no ABR signal can be detected at the current acoustic intensity, and the storage module saving the current iteration number by
If the maximum number of iterations has not been reached, the adaptive averaging calculation module further commands the control module to add new data, and after adding the newly acquired data to the original data, the adaptive averaging is performed at the current acoustic intensity. repeatedly performing the modulus computation process;
10. Auditory brainstem according to claim 8 or 9, characterized in that k is a preset fixed value or a value obtained after calculating all data points with a preset ratio. Automatic measurement device for reaction.
前記適応平均法演算モジュールが、制御モジュールにデータを収集する指令を出し、入力モジュールによって入力されるデータタイプが複数回の平均記録である場合、適応平均法演算モジュールは、今回の反復で新たに加えられる平均曲線avgA加を、前回の反復の平均曲線avgA旧、と組み合わせ、重みに従って現在の平均曲線avgA新を計算することと、
Figure 2020220564000002
ここでは、avgA新、avgA旧、avgA加は、それぞれQグループのABR測定平均データであり、Mは、現在の反復回数であり、適応平均法演算モジュールは、Qグループの現在の平均曲線avgA新に対して、二グループ間毎の相互相関関数を計算し、かつ相互相関関数の最大値に対応するQグループのタイムラグを取得することと、
前記Q回の並行判断において、各グループのタイムラグの絶対値をそれぞれk値と比較し、タイムラグずれがk個のデータポイント内にあるか否かを判断し、タイムラグの絶対値がkよりも小さい場合、タイムラグずれがk個のデータポイント内にあることを意味し、安定したタイムロック信号が検出されていることを表し、タイムラグの絶対値がk値以上の場合、タイムラグずれがk個のデータポイント内にないことを意味し、安定したタイムロック信号が検出されていないことを表すことと、
前記適応平均法演算モジュールが、Q回の並行判断のうちの各回にいずれも安定したタイムロック信号が検出されているか否かを判断し、Q回のうちの各回にいずれも安定したタイムロック信号が検出されている場合、制御モジュールにデータの収集を停止する指令を出し、現在の音響強度でABR信号が検出されていることを示す戻り情報を出し、かつ記憶モジュールによって現在の反復回数を保存することと、
前記Q回のうちの各回に必ずしも安定したタイムロック信号が検出されていない場合、さらに現在の音響強度での最大反復回数に達しているか否かを判断し、最大反復回数に達している場合、制御モジュールにデータの収集を停止する指令を出し、現在の音響強度でABR信号が検出されていることを示す戻り情報を出し、かつ記憶モジュールによって現在の反復回数を保存することと、
最大反復回数に達していない場合、さらに現在収集されたABR測定平均データに新たに収集されたABR測定平均データを加え、現在の音響強度で適応平均法の演算プロセスを繰り返し実行することと、
kが、予め設定された固定値であり、または、全てのデータポイントを予め設定された比率で計算した後に得られる値であることと、を特徴とする請求項8または9に記載の聴性脳幹反応の自動測定装置。
If the adaptive averaging computation module commands the control module to collect data and the data type input by the input module is multiple average records, the adaptive averaging computation module will newly in this iteration combining the added average curve avgA addition with the average curve avgA old of the previous iteration and calculating the current average curve avgA new according to the weight;
Figure 2020220564000002
Here, avgA new, avgA old, and avgA addition are the ABR measurement average data of the Q group respectively, M is the current iteration number, and the adaptive averaging calculation module is the current average curve avgA new of the Q group. , calculating the cross-correlation function between each two groups, and obtaining the time lag of the Q group corresponding to the maximum value of the cross-correlation function;
In the Q parallel determinations, the absolute value of the time lag of each group is respectively compared with the k value, and it is determined whether the time lag deviation is within k data points, and the absolute value of the time lag is less than k. If the time lag is within k data points, it means that a stable time lock signal has been detected. means not within a point, indicating that no stable timelock signal has been detected;
The adaptive averaging calculation module determines whether a stable time-locked signal is detected in each of Q parallel determinations, and determines whether a stable time-locked signal is detected in each of Q parallel determinations. is detected, command the control module to stop collecting data, provide return information indicating that an ABR signal is being detected at the current acoustic intensity, and store the current number of iterations by the storage module. and
If a stable time lock signal is not necessarily detected in each of the Q times, further determine whether the maximum number of repetitions at the current sound intensity has been reached, and if the maximum number of repetitions has been reached, commanding the control module to stop collecting data, issuing return information indicating that an ABR signal is being detected at the current acoustic intensity, and storing the current number of iterations with the storage module;
if the maximum number of iterations has not been reached, adding newly collected ABR measured average data to the currently collected ABR measured average data, and repeatedly executing the adaptive averaging calculation process at the current acoustic intensity;
10. Auditory brainstem according to claim 8 or 9, characterized in that k is a preset fixed value or a value obtained after calculating all data points with a preset ratio. Automatic measurement device for reaction.
JP2021562859A 2019-04-30 2019-09-17 Apparatus and method for automatic measurement of auditory brainstem response Pending JP2022534851A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2023178139A JP2023179702A (en) 2019-04-30 2023-10-16 Automatic measuring device and method for auditory brainstem response

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201910362550.0A CN109998539B (en) 2019-04-30 2019-04-30 Automatic testing arrangement of auditory brainstem response
CN201910362550.0 2019-04-30
PCT/CN2019/106154 WO2020220564A1 (en) 2019-04-30 2019-09-17 Auditory brainstem response automatic test apparatus and method

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2023178139A Division JP2023179702A (en) 2019-04-30 2023-10-16 Automatic measuring device and method for auditory brainstem response

Publications (2)

Publication Number Publication Date
JP2022534851A JP2022534851A (en) 2022-08-04
JPWO2020220564A5 true JPWO2020220564A5 (en) 2023-04-24

Family

ID=67175413

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2021562859A Pending JP2022534851A (en) 2019-04-30 2019-09-17 Apparatus and method for automatic measurement of auditory brainstem response
JP2023178139A Pending JP2023179702A (en) 2019-04-30 2023-10-16 Automatic measuring device and method for auditory brainstem response

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2023178139A Pending JP2023179702A (en) 2019-04-30 2023-10-16 Automatic measuring device and method for auditory brainstem response

Country Status (5)

Country Link
US (1) US20220183631A1 (en)
JP (2) JP2022534851A (en)
CN (1) CN109998539B (en)
DE (1) DE112019007067T5 (en)
WO (1) WO2020220564A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109998539B (en) * 2019-04-30 2020-05-12 上海交通大学医学院附属第九人民医院 Automatic testing arrangement of auditory brainstem response

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5697379A (en) * 1995-06-21 1997-12-16 Neely; Stephen T. Method and apparatus for objective and automated analysis of auditory brainstem response to determine hearing capacity
DE19548982A1 (en) * 1995-12-28 1997-07-03 Pilot Blankenfelde Medizinisch Process for automatic hearing threshold determination, especially in newborns and small children
US6200273B1 (en) * 1999-04-26 2001-03-13 House Ear Institute Power-optimized cumulative, sequential statistical method for detection of auditory evoked potentials
US6196977B1 (en) * 1999-04-26 2001-03-06 House Ear Institute Method for detection on auditory evoked potentials using a point optimized variance ratio
CN100341462C (en) * 2005-10-20 2007-10-10 上海交通大学 Double ear alternative stimulated brainstem auditory evoked potential(BAEP) tracing device
US8825149B2 (en) * 2006-05-11 2014-09-02 Northwestern University Systems and methods for measuring complex auditory brainstem response
KR100839109B1 (en) * 2006-09-20 2008-06-19 [주]이어로직코리아 The Method and Device for Objective Automated Audiometry
CN106725515A (en) * 2015-11-25 2017-05-31 中国科学院声学研究所 A kind of quick universal newborn hearing screening assessment system
CN109998539B (en) * 2019-04-30 2020-05-12 上海交通大学医学院附属第九人民医院 Automatic testing arrangement of auditory brainstem response

Similar Documents

Publication Publication Date Title
Sadhukhan et al. R-peak detection algorithm for ECG using double difference and RR interval processing
JP5019457B2 (en) Detection of rapid changes in blood pressure
CN110037668B (en) System for judging age, health state and malignant arrhythmia identification by combining pulse signal time-space domain with model
US20110201961A1 (en) Morphological clustering and analysis of intracranial pressure pulses (mocaip)
KR20180052943A (en) Atrial fibrillation discriminating device and atrial fibrillation discrimination method using a neural network
WO2013053290A1 (en) Device and method for improving identification accuracy of fetal heart rate deceleration
CN110123304B (en) Dynamic electrocardio noise filtering method based on multi-template matching and correlation coefficient matrix
CN107622259B (en) T wave detection method, electrocardiogram data analysis method and device
CN105686846A (en) Heart sound signal characteristic automatic extraction method
KR20130085632A (en) R-peak detection method of ecg signal using adaptive local threshold
EP3892190A1 (en) Method for analyzing electrocardiography signal
Bruun et al. Automatic atrial fibrillation detection: A novel approach using discrete wavelet transform and heart rate variability
EP1372473A1 (en) Apparatus and method for detection and quantification of oscillatory signals
AU2002248091A1 (en) Apparatus and method for detection and quantification of oscillatory signals
CN111839494A (en) Heart rate monitoring method and system
CN107753012B (en) Magnetocardiogram signal quality evaluation method, system and server
Zhou et al. A hybrid expert system for individualized quantification of electrical status epilepticus during sleep using biogeography-based optimization
JP6957011B2 (en) Sleep stage determination device, sleep stage determination method and program
Wan et al. A novel atrial fibrillation automatic detection algorithm based on ensemble learning and multi-feature discrimination
CN114052741A (en) Electrocardiogram recognition method based on multi-scale autoregressive model
JPWO2020220564A5 (en)
JP2023179702A (en) Automatic measuring device and method for auditory brainstem response
CN109846473B (en) Method for detecting single-lead 10-second electrocardiogram noise interference degree
JP2009507260A (en) Signal recognition method using a microcontroller
Moukadem et al. Phonocardiogram signal processing module for auto-diagnosis and telemedicine applications