JP5152800B2 - Noise suppression evaluation apparatus and program - Google Patents

Noise suppression evaluation apparatus and program Download PDF

Info

Publication number
JP5152800B2
JP5152800B2 JP2008179444A JP2008179444A JP5152800B2 JP 5152800 B2 JP5152800 B2 JP 5152800B2 JP 2008179444 A JP2008179444 A JP 2008179444A JP 2008179444 A JP2008179444 A JP 2008179444A JP 5152800 B2 JP5152800 B2 JP 5152800B2
Authority
JP
Japan
Prior art keywords
noise
kurtosis
suppression
index value
acoustic signal
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.)
Expired - Fee Related
Application number
JP2008179444A
Other languages
Japanese (ja)
Other versions
JP2010020013A (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.)
Yamaha Corp
Original Assignee
Yamaha Corp
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 Yamaha Corp filed Critical Yamaha Corp
Priority to JP2008179444A priority Critical patent/JP5152800B2/en
Priority to EP09164896.4A priority patent/EP2144233A3/en
Priority to US12/499,734 priority patent/US20100008520A1/en
Publication of JP2010020013A publication Critical patent/JP2010020013A/en
Application granted granted Critical
Publication of JP5152800B2 publication Critical patent/JP5152800B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、音響信号の雑音成分を抑圧する処理を評価する技術に関する。   The present invention relates to a technique for evaluating a process for suppressing a noise component of an acoustic signal.

信号成分（目的音の成分）と雑音成分とが重畳された音響信号について雑音成分を抑圧する技術が従来から提案されている。例えば非特許文献１および非特許文献２には、周波数領域で音響信号から雑音成分を抑圧するスペクトラルサブトラクション（Spectral Subtraction：SS）技術が開示されている。   Conventionally, a technique for suppressing a noise component of an acoustic signal in which a signal component (target sound component) and a noise component are superimposed has been proposed. For example, Non-Patent Document 1 and Non-Patent Document 2 disclose a spectral subtraction (SS) technique for suppressing a noise component from an acoustic signal in a frequency domain.

しかし、非特許文献１および非特許文献２のように周波数領域で音響信号から雑音成分を抑圧する方法においては、雑音成分の抑圧後に時間軸上および周波数軸上に分散的に雑音成分が残存し、耳障りなミュージカルノイズ（バーディノイズ）として受聴者に知覚されるという問題がある。そこで、非特許文献３には、雑音成分の抑圧に起因したミュージカルノイズを発生後に除去する技術が開示されている。
Steven F. Boll, "Suppression of Acoustic Noise in Speech Using Spectral Subtraction", IEEE TRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING, Vol. ASSP-27, No. 2, April 1979 Yariv Ephraim, David Malah, "Speech Enhancement Using a Minimum Mean-Square Error Short-Time Spectral Amplitude Estimator", IEEE TRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING, Vol. ASSP-32, No. 6, December 1984 阿部友実，松本光春，橋本周司，"時間−周波数Ｍ変換によるミュージカルノイズ除去"，日本音響学会講演論文集，3-6-9，p.727−p.730，2008年３月 However, in the method of suppressing the noise component from the acoustic signal in the frequency domain as in Non-Patent Document 1 and Non-Patent Document 2, the noise component remains dispersively on the time axis and the frequency axis after the noise component is suppressed. There is a problem that it is perceived by the listener as harsh musical noise (birdie noise). Therefore, Non-Patent Document 3 discloses a technique for removing musical noise resulting from noise component suppression after generation.
Steven F. Boll, "Suppression of Acoustic Noise in Speech Using Spectral Subtraction", IEEE TRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING, Vol. ASSP-27, No. 2, April 1979 Yariv Ephraim, David Malah, "Speech Enhancement Using a Minimum Mean-Square Error Short-Time Spectral Amplitude Estimator", IEEE TRANSACTIONS ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING, Vol. ASSP-32, No. 6, December 1984 Tomomi Abe, Mitsuharu Matsumoto, Shuji Hashimoto, "Musical noise removal by time-frequency M conversion", Proceedings of the Acoustical Society of Japan, 3-6-9, p.727-p.730, March 2008

ところで、雑音の抑圧後におけるミュージカルノイズの発生の程度を定量的に評価できれば、例えばミュージカルノイズを適切に除去できるように雑音成分の抑圧の度合を可変に制御するといった構成も実現され得る。しかし、ミュージカルノイズの発生の程度を定量的に評価する方法については非特許文献１や非特許文献２には何ら開示されていない。非特許文献３は、発生後のミュージカルノイズの除去を開示するに過ぎず、非特許文献１や非特許文献２と同様にミュージカルノイズの定量的な評価については何ら開示しない。以上の事情を考慮して、本発明は、ミュージカルノイズの発生の程度の定量的な指標を提供することを目的とする。   By the way, if the degree of occurrence of musical noise after noise suppression can be quantitatively evaluated, for example, a configuration of variably controlling the degree of noise component suppression so that musical noise can be appropriately removed can be realized. However, NPL 1 and NPL 2 do not disclose any method for quantitatively evaluating the degree of occurrence of musical noise. Non-Patent Document 3 merely discloses the removal of musical noise after generation, and does not disclose any quantitative evaluation of musical noise as in Non-Patent Document 1 and Non-Patent Document 2. In view of the above circumstances, an object of the present invention is to provide a quantitative indicator of the degree of occurrence of musical noise.

以上の課題を解決するために、本発明に係る雑音抑圧評価装置は、雑音成分の抑圧前または抑圧後（すなわち抑圧前および抑圧後の少なくとも一方）の音響信号の強度の度数分布における尖度に応じて変化する雑音指標値を、周波数領域における雑音成分の抑圧後にミュージカルノイズが発生する程度の指標として算定する指標算定手段を具備する。上の構成においては、雑音成分の抑圧前または抑圧後の音響信号の強度の度数分布における尖度に応じて変化する雑音指標値を、雑音成分の抑圧後の音響信号におけるミュージカルノイズの発生の度合の定量的な指標として利用できる。   In order to solve the above-described problems, the noise suppression evaluation apparatus according to the present invention has a kurtosis in the frequency distribution of the intensity of an acoustic signal before or after suppression of noise components (that is, at least one before and after suppression). There is provided an index calculating means for calculating a noise index value that changes accordingly as an index of the degree to which musical noise is generated after suppression of noise components in the frequency domain. In the above configuration, the noise index value that changes according to the kurtosis in the frequency distribution of the intensity of the acoustic signal before or after suppression of the noise component is set to the degree of occurrence of musical noise in the acoustic signal after suppression of the noise component. It can be used as a quantitative indicator.

本発明の好適な態様において、指標算定手段は、雑音成分の抑圧の度合を示す抑圧係数と尖度に応じた尖度指標値との関係を特定する相関特定手段と、相関特定手段が特定した関係において尖度指標値が所定値に接近するときの抑圧係数を雑音指標値として決定する指標決定手段とを含む。以上の態様においては、ミュージカルノイズの発生を所定の程度に制御するために必要な雑音抑圧の度合という観点から、雑音成分の抑圧後の音響信号におけるミュージカルノイズの発生の程度が表現される。なお、尖度指標値の目標となる所定値は、固定値および可変値の何れでもよい。   In a preferred aspect of the present invention, the index calculating means specifies the correlation specifying means for specifying the relationship between the suppression coefficient indicating the degree of suppression of the noise component and the kurtosis index value corresponding to the kurtosis, and the correlation specifying means specifies And an index determination unit that determines a suppression coefficient when the kurtosis index value approaches a predetermined value as a noise index value. In the above aspect, from the viewpoint of the degree of noise suppression necessary for controlling the generation of musical noise to a predetermined level, the degree of occurrence of musical noise in the acoustic signal after suppression of the noise component is expressed. Note that the predetermined value that is the target of the kurtosis index value may be either a fixed value or a variable value.

本発明の好適な態様において、指標算定手段は、雑音成分の抑圧前の音響信号の強度の度数分布における第１尖度を算定する第１尖度算定手段と、雑音成分の抑圧後の音響信号の強度の度数分布における第２尖度を算定する第２尖度算定手段と、第１尖度と第２尖度とから雑音指標値を算定する算定手段とを含む。以上の態様においては、第１尖度および第２尖度の双方に応じて雑音指標値が算定される（さらには雑音抑圧手段による抑圧の度合が制御される）から、例えば第１尖度および第２尖度の一方のみから雑音指標値を算定する構成と比較して、ミュージカルノイズの発生の度合を正確に示す雑音指標値が提供される。   In a preferred aspect of the present invention, the index calculating means includes first kurtosis calculating means for calculating the first kurtosis in the frequency distribution of the intensity of the acoustic signal before suppression of the noise component, and the acoustic signal after suppression of the noise component. The second kurtosis calculating means for calculating the second kurtosis in the intensity distribution of the intensity of the light and the calculating means for calculating the noise index value from the first kurtosis and the second kurtosis. In the above aspect, since the noise index value is calculated according to both the first kurtosis and the second kurtosis (and the degree of suppression by the noise suppression unit is controlled), for example, the first kurtosis and Compared with the configuration in which the noise index value is calculated from only one of the second kurtosis, a noise index value that accurately indicates the degree of occurrence of musical noise is provided.

以上の各態様においては、雑音成分の抑圧前の音響信号の第１尖度が小さいほど、雑音指標値が示すミュージカルノイズの発生の程度が増加するように、指標算定手段が雑音指標値を算定する構成（第２尖度の使用が必須でない構成）や、雑音成分の抑圧後の音響信号の第２尖度が小さいほど、雑音指標値が示すミュージカルノイズの発生の程度が減少するように、指標算定手段が雑音指標値を算定する構成（第１尖度の使用が必須でない構成）が好適である。なお、第２尖度は、雑音抑圧手段による実際の処理後の音響信号から算定されるほか、雑音抑圧手段の動作を模擬する（例えば式(16)の演算を実行する）ことで抑圧前の音響信号から算定（推定）される。   In each of the above aspects, the index calculation means calculates the noise index value so that the degree of occurrence of musical noise indicated by the noise index value increases as the first kurtosis of the acoustic signal before suppression of the noise component decreases. In order to reduce the degree of occurrence of musical noise indicated by the noise index value, the smaller the second kurtosis of the sound signal after suppressing the noise component (the configuration in which the use of the second kurtosis is not essential) or the noise component is suppressed, A configuration in which the index calculation means calculates the noise index value (a configuration in which the use of the first kurtosis is not essential) is preferable. Note that the second kurtosis is calculated from the acoustic signal actually processed by the noise suppression means, and also simulates the operation of the noise suppression means (for example, by executing the calculation of equation (16)). Calculated (estimated) from the acoustic signal.

また、雑音成分の抑圧の前後における尖度の変化度がミュージカルノイズの発生の度合に最も顕著に反映されるという傾向を考慮すると、雑音成分の抑圧前の音響信号の第１尖度に対する抑圧後の音響信号の第２尖度の相対比が大きいほど、雑音指標値が示すミュージカルノイズの発生の程度が増加するように、指標算定手段が第１尖度および第２尖度に応じた雑音指標値を算定する構成が好適である。特に、第１尖度に対する第２尖度の相対比の対数値がミュージカルノイズの発生の度合に対して高い相関を示すという傾向を考慮すると、第１尖度に対する第２尖度の相対比の対数値が大きいほど、雑音指標値が示すミュージカルノイズの発生の程度が増加するように、指標算定手段が当該対数値に応じた雑音指標値を算定する構成が好適である。   Further, in consideration of the tendency that the degree of change in kurtosis before and after noise component suppression is most significantly reflected in the degree of occurrence of musical noise, after suppression of the first kurtosis of the acoustic signal before noise component suppression. The noise index corresponding to the first kurtosis and the second kurtosis is calculated by the index calculation means so that the degree of occurrence of musical noise indicated by the noise index value increases as the relative ratio of the second kurtosis of the acoustic signal increases. A configuration for calculating the value is preferable. In particular, considering the tendency that the logarithm of the relative ratio of the second kurtosis to the first kurtosis shows a high correlation with the degree of occurrence of musical noise, the relative ratio of the second kurtosis to the first kurtosis A configuration in which the index calculation means calculates the noise index value corresponding to the logarithmic value is preferable so that the degree of occurrence of musical noise indicated by the noise index value increases as the logarithmic value increases.

なお、指標算定手段が算定した雑音指標値は、雑音抑圧装置における雑音成分の抑圧に利用される。すなわち、雑音抑圧装置は、以上の各態様に係る雑音抑圧評価装置（指標算定手段）と、音響信号の雑音成分を周波数領域で抑圧する雑音抑圧手段と、雑音抑圧手段による雑音成分の抑圧の度合を雑音指標値に応じて可変に制御する抑圧制御手段とを具備する。以上の構成においては、雑音抑圧手段による雑音成分の抑圧の度合が雑音指標値に応じて可変に制御されるから、雑音成分の抑圧の度合が固定された従来の技術と比較して、ミュージカルノイズの発生を有効に制御（典型的には抑制）しながら雑音成分を抑圧することが可能である。例えば、雑音指標値が示すミュージカルノイズの発生の度合が低いほど雑音成分の抑圧の度合が増加するように、雑音抑圧手段による抑圧の度合を抑圧制御手段が雑音指標値に応じて制御する構成においては、ミュージカルノイズの発生を有効に抑制しながら雑音成分を抑圧することが可能である。   Note that the noise index value calculated by the index calculation means is used for noise component suppression in the noise suppression device. That is, the noise suppression device includes a noise suppression evaluation device (index calculation unit) according to each of the above aspects, a noise suppression unit that suppresses a noise component of an acoustic signal in a frequency domain, and a degree of noise component suppression by the noise suppression unit. Suppression control means for variably controlling according to the noise index value. In the above configuration, since the degree of suppression of the noise component by the noise suppression unit is variably controlled according to the noise index value, compared with the conventional technique in which the degree of suppression of the noise component is fixed, It is possible to suppress the noise component while effectively controlling (typically suppressing) the occurrence of. For example, in a configuration in which the suppression control unit controls the degree of suppression by the noise suppression unit according to the noise index value so that the degree of suppression of the noise component increases as the degree of occurrence of musical noise indicated by the noise index value decreases. Can suppress noise components while effectively suppressing the generation of musical noise.

