The Journal of the Acoustical Society of Korea. 30 September 2016. 403-409
https://doi.org/10.7776/ASK.2016.35.5.403

ABSTRACT


MAIN

  • I. 서 론

  • II. 음향 벡터 센서 신호모델

  • III. 음향 벡터 센서 배열의 뒤틀림 보상기법

  •   3.1 음향 벡터 센서를 이용한 지연합 빔형성

  •   3.2 음향 벡터 센서 배열의 뒤틀림 보상

  • IV. 모의실험

  • V. 결 론

I. 서  론

음원의 방향을 추정하기 위해서 일반적으로 센서 배열에 도달하는 음원의 시간지연 차이를 이용한다. 이때 선배열 센서의 경우, 음원의 방향이 방위각과 고각으로 구분되지 않고 원추각의 형태로 결정되어지므로 좌 우 구분 모호성이 발생한다.[1] 이러한 문제를 해결하기 위한 다양한 방법들이 제안되었으며, 운용적인 방법으로 TMA(Target Motion Analysis)가 있다.[2] 하지만 TMA는 긴 관측시간이 필요로 하는 단점이 있기 때문에 단시간 내에 좌 우 구분 모호성을 해결할 수 있는 방법에 대한 연구가 진행되었다. 이러한 연구로써 센서 배열의 구조적 특성을 이용하는 이중배열(twin arrays)[3,4]과 삼중배열(triplet arrays)[5,6] 방법이 제안되었다. 이중배열은 선배열센서 두 개가 나란히 배치된 구조를 가지며, 삼중배열센서는 세 개의 전방향 센서가 근접한 삼각형의 구조를 가지며 이러한 형태가 배열화 되어있는 센서 배열이다. 이때 근접한 센서사이의 거리 때문에 수신된 잡음신호들이 강한 일관성이 발생한다. 이는 삼중배열센서를 이용한 카디오이드 빔형성의 성능 저하에 영향을 미치며, 이를 해결하기 위하여 센서간격을 넓힐 필요가 있다. 이는 센서 배열의 직경이 증가하게 되며, 플랫폼에서 긴 센서 배열을 운용함에 있어 어려움을 발생시킨다. 따라서 음향 벡터 센서를 이용한 단일배열센서를 이용하는 방법이 제안되었다.

벡터 센서[7,8]는 하나의 전방향 센서와 세 개의 지향성 센서로 구성되며, 세 개의 지향성 센서는 서로 수직(http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC98FC.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC997A.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC999B.gif)방향의 속도 성분을 측정하도록 배치되어 있다. 이러한 지향성 센서의 추가적인 정보는 신호 대 잡음비(Signal to Noise Ratio, SNR)의 성능향상뿐만 아니라 음원의 방향성분을 통해 정확한 음원의 방향추정에 유리하다. 정확한 신호의 방향추정은 정확한 축의 정보를 알고 있을 때 가능하다. 하지만 배열화된 음향 벡터 센서들은 운용중에 쉽게 뒤틀림(roll)이 발생할 수 있으며, 이것은 음향 벡터 센서의 축을 틀어지게 하는데 영향을 끼친다. 이러한 축의 뒤틀림은 잘못된 방향의 신호정보를 인가하여 빔형성 기법의 성능저하를 나타낸다. 따라서 본 논문에서는 음향 벡터 센서 축의 뒤틀림에 따른 좌 우 구분의 성능저하를 분석하고 보상을 통한 빔형성 기법의 성능향상을 확인하였다.

I장의 서론에 이어 II장에서는 음향 벡터 센서의 신호모델을 설명하고 III장에서는 센서 배열의 뒤틀림에 따른 성능분석과 보상방법을 제안하였다. IV장에서 시뮬레이션을 통하여 보상에 대한 성능향상을 확인하였으며 V장에서 결론을 맺는다.

II. 음향 벡터 센서 신호모델

Fig. 1은 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC99BB.gif개의 음향 벡터 센서가 등간격으로 배치되었을 때 원거리 음원으로부터 신호를 수신하는 상황을 나타낸 것이다. 음원의 위치는 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC99DB.gif축 기준으로 반시계방향으로 방위 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9A78.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9AA8.gif축 기준으로 고각 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9B36.gif에 위치하며, 배열의 중심에서부터 측정된다. 선 배열 센서의 특성상 방위각과 고각의 성분을 분리할 수 없기 때문에 표적의 방위는 원추각(http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9B66.gif)으로 나타난다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9B96.gif번째 센서에 수신되는 신호 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9C33.gif은 Eq.(1)과 같다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F1.jpg

