KR101065243B1 - 주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치 - Google Patents

주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치 Download PDF

Info

Publication number
KR101065243B1
KR101065243B1 KR1020100058740A KR20100058740A KR101065243B1 KR 101065243 B1 KR101065243 B1 KR 101065243B1 KR 1020100058740 A KR1020100058740 A KR 1020100058740A KR 20100058740 A KR20100058740 A KR 20100058740A KR 101065243 B1 KR101065243 B1 KR 101065243B1
Authority
KR
South Korea
Prior art keywords
frequency band
sound source
band
source level
instantaneous noise
Prior art date
Application number
KR1020100058740A
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 KR1020100058740A priority Critical patent/KR101065243B1/ko
Application granted granted Critical
Publication of KR101065243B1 publication Critical patent/KR101065243B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H15/00—Measuring mechanical or acoustic impedance
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

본 발명은 주파수 대역별 순간 소음에 대한 음원 준위를 산출하는 방법에 있어서, 센서로부터 측정된 순간 소음 신호에 대해 저주파 필터를 이용하여 저주파 대역 필터링을 수행하는 단계; 비주기적인 임펄스성 신호의 실효값을 산출하기 위한 지수 가중치 함수를 설계하는 단계; 상기 저주파 대역 필터링을 수행한 신호와 상기 지수 가중치 함수를 이용하여 컨볼루션(convolution)을 수행하는 단계; 상기 컨볼루션을 수행한 후, 얻어진 데이터에 대해 다중 대역 필터링를 수행하여 다수의 대역별 데이터를 얻는 단계; 다수의 각 대역별 소음 신호에 대해 실효값을 산출하는 단계; 상기 산출된 실효값으로부터 각 대역별 첨두값을 구하는 단계; 및 상기 각 대역별 첨두값에 매개변수를 이용하여 보상하는 단계를 포함하여 이루어지는 것을 특징으로 한다.

Description