また、以上の各態様に係る雑音抑圧評価装置は、雑音の抑圧に専用されるＤＳＰ（Digital Signal Processor）などのハードウェア（電子回路）によって実現されるほか、ＣＰＵ（Central Processing Unit）などの汎用の演算処理装置とプログラムとの協働によっても実現される。本発明に係るプログラムは、雑音成分の抑圧前または抑圧後の音響信号の強度の度数分布における尖度に応じて変化する雑音指標値を、周波数領域における雑音成分の抑圧後にミュージカルノイズが発生する程度の指標として算定する指標算定処理をコンピュータに実行させる。以上のプログラムによれば、本発明の各態様に係る雑音抑圧評価装置と同様の作用および効果が奏される。本発明のプログラムは、コンピュータが読取可能な記録媒体に格納された形態で利用者に提供されてコンピュータにインストールされるほか、通信網を介した配信の形態でサーバ装置から提供されてコンピュータにインストールされる。   The noise suppression evaluation apparatus according to each of the above aspects is realized by hardware (electronic circuit) such as a DSP (Digital Signal Processor) dedicated to noise suppression, and a general purpose such as a CPU (Central Processing Unit). This is also realized by cooperation between the arithmetic processing unit and the program. The program according to the present invention uses a noise index value that changes according to the kurtosis in the frequency distribution of the intensity of the acoustic signal before or after suppressing the noise component, to the extent that musical noise is generated after suppressing the noise component in the frequency domain. Causes the computer to execute an index calculation process for calculating the index of According to the above program, operations and effects similar to those of the noise suppression evaluation apparatus according to each aspect of the present invention are exhibited. The program of the present invention is provided to a user in a form stored in a computer-readable recording medium and installed in the computer, or provided from a server device in a form of distribution via a communication network and installed in the computer. Is done.

＜Ａ：第１実施形態＞
図１は、本発明の第１実施形態に係る雑音抑圧装置のブロック図である。雑音抑圧装置１００には、音響の波形を表す時間領域の音響信号ＶINが供給される。音響信号ＶINの供給元（図示略）は、例えば、周囲の音響に応じた音響信号ＶINを生成する収音機器や、記録媒体から音響信号ＶINを取得して出力する再生装置である。音響信号ＶINには信号成分ｓと雑音成分ｎとが混在する（ＶIN＝ｓ＋ｎ）。雑音抑圧装置１００は、音響信号ＶINの雑音成分ｎを抑圧することで音響信号ＶOUT（理想的にはＶOUT＝ｓ）を生成して出力する。音響信号ＶOUTは、例えばスピーカ装置やヘッドホンなどの放音装置（図示略）に供給されて音波として再生される。 <A: First Embodiment>
FIG. 1 is a block diagram of a noise suppression apparatus according to the first embodiment of the present invention. The noise suppression apparatus 100 is supplied with a time-domain acoustic signal VIN representing an acoustic waveform. The supply source (not shown) of the acoustic signal VIN is, for example, a sound collecting device that generates an acoustic signal VIN corresponding to surrounding sounds, or a playback device that acquires and outputs the acoustic signal VIN from a recording medium. The acoustic signal VIN includes a signal component s and a noise component n (VIN = s + n). The noise suppression device 100 generates and outputs an acoustic signal VOUT (ideally VOUT = s) by suppressing the noise component n of the acoustic signal VIN. The acoustic signal VOUT is supplied to a sound emitting device (not shown) such as a speaker device or headphones and reproduced as a sound wave.

雑音抑圧装置１００は、演算処理装置１２と記憶装置１４とを含むコンピュータシステムで実現される。記憶装置１４は、音響信号ＶINから音響信号ＶOUTを生成するためのプログラムや各種のデータを記憶する。半導体記憶装置や磁気記憶装置などの公知の記録媒体が記憶装置１４として任意に採用される。   The noise suppression device 100 is realized by a computer system including an arithmetic processing device 12 and a storage device 14. The storage device 14 stores a program and various data for generating the acoustic signal VOUT from the acoustic signal VIN. A known recording medium such as a semiconductor storage device or a magnetic storage device is arbitrarily adopted as the storage device 14.

演算処理装置１２は、記憶装置１４に格納されたプログラムを実行することで複数の要素（周波数分析部２２，雑音推定部２４，雑音抑圧部２６，波形合成部２８，指標算定部３２，ＳＮ比算定部３４，抑圧制御部３６）として機能する。なお、音響信号ＶINの処理に専用される電子回路（ＤＳＰ）が演算処理装置１２の各要素を実現する構成や、演算処理装置１２の各要素を複数の集積回路に分散的に搭載した構成も採用される。   The arithmetic processing device 12 executes a program stored in the storage device 14 to execute a plurality of elements (frequency analysis unit 22, noise estimation unit 24, noise suppression unit 26, waveform synthesis unit 28, index calculation unit 32, SN ratio). The calculation unit 34 and the suppression control unit 36) function. In addition, a configuration in which an electronic circuit (DSP) dedicated to processing of the acoustic signal VIN realizes each element of the arithmetic processing device 12 or a configuration in which each element of the arithmetic processing device 12 is mounted on a plurality of integrated circuits in a distributed manner. Adopted.

図１の周波数分析部２２は、図２に示すように、音響信号ＶINを時間軸上で区分した複数のフレームＦRの各々についてフーリエ変換を実行することで各フレームＦRの周波数スペクトルＸm(ejω)（図１や図２では「Ｘ」と簡略化されている）を算定する。第ｍ番目のフレームＦRの周波数スペクトルＸm(ejω)は、信号成分ｓの周波数スペクトルＳm(ejω)と雑音成分ｎの周波数スペクトルＮm(ejω)との加算に相当する（式(1)）。
Ｘm(ejω)＝Ｓm(ejω)＋Ｎm(ejω) ……(1) As shown in FIG. 2, the frequency analysis unit 22 in FIG. 1 performs a Fourier transform on each of a plurality of frames FR obtained by dividing the acoustic signal VIN on the time axis, so that the frequency spectrum Xm (e jω of each frame FR is obtained. ) (Simplified as “X” in FIGS. 1 and 2). Frequency spectrum Xm of the m-th frame FR (e j [omega]) is equivalent to the addition of the frequency spectrum Nm (e j [omega]) of the frequency spectrum Sm (e j [omega]) and the noise component n of the signal component s (Equation (1) ).
Xm (e jω ) = Sm (e jω ) + Nm (e jω ) (1)

図１の雑音推定部２４は、音響信号ＶINに重畳された雑音成分ｎの周波数スペクトル（以下「推定雑音スペクトル」という）ψm(ejω)を音響信号ＶINの複数のフレームＦRの各々について推定する。図１に示すように、雑音推定部２４は判定部２４２と推定部２４４とで構成される。判定部２４２は、各フレームＦRにおける信号成分ｓの有無を周波数スペクトルＸm(ejω)に応じて判定する。信号成分ｓの有無の判定には公知の技術が任意に採用される。 1 estimates a frequency spectrum (hereinafter referred to as “estimated noise spectrum”) ψm (e jω ) of a noise component n superimposed on the acoustic signal VIN for each of a plurality of frames FR of the acoustic signal VIN. . As shown in FIG. 1, the noise estimation unit 24 includes a determination unit 242 and an estimation unit 244. The determination unit 242 determines the presence / absence of the signal component s in each frame FR according to the frequency spectrum Xm (e jω ). A known technique is arbitrarily employed to determine the presence or absence of the signal component s.

推定部２４４は、判定部２４２による判定の結果を利用して推定雑音スペクトルψm(ejω)を算定する。さらに詳述すると、推定部２４４は、信号成分ｓを含まない（または少ない）と判定部２４２が判定した区間（以下「雑音区間」という）内の各フレームＦRについて周波数スペクトルＸm(ejω)を平均することで推定雑音スペクトルψm(ejω)を算定する。雑音区間では周波数スペクトルＸm(ejω)が周波数スペクトルＮm(ejω)に略一致するから、推定雑音スペクトルψm(ejω)は以下の式(2)で周波数スペクトルＮm(ejω)から算定される。式(2)の演算子Ｅは、期待値（平均値）の算定を意味する。
ψm(ejω)＝Ｅ｛|Ｎm(ejω)|2｝ ……(2) The estimation unit 244 calculates the estimated noise spectrum ψm (e jω ) using the result of determination by the determination unit 242. More specifically, the estimation unit 244 calculates the frequency spectrum Xm (e jω ) for each frame FR in the section (hereinafter referred to as “noise section”) determined by the determination unit 242 as not including (or reducing) the signal component s. The estimated noise spectrum ψm (e jω ) is calculated by averaging. In the noise interval, the frequency spectrum Xm (e jω ) substantially matches the frequency spectrum Nm (e jω ), so the estimated noise spectrum ψm (e jω ) is calculated from the frequency spectrum Nm (e jω ) using the following equation (2). The The operator E in the equation (2) means calculation of an expected value (average value).
ψm (e jω ) = E {| Nm (e jω ) | 2 } (2)

また、推定部２４４は、判定部２４２が信号成分ｓを含むと判定したフレームＦRについては、直前の推定雑音スペクトルψm-1(ejω)と同じ推定雑音スペクトルψm(ejω)を設定する（ψm(ejω)＝ψm-1(ejω)）。以上のように推定雑音スペクトルψm(ejω)はフレームＦR毎に順次に更新される。なお、推定雑音スペクトルψm(ejω)の推定には公知の技術が任意に採用される。 Further, the estimation unit 244 sets the same estimated noise spectrum ψm (e jω ) as the immediately preceding estimated noise spectrum ψm-1 (e jω ) for the frame FR determined by the determination unit 242 to include the signal component s ( ψm (e jω ) = ψm-1 (e jω )). As described above, the estimated noise spectrum ψm (e jω ) is sequentially updated for each frame FR. A known technique is arbitrarily employed for estimating the estimated noise spectrum ψm (e jω ).

雑音抑圧部２６は、音響信号ＶINの雑音成分ｎ（周波数スペクトルＮm(ejω)）を周波数領域で抑圧する。さらに詳述すると、雑音抑圧部２６は、周波数分析部２２が順次に算定する周波数スペクトルＸm(ejω)から推定雑音スペクトルψm(ejω)を減算することで周波数スペクトルＹm(ejω)（図１では「Ｙ」と簡略化されている）を算定する（Spectral Subtraction）。 The noise suppression unit 26 suppresses the noise component n (frequency spectrum Nm (e jω )) of the acoustic signal VIN in the frequency domain. More specifically, the noise suppression unit 26 subtracts the estimated noise spectrum ψm (e jω ) from the frequency spectrum Xm (e jω ) sequentially calculated by the frequency analysis unit 22 to thereby obtain a frequency spectrum Ym (e jω ) (FIG. 1 is simplified as “Y” (Spectral Subtraction).

雑音抑圧部２６が算定する周波数スペクトルＹm(ejω)は、式(3)に示すように、推定雑音スペクトルψm(ejω)に応じて算定されたパワースペクトルＰmの平方根に周波数スペクトルＸm(ejω)の位相成分ｅjθx(ejω)を付加することで算定される。
Ｙm(ejω)＝（Ｐm）１/２・ｅjθx(ejω) ……(3)
式(3)のパワースペクトルＰmは以下の式(4a)および式(4b)で算定される。

Figure 0005152800

すなわち、パワースペクトルＰmのうち周波数スペクトルＸm(ejω)の強度の自乗|Ｘm(ejω)|２が推定雑音スペクトルψm(ejω)と係数αmとの乗算値（αm・ψm(ejω)）を上回る周波数帯域の成分は、式(4a)に示すように、当該乗算値（αm・ψm(ejω)）を周波数スペクトルＸm(ejω)の強度の自乗|Ｘm(ejω)|２から減算することで算定される。一方、パワースペクトルＰmのうち周波数スペクトルＸm(ejω)の強度の自乗|Ｘm(ejω)|２が推定雑音スペクトルψm(ejω)と係数αmとの乗算値（αm・ψm(ejω)）を下回る周波数帯域の成分は、式(4b)に示すように、推定雑音スペクトルψm(ejω)と係数（フロアリング係数）βmとの乗算値（βm・ψm(ejω)）に設定される。係数αmや係数βmの詳細は後述する。 The frequency spectrum Ym (e jω ) calculated by the noise suppression unit 26 is the frequency spectrum Xm (e) at the square root of the power spectrum Pm calculated according to the estimated noise spectrum ψm (e jω ), as shown in Equation (3). It is calculated by adding a phase component e jθx (ejω) of jω ).
Ym (e jω ) = (Pm) 1/2 · e jθx (ejω) (3)
The power spectrum Pm of the equation (3) is calculated by the following equations (4a) and (4b).
Figure 0005152800
That is, the square of the intensity of the frequency spectrum Xm (e jω ) | Xm (e jω ) | 2 of the power spectrum Pm is a product of the estimated noise spectrum ψm (e jω ) and the coefficient αm (αm · ψm (e jω ) ) component of the frequency band above the formula (as shown in 4a), the multiplication value (αm · ψm (e jω) square of the intensity of) the frequency spectrum Xm (e jω) | Xm ( e jω) | 2 Calculated by subtracting from On the other hand, the square of the intensity of the frequency spectrum Xm (e jω ) | Xm (e jω ) | 2 of the power spectrum Pm is the product of the estimated noise spectrum ψm (e jω ) and the coefficient αm (αm · ψm (e jω ) )) Is set to the product of the estimated noise spectrum ψm (e jω ) and the coefficient (flooring coefficient) βm (βm · ψm (e jω )) as shown in equation (4b). The Details of the coefficient αm and the coefficient βm will be described later.

