KR101323592B1 - 대역 분할 처리를 통한 데몬 처리 성능 향상 방법 - Google Patents

대역 분할 처리를 통한 데몬 처리 성능 향상 방법 Download PDF

Info

Publication number
KR101323592B1
KR101323592B1 KR1020130088729A KR20130088729A KR101323592B1 KR 101323592 B1 KR101323592 B1 KR 101323592B1 KR 1020130088729 A KR1020130088729 A KR 1020130088729A KR 20130088729 A KR20130088729 A KR 20130088729A KR 101323592 B1 KR101323592 B1 KR 101323592B1
Authority
KR
South Korea
Prior art keywords
signal
band
noise
spectrum
band signal
Prior art date
Application number
KR1020130088729A
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 KR1020130088729A priority Critical patent/KR101323592B1/ko
Application granted granted Critical
Publication of KR101323592B1 publication Critical patent/KR101323592B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523—Details of pulse systems
    • G01S7/526—Receivers
    • G01S7/527—Extracting wanted echo signals
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/74—Systems using reradiation of acoustic waves, e.g. IFF, i.e. identification of friend or foe
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

본 발명은 데몬 처리에 관한 것으로서, 더 상세하게는 함정의 방사 소음을 이용하여 함정을 식별하는 수동 소나에 적용되어 적은 계산량으로 정확한 데몬 처리가 가능한 데몬 처리 성능 향상 방법에 대한 것이다.
본 발명의 데몬 처리 기법에 따르면, 분할된 대역이 서로 독립적이므로 일반적인 기법에 비해 적은 계산량으로 필요한 성능을 나타내는 것이 가능하다.

Description

