JP6613728B2 - Noise suppression device, program and method - Google Patents

Noise suppression device, program and method Download PDF

Info

Publication number
JP6613728B2
JP6613728B2 JP2015169943A JP2015169943A JP6613728B2 JP 6613728 B2 JP6613728 B2 JP 6613728B2 JP 2015169943 A JP2015169943 A JP 2015169943A JP 2015169943 A JP2015169943 A JP 2015169943A JP 6613728 B2 JP6613728 B2 JP 6613728B2
Authority
JP
Japan
Prior art keywords
channel
suppression
noise
channels
input
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2015169943A
Other languages
Japanese (ja)
Other versions
JP2017049277A (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.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry Co Ltd
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 Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP2015169943A priority Critical patent/JP6613728B2/en
Publication of JP2017049277A publication Critical patent/JP2017049277A/en
Application granted granted Critical
Publication of JP6613728B2 publication Critical patent/JP6613728B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

この発明は、雑音抑圧装置、プログラム及び方法に関し、例えば、共通の音響空間の異なった位置に配置された複数のマイクから入力された複数チャネルの入力信号の雑音抑圧に適用し得る。   The present invention relates to a noise suppression device, a program, and a method, and can be applied to noise suppression of input signals of a plurality of channels input from a plurality of microphones arranged at different positions in a common acoustic space, for example.

自然環境において雑音はいたる所に存在するため、一般に実世界で音声を収録すると、観測信号には種々の発信元からの雑音が混入してしまう。それらの雑音は、人が聴くにしても音声の了解性を低下させる。また、音声認識装置や音像定位装置、音源分離装置等の音声処理装置に雑音を含む信号が入力されると、音声処理の性能(例えば音声認識率、音像の定位精度、音源分離音の音質)を低下させる。そのため、入力信号に混入した雑音成分(非目的音成分)を抑圧して音声成分(目的音成分)を強調する技術の需要は高く、これまでに様々な雑音抑圧方法(音声強調方法と呼ばれることもある)が開発されてきた。   Since noise exists everywhere in the natural environment, generally, when recording voice in the real world, noise from various sources is mixed in the observed signal. Those noises reduce the intelligibility of speech even when a person listens. In addition, when a signal including noise is input to a speech processing device such as a speech recognition device, a sound image localization device, or a sound source separation device, the performance of the speech processing (for example, speech recognition rate, accuracy of sound image localization, sound quality of sound source separation sound) Reduce. For this reason, there is a great demand for technology that enhances the speech component (target sound component) by suppressing the noise component (non-target sound component) mixed in the input signal, and so far various noise suppression methods (speech enhancement methods) Has also been developed.

ここで、本明細書で使用する用語を以下のように定義する。以下では、雑音抑圧装置に入力される信号を「入力信号」と呼び、当該入力信号を周波数解析した(周波数領域に変換した)ものを「入力スペクトル」と呼ぶものとする。また、以下では、入力スペクトルの各要素の絶対値をとったものを「入力振幅スペクトル」と呼び、入力スペクトルの各要素のパワー(例えば絶対値の2乗)を算出したものを「入力パワースペクトル」と呼ぶものとする。さらに、以下では、雑音抑圧装置によって雑音を抑圧されたスペクトルを「抑圧後スペクトル」と呼び、抑圧後スペクトルの各要素の絶対値をとったものを「抑圧後振幅スペクトル」と呼び、抑圧後スペクトルの各要素のパワーを算出したものを「抑圧後パワースペクトル」と呼ぶものとする。さらにまた、以下では、抑圧後スペクトルから復元された信号波形を「抑圧後信号」と呼ぶものとする。また、以下では、入力信号に含まれる雑音成分のスペクトルを、「雑音振幅スペクトル」と呼び、雑音振幅スペクトルの各要素のパワーを算出したものを「雑音パワースペクトル」と呼ぶものとする。   Here, terms used in the present specification are defined as follows. Hereinafter, a signal input to the noise suppression device is referred to as an “input signal”, and a frequency analysis (converted into the frequency domain) of the input signal is referred to as an “input spectrum”. In the following, the absolute value of each element of the input spectrum is referred to as “input amplitude spectrum”, and the power of each element of the input spectrum (for example, the square of the absolute value) is calculated as “input power spectrum”. ". Furthermore, in the following, the spectrum in which noise is suppressed by the noise suppression device is referred to as “post-suppression spectrum”, and the absolute value of each element of the post-suppression spectrum is referred to as “post-suppression amplitude spectrum”. The power of each element is calculated as a “suppressed power spectrum”. Furthermore, hereinafter, a signal waveform restored from the post-suppression spectrum is referred to as a “post-suppression signal”. Hereinafter, the spectrum of the noise component included in the input signal is referred to as “noise amplitude spectrum”, and the power of each element of the noise amplitude spectrum is calculated as “noise power spectrum”.

ところで、従来、音像定位や音源分離のような複数チャネルの入力信号を用いる多チャネル信号処理に対して、チャネルごとに雑音抑圧を行う方法には、二つの問題がある。第1の問題点は、独立に算出された抑圧ゲインを各チャネルに適用すると、チャネル間のパワーバランスが崩れ、多チャネル信号処理の性能が低下するという問題である。第2の問題点は、チャネルごとに抑圧ゲインを算出するため、チャネル数の増加に伴って処理量が増加するという問題である。   By the way, conventionally, there are two problems in the method of performing noise suppression for each channel for multi-channel signal processing using input signals of a plurality of channels such as sound image localization and sound source separation. The first problem is that when the independently calculated suppression gain is applied to each channel, the power balance between the channels is lost and the performance of multichannel signal processing is degraded. The second problem is that the amount of processing increases as the number of channels increases because the suppression gain is calculated for each channel.

以上のような問題に対応するための従来技術として特許文献1の記載技術がある。特許文献1の記載技術では、複数チャネルの入力信号から1つの代表信号を生成し、該代表信号に基づいて算出された抑圧ゲインを各チャネルに共通に適用することで、上記の性能低下問題と処理量増加問題を解決しようとしている。特許文献1に記載された装置において、代表信号は、入力パワースペクトルの全チャネルの平均、重み付き加算、最大パワーの選択、最小パワーの選択によって生成される。特許文献1に記載された装置では、このようにして得られた代表信号に基づいて抑圧ゲインは算出される。特許文献1に記載された装置では、抑圧ゲインの算出は一度しか行われないため、チャネルごとに抑圧ゲインを算出するのに比べて処理量を削減できる。特許文献1に記載された装置では、得られた抑圧ゲインは、各チャネルに共通に適用されるので、チャネル間のパワーバランスは保たれる。   As a conventional technique for dealing with the above problems, there is a technique described in Patent Document 1. In the technique described in Patent Document 1, one representative signal is generated from input signals of a plurality of channels, and the suppression gain calculated based on the representative signal is commonly applied to each channel. I'm trying to solve the problem of increased throughput. In the apparatus described in Patent Document 1, the representative signal is generated by averaging all channels of the input power spectrum, weighted addition, selection of maximum power, and selection of minimum power. In the apparatus described in Patent Document 1, the suppression gain is calculated based on the representative signal obtained in this way. In the apparatus described in Patent Document 1, since the suppression gain is calculated only once, the amount of processing can be reduced compared to calculating the suppression gain for each channel. In the apparatus described in Patent Document 1, the obtained suppression gain is commonly applied to each channel, so that the power balance between the channels is maintained.

特許第5435204号公報Japanese Patent No. 5435204

しかしながら、特許文献1に記載された装置では、代表信号を生成するための処理量が、チャネル数に応じてO(n)の処理量オーダ(チャネル数に比例する量の処理)が必要である。例えば、チャネル数をM、入力パワースペクトルの要素数をKとする。このとき、特許文献1に記載された装置を用いて、重み付き加算によって代表信号を得る場合、M*K回の乗算が必要である。また、特許文献1に記載された装置において、最大パワーを選択することによって代表信号を得る場合、(M−1)*K回の比較が必要である。   However, in the apparatus described in Patent Document 1, the processing amount for generating the representative signal requires an O (n) processing amount order (processing in proportion to the number of channels) according to the number of channels. . For example, M is the number of channels and K is the number of elements of the input power spectrum. At this time, when a representative signal is obtained by weighted addition using the device described in Patent Document 1, M * K multiplications are required. Moreover, in the apparatus described in Patent Document 1, when obtaining a representative signal by selecting the maximum power, (M−1) * K comparisons are necessary.

そのため、複数チャネルの入力信号について雑音抑圧処理を行う際の処理を効率的に行うことができる雑音抑圧装置、プログラム及び方法が望まれている。   Therefore, a noise suppression device, a program, and a method that can efficiently perform processing when performing noise suppression processing on a plurality of channels of input signals are desired.

第1の本発明は、複数チャネルの入力信号のそれぞれに対応する入力スペクトルから雑音を抑圧する雑音抑圧装置において、(1)所定の単位時間ごとに、所定の方式により前記複数チャネルから1つのチャネルを選択するチャネル選択手段と、(2)選択されたチャネルの入力スペクトルを用いて当該入力スペクトルに含まれる雑音成分を抑圧するための抑圧ゲインを算出する抑圧ゲイン算出手段と、(3)前記抑圧ゲインを用いて、前記複数チャネルの入力スペクトルのそれぞれについて雑音成分を抑圧する雑音抑圧手段とを有すことを特徴とする。 According to a first aspect of the present invention, there is provided a noise suppression apparatus that suppresses noise from an input spectrum corresponding to each of a plurality of channel input signals. (1) One channel from the plurality of channels according to a predetermined method every predetermined unit time. (2) a suppression gain calculation means for calculating a suppression gain for suppressing a noise component included in the input spectrum using the input spectrum of the selected channel, and (3) the suppression Noise suppression means for suppressing a noise component for each of the input spectra of the plurality of channels using a gain is provided.

第2の本発明の雑音抑圧プログラムは、複数チャネルの入力信号のそれぞれに対応する入力スペクトルから雑音を抑圧する雑音抑圧装置に搭載されたコンピュータを、(1)所定の単位時間ごとに、所定の方式により前記複数チャネルから1つのチャネルを選択するチャネル選択手段と、(2)選択されたチャネルの入力スペクトルを用いて当該入力スペクトルに含まれる雑音成分を抑圧するための抑圧ゲインを算出する抑圧ゲイン算出手段と、(3)前記抑圧ゲインを用いて、前記複数チャネルの入力スペクトルのそれぞれについて雑音成分を抑圧する雑音抑圧手段として機能させることを特徴とする。 A noise suppression program according to a second aspect of the present invention provides a computer mounted on a noise suppression apparatus that suppresses noise from an input spectrum corresponding to each of a plurality of channels of input signals . suppression gain calculating the channel selection means for selecting one channel from the plurality of channels by method a suppression gain for suppressing a noise component contained in the input spectrum by using the input spectrum of the channel selected (2) And (3) functioning as noise suppression means for suppressing noise components for each of the input spectra of the plurality of channels using the suppression gain.

第3の本発明は、複数チャネルの入力信号のそれぞれに対応する入力スペクトルから雑音を抑圧する雑音抑圧方法において、(1)チャネル選択手段、抑圧ゲイン算出手段、及び雑音抑圧手段を備え、(2)前記チャネル選択手段は、所定の単位時間ごとに、所定の方式により前記複数チャネルから1つのチャネルを選択し、(3)前記抑圧ゲイン算出手段は、選択されたチャネルの入力スペクトルを用いて当該入力スペクトルに含まれる雑音成分を抑圧するための抑圧ゲインを算出し、(4)前記雑音抑圧手段は、前記抑圧ゲインを用いて、前記複数チャネルの入力スペクトルのそれぞれについて雑音成分を抑圧することを特徴とする。 According to a third aspect of the present invention, in a noise suppression method for suppressing noise from an input spectrum corresponding to each of a plurality of channel input signals, (1) a channel selection unit, a suppression gain calculation unit, and a noise suppression unit are provided. ) The channel selection means selects one channel from the plurality of channels by a predetermined method every predetermined unit time. (3) The suppression gain calculation means uses the input spectrum of the selected channel to Calculating a suppression gain for suppressing a noise component included in the input spectrum, and (4) the noise suppression means suppresses the noise component for each of the input spectra of the plurality of channels using the suppression gain. Features.