図１の波形合成部２８は、雑音抑圧部２６がフレームＦR毎に算定した周波数スペクトルＹm(ejω)から時間領域の音響信号ＶOUTを合成する。さらに詳述すると、波形合成部２８は、周波数スペクトルＹm(ejω)に対する逆フーリエ変換で算定した時間領域の信号を複数のフレームＦRについて時間軸上で重複させて加算することで音響信号ＶOUTを算定する。 The waveform synthesizer 28 in FIG. 1 synthesizes the time-domain acoustic signal VOUT from the frequency spectrum Ym (e jω ) calculated by the noise suppressor 26 for each frame FR. More specifically, the waveform synthesizing unit 28 adds the acoustic signal VOUT by overlapping the time domain signals calculated by the inverse Fourier transform on the frequency spectrum Ym (e jω ) on the time axis for a plurality of frames FR. Calculate.

以上のように周波数スペクトルＸm(ejω)から推定雑音スペクトルψm(ejω)（αm・ψm(ejω)）を減算することで雑音成分ｎを抑圧した音響信号ＶOUTには、時間軸上および周波数軸上に分散的にミュージカルノイズが点在する場合がある。図１の指標算定部３２は、音響信号ＶOUTにおけるミュージカルノイズの発生の度合の定量的な指標となる雑音指標値σmをフレームＦR毎に算定する。雑音指標値σmの詳細については後述する。 As described above, the acoustic signal VOUT in which the noise component n is suppressed by subtracting the estimated noise spectrum ψm (e jω ) (αm · ψm (e jω )) from the frequency spectrum Xm (e jω ) There are cases where musical noise is scattered in a distributed manner on the frequency axis. The index calculation unit 32 in FIG. 1 calculates a noise index value σm, which is a quantitative index of the degree of occurrence of musical noise in the acoustic signal VOUT, for each frame FR. Details of the noise index value σm will be described later.

ＳＮ比算定部３４は、音響信号ＶINのＳＮ比ξmをフレームＦR毎に算定する。さらに詳述すると、ＳＮ比算定部３４は、第ｍ番目のフレームＦRの推定雑音スペクトルψm(ejω)の強度|ψm(ejω)|と直前（第(m-1)番目）のフレームＦRの周波数スペクトルＹm(ejω)の強度の自乗|Ｙm(ejω)|２との相対比を、第ｍ番目のフレームＦRのＳＮ比ξmとして算定する（ξm＝|Ｙm(ejω)|２／|ψm(ejω)|）。ただし、ＳＮ比ξmを算定する方法は任意である。また、ＳＮ比ξmの更新の周期はフレームＦRに限定されない。 The SN ratio calculation unit 34 calculates the SN ratio ξm of the acoustic signal VIN for each frame FR. More specifically, the S / N ratio calculation unit 34 determines the intensity | ψm (e jω ) | of the estimated noise spectrum ψm (e jω ) of the mth frame FR and the immediately preceding ((m−1) th) frame FR. of the square of the intensity of the frequency spectrum Ym (e jω) | Ym ( e jω) | relative ratio of 2, calculates the SN ratio Kushiemu of the m-th frame FR (ξm = | Ym (e jω) | 2 / | Ψm (e jω ) |). However, the method for calculating the SN ratio ξm is arbitrary. Further, the cycle of updating the SN ratio ξm is not limited to the frame FR.

抑圧制御部３６は、雑音抑圧部２６が音響信号ＶIN（周波数スペクトルＸm(ejω)）から雑音成分ｎ（推定雑音スペクトルψm(ejω)）を抑圧する度合を、指標算定部３２が算定した雑音指標値σmとＳＮ比算定部３４が算定したＳＮ比ξmとに応じて可変に（適応的に）制御することで、雑音抑圧部２６による処理後の音響信号ＶOUTにおけるミュージカルノイズの発生を抑制する。 In the suppression control unit 36, the index calculation unit 32 calculates the degree to which the noise suppression unit 26 suppresses the noise component n (estimated noise spectrum ψm (e jω )) from the acoustic signal VIN (frequency spectrum Xm (e jω )). By variably (adaptively) controlling the noise index value σm and the SN ratio ξm calculated by the SN ratio calculation unit 34, the generation of musical noise in the acoustic signal VOUT processed by the noise suppression unit 26 is suppressed. To do.

次に、雑音指標値σmの算定について説明する。図３の部分(A)は、音響信号ＶIN（信号成分ｓが少ないと判定部２４２が判定した雑音区間）における強度の度数分布（強度を確率変数とする確率密度関数）である。図３の部分(A)に示すように、音響信号ＶINの強度は、強度がゼロから増加するほど度数が減少するように非線形に分布する。   Next, calculation of the noise index value σm will be described. Part (A) in FIG. 3 is a frequency distribution of the intensity (a probability density function having the intensity as a random variable) in the acoustic signal VIN (the noise section in which the determination unit 242 determines that the signal component s is small). As shown in part (A) of FIG. 3, the intensity of the acoustic signal VIN is non-linearly distributed so that the frequency decreases as the intensity increases from zero.

図３の部分(B)に図示した範囲ＡSSは、雑音抑圧部２６が音響信号ＶIN（周波数スペクトルＸm(ejω)）から減算する成分（αm・ψm(ejω)）の強度に相当する。雑音抑圧部２６による抑圧で範囲ＡSS内の強度はゼロに変化するから、雑音成分ｎが抑圧された音響信号ＶOUTの強度の度数分布（図３の部分(C)）においては、ゼロに近い強度の度数が、雑音成分ｎの抑圧前の度数分布（図３の部分(A)）と比較して増加する。すなわち、強度がゼロの近傍となる範囲における度数分布の傾斜は雑音抑圧部２６による抑圧後に急峻な形状に変化する。度数分布の形状（傾斜の急峻度）の尺度として尖度（kurtosis）を導入すると、音響信号ＶOUTの第ｍ番目のフレームＦRの尖度ＫSSm（図３の部分(C)）は、雑音抑圧部２６が音響信号ＶINから雑音成分ｎを抑圧することで、抑圧前の音響信号ＶINの第ｍ番目のフレームＦRの尖度ＫXm（図３の部分(A)）から増加する（ＫSSm＞ＫXm）。尖度κは、ｎ次のモーメントμnから以下の式(5)で算定される高次統計量である。

Figure 0005152800
3 corresponds to the intensity of the component (αm · ψm (e jω )) that the noise suppression unit 26 subtracts from the acoustic signal VIN (frequency spectrum Xm (e jω )). Since the intensity within the range ASS changes to zero due to the suppression by the noise suppression unit 26, the intensity distribution of the intensity of the acoustic signal VOUT in which the noise component n is suppressed (part (C) in FIG. 3) is an intensity close to zero. Is increased compared to the frequency distribution before suppression of the noise component n (part (A) in FIG. 3). That is, the slope of the frequency distribution in a range where the intensity is near zero changes to a steep shape after being suppressed by the noise suppression unit 26. When kurtosis is introduced as a measure of the shape of the frequency distribution (steepness of inclination), the kurtosis KSSm (part (C) in FIG. 3) of the mth frame FR of the acoustic signal VOUT is the noise suppression unit. 26 suppresses the noise component n from the acoustic signal VIN, thereby increasing the kurtosis KXm (part (A) in FIG. 3) of the mth frame FR of the acoustic signal VIN before suppression (KSSm> KXm). The kurtosis κ is a higher-order statistic calculated from the n-th moment μn by the following equation (5).
Figure 0005152800

ミュージカルノイズには、高頻度でゼロの付近の強度になるという傾向がある。したがって、度数分布にて強度がゼロとなる度数が雑音成分ｎの抑圧の前後で増加するほど、雑音成分ｎの抑圧に起因して発生したミュージカルノイズが多いと評価できる。すなわち、雑音成分ｎの抑圧の前後における尖度κの変化（ＫXm→ＫSSm）の程度が大きいほど、雑音成分ｎの抑圧に起因して発生するミュージカルノイズは多い。例えば、音響信号ＶINの尖度ＫXmが同等である複数の場合を想定すると、尖度ＫSSmが大きい場合ほど、雑音成分ｎの抑圧後のミュージカルノイズは多いと評価される。また、例えば雑音成分ｎの抑圧後の尖度ＫSSmが同等である複数の場合を想定すると、雑音成分ｎの抑圧前の音響信号ＶINの尖度ＫXmが小さい場合ほど、雑音成分ｎの抑圧後のミュージカルノイズは多いと評価される。   There is a tendency for musical noise to have an intensity around zero with high frequency. Therefore, it can be evaluated that the more the frequency at which the intensity becomes zero in the frequency distribution increases before and after the suppression of the noise component n, the more musical noise is generated due to the suppression of the noise component n. That is, the greater the degree of change in kurtosis κ (KXm → KSSm) before and after suppression of noise component n, the more musical noise is generated due to suppression of noise component n. For example, assuming a plurality of cases where the kurtosis KXm of the acoustic signal VIN is equivalent, it is evaluated that the greater the kurtosis KSSm, the more musical noise after the suppression of the noise component n. Further, for example, assuming a plurality of cases where the kurtosis KSSm after suppression of the noise component n is equal, the smaller the kurtosis KXm of the acoustic signal VIN before suppression of the noise component n, the smaller the kurtosis KSSm after suppression of the noise component n. It is evaluated that there are many musical noises.

以上の傾向から、図１の指標算定部３２は、雑音成分ｎの抑圧前の音響信号ＶINの尖度ＫXmと雑音成分ｎの抑圧後の尖度ＫSSmとに応じて、第ｍ番目のフレームＦRにおけるミュージカルノイズの程度の指標となる雑音指標値σmを算定する。雑音指標値σmの算定には、音響信号ＶINから抽出されたＭ個の強度ｘi（ｘ1〜ｘM）の集合ｇXが使用される。図２に示すように、集合ｇXは、第ｍ番目のフレームＦRを最後として時間的に手前にあるｎt個のフレームＦRの各々の周波数スペクトルＸのｎf個の強度ｘiの集合（Ｍ＝ｎt×ｎf）としてフレームＦR毎に順次に特定される。雑音指標値σmの算定に使用される演算式の導出を以下に例示する。   From the above tendency, the index calculation unit 32 in FIG. 1 determines the m-th frame FR according to the kurtosis KXm of the acoustic signal VIN before suppression of the noise component n and the kurtosis KSSm of the noise component n after suppression. A noise index value σm, which is an index of the degree of musical noise at, is calculated. For calculating the noise index value σm, a set gX of M intensities xi (x1 to xM) extracted from the acoustic signal VIN is used. As shown in FIG. 2, the set gX is a set of nf intensities xi of the frequency spectrum X of each of the nt frames FR that are temporally preceding the mth frame FR (M = nt × nf) is sequentially specified for each frame FR. Derivation of an arithmetic expression used for calculating the noise index value σm will be exemplified below.

まず、雑音成分ｎの抑圧前の尖度ＫXmの算定について説明する。音響信号ＶINの強度（集合ｇXのＭ個の強度ｘ1〜ｘM）の度数分布は、以下の式(6)の関数Ｇa(ｘ;ｋ,θ)で近似される。

Figure 0005152800

式(6)における係数Ｃは、ガンマ関数Γ(k)を利用して以下のように定義される。
Figure 0005152800
 First, calculation of the kurtosis KXm before suppression of the noise component n will be described. The frequency distribution of the intensity of the acoustic signal VIN (M intensities x1 to xM of the set gX) is approximated by a function Ga (x; k, θ) of the following equation (6).
Figure 0005152800
The coefficient C in Equation (6) is defined as follows using the gamma function Γ (k).
Figure 0005152800

２次のモーメントμ2の定義式における関数Ｐ(x)に式(6)の関数Ｇa(ｘ;ｋ,θ)を代入することで以下の式(7)が導出される。

Figure 0005152800
By substituting the function Ga (x; k, θ) of the equation (6) into the function P (x) in the definition equation of the second-order moment μ2, the following equation (7) is derived.
Figure 0005152800

２次のモーメントμ2の導出と同様に、４次のモーメントμ4の定義式における関数Ｐ(x)に式(6)の関数Ｇa(ｘ;ｋ,θ)を代入することで以下の式(8)が導出される。

Figure 0005152800
Similar to the derivation of the second moment μ2, by substituting the function Ga (x; k, θ) of the equation (6) into the function P (x) in the definition equation of the fourth moment μ4, the following equation (8 ) Is derived.
Figure 0005152800

式(7)の２次のモーメントμ2と式(8)の４次のモーメントμ4とを式(5)に代入することで、雑音成分ｎの抑圧前の音響信号ＶINの尖度ＫXmは以下の式(9)のように定義される。

Figure 0005152800