대역 분할 처리를 통한 데몬 처리 성능 향상 방법{Method for enhancing DEMON processing performance using multiband processing}
본 발명은 데몬 처리에 관한 것으로서, 더 상세하게는 함정의 방사 소음을 이용하여 함정을 식별하는 수동 소나에 적용되어 적은 계산량으로 정확한 데몬 처리가 가능한 데몬 처리 성능 향상 방법에 대한 것이다.
수중에서 원거리까지 전달 가능한 음파는 함정의 탐지 및 식별하기 위한 효과적인 수단으로 각종 소나에 활용되고 있다. 그 중 수동 소나의 경우 표적 식별을 위해서 방사 소음 신호의 분석이 필수적이며, 특히 함정의 프로펠러 축 회전 수, 날개 회전 수, 축 수 등의 정보를 확인하기 위해서 데몬(DEMON: Demodulation Envelop Modulated On Noise) 처리를 통한 신호 분석이 수행된다.
기존의 데몬 처리는 진폭 변조된 프로펠러의 소음 신호를 복조하여 변조 신호를 추출하는 신호 처리 기법으로 그 성능의 향상을 위해 많은 연구가 진행되고 있다.
이들 연구 기법들은 신호대 잡음비(signal to noise ratio)를 높이기 위해서 시간 영역에 대해 음향 신호를 분할하여 처리하는 방법을 사용하였다. 즉 일정 시간동안 음향 신호를 수집하여 처리 구간(section)을 구분한다. 그 다음 각각의 구간에 대한 분할 처리 후 평균을 취해 잡음의 분산(variance)을 감소시켜 신호대 잡음비를 향상시켰다. 또한 일반적으로 처리 구간은 음향 신호 수집 시간을 단축하기 위해 시간 영역에서 중첩(overlap)하여 사용하였다.
그러나 이러한 중첩 기법은 중첩 비율이 너무 높으면 각 처리 구간 음향 신호가 서로 독립적이지 않아서 신호대 잡음비 성능 향상이 미미하고, 반대로 중첩 비율이 너무 낮으면 음향 신호 수집에 많은 시간이 소요되는 단점이 있다.
또한, 기존의 데몬 신호 처리 기법은 신호 중첩을 이용하므로 신호대 잡음비 향상을 가져올 수 있지만 처리 음향 신호들의 중첩으로 인해 상관관계를 가져서 성능향상이 제한적이었다. 또한, 필요로 하는 성능 향상을 위해서는 많은 계산량이 소요 되었다.
1. 한국등록특허번호 제10-1104200호 2. 한국공개특허번호 제10-2012-0104025호
1. 황수복, 김진석, 이철묵 "토널 신호 간섭에 강인한 데몬 처리 기법", 한국음향학회지 31권 6호, pp 384-390, 2012년.
본 발명은 위 배경기술에 따른 문제점을 해소하기 위해 제안된 것으로서, 적은 계산량으로 필요한 성능을 나타낼 수 있는 대역 분할 처리를 통한 데몬 처리 성능 향상 방법을 제공하는데 그 목적이 있다.
또한, 본 발명은 함정을 식별하는 수동 소나에 적용되어 적은 계산량으로 정확한 데몬 처리가 가능한 대역 분할 처리를 통한 데몬 처리 성능 향상 방법을 제공하는데 다른 목적이 있다.
본 발명은 위에서 제시된 과제를 달성하기 위해, 적은 계산량으로 필요한 성능을 나타낼 수 있는 대역 분할 처리를 통한 데몬 처리 성능 향상 방법을 제공한다.
상기 데몬 처리 성능 향상 방법은,
프로펠러 소음 신호가 입력되는 신호 입력 단계;
다수의 대역 통과 필터부를 이용하여 분할 대역 신호를 생성하는 분할 대역 신호 생성 단계;
생성된 분할 대역 신호에 대하여 제곱을 수행하는 제곱 수행 단계;
다수의 저역 통과 필터부를 이용하여 제곱 수행된 분할 대역 신호로부터 고주파 신호를 제거하는 고주파 신호 제거 단계;
고주파 신호가 제거된 분할 대역 신호로부터 DC(Direct Current) 신호를 제거하는 DC 신호 제거 단계;
DC 신호가 제거된 분할 대역 신호를 FFT(Fast Fourier Transform) 변환하는 FFT 변환 단계;
FFT 변환된 분할 대역 신호를 배경 소음이 제거되도록 소음 규준화를 수행하여 스펙트럼을 생성하는 스펙트럼 생성 단계; 및
생성된 스펙트럼을 평균화하여 평균 스펙트럼을 생성하는 평균 스펙트럼 생성 단계;를 포함하는 것을 특징으로 할 수 있다.
이때, 상기 다수의 대역 통과 필터부는 주파수 영역에서 중첩되지 않는 것을 특징으로 할 수 있다.
이때, 상기 분할 대역 신호는 서로 독립적인 것을 특징으로 할 수 있다.
또한, 상기 분할 대역 신호는, 수학식
Figure 112013067861613-pat00001
(여기서, 연산자
Figure 112013067861613-pat00002
는 합성곱(convolution)을 의미하고,
Figure 112013067861613-pat00003
,
Figure 112013067861613-pat00004
는 주파수 영역에서 서로 중첩되지 않는 대역 통과 필터의 임펄스 응답(impulse response)을 나타내며,
Figure 112013067861613-pat00005
을 나타낸다)으로 표시되는 것으로 특징으로 할 수 있다.
또한, 상기 평균 스펙트럼은, 수학식
Figure 112013067861613-pat00006
(여기서,
Figure 112013067861613-pat00007
는 소음 규준화된 이후의 스펙트럼을 나타내며,
Figure 112013067861613-pat00008
,
Figure 112013067861613-pat00009
이다)에 의해 산출되는 것을 특징으로 할 수 있다.
다른 한편으로 본 발명의 다른 실시예는, 프로펠러 소음 신호가 입력됨에 따라 분할 대역 신호를 생성하는 대역 통과 필터 블럭; 생성된 분할 대역 신호에 대하여 제곱을 수행하는 제곱 수행 블럭; 제곱 수행된 분할 대역 신호로부터 고주파 신호를 제거하는 저역 통과 필터 블럭; 고주파 신호가 제거된 분할 대역 신호로부터 DC(Direct Current) 신호를 제거하는 노치 필터 블럭; DC 신호가 제거된 분할 대역 신호를 FFT(Fast Fourier Transform) 변환하는 FFT 변환 블럭; FFT 변환된 분할 대역 신호를 배경 소음이 제거되도록 소음 규준화를 수행하여 스펙트럼을 생성하는 소음 규준화 블럭; 및 생성된 스펙트럼을 평균화하여 평균 스펙트럼을 생성하는 평균화부;를 포함하는 것을 특징으로 하는 대역 분할 처리를 통한 데몬 처리 성능 향상 장치를 제공한다.
또한, 상기 대역 통과 필터 블럭, 제곱 수행 블럭, 저역 통과 필터 블럭, 노치 필터 블럭, FFT 변환 블럭, 및 소음 규준화 블럭은 각각 병렬로 연결되는, 다수의 대역 통과 필터부, 다수의 제곱부, 다수의 저역 통과 필터부, 다수의 노치 필터부, 다수의 FFT 변환부, 및 다수의 소음 규준화부를 갖는 것을 특징으로 할 수 있다.
일반적인 데몬 신호 처리 기법은 신호 중첩을 이용하여 신호 대 잡음비 향상을 가져올 수 있지만 처리 음향 신호들의 중첩으로 인해 상관관계를 가져서 성능향상이 제한적이었다. 또한 필요로 한 성능 향상을 위해서는 많은 계산량이 소요되었다.
그러나, 본 발명의 데몬 처리 기법에 따르면, 분할된 대역이 서로 독립적이므로 일반적인 기법에 비해 적은 계산량으로 필요한 성능을 나타내는 것이 가능하다.
또한, 본 발명의 다른 효과로서는 함정의 방사 소음을 이용하여 함정을 식별하는 수동 소나에 적용되어 적은 계산량으로 정확한 데몬 처리가 가능하게 할 것으로 기대된다는 점을 들 수 있다.
도 1은 일반적인 데몬 처리 기법의 블록도이다.
도 2는 도 1에 중첩 기법을 적용하는 개념도이다.
도 3은 본 발명의 일실시예에 따른 대역 분할에 의한 데몬 성능 향상 장치의 구성도를 나타낸 것이다.
도 4는 본 발명의 일실시예에 따른 대역 분할에 의한 데몬 성능 향상 과정을 보여주는 흐름도이다.
본 발명은 다양한 변경을 가할 수 있고 여러 가지 실시예를 가질 수 있는바, 특정 실시예들을 도면에 예시하고 상세한 설명에 구체적으로 설명하고자 한다. 그러나 이는 본 발명을 특정한 실시 형태에 대해 한정하려는 것이 아니며, 본 발명의 사상 및 기술 범위에 포함되는 모든 변경, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다.
각 도면을 설명하면서 유사한 참조부호를 유사한 구성요소에 대해 사용한다.
제 1, 제 2등의 용어는 다양한 구성요소들을 설명하는데 사용될 수 있지만, 상기 구성요소들은 상기 용어들에 의해 한정되어서는 안 된다. 상기 용어들은 하나의 구성요소를 다른 구성요소로부터 구별하는 목적으로만 사용된다.
예를 들어, 본 발명의 권리 범위를 벗어나지 않으면서 제 1 구성요소는 제 2 구성요소로 명명될 수 있고, 유사하게 제 2 구성요소도 제 1 구성요소로 명명될 수 있다. "및/또는" 이라는 용어는 복수의 관련된 기재된 항목들의 조합 또는 복수의 관련된 기재된 항목들 중의 어느 항목을 포함한다.
다르게 정의되지 않는 한, 기술적이거나 과학적인 용어를 포함해서 여기서 사용되는 모든 용어들은 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에 의해 일반적으로 이해되는 것과 동일한 의미가 있다.
일반적으로 사용되는 사전에 정의되어 있는 것과 같은 용어들은 관련 기술의 문맥상 가지는 의미와 일치하는 의미를 가지는 것으로 해석되어야 하며, 본 출원에서 명백하게 정의하지 않는 한, 이상적이거나 과도하게 형식적인 의미로 해석되지 않아야 한다.
이하 첨부된 도면을 참조하여 본 발명의 일실시예에 따른 대역 분할 처리를 통한 데몬 처리 성능 향상 장치 및 방법을 상세하게 설명하기로 한다.
일반적으로 프로펠러(미도시)에 의한 소음 신호
Figure 112013067861613-pat00010
은 프로펠러(미도시)의 회전으로 발생하는 하모닉 변조 신호와 광대역 캐비테이션 소음 신호와의 곱으로서 다음식과 같이 표현된다.
Figure 112013067861613-pat00011
여기서,
Figure 112013067861613-pat00012
은 수신된 음향 신호의 개수이고,
Figure 112013067861613-pat00013
은 변조 지수,
Figure 112013067861613-pat00014
은 변조 주파수,
Figure 112013067861613-pat00015
는 샘플링(sampling) 주파수이다.
Figure 112013067861613-pat00016
은 광대역 캐비테이션 소음 신호이고
Figure 112013067861613-pat00017
은 광대역 소음이며
Figure 112013067861613-pat00018
과
Figure 112013067861613-pat00019
은 서로 독립(independent)인 특성을 가진다.
일반적으로 변조 지수
Figure 112013067861613-pat00020
는 0.1 ~ 0.5의 값을 가지며, 변조 주파수
Figure 112013067861613-pat00021
은 10 Hz ~ 100Hz의 값을 가진다. 그리고 광대역 캐비테이션 소음 신호
Figure 112013067861613-pat00022
과 광대역 소음
Figure 112013067861613-pat00023
의 주파수 대역 폭(bandwidth)은 수kHz 이상이며, 따라서 샘플링 주파수
Figure 112013067861613-pat00024
도 수 kHz이상의 값을 가진다.
도 1은 일반적인 데몬 처리 기법의 블록도이다. 도 1을 참조하면, 데몬 프로세싱 시스템(100)에서, 프로펠러 소음 신호
Figure 112013067861613-pat00025
은 대역 통과 필터부(bandpass filter)(110)를 적용하여 캐비테이션 신호가 존재하는 주파수 대역을 선택한다. 그 후 선택된 주파수 대역에서 변조 신호의 복조를 위해 제곱부(120)에 의해 제곱(square)이 수행된다.
일반적으로 변조 신호는 저주파 대역에 존재하므로 저주파 대역 통과 필터부(lowpass filter)(130) 및 노치 필터부(notch filter)(140)를 이용하여 고주파 신호 및 DC 신호를 제거한다.
다음으로 시간 영역(time domain) 신호를 주파수 영역(frequency domain)으로 변환하기 위해 FFT(Fast Fourier Transform)부(150)를 이용하여 FFT를 수행하고 마지막으로 배경 소음을 제거하기 위해 소음 규준화부(noise normalization)(160)를 이용하여 소음 규준화를 수행한 후 스펙트럼(spectrum)
Figure 112013067861613-pat00026
을 생성한다.
일반적으로 데몬 처리 기법들은 스펙트럼의 신호대 잡음비를 높이기 위해 Welch 기법, 즉 중첩 기법을 사용한다. 즉 다음식과 같이
Figure 112013067861613-pat00027
을
Figure 112013067861613-pat00028
개의 구간으로 나누어 처리하여 분산을 감소시켜 신호대 잡음비를 향상시킨다. 이러한 중첩 기법을 도 1의 블럭도에 적용하면 도 2와 같이 나타낼 수 있다.
Figure 112013067861613-pat00029
여기서,
Figure 112013067861613-pat00030
은 구간의 크기, 즉 FFT(Fast Fourier Transform) 크기이며,
Figure 112013067861613-pat00031
는 중첩된 구간의 시작 위치이다. 만약
Figure 112013067861613-pat00032
이면 구간들이 서로 중첩되고,
Figure 112013067861613-pat00033
이면 구간들은 중첩되지 않는다. 그리고 중첩 비율
Figure 112013067861613-pat00034
은
Figure 112013067861613-pat00035
이다.
도 2에서 "데몬 처리부(100a 내지 100K)"는 도 1의 "데몬 처리 블럭(100) "과 동일하며, 평균화부(210)는 각 데몬 처리부(100a 내지 100K)의 출력 신호
Figure 112013067861613-pat00036
에 대해 평균을 취하여 스펙트럼
Figure 112013067861613-pat00037
을 생성한다. 이것은 다음식과 같이 나타낼 수 있다.
Figure 112013067861613-pat00038
이러한 시간 영역에 대한 중첩 기법에 의한 출력 스펙트럼
Figure 112013067861613-pat00039
의 분산은 개별 스펙트럼
Figure 112013067861613-pat00040
,
Figure 112013067861613-pat00041
의 분산과 다음식과 같은 관계를 가진다.
Figure 112013067861613-pat00042
여기서,
Figure 112013067861613-pat00043
은
Figure 112013067861613-pat00044
의 분산이며, 등호는
Figure 112013067861613-pat00045
,
Figure 112013067861613-pat00046
가 서로 독립(independent)일 때 성립하고, 부등호는
Figure 112013067861613-pat00047
가 서로 상관관계(correlation)가 있을 때를 나타낸다.
그러므로 시간 영역에 대한 중첩 기법은 개별 스펙트럼
Figure 112013067861613-pat00048
,
Figure 112013067861613-pat00049
가 서로 독립이 아니므로, 스펙트럼
Figure 112013067861613-pat00050
의 분산은
Figure 112013067861613-pat00051
의 평균에 의한 감소 효과가 제한적이다.
하모닉 변조 신호와 광대역 캐비테이션 소음 신호와의 곱을 표현한 식에서 보는 바와 같이 프로펠러(미도시)의 의한 소음 신호는 날개의 회전으로 발생하는 변조 신호를 광대역 캐비테이션 소음 신호가 캐리어(carrier)의 역할을 수행하여 발생한다.
이는 무수히 많은 캐리어가 존재하므로 그 중 일부 캐리어 성분만으로도 변조 신호를 복조할 수 있음을 의미한다. 다시 말해 광대역 캐비테이션 소음 신호에서 일부 주파수 대역 신호를 이용하여 날개의 회전으로 발생하는 변조 신호를 복조할 수 있다. 따라서 본 발명의 일실시예에서는 이러한 성질을 이용하여 대역 분할에 의한 데몬 성능 향상 기법을 제안한다.
이를 자세히 살펴보면, 수학식 1에서 나타낸 프로펠러에 의한 소음 신호
Figure 112013067861613-pat00052
은 다음식과 같이
Figure 112013067861613-pat00053
개의 스펙트럼으로 분할(spectrally- disjoint)된 대역 신호
Figure 112013067861613-pat00054
으로 나타낼 수 있다.
Figure 112013067861613-pat00055
여기서, 연산자 ＊는 합성곱(convolution)을 의미하고,
Figure 112013067861613-pat00056
,
Figure 112013067861613-pat00057
는 주파수 영역에서 서로 중첩되지 않는 대역 통과 필터의 임펄스 응답(impulse response)이다. 그리고
Figure 112013067861613-pat00058
은 다음식과 같이 정의된다.
Figure 112013067861613-pat00059
따라서
Figure 112013067861613-pat00060
을 시간 영역에 대해 중첩되지 않은
Figure 112013067861613-pat00061
개의 대역으로 나눈
Figure 112013067861613-pat00062
,
Figure 112013067861613-pat00063
,
Figure 112013067861613-pat00064
에 대해 데몬 처리를 수행하여 평균을 취하면
Figure 112013067861613-pat00065
는 서로 독립이므로 분산을 감소시켜 신호대 잡음비를 향상시킬 수 있다.
도 3은 본 발명의 일실시예에 따른 대역 분할에 의한 데몬 성능 향상 장치의 구성도를 나타낸 것이다. 도 3을 참조하면, 데몬 성능 향상 장치는, 프로펠러 소음 신호가 입력됨에 따라 분할 대역 신호를 생성하는 대역 통과 필터 블럭(310)과, 생성된 분할 대역 신호에 대하여 제곱을 수행하는 제곱 수행 블럭(320)과, 제곱 수행된 분할 대역 신호로부터 고주파 신호를 제거하는 저역 통과 필터 블럭(330)과, 고주파 신호가 제거된 분할 대역 신호로부터 DC(Direct Current) 신호를 제거하는 노치 필터 블럭(340)과, DC 신호가 제거된 분할 대역 신호를 FFT(Fast Fourier Transform) 변환하는 FFT 변환 블럭(350)과, FFT 변환된 분할 대역 신호를 배경 소음이 제거되도록 소음 규준화를 수행하여 스펙트럼을 생성하는 소음 규준화 블럭(360)과, 생성된 스펙트럼을 평균화하여 평균 스펙트럼을 생성하는 평균화부(370) 등을 포함하여 구성된다.
여기서, 상기 대역 통과 필터 블럭(310)은 다수의 대역 통과 필터부(311)가 병렬로 구성되며, 제곱 수행 블럭(320)은 다수의 제곱부(321)가 병렬로 구성되며, 저역 통과 필터 블럭(330)은 다수의 저역 통과 필터부(331)가 병렬로 구성되며, 노치 필터 블럭(340)은 다수의 노치 필터부(341)가 병렬로 구성되며, FFT 변환 블럭(350)은 다수의 FFT 변환부(351)가 병렬로 구성되며, 소음 규준화 블럭(360)은 다수의 소음 규준화부(361)가 병렬로 구성된다.
따라서, 프로펠러 소음 신호
Figure 112013067861613-pat00066
,
Figure 112013067861613-pat00067
을 주파수 영역에서 중첩되지 않는
Figure 112013067861613-pat00068
개의 대역 통과 필터부(310)를 적용하여 수학식 5에서 나타낸 분할 대역 신호
Figure 112013067861613-pat00069
을 생성한다.
그 후 선택된 주파수 대역에서 변조 신호의 복조를 위해
Figure 112013067861613-pat00070
개의 제곱부(320)를 이용하여 제곱을 수행하고 저주파 대역 통과 필터부(330)와 노치 필터부(340)를 이용하여 고주파 신호 및 DC 신호를 제거한다.
다음으로 시간 영역 신호를 주파수 영역으로 변환하기 위해 FFT(Fast Fourier Transform)부(350)를 이용하여 FFT를 수행하고 소음 규준화부(360)를 통해 배경 소음을 제거하기 위해 소음 규준화를 수행한 후 스펙트럼
Figure 112013067861613-pat00071
,
Figure 112013067861613-pat00072
,
Figure 112013067861613-pat00073
을 생성한다. 최종적으로 분산을 감소시켜 신호대 잡음비를 높이기 위해 다음식과 같이 평균화부(310)를 이용하여
Figure 112013067861613-pat00074
의 평균을 취하여 스펙트럼
Figure 112013067861613-pat00075
을 생성한다.
Figure 112013067861613-pat00076
도 4는 본 발명의 일실시예에 따른 대역 분할에 의한 데몬 성능 향상 과정을 보여주는 흐름도이다. 도 4를 참조하면, 프로펠러 소음 신호가 입력되면, 다수의 대역 통과 필터부를 이용하여 분할 대역 신호가 생성된다(단계 S400,S410).
분할 대역 신호가 생성되면 생성된 분할 대역 신호에 대하여 제곱을 수행한다(단계 S420). 부연하면, 선택된 주파수 대역에서 변조 신호의 복조를 위해 제곱을 수행하게 된다.
제곱 수행이 이루어지면, 다수의 저역 통과 필터부를 이용하여 제곱 수행된 분할 대역 신호로부터 고주파 신호를 제거하고, 또한 고주파 신호가 제거된 분할 대역 신호로부터 DC(Direct Current) 신호를 제거한다(단계 S430).
DC 신호 및/또는 고주파 신호가 제거되면 이러한 신호가 제거된 분할 대역 신호를 FFT(Fast Fourier Transform) 변환한다(단계 S440).
FFT 변환된 분할 대역 신호를 배경 소음이 제거되도록 소음 규준화를 수행하여 스펙트럼을 생성하고, 생성된 스펙트럼을 평균화하여 평균 스펙트럼을 생성한다(단계 S450,S460).
100: 데몬 프로세싱 시스템
110: 대역 통과 필터부
120: 제곱부
130: 저역 통과 필터부
140: 노치 필터부
150: FFT(Fast Fourier Transform)
160: 소음 규준화부
210: 평균화부
310: 대역 통과 필터 블럭 311: 대역 통과 필터부
320: 제곱 수행 블럭 321: 제곱부
330; 저역 통과 필터 블럭 331: 저역 통과 필터부
340: 노치 필터 블럭 341: 노치 필터부
350: FFT 변환 블럭 351: FFT 변환부
360: 소음 규준화 블럭 361: 소음 규준화부
370: 평균화부