Fig. 1. Configuration of acoustic vector sensor arraywith source.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9DAB.gif, (1)

여기서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9DDB.gif이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9E1A.gif은 센서 개수이다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9E98.gifhttp://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9ED8.gif이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9F17.gif은 음원의 개수이다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9F47.gif는 전방향 센서신호의 출력이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PIC9FA6.gif는 지향성 센서신호의 출력이다. 표적의 방향벡터 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA005.gif은 Eq.(2)와 같이 나타난다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA092.gifhttp://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA0C2.gif번째 센서에서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA0F2.gif번째 표적신호의 시간지연을 의미한다. 이때 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA122.gifhttp://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA181.gif번째 센서의 위치벡터이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA1C0.gif는 수중에서 음속을 의미한다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA24E.gifhttp://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA26E.gif번째 표적신호이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA3B7.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA454.gif는 전방향 센서와 지향성 센서에 수신되는 백색 가우시안 잡음(white Gaussian noise)이다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA540.gif, (2)

여기서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA5BE.gifhttp://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA5DE.gif번째 표적의 방위각이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA61D.gifhttp://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA63E.gif번째 표적의 고각을 나타낸다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA66E.gif는 전치행렬(transpose)를 의미한다.

III. 음향 벡터 센서 배열의 뒤틀림 보상기법

3.1 음향 벡터 센서를 이용한 지연합 빔형성

음향 벡터 센서를 이용한 지연합 빔형성(delay-sum beamforming) 기법은 Fig. 2와 같이 2단계로 나누어 수행한다. 블록도에서와 같이 첫 번째로 각각의 음향 벡터 센서에서 카디오이드 빔형성을 통해 특정방위의 신호를 얻는다. 두 번째로 첫 번째에서 얻은 신호를 이용하여 지연합 빔형성을 통해 최종 빔출력을 얻는다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F2.jpg

Fig. 2 Block-diagram of acoustic vector sensor arrays beamforming.

Step. 1 각각의 음향 벡터 센서에서 카디오이드 빔형성을 통해 특정방위의 신호를 얻는다.

Step. 2 Step.1에서 얻은 신호를 지연합 빔형성을 통해 최종 빔출력을 얻는다.

각각의 음향 벡터 센서에서 카디오이드 빔형성 기법을 통한 빔출력신호는 Eq.(3)과 같다.[4,5]

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA6DC.gif, (3)

여기서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA70C.gif이며 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA72C.gif은 센서 개수이다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA7AA.gif는 카디오이드 빔형성 기법에서 가중벡터를 의미하며, 관측하는 방향의 방위각이 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA7DA.gif, 고각이 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA7FA.gif인 경우 Eq.(4)와 같이 나타낼 수 있다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA859.gif. (4)

각각의 음향 벡터 센서에서 카디오이드 빔출력신호를 이용하여 지연합빔형성 기법을 적용하여 얻은 최종 빔출력 신호는 Eq.(5)와 같이 나타낼 수 있다.[9]

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA906.gif, (5)

여기서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICA926.gif는 선배열센서에서 표적의 방위각을 의미한다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAA11.gif는 카디오이드 빔출력 신호로 구성된 벡터이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAA41.gif는 선배열 지연합 빔형성에서 가중벡터를 나타내고 Eq.(6)과 같이 표현된다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAAA0.gif, (6)

여기서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAAEF.gif이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAB0F.gif는 센서사이의 거리, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAB30.gif는 조향각이다.

빔형성기의 빔출력 에너지 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICABBD.gif는 다음과 같다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAC6A.gif, (7)

여기서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICACF8.gif은 크기가 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAD18.gif인 센서 출력의 공분산행렬이다.

3.2 음향 벡터 센서 배열의 뒤틀림 보상

음향 벡터 센서에서 정확한 신호의 방향추정은 축의 정보를 정확히 알고 있다는 가정하에 이루어진다. 배열화된 음향 벡터 센서들은 운용중에 뒤틀림이 쉽게 발생할 수 있으며, 이로 인하여 빔형성 기법의 성능저하를 발생시킨다. 따라서 적절한 보상을 통하여 성능복원이 필요하다. 보상은 배열센서 사이에 뒤틀림의 정도를 확인할 수 있는 롤 센서의 정보를 이용한다. 센서 배열의 구조상 튜브안에 센서를 고정시키기 때문에 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAD86.gif축의 변화는 센서 배열의 형태에 주로 영향을 받으며, 뒤틀림은 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICADA7.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICADD7.gif축의 변화가 크게 작용한다. 따라서 Fig. 1에서 음향 벡터 센서 배열의 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICADF7.gif축은 고정되어 있고 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAE17.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAE37.gif축의 변화를 뒤틀림이라 가정한다. http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAE58.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAE78.gif축의 변화에 따른 보상은 음향 벡터 센서의 가중벡터에 롤 센서에서 얻은 뒤틀림정도를 이용하여 이루어지며, 보상된 가중벡터 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICAEF6.gif는 Eq.(8)과 같다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB04F.gif, (8)

여기서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB06F.gif는 롤 센서에서 얻은 뒤틀림 각도정보이며, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB09F.gif는 롤 센서의 오차를 의미한다.

Fig. 3은 센서 배열의 뒤틀림 정도에 따른 빔출력에너지 결과를 나타낸 것이다. 실선은 뒤틀림이 발생하지 않은 경우의 결과이며, 70°에 위치한 표적이 좌 우 구분되어 나타나는 결과를 확인할 수 있다. 대시-닷 선은 뒤틀림이 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB0DE.gif, 대시 선은 뒤틀림이 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB14D.gif, 점선은 뒤틀림이 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB17C.gif만큼 발생했을 때의 결과이다. 뒤틀림은 첫 번째 센서부터 마지막 센서까지 선형적으로 뒤틀림이 발생했다고 가정하였다. 즉, 뒤틀림이 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB1AC.gif만큼 발생하면 센서사이의 뒤틀림 각도 차이는 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB1DC.gif만큼 변화하게 되고 마지막 센서는 첫 번째 센서와 비교하여 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB21C.gif만큼 회전하게 된다. 따라서 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB24C.gif축과 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB2CA.gif축이 바뀐 신호가 수신되게 되며, 뒤틀림이 더 크게 발생할수록 좌 우 구분 성능이 저하되는 것을 알 수 있다. 실제 센서 배열에서 뒤틀림은 선형적으로 발생하지 않는다. 따라서 롤 센서가 음향 벡터 센서마다 위치하게 되어 정확한 뒤틀림 보상이 된다면 가장 효과적이지만 롤 센서의 사용이 제한적이라면 보상효과가 제한적으로 나타날 수 있다. 따라서 본 논문에서는 시뮬레이션을 통해 뒤틀림 보상에 따른 성능분석을 통해 좌 우 구분성능을 확인한다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F3.jpg

Fig. 3. Beam power without roll compensation.

IV. 모의실험

본 장에서는 음향 벡터 센서 배열의 뒤틀림에 따른 빔형성 성능분석과 뒤틀림 보상에 따른 성능 복원에 대한 분석을 수행하였다. 실험에서 음향 벡터 센서의 개수는 32개의 채널을 이용하였으며, 센서간 거리는 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB309.gif이며, 설계주파수는 1000 Hz로 설정하였다. Fig. 1과 같이 ULA(Uniform Line Array)를 가정하였으며, 뒤틀림은 3.2절에서 언급한 바와 같이 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB339.gif, http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB359.gif축에 회전으로 발생한다고 가정한다.

Fig. 4는 Fig. 3의 결과와 동일한 조건에서 뒤틀림 보상을 시행한 빔출력 결과이다. 보상을 위한 롤 센서는 음향 벡터 센서의 개수와 동일하게 32개를 이용하였으며, 음향 벡터 센서와 동일한 위치에 있다고 가정하였다. 또한 롤 센서 오차(http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB379.gif)는 가우시안 분포를 가정하였으며, 표준편차는 http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/PICB39A.gif로 설정하였다. 롤 센서를 음향 벡터 센서와 동일한 위치에서 동일한 개수만큼 사용하였기 때문에 뒤틀림 보상이 잘 수행되며, Fig. 4에서 보여 주는 것과 같이 우수한 보상성능 나타낸다. 그러나 실제 상황에서 보상은 이보다 더 열악한 환경에서 이루어진다. 본 모의실험에서 롤 센서의 측정 오차 및 사용되는 개수에 따른 성능 차이를 분석한다. 그리고 음향 벡터 센서 배열의 뒤틀림이 임의적으로 발생한다고 가정하였을 때 사용되는 롤 센서의 개수에 따른 성능 차이를 분석한다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F4.jpg

Fig. 4. Beam power with roll compensation.

Fig. 5는 음원이 위치한 방위각의 좌 우 대칭점의 빔 출력차이를 롤 센서의 오차에 따라 나타낸 결과이다. 롤 센서의 개수나 위치 조건은 앞의 Fig. 4와 동일하다. x축의 롤각도는 앞의 Fig. 3에서 설명된 센서 배열의 처음과 끝의 롤각도 전체변화량을 의미하며, 이번 실험에서도 센서 배열에서 롤이 선형적으로 변화한다고 가정한다. y축의 좌 우 빔출력 차이 값은 표적이 70°에 위치할 때 -70°의 빔출력 값과의 차이를 의미한다. 삼각형으로 마킹된 실선은 뒤틀림 보상이 없을 때의 결과를 나타내며, 비틀림이 점점 많이 발생할수록 좌 우 구분성능 저하가 크게 나타나는 것을 확인할 수 있다. 롤각도가 270°를 넘어갈 때 좌 우 구분성능이 올라가는 형태는 전체 배열에서 음원방향을 조향하는 벡터 센서의 수가 증가하여 나타나는 현상이다. 즉, 롤각도가 270°이면 전체배열의 1/3은 음원방향을 나머지 2/3은 반대방향을 조향하게되며, 롤각도가 360°이면 다시 전체배열의 1/2는 음원방향을 나머지 1/2는 반대방향을 조향하게 되어서 나타나는 현상이다. 대시-닷선(원형, 별표, 사각형, 다이아몬드형)은 뒤틀림 보상이 이루어진 결과를 나타내며, 보상 시에 롤 센서의 측정 오차에 대한 표준편차 5°, 10°, 20°, 40°일 경우를 나타낸 것이다. 각각의 경우에 대해서 몬테카를로 시뮬레이션 200회의 수행결과를 나타낸 것이다. 롤 센서의 측정 오차가 클수록 좌우 구분 성능이 저하됨을 확인할 수 있으나, 성능 저하의 정도가 크지 않음을 확인할 수 있다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F5.jpg

Fig. 5. Difference of port/starboard power for compen-sation error.

Fig. 6은 센서 배열의 양 끝에 두 개의 롤 센서가 위치해있다고 가정하고 센서 배열의 뒤틀림에 따른 빔출력을 나타낸 것이다. 실험조건은 Fig. 4와 동일하다. 단 센서 배열의 뒤틀림은 0° ~ 270° 만큼 선형적으로 발생시켰다. 따라서 양끝의 롤 센서에는 0°와 270°의 값이 입력되게 된다. 하지만 센서 배열의 변화가 시계방향인지 반시계방향인지는 롤 센서만의 정보로 확인이 불가능하다. 실선은 센서 배열의 뒤틀림 방향과 롤 센서정보를 이용한 뒤틀림 보상방향이 일치하는 경우 적절한 보상으로 인하여 좌 우 구분이 잘 되었으며, 대시-닷선은 뒤틀림 방향과 보상방향이 일치하지 않아서 보상이 이루어지지 않는 결과를 나타내었다. 따라서 잘못된 방향으로 보상이 되었을 때 Fig. 6과 같이 보상의 성능이 심각하게 저하된다. 따라서 적절한 개수의 롤 센서를 센서 배열에 배치해야 보상성능을 만족시킬 수 있다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F6.jpg

Fig. 6. Beam power for roll-direction.

Fig. 7은 실제 음향 벡터 센서 배열의 뒤틀림의 방향을 예측할 수 없기 때문에 음향 벡터 센서 배열의 뒤틀림 변화가 임의적으로 발생하였 때 롤 센서의 개수에 따른 보상성능을 확인한 결과이다. 음향 벡터 센서 배열의 뒤틀림을 다음과 같이 정의하였다. 각각의 음향 벡터 센서 사이의 뒤틀림 변화량을 가우시안분포로 가정하였으며, 표준편차가 10°, 30°, 50°에 대하여 실험하였다. 롤 센서의 배치는 음향 벡터 센서 배열의 전체길이에 등간격으로 배치하는 것을 가정하였다. 나머지 조건은 Fig. 3과 동일하다. y축의 좌 우 빔출력 차이는 앞의 Fig. 5에서 설명한 내용과 동일하게 수행되었다. 실선은 음향 벡터 센서 배열에서 음향 벡터 센서 사이의 뒤틀림 변화량이 표준편차가 10°인 경우에 대한 결과이다. 음향 벡터 센서 사이의 뒤틀림변화가 크지 않아 롤 센서의 개수에 큰 영향을 받지 않았다. 하지만 변화량의 표준편차가 30°, 50°인 경우는 센서 배열의 뒤틀림 변화가 커졌기 때문에 롤 센서의 개수에 따른 성능차이가 크게 나타났다. 실험결과 롤 센서의 개수가 증가함에 따라 보상성능이 증가하였으며, 좌 우 구분성능을 만족하기 위해서 설계사양에 맞는 롤 센서의 사용이 필요할 것으로 판단된다. Fig. 8은 Fig. 7의 모의실험에서 음향 벡터 센서 사이의 뒤틀림 변화량의 표준편차가 30°인 경우에 대해 음원의 방위각에 따른 좌 우 구분성능을 나타낸 것이다. 음원의 방향이 음향 벡터 센서 배열의 0°, 180° end-fire 근처에서는 좌 우 구분성능이 저하되어 보상의 성능여부를 확인하기 어렵다. 90° broadside 근처에서는 롤 센서의 개수에 따라 성능이 Fig. 7과 같이 나타난다.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F7.jpg

Fig. 7. Difference of port/starboard power for a number of roll sensors.

http://static.apub.kr/journalsite/sites/ask/2016-035-05/N0660350510/images/Figure_ASK_35_05_10_F8.jpg

Fig. 8. Difference of port/starboard power for source bearing.

V. 결  론

본 논문에서는 음향 벡터 센서 배열의 뒤틀림에 따른 빔형성 결과에서 좌 우 구분 모호성에 대하여 분석하였다. 그리고 롤 센서의 오차, 센서 배열의 비선형변화에서 센서개수에 따른 좌 우 구분 성능에 대한 영향을 분석하였다. 음향 벡터 센서 배열의 뒤틀림에 대한 적절한 보상을 통한 음향 벡터 센서 배열의 좌 우 구분 성능향상을 확인하였다. 롤 센서의 추정오차에 대한 영향은 좌 우 구분성능에 크게 영향을 미치지 않았으며, 비선형으로 변화하는 센서 배열에서 롤 센서의 개수에 따른 보상성능은 크게 차이가 났다. 따라서 적절한 롤 센서를 이용하여 센서 배열을 구성하는 것이 필요하다. 음향 벡터 센서 배열을 운용함에 있어 뒤틀림에 따른 좌 우 구분 성능저하가 심각하게 나타나기 때문에 적절한 뒤틀림 보상에 대한 고려가 필요함을 확인하였다.

Acknowledgements

이 논문은 2016년도 한화탈레스(주)의 재원을 지원 받아 수행된 연구이다(계약번호 : C-16-001).

References

1
1.R. J. Urick, Principles of Underwater Sound (McGrawHill, New York, 1983), pp. 54-60.
2
2.S. C. Nardone, A. Lindgren, and Kai Gong, “Fundamental properties and performance of conventional bearings-only target motion analysis,” IEEE Transactions on automatic control 29, 775-787 (1984).
3
3.T. Warhonowicz, H. Schmidt-Schierhorn, and H. Hostermann, “Port/Starboard discrimination performance by twin line array for a LFAS sonar system,” in Proc. UDT, Europe, 398 (1999).
4
4.J. P. Feuillet, W. S. Allensworth, and B. K. Newhall, “Nonambigous beamforming for a high resolution twin-line array,” J. Acoust. Soc. Am. 97, 3292 (1995).
5
5.G. W. M. van Mierlo, S. P. Beerens, R. Been, Y. Doisy, and E. Trouve, “Port-Starboard discrimination on hydrophone triplets in active and passive towed arrays,” in Proc. UDT, Pacific, 63-68 (2000).
6
6.Y. W. Ryu, E. C. Jeong, and S. Y. Chun, “Analysis on noise correlation of triplet line array sensors in the east sea” (in Korean), J. Acoust. Soc. Kr. 34, 479-486 (2015).
7
7.A. Nehorai and E. Paldi, “Acoustic vector sensor array processing,” IEEE Trans. signal Process, 42, 2481-2491  (1994).
8
8.J. Cao, “Survey on acoustic vector sensor and its applications in signal processing,” IEEE Control Conference (CCC), Chinese, (2014).
9
9.R. O. Nielsen, Sonar Signal Processing (Norwood, MA: Artech House, 1991), pp. 51-69.
페이지 상단으로 이동하기