本発明によれば、複数チャネルの入力信号について雑音抑圧処理を行う際の処理を効率的に行う雑音抑圧装置を提供することができる。   ADVANTAGE OF THE INVENTION According to this invention, the noise suppression apparatus which performs efficiently the process at the time of performing a noise suppression process about the input signal of multiple channels can be provided.

第1の実施形態に係る雑音抑圧装置の機能的構成について示したブロック図である。It is the block diagram shown about the functional structure of the noise suppression apparatus which concerns on 1st Embodiment. 第1の実施形態に係る帯域別抑圧手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the suppression means classified by band which concerns on 1st Embodiment. 第1の実施形態に係るチャネル選択手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the channel selection means which concerns on 1st Embodiment. 第1の実施形態に係るチャネル決定手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the channel determination means which concerns on 1st Embodiment. 第1の実施形態に係る抑圧ゲイン算出手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the suppression gain calculation means which concerns on 1st Embodiment. 第2の実施形態に係るチャネル選択手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the channel selection means which concerns on 2nd Embodiment. 第2の実施形態に係るチャネル決定手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the channel determination means which concerns on 2nd Embodiment. 第3の実施形態に係るチャネル選択手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the channel selection means which concerns on 3rd Embodiment. 第3の実施形態に係るチャネル決定手段の内部構成について示したブロック図である。It is the block diagram shown about the internal structure of the channel determination means which concerns on 3rd Embodiment. 第3の実施形態に係る閾値セットのデータ構造について示した図である。It is the figure shown about the data structure of the threshold value set which concerns on 3rd Embodiment. 第3の実施形態に係るチャネル決定手段の動作について示したフロー図である。It is the flowchart shown about operation | movement of the channel determination means which concerns on 3rd Embodiment.

(A)第1の実施形態
以下、本発明による雑音抑圧装置、プログラム及び方法の第1の実施形態を、図面を参照しながら詳述する。
(A−1)第1の実施形態の構成
図1は、この実施形態の雑音抑圧装置100の機能的構成を示すブロック図である。なお、図1において括弧内の符号は後述する第2及び第3の実施形態でのみ用いられる符号である。
(A) First Embodiment Hereinafter, a first embodiment of a noise suppression device, a program, and a method according to the present invention will be described in detail with reference to the drawings.
(A-1) Configuration of First Embodiment FIG. 1 is a block diagram showing a functional configuration of a noise suppression device 100 of this embodiment. In FIG. 1, reference numerals in parentheses are used only in second and third embodiments described later.

雑音抑圧装置100は、M個のチャネルの入力信号is1〜isM(時間領域の信号)の入力を受付ける。以下では、入力信号is1〜isMには、それぞれチャネルを識別するためのチャネル識別子としてチャネル番号が付与されているものとする。例えば、入力信号is1のチャネル番号は1、入力信号is2のチャネル番号は2、…、入力信号isMのチャネル番号はMとなる。入力信号is1〜isMは、例えば、それぞれ異なる図示しないマイクで捕捉された音声信号(音響信号)である。入力信号is1〜isMに係る図示しないマイクは、例えば、同じ室内等の限られた範囲内に配置されているものとする。   The noise suppression apparatus 100 receives input of input signals is1 to isM (time domain signals) of M channels. In the following, it is assumed that a channel number is assigned to each of the input signals is1 to isM as a channel identifier for identifying a channel. For example, the channel number of the input signal is1 is 1, the channel number of the input signal is2 is 2,..., And the channel number of the input signal isM is M. The input signals is1 to isM are audio signals (acoustic signals) captured by different microphones (not shown), for example. Assume that microphones (not shown) related to the input signals is1 to isM are arranged in a limited range such as in the same room.

そして、雑音抑圧装置100は、雑音抑圧の処理後の出力として、入力信号is1〜isのそれぞれに対応する抑圧後信号os1〜osM(時間領域の信号)を出力する。なお、この実施形態では、雑音抑圧装置100は時間領域の信号の入力を受付けて、時間領域の信号の出力を行うものとして説明するが、入出力する信号の形式は限定されないものであり、例えば、周波数領域(スペクトル)の信号入力や、周波数領域の信号出力を行う構成としてもよい。   Then, the noise suppression device 100 outputs post-suppression signals os1 to osM (time domain signals) corresponding to the input signals is1 to is as outputs after the noise suppression processing. In this embodiment, the noise suppression apparatus 100 is described as receiving a time-domain signal and outputting a time-domain signal. However, the format of the input / output signal is not limited. For example, Also, a configuration may be adopted in which signal input in the frequency domain (spectrum) or signal output in the frequency domain is performed.

雑音抑圧装置100は、入力信号is1〜isMのそれぞれに対応する周波数解析手段101(101−1〜101−M)と、帯域別抑圧手段102と、抑圧後信号os1〜osMのそれぞれを出力する波形復元手段103(103−1〜103−M)とを有している。   The noise suppression apparatus 100 outputs frequency analysis means 101 (101-1 to 101-M), band-specific suppression means 102, and post-suppression signals os1 to osM corresponding to the input signals is1 to isM, respectively. And restoring means 103 (103-1 to 103-M).

雑音抑圧装置100は、例えば、プロセッサ及びメモリを有するコンピュータ(プログラムの実施構成)にコンピュータプログラム(実施形態に係る雑音抑圧プログラム)をインストールすることにより構築してもよい。また、雑音抑圧装置100は、一部又は全部をハードウェア(例えば、専用の半導体チップや電気回路等)を用いて実現するようにしてもよい。   The noise suppression device 100 may be constructed, for example, by installing a computer program (a noise suppression program according to the embodiment) in a computer having a processor and a memory (program implementation configuration). The noise suppression apparatus 100 may be realized partially or entirely using hardware (for example, a dedicated semiconductor chip or an electric circuit).

周波数解析手段101−1〜101−Mは、それぞれ入力信号is1〜isMを時間領域の信号から周波数領域の信号に変換する処理を行う。以下では、入力信号is1〜isMを周波数領域に変換した信号(周波数解析手段101−1〜101−Mが出力する信号)を、それぞれ入力スペクトルIS1〜ISMと表すものとする。例えば、周波数解析手段101−1は、入力信号is1を周波数変換して入力スペクトルIS1を生成して出力する。各周波数解析手段101が行う変換処理の具体的方法については限定されないものであるが、例えば高速フーリエ変換(Fast Fourier Transform:FFT)やウェーブレット変換やフィルタパンク等を適用するようにしてもよい。この実施形態では、各周波数解析手段101は、FFTを用いた変換処理を行うものとして説明する。なお、以下では、各入力スペクトルIS1〜ISMは、複素数で表されるデータであるものとして説明する。   The frequency analysis units 101-1 to 101-M perform processing for converting the input signals is1 to isM from time domain signals to frequency domain signals, respectively. Hereinafter, signals obtained by converting the input signals is1 to isM into the frequency domain (signals output from the frequency analysis means 101-1 to 101-M) are represented as input spectra IS1 to ISM, respectively. For example, the frequency analysis unit 101-1 converts the frequency of the input signal is 1 to generate and output the input spectrum IS 1. Although the specific method of the conversion process performed by each frequency analysis unit 101 is not limited, for example, fast Fourier transform (FFT), wavelet transform, filter puncture, or the like may be applied. In this embodiment, each frequency analysis unit 101 will be described as performing a conversion process using FFT. In the following description, each input spectrum IS1 to ISM is assumed to be data represented by complex numbers.

帯域別抑圧手段102は、周波数領域に変換された入力スペクトルIS1〜ISMについて、帯域別に雑音成分を抑圧する処理を行い、雑音成分抑圧後の信号(周波数領域の信号)として抑圧後スペクトルOS1〜OSMを生成する。帯域別抑圧手段102は、抑圧後スペクトルOS1〜OSMを、それぞれ波形復元手段103−1〜103−Mに供給する。帯域別抑圧手段102は、例えば、入力スペクトルIS1について雑音抑圧する処理を行って抑圧後スペクトルOS1を生成し、波形復元手段103−1に供給する。   The band-by-band suppression unit 102 performs a process of suppressing noise components for each band for the input spectra IS1 to ISM converted to the frequency domain, and the post-suppression spectrums OS1 to OSM are obtained as noise component suppressed signals (frequency domain signals). Is generated. The band suppression unit 102 supplies the suppressed spectra OS1 to OSM to the waveform restoration units 103-1 to 103-M, respectively. For example, the band-by-band suppression unit 102 performs a process of noise suppression on the input spectrum IS1 to generate a post-suppression spectrum OS1 and supplies it to the waveform restoration unit 103-1.

波形復元手段103−1〜103−Mは、それぞれ帯域別抑圧手段102から供給された抑圧後スペクトルOS1〜OSMを時間領域の信号に復元し、抑圧後信号os1〜osMとして出力する。例えば、波形復元手段103−1は、抑圧後スペクトルOS1について時間領域に変換して、抑圧後信号os1を出力する。各波形復元手段103は、周波数解析手段101−1〜101−Mと対を成す方法により各抑圧後スペクトルを時間領域の信号に変換する。例えば、周波数解析手段101−1〜101−MでFFTを用いた周波数変換が行われている場合、波形復元手段103−1〜103−Mは、逆高速フーリエ変換(Inverse FFT:IFFT)を用いて、各抑圧後スペクトルを時間領域の信号に復元する。   The waveform restoration means 103-1 to 103 -M restore the post-suppression spectra OS <b> 1 to OSM supplied from the band-by-band suppression means 102, respectively, to time domain signals and output them as post-suppression signals os <b> 1 to osM. For example, the waveform restoration unit 103-1 converts the post-suppression spectrum OS <b> 1 into the time domain and outputs the post-suppression signal os <b> 1. Each waveform restoration unit 103 converts each suppressed spectrum into a signal in the time domain by a method paired with the frequency analysis units 101-1 to 101 -M. For example, when frequency conversion using FFT is performed by the frequency analysis units 101-1 to 101-M, the waveform restoration units 103-1 to 103-M use inverse fast Fourier transform (Inverse FFT: IFFT). Thus, each suppressed spectrum is restored to a time domain signal.

次に、帯域別抑圧手段102の内部構成について図2を用いて説明する。   Next, the internal configuration of the band-by-band suppression unit 102 will be described with reference to FIG.

図2は、帯域別抑圧手段102の内部構成について示したブロック図である。   FIG. 2 is a block diagram showing the internal configuration of the band-by-band suppression unit 102.

図2に示すように、帯域別抑圧手段102は、チャネル選択手段201、抑圧ゲイン算出手段202、及び乗算手段203(203−1〜203−M)を有している。   As shown in FIG. 2, the band-by-band suppression unit 102 includes a channel selection unit 201, a suppression gain calculation unit 202, and a multiplication unit 203 (203-1 to 203-M).

帯域別抑圧手段102に供給された入力スペクトルIS1〜ISMは、すべてチャネル選択手段201に供給されると共に、それぞれの乗算手段203−1〜203−Mにも供給される。例えば、入力スペクトルIS1は、乗算手段203−1に供給されることになる。   All of the input spectra IS1 to ISM supplied to the band-by-band suppression unit 102 are supplied to the channel selection unit 201 and also supplied to the multiplication units 203-1 to 203-M. For example, the input spectrum IS1 is supplied to the multiplication unit 203-1.

チャネル選択手段201は、入力スペクトルIS1〜ISMのうちいずれかを、所定の方式により選択して出力する処理を行う。以下では、チャネル選択手段201で選択された入力スペクトル(入力スペクトルIS1〜ISMのうち選択されたもの)を選択入力スペクトルと呼ぶものとする。チャネル選択手段201は、選択入力スペクトルを、抑圧ゲイン算出手段202に供給する。   The channel selection means 201 performs a process of selecting and outputting any one of the input spectra IS1 to ISM by a predetermined method. Hereinafter, the input spectrum selected by the channel selection unit 201 (selected one of the input spectra IS1 to ISM) is referred to as a selected input spectrum. The channel selection unit 201 supplies the selected input spectrum to the suppression gain calculation unit 202.

