Otoacoustic emission signal detection method, device, electronic device and storage medium

Abstract

The present invention relates to an otoacoustic emission signal detection method, device, electronic device and storage medium. The method comprises: obtaining otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweep. An online whole-segment signal artifact removal algorithm is called to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal respectively, and the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal are respectively subjected to filtering tracking and two-tone suppression superposition to obtain a re-collected signal. An offline segmented signal artifact removal algorithm is called to perform artifact removal on the re-collected signal, and the re-collected signal processed by the offline segmented signal artifact removal algorithm is subjected to segmented signal stacking averaging operation and filtering tracking to extract the re-collected signal of the stimulation frequency otoacoustic emission signal. The re-collected signal is processed by FFT operation to obtain an amplitude signal and a phase signal corresponding to the otoacoustic emission signal with high frequency resolution.

Description

Technical Field

[0001] The present invention relates to the technical field of otoacoustic emissions, and in particular to an otoacoustic emission signal detection method, device, electronic equipment and storage medium. Background Art

[0002] Otoacoustic emission is a weak sound signal generated by the active mechanism of the inner ear cochlea, which is conducted through the middle ear and released into the external auditory canal and can be captured by a high-sensitivity microphone. Otoacoustic emission signals can be divided into two categories according to the cause of their generation: spontaneous otoacoustic emission and evoked otoacoustic emission. Evoked otoacoustic emission can be divided into distortion product otoacoustic emission, transient evoked otoacoustic emission and stimulation frequency otoacoustic emission according to the type of evoked sound. Among them, distortion product otoacoustic emission (DPOAE) is an otoacoustic emission signal recorded in the external auditory canal when two pure tones of different frequencies (F1 and F2) stimulate the cochlea at the same time, which is induced by the cochlea to actively generate the otoacoustic emission signal. DPOAE refers to an otoacoustic emission signal with a frequency of 2F1-F2 that is different from the frequency of the pure tones of stimulation (F1 and F2) in the evoked response. Stimulation frequency otoacoustic emission (SFOAE) is a single-frequency sound signal recorded in the external auditory canal that is the same as the test tone frequency and is actively emitted by the inner ear cochlea.

[0003] At present, the traditional otoacoustic emission (OAEs) signal acquisition paradigm usually uses pure tone as the test sound. When collecting OAEs signals of multiple frequencies, it is necessary to modify the frequency of the test sound and repeat the paradigm. When OAEs with high frequency resolution need to be collected, this method becomes more time-consuming. For example, in the prior art, a hearing sensitivity detection system based on stimulation frequency otoacoustic emissions is proposed, which integrates the function of measuring auditory frequency sensitivity based on SFOAE suppression tuning curve, but the technology does not specifically describe the rapid detection of multiple otoacoustic emission signals with high frequency resolution. In addition, a stimulation frequency otoacoustic emission tuning curve detection and calibration system is also proposed. The technology determines the SFOAE suppression tuning curve of the subject based on the up and down tracking method. The specific method is: for a specific suppression sound frequency, the suppression sound intensity is continuously increased or decreased until the measured SFOAE intensity reaches the suppression criterion intensity and stops. The current suppression sound intensity is the threshold intensity extracted at the suppression sound frequency. The advantage of this technology is that the threshold strength at each suppression frequency is obtained by the computer searching up and down according to the test results and finally converging, but this technology can only test the results at a limited number of suppression sound frequency points, and the frequency resolution is highly limited. It is time-consuming to test the results of higher frequency resolutions, and the test efficiency is low.

[0004] In summary, existing otoacoustic emission signal detection methods are time-consuming and have low detection efficiency when detecting otoacoustic emission signals with high frequency resolution. Summary of the invention

[0005] Based on this, it is necessary to provide a method, device, electronic device and storage medium for detecting otoacoustic emission signals with high frequency resolution, which has short detection time and high detection efficiency, in order to address the above technical problems.

[0006] The present invention provides a method for detecting otoacoustic emission signals, the method comprising:

[0007] Acquiring otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweep, wherein the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals;

[0008] Calling the online whole-segment signal artifact removal algorithm to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and performing filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal, respectively, to obtain a first acquisition signal and a second acquisition signal;

[0009] Calling an offline segmented signal artifact removal algorithm to remove artifacts from the first and second acquisition signals, and performing a segmented signal stacking average operation on the first and second acquisition signals processed by the offline segmented signal artifact removal algorithm, so as to extract a third acquisition signal by filtering and tracking the second acquisition signal after the stacking average operation;

[0010] The first sampled signal and the third sampled signal after the stacked average operation are processed by FFT operation to obtain amplitude signals and phase signals corresponding to the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal with high frequency resolution, respectively;

[0011] Among them, the first acquisition signal is the acquisition signal of the distortion product otoacoustic emission signal, the second acquisition signal includes the acquisition signal of the stimulation frequency otoacoustic emission signal and the stimulation frequency otoacoustic emission signal induced by the inhibitory sound, and the third acquisition signal is the acquisition signal of the stimulation frequency otoacoustic emission signal.

[0012] In one embodiment, the step of acquiring otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweeping comprises:

[0013] Obtaining basic parameters of the frequency sweep, the basic parameters at least including a frequency sweep range, a frequency sweep mode, and a frequency sweep rate, and determining the relationship between the instantaneous frequency of the frequency sweep tone and the frequency sweep time according to the frequency sweep rate to obtain an instantaneous frequency expression of the logarithmic frequency sweep;

[0014] According to the derivative relationship between the instantaneous frequency and the phase of the swept frequency tone, an instantaneous phase expression is obtained, and based on the instantaneous phase expression, a time domain expression of the swept frequency tone is calculated;

[0015] Among them, the instantaneous phase expression is: The time domain expression is

[0016] In the formula, is the instantaneous phase of the swept tone at time t, p(t) is the time domain of the swept tone at time t, and A(t) is the amplitude of the swept tone at time t.

