KR101847638B1 - 음향 분석 방법 및 이를 이용한 비파괴 검사 장치 - Google Patents

음향 분석 방법 및 이를 이용한 비파괴 검사 장치 Download PDF

Info

Publication number
KR101847638B1
KR101847638B1 KR1020160156741A KR20160156741A KR101847638B1 KR 101847638 B1 KR101847638 B1 KR 101847638B1 KR 1020160156741 A KR1020160156741 A KR 1020160156741A KR 20160156741 A KR20160156741 A KR 20160156741A KR 101847638 B1 KR101847638 B1 KR 101847638B1
Authority
KR
South Korea
Prior art keywords
short
value
frame
sound
term
Prior art date
Application number
KR1020160156741A
Other languages
English (en)
Inventor
김영섭
김정호
Original Assignee
주식회사 시그널웍스
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 주식회사 시그널웍스 filed Critical 주식회사 시그널웍스
Priority to KR1020160156741A priority Critical patent/KR101847638B1/ko
Priority to PCT/KR2017/013422 priority patent/WO2018097620A1/ko
Application granted granted Critical
Publication of KR101847638B1 publication Critical patent/KR101847638B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/045—Analysing solids by imparting shocks to the workpiece and detecting the vibrations or the acoustic waves caused by the shocks
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01Q—SCANNING-PROBE TECHNIQUES OR APPARATUS; APPLICATIONS OF SCANNING-PROBE TECHNIQUES, e.g. SCANNING PROBE MICROSCOPY [SPM]
    • G01Q10/00—Scanning or positioning arrangements, i.e. arrangements for actively controlling the movement or position of the probe
    • G01Q10/04—Fine scanning or positioning
    • G01Q10/06—Circuits or algorithms therefor
    • G06F19/00—

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

음향 분석 방법 및 이를 이용한 비파괴 검사 장치가 개시된다. 본 발명의 일 실시예에 따른 비파괴 검사 장치는, 피 검사물에 대한 타격에 의해 발생하는 음향 신호를 획득하는 음향 신호 획득부, 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하는 전처리부, 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하는 단구간 프레임 에너지 산출부, 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하는 스펙트럼 밀도 산출부, 각 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하는 판별 추정값 산출부, 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정하는 임계값 설정부, 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 임계값을 기반으로 음향 신호에서 타격음 구간을 결정하는 타격음 구간 결정부, 및 기 결정된 타격음 구간에 대응하는 음향 신호를 기반으로 피 검사물의 결함 여부를 분석하는 결함 여부 분석부를 포함한다.

Description