チャネル選択手段201は、所定の順序で選択入力スペクトルとするチャネル番号を切替える処理を行う。例えば、チャネル選択手段201は、1単位時間(例えば、抑圧ゲイン算出手段202の処理単位時間)ごとに、選択入力スペクトルとするチャネル番号を切替える。このとき、チャネル選択手段201がチャネル番号を切替える順序は限定されないものである。この実施形態では、チャネル選択手段201は、チャネル番号順に循環的にチャネル番号を選択(選択入力スペクトルとするチャネル番号を選択)するものとする。言い換えると、チャネル選択手段201は、チャネル番号順にラウンドロビンで選択入力スペクトルとするチャネル番号を選択するものとする。すなわち、この実施形態のチャネル選択手段201は、「1、2、…、M、1、2、…M、1、…」という順序で選択入力スペクトルとするチャネル番号を選択するものとする。   The channel selection unit 201 performs processing for switching the channel number to be the selected input spectrum in a predetermined order. For example, the channel selection unit 201 switches the channel number to be the selected input spectrum every unit time (for example, the processing unit time of the suppression gain calculation unit 202). At this time, the order in which the channel selection unit 201 switches the channel numbers is not limited. In this embodiment, the channel selection means 201 shall select a channel number cyclically in order of a channel number (select a channel number used as a selection input spectrum). In other words, it is assumed that the channel selection unit 201 selects a channel number as a selection input spectrum in round-robin order in the order of channel numbers. That is, it is assumed that the channel selection unit 201 of this embodiment selects channel numbers to be selected input spectra in the order of “1, 2,..., M, 1, 2,.

なお、チャネル選択手段201が、選択入力スペクトルとするチャネル番号を選択する順序のアルゴリズムについては限定されないものである。例えば、チャネル番号順に循環的に選択入力スペクトルとするチャネル番号を選択する例以外にも、後述する第2の実施形態のように予め設定された順序としたり、後述する第3の実施形態のように疑似乱数に基づく順序(ランダムな順序)で選択するようにしてもよい。   Note that the algorithm of the order in which the channel selection unit 201 selects the channel number to be the selected input spectrum is not limited. For example, in addition to the example of selecting a channel number as a cyclically selected input spectrum in the order of channel numbers, the order may be set in advance as in the second embodiment described later, or as in the third embodiment described later. May be selected in an order based on pseudo-random numbers (random order).

抑圧ゲイン算出手段202は、選択入力スペクトルに基づいて、雑音成分のスペクトルを推定して取得する。以下では、抑圧ゲイン算出手段202が推定した雑音スペクトルを「抑圧ゲイン」とも呼ぶものとする。抑圧ゲイン算出手段202が抑圧ゲインを算出する処理の詳細(詳細構成)については後述する。抑圧ゲイン算出手段202は、算出した抑圧ゲインを、それぞれの乗算手段203−1〜203−Mに供給する。   The suppression gain calculation unit 202 estimates and acquires the noise component spectrum based on the selected input spectrum. Hereinafter, the noise spectrum estimated by the suppression gain calculation unit 202 is also referred to as “suppression gain”. Details (detailed configuration) of the processing in which the suppression gain calculation unit 202 calculates the suppression gain will be described later. The suppression gain calculation unit 202 supplies the calculated suppression gain to each multiplication unit 203-1 to 203-M.

乗算手段203−1〜203−Mは、それぞれ入力スペクトルIS1〜ISMに対して共通の抑圧ゲイン(抑圧ゲイン算出手段202から供給された抑圧ゲイン)を乗じる処理を行う。これにより、乗算手段203−1〜203−Mは、それぞれ入力スペクトルIS1〜ISMに含まれる雑音成分を抑圧して抑圧後スペクトルOS1〜OSMを得る。   Multiplication means 203-1 to 203-M perform a process of multiplying the input spectra IS1 to ISM by a common suppression gain (a suppression gain supplied from the suppression gain calculation means 202), respectively. Thereby, the multiplication means 203-1 to 203-M suppress the noise components included in the input spectra IS1 to ISM, respectively, and obtain post-suppression spectra OS1 to OSM.

そして、乗算手段203−1〜203−Mは、それぞれ抑圧後スペクトルOS1〜OSMを出力する。例えば、乗算手段203−1は、入力スペクトルIS1に抑圧ゲインを乗じて、抑圧後スペクトルOS1を取得して出力する処理を行うことになる。なお、抑圧後スペクトルOS1〜OSMは、それぞれ波形復元手段103−1〜103−Mに供給されることになる。すなわち、乗算手段203−1〜203−Mは、入力スペクトルIS1〜ISMのそれぞれについて共通する抑圧ゲインを用いた雑音抑圧を行う雑音抑圧手段として機能することになる。   Then, the multiplication means 203-1 to 203-M output the suppressed spectra OS1 to OSM, respectively. For example, the multiplying unit 203-1 performs a process of multiplying the input spectrum IS1 by a suppression gain and acquiring and outputting the suppressed spectrum OS1. Note that the suppressed spectra OS1 to OSM are supplied to the waveform restoration units 103-1 to 103 -M, respectively. That is, the multiplication means 203-1 to 203-M function as noise suppression means for performing noise suppression using a suppression gain common to each of the input spectra IS1 to ISM.

次に、チャネル選択手段201の内部構成について図3を用いて説明する。   Next, the internal configuration of the channel selection unit 201 will be described with reference to FIG.

図3は、チャネル選択手段201の内部構成について示したブロック図である。   FIG. 3 is a block diagram showing the internal configuration of the channel selection means 201.

図3に示すように、チャネル選択手段201は、フレーム計数手段301、チャネル決定手段302、チャネル記憶手段303、及び入力スペクトル選択手段304を有している。   As shown in FIG. 3, the channel selection unit 201 includes a frame counting unit 301, a channel determination unit 302, a channel storage unit 303, and an input spectrum selection unit 304.

フレーム計数手段301は、所定の単位時間(例えば、抑圧ゲイン算出手段202の処理単位時間)ごとに処理フレーム番号を1ずつインクリメントし、得られたフレーム番号をチャネル決定手段302に与える。   The frame counting unit 301 increments the processing frame number by one for every predetermined unit time (for example, the processing unit time of the suppression gain calculation unit 202), and gives the obtained frame number to the channel determination unit 302.

チャネル決定手段302は、与えられたフレーム番号と、チャネル記憶手段303から供給されたチャネル番号(1単位時間前のチャネル番号)とに基づいて、次のチャネル番号(次に、選択入力スペクトルとして選択するチャネル番号)を決定する。そして、チャネル決定手段302は、決定したチャネル番号を選択チャネル番号として入力スペクトル選択手段304に与える。チャネル決定手段302の詳細構成については後述する。   The channel determination unit 302 selects the next channel number (next as the selected input spectrum) based on the given frame number and the channel number (channel number one unit time before) supplied from the channel storage unit 303. Channel number) to be determined. Then, the channel determination unit 302 gives the determined channel number to the input spectrum selection unit 304 as the selected channel number. The detailed configuration of the channel determination unit 302 will be described later.

チャネル記憶手段303は、チャネル決定手段302から与えられたチャネル番号を記憶し、記憶していたチャネル番号を1単位時間後にチャネル決定手段302に与える。   The channel storage means 303 stores the channel number given from the channel determination means 302 and gives the stored channel number to the channel determination means 302 after one unit time.

入力スペクトル選択手段304は、チャネル決定手段302から与えられたチャネル番号に対応する入力スペクトルを選択入力スペクトルとして選択する。例えば、チャネル決定手段302から与えられたチャネル番号がmだった場合、選択入力スペクトルはチャネルmの入力スベクトル(入力スペクトルISm)を選択入力スペクトルとして取り扱う。   The input spectrum selection unit 304 selects an input spectrum corresponding to the channel number given from the channel determination unit 302 as a selected input spectrum. For example, when the channel number given from the channel determination means 302 is m, the selected input spectrum handles the input vector of the channel m (input spectrum ISm) as the selected input spectrum.

次に、チャネル決定手段302の内部構成について、図4を用いて説明する。   Next, the internal configuration of the channel determination unit 302 will be described with reference to FIG.

図4は、チャネル決定手段302の内部構成について示したブロック図である。   FIG. 4 is a block diagram showing the internal configuration of the channel determination unit 302.

図4に示すように、チャネル決定手段302は、チャネル計数手段401、定数1値供給手段402、及び初期チャネルスイッチ手段403を有している。   As shown in FIG. 4, the channel determination unit 302 includes a channel counting unit 401, a constant 1 value supply unit 402, and an initial channel switch unit 403.

チャネル決定手段302に与えられたフレーム番号は初期チャネルスイッチ手段403に与えられる。また、チャネル決定手段302に与えられた1単位時間前のチャネル番号はチャネル計数手段401に与えられる。   The frame number given to the channel determining means 302 is given to the initial channel switching means 403. Further, the channel number one unit time before given to the channel determining means 302 is given to the channel counting means 401.

チャネル計数手段401は、与えられた1単位時間前のチャネル番号を更新し、得られた暫定選択チャネル番号を初期チャネルスイッチ手段403に与える。チャネル計数手段401におけるチャネル番号の更新は、以下のように行われる。単位時間前のチャネル番号が1以上かつ(M−1)以下である場合には、チャネル計数手段401は、当該単位時間前のチャネル番号を1だけインクリメントした値を暫定選択チャネル番号とする。そして、チャネル計数手段401は、単位時間前のチャネル番号がMである場合には、暫定選択チャネル番号を初期値の1とする。   The channel counting means 401 updates the given channel number one unit time ago, and gives the obtained provisional selection channel number to the initial channel switching means 403. The channel number is updated in the channel counting means 401 as follows. When the channel number before the unit time is 1 or more and (M−1) or less, the channel counting means 401 sets a value obtained by incrementing the channel number before the unit time by 1 as the provisional selection channel number. Then, when the channel number before the unit time is M, the channel counting unit 401 sets the temporary selection channel number to 1 as an initial value.

定数1値供給手段402は、常に値「1」を初期チャネルスイッチ手段403に与える。値「1」は、チャネル番号の最初の値(初期値)である。なお、第1の実施形態では初期値を値「1」としたが、チャネル番号を表す1以上M以下の整数であればどの値を初期値として適用しても良い。   The constant 1 value supply means 402 always gives the value “1” to the initial channel switch means 403. The value “1” is the first value (initial value) of the channel number. In the first embodiment, the initial value is the value “1”, but any value may be applied as the initial value as long as it is an integer of 1 or more and M or less representing the channel number.

初期チャネルスイッチ手段403は、フレーム番号が初期フレーム番号であれば定数1値供給手段402から与えられた値「1」を選択チャネル番号とし、そうでなければチャネル計数手段から与えられた暫定選択チャネル番号を選択チャネル番号とする。初期フレーム番号は実装によって異なり、例えば0又は1とするようにしてもよい。また、初期動作を安定化させるために、初期チャネルスイッチ手段403は、フレーム番号が小さな整数値(例えば4)以下であるときに定数1値供給手段402から与えられた値「1」を選択チャネル番号としても良い。   If the frame number is the initial frame number, the initial channel switch means 403 uses the value “1” given from the constant 1 value supply means 402 as the selected channel number, otherwise the provisional selection channel given from the channel counting means. The number is the selected channel number. The initial frame number varies depending on the implementation, and may be set to 0 or 1, for example. In order to stabilize the initial operation, the initial channel switch means 403 selects the value “1” given from the constant 1 value supply means 402 when the frame number is a small integer value (for example, 4) or less as the selected channel. It is good also as a number.

次に、抑圧ゲイン算出手段202の内部構成について、図5を用いて説明する。   Next, the internal configuration of the suppression gain calculation unit 202 will be described with reference to FIG.

図5は、抑圧ゲイン算出手段202の内部構成について示したブロック図である。   FIG. 5 is a block diagram showing the internal configuration of the suppression gain calculation unit 202.

この実施形態では、抑圧ゲイン算出手段202は、事後SNRに基づいて判定される音声区間に基づいて推定された雑音パワースペクトルを用い、更にDecision−Directed法によって推定した事前SNRに基づくウィナーフィルタ(WF)で抑圧ゲインを算出(推定)するものとして説明する。なお、抑圧ゲイン算出手段202において抑圧ゲインを算出(推定)する具体的な算出方法は限定されないものである。   In this embodiment, the suppression gain calculation means 202 uses the noise power spectrum estimated based on the speech section determined based on the posterior SNR, and further uses the Wiener filter (WF) based on the a priori SNR estimated by the Decision-Directed method. ) To calculate (estimate) the suppression gain. A specific calculation method for calculating (estimating) the suppression gain in the suppression gain calculation unit 202 is not limited.