[0017] In one embodiment, the calling of the online whole signal artifact removal algorithm performs artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal respectively, and performs filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after the artifact removal respectively, including:

[0018] Acquire a first sweep frequency test paradigm of the distortion product otoacoustic emission signal, wherein the first sweep frequency test paradigm includes a first single-segment test sound paradigm and a second single-segment test sound paradigm respectively corresponding to a pair of sweep frequency stimulation sounds having a fixed ratio of instantaneous frequencies;

[0019] The online whole-segment signal artifact removal algorithm is called to perform filtering and tracking on the first frequency sweep test paradigm to remove the stimulus sound echo artifact in the first frequency sweep test paradigm and obtain the first echo signal.

[0020] In one embodiment, the calling of the online whole signal artifact removal algorithm performs artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, respectively, and performs filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after the artifact removal, respectively, and further includes:

[0021] Acquire a second frequency sweep test paradigm of the stimulation frequency otoacoustic emission signal, wherein the second frequency sweep test paradigm includes at least four test sound paradigms based on two-tone suppression, namely, a first segmented test sound paradigm and a second segmented test sound paradigm corresponding to stimulation sounds with opposite phases, and a third segmented test sound paradigm and a fourth segmented test sound paradigm corresponding to stimulation sounds with opposite phases and suppression sounds with the same phase;

[0022] Calling the online whole-segment signal artifact removal algorithm to perform dual-tone suppression superposition on the second frequency sweep test paradigm to remove the echo signal artifacts of the stimulation sound and the suppression sound in the second frequency sweep test paradigm, and obtaining the second echo signal;

[0023] In one embodiment, the method further comprises:

[0024] When the recovered signal reaches the preset stacking average, the reference signal corresponding to each recovered signal and the power of each reference signal are calculated to obtain the median and interquartile range of all powers;

[0025] When the power is not within the set range, the recovery signal corresponding to the power not within the set range is marked as an outlier, and the entire recovery signal corresponding to the outlier is deleted;

[0026] The setting range is (μ-1.5σ, μ+1.5σ), where μ is the median and σ is the interquartile range.

[0027] In one embodiment, the calling of the offline segmented signal artifact removal algorithm to remove artifacts from the first and second acquisition signals, and performing a segmented signal stacking average operation on the first and second acquisition signals processed by the offline segmented signal artifact removal algorithm, so as to extract the third acquisition signal by filtering and tracking the second acquisition signal after the stacking average operation, includes:

[0028] Calling the offline segmented signal artifact removal algorithm to remove short-term noise in the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and obtaining the first and second recollection signals without artifact interference by superposition and averaging operation of the segmented signals;

[0029] The second echo signal without artifact interference is filtered and tracked to filter out the otoacoustic emission signal of the stimulation frequency induced by the suppression sound in the second echo signal, so as to obtain the third echo signal.

[0030] In one embodiment, the method further comprises:

[0031] Obtaining instantaneous filtering frequencies of the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and determining pole frequency, zero frequency and instantaneous gain;

[0032] Based on the pole frequency, the zero frequency and the instantaneous gain, an expression of a filter is obtained to construct a tracking filter for filtering and tracking otoacoustic emission signals;

[0033] The pole frequency is the instantaneous filtering frequency corresponding to the retained otoacoustic emission signal, and the zero frequency is the instantaneous filtering frequency corresponding to the otoacoustic emission signal to be filtered out.

[0034] The present invention also provides an otoacoustic emission signal detection device, the device comprising:

[0035] A signal acquisition module, used to acquire otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweep, wherein the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals;

[0036] The first algorithm processing module is used to call the online whole-segment signal artifact removal algorithm to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and perform filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after the artifact removal, to obtain a first recollection signal and a second recollection signal respectively;

[0037] The second algorithm processing module is used to call the offline segmented signal artifact removal algorithm to remove artifacts from the first and second acquisition signals, and perform segmented signal stacking and averaging operations on the first and second acquisition signals processed by the offline segmented signal artifact removal algorithm, so as to extract the third acquisition signal by filtering and tracking the second acquisition signal after the stacking and averaging operations;

[0038] An FFT operation module is used to process the first acquisition signal and the third acquisition signal after the stacked average operation by FFT operation to obtain amplitude signals and phase signals corresponding to the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal with high frequency resolution, respectively;

[0039] Among them, the first acquisition signal is the acquisition signal of the distortion product otoacoustic emission signal, the second acquisition signal includes the acquisition signal of the stimulation frequency otoacoustic emission signal and the stimulation frequency otoacoustic emission signal induced by the inhibitory sound, and the third acquisition signal is the acquisition signal of the stimulation frequency otoacoustic emission signal.

[0040] The present invention also provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor implements any one of the above-described methods for detecting otoacoustic emission signals when executing the computer program.

[0041] The present invention also provides a computer storage medium storing a computer program, wherein when the computer program is executed by a processor, the otoacoustic emission signal detection method as described above is implemented.

[0042] The above-mentioned otoacoustic emission signal detection method, device, electronic device and storage medium obtain distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweeping, and call the online whole-segment signal artifact removal algorithm to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal respectively, and filter and track the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal and perform two-tone suppression superposition to obtain the corresponding acquisition signal, that is, the first acquisition signal and the second acquisition signal. Subsequently, call the offline segmented signal artifact removal algorithm to perform artifact removal on the acquisition signal, and perform segmented signal stacking and averaging operation on the acquisition signal processed by the offline segmented signal artifact removal algorithm, and then further filter and track the second acquisition signal after the stacking and averaging operation to extract the acquisition signal of the stimulation frequency otoacoustic emission signal. Finally, the collected signals of the distortion product otoacoustic emission signals after the stacked average operation and the collected signals of the extracted stimulation frequency otoacoustic emission signals are processed by FFT operation to obtain the amplitude signal and phase signal corresponding to the distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals with high frequency resolution. This method uses exponential frequency sweeping to obtain OAEs signals (otoacoustic emission signals) at multiple frequencies in one frequency sweeping test, and by jointly applying the artifact removal algorithm of the online whole signal and the artifact removal algorithm of the offline segmented signal, it realizes the rapid detection of multiple otoacoustic emission signals with high frequency resolution, reduces the detection time, and improves the detection efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0044] Figure 1 This is a schematic diagram of a flow chart of the otoacoustic emission signal detection method provided by the present invention;