式(6)に付記した変数ｋおよび変数γの定義式から理解されるように、式(9)における変数ｋ（変数ｋを定義する変数γ）の算定には、集合ｇXの強度ｘ1〜ｘMが使用される。 By substituting the second-order moment μ2 of equation (7) and the fourth-order moment μ4 of equation (8) into equation (5), the kurtosis KXm of the acoustic signal VIN before suppression of the noise component n is It is defined as equation (9).
Figure 0005152800
As understood from the definition formulas of the variable k and the variable γ appended to the equation (6), the calculation of the variable k (the variable γ defining the variable k) in the equation (9) is performed with the intensities x1 to xM of the set gX. Is used.

次に、雑音成分ｎの抑圧後の尖度ＫSSmの算定について説明する。ガンマ関数Γ(k)の平均ｋ・θを正規化することで以下の式(10)が導出される。

Figure 0005152800

いま、推定雑音スペクトルψm(ejω)のＡ倍（Ａ・ψm(ejω)）を雑音抑圧部２６が周波数スペクトルＸm(ejω)から減算する場合（式(4a)の係数αmを係数Ａに設定した場合）を想定すると、雑音成分ｎ（推定雑音スペクトルψm(ejω)）の抑圧後の音響信号ＶOUTの強度の度数分布を近似する関数Ｇb(ｘ;ｋ,θ)は、関数Ｇa(ｘ;ｋ,θ)の定義式(6)の強度ｘを強度(ｘ＋Ａ)に置換した以下の式(11)のように推定される。
Figure 0005152800
 Next, calculation of the kurtosis KSSm after suppression of the noise component n will be described. The following equation (10) is derived by normalizing the average k · θ of the gamma function Γ (k).
Figure 0005152800
Now, when the noise suppression unit 26 subtracts A times (A · ψm (e jω )) of the estimated noise spectrum ψm (e jω ) from the frequency spectrum Xm (e jω ) (the coefficient αm in equation (4a) is the coefficient A Assuming that the noise component n (estimated noise spectrum ψm (e jω )) is suppressed, the function Gb (x; k, θ) approximating the frequency distribution of the intensity of the acoustic signal VOUT is It is estimated as the following formula (11) in which the strength x in the definition formula (6) of (x; k, θ) is replaced with the strength (x + A).
Figure 0005152800

式(8)と同様に、４次のモーメントμ4の定義式における関数Ｐ(x)に式(11)の関数Ｇb(ｘ;ｋ,θ)を代入することで式(12)が導出される。

Figure 0005152800
Similar to the equation (8), the equation (12) is derived by substituting the function Gb (x; k, θ) of the equation (11) into the function P (x) in the definition equation of the fourth-order moment μ4. .
Figure 0005152800

式(12)の(ｘ＋Ａ)k-1は以下の式(13)のようにテイラー展開される。

Figure 0005152800
(X + A) k-1 in equation (12) is Taylor-expanded as in equation (13) below.
Figure 0005152800

式(13)の高次項を便宜的に無視したうえで式(12)に代入すると、４次のモーメントμ4を近似する以下の式(14)が導出される。

Figure 0005152800
If the higher order term of the equation (13) is ignored for convenience and substituted into the equation (12), the following equation (14) approximating the fourth-order moment μ4 is derived.
Figure 0005152800

２次のモーメントμ2についても同様に、定義式の関数Ｐ(x)（式(7)）に式(11)の関数Ｇb(ｘ;ｋ,θ)を代入したうえで式(13)の高次項を無視することで、以下の式(15)が導出される。

Figure 0005152800
Similarly for the second-order moment μ2, the function Gb (x; k, θ) of the equation (11) is substituted into the function P (x) (the equation (7)) of the defining equation, and then the high of the equation (13) By ignoring the next term, the following equation (15) is derived.
Figure 0005152800

式(14)の４次のモーメントμ4と式(15)の２次のモーメントμ2とを式(5)に代入することで、雑音成分ｎの抑圧後の尖度ＫSSmを変数ｋと係数（以下「抑圧係数」という）Ａとで定義する以下の式(16)が導出される。なお、式(16)の導出には式(10)を使用した。

Figure 0005152800
By substituting the fourth-order moment μ4 in equation (14) and the second-order moment μ2 in equation (15) into equation (5), the kurtosis KSSm after suppression of noise component n is expressed as variable k and coefficient The following expression (16) defined by A) (referred to as “suppression coefficient”) is derived. Note that the equation (10) was used to derive the equation (16).
Figure 0005152800

図４は、指標算定部３２の具体的な構成を示すブロック図である。図４に示すように、指標算定部３２は、相関特定部４２と指標決定部４４とを含んで構成される。相関特定部４２は、雑音成分ｎ（推定雑音スペクトルψm(ejω)）の抑圧の度合を表す抑圧係数Ａと尖度ＫXmおよび尖度ＫSSmに応じた尖度指標値Ｒmとの関係を特定する。 FIG. 4 is a block diagram illustrating a specific configuration of the index calculation unit 32. As shown in FIG. 4, the index calculating unit 32 includes a correlation specifying unit 42 and an index determining unit 44. The correlation specifying unit 42 specifies the relationship between the suppression coefficient A indicating the degree of suppression of the noise component n (estimated noise spectrum ψm (e jω )), the kurtosis KXm, and the kurtosis index value Rm corresponding to the kurtosis KSSm. .

尖度指標値Ｒmは、以下の式(17a)に示すように、尖度ＫXmに対する尖度ＫSSmの相対比ＫRm（ＫRm＝ＫSSm／ＫXm）を変数とする関数Ｆaで定義される。関数Ｆaは、相対比ＫRmに対して尖度指標値Ｒmが単調増加するように尖度指標値Ｒmと相対比ＫRmとの関係を定義する。
Ｒm＝Ｆa(ＫRm)＝Ｆa(ＫSSm／ＫXm) ……(17a) The kurtosis index value Rm is defined by a function Fa using a relative ratio KRm (KRm = KSSm / KXm) of the kurtosis KSSm to the kurtosis KXm as shown in the following equation (17a). The function Fa defines the relationship between the kurtosis index value Rm and the relative ratio KRm so that the kurtosis index value Rm monotonously increases with respect to the relative ratio KRm.
Rm = Fa (KRm) = Fa (KSSm / KXm) (17a)

図３を参照して説明したように、雑音成分ｎの抑圧後のミュージカルノイズが多いほど相対比ＫRmは大きくなる（尖度ＫXmから尖度ＫSSmへの変化度が大きくなる）から、尖度指標値Ｒmが大きいほど、雑音成分ｎの抑圧後のミュージカルノイズが多いと評価できる。換言すると、尖度指標値Ｒmが小さいほど（すなわち尖度ＫXmから尖度ＫSSmへの変化度が小さいほど）、雑音成分ｎの抑圧後のミュージカルノイズは少ないと評価できる。   As described with reference to FIG. 3, the relative ratio KRm increases (the degree of change from the kurtosis KXm to the kurtosis KSSm increases) as the musical noise after suppression of the noise component n increases. It can be evaluated that the larger the value Rm, the more musical noise after the noise component n is suppressed. In other words, it can be evaluated that the smaller the kurtosis index value Rm (ie, the smaller the degree of change from the kurtosis KXm to the kurtosis KSSm), the less the musical noise after suppression of the noise component n.

式(9)に示すように尖度ＫXmは変数ｋの関数であり、式(16)に示すように尖度ＫSSmは変数ｋおよび抑圧係数Ａの関数であるから、関数Ｆaは、以下の式(17b)に示すように、尖度指標値Ｒmと変数ｋおよび抑圧係数Ａとの関係を定義する関数に相当する。
Ｒm＝Ｆa(ＫSSm／ＫXm)＝Ｆa(Ａ,ｋ) ……(17b) The kurtosis KXm is a function of the variable k as shown in the equation (9), and the kurtosis KSSm is a function of the variable k and the suppression coefficient A as shown in the equation (16). As shown in (17b), this corresponds to a function that defines the relationship between the kurtosis index value Rm, the variable k, and the suppression coefficient A.
Rm = Fa (KSSm / KXm) = Fa (A, k) (17b)

第ｍ番目の１個のフレームＦRに着目すると、変数ｋは、集合ｇX（第ｍ番目を最後とするｎt個のフレームＦR）のＭ個の強度ｘ1〜ｘMから算定される固定値である。したがって、尖度指標値Ｒmは、以下の式(17c)に示すように、抑圧係数Ａを変数とした関数Ｆaで定義される。
Ｒm＝Ｆa(Ａ) ……(17c)
図４の相関特定部４２は、集合ｇXのＭ個の強度ｘ1〜ｘMから算定した変数ｋを式(17b)に代入することで、尖度指標値Ｒmと抑圧係数Ａとの関係を定義する式(17c)の関数Ｆaを特定する。変数ｋはフレームＦR毎に変化するから、相関特定部４２は、フレームＦR毎に関数Ｆaを特定する。 Focusing on the mth frame FR, the variable k is a fixed value calculated from M intensities x1 to xM of the set gX (nt frames FR with mth as the last). Therefore, the kurtosis index value Rm is defined by a function Fa using the suppression coefficient A as a variable, as shown in the following equation (17c).
Rm = Fa (A) (17c)
The correlation specifying unit 42 in FIG. 4 defines the relationship between the kurtosis index value Rm and the suppression coefficient A by substituting the variable k calculated from the M intensities x1 to xM of the set gX into the equation (17b). The function Fa in the equation (17c) is specified. Since the variable k changes for each frame FR, the correlation specifying unit 42 specifies the function Fa for each frame FR.

図４の指標決定部４４は、相関特定部４２が特定した関数Ｆaで定義される尖度指標値Ｒmが所定値Ｒrefに合致するときの抑圧係数Ａを雑音指標値σmとして決定する。すなわち、指標決定部４４は、以下の式(18)の演算を実行することでフレームＦR毎に雑音指標値σmを算定する。式(18)における演算子Ｆa-1は、関数Ｆaの逆写像を意味する。
σm＝Ｆa-1(Ｒref) ……(18) The index determining unit 44 in FIG. 4 determines the suppression coefficient A when the kurtosis index value Rm defined by the function Fa specified by the correlation specifying unit 42 matches the predetermined value Rref as the noise index value σm. That is, the index determination unit 44 calculates the noise index value σm for each frame FR by executing the following expression (18). The operator Fa -1 in equation (18) means the inverse mapping of the function Fa.
σm = Fa -1 (Rref) (18)

以上のように、雑音指標値σmは、雑音抑圧部２６による雑音成分ｎの抑圧後に発生するミュージカルノイズを所定の度合に制御する（すなわち尖度指標値Ｒmを所定値Ｒrefに制御する）ための係数αm（式(4a)）の数値に相当する。また、尖度指標値Ｒmが増加するほど雑音指標値σmは大きい数値となるから、雑音指標値σmには、尖度指標値Ｒmと同様に、抑圧係数Ａのもとで音響信号ＶINの雑音成分ｎを抑圧した場合に発生するミュージカルノイズの程度の指標としての意義もある。すなわち、音響信号ＶINは、雑音指標値σmが大きいほどミュージカルノイズを発生させ易い特性と評価され、雑音指標値σmが小さいほどミュージカルノイズを発生させ難い特性と評価される。以上のように、尖度指標値Ｒmおよび雑音指標値σmは、抑圧係数Ａのもとで雑音成分ｎを抑圧したときに音響信号ＶOUTに発生するミュージカルノイズの程度を定量的に表す指標として機能する。   As described above, the noise index value σm is used to control the musical noise generated after the noise suppression unit 26 suppresses the noise component n to a predetermined degree (that is, to control the kurtosis index value Rm to the predetermined value Rref). This corresponds to the numerical value of the coefficient αm (formula (4a)). Since the noise index value σm increases as the kurtosis index value Rm increases, the noise index value σm includes the noise of the acoustic signal VIN under the suppression coefficient A as in the kurtosis index value Rm. There is also a significance as an index of the degree of musical noise generated when the component n is suppressed. That is, the acoustic signal VIN is evaluated as a characteristic that is more likely to generate musical noise as the noise index value σm is larger, and is evaluated as a characteristic that is less likely to generate musical noise as the noise index value σm is smaller. As described above, the kurtosis index value Rm and the noise index value σm function as indexes that quantitatively represent the degree of musical noise generated in the acoustic signal VOUT when the noise component n is suppressed under the suppression coefficient A. To do.

図１の抑圧制御部３６は、雑音抑圧部２６が雑音成分ｎの抑圧（式(4a)および式(4b)）に使用する係数αmおよび係数βmを、指標算定部３２が算定した雑音指標値σmとＳＮ比算定部３４が算定したＳＮ比ξmとに応じて可変に制御する。抑圧制御部３６の具体的な動作を以下に説明する。   The suppression control unit 36 in FIG. 1 includes the noise index value calculated by the index calculation unit 32 using the coefficient αm and the coefficient βm used by the noise suppression unit 26 for the suppression of the noise component n (equations (4a) and (4b)). Control is variably performed according to σm and the SN ratio ξm calculated by the SN ratio calculation unit 34. A specific operation of the suppression control unit 36 will be described below.

