Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7675989B2 - GENERATION DEVICE, GENERATION METHOD, AND PROGRAM - Google Patents
[go: Go Back, main page]

JP7675989B2 - GENERATION DEVICE, GENERATION METHOD, AND PROGRAM - Google Patents

GENERATION DEVICE, GENERATION METHOD, AND PROGRAM Download PDF

Info

Publication number
JP7675989B2
JP7675989B2 JP2022025303A JP2022025303A JP7675989B2 JP 7675989 B2 JP7675989 B2 JP 7675989B2 JP 2022025303 A JP2022025303 A JP 2022025303A JP 2022025303 A JP2022025303 A JP 2022025303A JP 7675989 B2 JP7675989 B2 JP 7675989B2
Authority
JP
Japan
Prior art keywords
signal
microphone
noise
sound
estimated
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
JP2022025303A
Other languages
Japanese (ja)
Other versions
JP2023121933A (en
Inventor
弘章 伊藤
伸 村田
記良 鎌土
直樹 信夫
陽一 羽田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
THE UNIVERSITY OF ELECTRO-COMUNICATINS
NTT Inc
NTT Inc USA
Original Assignee
THE UNIVERSITY OF ELECTRO-COMUNICATINS
Nippon Telegraph and Telephone Corp
NTT Inc USA
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 THE UNIVERSITY OF ELECTRO-COMUNICATINS, Nippon Telegraph and Telephone Corp, NTT Inc USA filed Critical THE UNIVERSITY OF ELECTRO-COMUNICATINS
Priority to JP2022025303A priority Critical patent/JP7675989B2/en
Publication of JP2023121933A publication Critical patent/JP2023121933A/en
Application granted granted Critical
Publication of JP7675989B2 publication Critical patent/JP7675989B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Description

本発明は、特定の位置での外部の騒音を抑圧する能動的騒音抑圧(ANC:Active Noise Control)の技術に関する。 The present invention relates to active noise control (ANC) technology that suppresses external noise at a specific position.

従来の能動的騒音抑圧技術として非特許文献1が知られている。能動的騒音抑制では、参照マイク、エラーマイク、キャンセルスピーカを一般的に用いる。図1は従来の騒音抑圧装置の構成例を示す。参照マイク91で騒音源の発する騒音を収音する。キャンセルスピーカ92は、抑圧信号生成装置90で生成されたキャンセル信号を再生して、騒音を相殺するキャンセル音を発する。さらに、エラーマイク93で騒音の消し残しを収音し、フィードバックする。抑圧信号生成装置90は、参照マイク91の収音信号とエラーマイク93の収音信号とを用いて、騒音の消し残しが小さくなるようにキャンセル信号を能動的に制御し、生成する。エラーマイク93の設置位置において、騒音の消し残しが小さくなるように、キャンセルスピーカ92がキャンセル音を発するため、キャンセル音はエラーマイク93の設置位置で最も効率よく騒音を抑圧する。そのため、エラーマイク93はユーザの耳元近くに設置される。 Non-Patent Document 1 is known as a conventional active noise suppression technology. In active noise suppression, a reference microphone, an error microphone, and a cancellation speaker are generally used. FIG. 1 shows an example of the configuration of a conventional noise suppression device. The reference microphone 91 picks up noise emitted by a noise source. The cancellation speaker 92 reproduces the cancellation signal generated by the suppression signal generation device 90 and emits a cancellation sound that cancels the noise. Furthermore, the error microphone 93 picks up the remaining noise and feeds it back. The suppression signal generation device 90 actively controls and generates a cancellation signal using the collected signal of the reference microphone 91 and the collected signal of the error microphone 93 so that the remaining noise is reduced. Since the cancellation speaker 92 emits the cancellation sound so that the remaining noise is reduced at the installation position of the error microphone 93, the cancellation sound most efficiently suppresses noise at the installation position of the error microphone 93. For this reason, the error microphone 93 is installed close to the user's ear.

梶川, 「アクティブノイズコントロールの最近の話題と応用」, 研究報告音楽情報科学(MUS), vol. 2015-MUS-107, no. 3, pp. 1-6, 5月 2015.Kajikawa, "Recent Topics and Applications of Active Noise Control", Research Report Music Information Science (MUS), vol. 2015-MUS-107, no. 3, pp. 1-6, May 2015.

しかしながら、実際の利用に際しては、エラーマイク93をユーザの耳元近くに設置できない場合もあり、エラーマイク93の設置位置とユーザの耳元との距離が大きくなると、前述の通り、エラーマイク93の設置位置で最も効率よく騒音を抑圧し、ユーザの耳元では騒音の消し残しが大きくなり、抑圧性能が低下し、ユーザが騒音抑圧の恩恵を十分に得られない場合がある。例えば、騒音源から耳元までの距離が100mmであり、エラーマイク93をユーザの耳元(0mm)に設置した場合の抑圧性能は-∞dBであり、エラーマイク93を騒音源と耳元との中間地点に設置した場合の抑圧性能は-7.38dBであることをシミュレーションにて確認した。図2は、従来技術の抑圧可能領域(スイートスポット)Sと所望のスイートスポットSとの違いを説明するための図である。 However, in actual use, there are cases where the error microphone 93 cannot be installed close to the user's ear, and if the distance between the installation position of the error microphone 93 and the user's ear is large, as described above, the noise is most efficiently suppressed at the installation position of the error microphone 93, and the noise remaining at the user's ear increases, the suppression performance decreases, and the user may not be able to fully benefit from the noise suppression. For example, it was confirmed by simulation that when the distance from the noise source to the ear is 100 mm and the error microphone 93 is installed at the user's ear (0 mm), the suppression performance is -∞ dB, and when the error microphone 93 is installed at the midpoint between the noise source and the ear, the suppression performance is -7.38 dB. Figure 2 is a diagram for explaining the difference between the suppressible region (sweet spot) S 1 of the conventional technology and the desired sweet spot S 2 .

本発明は、実際のエラーマイクの配置位置で収音した収音信号から、ユーザの耳元で収音した場合に得られる収音信号を推定し、キャンセル信号を能動的に制御するために、実際のエラーマイクの配置位置で収音した収音信号の代わりに推定により得られた収音信号を使用することで、実際のエラーマイクから離れたユーザの耳元でも高い抑圧性能を実現する生成装置、生成方法、そのプログラムを提供することを目的とする。 The present invention aims to provide a generation device, generation method, and program that estimates the sound signal that would be obtained if sound were picked up near the user's ear from the sound signal picked up at the actual error microphone placement position, and uses the estimated sound signal instead of the sound signal picked up at the actual error microphone placement position to actively control the cancellation signal, thereby achieving high suppression performance even near the user's ear away from the actual error microphone.

上記の課題を解決するために、本発明の一態様によれば、生成装置は、アクティブノイズコントロールに用いるキャンセル信号を生成する。生成装置は、騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成する。 In order to solve the above problem, according to one aspect of the present invention, a generating device generates a cancellation signal to be used in active noise control. The generating device generates the cancellation signal so that the point at which the amount of noise suppression is maximized is located closer to the user than the installation position of the error microphone.

本発明によれば、ユーザの耳元にエラーマイクを配置できない場合に従来よりも高い抑圧性能を実現することができるという効果を奏する。 The present invention has the effect of achieving higher suppression performance than conventional techniques when an error microphone cannot be placed close to the user's ear.

従来のアクティブノイズコントロールを説明するための図。FIG. 1 is a diagram for explaining a conventional active noise control. 従来技術の抑圧可能領域を説明するための図。FIG. 1 is a diagram for explaining a suppressible region of the conventional technique. 第一実施形態に係る騒音抑圧システムの機能ブロック図。1 is a functional block diagram of a noise suppression system according to a first embodiment. 第一実施形態に係る騒音抑圧システムの処理フローの例を示す図。FIG. 2 is a diagram showing an example of a processing flow of the noise suppression system according to the first embodiment. 推定方法1を説明するための図。FIG. 13 is a diagram for explaining an estimation method 1. 推定方法2を説明するための図。FIG. 13 is a diagram for explaining an estimation method 2. 推定方法3を説明するための図。FIG. 13 is a diagram for explaining estimation method 3. 第一実施形態のシミュレーション結果を示す図。FIG. 4 is a diagram showing a simulation result of the first embodiment. 本手法を適用するコンピュータの構成例を示す図。FIG. 13 is a diagram showing an example of the configuration of a computer to which the present technique is applied.

以下、本発明の実施形態について、説明する。なお、以下の説明に用いる図面では、同じ機能を持つ構成部や同じ処理を行うステップには同一の符号を記し、重複説明を省略する。以下の説明において、テキスト中で使用する記号「^」「-」等は、本来直後の文字の真上に記載されるべきものであるが、テキスト記法の制限により、当該文字の直前に記載する。式中においてはこれらの記号は本来の位置に記述している。また、ベクトルや行列の各要素単位で行われる処理は、特に断りが無い限り、そのベクトルやその行列の全ての要素に対して適用されるものとする。 Hereinafter, an embodiment of the present invention will be described. In the drawings used in the following description, components having the same functions and steps performing the same processing are denoted with the same reference numerals, and duplicated description will be omitted. In the following description, symbols such as "^" and " - " used in the text should be written directly above the character immediately following them, but due to limitations in text notation, they are written immediately before the character in question. In the formulas, these symbols are written in their original positions. Furthermore, unless otherwise specified, processing performed on each element of a vector or matrix is assumed to be applied to all elements of that vector or matrix.

<第一実施形態のポイント>
本実施形態では、耳元から離れた位置に設置されたエラーマイクの収音信号から耳元での観測音圧を推定する。例えば、実際のエラーマイクの収音信号から、耳元に配置された仮想的なエラーマイクの収音信号を推定し、ANCにおいて、仮想的なエラーマイクの収音信号を従来のエラーマイクの収音信号として用いる。このような構成とすることで、スイートスポットの位置をエラーマイクの設置位置から仮想的なエラーマイクの位置に変更し、耳元での消し残りをキャンセルする音を出すことができる。
<Key Points of the First Embodiment>
In this embodiment, the observed sound pressure at the ear is estimated from the pickup signal of the error microphone installed at a position away from the ear. For example, the pickup signal of a virtual error microphone placed at the ear is estimated from the pickup signal of the actual error microphone, and the pickup signal of the virtual error microphone is used as the pickup signal of the conventional error microphone in the ANC. With this configuration, the position of the sweet spot can be changed from the installation position of the error microphone to the position of the virtual error microphone, and a sound that cancels the residual sound at the ear can be produced.

仮想的なエラーマイクの収音信号を推定する方法としては様々な方法が考えられる。例えば、実際のエラーマイクの設置位置から耳元までの距離減衰、位相遅延を考慮して音圧を推定する。また、例えば、球面上に配置した実際のエラーマイクから、球面調和関数を用いて耳元の音圧を推定する。 There are various methods for estimating the sound pickup signal of a virtual error microphone. For example, the sound pressure is estimated by taking into account the distance attenuation and phase delay from the actual error microphone's installation position to the ear. In addition, for example, the sound pressure at the ear is estimated using spherical harmonics from an actual error microphone placed on a sphere.

<第一実施形態>
図3は第一実施形態に係る騒音抑圧システムの機能ブロック図を、図4はその処理フローを示す。
First Embodiment
FIG. 3 is a functional block diagram of the noise suppression system according to the first embodiment, and FIG. 4 shows the processing flow thereof.

騒音抑圧システムは、参照マイク91、キャンセルスピーカ92、エラーマイク93、抑圧信号生成部110および音圧推定部120を含む。抑圧信号生成部110および音圧推定部120からなる装置を抑圧信号生成装置ともいう。 The noise suppression system includes a reference microphone 91, a cancellation speaker 92, an error microphone 93, a suppression signal generation unit 110, and a sound pressure estimation unit 120. The device consisting of the suppression signal generation unit 110 and the sound pressure estimation unit 120 is also called a suppression signal generation device.

抑圧信号生成装置は、参照マイク91の収音信号x(r)と、エラーマイク93の収音信号x(e)とを入力とし、騒音の抑圧量が最大となる地点が、エラーマイク93の設置位置よりもユーザ側に位置するようにキャンセル信号(以下、「抑圧信号」ともいう)yを生成して、キャンセルスピーカ92に出力する。 The suppression signal generating device receives the pickup signal x(r) of the reference microphone 91 and the pickup signal x(e) of the error microphone 93 as input, generates a cancellation signal (hereinafter also referred to as the "suppression signal") y so that the point at which the amount of noise suppression is maximum is located closer to the user than the installation position of the error microphone 93, and outputs the cancellation signal to the cancellation speaker 92.

抑圧信号生成装置は、例えば、中央演算処理装置(CPU: Central Processing Unit)、主記憶装置(RAM: Random Access Memory)などを有する公知又は専用のコンピュータに特別なプログラムが読み込まれて構成された特別な装置である。抑圧信号生成装置は、例えば、中央演算処理装置の制御のもとで各処理を実行する。抑圧信号生成装置に入力されたデータや各処理で得られたデータは、例えば、主記憶装置に格納され、主記憶装置に格納されたデータは必要に応じて中央演算処理装置へ読み出されて他の処理に利用される。抑圧信号生成装置の各処理部は、少なくとも一部が集積回路等のハードウェアによって構成されていてもよい。抑圧信号生成装置が備える各記憶部は、例えば、RAM(Random Access Memory)などの主記憶装置、またはリレーショナルデータベースやキーバリューストアなどのミドルウェアにより構成することができる。ただし、各記憶部は、必ずしも抑圧信号生成装置がその内部に備える必要はなく、ハードディスクや光ディスクもしくはフラッシュメモリ(Flash Memory)のような半導体メモリ素子により構成される補助記憶装置により構成し、抑圧信号生成装置の外部に備える構成としてもよい。 The suppression signal generating device is a special device configured by loading a special program into a publicly known or dedicated computer having, for example, a central processing unit (CPU), a main memory (RAM), etc. The suppression signal generating device executes each process under the control of the central processing unit, for example. Data input to the suppression signal generating device and data obtained by each process are stored in, for example, the main memory, and the data stored in the main memory is read out to the central processing unit as necessary and used for other processes. At least a part of each processing unit of the suppression signal generating device may be configured by hardware such as an integrated circuit. Each storage unit of the suppression signal generating device may be configured by, for example, a main storage device such as a RAM (Random Access Memory), or middleware such as a relational database or a key-value store. However, each storage unit does not necessarily have to be provided inside the suppression signal generating device, and may be configured by an auxiliary storage device configured by a semiconductor memory element such as a hard disk, optical disk, or flash memory, and may be configured to be provided outside the suppression signal generating device.

以下、各部について説明する。 Each part is explained below.

<参照マイク91>
参照マイク91は、抑圧対象の音を収音し(S91)、収音信号x(r)を出力する。参照マイク91で収音した抑圧対象の音を、以下「騒音」と記載する。
<Reference Microphone 91>
The reference microphone 91 picks up the sound to be suppressed (S91) and outputs a picked-up sound signal x(r). The sound to be suppressed picked up by the reference microphone 91 is hereinafter referred to as "noise."

<キャンセルスピーカ92>
キャンセルスピーカ92は、キャンセル信号yを入力とし、キャンセル信号yを再生する(S92)。キャンセルスピーカ92から再生される再生音と抑圧対象の騒音とが完全な逆位相となる場合、再生音と抑圧対象の騒音とが重なる、すなわち、音波同士が重畳する、と波が打ち消し合うため、騒音が抑圧される。
<Cancellation Speaker 92>
The cancellation speaker 92 receives the cancellation signal y as input and reproduces the cancellation signal y (S92). When the playback sound reproduced from the cancellation speaker 92 and the noise to be suppressed are in completely opposite phase, the playback sound and the noise to be suppressed overlap, that is, the sound waves overlap, and the waves cancel each other out, thereby suppressing the noise.

<エラーマイク93>
エラーマイク93は、騒音の消し残しを含む、キャンセルスピーカ92から再生される再生音で抑圧されなかった音を収音し(S93)、収音信号x(e)を出力する。エラーマイク93は、観測点(例えば、ユーザの耳元)よりも騒音源に近い位置に配置される。例えば、エラーマイク93は、図5のように耳元よりも騒音源に0.05m近い位置に配置される。
<Error Microphone 93>
The error microphone 93 picks up sounds that are not suppressed by the sound reproduced from the cancellation speaker 92, including residual noise (S93), and outputs a picked-up sound signal x(e). The error microphone 93 is placed at a position closer to the noise source than the observation point (e.g., the user's ear). For example, the error microphone 93 is placed at a position 0.05 m closer to the noise source than the user's ear, as shown in FIG. 5.

<音圧推定部120>
音圧推定部120は、エラーマイク93の出力信号(収音信号)x(e)を入力とし、エラーマイク93よりも観測点に近い位置にマイク130を設置した場合に収音されると推定される信号である、推定収音信号x(v)を算出し、出力する。すなわち、音圧推定部120は、キャンセルスピーカ92から再生される再生音で抑圧されなかった音がマイク130の設置位置で収音される場合に得られる収音信号を推定し(S120)、推定した収音信号を推定収音信号x(v)として出力する。以下、推定収音信号x(v)の推定方法を3つ例示する。ここで、マイク130は実際には設置せず仮想的に設置されるものであり、以下仮想マイク130と記載する。
<Sound Pressure Estimation Unit 120>
The sound pressure estimation unit 120 receives the output signal (sound pickup signal) x(e) of the error microphone 93 as input, calculates an estimated sound pickup signal x(v) which is a signal estimated to be picked up when the microphone 130 is installed at a position closer to the observation point than the error microphone 93, and outputs the signal. That is, the sound pressure estimation unit 120 estimates a sound pickup signal obtained when a sound that is not suppressed by the sound reproduced from the cancellation speaker 92 is picked up at the installation position of the microphone 130 (S120), and outputs the estimated sound pickup signal as an estimated sound pickup signal x(v). Three examples of estimation methods for the estimated sound pickup signal x(v) are given below. Here, the microphone 130 is not actually installed but is installed virtually, and is hereinafter referred to as a virtual microphone 130.

(推定方法1)
本推定方法では、エラーマイク93と仮想マイク130から、距離減衰と位相遅延に基づき、実際のエラーマイク93の収音信号x(e)から仮想マイク130の収音信号x(v)を推定する。図5は、騒音源とエラーマイク93、仮想マイク130の位置関係を説明するための図である。
(Estimation method 1)
In this estimation method, a pickup signal x(v) of the virtual microphone 130 is estimated from an actual pickup signal x(e) of the error microphone 93 based on distance attenuation and phase delay from the error microphone 93 and the virtual microphone 130. Fig. 5 is a diagram for explaining the positional relationship between the noise source, the error microphone 93, and the virtual microphone 130.

本推定方法では、騒音源の位置を仮定し、騒音源からエラーマイク93および仮想マイク130に騒音が平面波で伝搬すると仮定する。騒音源からエラーマイク93までの伝達関数と、騒音源から観測点(仮想マイク130の位置)までの伝達関数から、エラーマイク93から観測点までの距離減衰および位相ずれを推定することで、仮想マイク130で収音される収音信号を推定する。音圧推定部120は、エラーマイク93の出力信号(収音信号)x(e)から、次式により、仮想マイク130の収音信号を推定し、推定収音信号x(v)=[^Gp1 ^Gp2]を出力する。 In this estimation method, the position of a noise source is assumed, and it is assumed that noise propagates as a plane wave from the noise source to the error microphone 93 and the virtual microphone 130. The sound signal picked up by the virtual microphone 130 is estimated by estimating the distance attenuation and phase shift from the error microphone 93 to the observation point from the transfer function from the noise source to the error microphone 93 and the transfer function from the noise source to the observation point (position of the virtual microphone 130). The sound pressure estimation unit 120 estimates the sound signal picked up by the virtual microphone 130 from the output signal (sound pickup signal) x(e) of the error microphone 93 by the following equation, and outputs the estimated sound pickup signal x(v)=[^G p1 ^G p2 ].

^Gpn=wnx(e) (n=1,2) (1)
ここで、ゲイン減衰のみを考慮するとwn
wn=|Gpn|/|Ge| (2)
であり、位相ずれのみを考慮するとwn
wn=exp((arg Gpn-arg Ge)j) (3)
である。式(2),(3)のGe、Gpnは、仮定した騒音源の位置と観測点から推定処理に先立ち予め算出しておく。例えば、仮定した騒音源の位置に騒音源用のスピーカを配置し、騒音源用のスピーカで所定の信号を再生し、エラーマイクの位置に配置したマイクで収音した収音信号からGeを求め、観測点の位置に配置したマイクで収音した収音信号からGpnを求めておく。
(推定方法2)
本推定方法では、頭部近傍に等間隔に配置した複数のエラーマイクの収音信号から、球面調和関数展開係数を利用して仮想的なエラーマイクの収音信号を推定する。図6は、実際のエラーマイクの位置関係を説明するための図である。
^G pn =w n x(e) (n=1,2) (1)
Here, if we only consider the gain attenuation, wn is
wn = | Gpn |/| Ge | (2)
If we only consider the phase shift, wn is
w n =exp((arg G pn -arg G e )j) (3)
G e and G pn in formulas (2) and (3) are calculated in advance from the assumed noise source position and observation point prior to estimation processing. For example, a loudspeaker for the noise source is placed at the assumed noise source position, a predetermined signal is played back by the loudspeaker for the noise source, G e is obtained from a signal picked up by a microphone placed at the error microphone position, and G pn is obtained from a signal picked up by a microphone placed at the observation point position.
(Estimation method 2)
In this estimation method, a signal picked up by a virtual error microphone is estimated from signals picked up by multiple error microphones arranged at equal intervals near the head, using spherical harmonic expansion coefficients. Fig. 6 is a diagram for explaining the positional relationship of the actual error microphones.

本推定方法では、半径reの球面上に等間隔にエラーマイクを配置し、半径rの球面上の音圧を推定する。例えば、中心からエラーマイクまでの距離をre=0.15mとし、(i)6個のエラーマイクを正六面体の各面の中心に配置する(図6の(i)参照)、(ii)12個のエラーマイクを正十二面体の各面の中心に配置する(図6の(ii)参照)ことで、等間隔にエラーマイクを配置することができる。例えば、中心から観測点(仮想マイク130の位置)までの距離をr=0.08mとして推定する。 In this estimation method, error microphones are arranged at equal intervals on a sphere of radius r e , and the sound pressure on the sphere of radius r is estimated. For example, the distance from the center to the error microphones is set to r e =0.15 m, and the error microphones can be arranged at equal intervals by (i) arranging six error microphones at the center of each face of a regular hexahedron (see (i) of FIG. 6 ), or (ii) arranging twelve error microphones at the center of each face of a regular dodecahedron (see (ii) of FIG. 6 ). For example, the distance from the center to the observation point (the position of the virtual microphone 130) is estimated as r=0.08 m.

球面調和関数展開を利用することで、ある球面上での観測音圧から、任意の球面上における観測音圧を推定することが可能である。 By using spherical harmonic expansion, it is possible to estimate the observed sound pressure on any sphere from the observed sound pressure on a certain sphere.

半径re上のL個のエラーマイクから音圧の観測値p(θ11),p(θ22),…,p(θLL)を得る。例えば、L個のエラーマイク93の収音信号x(e)=[p(θ11),p(θ22),…,p(θLL)]とする。 Observation values of sound pressure p( θ1 , φ1 ),p( θ2 , φ2 ),...,p( θL , φL ) are obtained from L error microphones on a radius r e . For example, the sound pickup signals x(e) of the L error microphones 93 are set to [p( θ1 , φ1 ),p( θ2 , φ2 ),...,p( θL , φL )].

音圧推定部120は、次式により、球面調和関数Ym n(・)に対する半径re上の音場係数Pnm(re)を求める。

Figure 0007675989000001

音圧推定部120は、求めた音場係数Pnm(re)を用いて、次式により、半径r上の音場係数Pnm(r)を求める。
Figure 0007675989000002

ただし、kは波数であり、aは音を反射する剛球の半径とし、jnはn次の球面ベッセル関数であり、jn'はjnの微分であり、hn (2)は第二種球ハンケル関数であり、hn'(2)はhn (2)の微分である。
音圧推定部120は、再合成により、観測点(r,θ,φ)における音圧の推定値^p(r,θ,φ)を得る。
Figure 0007675989000003

なお、推定収音信号x(v)=^p(r,θ,φ)とする。 The sound pressure estimation unit 120 obtains a sound field coefficient P nm (r e ) on a radius r e for the spherical harmonic function Y m n (·) by the following equation.
Figure 0007675989000001

The sound pressure estimation unit 120 uses the obtained sound field coefficient P nm (r e ) to obtain a sound field coefficient P nm (r) on a radius r according to the following equation.
Figure 0007675989000002

where k is the wave number, a is the radius of the hard sphere that reflects sound, j n is the nth order spherical Bessel function, j n ' is the derivative of j n , h n (2) is the second kind of spherical Hankel function, and h n ' (2) is the derivative of h n (2) .
The sound pressure estimation unit 120 obtains an estimate value ^p(r, θ, φ) of the sound pressure at the observation point (r, θ, φ) by resynthesis.
Figure 0007675989000003

The estimated picked-up signal is x(v)=^p(r,θ,φ).

以下、式(5)の導出について説明する。 The derivation of equation (5) is explained below.

騒音源を点音源とし、半径aの剛球での反射を考慮したとき、点(r,θ,φ)における音圧は、

Figure 0007675989000004

である。なお、Bnmは騒音源の座標と信号で定まる係数である。球面調和関数展開は、
Figure 0007675989000005

であり、音場係数Pnm(r)、音場係数Pnm(re)は、次式で表される。
Figure 0007675989000006

Figure 0007675989000007

式(10)をBnm=…の形に変形して(9)に代入すると、
Figure 0007675989000008

Figure 0007675989000009

となる。
Figure 0007675989000010

なお、球面調和関数展開における最大次数Nは以下の制約を受ける。 When the noise source is a point source and reflection from a hard sphere of radius a is taken into consideration, the sound pressure at point (r, θ, φ) is
Figure 0007675989000004

Note that B nm is a coefficient determined by the coordinates of the noise source and the signal. The spherical harmonic expansion is
Figure 0007675989000005

and the sound field coefficients P nm (r) and P nm (r e ) are expressed by the following equations.
Figure 0007675989000006

Figure 0007675989000007

Transforming equation (10) into the form B nm = ... and substituting it into (9), we get
Figure 0007675989000008

Figure 0007675989000009

It becomes.
Figure 0007675989000010

The maximum order N in the spherical harmonic expansion is subject to the following constraints:

(N+1)2<L
ここで、球面調和関数Ym n(・)の各モードに対応するだけのスピーカ数が必要である。L=6ならばN=1であり、L=12ならばN=2である。
(N+1) 2 <L
Here, the number of speakers required corresponds to each mode of the spherical harmonic function Y m n (·). If L=6, then N=1, and if L=12, then N=2.

さらに、空間エイリアシングが起こらない条件として、Nは以下の制約を受ける。 In addition, to prevent spatial aliasing, N is subject to the following constraints:

kr<N
頭部と仮想的なエラーマイク間の距離が制限される。例えば、周波数300Hzとして、N=1のとき推定可能領域は頭部との距離r=0.18m以内に制限される。
kr<N
The distance between the head and the virtual error microphone is limited. For example, when the frequency is 300 Hz and N=1, the estimation area is limited to within the distance r=0.18 m from the head.

(推定方法3)
本推定方法では、頭部近傍に非等間隔に配置した複数のエラーマイクの収音信号から、最小二乗法で推定された球面調和関数展開係数を利用して仮想的なエラーマイクの収音信号を推定する。図7は、エラーマイクの位置関係を説明するための図である。例えば、頭部後方4点(方位角4つ(0°,30°,150°,180°)×仰角0°)または頭部後方12点(方位角4つ(0°,30°,150°,180°)×仰角3つ(-30°,0°,30°))にエラーマイクを配置する。エラーマイクの設置半径は、設置する環境や抑圧したい音の種類により決定すればよい。例えば、走行音を抑圧したい鉄道などに用いる場合、座席サイズを考慮して、エラーマイクの設置半径は0.13mとする。
(Estimation method 3)
In this estimation method, the signal picked up by a virtual error microphone is estimated from the signal picked up by a plurality of error microphones arranged at non-uniform intervals near the head, using spherical harmonic expansion coefficients estimated by the least squares method. FIG. 7 is a diagram for explaining the positional relationship of the error microphones. For example, error microphones are arranged at four points behind the head (four azimuth angles (0°, 30°, 150°, 180°) × elevation angle 0°) or at 12 points behind the head (four azimuth angles (0°, 30°, 150°, 180°) × three elevation angles (-30°, 0°, 30°)). The installation radius of the error microphones may be determined according to the environment in which they are installed and the type of sound to be suppressed. For example, when used in a train where running noise is to be suppressed, the installation radius of the error microphones is set to 0.13 m, taking into account the seat size.

本推定方法では、半径reの球面上に非等間隔にエラーマイクを配置し、半径rの球面上の音圧を推定する。本推定方法では、球面調和関数展開を直接利用できないので、最小二乗法により球面調和関数展開係数を推定し、半径rの球面上の音圧を得る。 In this estimation method, error microphones are placed at non-uniform intervals on a sphere of radius r e , and the sound pressure on the sphere of radius r is estimated. Since this estimation method cannot directly use spherical harmonic function expansion, the spherical harmonic function expansion coefficients are estimated by the least squares method to obtain the sound pressure on the sphere of radius r.

半径re上のL個のエラーマイクから音圧の観測値p(re11),p(re22),…,p(reLL)を得る。エラーマイク93の収音信号x(e)=-p=[p(re11),p(re22),…,p(reLL)]Tとする。-pは以下のように表される。




Figure 0007675989000011

Figure 0007675989000012

なお、ωi=ii)とし、i=1,2,…,Lとする。 Observation values of sound pressure p(r e11 ),p(r e22 ),...,p(r eLL ) are obtained from L error microphones on a radius of r e . The sound pickup signal of the error microphone 93 is defined as x(e)= - p=[p(r e11 ),p(r e22 ),...,p(r eLL )] T. -p is expressed as follows.




Figure 0007675989000011

Figure 0007675989000012

Here, ω i =i , φ i ), where i = 1, 2, ..., L.