주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치{A METHOD AND APPARATUS FOR CALCULATION SOURCE LEVEL OF TRANSIENT NOISE THROUGH UNDERWATER AND STRUCTURE-BORNE BY FREQUENCY BAND}
본 발명은 선박, 자동차 및 기계류에서 발생할 수 있는 수중음향 및 선체 구조전달 순간 소음의 진단 및 평가를 위해 주파수 대역별 순간 소음 음원 준위 산출 방법에 관한 것이다.
함정의 특수성능 시험평가 중 수중방사소음 시험평가는 크게 3가지로 개별장비 평가, 항해 중 종합 평가, 순간 소음 평가로 분류된다. 이 중 순간 소음 시험평가는 최근 수중 함정이 건조되면서 처음으로 시험평가 수행이 요구되었다. 그리고 함정이나 선박 건조시 순간 소음에 대한 진단 목적으로 선체 구조전달 순간 소음을 사전에 평가함으로써, 함 건조과정에서 미리 순간 소음을 줄일 수 있는 대책이 요구되고 있다. 함정에서의 순간 소음이란 어뢰 발사관, 기동 방향을 조정하는 타기(수평타, 수직타), 견시 및 공격 목적으로 잠망경 등이 작동할 때 짧은 시간 동안(수초에서 수십초 이내) 발생하는 소음이다.
순간 소음의 음원 준위 평가는 일반적으로 적용되는 스펙트럼 분석 기법(정상상태 신호 분석)이 아니라, 비주기적이며 무작위로 발생하는 임펄스성 소음을 분석하기 위해 지수 가중치 함수를 두어 실효값을 구하는 방법이 효과적이다.
본 발명에서는 선박뿐만 아니라 일반 기계류, 자동차 등에서 발생할 수 있는 순간 소음의 주파수 대역별 음원 준위 산출 방법을 제안하여 해결하고자 함에 있다.
본 발명은 주파수 대역별 순간 소음에 대한 음원 준위를 산출하는 방법에 있어서, 센서로부터 측정된 순간 소음 신호에 대해 저주파 필터를 이용하여 저주파 대역 필터링을 수행하는 단계; 비주기적인 임펄스성 신호의 실효값을 산출하기 위한 지수 가중치 함수를 설계하는 단계; 상기 저주파 대역 필터링을 수행한 신호와 상기 지수 가중치 함수를 이용하여 컨볼루션(convolution)을 수행하는 단계; 상기 컨볼루션을 수행한 후, 얻어진 데이터에 대해 다중 대역 필터링를 수행하여 다수의 대역별 데이터를 얻는 단계; 다수의 각 대역별 소음 신호에 대해 실효값을 산출하는 단계; 상기 산출된 실효값으로부터 각 대역별 첨두값을 구하는 단계; 및 상기 각 대역별 첨두값에 매개변수를 이용하여 보상하는 단계를 포함하여 이루어진다.
또한, 가속도 센서를 활용한 선체 구조전달 순간소음 음원준위는 상기 가속도 센서로부터 측정된 가속도 신호를 속도 신호로 변환하기 위해 상기 다수의 대역별 첨두값에 대해 적분을 수행하는 단계를 더 포함하여 이루어지는 것을 특징으로 한다.
또한, 상기 수중음향 신호의 각 대역별 첨두값은 1uPa을 기준값(Reference Value)으로 하여 상기 주파수 대역별 순간 소음에 대한 음원 준위를 산출하는 것을 특징으로 한다.
또한, 상기 선체 구조전달 신호의 각 대역별 첨두값에 대한 적분값은 5E-8 m/s을 기준값으로 하여 음원 준위를 산출하는 것을 특징으로 한다.
또한, 상기 매개변수 보상은 상기 주파수 대역별 순간 소음에 대한 음원 준위 산출을 위해 증폭기 이득, 센서 수신감도, 전달손실 등의 값들을 정의한 후, 상기 주파수 대역별 첨두값에 적용하는 것을 특징으로 한다.
또한, 주파수 대역별 순간 소음에 대한 음원 준위를 산출하기 위한 장치의 유저 인터페이스(User Interface)에 있어서, 주파수 대역별 순간 소음에 대한 음원 준위 산출을 위해 필요한 파라미터값들을 설정하기 위한 파라미터값 설정부; 다중 대역 필터링을 수행하기 위해, 상기 다중 대역 필터링에 사용되는 필터를 선택하기 위한 필터 선택부; 상기 주파수 대역별 순간 소음에 대한 음원 준위의 기준값과 산출값에 대한 정보를 표시하기 위한 정보 표시부; 및 순간 소음에 대한 음원 준위 산출, 자료 저장, 수치 초기화, 프로그램을 종료하기 위해 구비되는 버튼부를 포함하여 이루어지는 것을 특징으로 한다.
또한, 상기 정보 표시부는 시간 표시부, 그래프 표시부, 소음 표시부, 데이터 표시부를 포함하여 구성되는 것을 특징으로 한다.
본 발명은 선박뿐만 아니라 일반 기계류, 자동차 등에서 발생할 수 있는 순간 소음의 주파수 대역별 음원 준위 산출 방법을 통해, 함정이나 선박 건조시 순간 소음에 대한 진단 목적으로 선체 구조전달 순간 소음을 사전에 평가함으로써, 함 건조과정에서 미리 순간 소음을 줄일 수 있는 효과가 있다.
도 1은 주파수 대역별 수중 순간 소음 음원 준위 산출 방법에 관한 흐름도.
도 2는 주파수 대역별 선체 구조전달 순간 소음 음원 준위 산출 방법에 관한 흐름도.
도 3은 주파수 대역별 순간 소음 음원 준위 산출을 위한 장치의 사용자 인터페이스(User Interface)를 나타낸 도
도 4는 주파수 대역별 순간 소음 음원 준위 산출을 위한 장치의 실제 사용자 인터페이스(User Interface)를 나타낸 도.
이하, 본 발명의 바람직한 실시 예를 첨부된 도면의 참조와 함께 상세히 설명한다. 그리고, 본 발명을 설명함에 있어서, 관련된 공지기능 혹은 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단될 경우 그 상세한 설명은 생략할 것이다. 그리고 후술되는 용어들은 본 발명에서의 기능을 고려하여 정의된 용어들로서 이는 사용자, 운용자의 의도 또는 관례 등에 따라 달라질 수 있다. 그러므로 그 정의는 본 명세서 전반에 걸친 내용을 토대로 내려져야 할 것이다.
이하, 선박, 자동차 및 기계류에서 발생할 수 있는 수중음향 및 선체 구조전달 순간 소음의 진단 및 평가를 위해 주파수 대역별 순간 소음 음원 준위 산출 방법에 관해 설명한다.
먼저, 도 1은 주파수 대역별 수중 순간 소음 음원 준위 산출 방법에 관한 흐름도를 나타낸다.
상기 도 1을 참조하면, 주파수 대역별 순간 소음에 대한 음원 준위를 산출하기 위해, 수중 음향 센서로부터 측정된 순간 소음 신호에 대해 저주파 대역 통과 필터(Lowpass Filter)를 이용하여 저주파 대역 필터링을 수행한다.(S100) 상기 수중 음향 센서로부터 획득한 소음 신호(P)에 대해 분석하고자 하는 주파수 대역(0~XkHz)까지 저주파 필터를 적용하여 소음 신호(PL)를 얻는다. 소음 신호(PL)는 하기 [수학식 1]과 같다.
Figure 112010039792690-pat00001
다음으로, 비주기적인 임펄스성 신호의 실효값을 산출하기 위한 지수 가중치 함수(Exponential Weighting Function)를 산출한다.(S110)
비주기적인 임펄스성 신호의 실효값을 얻기 위해, 하기 [수학식 2]와 같이 지수 가중치 함수(A)를 설계하여 적용한다.
Figure 112010039792690-pat00002
여기서 T는 샘플링 시간(sampling time)이며, τ는 적분시간이다. B'는 지수 가중치 함수이며, C는 측정된 신호의 에너지 평가시 B'에 의해 원 신호의 에너지보다 초과 추정되는 것을 방지하기 위해 적분시간과 샘플링 시간과의 관계로 정의된 별도의 가중치이다. 지수 가중치 함수 A는 원 신호의 에너지보다 초과 추정이 되지 않기 위해 A(n)의 합이 1이 되도록 N의 값을 결정한다.
다음으로, 저주파 대역 필터링을 수행한 신호와 상기 지수 가중치 함수를 이용하여 컨볼루션(convolution)을 수행한다.(S120)
저주파 필터를 통과한 PL과 지수 가중치 함수 A를 가지고 하기 [수학식 3]과 같이 컨볼루션을 수행한다.
Figure 112010039792690-pat00003
다음으로, 상기 컨볼루션을 수행한 후, 얻어진 데이터에 대해 다중 대역 필터링(Multiple Bandpass Filtering)을 수행하여 다수의 대역별 데이터를 얻는다.(S130)
컨볼루션을 취해 얻어진 데이터(PC)에 대해 다중 대역통과 필터링을 수행하여 각 대역별 데이터(PCB)를 획득한다. 이는 M개의 B Hz 주파수 대역폭을 설정하여 수행하며, 각 시험평가 조건에 따라 달리 설정한다. M개의 대역에 대한 디지털 필터 설계는 실제 신호처리 샘플링 주파수 및 필터 종류에 따라 달라지며, 필터 종류는 Butterworth, Chebyschef II 등 여러 가지 필터를 적용할 수 있다.
다음으로, 다수의 각 대역별 소음 신호에 대해 실효값을 산출한다.(S140)
대역통과 필터를 거친 M개 대역의 소음 신호(PCB)에 대해 각 대역별 실효값을 취하기 위해 하기 [수학식 4]와 같이 수행한다.
Figure 112010039792690-pat00004
다음으로, 상기 산출된 실효값으로부터 각 대역별 첨두값을 구한다.(S150)
첨두값은 하기 [수학식 5]와 같이, M개의 대역별 실효값으로부터 각 대역별 첨두값을 획득하여 구한다.
Figure 112010039792690-pat00005
다음으로, 상기 각 대역별 첨두값에 매개변수(Parameter Compensation)를 이용하여 보상을 수행한다.(S160) 상기 매개변수 보상은 상기 주파수 대역별 순간 소음에 대한 음원 준위 산출을 위해 증폭기 이득, 센서 수신감도, 전달손실 등의 값들을 정의한 후, 상기 주파수 대역별 첨두값에 적용한다.
다음으로, 상기 매개변수 보상을 수행한 값을 이용하여, 수중에서의 주파수 대역별 순간 소음 음원 준위(Source Level)를 산출한다.(S170)
상기 주파수 대역별 순간 소음 음원 준위는 하기 [수학식 6]와 같이,주파수 대역별 첨두값에 매개변수 보상을 하여 최종적으로 얻어진다. 즉, 각 대역별 첨두값에 1uPa을 기준값(Reference Value)으로 하여 상기 주파수 대역별 순간 소음에 대한 음원 준위를 산출하게 된다. 여기서, 산출된 음원 준위의 단위는 dB(decibel)로 표현된다.
Figure 112010039792690-pat00006
Figure 112010039792690-pat00007
Figure 112010039792690-pat00008

