JP7675989B2 - GENERATION DEVICE, GENERATION METHOD, AND PROGRAM - Google Patents
GENERATION DEVICE, GENERATION METHOD, AND PROGRAM Download PDFInfo
- 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
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
しかしながら、実際の利用に際しては、エラーマイク93をユーザの耳元近くに設置できない場合もあり、エラーマイク93の設置位置とユーザの耳元との距離が大きくなると、前述の通り、エラーマイク93の設置位置で最も効率よく騒音を抑圧し、ユーザの耳元では騒音の消し残しが大きくなり、抑圧性能が低下し、ユーザが騒音抑圧の恩恵を十分に得られない場合がある。例えば、騒音源から耳元までの距離が100mmであり、エラーマイク93をユーザの耳元(0mm)に設置した場合の抑圧性能は-∞dBであり、エラーマイク93を騒音源と耳元との中間地点に設置した場合の抑圧性能は-7.38dBであることをシミュレーションにて確認した。図2は、従来技術の抑圧可能領域(スイートスポット)S1と所望のスイートスポットS2との違いを説明するための図である。
However, in actual use, there are cases where the
本発明は、実際のエラーマイクの配置位置で収音した収音信号から、ユーザの耳元で収音した場合に得られる収音信号を推定し、キャンセル信号を能動的に制御するために、実際のエラーマイクの配置位置で収音した収音信号の代わりに推定により得られた収音信号を使用することで、実際のエラーマイクから離れたユーザの耳元でも高い抑圧性能を実現する生成装置、生成方法、そのプログラムを提供することを目的とする。 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.
以下、本発明の実施形態について、説明する。なお、以下の説明に用いる図面では、同じ機能を持つ構成部や同じ処理を行うステップには同一の符号を記し、重複説明を省略する。以下の説明において、テキスト中で使用する記号「^」「-」等は、本来直後の文字の真上に記載されるべきものであるが、テキスト記法の制限により、当該文字の直前に記載する。式中においてはこれらの記号は本来の位置に記述している。また、ベクトルや行列の各要素単位で行われる処理は、特に断りが無い限り、そのベクトルやその行列の全ての要素に対して適用されるものとする。 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
抑圧信号生成装置は、参照マイク91の収音信号x(r)と、エラーマイク93の収音信号x(e)とを入力とし、騒音の抑圧量が最大となる地点が、エラーマイク93の設置位置よりもユーザ側に位置するようにキャンセル信号(以下、「抑圧信号」ともいう)yを生成して、キャンセルスピーカ92に出力する。
The suppression signal generating device receives the pickup signal x(r) of the
抑圧信号生成装置は、例えば、中央演算処理装置(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で収音した抑圧対象の音を、以下「騒音」と記載する。
<
The
<キャンセルスピーカ92>
キャンセルスピーカ92は、キャンセル信号yを入力とし、キャンセル信号yを再生する(S92)。キャンセルスピーカ92から再生される再生音と抑圧対象の騒音とが完全な逆位相となる場合、再生音と抑圧対象の騒音とが重なる、すなわち、音波同士が重畳する、と波が打ち消し合うため、騒音が抑圧される。
<
The
<エラーマイク93>
エラーマイク93は、騒音の消し残しを含む、キャンセルスピーカ92から再生される再生音で抑圧されなかった音を収音し(S93)、収音信号x(e)を出力する。エラーマイク93は、観測点(例えば、ユーザの耳元)よりも騒音源に近い位置に配置される。例えば、エラーマイク93は、図5のように耳元よりも騒音源に0.05m近い位置に配置される。
<Error Microphone 93>
The
<音圧推定部120>
音圧推定部120は、エラーマイク93の出力信号(収音信号)x(e)を入力とし、エラーマイク93よりも観測点に近い位置にマイク130を設置した場合に収音されると推定される信号である、推定収音信号x(v)を算出し、出力する。すなわち、音圧推定部120は、キャンセルスピーカ92から再生される再生音で抑圧されなかった音がマイク130の設置位置で収音される場合に得られる収音信号を推定し(S120)、推定した収音信号を推定収音信号x(v)として出力する。以下、推定収音信号x(v)の推定方法を3つ例示する。ここで、マイク130は実際には設置せず仮想的に設置されるものであり、以下仮想マイク130と記載する。
<Sound
The sound
(推定方法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
本推定方法では、騒音源の位置を仮定し、騒音源からエラーマイク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
^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(θ1,φ1),p(θ2,φ2),…,p(θL,φL)を得る。例えば、L個のエラーマイク93の収音信号x(e)=[p(θ1,φ1),p(θ2,φ2),…,p(θL,φL)]とする。
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
音圧推定部120は、次式により、球面調和関数Ym
n(・)に対する半径re上の音場係数Pnm(re)を求める。
音圧推定部120は、求めた音場係数Pnm(re)を用いて、次式により、半径r上の音場係数Pnm(r)を求める。
ただし、kは波数であり、aは音を反射する剛球の半径とし、jnはn次の球面ベッセル関数であり、jn'はjnの微分であり、hn
(2)は第二種球ハンケル関数であり、hn'(2)はhn
(2)の微分である。
音圧推定部120は、再合成により、観測点(r,θ,φ)における音圧の推定値^p(r,θ,φ)を得る。
なお、推定収音信号x(v)=^p(r,θ,φ)とする。
The sound
The sound
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
The estimated picked-up signal is x(v)=^p(r,θ,φ).
以下、式(5)の導出について説明する。 The derivation of equation (5) is explained below.
騒音源を点音源とし、半径aの剛球での反射を考慮したとき、点(r,θ,φ)における音圧は、
である。なお、Bnmは騒音源の座標と信号で定まる係数である。球面調和関数展開は、
であり、音場係数Pnm(r)、音場係数Pnm(re)は、次式で表される。
式(10)をBnm=…の形に変形して(9)に代入すると、
となる。
なお、球面調和関数展開における最大次数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
Note that B nm is a coefficient determined by the coordinates of the noise source and the signal. The spherical harmonic expansion is
and the sound field coefficients P nm (r) and P nm (r e ) are expressed by the following equations.
Transforming equation (10) into the form B nm = ... and substituting it into (9), we get
It becomes.
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(re,θ1,φ1),p(re,θ2,φ2),…,p(re,θL,φL)を得る。エラーマイク93の収音信号x(e)=-p=[p(re,θ1,φ1),p(re,θ2,φ2),…,p(re,θL,φL)]Tとする。-pは以下のように表される。
なお、ωi=(θi,φi)とし、i=1,2,…,Lとする。
Observation values of sound pressure p(r e ,θ 1 ,φ 1 ),p(r e ,θ 2 ,φ 2 ),...,p(r e ,θ L ,φ L ) are obtained from L error microphones on a radius of r e . The sound pickup signal of the
Here, ω i = (θ i , φ i ), where i = 1, 2, ..., L.
音圧推定部120は、絶対値の二乗誤差が最小となる-P(re)を解として求める。
推定したい点が半径r上にあるとしたとき、展開係数Pnm(r)および推定値の計算は推定方法2と同様である。つまり、
音圧推定部120は、次式により、球面調和関数Ym
n(・)に対する半径re上の音場係数Pnm(re)を求める。
音圧推定部120は、求めた音場係数Pnm(re)を用いて、次式により、半径r上の音場係数Pnm(r)を求める。
音圧推定部120は、再合成により、観測点(r,θ,φ)における音圧の推定値^p(r,θ,φ)を得る。
なお、推定収音信号x(v)=^p(r,θ,φ)とする。
The sound
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
The sound
The sound
The estimated picked-up signal is x(v)=^p(r,θ,φ).
<抑圧信号生成部110>
抑圧信号生成部110は、収音信号x(r)と推定収音信号x(v)とを入力とし、仮想マイク130の設置位置における騒音を抑圧するためのキャンセル信号yを生成し(S110)、出力する。
<Suppression
The
キャンセル信号の生成方法としては、従来技術を用いることができる。例えば、非特許文献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
<効果>
以上の構成により、ユーザの耳元にエラーマイクを配置できない場合に従来よりも高い抑圧性能を実現することができる。図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
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。 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.
騒音の抑圧量が最大となる地点が、エラーマイクの設置位置よりもユーザ側に位置するようにキャンセル信号を生成し、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.
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)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2024117876A (en) * | 2023-02-20 | 2024-08-30 | 日本電信電話株式会社 | GENERATION DEVICE, GENERATION METHOD, AND PROGRAM |
Citations (4)
| 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)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7466998B2 (en) * | 2020-07-03 | 2024-04-15 | アルプスアルパイン株式会社 | Active Noise Control System |
-
2022
- 2022-02-22 JP JP2022025303A patent/JP7675989B2/en active Active
Patent Citations (4)
| 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 |