音圧推定部120は、絶対値の二乗誤差が最小となる-P(re)を解として求める。

Figure 0007675989000013

推定したい点が半径r上にあるとしたとき、展開係数Pnm(r)および推定値の計算は推定方法2と同様である。つまり、
音圧推定部120は、次式により、球面調和関数Ym n(・)に対する半径re上の音場係数Pnm(re)を求める。
Figure 0007675989000014

音圧推定部120は、求めた音場係数Pnm(re)を用いて、次式により、半径r上の音場係数Pnm(r)を求める。
Figure 0007675989000015

音圧推定部120は、再合成により、観測点(r,θ,φ)における音圧の推定値^p(r,θ,φ)を得る。
Figure 0007675989000016

なお、推定収音信号x(v)=^p(r,θ,φ)とする。 The sound pressure estimation unit 120 obtains, as a solution, −P (r e ) that minimizes the absolute square error.
Figure 0007675989000013

When the point to be estimated is on the radius r, the calculation of the expansion coefficients P nm (r) and the estimated value is the same as in Estimation Method 2. That is,
The sound pressure estimation unit 120 obtains a sound field coefficient P nm (r e ) on a radius r e for the spherical harmonic function Y m n (·) by the following equation.
Figure 0007675989000014

The sound pressure estimation unit 120 uses the obtained sound field coefficient P nm (r e ) to obtain a sound field coefficient P nm (r) on a radius r according to the following equation.
Figure 0007675989000015