[0045] Figure 2 A schematic diagram of a process of extracting OAEs spectrum with high frequency resolution in a method for detecting otoacoustic emission signals in a specific embodiment of the present invention;

[0046] Figure 3 A schematic flow chart of an online artifact (AR) removal algorithm of an otoacoustic emission signal detection method in a specific embodiment of the present invention;

[0047] Figure 4A schematic flow chart of an offline artifact (AR) removal algorithm of a method for detecting otoacoustic emission signals in a specific embodiment of the present invention;

[0048] Figure 5 The second flow chart of the otoacoustic emission signal detection method provided by the present invention;

[0049] Figure 6 The third flow chart of the otoacoustic emission signal detection method provided by the present invention;

[0050] Figure 7 The fourth flow chart of the otoacoustic emission signal detection method provided by the present invention;

[0051] Figure 8 The fifth flow chart of the otoacoustic emission signal detection method provided by the present invention;

[0052] Fig. 9 The sixth flow chart of the otoacoustic emission signal detection method provided by the present invention;

[0053] Fig.10 Schematic diagram of the seventh flow chart of the otoacoustic emission signal detection method provided by the present invention;

[0054] Fig.11 A schematic diagram of the structure of the otoacoustic emission signal detection device provided by the present invention;

[0055] Fig.12 This is a diagram of the internal structure of the computer device provided by the present invention. DETAILED DESCRIPTION

[0056] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0057] Combine the following Figure 1-Figure 12 The present invention describes an otoacoustic emission signal detection method, device, electronic device and storage medium.

[0058] like Figure 1 As shown, in one embodiment, a method for detecting otoacoustic emission signals comprises the following steps:

[0059] Step S110, obtaining otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweep, wherein the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals.

[0060] Specifically, the server acquires otoacoustic emission signals at multiple frequencies in a single frequency sweep test through exponential frequency sweep, and the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals.

[0061] Combination Figure 2 As shown, in a specific embodiment, the otoacoustic emission signal detection method provided by the present invention uses a swept tone whose frequency changes continuously with time instead of a traditional pure tone as a test sound, and adopts an exponential frequency sweep in the hope of obtaining OAEs signals (otoacoustic emission signals) at multiple frequencies in one swept frequency test.

[0062] Step S120, calling the online whole signal artifact removal algorithm to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and performing filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal, to obtain the first re-collected signal and the second re-collected signal, respectively.

[0063] Specifically, the server calls the online whole signal artifact removal algorithm to remove artifacts from the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and performs filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal, to obtain corresponding re-acquisition signals, namely the first re-acquisition signal and the second re-acquisition signal.

[0064] Combination Figure 2 As shown, in a specific embodiment, the method for detecting otoacoustic emission signals provided by the present invention comprises the following steps for extracting DPOAEs (distortion product otoacoustic emission signals): 1) online whole-segment signal artifact removal; 2) tracking and filtering of swept frequency signals; 3) offline segmented signal artifact removal; 4) stacking and averaging of segmented signals; 5) finally obtaining the amplitude and phase signals of DPOAEs signals with high frequency resolution through FFT operation. The steps for extracting SFOAEs (stimulation frequency otoacoustic emission signals) are: 1) online whole-segment signal artifact removal; 2) dual-tone suppression and superposition; 3) offline segmented signal artifact removal; 4) stacking and averaging of segmented signals, 5) tracking and filtering of swept frequency signals; 6) finally obtaining the amplitude and phase signals of SFOAE signals with high frequency resolution through FFT operation.

[0065] In this embodiment, the synthesis of the sweep tone is generally divided into three steps. The first step is to determine the basic parameters of the sweep tone, namely, the sweep tone range, the sweep tone mode, and the sweep tone rate. Assume that the total duration of the sweep tone is T, and the sweep tone range is the frequency f0 at the start time and the frequency f at the end time. T , based on the exponential sweep method, the sweep rate b is obtainede , and then determine the instantaneous frequency f of the swept tone e (t), f e (t) is the instantaneous rate at time t. The third step is to find the instantaneous phase expression based on the derivative relationship between frequency and phase: Then the time domain expression of the swept frequency tone is calculated as

[0066]

[0067] In the formula, is the instantaneous phase of the swept tone at time t, p(t) is the time domain of the swept tone at time t, and A(t) is the amplitude of the swept tone at time t

[0068] In this embodiment, the sweep frequency test paradigm of DPOAEs is a single-segment test sound paradigm including a pair of sweep frequency stimulation sounds p1 and p2, and the instantaneous frequencies of p1 and p2 are in a fixed ratio (usually f2 / f1 is 1.22), and the instantaneous frequency of the induced DPOAEs signal is usually f DPOAE (t) = 2f1(t) - f2(t), therefore, in the echo signal r, the DPOAEs signal can be obtained by filtering and removing the stimulus echo artifacts r1 and r2. In contrast, the frequency sweep test paradigm of SFOAEs is more complex, that is, a four-segment test sound paradigm based on two-tone suppression. The four-segment test sound is represented as p A 、p B 、p C and p D , segments A and B only contain stimuli with opposite phases. and The C and D segments contain stimulating sounds with opposite phases and suppressing sounds with the same phase, respectively. Both the stimulation and inhibition sounds in the paradigm are swept frequency tones, and their frequencies are also in a fixed ratio.

[0069] Combination Figure 3 As shown, in this embodiment, during the online artifact removal process of DPOAEs and SFOAEs, first, when the acquisition signal is collected after the stimulation paradigm is applied, the online artifact rejection (AR) algorithm will be activated to delete the entire acquisition signal containing long-term interference. Secondly, after collecting a preset number of acquisition signals through the online AR, DPOAEs and SFOAEs are respectively superimposed by tracking filtering and two-tone suppression of the SFOAEs signal expression to obtain their own original OAEs time domain signals.