例えば、抑圧ゲイン算出手段202は、雑音パワースペクトルの推定には重み付き雑音推定(参考文献1「IEICE Trans.on Fundam.,Vol.E85−A、No.7,PP.1710―1718,Jul.2002」参照)や、MAP推定に基づく方法(参考文献2「電子情報通信学会技術報告,Vol.113,No.503,EA2013−121,PP.7−12,2014年3月」参照)等を用いた抑圧ゲインの算出(推定)処理を行うようにしてもよい、また、抑圧ゲイン算出手段202は、スペクトル減算法(参考文献3「IEEE Trans.,Vol.ASSP−27,No.2,pp.113−120,Apr.1979」)や、MMSESTSA法(参考文献4「IEEE Trans.,Vol.ASSP−32,No.6,PP.1109−1121,Dec.1984」)を用いた抑圧ゲインの算出(推定)処理を行うようにしてもよい。   For example, the suppression gain calculation means 202 uses weighted noise estimation (Reference 1 “IEICE Trans. On Fundam., Vol. E85-A, No. 7, PP. 1710-1718, Jul. 2002 ”) and a method based on MAP estimation (see Reference 2“ Technical Report of the Institute of Electronics, Information and Communication Engineers, Vol. 113, No. 503, EA 2013-121, PP. 7-12, March 2014 ”), etc. The used suppression gain calculation (estimation) processing may be performed, and the suppression gain calculation means 202 may be a spectral subtraction method (reference document 3 “IEEE Trans., Vol. ASSP-27, No. 2, pp. 113-120, Apr. 1979 ") and the MMSESTSA method (reference document 4" IEEE Trans., Vol. ASSP-32, No. 6, PP. 1109-1121, Dec. 1984 ") may be used to perform suppression gain calculation (estimation) processing.

図5に示すように、抑圧ゲイン算出手段202は、パワー算出手段501、雑音平滑化手段502、雑音記憶手段503、事後SNR算出手段504、事後SNR記憶手段505、閾値判定手段506、閾値供給手段507、事前SNR推定手段508、抑圧ゲイン決定手段509、及び抑圧ゲイン記憶手段510を有している。   As shown in FIG. 5, the suppression gain calculation unit 202 includes a power calculation unit 501, a noise smoothing unit 502, a noise storage unit 503, a post-SNR calculation unit 504, a post-SNR storage unit 505, a threshold determination unit 506, and a threshold supply unit. 507, an a priori SNR estimation unit 508, a suppression gain determination unit 509, and a suppression gain storage unit 510.

抑圧ゲイン算出手段202は、与えられた選択入力スペクトルをパワー算出手段501に与える。パワー算出手段501は、与えられた選択入力スペクトルの各要素の絶対値の2乗を要素ごとに計算し、得られた入力パワースペクトルを雑音平滑化手段502と事後SNR算出手段に与える。   The suppression gain calculation unit 202 gives the given selection input spectrum to the power calculation unit 501. The power calculation unit 501 calculates the square of the absolute value of each element of the given selected input spectrum for each element, and gives the obtained input power spectrum to the noise smoothing unit 502 and the post SNR calculation unit.

雑音平滑化手段502は、雑音記憶手段503から与えられる単位時間前の雑音パワースペクトルN’(k)と閾値判定手段506から与えられる音声区間判定結果Vに基づいて以下の式(1)によって雑音パワーを推定し、得られた雑音パワースペクトルN(k)を雑音記憶手段503と事後SNR算出手段504に与える。以下の(1)式において、k(k=0,...,K−1)は要素(周波数帯域)のインデックス、τは時定数、X(k)は入力パワースペクトルを示している。また、以下の(1)式では、時定数τは0.8程度とすることが好適である。

Figure 0006613728
Based on the noise power spectrum N ′ (k) before unit time given from the noise storage means 503 and the speech section judgment result V given from the threshold judgment means 506, the noise smoothing means 502 performs noise according to the following equation (1). The power is estimated, and the obtained noise power spectrum N (k) is given to the noise storage means 503 and the posterior SNR calculation means 504. In the following equation (1), k (k = 0,..., K−1) is an index of an element (frequency band), τ is a time constant, and X (k) is an input power spectrum. In the following equation (1), the time constant τ is preferably about 0.8.
Figure 0006613728

雑音記憶手段503は、与えられた雑音パワースペクトルを記憶し、記憶していた雑音パワースペクトルを1単位時間後に雑音平滑化手段502に与える。   The noise storage means 503 stores the given noise power spectrum, and gives the stored noise power spectrum to the noise smoothing means 502 after one unit time.

事後SNR算出手段504は、要素ごとに与えられた入力パワースペクトルを与えられた雑音パワースペクトルで除し、得られた事後SNRを事後SNR記憶手段505及び事前SNR推定手段508に与える。   The posterior SNR calculation means 504 divides the input power spectrum given for each element by the given noise power spectrum, and gives the obtained posterior SNR to the posterior SNR storage means 505 and the prior SNR estimation means 508.

事後SNR記憶手段505は、与えられた事後SNRを記憶し、記憶していた事後SNRを1単位時間後に閾値判定手段506及び事前SNR推定手段508に与える。   The posterior SNR storage unit 505 stores the given posterior SNR, and provides the stored posterior SNR to the threshold determination unit 506 and the prior SNR estimation unit 508 after one unit time.

閾値判定手段506は、事後SNR記憶手段505から与えられた単位時間前の事後SNRが、閾値供給手段507から与えられた事後SNRの閾値(詳細については後述する)より大きければ音声区間判定結果VにTrue(真値)をセットし、そうでなければVにFalse(偽値)をセットする。そして、閾値判定手段506は、得られた音声区間判定結果Vを雑音平滑化手段502に与える。   The threshold determination unit 506 determines the voice segment determination result V if the posterior SNR before the unit time given from the posterior SNR storage unit 505 is larger than the posterior SNR threshold given from the threshold supply unit 507 (details will be described later). True (true value) is set to V. Otherwise, False (false value) is set to V. Then, the threshold determination unit 506 gives the obtained speech segment determination result V to the noise smoothing unit 502.

閾値供給手段507は、事前に定められた事後SNR閾値を閾値判定手段506に与える。事後SNR閾値と比較される事後SNRは、1より大きければ入力パワーが雑音パワーより大きく、1より小さければ入力パワーが雑音パワーより小さいことを表しているが、入力パワーは単位時間ごとに細かく変化するのに対して雑音パワーはそれほど変化しない。そこで、入力パワーが雑音パワーより多少大きくても雑音平滑化手段502が雑音パワーを更新するようにするために、事後SNR閾値には2.0程度の値とすることが好適である。   The threshold supply unit 507 gives a predetermined posterior SNR threshold to the threshold determination unit 506. The posterior SNR compared with the posterior SNR threshold indicates that the input power is greater than the noise power if it is greater than 1, and the input power is smaller than the noise power if it is less than 1, the input power varies finely per unit time. However, the noise power does not change that much. Therefore, in order for the noise smoothing means 502 to update the noise power even if the input power is slightly larger than the noise power, it is preferable to set the posterior SNR threshold to a value of about 2.0.

事前SNR推定手段508は、与えられた事後SNR、単位時間前の事後SNR及び後述する単位時間前の抑圧ゲインに基づいて事前SNRを算出(推定)し、得られた事前SNRを抑圧ゲイン決定手段509に与える。事後SNRをγ(k)、単位時間前の事後SNRをγ’(k)、単位時間前の抑圧ゲインをG’(k)、事前SNRをξ(k)とすると、事前SNRの推定はDecision−Directed法と呼ばれる以下の(2)式によって行うことができる。以下の(2)式において、αは忘却係数であり、0.99程度の値とすることが好適である。   The a priori SNR estimation unit 508 calculates (estimates) the a priori SNR based on the given a posteriori SNR, the a posteriori SNR before unit time, and a suppression gain before unit time described later, and uses the obtained a priori SNR as a suppression gain determination unit. 509. Assuming that the posterior SNR is γ (k), the posterior SNR before unit time is γ ′ (k), the suppression gain before unit time is G ′ (k), and the prior SNR is ξ (k), the prior SNR is estimated as Decision It can be performed by the following equation (2) called the -Directed method. In the following formula (2), α is a forgetting factor, and is preferably set to a value of about 0.99.

抑圧ゲイン決定手段509は、事前SNRξ(k)に基づいて以下の(3)式によって抑圧ゲインG(k)を算出し、得られた抑圧ゲインを抑圧ゲイン記憶手段510に与える。当該抑圧ゲインは抑圧ゲイン算出手段202の出力として出力される。

Figure 0006613728
The suppression gain determination unit 509 calculates the suppression gain G (k) by the following equation (3) based on the prior SNRξ (k), and gives the obtained suppression gain to the suppression gain storage unit 510. The suppression gain is output as an output of the suppression gain calculation unit 202.
Figure 0006613728

抑圧ゲイン記憶手段510は、与えられた抑圧ゲインG(k)を記憶し、記憶した抑圧ゲインを1単位時間後に事前SNR推定手段508に与える。   The suppression gain storage unit 510 stores the given suppression gain G (k), and provides the stored suppression gain to the prior SNR estimation unit 508 after one unit time.

(A−2)第1の実施形態の動作
次に、以上のような構成を有する第1の実施形態の雑音抑圧装置100の動作(この実施形態の雑音抑圧方法)について図1〜図4を用いて説明する。
(A-2) Operation of the First Embodiment Next, the operation of the noise suppression device 100 of the first embodiment having the above configuration (the noise suppression method of this embodiment) will be described with reference to FIGS. It explains using.

周波数解析手段101−1〜101−Mは、それぞれ入力信号is1〜isMを時間領域の信号から周波数領域の信号に変換する処理を行う。   The frequency analysis units 101-1 to 101-M perform processing for converting the input signals is1 to isM from time domain signals to frequency domain signals, respectively.

そして、帯域別抑圧手段102は、周波数領域に変換された入力スペクトルIS1〜ISMについて、帯域別に雑音成分抑圧の処理を行い、雑音成分抑圧後の信号(周波数領域の信号)として抑圧後スペクトルOS1〜OSMを生成する。   Then, the band-by-band suppression unit 102 performs noise component suppression processing for each band on the input spectra IS1 to ISM converted into the frequency domain, and the post-suppression spectrum OS1 to OS1 as signals after noise component suppression (frequency domain signals). Create an OSM.

帯域別抑圧手段102では、チャネル選択手段201により、チャネル番号順に循環的に選択入力スペクトルに適用するチャネル番号が選択される。そして、抑圧ゲイン算出手段202では、チャネル選択手段201により選択された選択入力スペクトルに基づいて抑圧ゲインが生成される。そして、帯域別抑圧手段102では、乗算手段203−1〜203−Mにより、入力スペクトルIS1〜ISMに対して共通の抑圧ゲイン(抑圧ゲイン算出手段202から供給された抑圧ゲイン)を乗じる処理が行われる。これにより、乗算手段203−1〜203−Mは、それぞれ入力スペクトルIS1〜ISMに含まれる雑音成分を抑圧して抑圧後スペクトルOS1〜OSMを出力する。   In the band-by-band suppression unit 102, the channel selection unit 201 selects channel numbers to be cyclically applied to the selected input spectrum in the order of channel numbers. The suppression gain calculation unit 202 generates a suppression gain based on the selected input spectrum selected by the channel selection unit 201. In the band-by-band suppression unit 102, the multiplication units 203-1 to 203-M perform a process of multiplying the input spectra IS1 to ISM by a common suppression gain (a suppression gain supplied from the suppression gain calculation unit 202). Is called. Thereby, the multiplication means 203-1 to 203-M suppress the noise components included in the input spectra IS1 to ISM, respectively, and output the suppressed spectra OS1 to OSM.

そして、波形復元手段103−1〜103−Mは、それぞれ帯域別抑圧手段102から供給された抑圧後スペクトルOS1〜OSMを時間領域の信号に変換し、抑圧後信号os1〜osMとして出力する。   Then, the waveform restoration units 103-1 to 103 -M convert the suppressed spectra OS <b> 1 to OSM supplied from the band-based suppression units 102 into time-domain signals, and output them as post-suppression signals os <b> 1 to osM.

(A−3)第1の実施形態の効果
第1の実施形態によれば、以下のような効果を奏することができる。
(A-3) Effects of First Embodiment According to the first embodiment, the following effects can be achieved.

第1の実施形態の雑音抑圧装置100では、入力信号のチャネル数が増えた場合でも、チャネル選択手段201で、選択入力スペクトルに適用するチャネル番号の選択処理の量は一定となる。したがって、第1の実施形態の雑音抑圧装置100では、入力信号のチャネル数が増えた場合でも、チャネル数に応じてO(n)のオーダとはならず、従来よりも少ない処理量で雑音抑制を行うことができる。   In the noise suppression apparatus 100 of the first embodiment, the amount of channel number selection processing applied to the selected input spectrum by the channel selection unit 201 is constant even when the number of channels of the input signal increases. Therefore, in the noise suppression apparatus 100 according to the first embodiment, even when the number of channels of the input signal is increased, the order of O (n) is not increased according to the number of channels, and noise suppression is performed with a smaller processing amount than in the past. It can be performed.

また、第1の実施形態の雑音抑圧装置100では、任意の選択方法(例えば、チャネル番号順に循環的に選択する方法)で選択入力スペクトルを切替えながら抑圧ゲインを生成することで、実質的に複数チャネルの入力パワースペクトルを平均的に反映した抑圧ゲインを生成して雑音抑圧を行うことが可能となる。   Further, in the noise suppression device 100 of the first embodiment, a plurality of suppression gains are generated by switching the selection input spectrum by an arbitrary selection method (for example, a method of cyclically selecting in the order of channel numbers). Noise suppression can be performed by generating a suppression gain that reflects the average input power spectrum of the channel.

以下に、雑音抑圧装置100において、任意の選択方法(例えば、チャネル番号順に循環的に選択する方法)で選択入力スペクトルを切替えながら抑圧ゲインを生成する場合でも、雑音抑圧性能の維持が可能となることの数学的な裏付けについて説明する。   Hereinafter, in the noise suppression apparatus 100, it is possible to maintain the noise suppression performance even when generating the suppression gain while switching the selected input spectrum by an arbitrary selection method (for example, a method of cyclically selecting in order of channel numbers). Explain the mathematical support of this.

チャネル選択手段201において、選択入力スペクトルに適用するチャネル番号を選択する際の選択方法としては、例えば、以下の3つの選択方法が挙げられる。   As the selection method for selecting a channel number to be applied to the selected input spectrum in the channel selection means 201, for example, the following three selection methods can be mentioned.

第1の選択方法は、チャネル選択手段201が、抑圧ゲイン算出手段202が抑圧ゲイン算出を実行する単位時間ごとに、チャネル番号順に循環的にチャネル番号を選択する方法(上述の第1の実施形態と同様の選択方法)である。第1の選択方法では、チャネル選択手段201は、上述の第1の実施形態と同様に、例えば、チャネル1、チャネル2、・・・と選択していき、最後のチャネルMが選択された次の単位時間後は再びチャネル1を選択する。   The first selection method is a method in which the channel selection unit 201 cyclically selects channel numbers in order of channel number for each unit time for which the suppression gain calculation unit 202 executes suppression gain calculation (the first embodiment described above). Is the same selection method). In the first selection method, the channel selection unit 201 selects, for example, channel 1, channel 2,..., And the next channel M is selected, as in the first embodiment. After the unit time, channel 1 is selected again.