Claims (5)

  1. 프로펠러 소음 신호가 입력되는 신호 입력 단계;
    다수의 대역 통과 필터부를 이용하여 분할 대역 신호를 생성하는 분할 대역 신호 생성 단계;
    생성된 분할 대역 신호에 대하여 제곱을 수행하는 제곱 수행 단계;
    다수의 저역 통과 필터부를 이용하여 제곱 수행된 분할 대역 신호로부터 고주파 신호를 제거하는 고주파 신호 제거 단계;
    고주파 신호가 제거된 분할 대역 신호로부터 DC(Direct Current) 신호를 제거하는 DC 신호 제거 단계;
    DC 신호가 제거된 분할 대역 신호를 FFT(Fast Fourier Transform) 변환하는 FFT 변환 단계;
    FFT 변환된 분할 대역 신호를 배경 소음이 제거되도록 소음 규준화를 수행하여 스펙트럼을 생성하는 스펙트럼 생성 단계; 및
    생성된 스펙트럼을 평균화하여 평균 스펙트럼을 생성하는 평균 스펙트럼 생성 단계;
    을 포함하는 것을 특징으로 하는 대역 분할 처리를 통한 데몬 처리 성능 향상 방법.
  2. 제 1 항에 있어서,
    상기 다수의 대역 통과 필터부는 주파수 영역에서 중첩되지 않는 것을 특징으로 하는 대역 분할 처리를 통한 데몬 처리 성능 향상 방법
  3. 제 1 항에 있어서,
    상기 분할 대역 신호는 서로 독립적인 것을 특징으로 하는 대역 분할 처리를 통한 데몬 처리 성능 향상 방법.
  4. 제 1 항에 있어서,
    상기 분할 대역 신호는, 수학식
    Figure 112013067861613-pat00077
    (여기서, 연산자
    Figure 112013067861613-pat00078
    는 합성곱(convolution)을 의미하고,
    Figure 112013067861613-pat00079
    ,
    Figure 112013067861613-pat00080
    는 주파수 영역에서 서로 중첩되지 않는 대역 통과 필터의 임펄스 응답(impulse response)을 나타내며,
    Figure 112013067861613-pat00081
    을 나타낸다)으로 표시되는 것으로 특징으로 하는 대역 분할 처리를 통한 데몬 처리 성능 향상 방법.
  5. 제 1 항에 있어서,
    상기 평균 스펙트럼은, 수학식
    Figure 112013067861613-pat00082
    (여기서,
    Figure 112013067861613-pat00083
    는 소음 규준화된 이후의 스펙트럼을 나타내며,
    Figure 112013067861613-pat00084
    ,
    Figure 112013067861613-pat00085
    이다)에 의해 산출되는 것을 특징으로 하는 대역 분할 처리를 통한 데몬 처리 성능 향상 방법.