The sound pressure estimation unit 120 obtains an estimate ^p(r, θ, φ) of the sound pressure at the observation point (r, θ, φ) by resynthesis.
Figure 0007675989000016

The estimated picked-up signal is x(v)=^p(r,θ,φ).

<抑圧信号生成部110>
抑圧信号生成部110は、収音信号x(r)と推定収音信号x(v)とを入力とし、仮想マイク130の設置位置における騒音を抑圧するためのキャンセル信号yを生成し(S110)、出力する。
<Suppression signal generating unit 110>
The suppression signal generator 110 receives the collected sound signal x(r) and the estimated collected sound signal x(v) as input, generates a cancellation signal y for suppressing noise at the installation position of the virtual microphone 130 (S110), and outputs the cancellation signal y.

キャンセル信号の生成方法としては、従来技術を用いることができる。例えば、非特許文献1の方法を用いることができる。本実施形態では、収音信号x(r)、推定収音信号x(v)とキャンセル信号yによってフィードフォワード型ANCを実現する。騒音源からの騒音とキャンセル信号yの再生音との干渉音を仮想マイク130で検出した際に得られるだろう収音信号を推定するとともに、騒音源からの騒音を参照マイク91で検出し、適応ディジタルフィルタによって実現されている騒音制御フィルタに入力することでキャンセル信号yを生成し、キャンセルスピーカ92で再生する。キャンセル信号yの再生音は、キャンセルスピーカ92から仮想マイク130までの一連の伝達系である二次経路を伝播すると仮定する。そして、仮想マイク130の入力が最小となるように騒音制御フィルタの係数を適応アルゴリズムにより更新する。騒音制御フィルタの係数の更新方法としては従来の更新方法を用いることができるため、説明を省略する。フィードフォワード型ANCにおいては、二次経路を推定した二次経路モデルが二次経路の影響を適応アルゴリズムにおいて補償するため利用される。 Conventional techniques can be used as a method for generating the cancellation signal. For example, the method described in Non-Patent Document 1 can be used. In this embodiment, a feedforward type ANC is realized by the collected sound signal x(r), the estimated collected sound signal x(v), and the cancellation signal y. The collected sound signal that would be obtained when the interference sound between the noise from the noise source and the reproduced sound of the cancellation signal y is detected by the virtual microphone 130 is estimated, and the noise from the noise source is detected by the reference microphone 91. The cancellation signal y is generated by inputting the noise from the noise source to a noise control filter realized by an adaptive digital filter, and is reproduced by the cancellation speaker 92. It is assumed that the reproduced sound of the cancellation signal y propagates through a secondary path, which is a series of transmission systems from the cancellation speaker 92 to the virtual microphone 130. Then, the coefficients of the noise control filter are updated by an adaptive algorithm so that the input to the virtual microphone 130 is minimized. Since a conventional update method can be used as a method for updating the coefficients of the noise control filter, a description will be omitted. In the feedforward type ANC, a secondary path model that estimates the secondary path is used to compensate for the influence of the secondary path in the adaptive algorithm.