第2の選択方法は、チャネル選択手段201が、抑圧ゲイン算出手段202が抑圧ゲイン算出を実行する単位時間ごとに、予め定められた順序(例えば、予め設定されたリストに基づく順序)で、チャネル番号を選択する方法(後述する第2の実施形態と同様の選択方法)である。   In the second selection method, the channel selection unit 201 uses a predetermined order (for example, an order based on a preset list) for each unit time for which the suppression gain calculation unit 202 executes the suppression gain calculation. This is a method of selecting a number (the same selection method as in the second embodiment described later).

第3の選択方法は、チャネル選択手段201が、疑似乱数に基づく順序(ランダムな順序)でチャネル番号を選択する方法(後述する第3の実施形態と同様の選択方法)である。具体的には、チャネル選択手段201は、疑似乱数を発生させて、事前に定めておいたチャネル数より1つだけ少ない個数の閾値と該疑似乱数とを比較し、当該疑似乱数を下回る閾値の数に対応するチャネル番号を選択するようにしてもよい。   The third selection method is a method in which the channel selection unit 201 selects channel numbers in an order based on pseudo-random numbers (random order) (selection method similar to the third embodiment described later). Specifically, the channel selection unit 201 generates a pseudo-random number, compares the pseudo-random number with a threshold value that is one less than the predetermined number of channels, and sets a threshold value lower than the pseudo-random number. A channel number corresponding to the number may be selected.

チャネル選択手段201が、上述の第1〜第3の選択方法で、チャネル番号1〜M(M=4)の選択を行った場合、その処理量オーダは、第1と第2の選択方法ではO(1)であり、第3の選択方法ではO(logn)あるから、雑音抑圧装置100において抑圧ゲインを算出するための処理オーダもO(1)またはO(logn)となる。すなわち、この場合、雑音抑圧装置100では入力信号のチャネル数に左右されない処理量での抑圧ゲイン算出を実現できる。 When the channel selection unit 201 selects channel numbers 1 to M (M = 4) by the above-described first to third selection methods, the processing amount order is the first and second selection methods. Since O (1) and O (log 2 n) in the third selection method, the processing order for calculating the suppression gain in the noise suppression apparatus 100 is also O (1) or O (log 2 n). . That is, in this case, the noise suppression apparatus 100 can realize suppression gain calculation with a processing amount that does not depend on the number of channels of the input signal.

上述の第1の選択方法と第2の選択方法は、長期的(例えば、数秒以上)に見れば確率的にチャネルを選択しているのと同じであるから、上述の第3の選択方法の特別な場合と見なすことができる。そこで、以下では、上述の第3の選択方法の例を用いて、雑音抑圧装置100において抑圧性能が維持できることの理由説明を行う。   Since the first selection method and the second selection method described above are the same as selecting a channel stochastically in the long term (for example, several seconds or more), the third selection method described above It can be regarded as a special case. Therefore, in the following, the reason why the suppression performance can be maintained in the noise suppression device 100 will be described using the example of the third selection method described above.

抑圧ゲイン算出手段202で、抑圧ゲインを算出するためには、先に雑音パワースペクトルを推定しなければいけない。入力信号が単一チャネルの場合、一般に、時刻tにおける雑音パワースペクトルのk番目の要素(「k番目の帯域」と同義)N(k,t)は、入力パワースペクトルの要素X(k,t)と雑音重み係数W(k,t)の積の期待値として、以下の(4)式によって推定される。以下の(4)式において、演算子E,{}は時刻tに関する期待値を意味している。   In order for the suppression gain calculation means 202 to calculate the suppression gain, the noise power spectrum must first be estimated. When the input signal is a single channel, generally, the k-th element (synonymous with “k-th band”) N (k, t) of the noise power spectrum at time t is the element X (k, t) of the input power spectrum. ) And the noise weighting factor W (k, t) as an expected value is estimated by the following equation (4). In the following expression (4), operators E and {} mean expected values for time t.

また、以下の(4)式において、雑音重み係数W(k,t)は、入力パワースペクトルの要素が完全に雑音なら0とし、多分に音声を含むなら1とするものとする。なお、以下の(4)式において、音声区間を検出して雑音重み係数を0又は1の離散値としても良いし、音声区聞か雑音区間か暖昧な場合には雑音重み係数を0以上1以下の実数値としても良い。
N(k,t)=E{W(k,t)・X(k,t)} …(4)
In the following equation (4), the noise weighting coefficient W (k, t) is 0 if the input power spectrum element is completely noisy, and is 1 if it contains a lot of speech. In the following equation (4), the speech interval may be detected and the noise weighting coefficient may be a discrete value of 0 or 1. If the speech interval or the noise interval is warm, the noise weighting coefficient is 0 or more and 1 The following real values may be used.
N (k, t) = E t {W (k, t) · X (k, t)} (4)

そして、ここでは、複数チャネルから代表信号を生成し、該代表信号によって雑音パワースペクトルを推定する方法を定式化する。平均による代表信号の生成は重み付き加算の特別な場合であるから、ここでは、重み付き加算によって代表信号を生成する方法を適用するものとする。そして、チャネルmの入力パワースペクトルX(k,t)に対するチャネル重み係数をCとすると、雑音パワースペクトルの推定式は、以下の(5)式となる。ただし、Cは以下の(6)を満たすものとする。

Figure 0006613728
Here, a method for generating a representative signal from a plurality of channels and estimating a noise power spectrum using the representative signal is formulated. Since generation of a representative signal by averaging is a special case of weighted addition, a method of generating a representative signal by weighted addition is applied here. Then, assuming that the channel weight coefficient for the input power spectrum X m (k, t) of the channel m is C m , the noise power spectrum estimation formula is the following formula (5). However, C m shall meet the following (6).
Figure 0006613728

上述の(6)式が成り立つとすれば、チャネル重み係数Cは確率分布と見なせる。そこで、ここでは、以下の(7)式を満たすようなチャネル選択因子B(t)を導入する。チャネル選択因子B(t)は、ある時刻tにおいて唯一つのチャネル番号のみが1となり、他のチャネル番号はすべて0となるが、1となるチャネルmは確率分布Cに従って無作為に選ばれるような因子であるものとする。このようなチャネル選択因子B(t)を導入するために、まず以下の(5)式の期待値と、チャネルに関する和の計算順序を交代して以下の(8)式のように変形する。そして、以下の(8)式の確率分布Cをチャネル選択因子B(t)の期待値で書き直すと、以下の(9)式となる。チャネル選択因子B(t)は、入力パワースペクトルX(k,t)と雑音重み係数W(k,t)とは独立なので、以下の(10)式のように変形でき、最後に再び期待値とチャネルに関する和の計算順序を交代することで以下の(11)式を得ることができる。 If the above equation (6) is satisfied, the channel weighting coefficient C m is regarded as a probability distribution. Therefore, here, a channel selection factor B m (t) that satisfies the following expression (7) is introduced. The channel selection factor B m (t) is such that only one channel number is 1 at a certain time t and all other channel numbers are 0, but the channel m that is 1 is randomly selected according to the probability distribution C m. It is assumed that these are factors. In order to introduce such a channel selection factor B m (t), first, the expected value of the following equation (5) and the calculation order of the sum related to the channel are changed and transformed as the following equation (8). . Then, when the probability distribution C m in the following equation (8) is rewritten with the expected value of the channel selection factor B m (t), the following equation (9) is obtained. Since the channel selection factor B m (t) is independent of the input power spectrum X m (k, t) and the noise weighting factor W (k, t), the channel selection factor B m (t) can be transformed as the following equation (10), and finally again The following equation (11) can be obtained by changing the calculation order of the sum of the expected value and the channel.

そうすると、上述の(5)式と以下の(11)式の左辺が等しいので、チャネル選択因子B(t)を入力パワースペクトルX(k,t)に乗じてチャネルに関する和を取ることは、入力パワースペクトルX(k,t)から確率分布Cに従って無作為に1つのチャネルを選択することを意味することが分かる。

Figure 0006613728
Then, since the left side of the above equation (5) is equal to the following equation (11), the channel selection factor B m (t) is multiplied by the input power spectrum X m (k, t) to obtain the sum relating to the channel. It can be seen that this means that one channel is selected at random according to the probability distribution C m from the input power spectrum X m (k, t).
Figure 0006613728

従って、入力パワースペクトルX(k,t)から確率分布Cに従って無作為に1つのチャネルを選択して雑音パワースペクトルを推定することは、チャネル重み係数Cに基づいた重み付き加算によって代表信号を生成して雑音パワースペクトルを推定することと同じであることが分かる。 Therefore, estimating a noise power spectrum by randomly selecting one channel from the input power spectrum X m (k, t) according to the probability distribution C m is represented by weighted addition based on the channel weight coefficient C m. It can be seen that this is the same as generating a signal and estimating the noise power spectrum.