KR1020130088729A 2013-07-26 2013-07-26 대역 분할 처리를 통한 데몬 처리 성능 향상 방법 KR101323592B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020130088729A KR101323592B1 (ko) 2013-07-26 2013-07-26 대역 분할 처리를 통한 데몬 처리 성능 향상 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020130088729A KR101323592B1 (ko) 2013-07-26 2013-07-26 대역 분할 처리를 통한 데몬 처리 성능 향상 방법

Publications (1)

Publication Number Publication Date
KR101323592B1 true KR101323592B1 (ko) 2013-10-30

Family

ID=49639596

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020130088729A KR101323592B1 (ko) 2013-07-26 2013-07-26 대역 분할 처리를 통한 데몬 처리 성능 향상 방법

Country Status (1)

Country Link
KR (1) KR101323592B1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106483520A (zh) * 2016-09-27 2017-03-08 哈尔滨工程大学 一种船舶辐射噪声调制系数估计方法
CN116400337A (zh) * 2023-06-08 2023-07-07 中国人民解放军国防科技大学 基于线段检测的舰船噪声调制线谱提取与轴频估计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08320364A (ja) * 1995-05-26 1996-12-03 Nec Corp パッシブソーナーにおける変調周波数検出装置
US20060007784A1 (en) 2003-09-29 2006-01-12 Donald Lerro Sonar system and process
KR20120065408A (ko) * 2009-09-29 2012-06-20 아틀라스 엘렉트로닉 게엠베하 진폭 변조 광대역 소음을 분석하는 방법 및 장치

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08320364A (ja) * 1995-05-26 1996-12-03 Nec Corp パッシブソーナーにおける変調周波数検出装置
US20060007784A1 (en) 2003-09-29 2006-01-12 Donald Lerro Sonar system and process
KR20120065408A (ko) * 2009-09-29 2012-06-20 아틀라스 엘렉트로닉 게엠베하 진폭 변조 광대역 소음을 분석하는 방법 및 장치

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106483520A (zh) * 2016-09-27 2017-03-08 哈尔滨工程大学 一种船舶辐射噪声调制系数估计方法
CN106483520B (zh) * 2016-09-27 2018-12-07 哈尔滨工程大学 一种船舶辐射噪声调制系数估计方法
CN116400337A (zh) * 2023-06-08 2023-07-07 中国人民解放军国防科技大学 基于线段检测的舰船噪声调制线谱提取与轴频估计方法
CN116400337B (zh) * 2023-06-08 2023-08-18 中国人民解放军国防科技大学 基于线段检测的舰船噪声调制线谱提取与轴频估计方法