[0070] In addition, online AR is performed during the signal acquisition process, and the entire segment of the acquisition signal that is determined to be an artifact is removed. This process is activated when the acquisition signal reaches the preset stacking average number N. First, each acquisition signal r is calculated.i Corresponding reference signal To temporarily exclude the energy of the recovered test sound, calculate the power P of each reference signal i , and calculate all the powers P i The median μ and interquartile range σ of the sampled signal r i Mark as outlier and delete the entire segment until all N collected signals are not judged as outliers.

[0071] Step S130, calling the offline segmented signal artifact removal algorithm to remove artifacts from the first and second acquisition signals, and performing a segmented signal stacking average operation on the first and second acquisition signals processed by the offline segmented signal artifact removal algorithm, so as to extract the third acquisition signal by filtering and tracking the second acquisition signal after the stacking average operation.

[0072] Among them, the first collected signal is the collected signal of the distortion product otoacoustic emission signal, the second collected signal includes the collected signal of the stimulation frequency otoacoustic emission signal and the stimulation frequency otoacoustic emission signal induced by the inhibitory sound, and the third collected signal is the collected signal of the stimulation frequency otoacoustic emission signal.

[0073] Specifically, the server calls the offline segmented signal artifact removal algorithm to collect the distortion product otoacoustic emission signal r DP (t) and the echo signal r including the stimulation frequency otoacoustic emission signal SF (t) and the stimulation frequency otoacoustic emission signal r induced by the inhibitory sound SF′ The artifact removal of the second echo signal of the first echo signal (t) is performed, and the segmented signal stacking and averaging operation is performed on the first echo signal and the second echo signal after the offline segmented signal artifact removal algorithm, and then the second echo signal after the stacking and averaging operation is further filtered and tracked to extract the echo signal r of the stimulation frequency otoacoustic emission signal. SF (t).

[0074] Combination Figure 4 As shown, in a specific embodiment, the otoacoustic emission signal detection method provided by the present invention further removes the short-term noise in the OAEs by an offline AR algorithm, and obtains an OAEs signal with relatively no artifact interference by segmented superposition and averaging. SFOAEs are filtered once more than DPOAEs after this step to filter out SFOAEs that may be induced by suppressive sound. Finally, the obtained OAEs time domain signal is subjected to time-frequency analysis to obtain its spectrum.

[0075] In this embodiment, offline AR is a segmented artifact removal algorithm on N original OAEs signals. After being extracted, they are first divided into M frames Then the signal fragments of each frame are treated as a group for artifact removal. The removal method is similar to the process of the online AR algorithm. The power of the signal fragment is used as the basis for judging whether it is an outlier in the group. The number of valid signals retained in each frame group is recorded. j Then, the average value of the group signal of each frame is calculated according to the effective number and effective fragment, and finally the averaged frame signals are spliced ​​to complete the offline AR algorithm. This offline AR algorithm mainly removes the short-term impulse noise introduced by electrical interference.

[0076] Step S140, processing the first and third echoed signals respectively through FFT operation to obtain amplitude signals and phase signals corresponding to the distortion product otoacoustic emission signals and the stimulation frequency otoacoustic emission signals with high frequency resolution.

[0077] Specifically, the server processes the sampled signal of the distortion product otoacoustic emission signal in step S130 (i.e., the first sampled signal after the stacked averaging operation) and the sampled signal of the stimulation frequency otoacoustic emission signal (i.e., the third sampled signal) through FFT operation to obtain amplitude signals and phase signals corresponding to the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal with high frequency resolution, respectively.

[0078] The above-mentioned otoacoustic emission signal detection method obtains distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals at multiple frequencies in a single frequency sweep test by means of exponential frequency sweep, and calls the online whole-segment signal artifact removal algorithm to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal respectively, and performs filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal, respectively, to obtain the corresponding recollection signal. Subsequently, the offline segmented signal artifact removal algorithm is called to perform artifact removal on the recollection signal, and the recollection signal processed by the offline segmented signal artifact removal algorithm is subjected to segmented signal stacking and averaging operation, and then the second recollection signal after the stacking and averaging operation is further filtered and tracked to extract the recollection signal of the stimulation frequency otoacoustic emission signal. Finally, the recollection signal of the distortion product otoacoustic emission signal and the recollection signal of the stimulation frequency otoacoustic emission signal are processed respectively by FFT operation to obtain the amplitude signal and phase signal corresponding to the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal with high frequency resolution. This method uses exponential frequency sweeping, which can obtain OAEs signals (otoacoustic emission signals) at multiple frequencies in a single frequency sweep test. By jointly applying the online artifact removal algorithm for the entire signal and the offline artifact removal algorithm for the segmented signal, it can achieve rapid detection of multiple otoacoustic emission signals with high frequency resolution, reduce detection time, and improve detection efficiency.

[0079] like Figure 5 As shown, in one embodiment, the otoacoustic emission signal detection method provided by the present invention obtains otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweep, and previously includes the following steps:

[0080] Step S510, obtaining basic parameters of the frequency sweep, which at least include the frequency sweep range, the frequency sweep mode and the frequency sweep rate, and determining the relationship between the instantaneous frequency of the sweep tone and the frequency sweep time according to the frequency sweep rate to obtain the instantaneous frequency expression of the logarithmic frequency sweep.

[0081] Specifically, the server obtains basic parameters of the frequency sweep, which at least include the frequency sweep range, the frequency sweep mode and the frequency sweep rate, and determines the relationship between the instantaneous frequency of the sweep tone and the frequency sweep time according to the frequency sweep rate to obtain the instantaneous frequency expression of the logarithmic frequency sweep.