抑圧制御部３６は、例えば以下の式(19)を演算することで雑音指標値σmとＳＮ比ξmとに応じた係数αmを算定する。式(19)の係数ａg1や係数ａg2は、例えば音響信号ＶOUTのミュージカルノイズが効果的に低減されるように実験的または統計的に設定された正数である。
αm＝ｇ(σm,ξm)＝−ａg1・σm２＋ａg2・ξm ……(19)
式(19)から理解されるように、雑音指標値σmが増加するほど係数αmは減少する。したがって、雑音抑圧部２６による雑音成分ｎの抑圧でミュージカルノイズが発生する可能性が高い（雑音指標値σmが大きい）ほど、雑音抑圧部２６による周波数スペクトルＸm(ejω)からの減算量（αm・ψm(ejω)）は低減される。 The suppression control unit 36 calculates the coefficient αm according to the noise index value σm and the SN ratio ξm, for example, by calculating the following equation (19). The coefficient ag1 and the coefficient ag2 in the equation (19) are positive numbers experimentally or statistically set so that, for example, the musical noise of the acoustic signal VOUT is effectively reduced.
αm = g (σm, ξm) = − ag1 · σm 2 + ag2 · ξm (19)
As understood from the equation (19), the coefficient αm decreases as the noise index value σm increases. Accordingly, the more likely that the musical noise is generated by the suppression of the noise component n by the noise suppression unit 26 (the noise index value σm is larger), the subtraction amount (αm) from the frequency spectrum Xm (e jω ) by the noise suppression unit 26. • ψm (e jω )) is reduced.

例えば、図５の部分(A)および部分(B)のように音響信号ＶINの尖度ＫXmが抑圧後の尖度ＫSSmに対して充分に小さい場合（例えば、音響信号ＶINの尖度ＫXmが正規分布と比較して低いサブガウス的な特性の場合）には雑音指標値σm（尖度指標値Ｒm）が大きい数値となるから、係数αmが小さい数値に設定されることで周波数スペクトルＸm(ejω)からの減算量（αm・ψm(ejω)は低減される。一方、図６の部分(A)および部分(B)のように音響信号ＶINの尖度ＫXmが大きい場合（例えば、音響信号ＶINの尖度ＫXmが正規分布と比較して高いスーパーガウス的な特性の場合）には雑音指標値σmが小さい数値となるから、係数αmが大きい数値に設定されることで周波数スペクトルＸm(ejω)からの減算量（αm・ψm(ejω)）は増加する。以上のように係数αmは雑音指標値σmに応じて可変に設定されるから、式(19)におけるＳＮ比ξmの作用を便宜的に無視すると、雑音抑圧部２６による実際の処理後の音響信号ＶOUTの尖度指標値Ｒmは所定値（目標値）Ｒrefに略合致する。 For example, when the kurtosis KXm of the acoustic signal VIN is sufficiently smaller than the kurtosis KSSm after the suppression as in the parts (A) and (B) in FIG. 5 (for example, the kurtosis KXm of the acoustic signal VIN is normal) Since the noise index value σm (kurtosis index value Rm) is a large value in the case of sub-Gaussian characteristics lower than the distribution, the frequency spectrum Xm (e jω is set by setting the coefficient αm to a small value. ) Is reduced (αm · ψm (e jω ). On the other hand, when the kurtosis KXm of the acoustic signal VIN is large as in the parts (A) and (B) of FIG. When the kurtosis KXm of VIN is a super Gaussian characteristic higher than that of the normal distribution), the noise index value σm is a small numerical value. Therefore, the frequency spectrum Xm (e subtraction amount from jω) (αm · ψm (e jω)) is increased. above coefficient .alpha.m as the noise index value σm Accordingly, if the action of the SN ratio ξm in the equation (19) is ignored for convenience, the kurtosis index value Rm of the acoustic signal VOUT actually processed by the noise suppression unit 26 is a predetermined value (target) Value) Rref substantially matches.

また、式(19)から理解されるように、ＳＮ比算定部３４の算定したＳＮ比ξmが増加するほど係数αmは増加する。したがって、音響信号ＶINのＳＮ比ξmが大きいほど（つまり、信号成分ｓの強度が雑音成分ｎの強度に対して高いほど）、雑音抑圧部２６による周波数スペクトルＸm(ejω)からの減算量（αm・ψm(ejω)）は増加する。 Further, as understood from the equation (19), the coefficient αm increases as the SN ratio ξm calculated by the SN ratio calculator 34 increases. Therefore, as the SN ratio ξm of the acoustic signal VIN is larger (that is, the intensity of the signal component s is higher than the intensity of the noise component n), the subtraction amount from the frequency spectrum Xm (e jω ) by the noise suppression unit 26 ( αm · ψm (e jω )) increases.

また、抑圧制御部３６は、例えば式(20)を演算することで雑音指標値σmとＳＮ比ξmとに応じた係数βmを算定する。式(20)の係数ａh1や係数ａh2は、式(19)の係数ａg1や係数ａg2と同様に、例えば音響信号ＶOUTのミュージカルノイズが有効に低減されるように実験的または統計的に設定された正数である。
βm＝ｈ（σm,ξm）＝−ａh1・σm２＋ａh2・ξm ……(20) In addition, the suppression control unit 36 calculates a coefficient βm corresponding to the noise index value σm and the SN ratio ξm, for example, by calculating Equation (20). The coefficients ah1 and ah2 in the equation (20) are set experimentally or statistically so that, for example, the musical noise of the acoustic signal VOUT is effectively reduced, like the coefficients ag1 and ag2 in the equation (19). It is a positive number.
βm = h (σm, ξm) = − ah1 · σm 2 + ah2 · ξm (20)

式(20)から理解されるように、雑音指標値σmが増加するほど係数βmは減少する。したがって、雑音抑圧部２６による雑音成分ｎの抑圧でミュージカルノイズが発生する度合が高い（雑音指標値σmが大きい）ほど、周波数スペクトルＸm(ejω)の強度|Ｘm(ejω)|２が推定雑音スペクトルψm(ejω)と係数αmとの乗算（αm・ψm(ejω)）を下回る周波数帯域の成分の強度（βm・ψm(ejω)）は減少する。さらに、ＳＮ比ξmが増加するほど係数βmは増加する。したがって、音響信号ＶINのＳＮ比ξmが小さいほど、周波数スペクトルＸm(ejω)の強度|Ｘm(ejω)|２が推定雑音スペクトルψm(ejω)と係数αmとの乗算（αm・ψm(ejω)）を下回る周波数帯域の成分の強度（βm・ψm(ejω)）は減少する。 As understood from the equation (20), the coefficient βm decreases as the noise index value σm increases. Accordingly, degree of musical noise suppression of the noise component n by the noise suppressor 26 generates a high (noise index value σm is large) as the intensity of the frequency spectrum Xm (e jω) | Xm ( e jω) | 2 estimate The intensity (βm · ψm (e jω )) of the frequency band component lower than the multiplication (αm · ψm (e jω )) of the noise spectrum ψm (e jω ) and the coefficient αm decreases. Furthermore, the coefficient βm increases as the SN ratio ξm increases. Therefore, as the SN ratio ξm of the sound signal VIN is small, the intensity of the frequency spectrum Xm (e jω) | Xm ( e jω) | 2 multiplies the estimated noise spectrum ψm (e jω) and the coefficient αm (αm · ψm ( e j [omega])) the intensity of the component of the frequency band below (βm · ψm (e jω) ) is reduced.

以上に説明したように、本形態においては、雑音抑圧部２６による雑音成分ｎの抑圧の度合（αm・ψm(ejω)）が雑音指標値σmに応じて可変に制御される。さらに具体的には、雑音指標値σmが大きいほど雑音抑圧部２６による抑圧の度合（減算量）は低減される。したがって、雑音成分ｎの抑圧の度合が固定された従来の技術と比較して、音響信号ＶINが収録された環境（音響信号ＶINの特性）に拘わらず、ミュージカルノイズの発生を有効に抑制しながら音響信号ＶINの雑音成分ｎを効果的に抑圧することが可能である。本形態による効果の具体例を以下に説明する。 As described above, in this embodiment, the degree of suppression of the noise component n (αm · ψm (e jω )) by the noise suppression unit 26 is variably controlled according to the noise index value σm. More specifically, the degree of suppression (subtraction amount) by the noise suppression unit 26 decreases as the noise index value σm increases. Therefore, as compared with the conventional technique in which the degree of suppression of the noise component n is fixed, the generation of musical noise is effectively suppressed regardless of the environment (characteristics of the acoustic signal VIN) in which the acoustic signal VIN is recorded. It is possible to effectively suppress the noise component n of the acoustic signal VIN. A specific example of the effect of this embodiment will be described below.

雑音成分ｎが充分に抑制されるように係数αmを高目の固定値に設定した構成では、雑音成分ｎを充分に抑制することは確かに可能であるが、例えば音響信号ＶINが図５の部分(A)の特性である場合（つまり、ミュージカルノイズが発生し易い場合）に、雑音成分ｎの過剰な抑圧に起因して音響信号ＶOUTに顕著なミュージカルノイズが発生し易いという問題がある。本形態においては、図５の部分(A)および部分(B)のように雑音指標値σmが高い場合（音響信号ＶOUTにミュージカルノイズが発生し易い場合）に雑音抑圧部２６による抑圧の度合が低減されるから、音響信号ＶOUTのミュージカルノイズは有効に抑制される。   In the configuration in which the coefficient αm is set to a high fixed value so that the noise component n is sufficiently suppressed, it is certainly possible to sufficiently suppress the noise component n. For example, the acoustic signal VIN is shown in FIG. In the case of the characteristic of the portion (A) (that is, when musical noise is likely to occur), there is a problem that remarkable musical noise is likely to occur in the acoustic signal VOUT due to excessive suppression of the noise component n. In this embodiment, when the noise index value σm is high as in the parts (A) and (B) of FIG. 5 (when musical noise is likely to occur in the acoustic signal VOUT), the degree of suppression by the noise suppression unit 26 is increased. As a result, the musical noise of the acoustic signal VOUT is effectively suppressed.

一方、ミュージカルノイズが適切に抑制されるように係数αmを低目の固定値に設定した構成では、音響信号ＶOUTのミュージカルノイズを抑制することは確かに可能である。しかし、音響信号ＶINが図６の部分(A)の特性である場合には、雑音成分ｎの抑圧の度合を増加させても音響信号ＶOUTにミュージカルノイズは発生し難い状況であるにも拘わらず、雑音成分ｎの抑圧の度合が制限される（抑圧が不足する）という問題がある。本形態においては、図６の部分(A)および部分(B)のように雑音指標値σmが低い場合（音響信号ＶOUTにミュージカルノイズが発生し難い場合）に雑音抑圧部２６による抑圧の度合が増加するから、雑音成分ｎは音響信号ＶOUTにて効果的に抑制される。   On the other hand, in the configuration in which the coefficient αm is set to a low fixed value so that the musical noise is appropriately suppressed, it is certainly possible to suppress the musical noise of the acoustic signal VOUT. However, if the acoustic signal VIN has the characteristic of part (A) in FIG. 6, even if the degree of suppression of the noise component n is increased, musical noise is hardly generated in the acoustic signal VOUT. There is a problem that the degree of suppression of the noise component n is limited (suppression is insufficient). In the present embodiment, when the noise index value σm is low as in the parts (A) and (B) in FIG. 6 (when it is difficult to generate musical noise in the acoustic signal VOUT), the degree of suppression by the noise suppression unit 26 is reduced. Since it increases, the noise component n is effectively suppressed in the acoustic signal VOUT.

ところで、音響信号ＶINのＳＮ比ξmが高い場合には、雑音成分ｎの抑圧の度合が高くても、音響信号ＶOUTのミュージカルノイズは受聴者に知覚され難いという傾向がある。本形態においては、音響信号ＶINのＳＮ比ξmに応じて雑音抑圧部２６による抑圧の度合（係数αm）が制御される。さらに詳述すると、ＳＮ比ξmが高いほど雑音抑圧部２６による抑圧の度合（係数αm）は増加する。したがって、ＳＮ比ξmが高いことでミュージカルノイズが知覚され難い環境では、ミュージカルノイズの抑制に優先して雑音成分ｎが有効に抑制されるという利点がある。換言すると、音響信号ＶINのＳＮ比ξmが低い場合に雑音抑圧部２６による抑圧の度合（係数αm）が低減されるから、ＳＮ比ξmが低いことでミュージカルノイズが特に知覚され易い環境ではミュージカルノイズが優先的に抑制される。もっとも、ＳＮ比算定部３４を省略した構成（雑音抑圧部２６の制御に雑音指標値σmのみが反映される構成）も本発明では採用される。   By the way, when the SN ratio ξm of the acoustic signal VIN is high, even if the degree of suppression of the noise component n is high, the musical noise of the acoustic signal VOUT tends not to be perceived by the listener. In this embodiment, the degree of suppression (coefficient αm) by the noise suppression unit 26 is controlled in accordance with the SN ratio ξm of the acoustic signal VIN. More specifically, the degree of suppression by the noise suppression unit 26 (coefficient αm) increases as the SN ratio ξm increases. Therefore, in an environment where musical noise is difficult to perceive due to a high SN ratio ξm, there is an advantage that noise component n is effectively suppressed in preference to suppression of musical noise. In other words, when the SN ratio ξm of the acoustic signal VIN is low, the degree of suppression (coefficient αm) by the noise suppression unit 26 is reduced. Therefore, in an environment where musical noise is particularly perceptible due to the low SN ratio ξm, the musical noise is reduced. Is preferentially suppressed. However, a configuration in which the SN ratio calculation unit 34 is omitted (a configuration in which only the noise index value σm is reflected in the control of the noise suppression unit 26) is also employed in the present invention.