상기 [수학식 6]에서 PSL은 산출된 각 대역별 순간소음 음원준위이고, Lw는 기준값 1uPa에 대해 상기 [수학식 5]에서 얻어진 첨두값의 준위이며, Lsens는 수중 음향센서의 수신감도로서 센서의 제작 성능에 따라 결정된다. Lamp는 증폭기의 이득이며, Ltl은 음파전달 손실로서 수심(D)과 전달거리(R)에 의해 결정된다.
도 2는 주파수 대역별 선체 구조전달 순간 소음 음원 준위 산출 방법에 관한 흐름도를 나타낸다.
상기 도 2를 참조하면, 주파수 대역별 순간 소음에 대한 음원 준위를 산출하기 위해, 가속도 센서로부터 측정된 순간 소음 신호에 대해 저주파 대역 통과 필터(Lowpass Filter)를 이용하여 저주파 대역 필터링을 수행한다.(S200) 상기 가속도 센서로부터 획득한 소음 신호(P)에 대해 분석하고자 하는 주파수 대역(0~XkHz)까지 저주파 필터를 적용하여 상기 [수학식 1]과 같이 소음 신호(PL)를 얻는다.
Figure 112010039792690-pat00009
다음으로, 비주기적인 임펄스성 신호의 실효값을 산출하기 위한 지수 가중치 함수(Exponential Weighting Function)를 산출한다.(S210)
비주기적인 임펄스성 신호의 실효값을 얻기 위해, 상기 [수학식 2]와 같이 지수 가중치 함수(A)를 설계하여 적용한다.
Figure 112010039792690-pat00010
여기서 T는 샘플링 시간(sampling time)이며, τ는 적분시간이다. B'는 지수 가중치 함수이며, C는 측정된 신호의 에너지 평가시 B'에 의해 원 신호의 에너지보다 초과 추정되는 것을 방지하기 위해 적분시간과 샘플링 시간과의 관계로 정의된 별도의 가중치이다. 지수 가중치 함수 A는 원 신호의 에너지보다 초과 추정이 되지 않기 위해 A(n)의 합이 1이 되도록 N의 값을 결정한다.
다음으로, 저주파 대역 필터링을 수행한 신호와 상기 지수 가중치 함수를 이용하여 컨볼루션(convolution)을 수행한다.(S220)
저주파 필터를 통과한 PL과 지수 가중치 함수 A를 가지고 상기 [수학식 3]과 같이 컨볼루션을 수행한다.
Figure 112010039792690-pat00011
다음으로, 상기 컨볼루션을 수행한 후, 얻어진 데이터에 대해 다중 대역 필터링(Multiple Bandpass Filtering)를 수행하여 다수의 대역별 데이터를 얻는다.(S230)
컨볼루션을 취해 얻어진 데이터(PC)에 대해 다중 대역통과 필터링를 수행하여 각 대역별 데이터(PCB)를 획득한다. 이는 M개의 B Hz 주파수 대역폭을 설정하여 수행하며, 각 시험평가 조건에 따라 달리 설정한다. M개의 대역에 대한 디지털 필터 설계는 실제 신호처리 샘플링 주파수 및 필터 종류에 따라 달라지며, 필터 종류는 Butterworth, Chebyschef II 등 여러 가지 필터를 적용할 수 있다.
다음으로, 다수의 각 대역별 소음 신호에 대해 실효값을 산출한다.(S240)
대역통과 필터를 거친 M개 대역의 소음 신호(PCB)에 대해 각 대역별 실효값을 취하기 위해 상기 [수학식 4]와 같이 수행한다.
Figure 112010039792690-pat00012