[0082] Step S520, obtaining an instantaneous phase expression according to the derivative relationship between the instantaneous frequency and the phase of the swept-frequency tone, and calculating a time-domain expression of the swept-frequency tone based on the instantaneous phase expression.

[0083] Specifically, the server obtains the instantaneous phase expression according to the derivative relationship between the instantaneous frequency and the swept frequency tone phase obtained in step S510, that is, And calculate the time domain expression of the swept tone based on the instantaneous phase expression

[0084] In the formula, is the instantaneous phase of the swept tone at time t, p(t) is the time domain of the swept tone at time t, and A(t) is the amplitude of the swept tone at time t.

[0085] like Figure 6 As shown, in one embodiment, the otoacoustic emission signal detection method provided by the present invention calls an online whole-segment signal artifact removal algorithm to remove artifacts from the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and performs filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal, and specifically includes the following steps:

[0086] Step S122, obtaining a first sweep frequency test paradigm of the distortion product otoacoustic emission signal, wherein the first sweep frequency test paradigm includes a first single-segment test sound paradigm and a second single-segment test sound paradigm respectively corresponding to a pair of sweep frequency stimulation sounds with a fixed ratio of instantaneous frequencies.

[0087] Specifically, the server obtains a first sweep frequency test paradigm of the distortion product otoacoustic emission signal, which includes a first single-segment test sound paradigm and a second single-segment test sound paradigm corresponding to a pair of sweep frequency stimulation sounds with fixed instantaneous frequencies. For example, the instantaneous frequencies of p1 and p2 are in a fixed ratio (usually f2 / f1 is 1.22).

[0088] Step S124, calling the online whole signal artifact removal algorithm to perform filtering and tracking on the first frequency sweep test paradigm to remove the stimulus sound echo artifacts in the first frequency sweep test paradigm and obtain a first echo signal.

[0089] Specifically, the server calls the online whole signal artifact removal algorithm to perform filtering and tracking on the first frequency sweep test paradigm to remove the stimulus sound echo artifact in the first frequency sweep test paradigm and obtain the echo signal of the distortion product otoacoustic emission signal.

[0090] like Figure 7 As shown, in one embodiment, the otoacoustic emission signal detection method provided by the present invention calls an online whole-segment signal artifact removal algorithm to remove artifacts from the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and performs filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal, and specifically further includes the following steps:

[0091] Step S126, obtaining a second frequency sweep test paradigm of the stimulation frequency otoacoustic emission signal, the second frequency sweep test paradigm including at least four test sound paradigms based on dual-tone suppression, namely, a first segmented test sound paradigm and a second segmented test sound paradigm corresponding to stimulation sounds with opposite phases, and a third segmented test sound paradigm and a fourth segmented test sound paradigm corresponding to stimulation sounds with opposite phases and suppression sounds with the same phases.

[0092] Specifically, the server obtains a second frequency sweep test paradigm of the stimulation frequency otoacoustic emission signal, wherein the second frequency sweep test paradigm includes at least four test sound paradigms based on dual-tone suppression, namely, a first segmented test sound paradigm and a second segmented test sound paradigm corresponding to stimulation sounds with opposite phases, and a third segmented test sound paradigm and a fourth segmented test sound paradigm corresponding to stimulation sounds with opposite phases and suppression sounds with the same phase, for example, A 、p B 、p C and p D , segments A and B only contain stimuli with opposite phases. and The C and D segments contain stimulating sounds with opposite phases and suppressing sounds with the same phase, respectively.

[0093] Step S128, calling the online whole-segment signal artifact removal algorithm to perform dual-tone suppression superposition on the second frequency sweep test paradigm to remove the echo signal artifacts of the stimulation sound and the suppression sound in the second frequency sweep test paradigm to obtain a second echo signal.

[0094] Specifically, the server calls the online whole signal artifact removal algorithm to perform dual-tone suppression superposition on the second sweep frequency test paradigm to remove the echo signal artifacts of the stimulation sound and the suppression sound in the second sweep frequency test paradigm, and obtains the echo signal of the stimulation frequency otoacoustic emission signal, that is, the second echo signal.

[0095] like Figure 8 As shown, in one embodiment, the otoacoustic emission signal detection method provided by the present invention further includes the following steps:

[0096] Step S810, when the recovered signal reaches a preset stacking average, the reference signal corresponding to each recovered signal and the power of each reference signal are calculated to obtain the median and interquartile range of all powers.

[0097] Specifically, when the recovered signal reaches a preset stacking average, the server calculates the reference signal corresponding to each recovered signal and the power of each reference signal to obtain the median μ and the interquartile range σ of all powers.

[0098] Step S820: When the power is not within the set range, the recovery signal corresponding to the power not within the set range is marked as an outlier, and the entire recovery signal corresponding to the outlier is deleted.

[0099] Specifically, when the power is not within the set range, the server will mark the acquisition signal corresponding to the power not within the set range as an outlier, and delete the entire acquisition signal corresponding to the outlier.

[0100] like Fig. 9 As shown, in one embodiment, the otoacoustic emission signal detection method provided by the present invention calls an offline segmented signal artifact removal algorithm to remove artifacts from a first acquisition signal and a second acquisition signal, and performs a segmented signal stacking average operation on the first acquisition signal and the second acquisition signal processed by the offline segmented signal artifact removal algorithm, so as to extract a third acquisition signal by filtering and tracking the second acquisition signal after the stacking average operation, specifically comprising the following steps:

[0101] Step S132, calling an offline segmented signal artifact removal algorithm to remove short-term noise in the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and obtaining a first echo signal and a second echo signal without artifact interference by superposition and averaging of the segmented signals.

[0102] Specifically, the server calls an offline segmented signal artifact removal algorithm to remove short-term noise in the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and obtains the first and second echo signals without artifact interference by superposition and averaging of the segmented signals.

[0103] Step S134, filtering and tracking the second echo signal without artifact interference to filter out the otoacoustic emission signal of the stimulation frequency induced by the suppression sound in the second echo signal, and obtaining a third echo signal.