Similar Documents

Publication Publication Date Title
Yang et al. Component extraction for non-stationary multi-component signal using parameterized de-chirping and band-pass filter
Millioz et al. Circularity of the STFT and spectral kurtosis for time-frequency segmentation in Gaussian environment
ES2807529T3 (es) Método para la determinación de señal de ruido y aparato del mismo
Wang et al. Application of the dual-tree complex wavelet transform in biomedical signal denoising
Sharma et al. Improved eigenvalue decomposition-based approach for reducing cross-terms in Wigner–Ville distribution
CN102840907A (zh) 早期故障状态下滚动轴承振动信号特征提取和分析方法
CN108458871A (zh) 一种基于改进经验小波变换的齿轮箱故障识别方法
CN111769810B (zh) 一种基于能量峭度谱的流体机械调制频率提取方法
US10061014B2 (en) Radar signal processing method and apparatus for compensating for in-phase/quadrature-phase (I/Q) imbalance
KR101323592B1 (ko) 대역 분할 처리를 통한 데몬 처리 성능 향상 방법
Zhang et al. Improved local cepstrum and its applications for gearbox and rolling bearing fault detection
RU2367970C2 (ru) Устройство обнаружения узкополосных шумовых гидроакустических сигналов на основе вычисления интегрального вейвлет-спектра
CN107991531A (zh) 基于变分模态分解和去除噪声的电网谐波检测方法
KR20170124279A (ko) 수중 방사 소음 측정시 외부표적 소음 제거가 가능한 신호처리 방법 및 장치
CN106840670B (zh) 一种基于多尺度短时光滑分析的周期瞬态信号检测方法
KR20220109953A (ko) 셉스트럼 기반 유성 기어박스의 고장 감지 장치 및 방법
CN103438983B (zh) 一种信号随机平均谱的数据处理方法
CN109660475B (zh) 一种非合作相位编码水声通信信号自主辨识方法
Liu et al. A new motor fault detection method using multiple window S-method time-frequency analysis
CN110909827A (zh) 一种适用于风机叶片声音信号的降噪方法
Zaitsev Elimination on power line interference from ECG signal using combined bidirectional narrow-band notch filter
Boashash et al. Time-frequency analysis applied to signaturing of underwater acoustic signals
CN111308426B (zh) 一种适用于单天线接收机的低信噪比周期调频信号检测与分离方法
Schasse et al. Source localization based on the doppler effect
Huang et al. Short-time kurtogram for bearing fault feature extraction under time-varying speed conditions

Legal Events

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

Payment date: 20161005

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20170904

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20181002

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20191002

Year of fee payment: 7