以上のように、雑音抑圧装置100において、上述の第1〜第3の選択方法のいずれかで選択入力スペクトルを切替えながら抑圧ゲインを生成する場合でも、雑音抑圧性能の維持が可能であることが、数学的にも裏付けられている。すなわち、雑音抑圧装置100では、抑圧ゲインの算出するための処理量オーダがO(n)未満に削減できることと、処理量が削減されても抑圧性能を維持できるといえる。   As described above, in the noise suppression apparatus 100, it is possible to maintain the noise suppression performance even when the suppression gain is generated while switching the selected input spectrum by any of the first to third selection methods described above. Mathematically supported. That is, in the noise suppression apparatus 100, it can be said that the processing amount order for calculating the suppression gain can be reduced to less than O (n) and that the suppression performance can be maintained even if the processing amount is reduced.

(B)第2の実施形態
以下、本発明による雑音抑圧装置、プログラム及び方法の第2の実施形態を、図面を参照しながら詳述する。
(B) Second Embodiment Hereinafter, a second embodiment of the noise suppression device, program and method according to the present invention will be described in detail with reference to the drawings.

(B−1)第2の実施形態の構成及び動作
第2の実施形態の雑音抑圧装置100Aの全体構成についても上述の図1を用いて説明することができる。以下では、第2の実施形態について第1の実施形態との差異を説明する。
(B-1) Configuration and Operation of Second Embodiment The overall configuration of the noise suppression device 100A of the second embodiment can also be described using FIG. Hereinafter, differences of the second embodiment from the first embodiment will be described.

第1の実施形態の雑音抑圧装置100では、上述の第1の選択方法(チャネル番号順に循環的にチャネル番号を選択する方法)で選択入力スペクトルに適用するチャネル番号を選択して抑圧ゲインを生成していた。これに対して、第2の実施形態の雑音抑圧装置100Aでは、上述の第2の選択方法(予め定められた順序でチャネル番号を選択する方法)で、択入力スペクトルISSに適用するチャネル番号を選択して抑圧ゲインを生成する。   In the noise suppression apparatus 100 according to the first embodiment, a suppression gain is generated by selecting a channel number to be applied to a selected input spectrum by the above-described first selection method (a method of cyclically selecting channel numbers in order of channel numbers). Was. On the other hand, in the noise suppression apparatus 100A of the second embodiment, the channel number applied to the selected input spectrum ISS is selected by the above-described second selection method (method of selecting channel numbers in a predetermined order). Select to generate suppression gain.

図1、図2に示すように、雑音抑圧装置100Aでは、帯域別抑圧手段102(チャネル選択手段201)が帯域別抑圧手段102A(チャネル選択手段201A)に置き換えられている点で第1の実施形態と異なっている。   As shown in FIGS. 1 and 2, the noise suppression apparatus 100A is the first implementation in that the band-by-band suppression unit 102 (channel selection unit 201) is replaced with the band-by-band suppression unit 102A (channel selection unit 201A). It is different from the form.

図6は、第2の実施形態に係るチャネル選択手段201Aの内部構成について示したブロック図である。   FIG. 6 is a block diagram showing an internal configuration of the channel selection unit 201A according to the second embodiment.

図6に示すように、チャネル選択手段201Aは、チャネル決定手段302がチャネル決定手段302Aに置き換えられ、さらにチャネル記憶手段303が除外されている点で第1の実施形態と異なっている。   As shown in FIG. 6, the channel selection unit 201A is different from the first embodiment in that the channel determination unit 302 is replaced with the channel determination unit 302A, and the channel storage unit 303 is excluded.

チャネル決定手段302Aは、フレーム計数手段301から与えられたフレーム番号に基づいて、チャネル番号を決定し、得られたチャネル番号を選択チャネル番号として入力スペクトル選択手段304に与える。   The channel determining unit 302A determines a channel number based on the frame number given from the frame counting unit 301, and gives the obtained channel number to the input spectrum selecting unit 304 as a selected channel number.

図7は、第2の実施形態のチャネル決定手段302Aの内部構成について示したブロック図である。   FIG. 7 is a block diagram showing the internal configuration of the channel determination unit 302A of the second embodiment.

図7に示すように、チャネル決定手段302Aは、チャネル順序供給手段601及びチャネル順序選択手段602を有している。   As shown in FIG. 7, the channel determination unit 302A includes a channel order supply unit 601 and a channel order selection unit 602.

チャネル決定手段302Aに与えられたフレーム番号はチャネル順序選択手段602に与えられる。   The frame number given to the channel determining means 302A is given to the channel order selecting means 602.

チャネル順序供給手段601は、予め設定されたチャネル番号の順序を記述したリスト(以下、「チャネル順序情報」と呼ぶ)を、チャネル順序選択手段602に与える。チャネル順序情報は、例えば、配列やベクトル、リスト等任意の形式ですることができるが、片方向循環リストによって定義するのが好適である。以下では、チャネル順序情報が片方向循環リストで定義されているものとする。また、チャネル順序情報の要素はチャネル番号1〜Mのいずれかに限定されるが、チャネル順序の要素数は任意である。ただし、チャネル順序情報に記述される要素数がM個未満の場合、一度も参照されないチャネルが生じてしまう。また、チャネル順序情報の要素数がちょうどM個の場合は、統計的には各チャネル番号の出現確率等は第1の実施形態と変わらない結果となる。したがって、チャネル順序情報に設定する要素数はM個より多く設定するのが好適である。例えばM=4の場合、チャネル順序情報を{2,2,2,3,3,1,4}と定義すると、チャネル番号2が最も重要視され、次いでチャネル番号3が重要視され、チャネル番号1、4は参照するが重要視されない、という重み付けとなる。また例えば、チャネル順序情報として、上述の{2,2,2,3,3,1,4}の順序を変更して{1,2,3,4,2,3,2}と設定してもよい。   The channel order supply unit 601 gives a list (hereinafter referred to as “channel order information”) describing a preset order of channel numbers to the channel order selection unit 602. The channel order information can be in an arbitrary format such as an array, a vector, or a list, for example, but is preferably defined by a one-way circulation list. In the following, it is assumed that the channel order information is defined by a one-way circulation list. The elements of the channel order information are limited to any one of channel numbers 1 to M, but the number of elements of the channel order is arbitrary. However, when the number of elements described in the channel order information is less than M, a channel that is never referenced is generated. In addition, when the number of elements of the channel order information is exactly M, the appearance probability of each channel number is statistically the same as in the first embodiment. Therefore, it is preferable to set more than M elements to be set in the channel order information. For example, when M = 4, if channel order information is defined as {2, 2, 2, 3, 3, 1, 4}, channel number 2 is most important, then channel number 3 is most important, and channel number 1 and 4 are weighted so that they are referenced but not regarded as important. Also, for example, as the channel order information, the order of {2, 2, 2, 3, 3, 1, 4} is changed and set to {1, 2, 3, 4, 2, 3, 2} Also good.

チャネル順序選択手段602は、フレーム番号に基づいてチャネル順序から1つのチャネルを選択して選択チャネル番号とする。チャネル順序選択手段602は、最初に選択するチャネル順序を初期チャネル順序として事前に定めておいて、フレーム番号が初期フレーム番号であれば該初期チャネル順序を選択チャネル番号とし、そうでなければ次のチャネル順序を選択チャネル番号とする。   The channel order selection means 602 selects one channel from the channel order based on the frame number and sets it as the selected channel number. The channel order selection means 602 determines in advance the channel order to be selected first as the initial channel order. If the frame number is the initial frame number, the initial channel order is set as the selected channel number. The channel order is the selected channel number.

フレーム計数手段301が供給する初期フレーム番号は実装によって異なるが、例えば0又は1とするようにしてもよい。   The initial frame number supplied by the frame counting means 301 varies depending on the implementation, but may be set to 0 or 1, for example.

(B−2)第2の実施形態の効果
第2の実施形態によれば以下のような効果を奏することができる。
(B-2) Effects of Second Embodiment According to the second embodiment, the following effects can be achieved.

第2の実施形態の雑音抑圧装置100Aでは、第1の実施形態と同様に、入力信号のチャネル数が増えた場合でも、チャネル数に応じてO(n)の処理量オーダとはならず、従来よりも少ない処理量で雑音抑制を行うことができる。   In the noise suppression device 100A of the second embodiment, as in the first embodiment, even when the number of channels of the input signal increases, the processing amount order of O (n) does not become an order according to the number of channels. Noise suppression can be performed with a smaller amount of processing than before.

また、第2の実施形態の雑音抑圧装置100Aでは、チャネル順序情報のリストの記述内容で、チャネル番号ごとに出現確率の重み付けを行うことができる。例えば、第2の実施形態の雑音抑圧装置100Aでは、複数チャネルの入力パワースペクトルを、事前に定めた重要視する度合にしたがって評価しながら雑音を抑庄できる。例えば、雑音の音源の位置に最も近いマイクで捕捉された入力信号に係るチャネル番号の出現頻度が多くなるようにチャネル順序情報のリストを記述することで、より精度の高い抑圧ゲインを取得することができる。   Further, in the noise suppression device 100A of the second embodiment, the appearance probability can be weighted for each channel number with the description content of the list of channel order information. For example, in the noise suppression device 100A of the second embodiment, it is possible to suppress noise while evaluating the input power spectra of a plurality of channels according to a predetermined importance degree. For example, a more accurate suppression gain can be obtained by describing a list of channel order information so that the appearance frequency of the channel number related to the input signal captured by the microphone closest to the position of the noise source is increased. Can do.

(C)第3の実施形態
以下、本発明による雑音抑圧装置、プログラム及び方法の第3の実施形態を、図面を参照しながら詳述する。
(C) Third Embodiment Hereinafter, a third embodiment of the noise suppression device, program and method according to the present invention will be described in detail with reference to the drawings.

(C−1)第3の実施形態の構成及び動作
第3の実施形態の雑音抑圧装置100Bの全体構成についても上述の図1を用いて説明することができる。以下では、第3の実施形態について第3の実施形態との差異を説明する。
(C-1) Configuration and Operation of Third Embodiment The overall configuration of the noise suppression device 100B of the third embodiment can also be described using FIG. 1 described above. Hereinafter, differences of the third embodiment from the third embodiment will be described.

第1の実施形態の雑音抑圧装置100では、上述の第1の選択方法(チャネル番号順に循環的にチャネル番号を選択する方法)で選択入力スペクトルに適用するチャネル番号を選択して抑圧ゲインを生成していた。これに対して、第3の実施形態の雑音抑圧装置100では、上述の第3の選択方法(疑似乱数に基づく順序でチャネル番号を選択する方法)で、選択入力スペクトルに適用するチャネル番号を選択して抑圧ゲインを生成する。   In the noise suppression apparatus 100 according to the first embodiment, a suppression gain is generated by selecting a channel number to be applied to a selected input spectrum by the above-described first selection method (a method of cyclically selecting channel numbers in order of channel numbers). Was. On the other hand, in the noise suppression apparatus 100 of the third embodiment, the channel number to be applied to the selected input spectrum is selected by the above-described third selection method (a method of selecting channel numbers in the order based on pseudo-random numbers). To generate a suppression gain.

図1、図2に示すように、雑音抑圧装置100Bでは、帯域別抑圧手段102(チャネル選択手段201)が帯域別抑圧手段102B(チャネル選択手段201B)に置き換えられている点で第1の実施形態と異なっている。   As shown in FIGS. 1 and 2, the noise suppression apparatus 100B is the first implementation in that the band-by-band suppression unit 102 (channel selection unit 201) is replaced with the band-by-band suppression unit 102B (channel selection unit 201B). It is different from the form.

図8は、第3の実施形態に係るチャネル選択手段201Bの内部構成について示したブロック図である。   FIG. 8 is a block diagram showing the internal configuration of the channel selection unit 201B according to the third embodiment.

図8に示すように、チャネル選択手段201Bは、チャネル決定手段302がチャネル決定手段302Bに置き換えられ、さらにフレーム計数手段301及びチャネル記憶手段303が省略されている点で第1の実施形態と異なっている。   As shown in FIG. 8, the channel selection unit 201B is different from the first embodiment in that the channel determination unit 302 is replaced with the channel determination unit 302B, and the frame counting unit 301 and the channel storage unit 303 are omitted. ing.

チャネル決定手段302Bは、疑似乱数に基づいてチャネル番号を決定し、得られたチャネル番号を選択チャネル番号として入力スペクトル選択手段304に与える。   The channel determination unit 302B determines a channel number based on the pseudo random number, and gives the obtained channel number to the input spectrum selection unit 304 as a selected channel number.

図9は、第3の実施形態のチャネル決定手段302Bの内部構成について示したブロック図である。   FIG. 9 is a block diagram showing the internal configuration of the channel determination unit 302B of the third embodiment.