[0104] Specifically, the server performs filtering and tracking on the second echo signal without artifact interference to filter out the otoacoustic emission signal of the stimulation frequency induced by the suppression sound in the second echo signal to obtain the third echo signal.

[0105] like Fig.10 As shown, in one embodiment, the otoacoustic emission signal detection method provided by the present invention further includes the following steps:

[0106] Step S1010, obtaining instantaneous filtering frequencies of the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and determining the pole frequency, the zero frequency and the instantaneous gain.

[0107] Specifically, the server obtains the instantaneous filtering frequency of the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and determines the pole frequency, the zero frequency and the instantaneous gain.

[0108] The pole frequency is the instantaneous filtering frequency corresponding to the retained otoacoustic emission signal, and the zero frequency is the instantaneous filtering frequency corresponding to the otoacoustic emission signal to be filtered out.

[0109] Step S1020: obtaining a filter expression based on the pole frequency, the zero frequency and the instantaneous gain to construct a tracking filter for filtering and tracking otoacoustic emission signals.

[0110] Specifically, the server obtains an expression of a filter based on the pole frequency, the zero frequency and the instantaneous gain determined in step S1010 to construct a tracking filter for filtering and tracking otoacoustic emission signals.

[0111] In a specific embodiment, the otoacoustic emission signal detection method provided by the present invention, the tracking filter is a filter designed for the swept frequency OAEs signal to filter out interfering sounds whose frequency varies with time. Its essence is a set of dynamic zeros and poles that vary with time. The construction of the tracking filter includes four steps. The first step is to determine the instantaneous filtering frequency, the second step is to set the gain of the filter, and the third step is to determine the expression of the filter based on the gain and frequency; the fourth step is to filter the signal. When filtering, the input signal needs to be converted into a corresponding analytical signal. After filtering, the real part of the filter output signal is taken as the actual output signal. Taking the tracking filtering of SFOAEs as an example, the poles of the filter need to be set to the instantaneous frequency f of the SFOAEs induced by the stimulus sound. SF (t), and the zero point of the filter needs to be set to the instantaneous frequency f of the SFOAEs induced by the suppressed sound that may not be completely eliminated by the two-tone suppression superposition SF′ (t). It should be noted that due to the existence of the filter bandwidth, the artifacts with frequencies within the bandwidth cannot be filtered out. In actual operation, the bandwidth range can be narrowed by reducing the sweep rate to filter out artifacts with frequencies different from the poles as much as possible.

[0112] The otoacoustic emission signal detection device provided by the present invention is described below. The otoacoustic emission signal detection device described below and the otoacoustic emission signal detection method described above can be referred to each other.

[0113] like Fig.11 As shown, in one embodiment, an otoacoustic emission signal detection device includes a signal acquisition module 1110 , a first algorithm processing module 1120 , a second algorithm processing module 1130 and an FFT operation module 1140 .

[0114] The signal acquisition module 1110 is used to acquire otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweeping, where the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals.

[0115] The first algorithm processing module 1120 is used to call the online whole signal artifact removal algorithm to remove artifacts from the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and to perform filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after the artifact removal, so as to obtain the first re-collection signal and the second re-collection signal, respectively.

[0116] The second algorithm processing module 1130 is used to call the offline segmented signal artifact removal algorithm to remove artifacts from the first acquisition signal and the second acquisition signal, and perform segmented signal stacking and averaging operations on the first acquisition signal and the second acquisition signal processed by the offline segmented signal artifact removal algorithm, so as to extract the third acquisition signal by filtering and tracking the second acquisition signal after the stacking and averaging operations.

[0117] The FFT operation module 1140 is used to process the first and third echoed signals after the stacked average operation through FFT operation to obtain amplitude signals and phase signals corresponding to the distortion product otoacoustic emission signals and the stimulation frequency otoacoustic emission signals with high frequency resolution, respectively.

[0118] Among them, the first collected signal is the collected signal of the distortion product otoacoustic emission signal, the second collected signal includes the collected signal of the stimulation frequency otoacoustic emission signal and the stimulation frequency otoacoustic emission signal induced by the inhibitory sound, and the third collected signal is the collected signal of the stimulation frequency otoacoustic emission signal.

[0119] In this embodiment, the otoacoustic emission signal detection device provided by the present invention further includes an expression calculation module for:

[0120] The basic parameters of the frequency sweep are obtained, which at least include the frequency sweep range, the frequency sweep mode and the frequency sweep rate, and the relationship between the instantaneous frequency of the frequency sweep tone and the frequency sweep time is determined according to the frequency sweep rate to obtain the instantaneous frequency expression of the logarithmic frequency sweep.

[0121] According to the derivative relationship between the instantaneous frequency and the phase of the swept frequency tone, an instantaneous phase expression is obtained, and a time domain expression of the swept frequency tone is calculated based on the instantaneous phase expression.

[0122] The instantaneous phase expression is: The time domain expression is

[0123] In the formula, is the instantaneous phase of the swept tone at time t, p(t) is the time domain of the swept tone at time t, and A(t) is the amplitude of the swept tone at time t.

[0124] In this embodiment, in the otoacoustic emission signal detection device provided by the present invention, the first algorithm processing module is specifically used for:

[0125] A first sweep frequency test paradigm of distortion product otoacoustic emission signals is obtained, wherein the first sweep frequency test paradigm includes a first single-segment test sound paradigm and a second single-segment test sound paradigm respectively corresponding to a pair of sweep frequency stimulation sounds with fixed instantaneous frequencies.

[0126] The online whole signal artifact removal algorithm is called to perform filtering and tracking on the first frequency sweep test paradigm to remove the stimulus sound echo artifact in the first frequency sweep test paradigm and obtain the first echo signal.

[0127] In this embodiment, in the otoacoustic emission signal detection device provided by the present invention, the first algorithm processing module is further configured to:

[0128] A second frequency sweep test paradigm is used to obtain the stimulation frequency otoacoustic emission signal. The second frequency sweep test paradigm includes at least four test sound paradigms based on two-tone suppression, namely, a first segmented test sound paradigm and a second segmented test sound paradigm corresponding to stimulation sounds with opposite phases, and a third segmented test sound paradigm and a fourth segmented test sound paradigm corresponding to stimulation sounds with opposite phases and suppression sounds with the same phase.

[0129] The online whole-segment signal artifact removal algorithm is called to perform dual-tone suppression superposition on the second frequency sweep test paradigm to remove the echo signal artifacts of the stimulation sound and the inhibition sound in the second frequency sweep test paradigm to obtain the second echo signal.

[0130] In this embodiment, the otoacoustic emission signal detection device provided by the present invention further includes an online whole-segment signal artifact removal module, which is used to:

[0131] When the recovered signals reach the preset stacking average, the reference signal corresponding to each recovered signal and the power of each reference signal are calculated to obtain the median and interquartile range of all powers.

[0132] When the power is not within the set range, the recovery signal corresponding to the power not within the set range is marked as an outlier, and the entire recovery signal corresponding to the outlier is deleted.

[0133] In this embodiment, in the otoacoustic emission signal detection device provided by the present invention, the second algorithm processing module is specifically used for:

[0134] The offline segmented signal artifact removal algorithm is called to remove the short-term noise in the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and the first and second echo signals without artifact interference are obtained by superposition and averaging of the segmented signals.

[0135] The second echo signal without artifact interference is filtered and tracked to filter out the otoacoustic emission signal of the stimulation frequency induced by the inhibitory sound in the second echo signal, so as to obtain the third echo signal.

[0136] In this embodiment, the otoacoustic emission signal detection device provided by the present invention further includes a tracking filter building module, which is used to:

[0137] The instantaneous filtering frequency of the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal is obtained, and the pole frequency, the zero frequency and the instantaneous gain are determined.

[0138] Based on the pole frequency, the zero frequency and the instantaneous gain, the expression of the filter is obtained to construct a tracking filter for filtering and tracking of otoacoustic emission signals.

[0139] The pole frequency is the instantaneous filtering frequency corresponding to the retained otoacoustic emission signal, and the zero frequency is the instantaneous filtering frequency corresponding to the otoacoustic emission signal to be filtered out.

[0140] Fig.12 An example of a physical structure diagram of an electronic device is shown. The electronic device may be a smart terminal, and its internal structure diagram may be as follows: Fig.12 As shown. The electronic device includes a processor, an internal memory and a network interface connected through a system bus. Among them, the processor of the electronic device is used to provide computing and control capabilities. The memory of the electronic device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the electronic device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, an otoacoustic emission signal detection method is implemented, and the method includes:

[0141] Obtaining otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweeping, wherein the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals;

[0142] The online whole-segment signal artifact removal algorithm is called to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal are subjected to filtering tracking and two-tone suppression superposition to obtain the first acquisition signal and the second acquisition signal respectively;

[0143] Calling the offline segmented signal artifact removal algorithm to remove artifacts from the first acquisition signal and the second acquisition signal, and performing segmented signal stacking and averaging operations on the first acquisition signal and the second acquisition signal processed by the offline segmented signal artifact removal algorithm, so as to extract the third acquisition signal by filtering and tracking the second acquisition signal after the stacking and averaging operations;

[0144] The first and third echo signals after the stacked average operation are processed by FFT operation to obtain the amplitude signal and phase signal corresponding to the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal with high frequency resolution, respectively;

[0145] Among them, the first collected signal is the collected signal of the distortion product otoacoustic emission signal, the second collected signal includes the collected signal of the stimulation frequency otoacoustic emission signal and the stimulation frequency otoacoustic emission signal induced by the inhibitory sound, and the third collected signal is the collected signal of the stimulation frequency otoacoustic emission signal.

[0146] Those skilled in the art will understand that Fig.12 The structure shown in the figure is only a block diagram of a part of the structure related to the scheme of the present invention, and does not constitute a limitation on the electronic device to which the scheme of the present invention is applied. The specific electronic device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0147] On the other hand, the present invention further provides a computer storage medium storing a computer program, which implements the above-mentioned otoacoustic emission signal detection method when executed by a processor.

[0148] In another aspect, a computer program product or a computer program is provided, the computer program product or the computer program comprising computer instructions, the computer instructions being stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium, and when the processor executes the computer instructions, the above-mentioned otoacoustic emission signal detection method is implemented.

[0149] Those of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be implemented by instructing related hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided by the present invention may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory may include random access memory (RAM) or external cache memory.

[0150] By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0151] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0152] The above-mentioned embodiments only express several implementation methods of the present invention, and the description thereof is relatively specific and detailed, but it cannot be understood as limiting the scope of the patent of the present invention. It should be pointed out that, for ordinary technicians in this field, several variations and improvements can be made without departing from the concept of the present invention, which all belong to the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention shall be subject to the attached claims.

Claims

1. A method for detecting otoacoustic emission signals, characterized in that: The method comprises: Acquiring otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweep, wherein the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals; Calling the online whole-segment signal artifact removal algorithm to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and performing filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after artifact removal, respectively, to obtain a first acquisition signal and a second acquisition signal; Calling an offline segmented signal artifact removal algorithm to remove artifacts from the first and second acquisition signals, and performing a segmented signal stacking average operation on the first and second acquisition signals processed by the offline segmented signal artifact removal algorithm, so as to extract a third acquisition signal by filtering and tracking the second acquisition signal after the stacking average operation; The first sampled signal and the third sampled signal after the stacked average operation are processed by FFT operation to obtain amplitude signals and phase signals corresponding to the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal with high frequency resolution, respectively; Among them, the first acquisition signal is the acquisition signal of the distortion product otoacoustic emission signal, the second acquisition signal includes the acquisition signal of the stimulation frequency otoacoustic emission signal and the stimulation frequency otoacoustic emission signal induced by the inhibitory sound, and the third acquisition signal is the acquisition signal of the stimulation frequency otoacoustic emission signal.