다음으로, 상기 산출된 실효값으로부터 각 대역별 첨두값을 구한다.(S250)
첨두값은 상기 [수학식 5]와 같이, M개의 대역별 실효값으로부터 각 대역별 첨두값을 획득하여 구한다.
Figure 112010039792690-pat00013
다음으로, 상기 가속도 센서로부터 측정된 가속도 신호를 속도 신호로 변환하기 위해 상기 다수의 대역별 첨두값에 대해 적분을 수행한다.(S260) 여기서, 적분을 수행하는 과정은 선체 구조 전달로부터의 순간 소음에 대한 음원 준위를 산출하는 방법에만 적용된다. 상기 [수학식 6]을 통해 적분을 수행한다.
Figure 112010039792690-pat00014
다음으로, 상기 각 대역별 첨두값에 대한 적분값에 매개변수(Parameter Compensation)를 이용하여 보상을 수행한다.(S270)
상기 매개변수 보상은 상기 주파수 대역별 순간 소음에 대한 음원 준위 산출을 위해 증폭기 이득, 센서 수신감도 등의 값들을 정의한 후, 상기 주파수 대역별 첨두값에 적용한다.
다음으로, 상기 매개변수 보상을 수행한 값을 이용하여, 선체 구조 전달에서의 주파수 대역별 순간 소음 음원 준위(Source Level)를 산출한다.(S280)
상기 주파수 대역별 순간 소음 음원 준위는 상기 [수학식 7]와 같이,주파수 대역별 첨두값에 매개변수 보상을 하여 최종적으로 얻어진다. 상기 각 대역별 첨두값에 대한 적분값은 5E-8 m/s을 기준값으로 하여 음원 준위를 산출하게 된다. 여기서, 산출된 음원 준위의 단위는 dB(decibel)로 표현된다.
Figure 112010039792690-pat00015
Figure 112010039792690-pat00016
상기 [수학식 7]에서 PSL은 산출된 각 대역별 순간소음 음원준위이고, Lws는 기준값 5E-8에 대해 상기 [수학식 6]에서 얻어진 적분값의 준위이며, Lsens는 가속도센서의 수신감도로서 센서의 제작 성능에 따라 결정된다. Lamp는 증폭기의 이득이다.
도 3은 주파수 대역별 수중 및 선체 구조전달 순간 소음 음원 준위 산출을 위한 장치의 사용자 화면 구성도를 나타낸다.
상기 도 3을 참조하면, 주파수 대역별 순간 소음에 대한 음원 준위를 산출하기 위한 장치의 사용자 인터페이스(User Interface)(300)는 파라미터값 설정부(310), 필터 선택부(320), 정보 표시부(330), 버튼부(340)로 구성되어 있다.
먼저, 파라미터값 설정부(310)는 주파수 대역별 순간 소음에 대한 음원 준위 산출을 위해 필요한 파라미터값들을 정의한다.
필터 선택부(320)는 음원 준위 산출 과정 중, 다중 대역 필터링을 수행하기 위해 다중 대역 통과 필터링시 적용할 필터를 선택한다.
정보 표시부(330)는 그래프 표시부, 시간 표시부, 데이터 표시부,분석 테이블, 소음 표시부 등으로 구성되어 있다. 이는 음원 준위 산출 결과물들에 대한 정보를 사용자 화면에 표시해준다.
그래부 표시부(Ref. Data File)는 순간 소음에 대한 음원 준위의 기준값 및 산출값을 사용자 화면에 그래프로 보여주게 된다.
시간 표시부(Marking File)는 음원 준위 산출 장치의 작동 시간 구간을 화면에 표시한다.
상기 시간 표시부에서는 상기 음원 준위 산출 장치의 시작시간과 종료시간 정보를 인식하여 상기 소음 표시부에 인식된 시간정보를 연동하여 분석 구간을 설정한다.
데이터 표시부는 문제 소음의 발생 시점을 알기 위해, 주파수 대역별로 시간 영역에 대한 최종 처리 결과를 화면에 보여준다.
분석 테이블(Analysis Table)은 주파수 대역별 순간 소음에 대한 음원 준위의 기준값과 산출값을 수치로 보여준다.
소음 표시부(Raw Data File)는 음원 준위 산출 장치가 작동하여 발생한 실제 소음을 표시 및 청음이 가능하도록 해준다. 상기 소음 표시부는 음원 준위 산출 장치에서 발생한 소음이 아닌, 외부 소음을 제거할 수 있는 기능이 있다.
버튼부(340)는 순간 소음에 대한 음원 준위 산출, 자료 저장, 수치 초기화, 프로그램 종료 등을 수행하기 위해 각 기능에 해당하는 다수의 버튼으로 구성되어 있다.
도 4a, 4b는 주파수 대역별 순간 소음 음원 준위 산출을 위한 장치의 실제 사용자 인터페이스(User Interface)를 나타낸 도면이다.
상기 도 4a/b에서도 도시된 것처럼, 실제 소음을 화면상에 나타내고 소음에 대한 청취도 가능하다. 또한, 주파수 대역별 음원 준위 기준값과 산출값을 수치로 확인할 수 있으며, 이를 그래프를 통해서 확인할 수도 있다.
300: 사용자 인터페이스(User Interface)
310: 파라미터값 설정부
320: 필터 선택부
330: 정보 표시부
340: 버튼부