図9に示すように、チャネル決定手段302Bは、乱数生成手段701、閾値セット供給手段702、及び閾値判定手段703を有している。   As illustrated in FIG. 9, the channel determination unit 302B includes a random number generation unit 701, a threshold set supply unit 702, and a threshold determination unit 703.

乱数生成手段701は、疑似乱数rを生成して閾値判定手段703に与える。乱数生成手段701が、疑似乱数の生成に用いる方法は限定されないものであるが、例えば、処理量の少ない線形合同法を適用することが望ましい。また、乱数生成手段701において、処理量やリソースが問題とならないのであれば、例えば、より周期が長く規則性の少ないメルセンヌ・ツイスタを適用することが望ましい。   The random number generation unit 701 generates a pseudo random number r and gives it to the threshold value determination unit 703. The method used by the random number generation unit 701 for generating pseudo-random numbers is not limited. For example, it is desirable to apply a linear congruential method with a small processing amount. Further, in the random number generation unit 701, for example, if the processing amount and the resource are not problematic, it is desirable to apply, for example, a Mersenne twister having a longer period and less regularity.

線形合同法の計算には整数値しか用いられないため、得られる素の疑似乱数は整数値(整数疑似乱数と呼ぶ)であるが、この実施形態の例では、乱数生成手段701は、整数疑似乱数を整数疑似乱数の最大値で除して0.0〜1.0の実数値に変換した疑似乱数rを出力するものとする。   Since only an integer value is used for the calculation of the linear congruential method, the prime pseudorandom number obtained is an integer value (referred to as an integer pseudorandom number). A pseudo-random number r obtained by dividing the random number by the maximum value of the integer pseudo-random number and converting it to a real value of 0.0 to 1.0 is output.

閾値セット供給手段702は、チャネル数Mより1少ない(M−1)個の、0.0より大きく1.0より小さい実数値の集合を事前に定義して、定義された集合を閾値セットR(i=1,…,M−1)として閾値判定手段703に与える。閾値セットの構成は、制限されるものではないが、配列を用いた平衡二分探索木とする(配列の要素がR2i<R<…<R2i+1を満たすようにする)のが好適である(図10)。図10の「配列」は、添え字iの順に閾値が配列に格納されている様を示す。同じく「データ構造」は、閾値セットが平衡二分探索木の構造を持っている様を示す。同じく「値」は、(M−1)個の閾値が範囲0.0〜1.0の数直線上のどこに位置するかのイメージを示す。図10では閾値が等間隔に位置しているが、閾値の値はこれに制限されるものではない。同じく「選択チャネル」は,数直線上の乱数の位置と、後述する閾値判定手段703によって選択される選択チャネル番号との関係を示す。 The threshold set supply unit 702 predefines a set of real values greater than 0.0 and less than 1.0 (M−1), which is 1 less than the number of channels M, and sets the defined set as the threshold set R. i (i = 1,..., M−1) is given to the threshold determination means 703. The configuration of the threshold set is not limited, but is preferably a balanced binary search tree using an array (so that the elements of the array satisfy R 2i <R i <... <R 2i + 1 ). (FIG. 10). “Array” in FIG. 10 indicates that threshold values are stored in the array in the order of the subscript i. Similarly, “data structure” indicates that the threshold set has a structure of an equilibrium binary search tree. Similarly, “value” indicates an image of where (M−1) threshold values are located on a number line in the range of 0.0 to 1.0. In FIG. 10, the threshold values are located at equal intervals, but the threshold value is not limited to this. Similarly, “selected channel” indicates the relationship between the position of the random number on the number line and the selected channel number selected by the threshold determination means 703 described later.

閾値判定手段703は、与えられた疑似乱数rと与えられた閾値セットRに基づいて定まるチャネル番号を、選択チャネル番号として設定する。選択チャネル番号mは、図11に示したフローに従って選択される。図11に従った閾値判定手段703の動作を説明する。はじめに、選択チャネル番号mにチャネル数Mを代入して初期化し(S1)、疑似乱数rに乱数生成手段701(図11では関数rand()に相当)が生成した値(0.0〜1.0の実数値)を代入し(S2)、閾値セットの添え字iに1を代入して初期化する(S3)。疑似乱数rと閾値Rとの比較を行い(S4)、r<Rであれば選択チャネル番号mに添え字iを代入し(S5)、iを2倍して更新し(S6)、r≧Rであればiを2倍して1を加えて更新する(S7)。そして、添え字iとチャネル数Mとの比較を行い(S8)、i<MであればS4〜S7を更新し、i≧Mであれば処理を終了して選択チャネル番号mを確定する。以上の方法によれば、疑似乱数と閾値を比較する回数は、高々logMを切り上げた整数の回数、すなわち処理オーダO(logn)となり、特許文献1のO(n)より少ない回数で実現できる。 The threshold determination unit 703 sets a channel number determined based on the given pseudorandom number r and the given threshold set R i as the selected channel number. The selected channel number m is selected according to the flow shown in FIG. The operation of the threshold determination unit 703 according to FIG. 11 will be described. First, initialization is performed by substituting the channel number M for the selected channel number m (S1), and a value (0.0 to 1....) Generated by the random number generation means 701 (corresponding to the function rand () in FIG. 11) for the pseudorandom number r. (Real value of 0) is substituted (S2), and 1 is substituted for the subscript i of the threshold set for initialization (S3). Compares the pseudo-random number r and the threshold value R i (S4), by substituting r <R i is long if the selected channel number m subscript i (S5), and update doubled i (S6), If r ≧ R i , i is doubled and 1 is added to update (S7). Then, the subscript i is compared with the number of channels M (S8). If i <M, S4 to S7 are updated, and if i ≧ M, the process is terminated and the selected channel number m is determined. According to the above method, the number of comparisons between the pseudo random number and the threshold value is an integer number obtained by rounding up log 2 M at most, that is, the processing order O (log 2 n), which is smaller than O (n) in Patent Document 1. Can be realized.

したがって、第3の実施形態のチャネル決定手段302Bでは、閾値配列Rの定め方(各チャネル番号に適用するレンジ)によって、どのチャネルをどの程度重要視するのか(各チャネル番号の出現率)を調整することができる。例えば、R=i/Mとすれば、各チャネル番号の出現率は均等となり、各チャネル番号の出現率等は統計的には第1の実施形態と同じ結果となる。また、重要視するチャネル番号に対応するレンジを広く調整することで、第2の実施形態と同様に、統計的には重要視するチャネル番号の出現率を上げること等の調整が可能となる。 Therefore, in the channel determination unit 302B of the third embodiment, which channel is to be regarded as important (appearance rate of each channel number) depending on how the threshold value array R i is determined (range applied to each channel number). Can be adjusted. For example, if R i = i / M, the appearance rates of the respective channel numbers are equal, and the appearance rates of the respective channel numbers are statistically the same as those in the first embodiment. Further, by adjusting the range corresponding to the channel number to be regarded as important, it is possible to make adjustments such as increasing the appearance rate of the channel numbers regarded as important statistically, as in the second embodiment.

(C−2)第3の実施形態の効果
第3の実施形態によれば以下のような効果を奏することができる。
(C-2) Effects of Third Embodiment According to the third embodiment, the following effects can be achieved.

第3の実施形態の雑音抑圧装置100Bでは、第1、第2の実施形態と同様に、入力信号のチャネル数が増えた場合でも、チャネル数に応じてO(n)の処理量オーダとはならず、従来よりも少ない処理量で雑音抑制を行うことができる。第3の実施形態の雑音抑圧装置100Bでは、単位時間ごとに、乱数の生成は一回、閾値判定は高々logMを切り上げた整数の回数しか実行されないため、処理量オーダは平衡二分探索木と同じO(logn)となる。 In the noise suppression device 100B of the third embodiment, as in the first and second embodiments, even when the number of channels of the input signal increases, the processing amount order of O (n) according to the number of channels In addition, noise suppression can be performed with a smaller processing amount than in the past. In the noise suppression device 100B according to the third embodiment, random number generation is executed once per unit time, and threshold determination is executed only by an integer number obtained by rounding up log 2 M at most, so that the processing amount order is a balanced binary search tree. The same O (log 2 n).

また、第3の実施形態では、チャネル決定手段302Bに設定する閾値配列R(i=1,…,M−1)の内容を調整(各チャネル番号に対応するレンジを調整)することで、複数チャネルの入力パワースペクトルを事前に定めた自由度の高い重要視する度合にしたがって評価しながら雑音を抑制することができる。 Further, in the third embodiment, by adjusting the contents of the threshold array R i (i = 1,..., M−1) set in the channel determining unit 302B (adjusting the range corresponding to each channel number), Noise can be suppressed while evaluating the input power spectrum of a plurality of channels according to a predetermined degree of importance with a high degree of freedom.

第1の実施形態および第2の実施形態では、チャネルを選択する順序を一意的に定めていた。第1の実施形態では、複数チャネルを平均的に参照することしかできない。第2の実施形態では、各チャネルを重要視する度合を自由に設定するのが容易ではない。例えば、第2の実施形態において、チャネル数M=2で、チャネル1を30%、チャネル2を70%重要視したい場合には、わずか2チャネルのために{1,1,1,2,2,2,2,2,2,2}等のような要素数が10個以上のチャネル順序を用意しなければならない。しかしながら、第3の実施形態の雑音抑圧装置100Bでは疑似乱数を発生させ、発生させた疑似乱数と事前に定めた閾値に基づいてチャネルを決定する。これにより、例えばチャネル数M=2で重要視する度合を13%と87%にするというような、第2の実施形態では難しい設定も容易に実現できる。すなわち、第3の実施形態では、第2の実施形態の場合よりも各チャネル番号の出現率を容易に調整することができる。   In the first embodiment and the second embodiment, the order of selecting channels is uniquely determined. In the first embodiment, only a plurality of channels can be referred to on average. In the second embodiment, it is not easy to freely set the degree of importance of each channel. For example, in the second embodiment, if the number of channels is M = 2, and channel 1 is 30% and channel 2 is 70% important, {1, 1, 1, 2, 2 for only 2 channels , 2, 2, 2, 2, 2} and the like, a channel order having 10 or more elements must be prepared. However, the noise suppression apparatus 100B of the third embodiment generates a pseudo random number, and determines a channel based on the generated pseudo random number and a predetermined threshold value. This makes it possible to easily realize settings that are difficult in the second embodiment, such as 13% and 87%, for example, when the number of channels M = 2. That is, in the third embodiment, the appearance rate of each channel number can be adjusted more easily than in the second embodiment.

(D)他の実施形態
本発明は、上記の各実施形態に限定されるものではなく、以下に例示するような変形実施形態も挙げることができる。
(D) Other Embodiments The present invention is not limited to the above-described embodiments, and may include modified embodiments as exemplified below.

(D−1)第2及び第3の実施形態において、ユーザが、チャネル順序情報及び閾値配列Rを変更(編集)可能とするようにしてもよい。 In (D-1) second and third embodiments, a user may be capable of changing the channel order information and threshold array R i (Edit).

また、第2の実施形態において、重要視する度合をユーザが途中で変更可能とする場合、重要視する度合にしたがってチャネル順序を自動生成するチャネル順序自動生成手段を設けるようにしてもよい。このチャネル順序自動生成手段は、例えば、チャネル順序の最大要素数を事前に定めておいて、ユーザが指定した重要視する度合ができるだけ正しく反映されるようにチャネル順序を構成することが望ましい。   In the second embodiment, when the user can change the degree of importance on the way, channel order automatic generation means for automatically generating the channel order according to the degree of importance may be provided. It is desirable that this automatic channel order generation means, for example, determines the maximum number of elements of the channel order in advance and configures the channel order so that the degree of importance designated by the user is reflected as accurately as possible.

さらに、第3の実施形態において、重要視する度合をユーザが途中で変更可能とする場合、閾値配列供給手段の代わりにユーザの操作を受けて閾値配列を自動で調整する閾値配列生成手段を構成に含めることが望ましい。   Furthermore, in the third embodiment, when the user can change the degree of importance on the way, a threshold value array generating unit configured to automatically adjust the threshold value array in response to a user operation instead of the threshold value array supplying unit is configured. It is desirable to include.