<効果>
以上の構成により、ユーザの耳元にエラーマイクを配置できない場合に従来よりも高い抑圧性能を実現することができる。図8は、第一実施形態のシミュレーション結果を示す。(A)は騒音を300Hzの平面波とし、(B)は騒音を100Hzの平面波とする。
<Effects>
With the above configuration, it is possible to achieve higher suppression performance than in the past when it is not possible to place an error microphone close to the user's ear. Figure 8 shows the simulation results of the first embodiment. In (A), the noise is a 300 Hz plane wave, and in (B), the noise is a 100 Hz plane wave.

<その他の変形例>
本発明は上記の実施形態及び変形例に限定されるものではない。例えば、上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。
<Other Modifications>
The present invention is not limited to the above-mentioned embodiment and modified examples. For example, the above-mentioned various processes may be executed not only in chronological order as described, but also in parallel or individually depending on the processing capacity of the device executing the processes or as necessary. In addition, appropriate modifications are possible within the scope of the present invention.

<プログラム及び記録媒体>
上述の各種の処理は、図9に示すコンピュータの記憶部2020に、上記方法の各ステップを実行させるプログラムを読み込ませ、制御部2010、入力部2030、出力部2040などに動作させることで実施できる。
<Program and recording medium>
The various processes described above can be implemented by loading a program that executes each step of the above method into the memory unit 2020 of the computer shown in Figure 9, and operating the control unit 2010, input unit 2030, output unit 2040, etc.

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。 The program describing this processing can be recorded on a computer-readable recording medium. Examples of computer-readable recording media include magnetic recording devices, optical disks, magneto-optical recording media, and semiconductor memories.