Claims (9)

  1. 주파수 대역별 순간 소음에 대한 음원 준위를 산출하는 방법에 있어서,
    센서로부터 측정된 순간 소음 신호에 대해 저주파 필터를 이용하여 저주파 대역 필터링을 수행하는 단계;
    비주기적인 임펄스성 신호의 실효값을 산출하기 위해 지수 가중치 함수를 설계하는 단계;
    상기 저주파 대역 필터링을 수행한 신호와 상기 지수 가중치 함수를 이용하여 컨볼루션(convolution)을 수행하는 단계;
    상기 컨볼루션을 수행한 후, 얻어진 데이터에 대해 다중 대역 필터링를 수행하여 다수의 대역별 데이터를 얻는 단계;
    다수의 각 대역별 소음 신호에 대해 실효값을 산출하는 단계;
    상기 산출된 실효값으로부터 각 대역별 첨두값을 구하는 단계; 및
    상기 각 대역별 첨두값에 매개변수를 이용하여 보상하는 단계를 포함하여 이루어지는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 방법.
  2. 제 1항에 있어서,
    상기 센서는 가속도 센서이며, 상기 가속도 센서로부터 측정된 가속도 신호를 속도 신호로 변환하기 위해 상기 다수의 대역별 첨두값에 대해 적분을 수행하는 단계를 더 포함하여 이루어지는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 방법.
  3. 제2항에 있어서,
    상기 적분을 수행하는 단계는 하기 [수학식 6]을 이용하여 수행하는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 방법.
    [수학식 6]
    Figure 112010039792690-pat00017

  4. 제 1항에 있어서,
    상기 각 대역별 첨두값은 1uPa을 기준값(Reference Value)으로 하여 상기 주파수 대역별 순간 소음에 대한 음원 준위를 산출하는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 방법.
  5. 제 2항에 있어서,
    상기 각 대역별 첨두값에 대한 적분값은 5E-8 m/s을 기준값으로 하여 음원 준위를 산출하는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 방법.
  6. 제 1항에 있어서,
    상기 매개변수 보상은 상기 주파수 대역별 순간 소음에 대한 음원 준위 산출을 위해 증폭기 이득, 센서 수신감도, 전달손실 등의 값들을 정의한 후, 상기 주파수 대역별 첨두값에 적용하는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 방법.
  7. 제 1항에 있어서,
    상기 가중치 함수는 하기 [수학식 2]를 이용하여 산출하는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 방법.
    [수학식 2]
    Figure 112010039792690-pat00018


    여기서 는 샘플링 시간(sampling time)이며, 는 적분시간이다. 은 지수 가중치 함수이며, 는 측정된 신호의 에너지 평가시 에 의해 원 신호의 에너지보다 초과 추정되는 것을 방지하기 위해 적분시간과 샘플링 시간과의 관계로 정의된 별도의 가중치이다.
  8. 주파수 대역별 순간 소음에 대한 음원 준위를 산출하기 위한 장치의 유저 인터페이스(User Interface)에 있어서,
    주파수 대역별 순간 소음에 대한 음원 준위 산출을 위해 필요한 파라미터값들을 설정하기 위한 파라미터값 설정부;
    다중 대역 필터링을 수행하기 위해, 상기 다중 대역 필터링에 사용되는 필터를 선택하기 위한 필터 선택부;
    상기 주파수 대역별 순간 소음에 대한 음원 준위의 기준값과 산출값에 대한 정보를 표시하기 위한 정보 표시부; 및
    순간 소음에 대한 음원 준위 산출, 자료 저장, 수치 초기화, 프로그램을 종료하기 위해 구비되는 버튼부를 포함하여 이루어지는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 장치.
  9. 제 8항에 있어서,
    상기 정보 표시부는 시간 표시부, 그래프 표시부, 소음 표시부, 데이터 표시부를 포함하여 구성되는 것을 특징으로 하는 주파수 대역별 순간 소음에 대한 음원 준위 산출 장치.