(D−2)第2の実施形態及び第3の実施形態において、複数チャネルの入力パワースペクトルに応じて重要視する度合を自動で変更するチャネル重要度適応手段を構成に含めても良い。当該チャネル重要度適応手段は、例えば、各チャネルの音声パワー指標値が最大となるチャネルをより重要視するように、重要視する度合を変更する構成とするようにしてもよい。音声パワー指標値としては、例えば音声パワーの一定期間の最大値を適用することが望ましい。音声パワー指標値は、ある時刻においては選択チャネル番号のチャネルでのみ更新されるようにしてもよい。当該チャネル重要度適応手段は、例えば、1単位時間前における音声パワー指標値に忘却係数β(例えばβ=0.99)を乗じた値と、選択チャネル番号の入力パワースペクトルの最大要素とのうち大きい方を、新しい音声パワー指標値として更新するようにしてもよい。当該チャネル重要度適応手段は、例えば、1単位時間前における複数チャネルの音声パワー指標値のうち、最大の音声パワー指標値を記憶しておいて、選択チャネル番号の音声パワー指標値が当該最大の音声パワー指標値より大きければ該選択チャネル番号がより重要視されるように、重要視する度合を変更するようにしてもよい。このようにすることで、処理量オーダをO(logn)に保ちながら重要視する度合を変更することが可能となる。 (D-2) In the second and third embodiments, channel importance adapting means for automatically changing the degree of importance according to the input power spectrum of a plurality of channels may be included in the configuration. For example, the channel importance level adaptation unit may be configured to change the degree of importance so that the channel having the maximum audio power index value of each channel is more important. As the audio power index value, for example, it is desirable to apply the maximum value of audio power for a certain period. The audio power index value may be updated only on the channel of the selected channel number at a certain time. For example, the channel importance adapting means includes a value obtained by multiplying a voice power index value one unit time before by a forgetting factor β (for example, β = 0.99) and a maximum element of the input power spectrum of the selected channel number. The larger one may be updated as a new voice power index value. The channel importance adapting means stores, for example, the maximum audio power index value among the audio power index values of a plurality of channels one unit time before, and the audio power index value of the selected channel number is the maximum The degree of importance may be changed so that the selected channel number is more important if the value is larger than the voice power index value. In this way, it is possible to change the degree of emphasis while maintaining the processing amount order at O (log 2 n).

101、101−1〜101−M…周波数解析手段、102…帯域別抑圧手段、103−1〜103−M…波形復元手段、201…チャネル選択手段、202…抑圧ゲイン算出手段、203−1〜203−M…乗算手段、301…フレーム計数手段、302…チャネル決定手段、303…チャネル記憶手段、304…入力スペクトル選択手段、401…チャネル計数手段、402…定数1値供給手段、403…初期チャネルスイッチ手段、501…パワー算出手段、502…雑音平滑化手段、503…雑音記憶手段、504…事後SNR算出手段、505…事後SNR記憶手段、506…閾値判定手段、507…閾値供給手段、508…事前SNR推定手段、509…抑圧ゲイン決定手段、510…抑圧ゲイン記憶手段。   101, 101-1 to 101-M ... frequency analysis means, 102 ... band-based suppression means, 103-1 to 103-M ... waveform restoration means, 201 ... channel selection means, 202 ... suppression gain calculation means, 203-1 203-M... Multiplying means 301... Frame counting means 302. Channel determining means 303... Channel storage means 304... Input spectrum selecting means 401. Switch means 501... Power calculation means 502... Noise smoothing means 503. Noise storage means 504 .. Post-SNR calculation means 505 .. Post-SNR storage means 506... Threshold determination means 507. Prior SNR estimation means, 509 ... suppression gain determination means, 510 ... suppression gain storage means.

Claims (7)

複数チャネルの入力信号のそれぞれに対応する入力スペクトルから雑音を抑圧する雑音抑圧装置において、
所定の単位時間ごとに、所定の方式により前記複数チャネルから1つのチャネルを選択するチャネル選択手段と、
選択されたチャネルの入力スペクトルを用いて当該入力スペクトルに含まれる雑音成分を抑圧するための抑圧ゲインを算出する抑圧ゲイン算出手段と、
前記抑圧ゲインを用いて、前記複数チャネルの入力スペクトルのそれぞれについて雑音成分を抑圧する雑音抑圧手段と
を有すことを特徴とする雑音抑圧装置。
In a noise suppression device that suppresses noise from an input spectrum corresponding to each of multiple channel input signals,
Channel selection means for selecting one channel from the plurality of channels by a predetermined method every predetermined unit time;
Suppression gain calculating means for calculating a suppression gain for suppressing a noise component included in the input spectrum using the input spectrum of the selected channel;
A noise suppression device comprising: noise suppression means for suppressing a noise component for each of the input spectra of the plurality of channels using the suppression gain.
前記複数チャネルを構成するチャネルには、それぞれチャネル番号が付与されており、
前記チャネル選択手段は、前記複数チャネルからチャネル番号順に循環的にチャネル選択を行うことを特徴とする請求項1に記載の雑音抑圧装置。
Channel numbers are assigned to the channels constituting the plurality of channels,
The noise suppression apparatus according to claim 1, wherein the channel selection means performs channel selection cyclically in the order of channel numbers from the plurality of channels.
前記チャネル選択手段は、前記複数チャネルからチャネルを選択する順序が記述されたチャネル順序情報に従った順序で、前記複数チャネルからいずれかのチャネルを選択することを特徴とする請求項1に記載の雑音抑圧装置。   2. The channel selection unit according to claim 1, wherein the channel selection unit selects one of the plurality of channels in an order according to channel order information in which an order of selecting a channel from the plurality of channels is described. Noise suppression device. 疑似乱数を生成する疑似乱数生成手段をさらに備え、
前記チャネル選択手段は、前記疑似乱数生成手段で生成された疑似乱数に基づいて、前記複数チャネルから、いずれかのチャネルを選択することを特徴とする請求項1に記載の雑音抑圧装置。
A pseudo random number generating means for generating a pseudo random number;
The noise suppression apparatus according to claim 1, wherein the channel selection unit selects one of the plurality of channels based on the pseudorandom number generated by the pseudorandom number generation unit.
前記複数チャネルのチャネル数より1少ない数の閾値で構成される閾値配列を保持する閾値配列手段をさらに備え、
前記チャネル選択手段は、前記疑似乱数生成手段で生成された疑似乱数と、前記閾値配列を構成する閾値との比較結果に基づいたチャネルを選択する
ことを特徴とする請求項4に記載の雑音抑圧装置。
Further comprising threshold value array means for holding a threshold value array composed of threshold values that are one less than the number of channels of the plurality of channels;
5. The noise suppression according to claim 4, wherein the channel selection unit selects a channel based on a comparison result between the pseudo-random number generated by the pseudo-random number generation unit and the threshold value constituting the threshold value array. apparatus.
複数チャネルの入力信号のそれぞれに対応する入力スペクトルから雑音を抑圧する雑音抑圧装置に搭載されたコンピュータを、
所定の単位時間ごとに、所定の方式により前記複数チャネルから1つのチャネルを選択するチャネル選択手段と、
選択されたチャネルの入力スペクトルを用いて当該入力スペクトルに含まれる雑音成分を抑圧するための抑圧ゲインを算出する抑圧ゲイン算出手段と、
前記抑圧ゲインを用いて、前記複数チャネルの入力スペクトルのそれぞれについて雑音成分を抑圧する雑音抑圧手段と
して機能させることを特徴とする雑音抑圧プログラム。
A computer mounted in a noise suppression device that suppresses noise from the input spectrum corresponding to each of the input signals of multiple channels,
Channel selection means for selecting one channel from the plurality of channels by a predetermined method every predetermined unit time;
Suppression gain calculating means for calculating a suppression gain for suppressing a noise component included in the input spectrum using the input spectrum of the selected channel;
A noise suppression program that functions as a noise suppression unit that suppresses a noise component for each of the input spectra of the plurality of channels using the suppression gain.
複数チャネルの入力信号のそれぞれに対応する入力スペクトルから雑音を抑圧する雑音抑圧方法において、
チャネル選択手段、抑圧ゲイン算出手段、及び雑音抑圧手段を備え、
前記チャネル選択手段は、所定の単位時間ごとに、所定の方式により前記複数チャネルから1つのチャネルを選択し、
前記抑圧ゲイン算出手段は、選択されたチャネルの入力スペクトルを用いて当該入力スペクトルに含まれる雑音成分を抑圧するための抑圧ゲインを算出し、
前記雑音抑圧手段は、前記抑圧ゲインを用いて、前記複数チャネルの入力スペクトルのそれぞれについて雑音成分を抑圧する
ことを特徴とする雑音抑圧方法。
In a noise suppression method for suppressing noise from an input spectrum corresponding to each of input signals of a plurality of channels,
A channel selection means, a suppression gain calculation means, and a noise suppression means,
The channel selection means selects one channel from the plurality of channels by a predetermined method every predetermined unit time,
The suppression gain calculation means calculates a suppression gain for suppressing a noise component included in the input spectrum using the input spectrum of the selected channel,
The noise suppression method, wherein the noise suppression means suppresses a noise component for each of the input spectra of the plurality of channels using the suppression gain.
JP2015169943A 2015-08-31 2015-08-31 Noise suppression device, program and method Active JP6613728B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015169943A JP6613728B2 (en) 2015-08-31 2015-08-31 Noise suppression device, program and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015169943A JP6613728B2 (en) 2015-08-31 2015-08-31 Noise suppression device, program and method

Publications (2)

Publication Number Publication Date
JP2017049277A JP2017049277A (en) 2017-03-09
JP6613728B2 true JP6613728B2 (en) 2019-12-04

Family

ID=58280303

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015169943A Active JP6613728B2 (en) 2015-08-31 2015-08-31 Noise suppression device, program and method

Country Status (1)

Country Link
JP (1) JP6613728B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109562611B (en) 2016-08-02 2021-05-25 日本聚丙烯株式会社 Decorative film and method for producing decorative molded body using the same

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007116585A (en) * 2005-10-24 2007-05-10 Matsushita Electric Ind Co Ltd Noise cancel device and noise cancel method
JP5435204B2 (en) * 2006-07-03 2014-03-05 日本電気株式会社 Noise suppression method, apparatus, and program
CN103814584B (en) * 2011-10-14 2017-02-15 富士通株式会社 Sound processing device and sound processing method

Also Published As

Publication number Publication date
JP2017049277A (en) 2017-03-09

Similar Documents

Publication Publication Date Title
Erdogan et al. Improved mvdr beamforming using single-channel mask prediction networks.
JP4469882B2 (en) Acoustic signal processing method and apparatus
KR101910679B1 (en) Noise adaptive beamforming for microphone arrays
CN109661705B (en) Sound source separation device and method, and program
CN111133511B (en) sound source separation system
RU2420813C2 (en) Speech quality enhancement with multiple sensors using speech status model
JP6400218B2 (en) Audio source isolation
JP5967571B2 (en) Acoustic signal processing apparatus, acoustic signal processing method, and acoustic signal processing program
JP2017503388A5 (en)
JP2011107603A (en) Speech recognition device, speech recognition method and program
JP6987075B2 (en) Audio source separation
JP6724905B2 (en) Signal processing device, signal processing method, and program
US10718742B2 (en) Hypothesis-based estimation of source signals from mixtures
CN105580074B (en) Signal processing system and method
US20160042746A1 (en) Noise suppressing device, noise suppressing method, and a non-transitory computer-readable recording medium storing noise suppressing program
JP6538624B2 (en) Signal processing apparatus, signal processing method and signal processing program
JP5034735B2 (en) Sound processing apparatus and program
KR101243897B1 (en) Blind Source separation method in reverberant environments based on estimation of time delay and attenuation of the signals
JP6613728B2 (en) Noise suppression device, program and method
Barber et al. End-to-end Alexa device arbitration
JP5263020B2 (en) Signal processing device
JP6707914B2 (en) Gain processing device and program, and acoustic signal processing device and program
JP6596833B2 (en) Noise suppression device and program, noise estimation device and program, and SNR estimation device and program
Zheng et al. Statistical analysis and improvement of coherent-to-diffuse power ratio estimators for dereverberation
US20210174820A1 (en) Signal processing apparatus, voice speech communication terminal, signal processing method, and signal processing program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180515

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190419

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190507

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190705

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20191021

R150 Certificate of patent or registration of utility model

Ref document number: 6613728

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150