また、このプログラムの流通は、例えば、そのプログラムを記録したDVD、CD-ROM等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。 The program may be distributed, for example, by selling, transferring, or lending portable recording media such as DVDs or CD-ROMs on which the program is recorded. Furthermore, the program may be distributed by storing the program in a storage device of a server computer and transferring the program from the server computer to other computers via a network.

このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるASP(Application Service Provider)型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの(コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等)を含むものとする。 A computer that executes such a program, for example, first stores in its own storage device the program recorded on a portable recording medium or the program transferred from a server computer. Then, when executing a process, the computer reads the program stored on its own recording medium and executes the process according to the read program. As another execution form of this program, the computer may read the program directly from the portable recording medium and execute the process according to the program, or may execute the process according to the received program each time a program is transferred from the server computer to this computer. In addition, the server computer may not transfer the program to this computer, but may execute the above-mentioned process by a so-called ASP (Application Service Provider) type service that realizes the processing function only by issuing an execution instruction and obtaining the results. Note that the program in this form includes information used for processing by an electronic computer that is equivalent to a program (such as data that is not a direct command to the computer but has properties that specify the processing of the computer).

また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。 In addition, in this embodiment, the device is configured by executing a specific program on a computer, but at least a portion of the processing may be realized by hardware.

Claims (7)