なお、音響信号ＶOUTのミュージカルノイズは、推定雑音スペクトルψm(ejω)の減算（式(4a)）を主要因として発生するから、ミュージカルノイズを低減するという観点からすると、推定雑音スペクトルψm(ejω)の減算に適用される係数αmを雑音指標値σmに応じて可変に制御する構成が重要である。したがって、係数βmは所定値に固定する（雑音指標値σmに依存しない）という構成も採用される。しかし、係数βmを固定した構成では、周波数スペクトルＹm(ejω)のうち式(4a)が適用される帯域と式(4b)が適用される帯域とで強度の相違が過大となり、音響信号ＶOUTの再生音が聴感上で不自然な特性となる可能性がある。本形態においては、係数βmが係数αmと同様に雑音指標値σmとＳＮ比ξmとに応じて可変に制御されるから、式(4a)が適用される帯域と式(4b)が適用される帯域との強度の相違が抑制される。したがって、係数βmを固定した構成と比較して、再生音が聴感上で自然な印象に知覚される音響信号ＶOUTを生成できるという利点がある。 Note that the musical noise of the acoustic signal VOUT is generated mainly by subtraction of the estimated noise spectrum ψm (e jω ) (equation (4a)). From the viewpoint of reducing the musical noise, the estimated noise spectrum ψm (e The configuration in which the coefficient αm applied to the subtraction of jω ) is variably controlled according to the noise index value σm is important. Therefore, a configuration in which the coefficient βm is fixed to a predetermined value (not depending on the noise index value σm) is also adopted. However, in the configuration in which the coefficient βm is fixed, the difference in intensity between the band to which the expression (4a) is applied and the band to which the expression (4b) is applied is excessive in the frequency spectrum Ym (e jω ), and the acoustic signal VOUT May cause unnatural characteristics in the sense of hearing. In the present embodiment, the coefficient βm is variably controlled in accordance with the noise index value σm and the SN ratio ξm in the same manner as the coefficient αm. Therefore, the band to which the expression (4a) is applied and the expression (4b) are applied. Difference in intensity from the band is suppressed. Therefore, there is an advantage that the acoustic signal VOUT in which the reproduced sound is perceived as a natural impression on hearing can be generated as compared with the configuration in which the coefficient βm is fixed.

＜Ｂ：第２実施形態＞
図７は、本発明の第２実施形態に係る指標算定部３２のブロック図である。図７に示すように、本形態の指標算定部３２は、第１尖度算定部５１と第２尖度算定部５２と算定部５４とを含んで構成される。なお、本形態のうち第１実施形態と共通する要素については、以上と同じ符号を付して各々の詳細な説明を適宜に省略する。 <B: Second Embodiment>
FIG. 7 is a block diagram of the index calculation unit 32 according to the second embodiment of the present invention. As shown in FIG. 7, the index calculation unit 32 of the present embodiment includes a first kurtosis calculation unit 51, a second kurtosis calculation unit 52, and a calculation unit 54. In addition, about the element which is common in 1st Embodiment among this form, the same code | symbol as the above is attached | subjected and each detailed description is abbreviate | omitted suitably.

図７の第１尖度算定部５１は、音響信号ＶINの各フレームＦRについて尖度ＫXmを算定する。例えば、第１尖度算定部５１は、周波数スペクトルＸm(ejω)の時系列から抽出された集合ｇXのＭ個の強度ｘ1〜ｘMについて式(9)の演算を実行することで音響信号ＶINのフレームＦR毎に尖度ＫXmを算定する。同様に、第２尖度算定部５２は、雑音抑圧部２６による雑音成分ｎの抑圧後の各フレームＦRについて尖度ＫSSmを算定する。例えば、第２尖度算定部５２は、雑音抑圧部２６による実際の処理後の周波数スペクトルＹm(ejω)の時系列から図２の方法で抽出したＭ個の強度ｘ1〜ｘMについて式(9)の演算を実行することで、音響信号ＶOUTのフレームＦR毎に尖度ＫSSmを算定する。 The first kurtosis calculation unit 51 in FIG. 7 calculates the kurtosis KXm for each frame FR of the acoustic signal VIN. For example, the first kurtosis calculation unit 51 performs the calculation of the expression (9) on the M intensities x1 to xM of the set gX extracted from the time series of the frequency spectrum Xm (e jω ) to thereby obtain the acoustic signal VIN. The kurtosis KXm is calculated for each frame FR. Similarly, the second kurtosis calculation unit 52 calculates the kurtosis KSSm for each frame FR after the noise component n is suppressed by the noise suppression unit 26. For example, the second kurtosis calculation unit 52 uses the formula (9) for M intensities x1 to xM extracted by the method of FIG. 2 from the time series of the frequency spectrum Ym (e jω ) actually processed by the noise suppression unit 26. ), The kurtosis KSSm is calculated for each frame FR of the acoustic signal VOUT.

図７の算定部５４は、第１尖度算定部５１が算定した尖度ＫXmと第２尖度算定部５２が算定した尖度ＫSSmとから雑音指標値σmを算定する。さらに詳述すると、算定部５４は、尖度ＫXmに対する尖度ＫSSmの相対比ＫRmを算定し、相対比ＫRmを関数Ｆbに代入することで雑音指標値σmを算定する（σm＝Ｆb(ＫRm)＝Ｆb(ＫSSm／ＫXm)）。   7 calculates the noise index value σm from the kurtosis KXm calculated by the first kurtosis calculation unit 51 and the kurtosis KSSm calculated by the second kurtosis calculation unit 52. More specifically, the calculating unit 54 calculates the relative ratio KRm of the kurtosis KSSm to the kurtosis KXm, and calculates the noise index value σm by substituting the relative ratio KRm into the function Fb (σm = Fb (KRm)). = Fb (KSSm / KXm)).

関数Ｆbは、相対比ＫRmに対して雑音指標値σmが単調増加するように相対比ＫRmと雑音指標値σmとの関係を定義する。したがって、第１実施形態と同様に、雑音指標値σmは、雑音成分ｎの抑圧に起因してミュージカルノイズが発生する程度を定量的に評価するための指標として機能する。例えば、指標算定部３２の算定した雑音指標値σmが大きいほど（相対比ＫRmが大きいほど）、雑音成分ｎの抑圧後のミュージカルノイズが多いと評価できる。   The function Fb defines the relationship between the relative ratio KRm and the noise index value σm so that the noise index value σm monotonously increases with respect to the relative ratio KRm. Therefore, as in the first embodiment, the noise index value σm functions as an index for quantitatively evaluating the degree to which musical noise is generated due to suppression of the noise component n. For example, it can be evaluated that the greater the noise index value σm calculated by the index calculation unit 32 (the greater the relative ratio KRm), the more musical noise after suppression of the noise component n.

抑圧制御部３６は、雑音抑圧部２６が第ｍ番目のフレームＦRの処理に適用する係数αmおよび係数βmを、直前（第(m-1)番目）のフレームＦRについて指標算定部３２が算定した雑音指標値σm-1に応じて可変に設定する。係数αmや係数βmの算定には第１実施形態と同様の方法（式(19)，式(20)）が使用される。したがって、本形態においても第１実施形態と同様の効果が実現される。   In the suppression control unit 36, the index calculation unit 32 calculates the coefficient αm and the coefficient βm that the noise suppression unit 26 applies to the processing of the mth frame FR for the immediately previous ((m−1) th) frame FR. It is variably set according to the noise index value σm-1. For the calculation of the coefficient αm and the coefficient βm, the same method as in the first embodiment (Equation (19), Equation (20)) is used. Therefore, also in this embodiment, the same effect as that of the first embodiment is realized.

なお、以上の形態においては係数αmや係数βmが直前のフレームＦRの雑音指標値σm-1から算定される。一方、第１実施形態においては、第ｍ番目のフレームＦRに適用される係数αmや係数βmが、当該フレームＦRの音響信号ＶINから算定された雑音指標値σmに応じて設定される。したがって、雑音成分ｎの抑制の度合を音響信号ＶINの特性の変化（音響信号ＶINが収録された環境の変化）に対して迅速に適応させるという観点からすると、第２実施形態よりも第１実施形態が好適である。   In the above embodiment, the coefficient αm and the coefficient βm are calculated from the noise index value σm−1 of the immediately preceding frame FR. On the other hand, in the first embodiment, the coefficient αm and coefficient βm applied to the m-th frame FR are set according to the noise index value σm calculated from the acoustic signal VIN of the frame FR. Therefore, from the viewpoint of quickly adapting the degree of suppression of the noise component n to changes in the characteristics of the acoustic signal VIN (changes in the environment in which the acoustic signal VIN is recorded), the first implementation is more effective than the second embodiment. The form is preferred.

もっとも、第２実施形態においても、第ｍ番目のフレームＦRから算定された雑音指標値σmを当該フレームＦRの雑音成分ｎの抑制に使用する構成が採用される。例えば、係数αmや係数βmを仮定的に所定値（初期値）に設定した状態で雑音抑圧部２６が第ｍ番目のフレームＦRについて雑音成分ｎを抑圧し、抑圧後の第ｍ番目のフレームＦRについて指標算定部３２が算定した雑音指標値σmを、当該フレームＦRの実際の雑音成分ｎの抑圧に適用される係数αmや係数βmの算定に抑圧制御部３６が適用する。   However, in the second embodiment, a configuration in which the noise index value σm calculated from the m-th frame FR is used for suppressing the noise component n of the frame FR is employed. For example, the noise suppression unit 26 suppresses the noise component n for the mth frame FR in a state where the coefficient αm and the coefficient βm are assumed to be predetermined values (initial values), and the mth frame FR after the suppression. The suppression control unit 36 applies the noise index value σm calculated by the index calculation unit 32 for the calculation of the coefficient αm and the coefficient βm applied to the suppression of the actual noise component n of the frame FR.

第１実施形態および第２実施形態の例示から理解されるように、本発明は、雑音成分ｎの抑圧の前後の尖度（ＫXm，ＫSSm）を実際には算定せずに雑音指標値σmを算定する構成（第１実施形態）、および、雑音成分ｎの抑圧の前後の尖度（ＫXm，ＫSSm）を実際に算定したうえで雑音指標値σmを算定する構成（第２実施形態）の双方を包含する。   As understood from the illustrations of the first and second embodiments, the present invention calculates the noise index value σm without actually calculating the kurtosis (KXm, KSSm) before and after the suppression of the noise component n. Both the configuration for calculating (first embodiment) and the configuration for calculating the noise index value σm after actually calculating the kurtosis (KXm, KSSm) before and after suppression of the noise component n (second embodiment) Is included.

＜Ｃ：変形例＞
以上に例示した各形態には様々な変形が加えられる。具体的な変形の態様を例示すれば以下の通りである。なお、以下の例示から２以上の態様を任意に選択して組合せてもよい。 <C: Modification>
Various modifications can be made to each of the forms exemplified above. An example of a specific modification is as follows. Two or more aspects may be arbitrarily selected from the following examples and combined.

（１）変形例１
尖度ＫXmまたは尖度ＫSSmと雑音指標値σm（第１実施形態の尖度指標値Ｒm）との関係（尖度ＫXmや尖度ＫSSmから雑音指標値σmや尖度指標値Ｒmを算定する方法）は本発明において任意である。例えば、音響信号ＶOUTにおけるミュージカルノイズの発生の度合が雑音成分ｎの抑圧の前後における尖度の変化度（ＫXm→ＫSSm）に反映されるという傾向を考慮すると、抑圧前の尖度ＫXmと抑圧後の尖度ＫSSmとの差分|ＫSSm−ＫXm|に応じて雑音指標値σmや尖度指標値Ｒmを算定する構成も採用される。また、相対比ＫRmと尖度指標値Ｒmとの関係（関数Ｆaの内容）や相対比ＫRmと雑音指標値σmとの関係（関数Ｆbの内容）は適宜に変更される。例えば、第１実施形態においては、相対比ＫRmを尖度指標値Ｒmとする構成（Ｒm＝ＫRm）や、相対比ＫRmに所定の係数を加減または乗除して尖度指標値Ｒmを算定する構成が採用される。同様に、第２実施形態においては、相対比ＫRmを雑音指標値σmとして出力する構成（ＫRm＝σm）や、相対比ＫRmに所定の係数を加減または乗除して雑音指標値σmを算定する構成が採用される。 (1) Modification 1
Relationship between kurtosis KXm or kurtosis KSSm and noise index value σm (kurtosis index value Rm of the first embodiment) (method of calculating noise index value σm and kurtosis index value Rm from kurtosis KXm and kurtosis KSSm ) Is optional in the present invention. For example, considering the tendency that the degree of occurrence of musical noise in the acoustic signal VOUT is reflected in the degree of change in kurtosis before and after suppression of the noise component n (KXm → KSSm), kurtosis KXm before suppression and after suppression A configuration for calculating the noise index value σm and the kurtosis index value Rm according to the difference | KSSm−KXm | Further, the relationship between the relative ratio KRm and the kurtosis index value Rm (contents of the function Fa) and the relationship between the relative ratio KRm and the noise index value σm (contents of the function Fb) are appropriately changed. For example, in the first embodiment, a configuration in which the relative ratio KRm is the kurtosis index value Rm (Rm = KRm), or a configuration in which the kurtosis index value Rm is calculated by adding or subtracting or multiplying the relative ratio KRm by a predetermined coefficient. Is adopted. Similarly, in the second embodiment, a configuration in which the relative ratio KRm is output as the noise index value σm (KRm = σm), or a configuration in which the noise index value σm is calculated by adding or subtracting a predetermined coefficient to or from the relative ratio KRm. Is adopted.

