Vibration source estimation device, vibration source estimation method, and vibration source estimation program

A vibration source estimation method, comprising:

• receiving data samples of a plurality of vibration signals generated at a vibration position, the data samples being associated with a plurality of frames;
• calculating a cross correlation function for each frame based on the data samples;
• calculate a background noise model that models background noise using the cross correlation function;
• calculate a noise-subtracted cross correlation function by subtracting the background noise, determined by the background noise model, from the cross correlation function;
• calculating a weight for each frame based on a signal-to-noise ratio calculated from the noise-subtracted cross correlation function;
• multiplying the noise-subtracted cross correlation function by the weight to be a weighted cross correlation function for each frame;
• calculating a sum of weighted cross correlation functions for the plurality of frames as a weighted-cross-correlation function;
• and
• estimating the vibration position based on the weighted-cross-correlation function.

A non-transitory computer-readable storage medium storing instructions that, when executed by a computer, cause the computer to perform a method, the method comprising:

• receiving data samples of a plurality of vibration signals generated at a vibration position, the data samples being associated with a plurality of frames;
• calculating a cross correlation function based on the data samples;
• calculating a background noise model that models background noise using the cross correlation function;
• calculating a noise-subtracted cross correlation function by subtracting the background noise, determined by the background noise model, from the cross correlation function;
• calculating a weight for each frame based on a signal-to-noise ratio calculated from the noise-subtracted cross correlation function;
• multiplying the noise-subtracted cross correlation function by the weight to be a weighted cross correlation function for each frame;
• calculating a sum of weighted cross correlation functions for the plurality of frames as a weighted-cross-correlation function;
• and
• estimating the vibration position based on the weighted-cross-correlation function.