アクティブノイズコントロールに用いるキャンセル信号を生成する生成装置であって、
騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成し、
前記エラーマイクよりも観測点に近い位置に仮想マイクを設置した場合に収音される収音信号を推定し、推定収音信号x(v)を得る音圧推定部と、
抑圧対象の騒音を収音した収音信号x(r)と前記推定収音信号x(v)とを用いて、前記仮想マイクの設置位置における騒音を抑圧するためのキャンセル信号を生成する抑圧信号生成部とを含み、
前記仮想マイクは、実際には設置されずに仮想的に設置されるマイクであり、
前記音圧推定部は、ユーザの頭部近傍に等間隔に配置した複数の前記エラーマイクの収音信号x(e)から、球面調和関数展開係数を利用して、前記推定収音信号x(v)を得る、
生成装置。
A generating device for generating a cancellation signal for use in active noise control, comprising:
A cancellation signal is generated so that the point where the amount of noise suppression is maximum is located closer to the user than the installation position of the error microphone;
a sound pressure estimation unit that estimates a sound pickup signal that is picked up when a virtual microphone is installed at a position closer to the observation point than the error microphone, and obtains an estimated sound pickup signal x(v);
a suppression signal generation unit that generates a cancellation signal for suppressing noise at an installation position of the virtual microphone by using a pickup signal x(r) that picks up the noise to be suppressed and the estimated pickup signal x(v),
The virtual microphone is a microphone that is virtually installed without actually being installed,
The sound pressure estimation unit obtains the estimated sound pickup signal x(v) from the sound pickup signals x(e) of the plurality of error microphones arranged at equal intervals near the user's head by using spherical harmonic function expansion coefficients.
generator.
アクティブノイズコントロールに用いるキャンセル信号を生成する生成装置であって、
騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成し、
前記エラーマイクよりも観測点に近い位置に仮想マイクを設置した場合に収音される収音信号を推定し、推定収音信号x(v)を得る音圧推定部と、
抑圧対象の騒音を収音した収音信号x(r)と前記推定収音信号x(v)とを用いて、前記仮想マイクの設置位置における騒音を抑圧するためのキャンセル信号を生成する抑圧信号生成部とを含み、
前記仮想マイクは、実際には設置されずに仮想的に設置されるマイクであり、
前記音圧推定部は、前記エラーマイクと前記仮想マイクの位置関係から、距離減衰と位相遅延に基づき、前記エラーマイクの収音信号x(e)から
x(v)=[^G p1 ^G p2 ]
^G pn =w n x(e) (n=1,2)
により、前記推定収音信号x(v)を得、仮定した騒音源の位置からエラーマイクの位置までの伝達特性をG e とし、仮定した騒音源の位置から観測点の位置までの伝達特性をG pn とし、ゲイン減衰のみを考慮するとw n
w n =|G pn |/|G e |
であり、位相ずれのみを考慮するとw n
w n =exp((arg G pn -arg G e )j)
である、
生成装置。
A generating device for generating a cancellation signal for use in active noise control, comprising:
A cancellation signal is generated so that the point where the amount of noise suppression is maximum is located closer to the user than the installation position of the error microphone;
a sound pressure estimation unit that estimates a sound pickup signal that is picked up when a virtual microphone is installed at a position closer to the observation point than the error microphone, and obtains an estimated sound pickup signal x(v);
a suppression signal generation unit that generates a cancellation signal for suppressing noise at an installation position of the virtual microphone by using a pickup signal x(r) that picks up the noise to be suppressed and the estimated pickup signal x(v),
The virtual microphone is a microphone that is virtually installed without actually being installed,
The sound pressure estimation unit estimates a sound pressure from the sound signal x(e) of the error microphone based on distance attenuation and phase delay from the positional relationship between the error microphone and the virtual microphone.
x(v)=[^G p1 ^G p2 ]
^G pn =w n x(e) (n=1,2)
The estimated sound pickup signal x(v) is obtained by the above. The transfer characteristic from the assumed noise source position to the error microphone position is Ge , and the transfer characteristic from the assumed noise source position to the observation point position is Gpn . When only the gain attenuation is taken into consideration, wn is
w n =|G pn |/|G e |
If we only consider the phase shift, wn is
w n =exp((arg G pn -arg G e )j)
That is,
generator.
アクティブノイズコントロールに用いるキャンセル信号を生成する生成装置であって、
騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成し、
前記エラーマイクよりも観測点に近い位置に仮想マイクを設置した場合に収音される収音信号を推定し、推定収音信号x(v)を得る音圧推定部と、
抑圧対象の騒音を収音した収音信号x(r)と前記推定収音信号x(v)とを用いて、前記仮想マイクの設置位置における騒音を抑圧するためのキャンセル信号を生成する抑圧信号生成部とを含み、
前記仮想マイクは、実際には設置されずに仮想的に設置されるマイクであり、
前記音圧推定部は、ユーザの頭部近傍に非等間隔に配置した複数の前記エラーマイクの収音信号x(e)から、最小二乗法で推定された球面調和関数展開係数を利用して、前記推定収音信号x(v)を得る、
生成装置。
A generating device for generating a cancellation signal for use in active noise control, comprising:
A cancellation signal is generated so that the point where the amount of noise suppression is maximum is located closer to the user than the installation position of the error microphone;
a sound pressure estimation unit that estimates a sound pickup signal that is picked up when a virtual microphone is installed at a position closer to the observation point than the error microphone, and obtains an estimated sound pickup signal x(v);
a suppression signal generation unit that generates a cancellation signal for suppressing noise at an installation position of the virtual microphone by using a pickup signal x(r) that picks up the noise to be suppressed and the estimated pickup signal x(v),
The virtual microphone is a microphone that is virtually installed without actually being installed,
The sound pressure estimation unit obtains the estimated sound pickup signal x(v) from the sound pickup signals x(e) of the plurality of error microphones arranged at non-equidistant intervals near the user's head by using spherical harmonic function expansion coefficients estimated by a least squares method.
generator.
アクティブノイズコントロールに用いるキャンセル信号を生成する生成方法であって、
騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成し、
前記エラーマイクよりも観測点に近い位置に仮想マイクを設置した場合に収音される収音信号を推定し、推定収音信号x(v)を得る音圧推定ステップと、
抑圧対象の騒音を収音した収音信号x(r)と前記推定収音信号x(v)とを用いて、前記仮想マイクの設置位置における騒音を抑圧するためのキャンセル信号を生成する抑圧信号生成ステップとを含み、
前記仮想マイクは、実際には設置されずに仮想的に設置されるマイクであり、
前記音圧推定ステップは、ユーザの頭部近傍に等間隔に配置した複数の前記エラーマイクの収音信号x(e)から、球面調和関数展開係数を利用して、前記推定収音信号x(v)を得る、
生成方法。
1. A method for generating a cancellation signal for use in active noise control, comprising:
A cancellation signal is generated so that the point where the amount of noise suppression is maximum is located closer to the user than the installation position of the error microphone;
a sound pressure estimation step of estimating a sound pickup signal that would be picked up if a virtual microphone were installed at a position closer to the observation point than the error microphone, and obtaining an estimated sound pickup signal x(v);
a suppression signal generation step of generating a cancellation signal for suppressing noise at an installation position of the virtual microphone by using a pickup signal x(r) that picks up the noise to be suppressed and the estimated pickup signal x(v),
The virtual microphone is a microphone that is virtually installed without actually being installed,
The sound pressure estimation step obtains the estimated sound pickup signal x(v) from the sound pickup signals x(e) of the plurality of error microphones arranged at equal intervals near the user's head by using spherical harmonic function expansion coefficients.
Generation method.
アクティブノイズコントロールに用いるキャンセル信号を生成する生成方法であって、
騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成し、
前記エラーマイクよりも観測点に近い位置に仮想マイクを設置した場合に収音される収音信号を推定し、推定収音信号x(v)を得る音圧推定ステップと、
抑圧対象の騒音を収音した収音信号x(r)と前記推定収音信号x(v)とを用いて、前記仮想マイクの設置位置における騒音を抑圧するためのキャンセル信号を生成する抑圧信号生成ステップとを含み、
前記仮想マイクは、実際には設置されずに仮想的に設置されるマイクであり、
前記音圧推定ステップは、前記エラーマイクと前記仮想マイクの位置関係から、距離減衰と位相遅延に基づき、前記エラーマイクの収音信号x(e)から
x(v)=[^G p1 ^G p2 ]
^G pn =w n x(e) (n=1,2)
により、前記推定収音信号x(v)を得、仮定した騒音源の位置からエラーマイクの位置までの伝達特性をG e とし、仮定した騒音源の位置から観測点の位置までの伝達特性をG pn とし、ゲイン減衰のみを考慮するとw n
w n =|G pn |/|G e |
であり、位相ずれのみを考慮するとw n
w n =exp((arg G pn -arg G e )j)
である、
生成方法。
1. A method for generating a cancellation signal for use in active noise control, comprising:
A cancellation signal is generated so that the point where the amount of noise suppression is maximum is located closer to the user than the installation position of the error microphone;
a sound pressure estimation step of estimating a sound pickup signal that would be picked up if a virtual microphone were installed at a position closer to the observation point than the error microphone, and obtaining an estimated sound pickup signal x(v);
a suppression signal generation step of generating a cancellation signal for suppressing noise at an installation position of the virtual microphone by using a pickup signal x(r) that picks up the noise to be suppressed and the estimated pickup signal x(v),
The virtual microphone is a microphone that is virtually installed without actually being installed,
The sound pressure estimation step is to estimate a sound pressure from the sound signal x(e) of the error microphone based on distance attenuation and phase delay from the positional relationship between the error microphone and the virtual microphone.
x(v)=[^G p1 ^G p2 ]
^G pn =w n x(e) (n=1,2)
The estimated sound pickup signal x(v) is obtained by the above. The transfer characteristic from the assumed noise source position to the error microphone position is Ge , and the transfer characteristic from the assumed noise source position to the observation point position is Gpn . When only the gain attenuation is taken into consideration, wn is
w n =|G pn |/|G e |
If we only consider the phase shift, wn is
w n =exp((arg G pn -arg G e )j)
That is,
Generation method.
アクティブノイズコントロールに用いるキャンセル信号を生成する生成方法であって、1. A method for generating a cancellation signal for use in active noise control, comprising:
騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成し、A cancellation signal is generated so that the point where the amount of noise suppression is maximum is located closer to the user than the installation position of the error microphone;
前記エラーマイクよりも観測点に近い位置に仮想マイクを設置した場合に収音される収音信号を推定し、推定収音信号x(v)を得る音圧推定ステップと、a sound pressure estimation step of estimating a sound pickup signal that would be picked up if a virtual microphone were installed at a position closer to the observation point than the error microphone, and obtaining an estimated sound pickup signal x(v);
抑圧対象の騒音を収音した収音信号x(r)と前記推定収音信号x(v)とを用いて、前記仮想マイクの設置位置における騒音を抑圧するためのキャンセル信号を生成する抑圧信号生成ステップとを含み、a suppression signal generation step of generating a cancellation signal for suppressing noise at an installation position of the virtual microphone by using a pickup signal x(r) that picks up the noise to be suppressed and the estimated pickup signal x(v),
前記仮想マイクは、実際には設置されずに仮想的に設置されるマイクであり、The virtual microphone is a microphone that is virtually installed without actually being installed,
前記音圧推定ステップは、ユーザの頭部近傍に非等間隔に配置した複数の前記エラーマイクの収音信号x(e)から、最小二乗法で推定された球面調和関数展開係数を利用して、前記推定収音信号x(v)を得る、The sound pressure estimation step obtains the estimated sound pickup signal x(v) from the sound pickup signals x(e) of the plurality of error microphones arranged at non-equidistant intervals near the user's head by using spherical harmonic function expansion coefficients estimated by a least squares method.
生成方法。 Generation method.
請求項1から請求項の何れかの生成装置としてコンピュータを機能させるためのプログラム。 A program for causing a computer to function as the generating device according to any one of claims 1 to 3 .
JP2022025303A 2022-02-22 2022-02-22 GENERATION DEVICE, GENERATION METHOD, AND PROGRAM Active JP7675989B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2022025303A JP7675989B2 (en) 2022-02-22 2022-02-22 GENERATION DEVICE, GENERATION METHOD, AND PROGRAM

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2022025303A JP7675989B2 (en) 2022-02-22 2022-02-22 GENERATION DEVICE, GENERATION METHOD, AND PROGRAM

Publications (2)

Publication Number Publication Date
JP2023121933A JP2023121933A (en) 2023-09-01
JP7675989B2 true JP7675989B2 (en) 2025-05-14

Family

ID=87799150

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022025303A Active JP7675989B2 (en) 2022-02-22 2022-02-22 GENERATION DEVICE, GENERATION METHOD, AND PROGRAM

Country Status (1)

Country Link
JP (1) JP7675989B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2024117876A (en) * 2023-02-20 2024-08-30 日本電信電話株式会社 GENERATION DEVICE, GENERATION METHOD, AND PROGRAM

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008164670A (en) 2006-12-27 2008-07-17 Sony Corp Audio output device, audio output method, noise reduction device, noise reduction method, noise reduction processing program, noise reduction audio output device, and noise reduction audio output method
JP2012168283A (en) 2011-02-10 2012-09-06 Tokai Rubber Ind Ltd Active vibration and noise dampener
WO2021100461A1 (en) 2019-11-18 2021-05-27 ソニーグループ株式会社 Signal processing device, method, and program
WO2022031279A1 (en) 2020-08-05 2022-02-10 Harman International Industries, Incorporated Occupancy based active noise cancellation systems

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7466998B2 (en) * 2020-07-03 2024-04-15 アルプスアルパイン株式会社 Active Noise Control System

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008164670A (en) 2006-12-27 2008-07-17 Sony Corp Audio output device, audio output method, noise reduction device, noise reduction method, noise reduction processing program, noise reduction audio output device, and noise reduction audio output method
JP2012168283A (en) 2011-02-10 2012-09-06 Tokai Rubber Ind Ltd Active vibration and noise dampener
WO2021100461A1 (en) 2019-11-18 2021-05-27 ソニーグループ株式会社 Signal processing device, method, and program
WO2022031279A1 (en) 2020-08-05 2022-02-10 Harman International Industries, Incorporated Occupancy based active noise cancellation systems