以上の各形態においては尖度ＫXmおよび尖度ＫSSmの相対比ＫRmと雑音指標値σm（あるいは第１実施形態の尖度指標値Ｒm）との相関に着目したが、雑音成分ｎの抑圧後のミュージカルノイズの発生の程度は、相対比ＫRmの対数値に対して特に顕著な相関を示すという傾向がある。したがって、相対比ＫRmの対数値から雑音指標値σmを算定する構成（すなわち、以上の各形態における相対比ＫRmを相対比ＫRMの対数値に置換した構成）も好適である。相対比ＫRmの対数値を利用した構成によれば、ミュージカルノイズの発生の程度を、より正確に雑音指標値σmから評価できるという利点がある。   In each of the above embodiments, attention is paid to the correlation between the relative ratio KRm of the kurtosis KXm and the kurtosis KSSm and the noise index value σm (or the kurtosis index value Rm of the first embodiment). The degree of occurrence of musical noise tends to show a particularly significant correlation with the logarithmic value of the relative ratio KRm. Therefore, a configuration in which the noise index value σm is calculated from the logarithmic value of the relative ratio KRm (that is, a configuration in which the relative ratio KRm in each of the above embodiments is replaced with a logarithmic value of the relative ratio KRM) is also suitable. According to the configuration using the logarithmic value of the relative ratio KRm, there is an advantage that the degree of occurrence of musical noise can be more accurately evaluated from the noise index value σm.

（２）変形例２
以上の各形態においては、音響信号ＶINの尖度ＫXmと雑音成分ｎの抑圧後の尖度ＫSSmとの双方を雑音指標値σmの算定に使用したが、尖度ＫXmおよび尖度ＫSSmの一方のみを使用して雑音指標値σmを算定する構成も採用される。例えば、雑音成分ｎの抑圧前の尖度ＫXmが低いほど音響信号ＶOUTにミュージカルノイズが発生し易いという傾向のみに着目すると、音響信号ＶINの尖度ＫXmが低いほど雑音指標値σmが増加するように指標算定部３２が雑音指標値σmを算定する構成（雑音指標値σmが尖度ＫSSmに依存しない構成）も好適である。 (2) Modification 2
In the above embodiments, both the kurtosis KXm of the acoustic signal VIN and the kurtosis KSSm after suppression of the noise component n are used for the calculation of the noise index value σm, but only one of the kurtosis KXm and the kurtosis KSSm is used. A configuration is also employed in which the noise index value σm is calculated using. For example, focusing only on the tendency that musical noise is more likely to occur in the acoustic signal VOUT as the kurtosis KXm before suppression of the noise component n is lower, the noise index value σm increases as the kurtosis KXm of the acoustic signal VIN is lower. Further, a configuration in which the index calculation unit 32 calculates the noise index value σm (a configuration in which the noise index value σm does not depend on the kurtosis KSSm) is also preferable.

また、雑音成分ｎの抑圧後の尖度ＫSSmが高いほど音響信号ＶOUTにミュージカルノイズが発生し易いという傾向のみに着目すると、尖度ＫSSmが高いほど雑音指標値σmが増加するように指標算定部３２が雑音指標値σmを算定する構成（雑音指標値σmが尖度ＫXmに依存しない構成）も採用される。もっとも、雑音成分ｎの抑圧の前後における尖度の変化度（ＫXm→ＫSSm）が音響信号ＶOUTにおけるミュージカルノイズの発生の度合に最も顕著に反映されるという傾向を考慮すると、尖度ＫXmおよび尖度ＫSSmの双方に応じて雑音指標値σmを算定する構成が好適であり、尖度ＫXmから尖度ＫSSmへの変化度（尖度ＫXmと尖度ＫSSmとの相対比や差分）に応じて雑音指標値σmを算定する構成が格別に好適である。   Focusing only on the tendency that musical noise is more likely to occur in the acoustic signal VOUT as the kurtosis KSSm after suppression of the noise component n is higher, the index calculation unit increases the noise index value σm as the kurtosis KSSm increases. A configuration in which the noise index value σm 32 calculates (a configuration in which the noise index value σm does not depend on the kurtosis KXm) is also employed. However, considering the tendency that the degree of change in kurtosis (KXm → KSSm) before and after suppression of the noise component n is most prominently reflected in the degree of occurrence of musical noise in the acoustic signal VOUT, kurtosis KXm and kurtosis are considered. A configuration that calculates the noise index value σm according to both KSSm is suitable, and the noise index according to the degree of change from the kurtosis KXm to the kurtosis KSSm (the relative ratio or difference between the kurtosis KXm and the kurtosis KSSm) A configuration for calculating the value σm is particularly suitable.

（３）変形例３
以上の各形態においては、雑音指標値σmが相対比ＫRmに対して単調増加する場合を例示したが、相対比ＫRmの増減（尖度ＫXmや尖度ＫSSmの増減）と雑音指標値σmの増減との関係は、雑音指標値σmに応じて雑音抑圧部２６を制御する具体的な方法に応じて適宜に変更される。例えば、式(19)とは逆に雑音指標値σmが大きいほど係数αmが増加するように係数αmを定義した構成においては、相対比ＫRmが増加するほど雑音指標値σmが減少するように相対比ＫRmから雑音指標値σmが算定される。すなわち、相対比ＫRmが増加するほど（尖度ＫXmが減少するほど、または、尖度ＫSSmが増加するほど）、雑音指標値σmが示すミュージカルノイズの発生の度合が増加する構成が本発明では好適に採用され、相対比ＫRmの増加に対して雑音指標値σmが数値として増加するか減少するかは不問である。 (3) Modification 3
In each of the above embodiments, the case where the noise index value σm monotonically increases with respect to the relative ratio KRm is exemplified, but the increase / decrease in the relative ratio KRm (increase / decrease in kurtosis KXm or kurtosis KSSm) and the increase / decrease in the noise index value σm. Is appropriately changed according to a specific method of controlling the noise suppression unit 26 according to the noise index value σm. For example, in the configuration in which the coefficient αm is defined such that the coefficient αm increases as the noise index value σm increases, the noise index value σm decreases as the relative ratio KRm increases. A noise index value σm is calculated from the ratio KRm. That is, a configuration in which the degree of occurrence of musical noise indicated by the noise index value σm increases as the relative ratio KRm increases (the kurtosis KXm decreases or the kurtosis KSSm increases) is preferable in the present invention. No matter whether the noise index value σm increases or decreases as a numerical value with respect to an increase in the relative ratio KRm.

また、本発明が適用される範囲は、雑音指標値σmが示すミュージカルノイズの発生の度合が低いほど雑音成分ｎの抑圧の度合を増加させる構成に限定されない。例えば、音響信号ＶOUTにミュージカルノイズを積極的に発生させる場合（例えば音響信号ＶOUTに発生するミュージカルノイズの特性を調査する場合や、雑音抑圧部２６による処理の良否をミュージカルノイズから判断する場合）には、雑音抑圧度σmが示すミュージカルノイズの発生の度合が高いほど雑音成分ｎの抑圧の度合を増加させるという構成が採用される可能性もある。   The range to which the present invention is applied is not limited to a configuration in which the degree of suppression of the noise component n is increased as the degree of occurrence of musical noise indicated by the noise index value σm is lower. For example, when musical noise is positively generated in the acoustic signal VOUT (for example, when characteristics of the musical noise generated in the acoustic signal VOUT are investigated, or when the quality of the processing by the noise suppression unit 26 is determined from the musical noise). May adopt a configuration in which the degree of suppression of the noise component n increases as the degree of occurrence of musical noise indicated by the degree of noise suppression σm increases.

（４）変形例４
以上の形態においてはフレームＦR毎に雑音指標値σm（尖度ＫXm，尖度ＫSSm，相対比ＫRm，尖度指標値Ｒm）を算定したが、指標算定部３２が雑音指標値σmを算定する周期は任意である。例えば、相前後するフレームＦRで雑音指標値σmが殆ど変化しないという前提に立てば、所定の個数毎に順次に選択したフレームＦRのみについて雑音指標値σmを算定する構成や、複数のフレームＦRにわたる雑音指標値σmの平均値が抑圧制御部３６に指示される構成（あるいは複数のフレームＦRにわたる相対比ＫRmの平均値から雑音指標値σmを算定する構成）も採用される。また、判定部２４２が検出した雑音区間（信号成分ｓが少ない区間）毎に指標算定部３２が雑音指標値σmを算定し、信号成分ｓを含む区間内の各フレームＦRの係数αmの算定に、直前の雑音区間での雑音指標値σmを使用する構成（雑音区間以外では雑音指標値σmを更新しない構成）も好適である。 (4) Modification 4
In the above embodiment, the noise index value σm (kurtosis KXm, kurtosis KSSm, relative ratio KRm, kurtosis index value Rm) is calculated for each frame FR, but the cycle in which the index calculation unit 32 calculates the noise index value σm. Is optional. For example, if it is assumed that the noise index value σm hardly changes in successive frames FR, a configuration in which the noise index value σm is calculated only for the frames FR sequentially selected for each predetermined number, or over a plurality of frames FR. A configuration in which the average value of the noise index value σm is instructed to the suppression control unit 36 (or a configuration in which the noise index value σm is calculated from the average value of the relative ratio KRm over a plurality of frames FR) is also employed. In addition, the index calculation unit 32 calculates the noise index value σm for each noise section (section where the signal component s is small) detected by the determination unit 242, and calculates the coefficient αm of each frame FR in the section including the signal component s. A configuration using the noise index value σm in the immediately preceding noise section (a configuration in which the noise index value σm is not updated outside the noise section) is also suitable.

（５）変形例５
雑音成分ｎの抑圧の前後の尖度（ＫXm，ＫSSm）を算定するための具体的な方法は以上の例示に限定されない。例えば、音響信号ＶINの強度の度数分布を所定の関数（例えば式(6)や式(11)）で近似する構成は本発明において必須ではなく、音響信号ＶIN（周波数スペクトルＸm(ejω)）から直接的に尖度ＫXmを算定する構成や音響信号ＶOUT（周波数スペクトルＹm(ejω)）から直接的に尖度ＫSSmを算定する構成も採用される。 (5) Modification 5
The specific method for calculating the kurtosis (KXm, KSSm) before and after suppression of the noise component n is not limited to the above examples. For example, a configuration for approximating the frequency distribution of the intensity of the acoustic signal VIN with a predetermined function (for example, Equation (6) or Equation (11)) is not essential in the present invention, and the acoustic signal VIN (frequency spectrum Xm (e jω )) A configuration for directly calculating the kurtosis KXm from the sound signal VOUT (frequency spectrum Ym (e jω )) is also employed.

（６）変形例６
第１実施形態においてはミュージカルノイズの発生の度合の目標となる所定値Ｒrefを固定値としたが、所定値Ｒrefを可変値とした構成も好適である。例えば、指標決定部４４は、利用者からの指示（入力装置に対する操作）に応じて所定値Ｒrefを可変に設定する。以上の構成によれば、音響信号ＶOUTにおけるミュージカルノイズの発生の度合（雑音成分ｎの抑圧およびミュージカルノイズの低減の何れを優先させるか）を、例えば利用者の好みに応じて調整できるという利点がある。 (6) Modification 6
In the first embodiment, the predetermined value Rref that is a target of the degree of occurrence of musical noise is a fixed value, but a configuration in which the predetermined value Rref is a variable value is also suitable. For example, the index determination unit 44 variably sets the predetermined value Rref in accordance with an instruction from the user (operation on the input device). According to the above configuration, there is an advantage that the degree of occurrence of musical noise in the acoustic signal VOUT (whether priority is given to suppression of the noise component n or reduction of the musical noise) can be adjusted according to, for example, user's preference. is there.

（７）変形例７
第２実施形態においては、雑音抑圧部２６による実際の処理後の音響信号ＶOUT（周波数スペクトルＹm(ejω)）を尖度ＫSSmの算定に使用したが、第２尖度算定部５２が式(16)の演算で尖度ＫSSmを算定する構成も採用される。式(16)の抑圧係数Ａとしては、過去のフレームＦRについて算定した係数αm（例えば直前のフレームＦRについて算定した係数αm-1）が使用される。以上の態様においては、音響信号ＶOUT（周波数スペクトルＹm(ejω)）が雑音指標値σmの算定に不要であるから、雑音抑圧部２６による処理を待たずに雑音指標値σmを算定できるという利点がある。 (7) Modification 7
In the second embodiment, the acoustic signal VOUT (frequency spectrum Ym (e jω )) actually processed by the noise suppression unit 26 is used for the calculation of the kurtosis KSSm, but the second kurtosis calculation unit 52 uses the formula ( A configuration for calculating the kurtosis KSSm by the calculation of 16) is also adopted. As the suppression coefficient A in equation (16), a coefficient αm calculated for the past frame FR (for example, a coefficient αm−1 calculated for the immediately preceding frame FR) is used. In the above embodiment, since the acoustic signal VOUT (frequency spectrum Ym (e jω )) is not necessary for the calculation of the noise index value σm, the noise index value σm can be calculated without waiting for the processing by the noise suppression unit 26. There is.

（８）変形例８
雑音抑圧部２６による雑音成分ｎの抑圧には公知の技術が任意に採用される。例えば、推定雑音スペクトルψm(ejω)に応じた１未満の係数を周波数スペクトルＸm(ejω)の各周波数の強度に乗算することで雑音成分ｎを抑圧する構成が採用される。抑圧制御部３６は、周波数スペクトルＸm(ejω)に乗算される係数を雑音指標値σmに応じて可変に制御する。 (8) Modification 8
A known technique is arbitrarily employed for the suppression of the noise component n by the noise suppression unit 26. For example, a configuration is adopted in which the noise component n is suppressed by multiplying the intensity of each frequency of the frequency spectrum Xm (e jω ) by a coefficient less than 1 corresponding to the estimated noise spectrum ψm (e jω ). The suppression controller 36 variably controls the coefficient multiplied by the frequency spectrum Xm (e jω ) according to the noise index value σm.

（９）変形例９
以上の形態においては、音響信号ＶINから雑音成分ｎを抑圧する雑音抑圧部２６を具備する雑音抑圧装置１００を例示したが、ミュージカルノイズの発生の度合（雑音成分ｎの抑圧の適否）を評価するための雑音指標値σmや尖度指標値Ｒmの算定に利用される装置（雑音抑圧評価装置）としても本発明は特定される。雑音抑圧評価装置では図１の雑音抑圧部２６や抑圧制御部３６が省略される。また、以上の形態では雑音指標値σmを雑音抑圧部２６の制御に使用したが、雑音抑圧評価装置が算定した雑音指標値σmや尖度指標値Ｒmを使用する方法は任意である（雑音抑圧部２６の制御に限定されない）。例えば、音響信号ＶINの特性（ミュージカルノイズの発生し易さ）を評価するための定量的な指標として雑音指標値σmが使用される。また、雑音抑圧評価装置の算定した雑音指標値σmが可搬型の記録媒体や通信網を介して別個の雑音抑圧装置に供給されたうえで雑音成分ｎの抑圧に使用される。 (9) Modification 9
In the above embodiment, the noise suppression apparatus 100 including the noise suppression unit 26 that suppresses the noise component n from the acoustic signal VIN is exemplified. However, the degree of occurrence of musical noise (appropriateness of suppression of the noise component n) is evaluated. The present invention is also specified as a device (noise suppression evaluation device) used for calculating the noise index value σm and the kurtosis index value Rm. In the noise suppression evaluation apparatus, the noise suppression unit 26 and the suppression control unit 36 in FIG. 1 are omitted. In the above embodiment, the noise index value σm is used for control of the noise suppression unit 26, but the method using the noise index value σm and the kurtosis index value Rm calculated by the noise suppression evaluation apparatus is arbitrary (noise suppression). It is not limited to the control of the unit 26). For example, the noise index value σm is used as a quantitative index for evaluating the characteristics of the acoustic signal VIN (ease of occurrence of musical noise). Further, the noise index value σm calculated by the noise suppression evaluation apparatus is supplied to a separate noise suppression apparatus via a portable recording medium or communication network, and is used for suppression of the noise component n.

本発明の第１実施形態に係る雑音抑圧装置のブロック図である。1 is a block diagram of a noise suppression device according to a first embodiment of the present invention. 音響信号の区分について説明するための概念図である。It is a conceptual diagram for demonstrating the division of an acoustic signal. 音響信号の強度の度数分布が雑音成分の抑圧で変化する様子を示す概念図である。It is a conceptual diagram which shows a mode that the frequency distribution of the intensity | strength of an acoustic signal changes by suppression of a noise component. 指標算定部のブロック図である。It is a block diagram of an index calculation unit. 尖度の相対比が大きい場合（雑音指標値が大きい場合）を例示する概念図である。It is a conceptual diagram which illustrates the case where the relative ratio of kurtosis is large (when the noise index value is large). 尖度の相対比が小さい場合（雑音指標値が小さい場合）を例示する概念図である。It is a conceptual diagram which illustrates the case where the relative ratio of kurtosis is small (when the noise index value is small). 本発明の第２実施形態における指標算定部のブロック図である。It is a block diagram of the parameter | index calculation part in 2nd Embodiment of this invention.

符号の説明Explanation of symbols

１００……雑音抑圧装置、１２……演算処理装置、１４……記憶装置、２２……周波数分析部、２４……雑音推定部、２４２……判定部、２４４……推定部、２６……雑音抑圧部、２８……波形合成部、３２……指標算定部、３４……ＳＮ比算定部、３６……抑圧制御部、４２……相関特定部、４４……指標決定部、５１……第１尖度算定部、５２……第２尖度算定部、５４……算定部、ＶIN，ＶOUT……音響信号、ＫXm……雑音抑圧前の尖度、ＫSSm……雑音抑圧後の尖度、σm……雑音指標値、αm，βm……係数、ψm（ψm(ejω)）……推定雑音スペクトル、Ｘm（Ｘm(ejω)），Ｙm（Ｙm(ejω)）……周波数スペクトル。 DESCRIPTION OF SYMBOLS 100 ... Noise suppression apparatus, 12 ... Arithmetic processing unit, 14 ... Memory | storage device, 22 ... Frequency analysis part, 24 ... Noise estimation part, 242 ... Determination part, 244 ... Estimation part, 26 ... Noise Suppression unit 28... Waveform synthesis unit 32... Index calculation unit 34... SN ratio calculation unit 36 .. Suppression control unit 42 .. Correlation specifying unit 44. 1 kurtosis calculation unit, 52 …… second kurtosis calculation unit, 54 …… calculation unit, VIN, VOUT …… acoustic signal, KXm …… kurtosis before noise suppression, KSSm …… kurtosis after noise suppression, σm: Noise index value, αm, βm: Coefficient, ψm (ψm (e jω )) ... Estimated noise spectrum, Xm (Xm (e jω )), Ym (Ym (e jω )): Frequency spectrum.

Claims (8)

雑音成分の抑圧前または抑圧後の音響信号の強度の度数分布における尖度に応じて変化する雑音指標値を、周波数領域における雑音成分の抑圧後にミュージカルノイズが発生する程度の指標として算定する指標算定手段
を具備する雑音抑圧評価装置。 Index calculation that calculates the noise index value that changes according to the kurtosis in the frequency distribution of the intensity of the acoustic signal before or after suppression of the noise component as an index of the degree to which musical noise occurs after suppression of the noise component in the frequency domain A noise suppression evaluation apparatus comprising: means. 前記指標算定手段は、
雑音成分の抑圧の度合を示す抑圧係数と前記尖度に応じた尖度指標値との関係を特定する相関特定手段と、
前記相関特定手段が特定した関係において前記尖度指標値が所定値に接近するときの前記抑圧係数を前記雑音指標値として決定する指標決定手段とを含む
請求項１の雑音抑圧評価装置。 The index calculation means is:
Correlation specifying means for specifying a relationship between a suppression coefficient indicating a degree of suppression of noise components and a kurtosis index value corresponding to the kurtosis;
The noise suppression evaluation apparatus according to claim 1, further comprising: an index determination unit that determines, as the noise index value, the suppression coefficient when the kurtosis index value approaches a predetermined value in the relationship specified by the correlation specifying unit. 前記指標算定手段は、
雑音成分の抑圧前の音響信号の強度の度数分布における第１尖度を算定する第１尖度算定手段と、
雑音成分の抑圧後の音響信号の強度の度数分布における第２尖度を算定する第２尖度算定手段と、
前記第１尖度と前記第２尖度とから前記雑音指標値を算定する算定手段とを含む
請求項１の雑音抑圧評価装置。 The index calculation means is:
First kurtosis calculation means for calculating the first kurtosis in the frequency distribution of the intensity of the acoustic signal before suppression of the noise component;
A second kurtosis calculating means for calculating a second kurtosis in the frequency distribution of the intensity of the acoustic signal after suppression of the noise component;
The noise suppression evaluation apparatus according to claim 1, further comprising calculation means for calculating the noise index value from the first kurtosis and the second kurtosis. 前記指標算定手段は、雑音成分の抑圧前の音響信号の第１尖度が小さいほど、前記雑音指標値が示すミュージカルノイズの発生の程度が増加するように、前記第１尖度に応じた前記雑音指標値を算定する
請求項１から請求項３の何れかの雑音抑圧評価装置。 The index calculation means is configured to increase the degree of occurrence of musical noise indicated by the noise index value as the first kurtosis of the acoustic signal before suppression of the noise component is smaller. The noise suppression evaluation apparatus according to any one of claims 1 to 3, wherein a noise index value is calculated. 前記指標算定手段は、雑音成分の抑圧後の音響信号の第２尖度が小さいほど、前記雑音指標値が示すミュージカルノイズの発生の程度が減少するように、前記第２尖度に応じた前記雑音指標値を算定する
請求項１から請求項３の何れかの雑音抑圧評価装置。 The index calculation means is configured to adjust the second kurtosis according to the second kurtosis such that the smaller the second kurtosis of the acoustic signal after suppression of the noise component is, the smaller the degree of occurrence of musical noise indicated by the noise index value is. The noise suppression evaluation apparatus according to any one of claims 1 to 3, wherein a noise index value is calculated. 前記指標算定手段は、雑音成分の抑圧前の音響信号の第１尖度に対する抑圧後の音響信号の第２尖度の相対比が大きいほど、前記雑音指標値が示すミュージカルノイズの発生の程度が増加するように、前記第１尖度および前記第２尖度に応じた前記雑音指標値を算定する
請求項１から請求項３の何れかの雑音抑圧評価装置。 As the relative ratio of the second kurtosis of the acoustic signal after suppression to the first kurtosis of the acoustic signal before suppression of the noise component increases, the index calculation means increases the degree of occurrence of musical noise indicated by the noise index value. The noise suppression evaluation apparatus according to any one of claims 1 to 3, wherein the noise index value corresponding to the first kurtosis and the second kurtosis is calculated so as to increase. 前記指標算定手段は、前記第１尖度に対する前記第２尖度の相対比の対数値が大きいほど、前記雑音指標値が示すミュージカルノイズの発生の程度が増加するように、前記対数値に応じた前記雑音指標値を算定する
請求項６の雑音抑圧評価装置。 The index calculation means according to the logarithmic value such that the greater the logarithmic value of the relative ratio of the second kurtosis to the first kurtosis, the greater the degree of occurrence of musical noise indicated by the noise index value. The noise suppression evaluation apparatus according to claim 6, wherein the noise index value is calculated. 雑音成分の抑圧前または抑圧後の音響信号の強度の度数分布における尖度に応じて変化する雑音指標値を、周波数領域における雑音成分の抑圧後にミュージカルノイズが発生する程度の指標として算定する指標算定処理
をコンピュータに実行させるプログラム。
Index calculation that calculates the noise index value that changes according to the kurtosis in the frequency distribution of the intensity of the acoustic signal before or after suppression of the noise component as an index of the degree to which musical noise occurs after suppression of the noise component in the frequency domain A program that causes a computer to execute processing.
JP2008179444A 2008-07-09 2008-07-09 Noise suppression evaluation apparatus and program Expired - Fee Related JP5152800B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2008179444A JP5152800B2 (en) 2008-07-09 2008-07-09 Noise suppression evaluation apparatus and program
EP09164896.4A EP2144233A3 (en) 2008-07-09 2009-07-08 Noise supression estimation device and noise supression device
US12/499,734 US20100008520A1 (en) 2008-07-09 2009-07-08 Noise Suppression Estimation Device and Noise Suppression Device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008179444A JP5152800B2 (en) 2008-07-09 2008-07-09 Noise suppression evaluation apparatus and program

Publications (2)

Publication Number Publication Date
JP2010020013A JP2010020013A (en) 2010-01-28
JP5152800B2 true JP5152800B2 (en) 2013-02-27

Family

ID=41704995

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008179444A Expired - Fee Related JP5152800B2 (en) 2008-07-09 2008-07-09 Noise suppression evaluation apparatus and program

Country Status (1)

Country Link
JP (1) JP5152800B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5609157B2 (en) * 2010-02-26 2014-10-22 ヤマハ株式会社 Coefficient setting device and noise suppression device
JP5633673B2 (en) * 2010-05-31 2014-12-03 ヤマハ株式会社 Noise suppression device and program
JP5728903B2 (en) * 2010-11-26 2015-06-03 ヤマハ株式会社 Sound processing apparatus and program
JP5621637B2 (en) * 2011-02-04 2014-11-12 ヤマハ株式会社 Sound processor

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2771542B1 (en) * 1997-11-21 2000-02-11 Sextant Avionique FREQUENTIAL FILTERING METHOD APPLIED TO NOISE NOISE OF SOUND SIGNALS USING A WIENER FILTER

Also Published As

Publication number Publication date
JP2010020013A (en) 2010-01-28

Similar Documents

Publication Publication Date Title
JP5641186B2 (en) Noise suppression device and program
JP5152799B2 (en) Noise suppression device and program
KR101737824B1 (en) Method and Apparatus for removing a noise signal from input signal in a noisy environment
JP6019969B2 (en) Sound processor
JP5187666B2 (en) Noise suppression device and program
JP2008257049A (en) Noise suppressing device and program
EP2144233A2 (en) Noise supression estimation device and noise supression device
JP5152800B2 (en) Noise suppression evaluation apparatus and program
JP5003419B2 (en) Sound processing apparatus and program
US9418677B2 (en) Noise suppressing device, noise suppressing method, and a non-transitory computer-readable recording medium storing noise suppressing program
JP5915281B2 (en) Sound processor
JP2010220087A (en) Sound processing apparatus and program
US10297272B2 (en) Signal processor
JP5609157B2 (en) Coefficient setting device and noise suppression device
JP5942388B2 (en) Noise suppression coefficient setting device, noise suppression device, and noise suppression coefficient setting method
JP5633673B2 (en) Noise suppression device and program
JP5728903B2 (en) Sound processing apparatus and program
JP5239359B2 (en) Howling suppression device
JP2013246418A (en) Noise suppression device, method, and program
JP3849679B2 (en) Noise removal method, noise removal apparatus, and program
JP2013250356A (en) Coefficient setting device and noise suppression device
JP5772723B2 (en) Acoustic processing apparatus and separation mask generating apparatus
JP2015169901A (en) Acoustic processing device
JP2013182161A (en) Acoustic processing device and program
JP2015169900A (en) Noise suppression device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110629

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20110630

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120730

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120828

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: 20121113

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121128

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20151214

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5152800

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees