KR20100071617A - 무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서 - Google Patents

무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서 Download PDF

Info

Publication number
KR20100071617A
KR20100071617A KR1020080130397A KR20080130397A KR20100071617A KR 20100071617 A KR20100071617 A KR 20100071617A KR 1020080130397 A KR1020080130397 A KR 1020080130397A KR 20080130397 A KR20080130397 A KR 20080130397A KR 20100071617 A KR20100071617 A KR 20100071617A
Authority
KR
South Korea
Prior art keywords
transfer function
sound
impulse response
infinite impulse
response filter
Prior art date
Application number
KR1020080130397A
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 KR1020080130397A priority Critical patent/KR20100071617A/ko
Publication of KR20100071617A publication Critical patent/KR20100071617A/ko

Links

Images

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/0248—Filters characterised by a particular frequency response or filtering method
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S1/00—Two-channel systems
    • H04S1/002—Non-adaptive circuits, e.g. manually adjustable or static, for enhancing the sound image or the spatial distribution
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30—Control circuits for electronic adaptation of the sound field
    • H04S7/301—Automatic calibration of stereophonic sound system, e.g. with test microphone
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04S—STEREOPHONIC SYSTEMS 
    • H04S2420/00—Techniques used stereophonic systems covered by H04S but not provided for in its groups
    • H04S2420/01—Enhancing the perception of the sound image or of the spatial distribution using head related transfer functions [HRTF's] or equivalents thereof, e.g. interaural time difference [ITD] or interaural level difference [ILD]

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Stereophonic System (AREA)

Abstract

본 발명은 측정한 머리전달함수의 임펄스 응답을 사용하여 입체음향을 구현하는 유한 임펄스 응답(FIR) 필터의 계산량이 많은 것을 보완하여 적은 계산량으로 유사한 음향효과를 낼 수 있도록 하기 위하여 무한 임펄스 응답(IIR) 필터 머리전달함수를 구하는 방법과 이를 활용한 입체음향 재생 장치에 관한 것이다. 발명은 모노 음원을 수신하는 제1 채널 및 제2 채널을 구비하고, 상기 제1 및 제2 채널들 각각은 제1 위치 조정 머리관련 전달함수를 구현한 제1 무한 임펄스 응답 필터 및 제1 위치 조정 머리관련 전달함수를 구현한 제1 무한 임펄스 응답 필터를 탑재한 것이며; 상기 무한 임펄스 응답 필터의 전달함수는
Figure 112008087512184-PAT00001
에 의해 구현되는, 입체음향 장치를 제공한다. 적합하게는, 무한 임펄스 응답 필터 계수는 32차일 수 있다.
입체음향 장치, 머리전달함수, 무한 임펄스 응답 필터.

Description

무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서{3D production device using IIR filter-based head-related transfer function, and DSP for use in said device}
본 발명은 측정한 머리전달함수의 임펄스 응답을 사용하여 입체음향을 구현하는 유한 임펄스 응답(FIR) 필터의 계산량이 많은 것을 보완하여 적은 계산량으로 유사한 음향효과를 낼 수 있도록 하기 위하여 무한 임펄스 응답(IIR) 필터 머리전달함수를 활용한 입체음향 장치 및 이에 사용될 수 있는 디지털 신호 프로세서에 관한 것이다.
이 분야에 장치는 화자와 운영자 간의 긴장을 해소하여 보다 자연스러운 방송 운영을 할 수 있도록 화자의 음성을 인식하여 화자의 원고 읽는 속도에 맞추어 자동으로 원고 표출속도를 조절하기 위해 사용될 수 있는 장치이다.
사람은 20 ∼ 20 kHz까지를 들을 수 있으며 두 귀를 가지고 있어 소리의 위치를 인지할 수 있다. 특히, 두 귀가 일정한 거리로 떨어져 있으며, 특별한 귀바퀴 모양을 가지고 있어 음원의 위치가 전후좌우는 물론이고 상하 어디에 있는 인지할 수 있다. 그리고, 이런 위치를 인지하는데 있어 저음보다는 고음에 대하여 민감하 게 인지하는 것으로 알려져 있다.
컴퓨터게임 혹은 가상공간에서는 시각적인 효과뿐만 아니라 음향효과도 가상의 공간에서 다양한 위치와 음장효과가 제공되어야 현실감을 가질 수 있다. 특히, 헤드셋과 같은 음향재생장치를 사용하는 경우에는 가상공간의 특정 객체의 소리가 특정한 위치에서 발생하는 것과 같은 효과를 제공하는 것이 필요하다. 그리고, 화상회의장치에서는 화면에 배치된 사람의 위치에 따라서 소리가 적절한 위치에서 발생하는 것으로 인지할 수 있어야 한다.
이러한 음향효과를 제공하기 위하여 머리전달함수(HRTF, head related transfer function)을 사용한다. 종래의 머리전달함수는 FIR 필터로 구성된다. 머리전달함수는 사람의 머리와 유사한 형태의 더미헤드(dummy head)의 양쪽 귀속에 콘덴서 마이크를 넣어 일정한 간격으로 전후좌우, 상하에 배치된 스피커에서 임펄스를 발생하여 측정한 임펄스 응답을 주로 사용한다. 따라서, 머리전달함수는 임펄스 응답을 계수로 하는 FIR 필터들로 구성된다.
종래의 FIR 필터 머리전달함수는 대개 512차의 계수를 갖는다. 이를 이용하여 공간감이 없는 모노의 음원으로부터 방향감을 가진 입체 음향을 재생하기 위해서는 512차의 계수와 컨벌루션을 하여야 하기 때문에 계산량이 많아 입체음향을 구현하는데 메모리가 많아야 하며 고성능의 프로세서를 사용하여야 한다.
따라서, 본 발명에서는 소요 메모리와 프로세서의 부담을 줄이기 위하여 계산 량이 줄이고자 FIR 필터 머리전달함수를 유사한 성능을 가진 IIR 필터로 구성된 머리전달함수를 구하는 방법에 관한 것이며 IIR 필터 머리전달함수를 이용하여 입체음향을 실현하는 장치에 관한 것이다.
또한, 이러한 류의 장치에 사용될 수 있는 전용의 디지털 신호 프로세서를 제공하는 것이다.
본 발명에 따른 목적은,
모노 음원을 수신하는 제1 채널 및 제2 채널을 구비하고,
상기 제1 및 제2 채널들 각각은 제1 위치 조정 머리관련 전달함수를 구현한 제1 무한 임펄스 응답 필터 및 제1 위치 조정 머리관련 전달함수를 구현한 제1 무한 임펄스 응답 필터를 탑재한 것이며;
상기 무한 임펄스 응답 필터의 전달함수는
Figure 112008087512184-PAT00002
에 의해 구현되는, 입체음향 장치에 의해 달성될 수 있다.
바람직하게, 상기 무한 임펄스 응답 필터 계수는 32차일 수 있다.
본 발명의 또 다른 면에 따라서,
마이크로부터 수신된 음을 디지털화된 신호로 제공하는 아날로그-디지털 변환기;
제1항의 상기 무한 임펄스 응답(IIR) 필터가 탑재되는 디지털 신호 프로세서;
음향 데이터 저장 메모리;
음원위치 IIR 필터 계수 저장 메모리; 및
스피커를 포함하는, 입체음향 장치가 제공된다.
또한, 본 발명의 또 다른 면에서 따라서, 본 발명의 목적은,
모노 음원으로부터 수신된 입력신호에 대해 입체음향을 생성하기 위해 머리관련 전달함수를 무한 임펄스 응답필터로 구현하도록 한 디지털 신호 프로세서를 포함하고,
상기 무한 임펄스 응답 필터의 전달함수는
Figure 112008087512184-PAT00003
에 의해 구현되는, 입체음향 재생 전용 디지털 신호 프로세서에 의해서 달성될 수 있다.
바람직하게, 상기 무한 임펄스 응답 필터 계수는 32차일 수 있다.
본 발명에 의하여 발명된 IIR필터 머리전달함수를 사용함으로써 기존의 FIR 필터 머리전달함수를 사용한 입체음향재생에 비하여 계산량이 작음으로 인해서 소요 메모리가 작아도 되며, 처리를 위한 프로세서의 성능이 좋을 필요가 없기 때문에 소요 비용을 줄일 수 있다.
또한, 하나의 컴퓨터에서 여러 가지 작업을 수행할 때, 입체음향 재생에 대한 부담을 줄일 수 있기 때문에 효과적이다.
본 발명은 종래의 FIR 필터 머리전달함수를 이용한 입체음향 구현에 있어 계산량이 많은 단점을 보완하기 위하여 FIR 필터계수 즉, 임펄스 응답 특성과 유사한 특성을 가지는 IIR 필터 머리전달함수를 구하는 방법이 본 발명에 따라 제공되며, 이를 이용하여 입체음향을 재생하는데 사용될 수 있다.
종래의 FIR 필터 머리전달함수를 이용한 입체음향 재생의 구성은 도 1과 같다. 입체감이 없는 모노(mono) 음원을 특정한 방향에 대한 FIR 필터 머리전달함수(
Figure 112008087512184-PAT00004
,
Figure 112008087512184-PAT00005
)를 거쳐서 각각 오른쪽(R)과 왼쪽 스피커(L)로 출력함으로써 헤드셋을 착용한 청취자는 특정한 방향에서 소리가 발생한 것으로 인지하게 된다. 종래의 FIR 필터 머리전달함수는 오른쪽과 왼쪽 스피커로 출력하기 위하여 1 샘플의 음원신호에 대하여 0각각 512개의 계수에 대하여 컨벌루션을 하기 때문에 512 차의 계수를 저장하기 위한 메모리와 512회의 곱셈과 511회의 덧셈을 기본적으로 수행하여야 한다.
본 발명에서 발명한 IIR 필터 머리전달함수를 이용한 입체음향 재생은 장치의 구성은 종래의 입체음향 재생 장치의 구성과 유사하지만 오른쪽과 왼쪽 스피커로의 출력신호를 생성하는 머리전달함수를 IIR 필터를 사용하는 것이 다르다. IIR 필터의 계수를 32차를 사용하는 경우에는 32의 곱셈과 31회의 덧셈이 소요된다.
계산량은 적으면서 FIR 필터와 그 특성이 유사한 IIR 필터의 계수를 계산하는 과정은 본 발명에 의해 제공되는 이에 대해 구체적으로 도출 방법을 첨부한 도면인 도 2를 사용하여 설명한다. 그 절차는 다음과 같다.
첫 번째 단계는 전달함수를 상태방정식으로 변환하는 과정이다. 특정한 방향에 대한 머리전달함수의 FIR 필터 전달함수는 식 (1)과 같다.
Figure 112008087512184-PAT00006
(1)
여기서, bi(i=0,...,n)은 FIR 필터의 계수이다. 식 (1)의 전달함수를 상태방정식으로 변환하면
Figure 112008087512184-PAT00007
0 (2)
Figure 112008087512184-PAT00008
(3)
이 된다. 여기서, s[n]은 상태변수벡터이며, y[n]은 상태방정식의 출력이다.
식(1)에 대하여 각 상태방정식 (2)와 (3)의 각 행렬과 벡터는 다음과 같다.
Figure 112008087512184-PAT00009
(4)
Figure 112008087512184-PAT00010
(5)
Figure 112008087512184-PAT00011
(6)
Figure 112008087512184-PAT00012
(7)
두 번째 단계는 상태방정식 (2), (3)의 상호 그라미안 행렬(cross-gramian matrix)을 구하는 과정이며, FIR 필터의 경우에는 식(9)과 같은 핸켈 행렬(Hankel matrix)가 된다.
Figure 112008087512184-PAT00013
(8)
Figure 112008087512184-PAT00014
은 식(1)의 FIR 필터 계수이다.
세 번째 단계는 식 (2)와 (3) 에 대한 유사행렬(similarity matrix)을 구한다. 유사행렬은 식 (9)와 식 (10)을 이용하여 식 (11)과 (12)를 구한다.
Figure 112008087512184-PAT00015
(9)
Figure 112008087512184-PAT00016
,
Figure 112008087512184-PAT00017
,
Figure 112008087512184-PAT00018
,
Figure 112008087512184-PAT00019
(10)
Figure 112008087512184-PAT00020
(11)
Figure 112008087512184-PAT00021
(12)
여기서,
Figure 112008087512184-PAT00022
(
Figure 112008087512184-PAT00023
는 대각 행렬)로 두고
Figure 112008087512184-PAT00024
를 대각 행렬로 변환하면
Figure 112008087512184-PAT00025
(13)
Figure 112008087512184-PAT00026
(14)
Figure 112008087512184-PAT00027
(15)
네 번째 단계는 헨켈 행렬을 유사변환된 행렬로 분해한다.
Figure 112008087512184-PAT00028
,
Figure 112008087512184-PAT00029
(16)
여기서,
Figure 112008087512184-PAT00030
는 단위행렬이다.
다섯 번째 단계는 감소시킨 헨켈 행렬을 구한다.
Figure 112008087512184-PAT00031
(17)
Figure 112008087512184-PAT00032
(18)
Figure 112008087512184-PAT00033
,
Figure 112008087512184-PAT00034
는 감소시킨 헨켈 행렬에서 사용한 부분 행렬이고,
Figure 112008087512184-PAT00035
,
Figure 112008087512184-PAT00036
는 무시할 부분 행렬이다.
여섯 번째 단계는 감소시킨 헨켈 행렬을 이용하여 IIR 전달함수를 구한다.
Figure 112008087512184-PAT00037
(19)
IIR 전달함수는 감소시킨 차수 k 에 의하여 결정된다.
도 3은 이러한 본 발명에 따른 IIR필터 전달함수가 구현되어 탑재된 디지털신 호처리기(DSP, digital signal processor)를 이용하여 IIR 필터 입체음향을 구현하기 위한 구성예를 나타낸다.
마이크와 프리앰프 그리고 AD 변환기는 음성 및 음향을 수음하여 디지털신호처리기로 입력된다. 음향데이터 저장 메모리는 디지털화된 음향 데이터를 저장한다. 디지털신호처리기의 음원은 음향데이터 저장 메모리에 저장된다. 음원위치 IIR 필터 계수 저장 메모리는 측정된 입체음향 FIR 필터계수를 식(19)로 계산된 IIR 필터 계수를 저장한다. 디지털신호처리기는 마이크 혹은 음향데이터 저장메모리에 저장된 음원 데이터에 대하여 음원 위치 IIR 필터계수 저장 메모리에 저장된 IIR 필터 계수와 필터 계산 처리를 통하여 입체음향신호를 생성하여 앰프를 거쳐 스피커로 출력되도록 한다.
도 4는 식 (19)로 계산된 IIR 필터 계수를 이용하여 마이크입력신호 혹은 음향 데이터 저장 메모리에 저장된 음향 데이터에 대하여 필터링을 하는 IIR 필터 구조를 예시적으로 도시한 것으로서, 이에 대해 이를테면 직접형 IIR 필터 2 구조(direct form IIR filter structure)가 사용될 수 있다. 그러나 당업자는 발명이 이것으로 제한되지 않고 다양한 많은 디지털 신호 프로세서에 구현될 수 있음을 알 것이다. 이를 위해서, 식 (19)의 전달함수는 다음 식으로 표현할 수 있고 이에 따라 당업자는 본 발명의 IIR필터 기반의 머리전달함수를 전용 프로세서에 구현할 수 있을 것이다.
Figure 112008087512184-PAT00038
(20)
도 4에는 식(19)에 상응하는 식(20)의 구현을 위해서 필터 계수(p0...pn, 및 -d1...-dn)를 위한 기억소자, 가산기, 지연소자, 등으로 구성될 수 있는 디지털 신호 프로세서가 개략적으로 도시되었다.
도 5는 음원이 정면에서 330°에서 있는 것과 같은 효과를 제공하는 FIR 필터 머리전달함수에 대하여 차수를 감소시킨 IIR 필터 머리전달함수를 비교한 것이다. IIR 필터차수가 k=32 즉, 32차인 것은 특성이 유사하다. 32차인 경우에 대하여 저음에서 특성이 다소 차이가 있지만 저음은 상기 언급한 바와 같이 음원의 방향을 결정하는데 영향을 주지 않는다. 차수가 더 작은 경우에는 특성이 512차의 FIR 필터와 차이가 있어 32차까지가 적절하다.
상술한 바와 같이 본 발명에 의하여 고안된 IIR필터 머리전달함수를 사용함으로써 기존의 FIR 필터 머리전달함수를 사용한 입체음향재생에 비하여 계산량이 작음으로 인해서 소요 메모리가 작아도 되며, 처리를 위한 프로세서의 성능이 좋을 필요가 없기 때문에 소요 비용을 줄일 수 있다. 그리고, 하나의 컴퓨터에서 여러 가지 작업을 수행할 때, 입체음향 재생에 대한 부담을 줄일 수 있기 때문에 효과적이다.
도 1은 머리전달함수를 이용한 입체음향 구현 방법
도 2는 FIR 필터계수로부터 IIR 필터 계수를 구하는 절차
도 3은 디지털신호처리기(DSP)를 이용한 IIR필터 입체음향 구현 구성도
도 4는 IIR 필터 구현을 위한 디지털 신호 처리기의 개략도
도 5는 FIR 필터와 IIR 필터 특성 비교도

Claims (6)

  1. 모노 음원을 수신하는 제1 채널 및 제2 채널을 구비하고,
    상기 제1 및 제2 채널들 각각은 제1 위치 조정 머리관련 전달함수를 구현한 제1 무한 임펄스 응답 필터 및 제1 위치 조정 머리관련 전달함수를 구현한 제1 무한 임펄스 응답 필터를 탑재한 것이며;
    상기 무한 임펄스 응답 필터의 전달함수는
    Figure 112008087512184-PAT00039
    에 의해 구현되는, 입체음향 장치.
  2. 제1항에 있어서, 상기 무한 임펄스 응답 필터 계수는 32차인, 입체음향 장치.
  3. 마이크로부터 수신된 음을 디지털화된 신호로 제공하는 아날로그-디지털 변환기;
    제1항의 상기 무한 임펄스 응답(IIR) 필터가 탑재되는 디지털 신호 프로세서;
    음향 데이터 저장 메모리;
    음원위치 IIR 필터 계수 저장 메모리; 및
    스피커를 포함하는, 입체음향 장치.
  4. 제3항에 있어서, 상기 무한 임펄스 응답 필터 계수는 32차인, 입체음향 장치.
  5. 모노 음원으로부터 수신된 입력신호에 대해 입체음향을 생성하기 위해 머리관련 전달함수를 무한 임펄스 응답필터로 구현하도록 한 디지털 신호 프로세서를 포함하고,
    상기 무한 임펄스 응답 필터의 전달함수는
    Figure 112008087512184-PAT00040
    에 의해 구현되는, 입체음향 재생 전용 디지털 신호 프로세서.
  6. 제5항에 있어서, 상기 무한 임펄스 응답 필터 계수는 32차인, 입체음향 재생 전용 디지털 신호 프로세서.
KR1020080130397A 2008-12-19 2008-12-19 무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서 KR20100071617A (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020080130397A KR20100071617A (ko) 2008-12-19 2008-12-19 무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020080130397A KR20100071617A (ko) 2008-12-19 2008-12-19 무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서

Publications (1)

Publication Number Publication Date
KR20100071617A true KR20100071617A (ko) 2010-06-29

Family

ID=42369088

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020080130397A KR20100071617A (ko) 2008-12-19 2008-12-19 무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서

Country Status (1)

Country Link
KR (1) KR20100071617A (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2549826A (en) * 2016-02-18 2017-11-01 Google Inc Signal processing methods and systems for rendering audio on virtual loudspeaker arrays

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2549826A (en) * 2016-02-18 2017-11-01 Google Inc Signal processing methods and systems for rendering audio on virtual loudspeaker arrays
US10142755B2 (en) 2016-02-18 2018-11-27 Google Llc Signal processing methods and systems for rendering audio on virtual loudspeaker arrays
GB2549826B (en) * 2016-02-18 2020-02-19 Google Llc Signal processing methods and systems for rendering audio on virtual loudspeaker arrays

Similar Documents

Publication Publication Date Title
JP4780119B2 (ja) 頭部伝達関数測定方法、頭部伝達関数畳み込み方法および頭部伝達関数畳み込み装置
KR100619082B1 (ko) 와이드 모노 사운드 재생 방법 및 시스템
KR100677119B1 (ko) 와이드 스테레오 재생 방법 및 그 장치
CN100574516C (zh) 对多声道声音模拟2声道虚拟声音的方法和装置
KR20080092912A (ko) 수음·재생 방법 및 장치
US9607622B2 (en) Audio-signal processing device, audio-signal processing method, program, and recording medium
CN102256192A (zh) 声音信号的分离
US7826630B2 (en) Sound image localization apparatus
JP6515720B2 (ja) 頭外定位処理装置、頭外定位処理方法、及びプログラム
CN104010264A (zh) 双声道音频信号处理的方法和装置
US8320590B2 (en) Device, method, program, and system for canceling crosstalk when reproducing sound through plurality of speakers arranged around listener
JPH08182100A (ja) 音像定位方法および音像定位装置
JP2008502200A (ja) ワイドステレオ再生方法及びその装置
JP2011211312A (ja) 音像定位処理装置及び音像定位処理方法
JP5754595B2 (ja) トランスオーラルシステム
JPH09327099A (ja) 音響再生装置
US8958585B2 (en) Sound image localization apparatus
CN112567766B (zh) 信号处理装置、信号处理方法和介质
JP2005167381A (ja) デジタル信号処理装置及びデジタル信号処理方法、並びにヘッドホン装置
KR20100071617A (ko) 무한 임펄스 응답 필터 머리전달함수를 이용한 입체음향 장치 및 이에 사용되는 디지털 신호 프로세서
JP5776223B2 (ja) 音像制御装置および音像制御方法
CN109923877B (zh) 对立体声音频信号进行加权的装置和方法
KR100641454B1 (ko) 오디오 시스템의 크로스토크 제거 장치
JP2011259299A (ja) 頭部伝達関数生成装置、頭部伝達関数生成方法及び音声信号処理装置
JP5163685B2 (ja) 頭部伝達関数測定方法、頭部伝達関数畳み込み方法および頭部伝達関数畳み込み装置

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
SUBM Surrender of laid-open application requested