KR1020100058740A 2010-06-21 2010-06-21 주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치 KR101065243B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020100058740A KR101065243B1 (ko) 2010-06-21 2010-06-21 주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100058740A KR101065243B1 (ko) 2010-06-21 2010-06-21 주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치

Publications (1)

Publication Number Publication Date
KR101065243B1 true KR101065243B1 (ko) 2011-09-19

Family

ID=44957448

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100058740A KR101065243B1 (ko) 2010-06-21 2010-06-21 주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치

Country Status (1)

Country Link
KR (1) KR101065243B1 (ko)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001074616A (ja) 1999-09-06 2001-03-23 Mitsubishi Electric Corp 回転機の異常診断装置
KR100658302B1 (ko) 2005-12-15 2006-12-14 넥센타이어 주식회사 순간 접지면의 주파수 파형을 고려한 타이어의소음저감방법

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001074616A (ja) 1999-09-06 2001-03-23 Mitsubishi Electric Corp 回転機の異常診断装置
KR100658302B1 (ko) 2005-12-15 2006-12-14 넥센타이어 주식회사 순간 접지면의 주파수 파형을 고려한 타이어의소음저감방법

Similar Documents

Publication Publication Date Title
CN105723458B (zh) 用于滤除噪声和还原声信号遭到衰减的高频分量的动态阈值方法、系统、计算机可读介质
US10026418B2 (en) Abnormal frame detection method and apparatus
US20110075860A1 (en) Sound source separation and display method, and system thereof
US20140311245A1 (en) Pipe inspection system and related methods
EP2326933B1 (en) Improvements in and relating to apparatus for the airborne acoustic inspection of pipes
KR101430321B1 (ko) 오디오 시스템의 지각 품질을 결정하기 위한 방법 및 시스템
EP2554124A1 (en) Ultrasound observation device, operation method of ultrasound observation device, and operation program of ultrasound observation device
Avanesians et al. Wave separation: Application for arrival time detection in ultrasonic signals
KR101648597B1 (ko) 파형의 정밀 측정
CN111736222A (zh) 单炮数据信噪比确定方法及装置
RU2620569C1 (ru) Способ измерения разборчивости речи
KR101065243B1 (ko) 주파수 대역별 수중 및 선체 구조전달 순간소음 음원준위 산출 방법 및 장치
CN110132405A (zh) 一种摘挂挡的噪音识别系统及方法
Vergara et al. Measurement of cement porosity by centroid frequency profiles of ultrasonic grain noise
JP6485958B2 (ja) 異音解析装置
RU2736188C9 (ru) Способ отображения гидроакустической информации
JP6926695B2 (ja) 振幅変調成分の分析方法及び変調信号分析システム
EP2958343B1 (en) Apparatus for testing directionality in hearing instruments
JP6541179B2 (ja) 信号処理装置
CN115424636B (zh) 一种基于语音传输指数的音质测量方法及装置
WO2017179342A1 (ja) 信号処理装置、レーダ装置、および、信号処理方法
JP6347660B2 (ja) 魚群探知機、単体魚検出方法、および単体魚検出プログラム
JP6201722B2 (ja) マルチパス評価装置、及びマルチパス評価方法
EP4362496A1 (en) System and method for switching a frequency response and directivity of microphone
JP4198690B2 (ja) ノイズ測定方法、ノイズ測定プログラム、およびノイズ測定装置

Legal Events

Date Code Title Description
A201 Request for examination
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20140901

Year of fee payment: 4

LAPS Lapse due to unpaid annual fee