음향 분석 방법 및 이를 이용한 비파괴 검사 장치{METHOD FOR ANALYZING ACOUSTIC AND APPARATUS FOR TESTING NON DESTRUCTIVE OF USING THE METHOD}
본 발명의 실시예는 음향 분석 기술에 관한 것으로, 보다 상세하게는 음향 신호를 이용한 비과피 검사 기술에 관한 것이다.
기계 부품 또는 제품 등의 결함 여부를 검사하기 위한 방법으로 다양한 비파괴 검사법(Non Destructive Testing : NDT)이 적용되고 있다. 그 중 음향 공진(Acoustic Resonance)을 이용한 비파괴 검사법은 신속하면서 간단하게 피검사물의 결함을 테스트 할 수 있어 많이 사용되고 있다.
음향 공진을 이용한 비파괴 검사법에는 피검사물을 타격하여 발생한 음향 신호를 분석하여 피검사물의 결함 여부를 판단하는 방식이 있다. 종래에는 인라인 벨트에 장착된 센서의 트리거 신호 및 제품 별 발생하는 타격음의 지속 시간을 고려한 설정값을 통해 타격음을 획득하는 방식을 사용하였다.
종래의 방식에 의하면, 다품종 소량 생산의 경우, 피 검사 제품이 자주 바뀌기 때문에 각 피 검사 제품에 따른 설정값을 계속 변경해주어야 하는 번거로움이 있게 된다. 그리고, 노이즈 환경에서는 타격음에 대해 오분류할 가능성이 높아 신뢰성이 저하되는 문제점이 있게 된다.
한국등록특허공보 제10-1180687호(2012.09.07)
본 발명의 실시예는 잡음 환경에서 유효음 구간을 정확히 추출할 수 있는 음향 분석 방법 및 이를 이용한 비파괴 검사 장치를 제공하기 위한 것이다.
본 발명의 일 실시예에 따른 비파괴 검사 장치는, 피 검사물에 대한 타격에 의해 발생하는 음향 신호를 획득하는 음향 신호 획득부; 상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하는 전처리부; 상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하는 단구간 프레임 에너지 산출부; 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하는 스펙트럼 밀도 산출부; 상기 각 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하는 판별 추정값 산출부; 상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정하는 임계값 설정부; 상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는 타격음 구간 결정부; 및 상기 결정된 타격음 구간에 대응하는 음향 신호를 기반으로 상기 피 검사물의 결함 여부를 분석하는 결함 여부 분석부를 포함한다.
상기 스펙트럼 밀도 산출부는, 하기의 수학식을 통해 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density : SD)를 산출할 수 있다.
Figure 112016114802737-pat00001
M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
Ni : 기 설정된 잡음 신호의 i번째 MFCC 계수의 값
상기 판별 추정값 산출부는, 하기의 수학식을 통해 상기 판별 추정값을 산출할 수 있다.
Figure 112016114802737-pat00002
Y : 해당 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값
M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
Wi : i번째 MFCC 계수에 대해 기 설정된 가중치
V : 바이어스(Bias)로서 기 설정된 상수값
상기 임계값 설정부는, 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출하고, 상기 산출한 이동 평균값을 기반으로 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 임계값을 각각 설정할 수 있다.
상기 타격음 구간 결정부는, 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값과 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값 및 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 잡음 성분이 포함되는 정도에 따라 각각 설정되는 파라미터 가중 팩터를 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정할 수 있다.
상기 임계값 설정부는, 하기의 수학식을 통해 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출할 수 있다.
Figure 112016114802737-pat00003
p1(n)(즉, j=1인 경우) : n번째 단구간 프레임의 에너지
p2(n)(즉, j=2인 경우) : n번째 단구간 프레임의 스펙트럼 밀도
p3(n)(즉, j=3인 경우) : n번째 단구간 프레임의 판별 추정값
Figure 112016114802737-pat00004
: n번째 단구간 프레임의 에너지에 대한 이동 평균값
Figure 112016114802737-pat00005
: n번째 단구간 프레임의 스펙트럼 밀도에 대한 이동 평균값
Figure 112016114802737-pat00006
: n번째 단구간 프레임의 판별 추정값에 대한 이동 평균값
α : 이동 평균값의 변동량을 결정하는 상수
상기 타격음 구간 결정부는, 하기의 수학식을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정할 수 있다.
Figure 112016114802737-pat00007
R(n) : n번째 단구간 프레임이 타격음을 포함하는지 여부를 판단하기 위한 함수
δj(n) : n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값과 n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값
N : 타격음 구간을 결정하기 위해 분석 대상으로 취하는 n번째 단구간 프레임과 인접한 단구간 프레임의 개수
m : n번째 단구간 프레임과 인접한 단구간 프레임의 번호
βj : 해당 단구간 프레임에 잡음 성분이 포함되는 정도에 따라 설정되는 파라미터 가중 팩터
상기 타격음 구간 결정부는, 상기 파라미터 가중 팩터(βj)를 하기의 수학식에 의해 설정할 수 있다.
Figure 112016114802737-pat00008
β1 : 단구간 프레임의 에너지에 대한 파라미터 가중 팩터
β2 : 단구간 프레임의 스펙트럼 밀도에 대한 파라미터 가중 팩터
β3 : 단구간 프레임의 판별 추정값에 대한 파라미터 가중 팩터
Figure 112016114802737-pat00009
: 단구간 프레임의 에너지에 대한 이동 평균값
Figure 112016114802737-pat00010
: 단구간 프레임의 판별 추정값에 대한 이동 평균값
r1 : β1과 β3의 척도를 일치시키기 위한 비례 상수
상기 타격음 구간 결정부는, 기 설정된 개수의 단구간 프레임에서 연속적으로 상기 함수(R(n))의 출력 값이 기 설정된 출력값 이상인 것으로 나오는 경우, 해당 단구간 프레임들의 구간을 상기 타격음 구간의 시작인 것으로 판단하고, 상기 타격음 구간의 시작 이후에 존재하는 단구간 프레임에서 상기 함수(R(n))의 출력 값이 기 설정된 출력값 미만인 것으로 나오는 경우, 해당 단구간 프레임을 상기 타격음 구간의 종점인 것으로 판단할 수 있다.
본 발명의 일 실시예에 따른 방법은, 하나 이상의 프로세서들, 및 상기 하나 이상의 프로세서들에 의해 실행되는 하나 이상의 프로그램들을 저장하는 메모리를 구비한 컴퓨팅 장치에서 수행되는 방법으로서, 피 검사물에 대한 타격에 의해 발생하는 음향 신호를 획득하는 단계; 상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하는 단계; 상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하는 단계; 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하는 단계; 상기 각 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하는 단계; 상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정하는 단계; 및 상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는 단계를 포함한다.
상기 스펙트럼 밀도를 산출하는 단계는, 하기의 수학식을 통해 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density : SD)를 산출할 수 있다.
Figure 112016114802737-pat00011
M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
Ni : 기 설정된 잡음 신호의 i번째 MFCC 계수의 값
상기 판별 추정값을 산출하는 단계는, 하기의 수학식을 통해 상기 판별 추정값을 산출할 수 있다.
Figure 112016114802737-pat00012
Y : 해당 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값
M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
Wi : i번째 MFCC 계수에 대해 기 설정된 가중치
상기 임계값을 설정하는 단계는, 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출하는 단계; 및 상기 산출한 이동 평균값을 기반으로 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 임계값을 각각 설정하는 단계를 포함할 수 있다.
상기 타격음 구간을 결정하는 단계는, 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값과 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값 및 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 잡음 성분이 포함되는 정도에 따라 각각 설정되는 파라미터 가중 팩터를 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정할 수 있다.
상기 이동 평균값을 산출하는 단계는, 하기의 수학식을 통해 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출할 수 있다.
Figure 112016114802737-pat00013
p1(n)(즉, j=1인 경우) : n번째 단구간 프레임의 에너지
p2(n)(즉, j=2인 경우) : n번째 단구간 프레임의 스펙트럼 밀도
p3(n)(즉, j=3인 경우) : n번째 단구간 프레임의 판별 추정값
Figure 112016114802737-pat00014
: n번째 단구간 프레임의 에너지에 대한 이동 평균값
Figure 112016114802737-pat00015
: n번째 단구간 프레임의 스펙트럼 밀도에 대한 이동 평균값
Figure 112016114802737-pat00016
: n번째 단구간 프레임의 판별 추정값에 대한 이동 평균값
α : 이동 평균값의 변동량을 결정하는 상수
상기 타격음 구간을 결정하는 단계는, 하기의 수학식을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정할 수 있다.
Figure 112016114802737-pat00017
R(n) : n번째 단구간 프레임이 타격음을 포함하는지 여부를 판단하기 위한 함수
δj(n) : n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값과 n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값
N : 타격음 구간을 결정하기 위해 분석 대상으로 취하는 n번째 단구간 프레임과 인접한 단구간 프레임의 개수
m : n번째 단구간 프레임과 인접한 단구간 프레임의 번호
βj : 해당 단구간 프레임에 잡음 성분이 포함되는 정도에 따라 설정되는 파라미터 가중 팩터
상기 타격음 구간을 결정하는 단계는, 상기 파라미터 가중 팩터(βj)를 하기의 수학식에 의해 설정할 수 있다.
Figure 112016114802737-pat00018
β1 : 단구간 프레임의 에너지에 대한 파라미터 가중 팩터
β2 : 단구간 프레임의 스펙트럼 밀도에 대한 파라미터 가중 팩터
β3 : 단구간 프레임의 판별 추정값에 대한 파라미터 가중 팩터
Figure 112016114802737-pat00019
: 단구간 프레임의 에너지에 대한 이동 평균값
Figure 112016114802737-pat00020
: 단구간 프레임의 판별 추정값에 대한 이동 평균값
r1 : β1과 β3의 척도를 일치시키기 위한 비례 상수
상기 타격음 구간을 결정하는 단계는, 기 설정된 개수의 단구간 프레임에서 연속적으로 상기 함수(R(n))의 출력 값이 기 설정된 출력값 이상인 것으로 나오는 경우, 해당 단구간 프레임들의 구간을 상기 타격음 구간의 시작인 것으로 판단하는 단계; 및 상기 타격음 구간의 시작 이후에 존재하는 단구간 프레임에서 상기 함수(R(n))의 출력 값이 기 설정된 출력값 미만인 것으로 나오는 경우, 해당 단구간 프레임을 상기 타격음 구간의 종점인 것으로 판단하는 단계를 포함할 수 있다.
본 발명의 다른 실시예에 따른 컴퓨팅 장치는, 하나 이상의 프로세서들; 메모리; 및 하나 이상의 프로그램들을 포함하고, 상기 하나 이상의 프로그램들은 상기 메모리에 저장되고, 상기 하나 이상의 프로세서들에 의해 실행되도록 구성되며, 상기 하나 이상의 프로그램들은, 음향 신호를 획득하기 위한 명령; 상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하기 위한 명령; 상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하기 위한 명령; 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하기 위한 명령; 상기 각 단구간 프레임이 유효음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하기 위한 명령; 상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 유효음 구간의 경계를 결정하기 위한 임계값을 각각 설정하기 위한 명령; 및 상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 유효음 구간을 결정하기 위한 명령을 포함한다.
본 발명의 실시예에 의하면, 피 검사물의 타격에 의한 음향 신호를 획득하고, 음향 신호에서 특성 파라미터(단구간 프레임 에너지, 단구간 프레임의 스펙트럼 밀도, 단구간 프레임의 판별 추정값)을 산출하고, 이를 기반으로 타격음 구간을 검출함으로써, 잡음 환경에서도 신뢰성 있는 타격음 구간을 검출 할 수 있게 된다. 그리고, 타격음 구간을 자동으로 검출함으로써, 피 검사물의 종류에 따라 타격음 지속 시간을 고려한 설정값 등을 별도로 설정할 필요가 없으며, 다품종 소량 생산 제품의 결함 여부를 신속하고 용이하게 검사할 수 있게 된다.
도 1은 본 발명의 일 실시예에 따른 비파괴 검사 장치의 구성을 나타낸 블록도
도 2는 본 발명의 일 실시예에 따른 비파괴 검사 방법을 나타낸 흐름도
도 3은 예시적인 실시예들에서 사용되기에 적합한 컴퓨팅 장치를 포함하는 컴퓨팅 환경을 예시하여 설명하기 위한 블록도
이하, 도면을 참조하여 본 발명의 구체적인 실시형태를 설명하기로 한다. 이하의 상세한 설명은 본 명세서에서 기술된 방법, 장치 및/또는 시스템에 대한 포괄적인 이해를 돕기 위해 제공된다. 그러나 이는 예시에 불과하며 본 발명은 이에 제한되지 않는다.
본 발명의 실시예들을 설명함에 있어서, 본 발명과 관련된 공지기술에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우에는 그 상세한 설명을 생략하기로 한다. 그리고, 후술되는 용어들은 본 발명에서의 기능을 고려하여 정의된 용어들로서 이는 사용자, 운용자의 의도 또는 관례 등에 따라 달라질 수 있다. 그러므로 그 정의는 본 명세서 전반에 걸친 내용을 토대로 내려져야 할 것이다. 상세한 설명에서 사용되는 용어는 단지 본 발명의 실시예들을 기술하기 위한 것이며, 결코 제한적이어서는 안 된다. 명확하게 달리 사용되지 않는 한, 단수 형태의 표현은 복수 형태의 의미를 포함한다. 본 설명에서, "포함" 또는 "구비"와 같은 표현은 어떤 특성들, 숫자들, 단계들, 동작들, 요소들, 이들의 일부 또는 조합을 가리키기 위한 것이며, 기술된 것 이외에 하나 또는 그 이상의 다른 특성, 숫자, 단계, 동작, 요소, 이들의 일부 또는 조합의 존재 또는 가능성을 배제하도록 해석되어서는 안 된다.
이하의 설명에 있어서, 신호 또는 정보의 "전송", "통신", "송신", "수신" 기타 이와 유사한 의미의 용어는 일 구성요소에서 다른 구성요소로 신호 또는 정보가 직접 전달되는 것뿐만이 아니라 다른 구성요소를 거쳐 전달되는 것도 포함한다. 특히 신호 또는 정보를 일 구성요소로 "전송" 또는 "송신"한다는 것은 그 신호 또는 정보의 최종 목적지를 지시하는 것이고 직접적인 목적지를 의미하는 것이 아니다. 이는 신호 또는 정보의 "수신"에 있어서도 동일하다. 또한 본 명세서에 있어서, 2 이상의 데이터 또는 정보가 "관련"된다는 것은 하나의 데이터(또는 정보)를 획득하면, 그에 기초하여 다른 데이터(또는 정보)의 적어도 일부를 획득할 수 있음을 의미한다.
또한, 제1, 제2 등의 용어는 다양한 구성 요소들을 설명하는데 사용될 수 있지만, 상기 구성 요소들은 상기 용어들에 의해 한정되어서는 안 된다. 상기 용어들은 하나의 구성 요소를 다른 구성 요소로부터 구별하는 목적으로 사용될 수 있다. 예를 들어, 본 발명의 권리 범위를 벗어나지 않으면서 제1 구성 요소는 제2 구성 요소로 명명될 수 있고, 유사하게 제2 구성 요소도 제1 구성 요소로 명명될 수 있다.
한편, 본 발명의 실시예는 본 명세서에서 기술한 방법들을 컴퓨터상에서 수행하기 위한 프로그램, 및 상기 프로그램을 포함하는 컴퓨터 판독 가능 기록매체를 포함할 수 있다. 상기 컴퓨터 판독 가능 기록매체는 프로그램 명령, 로컬 데이터 파일, 로컬 데이터 구조 등을 단독으로 또는 조합하여 포함할 수 있다. 상기 매체는 본 발명을 위하여 특별히 설계되고 구성된 것들이거나, 또는 컴퓨터 소프트웨어 분야에서 통상적으로 사용 가능한 것일 수 있다. 컴퓨터 판독 가능 기록매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체, CD-ROM, DVD와 같은 광 기록 매체, 및 롬, 램, 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. 상기 프로그램의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 코드뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 코드를 포함할 수 있다.
도 1은 본 발명의 일 실시예에 따른 비파괴 검사 장치를 나타낸 블록도이다.
도 1을 참조하면, 비파괴 검사 장치(100)는 타격부(102), 음향 신호 획득부(104), 전처리부(106), 단구간 프레임 에너지 산출부(108), 스펙트럼 밀도 산출부(110), 판별 추정값 산출부(112), 임계값 설정부(114), 타격음 구간 결정부(116), 및 결함 여부 분석부(118)를 포함할 수 있다.
타격부(102)는 피 검사물(미도시)을 타격하는 역할을 한다. 예시적인 실시예에서, 타격부(102)는 미리 설정된 주파수의 진동을 발생시키는 가진기(加振機)(102a) 및 가진기에서 발생되는 진동을 전달받아 피 검사물(미도시)을 타격하는 타격기(102b)를 포함할 수 있다.
음향 신호 획득부(104)는 타격부(102)가 피 검사물(미도시)을 타격하였을 때 발생하는 음향 신호를 획득하는 역할을 한다. 음향 신호 획득부(104)는 마이크로폰(Microphone)을 포함하여 구현될 수 있다. 음향 신호 획득부(104)에서 획득한 음향 신호는 아날로그 신호일 수 있다.
전처리부(106)는 음향 신호 획득부(104)가 획득한 음향 신호를 샘플링하여 디지털 신호인 오디오 데이터로 변환하고, 변환된 오디오 데이터를 기 설정된 단구간 길이(예를 들어, 5 ~ 60ms)를 갖는 복수 개의 프레임(이하, 단구간 프레임으로 지칭할 수 있음)으로 구분할 수 있다. 여기서, 각 단구간 프레임은 이웃하는 단구간 프레임과 일정 범위가 겹쳐지게(즉, 오버랩(Overlap) 되게) 설정될 수 있다. 예를 들어, 단구간 프레임이 400개의 샘플을 포함하는 경우, 첫 번째 단구간 프레임은 0 ~ 400 번째 샘플까지로 하고, 두 번째 단구간 프레임은 160 ~ 560번째 샘플까지로 설정할 수 있다. 이와 같이, 음향 신호(즉, 오디오 신호)를 샘플링하여 기 설정된 단구간 길이의 복수 개의 프레임으로 구분하여 분석하는 방법을 단구간 분석 방법이라 한다.
단구간 프레임 에너지 산출부(108)는 음향 신호 획득부(104)가 획득한 음향 신호의 각 단구간 프레임 에너지를 산출할 수 있다. 즉, 단구간 프레임 에너지 산출부(108)는 전처리부(106)에서 구분된 각 단구간 프레임의 에너지를 산출할 수 있다. 단구간 프레임 에너지는 해당 단구간에 속하는 음향 신호가 내포하고 있는 에너지의 합을 의미할 수 있다. 단구간 프레임 에너지 산출부(108)는 하기 수학식 1을 통해 각 단구간 프레임 에너지를 산출할 수 있다.
(수학식 1)
Figure 112016114802737-pat00021
여기서, En은 n번째 단구간 프레임의 에너지를 나타내고, Xn(k)는 n번째 단구간 프레임의 k번째 샘플링에 대응하는 음향 신호를 나타내며, N은 해당 단구간 프레임에 포함되는 샘플링의 총 개수를 나타낼 수 있다.
스펙트럼 밀도 산출부(110)는 전처리부(106)에서 구분된 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density : SD)를 산출할 수 있다. 여기서, 스펙트럼 밀도는 해당 단구간 프레임에 대응하는 음향 신호에서 잡음 성분의 추정치를 차감하여 산출할 수 있다. 보다 상세하게는, 소정 단구간 프레임의 스펙트럼 밀도는 해당 단구간 프레임의 MFCC(Mel Frequency Cepstral Coefficient)에서 잡음의 추정치의 MFCC를 차감하여 산출할 수 있다. 여기서, 잡음의 추정치는 일반적인 잡음 신호(예를 들어, 화이트 노이즈 등)에서 MFCC 값을 구하여 산출할 수 있다.
스펙트럼 밀도 산출부(110)는 각 단구간 프레임의 스펙트럼 밀도(SD)를 하기의 수학식 2를 통해 산출할 수 있다.
(수학식 2)
Figure 112016114802737-pat00022
여기서, M은 MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수를 나타내고, Xi는 해당 단구간 프레임의 i번째 MFCC 계수의 값을 나타내며, Ni는 기 설정된 잡음 신호의 i번째 MFCC 계수의 값을 나타낼 수 있다. 수학식 2에 의하면, 해당 단구간 프레임의 음향 신호에서 잡음으로 추정되는 성분이 차감되는 효과를 얻을 수 있다.
단구간 프레임의 MFCC(Mel Frequency Cepstral Coefficient)를 구하는 방법을 간략히 살펴보면, 단구간 프레임을 FFT(Fast Fourier Transform)을 통해 주파수 영역으로 변환하고, 단구간 프레임의 주파수 영역을 복수 개의 필터 뱅크(Filter Bank)(예를 들어, Mel-spaced filter bank)에 적용하여 각 필터 뱅크에서의 에너지를 구한다. 여기서, 필터 뱅크의 개수는 20 ~ 40개가 적용될 수 있으나, 표준적으로는 26개가 사용될 수 있다. 그리고, 각 필터 뱅크에서의 에너지에 로그를 취한 후 DCT(Discrete Cosine Transform)을 적용하면 해당 단구간 프레임에 대한 MFCC(Mel Frequency Cepstral Coefficient)를 산출할 수 있다. MFCC를 산출하는 방법은 이미 공지된 기술이므로 이에 대한 보다 자세한 설명은 생략하기로 한다.
판별 추정값 산출부(112)는 전처리부(106)에서 구분된 각 단구간 프레임이 타격음 구간에 포함되는지와 관련한 판별 추정값을 산출할 수 있다. 즉, 판별 추정값 산출부(110)는 각 단구간 프레임이 타격음 구간에 포함되는지 아니면 잡음 또는 무음 구간에 포함되는지 여부에 따라 서로 다른 판별 추정값을 출력할 수 있다. 판별 추정값 산출부(112)는 각 단구간 프레임이 타격음 구간에 포함되는지와 관련하여 하기의 수학식 3에 의한 선형 판별 함수를 통해 판별 추정값을 산출할 수 있다.
(수학식 3)
Figure 112016114802737-pat00023
여기서, Y는 해당 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값(즉, 출력값)을 나타내고, Xi는 해당 단구간 프레임의 i번째 MFCC 계수의 값을 나타내며, Wi는 i번째 MFCC 계수에 대한 가중치(선형 판별 함수의 가중치)를 나타내며, V는 바이어스(Bias)로서 상수값을 나타낼 수 있다.
수학식 3에서 W는 가중치 벡터로서, 가중치 벡터(W)는 입력 벡터인 X가 타격음 구간에 포함되면 출력 신호(Y)가 1이고, 입력 벡터인 X가 잡음 구간 또는 무음 구간에 포함되면 출력 신호(Y)가 0이 되도록 미리 학습된 값일 수 있다. 여기서, 가중치 벡터(W)는 하기의 수학식 4와 같은 선형 회귀법에 의해 구할 수 있다.
(수학식 4)
Figure 112016114802737-pat00024
여기서, T는 전이 행렬(Transition Matrix)를 나타낼 수 있다.
임계값 설정부(114)는 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정할 수 있다. 즉, 음향 신호 획득부(104)에서 획득한 음향 신호에는 타격음 구간뿐만 아니라 잡음 구간 또는 무음 구간이 포함되어 있으므로, 타격음 구간과 그 이외의 구간(즉, 잡음 구간 또는 무음 구간)을 구분하기 위한 경계를 결정하기 위해 임계값을 설정할 수 있다.
구체적으로, 임계값 설정부(114)는 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값을 각각 산출할 수 있다. 임계값 설정부(114)는 하기의 수학식 5를 통해 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값을 각각 산출할 수 있다.
(수학식 5)
Figure 112016114802737-pat00025
여기서, p1(n)(즉, j=1인 경우)은 n번째 단구간 프레임의 에너지를 나타내고, p2(n)(즉, j=2인 경우)은 n번째 단구간 프레임의 스펙트럼 밀도를 나타내며, p3(n)(즉, j=3인 경우)은 n번째 단구간 프레임의 판별 추정값을 나타내고, α는 이동 평균값의 변동량을 결정하는 상수를 나타낼 수 있다.
그러면,
Figure 112016114802737-pat00026
은 n번째 단구간 프레임의 에너지에 대한 이동 평균값을 나타내고,
Figure 112016114802737-pat00027
는 n번째 단구간 프레임의 스펙트럼 밀도에 대한 이동 평균값을 나타내며,
Figure 112016114802737-pat00028
는 n번째 단구간 프레임의 판별 추정값에 대한 이동 평균값을 나타내게 된다.
그리고, 임계값 설정부(114)는 하기의 수학식 6을 통해 각 단구간 프레임의 에너지에 대한 임계값(T1(n)), 단구간 프레임의 스펙트럼 밀도에 대한 임계값(T2(n)), 및 단구간 프레임의 판별 추정값에 대해 임계값(T3(n))을 설정할 수 있다.
(수학식 6)
Figure 112016114802737-pat00029
여기서, a는 임계값 설정을 위한 가중치를 나타내고, b는 바이어스(Bias)로서 상수값을 나타낼 수 있다. a 및 b는 수학식 4와 유사한 방식의 선형 회귀법을 통해 구할 수 있다.
타격음 구간 결정부(116)는 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값과 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 비교하여 나온 값(이하, 비교 결과값이라 지칭할 수 있음)과 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 각각 설정된 파라미터 가중 팩터를 기반으로 타격음 구간을 결정할 수 있다.
구체적으로, 타격음 구간 결정부(116)는 하기의 수학식 7을 통해 음향 신호 획득부(104)가 획득한 음향 신호에서 타격음 구간을 결정할 수 있다. 예시적인 실시예에서, 타격음 구간 결정부(116)는 수학식 7에서 R(n)의 출력값이 0 이상의 값을 갖는 단구간 프레임이 연속적으로 기 설정된 개수 이상 나타나는지 여부에 따라 타격음 구간을 결정할 수 있다.
(수학식 7)
Figure 112016114802737-pat00030
여기서, R(n)은 n번째 단구간 프레임이 타격음을 포함하는지 여부를 판단하기 위한 함수(이하, 타격음 판단 함수라 지칭할 수 있음)이다. δj(n)은 상기 비교 결과값을 나타낸 것으로,
Figure 112016114802737-pat00031
이면 n번째 단구간 프레임이 타격음을 포함하는 것으로 잠정적인 판단을 하여 δj(n) = 1로 하고,
Figure 112016114802737-pat00032
이면 n번째 단구간 프레임이 타격음을 포함하지 않는 것으로 잠정적인 판단을 하여 δj(n) = -1로 할 수 있다. N은 타격음 구간을 결정하기 위해 분석 대상으로 취하는 n번째 단구간 프레임과 인접한 단구간 프레임의 개수를 나타낼 수 있다. m은 상기 n번째 단구간 프레임과 인접한 단구간 프레임의 번호를 나타낼 수 있다.
βj(j = 1, 2, 3)는 해당 단구간 프레임에 잡음 성분이 포함되는 정도에 따라 설정되는 파라미터 가중 팩터로서, β1은 해당 단구간 프레임의 에너지에 대한 파라미터 가중 팩터이고, β2은 해당 단구간 프레임의 스펙트럼 밀도에 대한 파라미터 가중 팩터이며, β3은 해당 단구간 프레임의 판별 추정값에 대한 파라미터 가중 팩터일 수 있다. βj(j = 1, 2, 3)는 아래의 수학식 8과 같이 설정될 수 있다.
(수학식 8)
Figure 112016114802737-pat00033
여기서, r1은 β1과 β3의 척도를 일치시키기 위한 비례 상수이다.
β1은 단구간 프레임의 에너지가 잡음 성분이 클수록 증가하기 때문에, 잡음 성분이 적은 경우에도 타격음 검출을 수행할 수 있도록
Figure 112016114802737-pat00034
(즉, 단구간 프레임의 에너지에 대한 이동 평균값)에 역수를 취하여 구할 수 있다. β2은 주변 잡음 성분을 고려하여 미리 계산된 적절한 상수로 설정할 수 있다. 예를 들어, β2은 0.1이 사용될 수 있다. β3은 0과 1사이의 값으로, 단구간 프레임의 판별 추정값이 잡음 성분이 클수록 감소하기 때문에, 잡음 성분이 큰 경우에도 타격음 검출을 수행할 수 있도록
Figure 112016114802737-pat00035
(즉, 단구간 프레임의 판별 추정값에 대한 이동 평균값) 및 최대값인 1과의 차이값을 사용할 수 있다.
한편, 타격음으로 의미 있는 구간은 최소 40ms 이상 지속되어야 한다. 예를 들어, 단구간 프레임의 시구간 길이가 8ms인 경우, 타격음 구간으로 결정하기 위해서는 최소 5개의 단구간 프레임이 연속적으로 타격음을 포함하는 것으로 판단되어야 한다.
따라서, 타격음 구간 결정부(116)는 기 설정된 개수의 단구간 프레임에서 연속적으로 R(n)의 출력값이 0 이상으로 판별되는 경우, 타격음 구간의 시작인 것으로 판단할 수 있다. 반면, 기 설정된 개수의 단구간 프레임에서 연속적으로 R(n)의 출력값이 0 이상이 되지 않는 경우, 해당 단구간 프레임들에 대해 잡음 또는 무음인 것으로 판단할 수 있다. 또한, 타격음 구간 결정부(116)는 타격음 구간이 시작된 이후, R(n)의 출력값이 0 미만으로 나오는 경우, 해당 단구간 프레임을 타격음 구간의 종점인 것으로 판단할 수 있다.
결함 여부 분석부(118)는 타격음 구간으로 결정된 음향 신호를 분석하여 피 검사물에 결함이 있는지 여부를 분석할 수 있다. 예시적인 실시예에서, 결함 여부 분석부(118)는 상기 음향 신호를 기 설정된 음향 신호 패턴과 비교하여 결함 여부를 판단할 수 있다.
본 발명의 실시예에 의하면, 피 검사물의 타격에 의한 음향 신호를 획득하고, 음향 신호에서 특성 파라미터(단구간 프레임 에너지, 단구간 프레임의 스펙트럼 밀도, 단구간 프레임의 판별 추정값)을 산출하고, 이를 기반으로 타격음 구간을 검출함으로써, 잡음 환경에서도 신뢰성 있는 타격음 구간을 검출 할 수 있게 된다. 그리고, 타격음 구간을 자동으로 검출함으로써, 피 검사물의 종류에 따라 타격음 지속 시간을 고려한 설정값 등을 별도로 설정할 필요가 없으며, 다품종 소량 생산 제품의 결함 여부를 신속하고 용이하게 검사할 수 있게 된다.
도 2는 본 발명의 일 실시예에 따른 비파괴 검사 방법을 설명하기 위한 흐름도이다. 도 2에 도시된 방법은 예를 들어, 전술한 비파과 검사 장치(100)에 의해 수행될 수 있다. 도시된 흐름도에서는 상기 방법을 복수 개의 단계로 나누어 기재하였으나, 적어도 일부의 단계들은 순서를 바꾸어 수행되거나, 다른 단계와 결합되어 함께 수행되거나, 생략되거나, 세부 단계들로 나뉘어 수행되거나, 또는 도시되지 않은 하나 이상의 단계가 부가되어 수행될 수 있다.
도 2를 참조하면, 비파괴 검사 장치(100)는 피 검사물(미도시)을 타격한 후(S 101), 타격에 따른 음향 신호를 획득한다(S 103).
다음으로, 비파괴 검사 장치(100)는 획득한 음향 신호를 샘플링하여 디지털 신호인 오디오 데이터로 변환하고, 변환된 오디오 데이터를 복수 개의 단구간 프레임으로 구분한다(S 105).
다음으로, 비파괴 검사 장치(100)는 각 단구간 프레임을 대상으로 단구간 프레임 에너지를 산출한다(S 107). 비파괴 검사 장치(100)는 상기 수학식 1을 통해 각 단구간 프레임의 에너지를 산출할 수 있다.
다음으로, 비파괴 검사 장치(100)는 각 단구간 프레임의 스펙트럼 밀도를 산출한다(S 109). 비파괴 검사 장치(100)는 상기 수학식 2를 통해 각 단구간 프레임의 스펙트럼 밀도를 산출할 수 있다.
다음으로, 비파괴 검사 장치(100)는 각 단구간 프레임이 타격음 구간에 포함되는지 여부에 대한 판별 추정값을 산출한다(S 111). 비파괴 검사 장치(100)는 상기 수학식 3 및 수학식 4를 통해 각 단구간 프레임의 판별 추정값을 산출할 수 있다.
다음으로, 비파괴 검사 장치(100)는 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정한다(S 113). 구체적으로, 비파괴 검사 장치(100)는 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값을 각각 산출한 후, 이를 이용하여 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 각각 설정할 수 있다. 비파괴 검사 장치(100)는 상기 수학식 5 및 수학식 6을 통해 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 각각 설정할 수 있다.
다음으로, 비파괴 검사 장치(100)는 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값과 각 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 비교하여 나온 값(즉, 비교 결과값)과 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 각각 설정된 파라미터 가중 팩터를 기반으로 타격음 구간을 결정한다(S 115). 비파괴 검사 장치(100)는 상기 수학식 7에서 R(n)의 출력값이 0 이상의 값을 갖는 단구간 프레임이 연속적으로 기 설정된 개수 이상 나타나는지 여부에 따라 타격음 구간을 결정할 수 있다.
다음으로, 비파괴 검사 장치(100)는 타격음 구간으로 결정된 음향 신호를 기반으로 피 검사물의 결함 여부를 분석한다(S 117).
여기서는, 상기 기술한 음향 분석 방법이 비파괴 검사 방법에 적용되는 것을 일 실시예로 설명하였으나, 본 발명이 이에 한정되는 것은 아니며 음향 신호를 이용한 다양한 기술 분야에 적용될 수 있음은 물론이다.
도 3은 예시적인 실시예들에서 사용되기에 적합한 컴퓨팅 장치를 포함하는 컴퓨팅 환경(10)을 예시하여 설명하기 위한 블록도이다. 도시된 실시예에서, 각 컴포넌트들은 이하에 기술된 것 이외에 상이한 기능 및 능력을 가질 수 있고, 이하에 기술되는 것 이외에도 추가적인 컴포넌트를 포함할 수 있다.
도시된 컴퓨팅 환경(10)은 컴퓨팅 장치(12)를 포함한다. 일 실시예에서, 컴퓨팅 장치(12)는 비파괴 검사를 위한 장치(예를 들어, 비파괴 검사 장치(100))일 수 있다.
컴퓨팅 장치(12)는 적어도 하나의 프로세서(14), 컴퓨터 판독 가능 저장 매체(16) 및 통신 버스(18)를 포함한다. 프로세서(14)는 컴퓨팅 장치(12)로 하여금 앞서 언급된 예시적인 실시예에 따라 동작하도록 할 수 있다. 예컨대, 프로세서(14)는 컴퓨터 판독 가능 저장 매체(16)에 저장된 하나 이상의 프로그램들을 실행할 수 있다. 상기 하나 이상의 프로그램들은 하나 이상의 컴퓨터 실행 가능 명령어를 포함할 수 있으며, 상기 컴퓨터 실행 가능 명령어는 프로세서(14)에 의해 실행되는 경우 컴퓨팅 장치(12)로 하여금 예시적인 실시예에 따른 동작들을 수행하도록 구성될 수 있다.
컴퓨터 판독 가능 저장 매체(16)는 컴퓨터 실행 가능 명령어 내지 프로그램 코드, 프로그램 데이터 및/또는 다른 적합한 형태의 정보를 저장하도록 구성된다. 컴퓨터 판독 가능 저장 매체(16)에 저장된 프로그램(20)은 프로세서(14)에 의해 실행 가능한 명령어의 집합을 포함한다. 일 실시예에서, 컴퓨터 판독 가능 저장 매체(16)는 메모리(랜덤 액세스 메모리와 같은 휘발성 메모리, 비휘발성 메모리, 또는 이들의 적절한 조합), 하나 이상의 자기 디스크 저장 디바이스들, 광학 디스크 저장 디바이스들, 플래시 메모리 디바이스들, 그 밖에 컴퓨팅 장치(12)에 의해 액세스되고 원하는 정보를 저장할 수 있는 다른 형태의 저장 매체, 또는 이들의 적합한 조합일 수 있다.
통신 버스(18)는 프로세서(14), 컴퓨터 판독 가능 저장 매체(16)를 포함하여 컴퓨팅 장치(12)의 다른 다양한 컴포넌트들을 상호 연결한다.
컴퓨팅 장치(12)는 또한 하나 이상의 입출력 장치(24)를 위한 인터페이스를 제공하는 하나 이상의 입출력 인터페이스(22) 및 하나 이상의 네트워크 통신 인터페이스(26)를 포함할 수 있다. 입출력 인터페이스(22) 및 네트워크 통신 인터페이스(26)는 통신 버스(18)에 연결된다. 입출력 장치(24)는 입출력 인터페이스(22)를 통해 컴퓨팅 장치(12)의 다른 컴포넌트들에 연결될 수 있다. 예시적인 입출력 장치(24)는 포인팅 장치(마우스 또는 트랙패드 등), 키보드, 터치 입력 장치(터치패드 또는 터치스크린 등), 음성 또는 소리 입력 장치, 다양한 종류의 센서 장치 및/또는 촬영 장치와 같은 입력 장치, 및/또는 디스플레이 장치, 프린터, 스피커 및/또는 네트워크 카드와 같은 출력 장치를 포함할 수 있다. 예시적인 입출력 장치(24)는 컴퓨팅 장치(12)를 구성하는 일 컴포넌트로서 컴퓨팅 장치(12)의 내부에 포함될 수도 있고, 컴퓨팅 장치(12)와는 구별되는 별개의 장치로 컴퓨팅 장치(12)와 연결될 수도 있다.
이상에서 본 발명의 대표적인 실시예들을 상세하게 설명하였으나, 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자는 상술한 실시예에 대하여 본 발명의 범주에서 벗어나지 않는 한도 내에서 다양한 변형이 가능함을 이해할 것이다. 그러므로 본 발명의 권리범위는 설명된 실시예에 국한되어 정해져서는 안 되며, 후술하는 특허청구범위뿐만 아니라 이 특허청구범위와 균등한 것들에 의해 정해져야 한다.
100 : 비파괴 검사 장치
102 : 타격부
104 : 음향 신호 획득부
106 : 전처리부
108 : 단구간 프레임 에너지 산출부
110 : 스펙트럼 밀도 산출부
112 : 판별 추정값 산출부
114 : 임계값 설정부
116 : 타격음 구간 결정부
118 : 결함 여부 분석부

