This invention relates to techniques for detecting anomalous sound of a machine from sound signals. For instance, it relates to techniques for generating an acoustic feature extractor for anomalous sound detection.
In a factory or the like, even a shutdown of an industrial device installed therein, such as a large-sized manufacturing or molding machine, due to a failure significantly hinders the factory's operation. It is thus necessary to routinely monitor the operational status of the device and immediately take measures in the event of an anomaly. One solution is to regularly dispatch a maintenance person to the field from an industrial device management service to check machine components for wear and the like. Since this requires enormous personnel and/or traveling costs and labor, however, it is difficult to implement this on all industrial devices and/or factories.
A solution to this is to install a microphone inside a machine so as to routinely monitor the operational sound of the machine. Through analysis of the operational sound, any occurrence of sound that is likely to be an anomaly (i.e., anomalous sound) is detected and an alert is raised, thereby solving the anomaly. However, setting the types of anomalous sound and/or methods for detecting them for every machine type or individual unit is even more expensive than manual monitoring. Thus, there is a need for automated design of rules for automatically detecting anomalous sound.
As a way to address this problem, anomalous sound detection based on statistical approaches is well known (see Non-patent Literature 1, for instance). Anomalous sound detection based on statistical approaches is generally classified into supervised anomalous sound detection and unsupervised anomalous sound detection. In supervised anomalous sound detection, a classifier is trained from training data for normal sound and anomalous sound; whereas in unsupervised anomalous sound detection, a classifier is trained only from training data for normal sound. For an industrial application, unsupervised anomalous sound detection is often used because training data for anomalous sound is difficult to collect.
A training/detection flow for unsupervised anomalous sound detection is as shown in
A challenge in unsupervised anomaly detection is design of an acoustic feature extractor F(·). In supervised anomalous sound detection, an acoustic feature that can correctly identify a target of determination is manually designed. For example, if it is known that the normal sound is a sinusoidal wave of 1000 Hz and the anomalous sound is a sinusoidal wave of 2000 Hz, a log power of a Mel filter bank (log-MFBO) is extracted per frame because the two sounds have different tone timbres. If the normal sound is steady engine noise and the anomalous sound is a “bumping” sound emitted by hitting of devices, a temporal difference of the power of a Mel filter bank (ΔMFBO) is extracted because the anomalous sound is a sporadic sound. “Deep training”, a type of supervised training, is said to be able to automatically design an acoustic feature from training data.
In unsupervised anomaly detection, however, it cannot be known anomalous sound with what kind of sound characteristics will occur. It is accordingly difficult to manually design the acoustic feature extractor and the use of deep training is also difficult. For example, if the anomalous sound is assumed to be a sinusoidal wave of 2000 Hz and log-MFBO is used as the acoustic feature due to the fact that the normal sound was a sinusoidal wave of 1000 Hz, anomalous sound like a “bumping” sound emitted by hitting of devices cannot be detected. The opposite is also true. This has forced use of, for example, Mel filter bank cepstrum coefficient (MFCC), which is a general-purpose acoustic feature, resulting in low accuracy of detection compared to supervised training.
An object of the present invention is to provide an anomalous sound detection training apparatus that is capable of generating an acoustic feature extractor for anomalous sound detection irrespective of whether training data for anomalous sound is available or not, as well as an acoustic feature extraction apparatus, an anomalous sound sampling apparatus, methods, and programs associated with the anomalous sound detection training apparatus.
An anomalous sound detection training apparatus according to an aspect of this invention includes: a first acoustic feature extraction unit that extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating unit that updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction unit that extracts an acoustic feature of anomalous sound based on simulated anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating unit that updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit. Processing by the first acoustic feature extraction unit, the normal sound model updating unit, the second acoustic feature extraction unit, and the acoustic feature extractor updating unit is repeatedly performed.
An anomalous sound detection training apparatus according to an aspect of this invention includes: a first acoustic feature extraction unit that extracts an acoustic feature of normal sound based on training data for normal sound by using an acoustic feature extractor; a normal sound model updating unit that updates a normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit; a second acoustic feature extraction unit that extracts an acoustic feature of anomalous sound based on training data for anomalous sound and extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor; and an acoustic feature extractor updating unit that updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit. Processing by the first acoustic feature extraction unit, the normal sound model updating unit, the second acoustic feature extraction unit, and the acoustic feature extractor updating unit is repeatedly performed.
It is possible to generate an acoustic feature extractor for anomalous sound detection irrespective of whether training data for anomalous sound is available or not.
[Overview of the Invention]
[[Structure of the Invention]]
A hierarchical structure of the present invention is shown in
(I) Irrespective of supervised/unsupervised, a “normal sound model” and an “acoustic feature extractor” are optimized simultaneously from training data.
(II) In the case of unsupervised, training data for anomalous sound is sampled.
(III) Irrespective of supervised/unsupervised, the “acoustic feature extractor” is represented as a matrix product and a logarithm operation, or as a multilayer perceptron, for example.
(IV) A training/detection algorithm that simultaneously implements (I) to (III) is provided.
The following will describe (I) to (III). Note that for the sake of description, they will be described in the order of (III), (I), (H). It will be followed by description on [Anomalous sound detection training apparatus and method] for performing them.
[[(III) Definition of Acoustic Feature Extractor]]
A series derived by Fourier-transformation of an observation signal and extraction of a power spectrum or amplitude spectrum is expressed as Xω,t∈R+Ω×T. Here, ω={1, 2, . . . , Ω} and t={1, 2, . . . , T} represent indices of frequency and time, respectively, Ω represents the number of frequency bins, and T represents the number of observed frames. R+ represents a set of non-negative real numbers.
Most cases of acoustic feature extraction can be represented as a matrix operation on this series with a logarithmic function. For example, log-MFBO can be described as spectrum Xt=(X1,t, . . . , XΩ,t)T of all frequencies in a certain frame multiplied by M Mel filter banks mm, of which the log is then taken.
Here, T represents transposition. Also, when ΔMFBO is defined as the temporal difference of Mel filter bank's output, it can be described as follows using a difference matrix D.
Here, EM is a unit matrix of M dimensions.
Herein, an acoustic feature extractor is defined with extension of this concept. First, past P1 frame through future P2 frame are taken from determination target frame t, and a vertical connection of them is defined as xt, where xt is represented as shown in
Here, Q=Ω×(P1+P2+1) holds. The modification to the first to second lines of formula (1′) is rewriting for the purpose of simplifying the notation; the contents and order of the elements are unchanged. An operation that multiplies xt by a matrix A∈RD×Q and takes the log thereof is defined as the acoustic feature extractor. That is, it will be as follows.
Thus, optimization of the acoustic feature extractor can be said to be optimization of matrix A.
[[(I) Simultaneous Optimization of Normal Sound Model and Acoustic Feature Extractor]]
First, a normal sound model (a probability density function followed by a feature amount extracted from normal sound) is defined as p(F(A,x)|z=0). Then, a negative logarithmic likelihood (the degree of anomaly) can be computed as below.
g(A,x)=−ln p(F(A,x)|z=0) (5)
Here, although the normal sound model may be any of various types, when it is a Gaussian mixture distribution, for example, it can be described as:
Here, K is the number of mixtures, and k, Σk, and wk are a mean vector, a covariance matrix, and a weight parameter for the kth distribution, respectively.
In anomalous sound detection, when the value of g(A,x) is greater than threshold ϕ, it is determined to be “anomalous”; when the value is smaller than threshold ϕ, it is determined to be “normal”. Although this determination can be described as a step function, for the simplicity of the subsequent calculations, it is described with a sigmoid function:
The sigmoid function 0≤σ{α}≤1 outputs a value closer to 0 as α is smaller and outputs a value closer to 1 as α is greater. That is, when σ{g(A,x−ϕ)} is equal to or greater than 0.5, it is determined to be anomalous; when it is equal to or smaller than 0.5, it is determined to be normal.
An aim of an anomalous sound detection algorithm is to output a great a {g(A,x−ϕ)} when the observation represents anomalous sound and to output a small σ{g(A,x−ϕ)} when it represents normal sound. That is, minimizing the following formula is the aim of optimization of anomalous sound detection.
L=−∫σ{g(A,x)−ϕ}(p(F(A,x)|z=1)−p(F(A,x)|z=0))dx (9)
Here, p(F(A,x)|z=1) is an anomalous sound model (a probability density function followed by a feature amount extracted from anomalous sound).
From the foregoing, simultaneous optimization of the normal sound model and the acoustic feature extractor should be performed so as to simultaneously satisfy both:
Optimization of the acoustic feature extraction can be said to be design of the matrix A such that formula (11) is optimized.
Optimization of formula (10) can be carried out with a maximum likelihood method, for example. Thus, consider the optimization of formula (11). The foremost challenge in calculating formula (9) is the integration operation. Thus, this integration is replaced with arithmetic mean of training data. First, an objective function is defined as:
L=−∫σ{g(A,x)−ϕ}p(F(A,x)|z=1)dx+∫σ{g(A,x)−ϕ}p(F(A,x)|z=0)dx (12)
Each of the terms can then be regarded as an expectation operation for σ{G(A,x−ϕ)}, thus they can be described as:
L=−E[σ{g(A,x)−ϕ}]x|z=1+E[σ{g(A,x)−ϕ}]x|z=0 (13)
Here, if E[·]x|z=0 is an expectation operation relating to normal sound, while E[·]x|z=1 is an expectation operation relating to anomalous sound.
Here, if replacing the expectation operations with arithmetic means, formula (11) can be approximated as follows:
Here, N is a time index of training data for anomalous sound. That is, the accuracy of anomalous sound detection will be maximized by calculating the arithmetic means of Ψ1, . . . , T using training data for normal sound, calculating the arithmetic means of Ψ1, . . . , T using training data for anomalous sound, and designing A so that the sum of the two arithmetic means is minimized. In the case of supervised anomalous sound detection where training data for anomalous sound has been obtained, this can be performed without modification. In the case of unsupervised anomalous sound detection, N pieces of anomalous sound data are sampled (simulated) in accordance with the sampling algorithm described below. In that case, N should be set about N=T.
Now that the objective function has been set with formula (15), it is then minimized. For the minimization, a gradient method may be used. When a steepest descent method is used as the gradient method, the matrix may be optimized for each dimension d as:
Here, λ is step size. For this optimization, AdaGrad, that is, an improved algorithm of the steepest descent method, or the like may be used, of course.
In a case where the normal sound model is represented as a Gaussian mixture distribution, a gradient vector ∂J/∂ad can be calculated follows using the chain rule of partial differentiation.
Here, each of the components can be calculated by:
[[(II) Sampling of Training Data for Anomalous Sound]]
For calculation of formula (15) in unsupervised anomalous sound detection, how training data for anomalous sound is sampled will be discussed.
The distribution used with anomalous sound model p(F(A,x)|z=1) is arbitrary. In unsupervised training, no previous knowledge on anomalous sound has been obtained. In order to detect every kind of anomalous sound, it is desirable to construct an anomalous sound model with as least assumptions as possible. That is, it is desirable to sample anomalous sound only from generally known findings about the physical characteristics or frequency of appearance of an amplitude spectrum and/or power spectrum.
In the following, an anomalous sound model is constructed with assumptions made only for the occurrence frequency, sparsity, and continuity in temporal frequency direction of amplitude values.
First, consider the magnitude of Xω,n.
A distribution that satisfies the characteristics is an exponential distribution, for example. An exponential distribution is a distribution such that:
In the distribution, the frequency of appearance of o decreases as the value of o increases. β is an average value. In the case of sampling Xω,n from an exponential distribution, it will be as follows.
Since β is a parameter relating to the volume of observed sound, it can be set from the training data for normal sound. For example, it can be set as:
If previous knowledge on stationary noise contained in observation signals or the like is available, this sampling may be modified as desired. For example, if an average spectrum Yω of stationary noise is known, sampling may be performed such that:
β here can be calculated by:
When sampled with formula (29), a volume average βω of anomalous sound will be constant independently of time. Depending on the distance to a microphone or the type of anomalous sound, the amplitude value, or the volume, might increase or decrease. Hence, it is desirable to vary the value of βω,n on a per-n basis. For representing this, formula (29) may be modified as follows, for example:
Here, θ is a parameter in unit of dB for correcting the maximum volume of anomalous sound, and Uniform(a,b) is a continuous uniform distribution of maximum a and minimum b. For example, θ may be set to about 0 to 12. As another way, similar effects could be provided by sampling θ from a normal distribution or the like and setting it always to, for example, 1, without performing sampling of bω.
It is also possible to focus on the sparsity of an amplitude spectrum (almost all of the elements being 0).
Here, Cat(π) is category distribution and 1Ω is a vector with a length of Ω whose elements are all 1. This makes only one spectrum bin become active at each time via zn, thereby driving the anomalous sound spectrum E to be sparse.
Also, sampling with formula (35) results in generation of a value of Eω,n that is completely independent from the temporal frequency direction. A sound signal typically changes smoothly in terms of time, so that the spectrum thereof smoothly changes in the temporal frequency direction. For achieving this, Fω,n may be smoothed in the temporal frequency direction.
While the smoothing can be done in various ways, it can be performed using a Gaussian filter as follows, for example:
Xω,n=Yω+conv2(Eω,n,G(F,T)) (38)
Here, conv2 is a convolution operation of two dimensions, and G(F,T) is a two-dimensional Gaussian filter with a size of F×T. Although the parameters should be controlled experimentally, they may be set at about F=11 and T=31, for example.
[Anomalous Sound Detection Training Apparatus and Method]
As exemplarily shown in
<Spectrum Calculation Unit 1>
Training data for normal sound is input to the spectrum calculation unit 1.
The training data for normal sound is sound signals collected by a microphone installed near (or alternatively inside) a device for which anomalous sound is to be detected. While the sampling frequency for collecting sound signals should be set so that the frequencies of normal sound can be recorded, it may be set at 16 kHz, for example.
The spectrum calculation unit 1 converts the training data into a frequency domain via an approach such as short-time Fourier transformation (STFT).
It then takes an absolute value of a complex spectrum to determine an amplitude spectrum Xω,t (step S1). Alternatively, Xω,t may be a power spectrum derived by determining the square root of the amplitude spectrum. Here, Fourier transformation length may be set at 32 ms and shift width may be set at 16 ms, for example.
In the following, Xω,t as either an amplitude spectrum or a power spectrum will be abbreviated as spectrum Xω,t. The spectrum Xω,t is output to the first acoustic feature extraction unit 3 and the second acoustic feature extraction unit 6.
If training data for anomalous sound is available, namely in the case of supervised training, the spectrum calculation unit 1 calculates the spectrum of anomalous sound based on the training data for anomalous sound. The spectrum of anomalous sound is the amplitude spectrum or the power spectrum of the anomalous sound, for example. The calculated spectrum of anomalous sound is output to the first acoustic feature extraction unit 3 and the second acoustic feature extraction unit 6.
<Initialization Unit 2>
Preset parameters for a training algorithm and preset parameters for an acoustic feature extractor are input to the initialization unit 2.
The various parameters may be set such that P1=3, P2=3, ε=min(Xω,t)/2, ϕ=100, and λ=5.0×10−3, for example.
The initialization unit 2 initializes the matrix A for the acoustic feature extractor and the normal sound model (step S2). After being initialized, the matrix A and the normal sound model are output to the first acoustic feature extraction unit 3.
The matrix A may be initialized with uniform random numbers, for example. While initialization of the normal sound model varies depending on the type of the probability density function, when a Gaussian mixture distribution is used, for example, the normal sound model can be initialized via a k-means algorithm or the like.
<First Acoustic Feature Extraction Unit 3>
To the first acoustic feature extraction unit 3, the acoustic feature function initialized by the initialization unit 2 and the training data for normal sound are input.
The first acoustic feature extraction unit 3 extracts an acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor (step S3). In this example, as the training data for normal sound, the first acoustic feature extraction unit 3 uses the spectrum Xω,t calculated by the spectrum calculation unit 1 based on the training data for normal sound to extract the acoustic feature of normal sound. The extracted acoustic feature is output to the normal sound model updating unit 4.
A first round of processing by the first acoustic feature extraction unit 3 is performed by using the acoustic feature extractor initialized by the initialization unit 2. Second and subsequent rounds of processing by the first acoustic feature extraction unit 3 are performed by using the acoustic feature extractor as updated by the acoustic feature extractor updating unit 7.
If the spectrum of anomalous sound has been calculated by the spectrum calculation unit 1 based on training data for anomalous sound, the first acoustic feature extraction unit 3 extracts the acoustic feature of anomalous sound based on this calculated spectrum of anomalous sound by using the acoustic feature extractor.
<Normal Sound Model Updating Unit 4>
The normal sound model updating unit 4 updates the normal sound model by using the acoustic feature extracted by the first acoustic feature extraction unit 3 (step S4).
While the way of update varies depending on the type of the probability density function for the normal sound model, it may be updated by maximum likelihood estimation, for example. In a case where a Gaussian mixture distribution is used as the normal sound model, one to several steps of update may be performed in an EM algorithm (see Reference Literature 1, for instance).
<Anomalous Sound Sampling Unit 5>
The anomalous sound sampling unit 5 samples simulated anomalous sound (step S5). The sampled anomalous sound is output to the second acoustic feature extraction unit 6.
The anomalous sound sampling unit 5 performs sampling of simulated anomalous sound based on formula (26), for example. The parameter β in formula (26) can be determined based on formula (27), for example. In this case, for Xω,t in formula (27), the spectrum Xω,t of normal sound calculated by the spectrum calculation unit 1 may be used.
If training data for anomalous sound is available, namely in the case of supervised training, sampling is not performed. That is, the subsequent processing may be performed using the training data for anomalous sound as the sampling result.
<Second Acoustic Feature Extraction Unit 6>
The second acoustic feature extraction unit 6 extracts the acoustic feature of anomalous sound based on simulated anomalous sound and extracts the acoustic feature of normal sound based on the training data for normal sound by using the acoustic feature extractor (step S6). The acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted are output to the acoustic feature extractor updating unit 7.
In this example, as the training data for normal sound, the second acoustic feature extraction unit 6 uses the spectrum Xω,t calculated by the spectrum calculation unit 1 based on the training data for normal sound to extract the acoustic feature of normal sound.
The first round of processing by the second acoustic feature extraction unit 6 is performed by using the acoustic feature extractor initialized by the initialization unit 2. The second and subsequent rounds of processing by the second acoustic feature extraction unit 6 are performed by using the acoustic feature extractor as updated by the acoustic feature extractor updating unit 7.
<Acoustic Feature Extractor Updating Unit 7>
The acoustic feature extractor updating unit 7 updates the acoustic feature extractor by using the acoustic feature of anomalous sound and the acoustic feature of normal sound that have been extracted by the second acoustic feature extraction unit 6 (step S7). The updated acoustic feature extractor is output to the first acoustic feature extraction unit 3 and the second acoustic feature extraction unit 6. Also, the acoustic feature extractor as finally updated after repeated control by the control unit 8 is output as a final result of training performed by the anomalous sound detection training apparatus and method.
The acoustic feature extractor updating unit 7 updates the acoustic feature extractor based on formula (16), for example. Prior to execution of this update, threshold ϕ may be updated according to the formula below using training data for normal sound x1, . . . , xT and sampled anomalous sound data x1, . . . xN. This stabilizes convergence.
<Control Unit 8>
The control unit 8 performs control so as to repeatedly perform processing by the first acoustic feature extraction unit 3, the normal sound model updating unit 4, the second acoustic feature extraction unit 6, and the acoustic feature extractor updating unit 7 (steps S81 and S82). This repeated processing is performed until the acoustic feature extractor and the normal sound model converge.
For example, the control unit 8 performs control so as to repeatedly perform the processing by the second acoustic feature extraction unit 6 and the acoustic feature extractor updating unit 7 (hereinafter referred to as a first repeated processing). Setting a first convergence determination condition as that the number of executions of processing has reached a certain number of times (for example, 50) or that an update norm for the matrix A has fallen below a certain value, the control unit 8 performs control so as to repeatedly perform the first repeated processing until the first convergence condition is satisfied, or in other words, until it is determined that convergence has been reached. The first convergence determination condition may also be a different condition.
The control unit 8 subsequently performs control so as to repeatedly perform the processing by the first acoustic feature extraction unit 3 and the normal sound model updating unit 4 and the first repeated processing (hereinafter referred to as a second repeated processing). Setting a second convergence condition as that the number of executions of processing has reached a certain number of times (for example, 200), the control unit 8 performs control so as to repeatedly perform the second repeated processing until the second convergence condition is satisfied, or in other words, until it is determined that convergence has been reached. The second convergence determination condition may also be a different condition.
[Anomalous Sound Detection Apparatus and Method]
As exemplarily shown in
<Spectrum Calculation Unit 11>
Operational sound of a machine for which anomalous sound is to be detected is collected through a microphone. The sampling rate used for this collection is similar to the one used for training. The collected sound signal is input to the spectrum calculation unit 11.
The spectrum calculation unit 11 obtains spectrum. Xω,t based on the collected sound signal in a similar manner to the spectrum calculation unit 1 (step S11). The obtained spectrum Xω,t is output to the acoustic feature extraction unit 12.
<Acoustic Feature Extraction Unit 12>
The acoustic feature extraction unit 12 uses the Obtained spectrum Xω,t to extract the acoustic feature of the collected sound signal based on the acoustic feature extractor that has been output by the anomalous sound detection training apparatus and method as the final result of training (step S12). The extracted acoustic feature is output to the degree-of-anomaly calculation unit 13.
<Degree-of-Anomaly Calculation Unit 3>
The degree-of-anomaly calculation unit 13 calculates a negative logarithmic likelihood g(A,x) by using the extracted acoustic feature and a normal sound model that has been output by the anomalous sound detection training apparatus and method as the final result of training (step S13). Hereinafter, a negative logarithmic likelihood will be referred to as the degree of anomaly. The calculated degree of anomaly is output to the determination unit 14.
<Determination Unit 14>
The determination unit 14 outputs “anomalous” if the degree of anomaly of a current frame is equal to or greater than threshold ϕ (step S14). Although the threshold should be adjusted in accordance with the machine and/or environment in question, it may be set at about 1500, for example.
As with voice activity detection, “hangover”, which suppresses detection errors with heuristic rules, may be used. While the hangover processing to be applied may be any of various types, the hangover processing should be set in accordance with the types of false detection of anomalous sound.
As an example, musical noise that occurs during noise suppression could be determined as a sporadic anomalous sound. Since sound like a sporadic hitting sound often exhibits change in spectrum shape for 100 ms or more, the degree of anomaly would remain equal to or greater than the threshold for ┌(100/the frame shift width for STFT)┐ frames continuously. However, because in musical noise an anomalous amplitude spectrum value only occurs in the relevant frames, there would be several frames at most where the degree of anomaly remains equal to or greater than the threshold continuously. Thus, a rule for anomaly determination may be set as “outputting “anomalous” if the degree of anomaly remains equal to or greater than the threshold for F1 frames or more continuously”, for example.
As another example, it is also conceivable that anomalous sound continues for a long period of time with the degree of anomaly slightly below the threshold due to a low volume of the anomalous sound. In such a case, a rule like “detecting a sound as anomalous if a total sum of its degree of anomaly over the last F2 frames is equal to or greater than ϕ1” may be added as a determination rule for continuous anomalous sound. Although ϕ1 should be determined by tuning, it may be set at about ϕ1=F2×(ϕ−250), for example.
By detecting anomalous sound from a large-sized manufacturing or molding machine installed in a factory or the like using the anomalous sound detection apparatus and method as described above, faster handling of a failure and/or failure prediction become possible. This can contribute to increased efficiency of an industry, for example, manufacturing industry in particular.
[Program and Recording Medium]
In the case of implementing the processing in the anomalous sound detection training apparatus or the anomalous sound detection apparatus with a computer, the processing details of the functions to be provided by the anomalous sound detection training apparatus or the anomalous sound detection apparatus are described by a program. By the computer then executing the program, the processing is implemented on the computer.
The program describing the processing details may be recorded in a computer-readable recording medium. The computer-readable recording medium may be any kind of medium, for example, a magnetic recording device, an optical disk, a magneto-optical recording medium, and a semiconductor memory.
In addition, individual processing means may be embodied by execution of a predetermined program on a computer or at least some of their processing details may be implemented in hardware.
[Modifications]
Using the optimization metrics for acoustic feature described above, a multilayer perceptron can also be trained unsupervised. Specifically, although the acoustic feature extractor was described using formula (3), it may instead be implemented with a multilayer perceptron. For example, when a simplest fully-connected multilayer perceptron is used, the acoustic feature extractor would be as shown below. In this manner, the acoustic feature extractor may also be represented by a multilayer perceptron.
F(xτ)=W(L)gτ(L-1)+b(L) (40)
gτ(k)=σθ{W(k)gτ(k-1)+b(k)} (41)
Here, L, (k), and b(k) represent the number of layers of a neural network, a weight matrix, and a bias vector, respectively. L is a predetermined positive integer equal to or greater than 2. σθ is an activation function and may be implemented as a sigmoid function, for example. Input to the input layer may be such that gτ(1)=xτ or gτ(1)=ln xτ, for example. gτ(k) is input to the kth layer of the multilayer perceptron. In this case, optimization of the acoustic feature extractor is optimization of Θ={(W(k),b(k)}, (k∈{2, . . . , L}).
For this optimization, a gradient method may be used as in formula (16). If a steepest descent method is used as the gradient method, it may be done such that:
For this optimization, AdaGrad, that is, an improved algorithm of the steepest descent method, or the like may be used, of course. The objective function J is given by
Thus, formulae (42) and (43) can be described as:
The gradients:
can be calculated as follows using the chain rule of partial differentiation from formula (17).
In a multilayer perceptron, gradients in the form of formulae (47) and (48) can be efficiently calculated by using error back propagation (see Reference Literature 2, for instance).
In addition to being executed chronologically in the order described, the processing described for the anomalous sound detection training apparatus or the anomalous sound detection apparatus may also be executed in parallel or individually depending on the processing ability of the apparatus executing the processing or on any necessity.
It goes without saying that other modifications may be made as desired without departing from the scope of this invention.
Number | Date | Country | Kind |
---|---|---|---|
JP2016-074401 | Apr 2016 | JP | national |
JP2017-004933 | Jan 2017 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2017/013738 | 3/31/2017 | WO |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2017/171051 | 10/5/2017 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
20060171625 | Jones | Aug 2006 | A1 |
20140046878 | Lecomte | Feb 2014 | A1 |
20140311245 | Horoshenkov | Oct 2014 | A1 |
20140328487 | Hiroe | Nov 2014 | A1 |
20150241584 | Aarre | Aug 2015 | A1 |
Number | Date | Country |
---|---|---|
8-271573 | Oct 1996 | JP |
2012-181280 | Sep 2012 | JP |
5767825 | Aug 2015 | JP |
Entry |
---|
Combined Chinese Office Action and Search Report dated Mar. 25, 2020 in Chinese Patent Application No. 201780020698 5 (with English translation), citing document AO therein, 24 pages. |
Extended European Search Report dated Mar. 10, 2020 in European Patent Application No. 17775582.4, 10 pages. |
Office Action dated Aug. 20, 2019 in corresponding Japanese Patent Application No. 2018-509685 (with English Translation), 8 pages. |
Partial Supplementary European Search Report dated Nov. 22, 2019, in Patent Application No. 17775582.4, citing documents AX-AY therein, 13 pages. |
Bonfigli, R. et al., “A Real-Time Implementation of an Acoustic Novelty Detector on the BeagleBoard-XM”, 2014 6th European Embedded Design In Education and Research Conference, XP032661946, Sep. 11, 2014, pp. 307-311. |
Ntalampiras, S. et al., “Probabilistic Novelty Detection for Acoustic Surveillance Under Real-World Conditions”, IEEE Transactions on Multimedia, XP011330160, vol. 13, No. 4, Aug. 1, 2011, pp. 713-719. |
Tara N. Sainath, et al., “Learning filter banks within a deep neural network framework”, ASRU 2013 IEEE Workshop, Jan. 2014, pp. 297-302. |
Luke Pfister, et al., “Learning sparsifying filter banks”, SPIE Proceedings, Wavelets and Sparsity XVI, vol. 9597, Aug. 2015, 11 pages. |
Alain Biem, et al., “An application of discriminative feature extraction to filter-bank-based speech recognition”, IEEE transactions on speech an audio processing, vol. 9, No. 2, Feb. 2001, pp. 96-110. |
Ide, et al., “Anomaly Detection and Change Detection”, Kodansha, Aug. 2015, pp. 6-7 and cover pages (with partial English language translation). |
Sadanori Konishi, “Introduction to Multivariate Analysis, Appendix C: EM algorithm”, Iwanami Shoten, Jan. 2010, pp. 294-298 and cover pages (with partial English language translation). |
Hideki Asoh, et al., “Deep Learning”, The Japanese Society for Artificial Intelligence, Kindai kagaku sha, Oct. 2015, p. 145 and cover pages (with partial English language translation). |
Number | Date | Country | |
---|---|---|---|
20190120719 A1 | Apr 2019 | US |