Also Published As

Publication number Publication date
JP2023121933A (en) 2023-09-01

Similar Documents

Publication Publication Date Title
JP7028238B2 (en) Signal processing equipment and methods, as well as programs
Zhang et al. A robust online secondary path modeling method with auxiliary noise power scheduling strategy and norm constraint manipulation
US10403259B2 (en) Multi-microphone feedforward active noise cancellation
Iwai et al. Multichannel feedforward active noise control system combined with noise source separation by microphone arrays
CN113421540A (en) Active noise reduction method, active noise reduction device and semi-in-ear active noise reduction earphone
JP7675989B2 (en) GENERATION DEVICE, GENERATION METHOD, AND PROGRAM
WO2021100461A1 (en) Signal processing device, method, and program
Wu et al. Active noise control pillow based on the combination of the fixed and adaptive feedback structures
JP7700882B2 (en) Noise suppression device, noise suppression system, noise suppression method, and program
CN120708588B (en) Noise reduction method, apparatus, computer device, readable storage medium, and program product
JP7782706B2 (en) Active noise control method, active noise control device, and program
JP7709142B2 (en) GENERATION DEVICE, GENERATION METHOD, AND PROGRAM
JP7847779B2 (en) Generation apparatus, generation method, and program
CN113132848B (en) Filter design method, device and in-ear active noise reduction headset
JP2024117876A (en) GENERATION DEVICE, GENERATION METHOD, AND PROGRAM
JP2023121744A (en) Stereophonic sound player
US20210375256A1 (en) Signal processing device and method, and program
JP7447993B2 (en) Elimination filter coefficient generation method, erasure filter coefficient generation device, program
CN113488016B (en) Coefficient determination method and device
JP2026029070A (en) Sound pressure estimation device, sound pressure estimation method, and program
JP2024158070A (en) Spatial active noise control device, spatial active noise control method, control filter calculation device, control filter calculation method, and program
Wang et al. Multi-Channel ANC with Adaptive Kernel Assisted on-Line Secondary Path Modeling
WO2023139753A1 (en) Noise suppression device, noise suppression method, and program
WO2025109672A1 (en) Cancel signal generation device
Wu et al. Effects of Reverberation on Active Noise Control Headrest Performance

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20220222

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240207

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240917

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20241112

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250106

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20250415

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250416

R150 Certificate of patent or registration of utility model

Ref document number: 7675989

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350