2. The otoacoustic emission signal detection method according to claim 1, characterized in that: The step of acquiring otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweeping comprises: Obtaining basic parameters of the frequency sweep, the basic parameters at least including a frequency sweep range, a frequency sweep mode, and a frequency sweep rate, and determining the relationship between the instantaneous frequency of the frequency sweep tone and the frequency sweep time according to the frequency sweep rate to obtain an instantaneous frequency expression of the logarithmic frequency sweep; According to the derivative relationship between the instantaneous frequency and the phase of the swept frequency tone, an instantaneous phase expression is obtained, and based on the instantaneous phase expression, a time domain expression of the swept frequency tone is calculated; Among them, the instantaneous phase expression is: , the time domain expression is ; In the formula, For sweep tone t The instantaneous phase of the moment, For sweep tone t The time domain of time, For sweep tone t The amplitude of the moment, Represents the inverse function of the instantaneous frequency at time t.

3. The otoacoustic emission signal detection method according to claim 2, characterized in that: The calling of the online whole signal artifact removal algorithm performs artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal respectively, and performs filtering tracking and dual-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after the artifact removal respectively, including: Acquire a first sweep frequency test paradigm of the distortion product otoacoustic emission signal, wherein the first sweep frequency test paradigm includes a first single-segment test sound paradigm and a second single-segment test sound paradigm respectively corresponding to a pair of sweep frequency stimulation sounds having a fixed ratio of instantaneous frequencies; The online whole-segment signal artifact removal algorithm is called to perform filtering and tracking on the first frequency sweep test paradigm to remove the stimulus sound echo artifact in the first frequency sweep test paradigm and obtain the first echo signal.

4. The method for detecting otoacoustic emission signals according to claim 3, characterized in that: The calling of the online whole signal artifact removal algorithm performs artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, respectively, and performs filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after the artifact removal, and further includes: Acquire a second frequency sweep test paradigm of the stimulation frequency otoacoustic emission signal, wherein the second frequency sweep test paradigm includes at least four test sound paradigms based on two-tone suppression, namely, a first segmented test sound paradigm and a second segmented test sound paradigm corresponding to stimulation sounds with opposite phases, and a third segmented test sound paradigm and a fourth segmented test sound paradigm corresponding to stimulation sounds with opposite phases and suppression sounds with the same phase; The online whole-segment signal artifact removal algorithm is called to perform dual-tone suppression superposition on the second frequency sweep test paradigm to remove the echo signal artifacts of the stimulation sound and the suppression sound in the second frequency sweep test paradigm to obtain the second echo signal.

5. The otoacoustic emission signal detection method according to claim 4, characterized in that: The method further comprises: When the recovered signal reaches the preset stacking average, the reference signal corresponding to each recovered signal and the power of each reference signal are calculated to obtain the median and interquartile range of all powers; When the power is not within the set range, the recovery signal corresponding to the power not within the set range is marked as an outlier, and the entire recovery signal corresponding to the outlier is deleted; Among them, the setting range is , where is the median, is the interquartile range.

6. The otoacoustic emission signal detection method according to claim 1, characterized in that: The method of calling the offline segmented signal artifact removal algorithm to remove artifacts from the first and second acquisition signals, and performing segmented signal stacking and averaging operations on the first and second acquisition signals processed by the offline segmented signal artifact removal algorithm, so as to extract the third acquisition signal by filtering and tracking the second acquisition signal after the stacking and averaging operations, includes: Calling the offline segmented signal artifact removal algorithm to remove short-term noise in the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and obtaining the first and second recollection signals without artifact interference by superposition and averaging operation of the segmented signals; The second echo signal without artifact interference is filtered and tracked to filter out the otoacoustic emission signal of the stimulation frequency induced by the suppression sound in the second echo signal, so as to obtain the third echo signal.

7. The otoacoustic emission signal detection method according to any one of claims 1 to 6, characterized in that: The method further comprises: Obtaining instantaneous filtering frequencies of the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and determining pole frequency, zero frequency and instantaneous gain; Based on the pole frequency, the zero frequency and the instantaneous gain, an expression of a filter is obtained to construct a tracking filter for filtering and tracking otoacoustic emission signals; The pole frequency is the instantaneous filtering frequency corresponding to the retained otoacoustic emission signal, and the zero frequency is the instantaneous filtering frequency corresponding to the otoacoustic emission signal to be filtered out.

8. An otoacoustic emission signal detection device, characterized in that: The device comprises: A signal acquisition module, used to acquire otoacoustic emission signals at multiple frequencies in a single frequency sweep test by exponential frequency sweep, wherein the otoacoustic emission signals include distortion product otoacoustic emission signals and stimulation frequency otoacoustic emission signals; The first algorithm processing module is used to call the online whole-segment signal artifact removal algorithm to perform artifact removal on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal, and perform filtering tracking and two-tone suppression superposition on the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal after the artifact removal, to obtain a first recollection signal and a second recollection signal respectively; The second algorithm processing module is used to call the offline segmented signal artifact removal algorithm to remove artifacts from the first and second acquisition signals, and perform segmented signal stacking and averaging operations on the first and second acquisition signals processed by the offline segmented signal artifact removal algorithm, so as to extract the third acquisition signal by filtering and tracking the second acquisition signal after the stacking and averaging operations; An FFT operation module is used to process the first acquisition signal and the third acquisition signal after the stacked average operation by FFT operation to obtain amplitude signals and phase signals corresponding to the distortion product otoacoustic emission signal and the stimulation frequency otoacoustic emission signal with high frequency resolution, respectively; Among them, the first acquisition signal is the acquisition signal of the distortion product otoacoustic emission signal, the second acquisition signal includes the acquisition signal of the stimulation frequency otoacoustic emission signal and the stimulation frequency otoacoustic emission signal induced by the inhibitory sound, and the third acquisition signal is the acquisition signal of the stimulation frequency otoacoustic emission signal.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images