Claims (19)

  1. 피 검사물에 대한 타격에 의해 발생하는 음향 신호를 획득하는 음향 신호 획득부;
    상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하는 전처리부;
    상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하는 단구간 프레임 에너지 산출부;
    상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하는 스펙트럼 밀도 산출부;
    상기 각 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하는 판별 추정값 산출부;
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정하는 임계값 설정부;
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는 타격음 구간 결정부; 및
    상기 결정된 타격음 구간에 대응하는 음향 신호를 기반으로 상기 피 검사물의 결함 여부를 분석하는 결함 여부 분석부를 포함하고,
    상기 스펙트럼 밀도 산출부는,
    하기의 수학식을 통해 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density : SD를 산출하는, 비파괴 검사 장치.
    Figure 112017092210556-pat00059

    M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
    Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
    Ni : 기 설정된 잡음 신호의 i번째 MFCC 계수의 값
  2. 삭제
  3. 피 검사물에 대한 타격에 의해 발생하는 음향 신호를 획득하는 음향 신호 획득부;
    상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하는 전처리부;
    상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하는 단구간 프레임 에너지 산출부;
    상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하는 스펙트럼 밀도 산출부;
    상기 각 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하는 판별 추정값 산출부;
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정하는 임계값 설정부;
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는 타격음 구간 결정부; 및
    상기 결정된 타격음 구간에 대응하는 음향 신호를 기반으로 상기 피 검사물의 결함 여부를 분석하는 결함 여부 분석부를 포함하고,
    상기 판별 추정값 산출부는,
    하기의 수학식을 통해 상기 판별 추정값을 산출하는, 비파괴 검사 장치.
    Figure 112017092210556-pat00037

    Y : 해당 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값
    M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
    Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
    Wi : i번째 MFCC 계수에 대해 기 설정된 가중치
    V : 바이어스(Bias)로서 기 설정된 상수값
  4. 청구항 1에 있어서,
    상기 임계값 설정부는,
    상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출하고, 상기 산출한 이동 평균값을 기반으로 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 임계값을 각각 설정하는, 비파괴 검사 장치.
  5. 청구항 4에 있어서,
    상기 타격음 구간 결정부는,
    상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값과 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값 및 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 잡음 성분이 포함되는 정도에 따라 각각 설정되는 파라미터 가중 팩터를 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는, 비파괴 검사 장치.
  6. 청구항 4에 있어서,
    상기 임계값 설정부는,
    하기의 수학식을 통해 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출하는, 비파괴 검사 장치.
    Figure 112016114802737-pat00038

    p1(n)(즉, j=1인 경우) : n번째 단구간 프레임의 에너지
    p2(n)(즉, j=2인 경우) : n번째 단구간 프레임의 스펙트럼 밀도
    p3(n)(즉, j=3인 경우) : n번째 단구간 프레임의 판별 추정값
    Figure 112016114802737-pat00039
    : n번째 단구간 프레임의 에너지에 대한 이동 평균값
    Figure 112016114802737-pat00040
    : n번째 단구간 프레임의 스펙트럼 밀도에 대한 이동 평균값
    Figure 112016114802737-pat00041
    : n번째 단구간 프레임의 판별 추정값에 대한 이동 평균값
    α : 이동 평균값의 변동량을 결정하는 상수
  7. 청구항 6에 있어서,
    상기 타격음 구간 결정부는,
    하기의 수학식을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는, 비파괴 검사 장치.
    Figure 112016114802737-pat00042

    R(n) : n번째 단구간 프레임이 타격음을 포함하는지 여부를 판단하기 위한 함수
    δj(n) : n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값과 n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값
    N : 타격음 구간을 결정하기 위해 분석 대상으로 취하는 n번째 단구간 프레임과 인접한 단구간 프레임의 개수
    m : n번째 단구간 프레임과 인접한 단구간 프레임의 번호
    βj : 해당 단구간 프레임에 잡음 성분이 포함되는 정도에 따라 설정되는 파라미터 가중 팩터
  8. 청구항 7에 있어서,
    상기 타격음 구간 결정부는,
    상기 파라미터 가중 팩터(βj)를 하기의 수학식에 의해 설정하는, 비파괴 검사 장치.
    Figure 112016114802737-pat00043

    β1 : 단구간 프레임의 에너지에 대한 파라미터 가중 팩터
    β2 : 단구간 프레임의 스펙트럼 밀도에 대한 파라미터 가중 팩터
    β3 : 단구간 프레임의 판별 추정값에 대한 파라미터 가중 팩터
    Figure 112016114802737-pat00044
    : 단구간 프레임의 에너지에 대한 이동 평균값
    Figure 112016114802737-pat00045
    : 단구간 프레임의 판별 추정값에 대한 이동 평균값
    r1 : β1과 β3의 척도를 일치시키기 위한 비례 상수
  9. 청구항 7에 있어서,
    상기 타격음 구간 결정부는,
    기 설정된 개수의 단구간 프레임에서 연속적으로 상기 함수(R(n))의 출력 값이 기 설정된 출력값 이상인 것으로 나오는 경우, 해당 단구간 프레임들의 구간을 상기 타격음 구간의 시작인 것으로 판단하고, 상기 타격음 구간의 시작 이후에 존재하는 단구간 프레임에서 상기 함수(R(n))의 출력 값이 기 설정된 출력값 미만인 것으로 나오는 경우, 해당 단구간 프레임을 상기 타격음 구간의 종점인 것으로 판단하는, 비파괴 검사 장치.
  10. 하나 이상의 프로세서들, 및
    상기 하나 이상의 프로세서들에 의해 실행되는 하나 이상의 프로그램들을 저장하는 메모리를 구비한 컴퓨팅 장치에서 수행되는 방법으로서,
    피 검사물에 대한 타격에 의해 발생하는 음향 신호를 획득하는 단계;
    상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하는 단계;
    상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하는 단계;
    상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하는 단계;
    상기 각 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하는 단계;
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정하는 단계; 및
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는 단계를 포함하고,
    상기 스펙트럼 밀도를 산출하는 단계는,
    하기의 수학식을 통해 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density : SD를 산출하는, 방법.
    Figure 112017092210556-pat00060

    M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
    Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
    Ni : 기 설정된 잡음 신호의 i번째 MFCC 계수의 값
  11. 삭제
  12. 하나 이상의 프로세서들, 및
    상기 하나 이상의 프로세서들에 의해 실행되는 하나 이상의 프로그램들을 저장하는 메모리를 구비한 컴퓨팅 장치에서 수행되는 방법으로서,
    피 검사물에 대한 타격에 의해 발생하는 음향 신호를 획득하는 단계;
    상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하는 단계;
    상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하는 단계;
    상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하는 단계;
    상기 각 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하는 단계;
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 타격음 구간의 경계를 결정하기 위한 임계값을 각각 설정하는 단계; 및
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는 단계를 포함하고,
    상기 판별 추정값을 산출하는 단계는,
    하기의 수학식을 통해 상기 판별 추정값을 산출하는, 방법.
    Figure 112017092210556-pat00047

    Y : 해당 단구간 프레임이 타격음 구간에 포함되는지 여부와 관련한 판별 추정값
    M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
    Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
    Wi : i번째 MFCC 계수에 대해 기 설정된 가중치
  13. 청구항 10에 있어서,
    상기 임계값을 설정하는 단계는,
    상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출하는 단계; 및
    상기 산출한 이동 평균값을 기반으로 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 임계값을 각각 설정하는 단계를 포함하는, 방법.
  14. 청구항 13에 있어서,
    상기 타격음 구간을 결정하는 단계는,
    상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값과 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값 및 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대해 잡음 성분이 포함되는 정도에 따라 각각 설정되는 파라미터 가중 팩터를 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는, 방법.
  15. 청구항 13에 있어서,
    상기 이동 평균값을 산출하는 단계는,
    하기의 수학식을 통해 상기 단구간 프레임 에너지, 상기 스펙트럼 밀도, 및 상기 판별 추정값에 대한 이동 평균값을 각각 산출하는, 방법.
    Figure 112016114802737-pat00048

    p1(n)(즉, j=1인 경우) : n번째 단구간 프레임의 에너지
    p2(n)(즉, j=2인 경우) : n번째 단구간 프레임의 스펙트럼 밀도
    p3(n)(즉, j=3인 경우) : n번째 단구간 프레임의 판별 추정값
    Figure 112016114802737-pat00049
    : n번째 단구간 프레임의 에너지에 대한 이동 평균값
    Figure 112016114802737-pat00050
    : n번째 단구간 프레임의 스펙트럼 밀도에 대한 이동 평균값
    Figure 112016114802737-pat00051
    : n번째 단구간 프레임의 판별 추정값에 대한 이동 평균값
    α : 이동 평균값의 변동량을 결정하는 상수
  16. 청구항 15에 있어서,
    상기 타격음 구간을 결정하는 단계는,
    하기의 수학식을 기반으로 상기 획득한 음향 신호에서 타격음 구간을 결정하는, 방법.
    Figure 112016114802737-pat00052

    R(n) : n번째 단구간 프레임이 타격음을 포함하는지 여부를 판단하기 위한 함수
    δj(n) : n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 이동 평균값과 n번째 단구간 프레임의 에너지, 스펙트럼 밀도, 및 판별 추정값에 대한 임계값을 비교하여 나온 비교 결과값
    N : 타격음 구간을 결정하기 위해 분석 대상으로 취하는 n번째 단구간 프레임과 인접한 단구간 프레임의 개수
    m : n번째 단구간 프레임과 인접한 단구간 프레임의 번호
    βj : 해당 단구간 프레임에 잡음 성분이 포함되는 정도에 따라 설정되는 파라미터 가중 팩터
  17. 청구항 16에 있어서,
    상기 타격음 구간을 결정하는 단계는,
    상기 파라미터 가중 팩터(βj)를 하기의 수학식에 의해 설정하는, 방법.
    Figure 112016114802737-pat00053

    β1 : 단구간 프레임의 에너지에 대한 파라미터 가중 팩터
    β2 : 단구간 프레임의 스펙트럼 밀도에 대한 파라미터 가중 팩터
    β3 : 단구간 프레임의 판별 추정값에 대한 파라미터 가중 팩터
    Figure 112016114802737-pat00054
    : 단구간 프레임의 에너지에 대한 이동 평균값
    Figure 112016114802737-pat00055
    : 단구간 프레임의 판별 추정값에 대한 이동 평균값
    r1 : β1과 β3의 척도를 일치시키기 위한 비례 상수
  18. 청구항 16에 있어서,
    상기 타격음 구간을 결정하는 단계는,
    기 설정된 개수의 단구간 프레임에서 연속적으로 상기 함수(R(n))의 출력 값이 기 설정된 출력값 이상인 것으로 나오는 경우, 해당 단구간 프레임들의 구간을 상기 타격음 구간의 시작인 것으로 판단하는 단계; 및
    상기 타격음 구간의 시작 이후에 존재하는 단구간 프레임에서 상기 함수(R(n))의 출력 값이 기 설정된 출력값 미만인 것으로 나오는 경우, 해당 단구간 프레임을 상기 타격음 구간의 종점인 것으로 판단하는 단계를 포함하는, 방법.
  19. 하나 이상의 프로세서들;
    메모리; 및
    하나 이상의 프로그램들을 포함하고,
    상기 하나 이상의 프로그램들은 상기 메모리에 저장되고, 상기 하나 이상의 프로세서들에 의해 실행되도록 구성되며,
    상기 하나 이상의 프로그램들은,
    음향 신호를 획득하기 위한 명령;
    상기 획득한 음향 신호를 샘플링하여 복수 개의 단구간 프레임으로 구분하기 위한 명령;
    상기 각 단구간 프레임에 대한 단구간 프레임 에너지를 산출하기 위한 명령;
    상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density)를 산출하기 위한 명령;
    상기 각 단구간 프레임이 유효음 구간에 포함되는지 여부와 관련한 판별 추정값을 산출하기 위한 명령;
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값에 대해 유효음 구간의 경계를 결정하기 위한 임계값을 각각 설정하기 위한 명령; 및
    상기 각 단구간 프레임에 대해 산출한 단구간 프레임 에너지, 스펙트럼 밀도, 및 판별 추정값과 상기 임계값을 기반으로 상기 획득한 음향 신호에서 유효음 구간을 결정하기 위한 명령을 포함하고,
    상기 스펙트럼 밀도를 산출하기 위한 명령은,
    하기의 수학식을 통해 상기 각 단구간 프레임의 스펙트럼 밀도(Spectrum Density : SD를 산출하는, 컴퓨팅 장치.
    Figure 112017092210556-pat00061

    M : MFCC(Mel Frequency Cepstral Coefficient)의 총 계수의 개수
    Xi : 해당 단구간 프레임의 i번째 MFCC 계수의 값
    Ni : 기 설정된 잡음 신호의 i번째 MFCC 계수의 값
KR1020160156741A 2016-11-23 2016-11-23 음향 분석 방법 및 이를 이용한 비파괴 검사 장치 KR101847638B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020160156741A KR101847638B1 (ko) 2016-11-23 2016-11-23 음향 분석 방법 및 이를 이용한 비파괴 검사 장치
PCT/KR2017/013422 WO2018097620A1 (ko) 2016-11-23 2017-11-23 이상 음원 감지 방법 및 이를 수행하기 위한 장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020160156741A KR101847638B1 (ko) 2016-11-23 2016-11-23 음향 분석 방법 및 이를 이용한 비파괴 검사 장치

Publications (1)

Publication Number Publication Date
KR101847638B1 true KR101847638B1 (ko) 2018-04-10

Family

ID=61974954

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020160156741A KR101847638B1 (ko) 2016-11-23 2016-11-23 음향 분석 방법 및 이를 이용한 비파괴 검사 장치

Country Status (2)

Country Link
KR (1) KR101847638B1 (ko)
WO (1) WO2018097620A1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20200125362A (ko) 2019-04-25 2020-11-04 울산대학교 산학협력단 음향신호를 이용한 비파괴 검사장치 및 검사방법
CN115964670A (zh) * 2022-12-17 2023-04-14 中国人民解放军32802部队 一种频谱异常检测方法

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111755025B (zh) * 2019-03-26 2024-02-23 苏州君林智能科技有限公司 一种基于音频特征的状态检测方法、装置及设备
CN112151055B (zh) * 2020-09-25 2024-04-30 北京猿力未来科技有限公司 音频处理方法及装置
CN113205815B (zh) * 2021-04-28 2023-02-28 维沃移动通信有限公司 语音处理方法和电子设备
CN114812796B (zh) * 2022-04-22 2024-05-31 中国人民解放军63811部队 设备状况评估方法以及相关设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007315863A (ja) * 2006-05-24 2007-12-06 Jtekt Corp アコースティックエミッション検出装置および制御装置
JP2009008823A (ja) * 2007-06-27 2009-01-15 Fujitsu Ltd 音響認識装置、音響認識方法、及び、音響認識プログラム
KR101657257B1 (ko) * 2015-06-25 2016-09-13 (주)세아스 절삭공구 용접 불량 검사 장치

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101180687B1 (ko) * 2010-11-12 2012-09-07 대진대학교 산학협력단 음향감지를 이용한 품질관리 시스템 및 그 제공방법
KR101637027B1 (ko) * 2014-08-28 2016-07-08 주식회사 아이티매직 사운드 신호에서 진단 신호를 추출하는 방법 및 진단 장치

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007315863A (ja) * 2006-05-24 2007-12-06 Jtekt Corp アコースティックエミッション検出装置および制御装置
JP2009008823A (ja) * 2007-06-27 2009-01-15 Fujitsu Ltd 音響認識装置、音響認識方法、及び、音響認識プログラム
KR101657257B1 (ko) * 2015-06-25 2016-09-13 (주)세아스 절삭공구 용접 불량 검사 장치

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20200125362A (ko) 2019-04-25 2020-11-04 울산대학교 산학협력단 음향신호를 이용한 비파괴 검사장치 및 검사방법
CN115964670A (zh) * 2022-12-17 2023-04-14 中国人民解放军32802部队 一种频谱异常检测方法

Also Published As

Publication number Publication date
WO2018097620A1 (ko) 2018-05-31

Similar Documents

Publication Publication Date Title
KR101847638B1 (ko) 음향 분석 방법 및 이를 이용한 비파괴 검사 장치
US9536523B2 (en) Method and system for identification of speech segments
CN106941008B (zh) 一种基于静音段的异源音频拼接篡改盲检测方法
US7957967B2 (en) Acoustic signal classification system
Skowronski et al. Acoustic detection and classification of microchiroptera using machine learning: lessons learned from automatic speech recognition
US20050192795A1 (en) Identification of the presence of speech in digital audio data
WO2011135001A1 (en) Assessing speech prosody
Venter et al. Automatic detection of African elephant (Loxodonta africana) infrasonic vocalisations from recordings
US9058821B2 (en) Computer-readable medium for recording audio signal processing estimating a selected frequency by comparison of voice and noise frame levels
Khanagha et al. Phonetic segmentation of speech signal using local singularity analysis
Zeng et al. Classifying watermelon ripeness by analysing acoustic signals using mobile devices
WO2010071428A1 (en) An apparatus, a method and a computer program for recognition of flow regimes in a multiphase fluid flowing in a conduit
Arsikere et al. Automatic estimation of the first three subglottal resonances from adults’ speech signals with application to speaker height estimation
KR101808810B1 (ko) 음성/무음성 구간 검출 방법 및 장치
CN106356076A (zh) 基于人工智能的语音活动性检测方法和装置
Singh et al. Linear Prediction Residual based Short-term Cepstral Features for Replay Attacks Detection.
Sadjadi et al. Robust front-end processing for speaker identification over extremely degraded communication channels
Ma et al. A percussion method with attention mechanism and feature aggregation for detecting internal cavities in timber
Xu et al. IPSO-VMD based signal feature extraction and internal defect detection of hardwood logs through acoustic impact test
US7292981B2 (en) Signal variation feature based confidence measure
Pal et al. Modified energy based method for word endpoints detection of continuous speech signal in real world environment
Hübschen et al. Bitrate and tandem detection for the amr-wb codec with application to network testing
Sharma et al. Speech recognition of Punjabi numerals using synergic HMM and DTW approach
US10482897B2 (en) Biological sound analyzing apparatus, biological sound analyzing method, computer program, and recording medium
JP2016080767A (ja) 周波数成分抽出装置、周波数成分抽出方法及び周波数成分抽出プログラム

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant