JP2012215715A - Noise estimation device and noise estimation program - Google Patents

Noise estimation device and noise estimation program Download PDF

Info

Publication number
JP2012215715A
JP2012215715A JP2011081276A JP2011081276A JP2012215715A JP 2012215715 A JP2012215715 A JP 2012215715A JP 2011081276 A JP2011081276 A JP 2011081276A JP 2011081276 A JP2011081276 A JP 2011081276A JP 2012215715 A JP2012215715 A JP 2012215715A
Authority
JP
Japan
Prior art keywords
noise
frequency characteristic
information
estimated
estimated noise
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.)
Granted
Application number
JP2011081276A
Other languages
Japanese (ja)
Other versions
JP5768450B2 (en
Inventor
Kaori Endo
香緒里 遠藤
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.)
Fujitsu Ltd
Original Assignee
Fujitsu 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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2011081276A priority Critical patent/JP5768450B2/en
Publication of JP2012215715A publication Critical patent/JP2012215715A/en
Application granted granted Critical
Publication of JP5768450B2 publication Critical patent/JP5768450B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Fittings On The Vehicle Exterior For Carrying Loads, And Devices For Holding Or Mounting Articles (AREA)

Abstract

PROBLEM TO BE SOLVED: To estimate noise with high estimation accuracy even for the noise which unexpectedly varies.SOLUTION: First noise estimation means (4) outputs a first estimated noise frequency characteristic indicating a noise frequency characteristic of noise included in input sound (6). Second noise estimation means (8) outputs a second estimated noise frequency characteristic of the noise specified by using prior information (12) on the noise relevant to the input sound as well as current information (10) to specify the noise. Adding means (14) outputs an estimated noise frequency characteristic in accordance with the first estimated noise frequency characteristic as well as the second estimated noise frequency characteristic.

Description

本発明は、騒音の周波数特性を推定する騒音推定装置及び騒音推定プログラムに関する。
The present invention relates to a noise estimation device and a noise estimation program for estimating frequency characteristics of noise.

騒音の周波数特性を推定する技術について、音声区間と非音声区間を判別し、非音声区間で騒音のスペクトルを学習し、騒音の周波数特性を推定することが知られている。   As a technique for estimating the frequency characteristics of noise, it is known to discriminate between speech and non-speech sections, learn a noise spectrum in the non-speech section, and estimate the noise frequency characteristics.

この周波数特性の推定技術について、特許文献1には音声期間では背景雑音推定で推定雑音信号の更新を停止させるか、または、雑音期間より遅い入力信号に対する追従速度で更新させることが記載されている。   Regarding this frequency characteristic estimation technique, Patent Document 1 describes that in an audio period, updating of an estimated noise signal is stopped by background noise estimation, or updating is performed at a tracking speed for an input signal that is slower than the noise period. .

また、特許文献2には、環境雑音が重畳した音声信号から、音声信号の電力、零交差率、パワースペクトルのピーク周波数、ピッチ周期などの音声特徴量を用いることにより、音声、非音声の判定を行うことが記載されている。   Further, Patent Document 2 discloses determination of speech and non-speech by using speech feature values such as power of the speech signal, zero-crossing rate, peak frequency of power spectrum, and pitch period from the speech signal on which environmental noise is superimposed. It is described to do.

特開平09−311698号公報JP 09-31698 A 特開2002−258881号公報JP 2002-258881 A

ところで、騒音の定常的な性質を捉えて周波数特性を推定した場合、突発的に騒音が変化した場合には実際の騒音と学習した騒音とが相違し、両者間の誤差が大きくなる。特に、音声が含まれる区間で騒音が変化した場合や、突発的な騒音の場合では、音の変化を学習することが難しく、誤差が顕著になる。   By the way, when the frequency characteristics are estimated by capturing the steady nature of the noise, if the noise suddenly changes, the actual noise and the learned noise are different, and the error between them increases. In particular, when the noise changes in a section including speech or sudden noise, it is difficult to learn the change of the sound, and the error becomes remarkable.

そこで、本開示の騒音推定装置および騒音推定方法の目的は、突発的に変化する騒音であっても、推定精度の高い騒音推定を行うことにある。
Therefore, an object of the noise estimation device and the noise estimation method of the present disclosure is to perform noise estimation with high estimation accuracy even for noise that changes suddenly.

上記目的を達成するため、本開示の騒音推定装置および騒音推定プログラムは、第1の騒音推定手段が入力音に含まれる騒音の騒音周波数特性を示す第1の推定騒音周波数特性を出力する。第2の騒音推定手段が前記入力音に関する騒音の事前情報と、前記騒音を特定する現在情報とを用いて特定された騒音の第2の推定騒音周波数特性を出力する。そして、加算手段は前記第1の推定騒音周波数特性と前記第2の推定騒音周波数特性とに応じた推定騒音周波数特性を出力する。
In order to achieve the above object, in the noise estimation device and the noise estimation program of the present disclosure, the first noise estimation unit outputs a first estimated noise frequency characteristic indicating the noise frequency characteristic of the noise included in the input sound. The second noise estimation means outputs a second estimated noise frequency characteristic of the noise specified using the prior information of the noise related to the input sound and the current information specifying the noise. Then, the adding means outputs an estimated noise frequency characteristic corresponding to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic.

本開示の騒音推定装置および騒音推定プログラムによれば、次のいずれかの効果が得られる。   According to the noise estimation device and the noise estimation program of the present disclosure, any of the following effects can be obtained.

(1) 入力音から得られた第1の騒音周波数特性と、入力音に関する騒音の事前情報とその騒音を特定する現在情報とから求めた第2の騒音周波数特性とに応じた推定騒音周波数特性を得ており、推定精度を高めることができる。   (1) Estimated noise frequency characteristic according to the first noise frequency characteristic obtained from the input sound and the second noise frequency characteristic obtained from the prior information of the noise related to the input sound and the current information specifying the noise Thus, the estimation accuracy can be improved.

(2) 周波数特性の推定に現在情報や事前情報を反映させているので、入力音に音声が含まれている場合や、トンネルなど突発的な騒音が変化する場合でも、騒音の推定精度を低下させることなく、騒音推定を行うことができる。
(2) Since current information and prior information are reflected in the estimation of frequency characteristics, the accuracy of noise estimation is reduced even if the input sound contains speech or sudden noise changes such as in tunnels. It is possible to perform noise estimation without making it.

そして、本発明の他の目的、特徴及び利点は、添付図面及び各実施の形態を参照することにより、一層明確になるであろう。
Other objects, features, and advantages of the present invention will become clearer with reference to the accompanying drawings and each embodiment.

第1の実施の形態に係る騒音推定装置および騒音推定プログラムの一例を示す図である。It is a figure which shows an example of the noise estimation apparatus and noise estimation program which concern on 1st Embodiment. 第2の実施の形態に係る騒音推定装置の構成例を示す図である。It is a figure which shows the structural example of the noise estimation apparatus which concerns on 2nd Embodiment. 特定騒音情報テーブルの構成例を示す図である。It is a figure which shows the structural example of a specific noise information table. 他の特定騒音情報テーブルの一例を示す図である。It is a figure which shows an example of another specific noise information table. 騒音推定の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of noise estimation. 第3の実施の形態に係る騒音推定装置の構成例を示す図である。It is a figure which shows the structural example of the noise estimation apparatus which concerns on 3rd Embodiment. 特定騒音情報テーブルの一例を示す図である。It is a figure which shows an example of a specific noise information table. 騒音推定の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of noise estimation. 第4の実施の形態に係る騒音推定装置の構成例を示す図である。It is a figure which shows the structural example of the noise estimation apparatus which concerns on 4th Embodiment. 特定音源騒音の周波数特性の更新処理手順の一例を示すフローチャートである。It is a flowchart which shows an example of the update process procedure of the frequency characteristic of specific sound source noise. 騒音推定の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of noise estimation. 第5の実施の形態に係る騒音推定装置の構成例を示す図である。It is a figure which shows the structural example of the noise estimation apparatus which concerns on 5th Embodiment. 特定音源騒音の周波数特性の更新処理手順の一例を示すフローチャートである。It is a flowchart which shows an example of the update process procedure of the frequency characteristic of specific sound source noise. 騒音推定の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of noise estimation. 第6の実施の形態に係る騒音推定装置の構成例を示す図である。It is a figure which shows the structural example of the noise estimation apparatus which concerns on 6th Embodiment. 騒音推定の処理手順を示すフローチャートである。It is a flowchart which shows the process sequence of noise estimation. コンピュータシステムの構成例を示す図である。It is a figure which shows the structural example of a computer system.

〔第1の実施の形態〕 [First Embodiment]

図1は第1の実施の形態に係る騒音推定装置および騒音推定プログラムの一例を示している。図1に示す構成は一例であり、斯かる構成に本発明が限定されるものではない。   FIG. 1 shows an example of a noise estimation device and a noise estimation program according to the first embodiment. The configuration shown in FIG. 1 is an example, and the present invention is not limited to such a configuration.

図1に示す騒音推定装置2には第1の騒音推定手段4が備えられ、この騒音推定手段4は入力音6に含まれる騒音の騒音周波数特性を示す第1の推定騒音周波数特性を出力する。一例としての入力音6には利用者の音声が含まれることを想定している。   The noise estimation apparatus 2 shown in FIG. 1 includes first noise estimation means 4, which outputs first estimated noise frequency characteristics indicating noise frequency characteristics of noise included in the input sound 6. . As an example, it is assumed that the input sound 6 includes a user's voice.

第2の騒音推定手段8は、現在情報10と事前情報12とを用いることにより、入力音6に対応する騒音周波数特性を算出し、第2の推定騒音周波数特性を出力する。現在情報10は日時、位置、エンジン回転数などの現在の特定騒音の情報である。また、事前情報12は、新幹線などの交通機関の経路やタイムテーブルなど、予め取得できる情報である。   The second noise estimation means 8 calculates the noise frequency characteristic corresponding to the input sound 6 by using the current information 10 and the prior information 12, and outputs the second estimated noise frequency characteristic. The current information 10 is current specific noise information such as date and time, position, and engine speed. The prior information 12 is information that can be acquired in advance such as a route of a transportation facility such as a bullet train or a time table.

そして、加算手段14はたとえば、第1の推定騒音周波数特性と第2の推定騒音周波数特性との加算により、第1の推定騒音周波数特性と第2の推定騒音周波数特性とに応じた推定騒音周波数特性を出力する。この加算手段14は、これらの推定騒音周波数特性を同一の比率で加算してもよいし、予め定めた寄与係数で重み付けした値で加算してもよい。   Then, the adding means 14 adds, for example, the first estimated noise frequency characteristic and the second estimated noise frequency characteristic, thereby estimating the estimated noise frequency according to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic. Output characteristics. The adding means 14 may add these estimated noise frequency characteristics at the same ratio, or may add the weighted values with a predetermined contribution coefficient.

斯かる構成では、第1の推定騒音周波数特性と第2の推定騒音周波数特性とに応じた騒音周波数特性が出力され、騒音が推定されるので、高い推定精度を得ることができる。入力音6として、周囲騒音が重畳した音声信号(重畳騒音)の周波数特性の推定に用いることができる。また、列車などの経路や運行時刻が定まっている環境で、予めエンジン音やトンネルの走行音などの騒音の周波数特性を事前情報として保持するので、騒音の推定周波数特性の推定性能を高めることができる。   In such a configuration, since the noise frequency characteristic corresponding to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic is output and noise is estimated, high estimation accuracy can be obtained. The input sound 6 can be used to estimate the frequency characteristics of a voice signal (superimposed noise) on which ambient noise is superimposed. In addition, since the frequency characteristics of noise such as engine noise and tunnel running noise are stored in advance as an advance information in an environment where the route and operating time of trains are fixed, it is possible to improve the estimation performance of the estimated frequency characteristics of noise. it can.

〔第2の実施の形態〕 [Second Embodiment]

図2は第2の実施の形態に係る騒音推定装置の一例を示している。図2に示す構成は一例であり、斯かる構成に本発明が限定されるものではない。図2において、図1と同一部分には同一符号を付してある。   FIG. 2 shows an example of a noise estimation apparatus according to the second embodiment. The configuration shown in FIG. 2 is an example, and the present invention is not limited to such a configuration. 2, the same parts as those in FIG. 1 are denoted by the same reference numerals.

図2に示す騒音推定装置2−1は、第1の騒音推定手段4−1と、事前情報保持手段16と、現存特定騒音情報算出手段18と、第2の騒音推定手段8と、第3の騒音推定手段20を備えている。   The noise estimation apparatus 2-1 shown in FIG. 2 includes a first noise estimation unit 4-1, a prior information holding unit 16, an existing specific noise information calculation unit 18, a second noise estimation unit 8, and a third Noise estimation means 20 is provided.

騒音推定手段4−1は、入力音から騒音の周波数特性を算出し、これを第1の推定騒音周波数特性として出力する。一例としての入力音には音声が含まれる入力音が含まれる。   The noise estimation means 4-1 calculates the frequency characteristic of noise from the input sound and outputs this as the first estimated noise frequency characteristic. As an example, the input sound includes an input sound including sound.

この騒音推定手段4−1の処理には、騒音推定の方法や音声判定の方法を用いることができる。騒音推定の方法は、時間軸上の入力信号を周波数軸上の信号に変換し、変換された周波数軸上の信号に基づいて、その入力信号における背景騒音を入力信号に追従しながら推定し、推定雑音信号を形成する方法、(たとえば、特許第03269969号)を利用すればよい。即ち、音声を含まない雑音(背景雑音)の信号スペクトル信号Xf,kと1分析フレーム前の推定雑音信号N' f,k−1を用いて、現分析フレームの推定雑音信号N' f,kを算出する方法を利用すればよい。この場合、Xf,kとN' f,k−1に係数を掛けてそれぞれをa・N' f,k−1及び(1−a)・Xf,kとし、a・N' f,k−1と(1−a)・Xf,kを加算して、現分析フレームの推定雑音信号N' f,kが算出される。また、音声判定の方法は、入力信号として環境騒音が重畳した音声信号が時系列順に与えられるフレーム毎に、当該フレームに音声が含まれているか否かを判定する方法(たとえば、特許第03849116号)を利用すればよい。即ち、音声信号の特徴量と、音声周波数の特徴量の高域成分を強調させた特徴量を算出して、音声信号の特徴量に加えて高域成分を強調させた特徴量を基に判定を行う方法を利用すれば良い。また、音声信号の特徴量と、音声の共振周期の規則性を算出して、音声信号の特徴量に加えて音声の共振周期の規則性を基に判定を行う方法を利用すれば良い。   For the processing of the noise estimation means 4-1, a noise estimation method or a voice determination method can be used. The noise estimation method converts an input signal on the time axis into a signal on the frequency axis, and estimates background noise in the input signal while following the input signal based on the converted signal on the frequency axis. A method of forming an estimated noise signal (for example, Japanese Patent No. 0326969) may be used. That is, using the signal spectrum signal Xf, k of noise (background noise) that does not include speech and the estimated noise signal N ′ f, k−1 one analysis frame before, the estimated noise signal N ′ f, k of the current analysis frame is used. A method for calculating the value may be used. In this case, Xf, k and N′f, k−1 are multiplied by a coefficient to give a · N′f, k−1 and (1-a) · Xf, k, respectively, and a · N′f, k− 1 and (1-a) · Xf, k are added to calculate the estimated noise signal N′f, k of the current analysis frame. In addition, a speech determination method is a method for determining whether or not speech is included in each frame in which a speech signal on which environmental noise is superimposed as an input signal is given in chronological order (for example, Japanese Patent No. 03849116). ). In other words, the feature quantity that emphasizes the high frequency component of the audio signal feature quantity and the audio frequency feature quantity is calculated, and the determination is made based on the feature quantity that emphasizes the high frequency component in addition to the feature quantity of the audio signal. You can use the method of doing. Further, a method may be used in which the feature amount of the audio signal and the regularity of the resonance period of the audio are calculated and the determination is made based on the regularity of the resonance cycle of the audio in addition to the feature amount of the audio signal.

事前情報保持手段16は、新幹線など、経路やタイムテーブルが定まっている場合に、予め騒音に関する情報を取得し、この情報を事前情報12として保持する。   The prior information holding means 16 acquires information about noise in advance and holds this information as the prior information 12 when a route or a time table is fixed, such as a bullet train.

事前情報の一例およびその取得は次のとおりである。ここで、周波数特性とは、各周波数毎のパワースペクトルの大きさを示す数値〔dB〕である。   An example of prior information and its acquisition are as follows. Here, the frequency characteristic is a numerical value [dB] indicating the magnitude of the power spectrum for each frequency.

1)基準騒音の収録 1) Recording of standard noise

この基準騒音の収録では、基準となるエンジン回転数での騒音の周波数特性を測定する。この場合、基準のエンジン回転数、周波数f、基準のエンジン回転数の条件の騒音の周波数特性は次のとおりである。
ref :基準のエンジン回転数
ref (f):周波数f、基準のエンジン回転数の条件の騒音の周波数特性
f:周波数(Hz)
In recording the reference noise, the frequency characteristic of the noise at the reference engine speed is measured. In this case, the frequency characteristics of noise under the conditions of the reference engine speed, the frequency f, and the reference engine speed are as follows.
r ref : reference engine speed n ref (f): frequency characteristic of noise under the condition of frequency f and reference engine speed f: frequency (Hz)

2)エンジン回転数と騒音スペクトルの変化の関係の算出 2) Calculation of relationship between engine speed and noise spectrum change

このエンジン回転数と騒音スペクトルの変化の関係の算出では、エンジン回転数を変えた際の騒音周波数特性を測定し、式(1) の係数β(f)を算出する。   In calculating the relationship between the engine speed and the change in the noise spectrum, the noise frequency characteristic when the engine speed is changed is measured, and the coefficient β (f) of equation (1) is calculated.

Figure 2012215715
但し、r:エンジン回転数、
α(f,r):エンジン回転数r,周波数fの場合の、エンジン騒音の周波数成分を算出するための係数、
β(f):周波数fのエンジン回転数と基準のエンジン回転数の差によって変化する周波数成分の大きさの係数。
Figure 2012215715
Where r: engine speed,
α (f, r): coefficient for calculating the frequency component of engine noise in the case of engine speed r and frequency f,
β (f): A coefficient of the magnitude of the frequency component that changes depending on the difference between the engine speed at the frequency f and the reference engine speed.

式(2) により、エンジン回転数がrの際のエンジン音の周波数特性ne (f,r)を算出する。 The frequency characteristic n e (f, r) of the engine sound when the engine speed is r is calculated from the equation (2).

Figure 2012215715
但し、ne (f,r):エンジン回転数がrの場合のエンジンによる騒音周波数特性〔dB〕。
Figure 2012215715
Where ne (f, r): noise frequency characteristic [dB] by the engine when the engine speed is r.

この場合、β(f)は、予め周波数成分毎にエンジン回転数の基準からのずれ(r−rref )とエンジン騒音のパワー成分の基準からのずれを測定し、算出する。 In this case, β (f) is calculated in advance by measuring a deviation (r−r ref ) from the reference of the engine speed and a deviation from the reference of the power component of the engine noise for each frequency component.

3) 騒音の変化要因と、変化した際の騒音の収録 3) Noise change factors and recording of noise when changed

この騒音の変化要因と、変化した際の騒音の収録では、トンネル、橋、対向車両など、騒音の変化要因毎に予め騒音を収録する。そして、変化要因による騒音の変化量を算出する。   In recording the noise change factors and the noises at the time of the change, noise is recorded in advance for each noise change factor such as a tunnel, a bridge, and an oncoming vehicle. Then, the amount of change in noise due to the change factor is calculated.

トンネルの周波数fの騒音スペクトル変化量t(i,f)〔dB〕を式(3) により算出する。   The amount of noise spectrum change t (i, f) [dB] at the frequency f of the tunnel is calculated by equation (3).

Figure 2012215715
但し、t(i,f):i番目のトンネルの周波数fの騒音スペクトル変化量〔dB〕、
obs (i,f,t):i番目のトンネルの周波数fの時刻tにおける騒音スペ クトル〔dB〕、
T:スペクトルを算出した時間フレームのインデックス、
N:スペクトルを平均する時間フレームの数。
Figure 2012215715
Where t (i, f): noise spectrum change amount [dB] of the frequency f of the i-th tunnel,
t obs (i, f, t): noise spectrum [dB] at time t at frequency f of the i-th tunnel,
T: Index of the time frame in which the spectrum was calculated,
N: The number of time frames that average the spectrum.

同様に、式(4) および式(5) から、橋の騒音b(i,f)と対向車両の騒音o(i,f)も、同様に算出する。b(i,f)はi番目の橋の周波数fの騒音スペクトル変化量〔dB〕である。また、o(i,f)はi番目の対向車両の周波数fの騒音スペクトル変化量〔dB〕である。bobs (i,f,t)は、i番目の橋の周波数fの時刻tにおける騒音スペクトル〔dB〕であり、oobs (i,f,t)は、i番目の対向車両の周波数fの時刻tにおける騒音スペクトル〔dB〕である。 Similarly, the bridge noise b (i, f) and the oncoming vehicle noise o (i, f) are calculated in the same manner from the equations (4) and (5). b (i, f) is a noise spectrum change amount [dB] of the frequency f of the i-th bridge. Further, o (i, f) is a noise spectrum change amount [dB] of the frequency f of the i-th oncoming vehicle. b obs (i, f, t) is the noise spectrum [dB] at the time t at the frequency f of the i-th bridge, and o obs (i, f, t) is the frequency f of the i-th oncoming vehicle. It is a noise spectrum [dB] at time t.

Figure 2012215715
Figure 2012215715

Figure 2012215715
Figure 2012215715

現存特定騒音情報算出手段18は、日時、位置、エンジン回転数などの現在情報が与えられた際にたとえば、予め特定騒音リストを定めておき、この特定騒音リストから現存するかどうかの判定情報を参照し、現存特定騒音情報を算出する。   The existing specific noise information calculation means 18 determines, for example, a specific noise list in advance when current information such as date and time, position, and engine speed is given, and determines whether or not it exists from this specific noise list. The existing specific noise information is calculated with reference.

この現存特定騒音情報算出手段18では、現在情報(日時、位置、エンジン回転数)の入力により、予め定めた特定騒音との関係付けにより、橋・トンネルの走行、対向車両の有無を調べる。たとえば、橋・トンネルの走行、対向車両とすれ違う日時・位置を登録しておき、現在の日時・位置を照らし合わせることで、現存特定騒音源情報を生成する。   The existing specific noise information calculating means 18 checks the presence / absence of an oncoming vehicle or a bridge / tunnel by inputting the current information (date and time, position, engine speed) and associating with the predetermined specific noise. For example, the existing specific noise source information is generated by registering the date and position of a bridge / tunnel running and a passing vehicle, and comparing the current date and position.

現存特定騒音源情報を図3および図4に示す。図3に示す特定騒音情報テーブル21において、項目2102には「基準のエンジン回転数」、「トンネルの数」、「橋の数」、「対向車両の数」が設定され、これらの値2104は運行条件や路線に応じて予め設定する。特定騒音情報テーブル21の「現在のエンジン回転数」について、現存特定騒音情報算出手段18により外部から情報を入力して更新する。   The existing specific noise source information is shown in FIGS. In the specific noise information table 21 shown in FIG. 3, “reference engine speed”, “number of tunnels”, “number of bridges”, and “number of oncoming vehicles” are set in the item 2102. It is set in advance according to operation conditions and routes. The “current engine speed” in the specific noise information table 21 is updated by inputting information from the outside by the existing specific noise information calculation means 18.

また、図4に示す特定騒音情報テーブル22では、項目2202、2204、2206が設定され、項目2202には「有無」、項目2204には特定騒音種、項目2206にはインデックスが設定される。項目2202の「有無」には、現在の日時・位置を照らし合わせて得た、現存する特定騒音源情報によって、情報を入力して更新する。「有無」の" 1" はその特定騒音が存在することを示し、" 0" はその特定騒音が存在しないことを示す。   In the specific noise information table 22 shown in FIG. 4, items 2202, 2204, and 2206 are set, “existence” is set in the item 2202, specific noise type is set in the item 2204, and an index is set in the item 2206. In the “presence / absence” of the item 2202, information is input and updated by existing specific noise source information obtained by comparing the current date / time and position. “1” of “presence / absence” indicates that the specific noise exists, and “0” indicates that the specific noise does not exist.

たとえば、現在、5番目のトンネルを走行中で、7番目のポイントで対向車とすれ違い中の場合には、5番目のインデックスのトンネルの有無の欄を" 1" 、7番目のインデックスの対向車両の有無の欄を“1”、その他は全て“0”とする。     For example, if you are currently driving through the 5th tunnel and are passing the oncoming vehicle at the 7th point, set the 5th index tunnel presence / absence column to "1" and the 7th index oncoming vehicle. The column of presence / absence of “1” is “1”, and all others are “0”.

騒音推定手段8は、事前情報12と現存特定騒音情報を用いて騒音の周波数特性を算出し、第2の推定騒音周波数特性として出力する。この騒音推定手段8では、つぎの処理を行う。   The noise estimation means 8 calculates the frequency characteristic of the noise using the prior information 12 and the existing specific noise information, and outputs it as the second estimated noise frequency characteristic. The noise estimation means 8 performs the following processing.

1)現存特定騒音情報を入力する。   1) Input existing specific noise information.

2)基準の周囲騒音情報を得る。すなわち、事前情報より、基準のエンジン回転数に対する騒音スペクトルnref (f)を得る。 2) Obtain standard ambient noise information. That is, the noise spectrum n ref (f) with respect to the reference engine speed is obtained from the prior information.

3)エンジン回転数による騒音補正をする。すなわち、現存特定騒音情報よりエンジン回転数(r)を取得し、エンジン回転数の変化による騒音周波数特性の補正を行う。これは、既述の式(2) によって算出する。   3) The noise is corrected by the engine speed. That is, the engine speed (r) is acquired from the existing specific noise information, and the noise frequency characteristic is corrected by the change in the engine speed. This is calculated by the equation (2) described above.

4)位置・時刻による騒音補正をする。すなわち、現存特定騒音情報より、現存する騒音情報を得る。トンネル、橋、対向車両など、騒音の変化要因による騒音周波数特性の補正を行う。これは、式(6) によって算出する。   4) Perform noise correction based on position and time. That is, the existing noise information is obtained from the existing specific noise information. Correction of noise frequency characteristics due to noise change factors such as tunnels, bridges, and oncoming vehicles. This is calculated by equation (6).

Figure 2012215715
但し、Nt:トンネルのインデックスの数、
Nb:橋のインデックスの数、
No:対向車両のインデックスの数、
Et(i):現在インデックスiのトンネルを走行中の際は“1”、その他は“0”である関数。
Eb(i):現在インデックスiの橋を走行中の際は" 1" 、その他は“0”である関数。
Eo(i):現在インデックスiの対向車両とすれ違い中の際は“1”、その他は“0”である関数。
Figure 2012215715
Where Nt: number of tunnel indexes,
Nb: number of bridge indexes,
No: Number of oncoming vehicle index,
Et (i): A function that is “1” when currently traveling through a tunnel of index i, and “0” otherwise.
Eb (i): A function that is “1” when the vehicle is currently traveling on a bridge of index i, and “0” in other cases.
Eo (i): A function that is “1” when passing the oncoming vehicle of the current index i, and “0” otherwise.

騒音推定手段20は、既述の加算手段14(第1の実施の形態)の一例である。この騒音推定手段20は、第1の推定騒音周波数特性と第2の推定騒音周波数特性とを用いて、予め定めた寄与係数で重み付け加算することにより第3の推定騒音周波数特性を算出し、これを推定騒音周波数特性として出力する。つまり、第1の推定騒音周波数特性と第2の推定騒音周波数特性とに応じた推定騒音周波数特性が出力される。   The noise estimation means 20 is an example of the above-described addition means 14 (first embodiment). The noise estimation means 20 calculates a third estimated noise frequency characteristic by weighting and adding with a predetermined contribution coefficient using the first estimated noise frequency characteristic and the second estimated noise frequency characteristic. Is output as an estimated noise frequency characteristic. That is, an estimated noise frequency characteristic corresponding to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic is output.

この騒音推定手段20は、騒音推定手段4−1で算出した騒音の周波数特性と、騒音推定手段10で算出した騒音の周波数特性を用いて式(7) によって算出する。   The noise estimation unit 20 calculates the noise frequency characteristic calculated by the noise estimation unit 4-1 and the noise frequency characteristic calculated by the noise estimation unit 10 according to Expression (7).

Figure 2012215715
但し、coeff2:寄与係数(0〜1の数)、
n1_pred(f):第1の推定騒音周波数特性、
n2_pred(f):第2の推定騒音周波数特性、
n3_pred(f):第3の推定騒音周波数特性。
寄与係数は、予め与えた値または後述の寄与係数算出手段22−1(図6)から得ればよい。
Figure 2012215715
Where coeff2: contribution coefficient (number of 0 to 1),
n1_pred (f): first estimated noise frequency characteristic,
n2_pred (f): second estimated noise frequency characteristic,
n3_pred (f): third estimated noise frequency characteristic.
The contribution coefficient may be obtained from a value given in advance or contribution coefficient calculation means 22-1 (FIG. 6) described later.

つぎに、この騒音推定の処理手順について、図5を参照する。図5は騒音推定の処理手順の一例を示している。   Next, FIG. 5 is referred to regarding the processing procedure of the noise estimation. FIG. 5 shows an example of a noise estimation processing procedure.

図5に示す処理手順は本開示の騒音推定プログラムの一例である。この処理手順では、既述の第1の騒音推定手段4−1が入力音から第1の推定騒音周波数特性を算出する(ステップS11)。現存特定騒音情報算出手段18は、日時、位置、現在のエンジン回転数を表す現在情報から現存特定騒音情報を算出する(ステップS12)。第2の騒音推定手段8が現存特定騒音情報と事前情報とを用いて第2の推定騒音周波数特性を算出する(ステップS13)。そして、第3の騒音推定手段20が既述の第1の推定騒音周波数特性と第2の推定騒音周波数特性とを用いて第3の推定騒音周波数特性を算出し、推定騒音周波数特性として出力する(ステップS14)。   The processing procedure illustrated in FIG. 5 is an example of the noise estimation program of the present disclosure. In this processing procedure, the first noise estimation means 4-1 described above calculates a first estimated noise frequency characteristic from the input sound (step S11). The existing specific noise information calculation means 18 calculates the existing specific noise information from the current information indicating the date and time, the position, and the current engine speed (step S12). The second noise estimation means 8 calculates the second estimated noise frequency characteristic using the existing specific noise information and the prior information (step S13). Then, the third noise estimation means 20 calculates the third estimated noise frequency characteristic using the first estimated noise frequency characteristic and the second estimated noise frequency characteristic described above, and outputs it as the estimated noise frequency characteristic. (Step S14).

この実施の形態の騒音推定装置2−1では、騒音推定手段4によって騒音の周波数特性(第1の推定騒音周波数特性)を算出し、騒音推定手段8によって事前情報と現在情報を用いて騒音の周波数特性(第2の推定騒音周波数特性)を算出し、騒音推定手段20で第1および第2の推定騒音周波数特性を用いている。   In the noise estimation apparatus 2-1 of this embodiment, the noise estimation means 4 calculates the noise frequency characteristic (first estimated noise frequency characteristic), and the noise estimation means 8 uses the prior information and the current information to calculate the noise level. A frequency characteristic (second estimated noise frequency characteristic) is calculated, and the noise estimation means 20 uses the first and second estimated noise frequency characteristics.

この実施の形態によれば、入力音に利用者の音声が含まれて正しく騒音を推定できない場合や、トンネルに出入りした際など、急に騒音が変化して追従が難しい際に、第2の推定騒音周波数特性を使用することで、騒音を推定できる。従来に比べて、高性能な推定騒音周波数特性を算出できる。   According to this embodiment, when the user's voice is included in the input sound and the noise cannot be estimated correctly, or when the noise suddenly changes and it is difficult to follow, such as when entering or exiting the tunnel, the second Noise can be estimated by using the estimated noise frequency characteristic. Compared with the prior art, it is possible to calculate a high-performance estimated noise frequency characteristic.

〔第3の実施の形態〕   [Third Embodiment]

図6は第3の実施の形態に係る騒音推定装置の一例を示している。図6に示す構成は一例であり、斯かる構成に本発明が限定されるものではない。図6において、図2と同一部分には同一符号を付してある。   FIG. 6 shows an example of a noise estimation apparatus according to the third embodiment. The configuration shown in FIG. 6 is an example, and the present invention is not limited to such a configuration. In FIG. 6, the same parts as those in FIG.

図6に示す騒音推定装置2−2は、第1の騒音推定手段4−1と、事前情報保持手段16と、現存特定騒音情報算出手段18と、第2の騒音推定手段8と、第3の騒音推定手段20と、寄与係数算出手段22−1を備えている。騒音推定手段4−1、事前情報保持手段16、現存特定騒音情報算出手段18および騒音推定手段8は第2の実施の形態と同様である。   The noise estimation apparatus 2-2 shown in FIG. 6 includes a first noise estimation unit 4-1, a prior information holding unit 16, an existing specific noise information calculation unit 18, a second noise estimation unit 8, and a third noise estimation unit 2-2. Noise estimation means 20 and contribution coefficient calculation means 22-1. The noise estimation means 4-1, the prior information holding means 16, the existing specific noise information calculation means 18 and the noise estimation means 8 are the same as those in the second embodiment.

寄与係数算出手段22−1は、現存特定騒音情報を用いて、現存する騒音の種類によって、第1の推定騒音周波数特性に対する第2の推定騒音周波数特性の寄与の大きさを示す寄与係数を算出する。   The contribution coefficient calculation means 22-1 uses the existing specific noise information to calculate a contribution coefficient indicating the magnitude of the contribution of the second estimated noise frequency characteristic to the first estimated noise frequency characteristic depending on the type of existing noise. To do.

この寄与係数算出手段22−1では、現存する特定騒音種の係数を算出する。この場合、特定騒音種の係数は、図7に示すように、予め特定騒音情報テーブル23として保持する。この特定騒音種の係数を用いて、第2の騒音周波数特性の係数を式(8) で算出する。   The contribution coefficient calculation means 22-1 calculates the coefficient of the existing specific noise type. In this case, the coefficient of the specific noise type is held in advance as the specific noise information table 23 as shown in FIG. The coefficient of the second noise frequency characteristic is calculated by the equation (8) using the coefficient of the specific noise type.

Figure 2012215715
但し、coeff2_org:予め設定された第2の騒音周波数特性の係数、
coeff2:第2の騒音の周波数特性の寄与係数(0〜1の数)、
coeff_n(k):特定音源種の係数(インデックスk)、
特定音源種の係数は、現存する場合は特定騒音情報テーブル21の値、現存しない場合は“1.0" とする。
Nk:特定騒音の数
Figure 2012215715
Where coeff2_org: coefficient of the second noise frequency characteristic set in advance,
coeff2: contribution coefficient (number of 0 to 1) of the frequency characteristic of the second noise,
coeff_n (k): coefficient of specific sound source type (index k),
The coefficient of the specific sound source type is set to the value of the specific noise information table 21 when existing, or “1.0” when not existing.
Nk: Number of specific noise

そして、騒音推定手段20は、第2の実施の形態と同様に、第1の推定騒音周波数特性と第2の推定騒音周波数特性を寄与係数で重み付け加算することにより推定騒音周波数特性を算出する。つまりこの実施の形態においても、第1の推定騒音周波数特性と第2の推定騒音周波数特性とに応じた推定騒音周波数特性が出力される。   Then, as in the second embodiment, the noise estimation means 20 calculates the estimated noise frequency characteristic by weighting and adding the first estimated noise frequency characteristic and the second estimated noise frequency characteristic with the contribution coefficient. That is, also in this embodiment, an estimated noise frequency characteristic corresponding to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic is output.

つぎに、この騒音推定の処理手順について、図8を参照する。図8は騒音推定の処理手順の一例を示している。   Next, the noise estimation processing procedure will be described with reference to FIG. FIG. 8 shows an example of a noise estimation processing procedure.

図8に示す処理手順は本開示の騒音推定プログラムの一例である。この処理手順では、第1の騒音推定手段4−1が入力音から第1の推定騒音周波数特性を算出する(ステップS21)。現存特定騒音情報算出手段18は現在情報から現存特定騒音情報を算出する(ステップS22)。第2の騒音推定手段8は現存特定騒音情報と事前情報とから第2の推定騒音周波数特性を算出する(ステップS23)。寄与係数算出手段22−1は現存特定騒音情報を入力し、寄与係数を算出し、出力する(ステップS24)。そして、第3の騒音推定手段20は第1の推定騒音周波数特性と第2の推定騒音周波数特性と寄与係数とから第3の推定騒音周波数特性を算出し、出力する(ステップS25)。   The processing procedure shown in FIG. 8 is an example of the noise estimation program of the present disclosure. In this processing procedure, the first noise estimation unit 4-1 calculates a first estimated noise frequency characteristic from the input sound (step S21). The existing specific noise information calculation means 18 calculates existing specific noise information from the current information (step S22). The second noise estimation means 8 calculates a second estimated noise frequency characteristic from the existing specific noise information and the prior information (step S23). The contribution coefficient calculation means 22-1 inputs the existing specific noise information, calculates a contribution coefficient, and outputs it (step S24). Then, the third noise estimation means 20 calculates and outputs a third estimated noise frequency characteristic from the first estimated noise frequency characteristic, the second estimated noise frequency characteristic, and the contribution coefficient (step S25).

この騒音推定装置2−2によれば、寄与係数算出手段22−1で、現存する特定騒音によって寄与係数を算出する。第1の実施の形態の効果に加え、音源の種類によってエンジン回転数が変化した際や、対向車両とのすれ違いなどの、急に騒音が変化する場合、第2の推定騒音周波数特性の寄与度を大きくすることができる。この結果、騒音変化の追従性を速めることができ、より騒音の推定性能を高めることができる。   According to this noise estimation apparatus 2-2, the contribution coefficient calculation means 22-1 calculates the contribution coefficient using the existing specific noise. In addition to the effects of the first embodiment, when the engine speed changes depending on the type of sound source, or when noise suddenly changes, such as when passing the oncoming vehicle, the contribution of the second estimated noise frequency characteristic Can be increased. As a result, the follow-up performance of the noise change can be accelerated, and the noise estimation performance can be further improved.

〔第4の実施の形態〕   [Fourth Embodiment]

図9は第4の実施の形態に係る騒音推定装置の一例を示している。図9に示す構成は一例であり、斯かる構成に本発明が限定されるものではない。図9において、図2と同一部分には同一符号を付してある。   FIG. 9 shows an example of a noise estimation apparatus according to the fourth embodiment. The configuration shown in FIG. 9 is an example, and the present invention is not limited to such a configuration. 9, the same parts as those in FIG. 2 are denoted by the same reference numerals.

図9に示す騒音推定装置2−3は、第1の騒音推定手段4−2と、事前情報保持手段16と、現存特定騒音情報算出手段18と、第2の騒音推定手段8と、第3の騒音推定手段20と、事前情報更新手段24を備えている。現存特定騒音情報算出手段18、騒音推定手段8および第3の騒音推定手段20は第2の実施の形態と同様である。   The noise estimation apparatus 2-3 shown in FIG. 9 includes a first noise estimation unit 4-2, a prior information holding unit 16, an existing specific noise information calculation unit 18, a second noise estimation unit 8, and a third noise estimation unit. Noise estimating means 20 and prior information updating means 24 are provided. The existing specific noise information calculation means 18, the noise estimation means 8, and the third noise estimation means 20 are the same as those in the second embodiment.

第1の騒音推定手段4−2は、この騒音の周波数特性(第1の推定騒音周波数特性)と、現在の周波数特性と、音声判定結果を算出する。第1の推定騒音周波数特性は、入力音から騒音の周波数特性を算出して得ることができる。入力音には音声が含まれる入力音が含まれる。音声判定結果の算出は、入力音に音声が含まれているか否かの判定である。この騒音推定手段4−2の処理には、騒音推定の方法や音声判定の方法を用いることができる。騒音推定の方法は、既述の時間軸上の入力信号を周波数軸上の信号に変換し、変換された周波数軸上の信号に基づいて、その入力信号における背景騒音を入力信号に追従しながら推定し、推定雑音信号を形成する方法(たとえば、特許第3269969号)を利用すればよい。即ち、音声を含まない雑音(背景雑音)の信号スペクトル信号Xf,kと1分析フレーム前の推定雑音信号N' f,k−1を用いて、現分析フレームの推定雑音信号N' f,kを算出する方法を利用すればよい。この場合、Xf,kとN' f,k−1に係数を掛けてそれぞれをa・N' f,k−1及び(1−a)・Xf,kとし、a・N' f,k−1と(1−a)・Xf,kを加算して、現分析フレームの推定雑音信号N' f,kが算出される。音声判定の方法は、既述の入力信号として環境騒音が重畳した音声信号が時系列順に与えられるフレーム毎に、当該フレームに音声が含まれているか否かを判定する方法(たとえば、特許第3849116号)を利用すればよい。即ち、音声信号の特徴量と、音声周波数の特徴量の高域成分を強調させた特徴量を算出して、音声信号の特徴量に加えて高域成分を強調させた特徴量を基に判定を行う方法を利用すれば良い。また、音声信号の特徴量と、音声の共振周期の規則性を算出して、音声信号の特徴量に加えて音声の共振周期の規則性を基に判定を行う方法を利用すれば良い。   The first noise estimation means 4-2 calculates the frequency characteristic of the noise (first estimated noise frequency characteristic), the current frequency characteristic, and the voice determination result. The first estimated noise frequency characteristic can be obtained by calculating the frequency characteristic of noise from the input sound. The input sound includes an input sound including sound. The calculation of the sound determination result is a determination of whether or not sound is included in the input sound. For the processing of the noise estimation unit 4-2, a noise estimation method or a voice determination method can be used. The noise estimation method converts the input signal on the time axis described above into a signal on the frequency axis, and follows the background noise in the input signal based on the converted signal on the frequency axis while following the input signal. A method of estimating and forming an estimated noise signal (for example, Japanese Patent No. 3269969) may be used. That is, using the signal spectrum signal Xf, k of noise (background noise) that does not include speech and the estimated noise signal N ′ f, k−1 one analysis frame before, the estimated noise signal N ′ f, k of the current analysis frame is used. A method for calculating the value may be used. In this case, Xf, k and N′f, k−1 are multiplied by a coefficient to give a · N′f, k−1 and (1-a) · Xf, k, respectively, and a · N′f, k− 1 and (1-a) · Xf, k are added to calculate the estimated noise signal N′f, k of the current analysis frame. The voice determination method is a method for determining whether or not a voice is included in the frame for each frame in which a voice signal on which environmental noise is superimposed as an input signal described above is given in chronological order (for example, Japanese Patent No. 3849116). No.) can be used. In other words, the feature quantity that emphasizes the high frequency component of the audio signal feature quantity and the audio frequency feature quantity is calculated, and the determination is made based on the feature quantity that emphasizes the high frequency component in addition to the feature quantity of the audio signal. You can use the method of doing. Further, a method may be used in which the feature amount of the audio signal and the regularity of the resonance period of the audio are calculated and the determination is made based on the regularity of the resonance cycle of the audio in addition to the feature amount of the audio signal.

事前情報更新手段24によって現存特定騒音情報と音声判定結果と現在の周波数特性を用いて、事前情報保持手段16の事前情報12を更新する。第2の騒音推定手段8は、事前情報更新手段24によって更新された事前情報12と、現存特定騒音情報とを用いて騒音の周波数特性(第2の推定騒音周波数特性)を算出する。そして、第3の騒音推定手段20は、第1の推定騒音周波数特性と第2の推定騒音周波数特性を、予め定めた寄与係数で重み付け加算することで推定騒音周波数特性を算出する。つまり、第3の騒音推定手段20では、第1の推定騒音周波数特性と第2の推定騒音周波数特性とに応じた推定騒音周波数特性出力が得られる。   The prior information updating unit 24 updates the prior information 12 of the prior information holding unit 16 using the existing specific noise information, the voice determination result, and the current frequency characteristic. The second noise estimation means 8 calculates noise frequency characteristics (second estimated noise frequency characteristics) using the prior information 12 updated by the prior information update means 24 and the existing specific noise information. Then, the third noise estimation unit 20 calculates the estimated noise frequency characteristic by weighting and adding the first estimated noise frequency characteristic and the second estimated noise frequency characteristic with a predetermined contribution coefficient. That is, the third noise estimation means 20 can obtain an estimated noise frequency characteristic output corresponding to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic.

つぎに、事前情報更新の処理について、図10を参照する。図10は事前情報更新の処理手順を示している。   Next, FIG. 10 will be referred to regarding the advance information update process. FIG. 10 shows a processing procedure for updating prior information.

この処理手順は、事前情報更新手段24で実行される。図10に示す処理手順では、非音声区間であるかを判定する(ステップS31)。非音声区間であれば(ステップS31のYES)、特定音源騒音種の数を算出する(ステップS32)。   This processing procedure is executed by the prior information update unit 24. In the processing procedure shown in FIG. 10, it is determined whether it is a non-voice section (step S31). If it is a non-voice section (YES in step S31), the number of specific sound source noise types is calculated (step S32).

算出された特定音源騒音種の数について、特定音源騒音種の数=0であるかを判定する(ステップS33)。算出した特定騒音種の数=0の場合、エンジン回転による騒音以外の特定騒音種がないことを表している。特定音源騒音種の数=0であれば(ステップS33のYES)、エンジン回転数から基準のエンジン回転数を減算した値の絶対値が閾値(TH)より小さいかすなわち、|エンジン回転数−基準のエンジン回転数|<THであるかを判定する(ステップS34)。   For the calculated number of specific sound source noise types, it is determined whether the number of specific sound source noise types = 0 (step S33). When the calculated number of specific noise types = 0, it indicates that there are no specific noise types other than noise caused by engine rotation. If the number of specific sound source noise types = 0 (YES in step S33), whether the absolute value of the value obtained by subtracting the reference engine speed from the engine speed is smaller than the threshold (TH), that is, | engine speed-reference It is determined whether or not the engine rotational speed | <TH (step S34).

|エンジン回転数−基準のエンジン回転数|<THであれば(ステップS34のYES)、基準のエンジン回転数の周波数特性を更新し(ステップS35)、この処理を終了する。ステップS34において、|エンジン回転数−基準のエンジン回転数|<THでなければ(ステップS34のNO)、この処理を終了する。   If | engine speed−reference engine speed | <TH (YES in step S34), the frequency characteristic of the reference engine speed is updated (step S35), and this process ends. In step S34, if | engine speed−reference engine speed | <TH is not satisfied (NO in step S34), the process is terminated.

基準のエンジン回転数の周波数特性の更新は式(9) を用いて行う。   Update the frequency characteristics of the standard engine speed using Equation (9).

Figure 2012215715
但し、coeff:0〜1の間の数(たとえば、0.9)、
nref(f):更新後の基準のエンジン回転数の周波数特性〔dB〕、
nref_prev(f):更新前の基準のエンジン回転数の周波数特性〔dB〕、
n_cur(f):現在の周波数特性〔dB〕。
Figure 2012215715
However, coeff: a number between 0 and 1 (for example, 0.9),
nref (f): frequency characteristics [dB] of the reference engine speed after the update,
nref_prev (f): frequency characteristic [dB] of the reference engine speed before update,
n_cur (f): Current frequency characteristic [dB].

また、ステップS33において、特定音源騒音種の数=0でなければ(ステップS33のNO)、特定音源騒音種の数=1であるかを判定する(ステップS36)。特定音源騒音種の数=1であれば(ステップS36のYES)、特定音源騒音の周波数特性を更新し(ステップS37)、この処理を終了する。特定音源騒音種の数=1でなければ(ステップS36のNO)、特定音源騒音の周波数特性を更新しないで、この処理を終了する。   In step S33, if the number of specific sound source noise types is not 0 (NO in step S33), it is determined whether the number of specific sound source noise types is 1 (step S36). If the number of specific sound source noise types = 1 (YES in step S36), the frequency characteristics of the specific sound source noise are updated (step S37), and this process ends. If the number of specific sound source noise types is not 1 (NO in step S36), the processing is terminated without updating the frequency characteristics of the specific sound source noise.

現存の特定騒音種の周波数特性は、次の式(10)、(11)または(12)を用いて現存の特定騒音種(1つ)を更新する。   For the frequency characteristics of the existing specific noise type, the existing specific noise type (one) is updated using the following formula (10), (11) or (12).

Figure 2012215715
Figure 2012215715

Figure 2012215715
Figure 2012215715

Figure 2012215715
但し、coeff:0〜1の間の数(たとえば、0.9)
t(i,f):i番目のトンネルの周波数fの騒音周波数特性〔dB〕
t_prev(i,f):更新前のi番目のトンネルの周波数fの騒音周波数特性〔dB〕
b(i,f):i番目の橋の周波数fの騒音周波数特性〔dB〕
b_prev(i,f):更新前のi番目の橋の周波数fの騒音周波数特性〔dB〕
o(i,f):i番目の対向車両の周波数fの騒音周波数特性〔dB〕
o_prev(i,f):更新前のi番目の対向車両の周波数fの騒音周波数特性〔dB〕
n_cur2(f):現在の周波数特性〔dB〕から、エンジン音の成分を除去した周波数fの騒音周波数特性〔dB〕
n_cur2(f)は、式(13)に示す通りである。
Figure 2012215715
However, coeff: a number between 0 and 1 (for example, 0.9)
t (i, f): noise frequency characteristic of the frequency f of the i-th tunnel [dB]
t_prev (i, f): noise frequency characteristic [dB] of the frequency f of the i-th tunnel before update
b (i, f): noise frequency characteristic of the frequency f of the i-th bridge [dB]
b_prev (i, f): noise frequency characteristic [dB] of the frequency f of the i-th bridge before update
o (i, f): noise frequency characteristic [dB] of the frequency f of the i-th oncoming vehicle
o_prev (i, f): noise frequency characteristic [dB] of the frequency f of the i-th oncoming vehicle before update
n_cur2 (f): Noise frequency characteristic [dB] of the frequency f obtained by removing the engine sound component from the current frequency characteristic [dB].
n_cur2 (f) is as shown in Expression (13).

Figure 2012215715
Figure 2012215715

つぎに、この騒音推定の処理手順について、図11を参照する。図11は騒音推定の処理手順の一例を示している。   Next, FIG. 11 is referred to regarding the processing procedure of the noise estimation. FIG. 11 shows an example of a processing procedure for noise estimation.

図11に示す処理手順は本開示の騒音推定プログラムの一例である。この処理手順では、第1の騒音推定手段4−2は入力音から現在の周波数特性を算出する(ステップS41)。第1の騒音推定手段4−2は入力音から第1の推定騒音周波数特性を算出する(ステップS42)。第1の騒音推定手段4−2は入力音の音声判定を行う(ステップS43)。   The processing procedure shown in FIG. 11 is an example of the noise estimation program of the present disclosure. In this processing procedure, the first noise estimation unit 4-2 calculates the current frequency characteristic from the input sound (step S41). The first noise estimation means 4-2 calculates a first estimated noise frequency characteristic from the input sound (step S42). The first noise estimation unit 4-2 performs voice determination of the input sound (step S43).

現存特定騒音情報算出手段18は現在情報から現存特定騒音情報を算出する(ステップS44)。事前情報更新手段24は現在の周波数特性と入力音の音声判定結果と現存特定騒音情報を入力し、事前情報保持手段16にある事前情報12を更新する(ステップS45)。   The existing specific noise information calculation means 18 calculates existing specific noise information from the current information (step S44). The prior information update unit 24 inputs the current frequency characteristics, the voice determination result of the input sound, and the existing specific noise information, and updates the prior information 12 in the prior information holding unit 16 (step S45).

第2の騒音推定手段8は現存特定騒音情報と事前情報12から第2の推定騒音周波数特性を算出する(ステップS46)。   The second noise estimation means 8 calculates a second estimated noise frequency characteristic from the existing specific noise information and the prior information 12 (step S46).

第3の騒音推定手段20は、第2の実施の形態と同様に、第1の推定騒音周波数特性と第2の推定騒音周波数特性を入力し、第3の推定騒音周波数特性を算出し、出力する(ステップS47)。   Similarly to the second embodiment, the third noise estimation means 20 inputs the first estimated noise frequency characteristic and the second estimated noise frequency characteristic, calculates the third estimated noise frequency characteristic, and outputs it. (Step S47).

この騒音推定装置2−3によれば、事前情報更新手段24は、予め保持した騒音と実際の騒音が異なる際に、実際の騒音に近づけるように事前情報12に保持されている騒音情報を更新できるため、より騒音の推定性能を高めることができる。   According to the noise estimation apparatus 2-3, the prior information update unit 24 updates the noise information held in the prior information 12 so as to approach the actual noise when the previously held noise differs from the actual noise. Therefore, the noise estimation performance can be further improved.

〔第5の実施の形態〕   [Fifth Embodiment]

図12は第5の実施の形態に係る騒音推定装置の一例を示している。図12に示す構成は一例であり、斯かる構成に本発明が限定されるものではない。   FIG. 12 shows an example of a noise estimation apparatus according to the fifth embodiment. The configuration shown in FIG. 12 is an example, and the present invention is not limited to such a configuration.

図12に示す騒音推定装置2−4は、第1の騒音推定手段4−2と、事前情報保持手段16と、現存特定騒音情報算出手段18と、第2の騒音推定手段8と、第3の騒音推定手段20と、寄与係数算出手段22−2と、信頼度算出手段26とを備えている。   The noise estimation device 2-4 shown in FIG. 12 includes a first noise estimation unit 4-2, a prior information holding unit 16, an existing specific noise information calculation unit 18, a second noise estimation unit 8, and a third noise estimation unit 4-2. Noise estimation means 20, contribution coefficient calculation means 22-2, and reliability calculation means 26.

第1の騒音推定手段4−2は、入力音から既述の方法で騒音の周波数特性(第1の推定騒音周波数特性)の算出、音声判定、現在の周波数特性の算出をする。音声判定は、入力音に音声が含まれるか否かの判定である。この第1の騒音推定手段4−2は、第4の実施の形態と同様であるので、具体的な処理についてはその説明を割愛する。   The first noise estimation means 4-2 calculates the frequency characteristics (first estimated noise frequency characteristics) of the noise, the voice determination, and the current frequency characteristics from the input sound by the method described above. The sound determination is a determination as to whether or not sound is included in the input sound. Since this first noise estimation means 4-2 is the same as that of the fourth embodiment, the description of the specific processing is omitted.

事前情報保持手段16は、新幹線など、経路やタイムテーブルが定まっている場合に、予め取得した騒音に関する事前情報12を保持している。   The prior information holding means 16 holds the prior information 12 related to noise acquired in advance when a route or a time table is fixed, such as a bullet train.

現存特定騒音情報算出手段18は、現在情報として日時、位置、エンジン回転数が与えられた際に、現存特定騒音情報を算出する。現存特定騒音情報は、予め定められた特定騒音リストについて、現存するかどうかの判定情報である。   The existing specific noise information calculating means 18 calculates the existing specific noise information when the date and time, the position, and the engine speed are given as the current information. The existing specific noise information is determination information as to whether or not a predetermined specific noise list exists.

信頼度算出手段26は、現存特定騒音情報と音声判定結果と現在の周波数特性を用いて事前情報の信頼度情報を算出する。信頼度情報は、予め定められた特定騒音リストについて、各々の特定騒音の信頼度を数値で表現された情報である。信頼度を表す数値は、たとえば、0:信頼度最小から1:信頼度最大の値である。   The reliability calculation means 26 calculates the reliability information of the prior information using the existing specific noise information, the voice determination result, and the current frequency characteristic. The reliability information is information in which the reliability of each specific noise is expressed by a numerical value for a predetermined specific noise list. The numerical value representing the reliability is, for example, a value from 0: minimum reliability to 1: maximum reliability.

第2の騒音推定手段8は、事前情報と現存特定騒音情報を用いて騒音の周波数特性(第2の推定騒音周波数特性)を算出する。   The second noise estimation means 8 calculates noise frequency characteristics (second estimated noise frequency characteristics) using the prior information and the existing specific noise information.

寄与係数算出手段22−2は、現存特定騒音情報と信頼度情報を用いて、寄与係数を算出する。この寄与係数は、第1の推定騒音周波数特性に対する第2の推定騒音周波数特性を寄与させる大きさを表わす情報である。   The contribution coefficient calculation means 22-2 calculates a contribution coefficient using the existing specific noise information and the reliability information. This contribution coefficient is information representing the magnitude of the contribution of the second estimated noise frequency characteristic to the first estimated noise frequency characteristic.

第3の騒音推定手段20は、第3の実施の形態と同様に、第1の推定騒音周波数特性と第2の推定騒音周波数特性を、寄与係数で重み付け加算し、推定騒音周波数特性を算出する。つまり、騒音推定手段20では第1の推定騒音周波数特性と第2の推定騒音周波数特性とに応じて推定騒音周波数特性出力が得られる。   Similarly to the third embodiment, the third noise estimation means 20 calculates the estimated noise frequency characteristic by weighting and adding the first estimated noise frequency characteristic and the second estimated noise frequency characteristic with the contribution coefficient. . That is, the noise estimation means 20 can obtain an estimated noise frequency characteristic output according to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic.

つぎに、信頼度算出手段26について、図13を参照する。図13は信頼度算出手段26の信頼度算出の処理手順の一例を示している。   Next, FIG. 13 is referred to for the reliability calculation means 26. FIG. 13 shows an example of a processing procedure of reliability calculation by the reliability calculation means 26.

この信頼度算出には、現存特定騒音情報、音声判定および現在の周波数特性を使用する。音声判定には第1の騒音推定手段4−2の判定結果を用いる。   For the reliability calculation, the existing specific noise information, voice determination, and current frequency characteristics are used. The determination result of the first noise estimation means 4-2 is used for the sound determination.

図13に示す処理手順では、非音声区間か否かを判定する(ステップS51)。第1の騒音推定手段4−2の音声判定結果が、非音声区間であれば(ステップS51のYES)、現存特定騒音情報の、特定騒音種の数を算出する(ステップS52)。音声判定結果が、非音声区間でなければ(ステップS51のNO)、この処理を終了する。   In the processing procedure shown in FIG. 13, it is determined whether or not it is a non-voice section (step S51). If the voice determination result of the first noise estimation unit 4-2 is a non-voice section (YES in step S51), the number of specific noise types in the existing specific noise information is calculated (step S52). If the speech determination result is not a non-speech interval (NO in step S51), this process ends.

算出した特定騒音種の数=0かを判定する(ステップS53)。算出した特定騒音種の数=0の場合、エンジン回転による騒音以外の特定騒音種がないことを表している。特定騒音種の数=0であれば(ステップS53のYES)、エンジン回転数から基準のエンジン回転数を減算した値の絶対値が閾値(TH)より小さいかすなわち、|エンジン回転数−基準のエンジン回転数|<THであるかを判定する(ステップS54)。   It is determined whether the calculated number of specific noise types = 0 (step S53). When the calculated number of specific noise types = 0, it indicates that there are no specific noise types other than noise caused by engine rotation. If the number of specific noise types = 0 (YES in step S53), whether the absolute value of the value obtained by subtracting the reference engine speed from the engine speed is smaller than the threshold (TH), that is, | engine speed-reference It is determined whether the engine speed | <TH (step S54).

|エンジン回転数−基準のエンジン回転数|<THであれば(ステップS54のYES)、基準のエンジン回転数の周波数特性の信頼度を算出し(ステップS55)、この処理を終了する。|エンジン回転数−基準のエンジン回転数|<THでなければ(ステップS54のNO)、基準のエンジン回転数の周波数特性の信頼度を算出することなく、この処理を終了する。   If | engine speed−reference engine speed | <TH (YES in step S54), the reliability of the frequency characteristic of the reference engine speed is calculated (step S55), and this process ends. If | engine speed−reference engine speed | <TH (NO in step S54), the process ends without calculating the reliability of the frequency characteristics of the reference engine speed.

〔基準のエンジン回転数の周波数特性の信頼度の算出〕 [Calculation of reliability of frequency characteristics of standard engine speed]

ステップS55における基準のエンジン回転数の周波数特性の信頼度は、式(14)及び(15)より算出する。   The reliability of the frequency characteristic of the reference engine speed in step S55 is calculated from equations (14) and (15).

Figure 2012215715
Figure 2012215715

Figure 2012215715
但し、nref_r(f):エンジン回転数が基準の場合のエンジンによる騒音周波数 特性の信頼度(0〜1)、
Δ(f):基準のエンジン回転数の騒音周波数特性の推定値と事前情 報の誤差〔dB〕、
TH3:小さな値(5〔dB〕など)。
Figure 2012215715
However, nref_r (f): reliability of the noise frequency characteristic by the engine when the engine speed is a reference (0 to 1),
Δ (f): Estimated value of noise frequency characteristic of standard engine speed and error [dB] of prior information,
TH3: A small value (5 [dB] or the like).

また、ステップS52で算出された特定騒音種の数が特定騒音種の数=0でなければ(ステップS53のNO)、特定騒音種の数=1であるかを判定する(ステップS56)。特定騒音種の数=1であれば(ステップS56のYES)、特定音源騒音の信頼度を算出し(ステップS57)、この処理を終了する。特定騒音種の数=1でなければ(ステップS56のNO)、特定音源騒音の信頼度を算出することなく、この処理を終了する。   If the number of specific noise types calculated in step S52 is not zero (NO in step S53), it is determined whether the number of specific noise types = 1 (step S56). If the number of specific noise types = 1 (YES in step S56), the reliability of specific sound source noise is calculated (step S57), and this process is terminated. If the number of specific noise types is not 1 (NO in step S56), the process ends without calculating the reliability of the specific sound source noise.

〔現存の特定騒音種の信頼度の算出〕 [Calculation of reliability of existing specific noise types]

ステップS57における、現存の特定騒音種の信頼度の算出は、式(16)、式(17)を用いて行う。   The calculation of the reliability of the existing specific noise type in step S57 is performed using Equation (16) and Equation (17).

Figure 2012215715
Figure 2012215715

Figure 2012215715
但し、nref_r(f):エンジン回転数が基準の場合のエンジンによる騒音周波数 特性の信頼度(0〜1)、
Δt(i,f):インデックスiのトンネルの騒音周波数特性の推定値と事 前情報の誤差〔dB〕。
Figure 2012215715
However, nref_r (f): reliability of the noise frequency characteristic by the engine when the engine speed is a reference (0 to 1),
Δt (i, f): an error [dB] between the estimated value of the noise frequency characteristic of the tunnel of index i and the prior information.

現存の特定騒音種が橋、対向車両の場合であれば、既述のトンネルの場合と同様に式(16)、式(17)を用いて信頼度を算出すればよい。     If the existing specific noise type is a bridge or an oncoming vehicle, the reliability may be calculated using Equations (16) and (17) as in the case of the tunnel described above.

このように事前情報に保持されている現存する特定騒音種の数と、現存する特定騒音種の信頼度を得る。   Thus, the number of existing specific noise types held in the prior information and the reliability of the existing specific noise types are obtained.

寄与係数算出手段22−2では、特定騒音種の信頼度を用いて、第2の騒音周波数特性の係数を式(18)により算出する。   In the contribution coefficient calculation means 22-2, the coefficient of the second noise frequency characteristic is calculated by the equation (18) using the reliability of the specific noise type.

Figure 2012215715
但し、coeff2=coeff2_org * rel =1.0
coeff2_org *rel ≦1
coeff2_org *rel >1
coeff2_org :予め設定された寄与係数
rel:現在の特定騒音種の信頼度
coeff2:第2の騒音の周波数特性の寄与係数(0〜1の数)
Nt:トンネルのインデックスの数
Nb:橋のインデックスの数
No:対向車両のインデックスの数
Et(i):現在インデックスiのトンネルを走行中の際は“1”、その他は“0”である関数。
Eb(i):現在インデックスiの橋を走行中の際は“1”、その他は“0”である関数。
Eo(i):現在インデックスiの対向車両とすれ違い中の際は“1”、その他は“0”である関数。
Ne:現存する特定騒音種の数
F:周波数の帯域数
nref_r(f):エンジン回転数が基準の場合のエンジンによる騒音スペクトルの信頼度(0〜1)
t_r(f):インデックスiのトンネル走行の場合の騒音スペクトルの信頼度(0〜1)
b_r(f):インデックスiの橋走行の場合の騒音スペクトルの信頼度(0〜1)
o_r(f):インデックスiの対向車両の場合の騒音スペクトルの信頼度(0〜1)
Figure 2012215715
However, coeff2 = coeff2_org * rel = 1.0
coeff2_org * rel ≦ 1
coeff2_org * rel> 1
coeff2_org: contribution coefficient set in advance
rel: The reliability of the current specific noise type
coeff2: contribution factor (number of 0 to 1) of the frequency characteristic of the second noise
Nt: Number of tunnel indexes
Nb: Number of bridge indexes
No: Number of oncoming vehicle indexes
Et (i): A function that is “1” when currently traveling through a tunnel of index i, and “0” otherwise.
Eb (i): A function that is “1” when the vehicle is currently running on the bridge of index i, and “0” in other cases.
Eo (i): A function that is “1” when passing the oncoming vehicle of the current index i, and “0” otherwise.
Ne: Number of existing specific noise types
F: Number of frequency bands
nref_r (f): reliability of the noise spectrum by the engine when the engine speed is a reference (0 to 1)
t_r (f): reliability of noise spectrum in the case of tunnel traveling at index i (0 to 1)
b_r (f): reliability of noise spectrum in the case of traveling on bridge with index i (0 to 1)
o_r (f): reliability of noise spectrum in the case of oncoming vehicle with index i (0 to 1)

つぎに、この騒音推定の処理手順について、図14を参照する。図14は騒音推定の処理手順の一例を示している。   Next, FIG. 14 is referred to regarding the processing procedure of this noise estimation. FIG. 14 shows an example of a processing procedure for noise estimation.

図14に示す処理手順は本開示の騒音推定プログラムの一例である。この処理手順では第1の騒音推定手段4−2が入力音から現在の周波数特性を算出する(ステップS61)。第1の騒音推定手段4−2は入力音から第1の推定騒音周波数特性を算出する(ステップS62)。第1の騒音推定手段4−2は入力音の音声判定を行う(ステップS63)。   The processing procedure illustrated in FIG. 14 is an example of the noise estimation program of the present disclosure. In this processing procedure, the first noise estimation means 4-2 calculates the current frequency characteristic from the input sound (step S61). The first noise estimation means 4-2 calculates a first estimated noise frequency characteristic from the input sound (step S62). The first noise estimation means 4-2 performs voice determination of the input sound (step S63).

現存特定騒音情報算出手段18が現在情報から現存特定騒音情報を算出する(ステップS64)。信頼度算出手段26が第1の推定騒音周波数特性と入力音の音声判定結果と現存特定騒音情報から事前情報の信頼度(信頼度情報)を算出する(ステップS65)。   The existing specific noise information calculation means 18 calculates existing specific noise information from the current information (step S64). The reliability calculation means 26 calculates the reliability (reliability information) of the prior information from the first estimated noise frequency characteristic, the voice determination result of the input sound, and the existing specific noise information (step S65).

寄与係数算出手段22−2は現存特定騒音情報と信頼度情報とから寄与係数を算出する(ステップS66)。第2の騒音推定手段8は現存特定騒音情報と事前情報12とから第2の推定騒音周波数特性を算出する(ステップS67)。そして、第3の騒音推定手段20は第1の推定騒音周波数特性と第2の推定騒音周波数特性と寄与係数を入力し、第3の推定騒音周波数特性を算出し、出力する(ステップS68)。   The contribution coefficient calculation means 22-2 calculates a contribution coefficient from the existing specific noise information and the reliability information (step S66). The second noise estimation means 8 calculates a second estimated noise frequency characteristic from the existing specific noise information and the prior information 12 (step S67). The third noise estimation means 20 receives the first estimated noise frequency characteristic, the second estimated noise frequency characteristic, and the contribution coefficient, calculates the third estimated noise frequency characteristic, and outputs it (step S68).

この騒音推定装置2−4によれば、信頼度算出手段26が予め保持した騒音と実際の騒音との差によって事前情報12の信頼度を算出し、寄与係数算出手段22−2が音源の信頼度によって、第2の推定騒音周波数特性の寄与係数を算出するので、騒音の推定性能をより高めることができる。   According to this noise estimation apparatus 2-4, the reliability of the prior information 12 is calculated based on the difference between the noise held in advance by the reliability calculation means 26 and the actual noise, and the contribution coefficient calculation means 22-2 is trusted by the sound source. Since the contribution coefficient of the second estimated noise frequency characteristic is calculated according to the degree, the noise estimation performance can be further improved.

〔第6の実施の形態〕   [Sixth Embodiment]

図15は第6の実施の形態に係る騒音推定装置の一例を示している。図15に示す構成は一例であり、斯かる構成に本発明が限定されるものではない。   FIG. 15 shows an example of a noise estimation apparatus according to the sixth embodiment. The configuration shown in FIG. 15 is an example, and the present invention is not limited to such a configuration.

図15に示す騒音推定装置2−5は、第1の騒音推定手段4−3と、事前情報保持手段16と、現存特定騒音情報算出手段18と、第2の騒音推定手段8と、第3の騒音推定手段20と、寄与係数算出手段22−3を備えている。   A noise estimation device 2-5 shown in FIG. 15 includes a first noise estimation unit 4-3, a prior information holding unit 16, an existing specific noise information calculation unit 18, a second noise estimation unit 8, and a third noise estimation unit 8. Noise estimation means 20 and contribution coefficient calculation means 22-3.

第1の騒音推定手段4−3は、入力音から既述の方法により騒音の周波数特性(第1の推定騒音周波数特性)の算出、音声判定を行う。   The first noise estimation means 4-3 calculates noise frequency characteristics (first estimated noise frequency characteristics) from the input sound according to the method described above, and performs voice determination.

この騒音推定手段4−3の処理には、騒音推定の方法や音声判定の方法を用いることができる。騒音推定の方法は、既述の時間軸上の入力信号を周波数軸上の信号に変換し、変換された周波数軸上の信号に基づいて、その入力信号における背景騒音を入力信号に追従しながら推定し、推定雑音信号を形成する方法(たとえば、特許第3269969号)を利用すればよい。即ち、音声を含まない雑音(背景雑音)の信号スペクトル信号Xf,kと1分析フレーム前の推定雑音信号N' f,k−1を用いて、現分析フレームの推定雑音信号N' f,kを算出する方法を利用すればよい。この場合、Xf,kとN' f,k−1に係数を掛けてそれぞれをa・N' f,k−1及び(1−a)・Xf,kとし、a・N' f,k−1と(1−a)・Xf,kを加算して、現分析フレームの推定雑音信号N' f,kが算出される。音声判定の方法は、既述の入力信号として環境騒音が重畳した音声信号が時系列順に与えられるフレーム毎に、当該フレームに音声が含まれているか否かを判定する方法(たとえば、特許第3849116号)を利用すればよい。即ち、音声信号の特徴量と、音声周波数の特徴量の高域成分を強調させた特徴量を算出して、音声信号の特徴量に加えて高域成分を強調させた特徴量を基に判定を行う方法を利用すれば良い。また、音声信号の特徴量と、音声の共振周期の規則性を算出して、音声信号の特徴量に加えて音声の共振周期の規則性を基に判定を行う方法を利用すれば良い。   For the processing of the noise estimation unit 4-3, a noise estimation method or a voice determination method can be used. The noise estimation method converts the input signal on the time axis described above into a signal on the frequency axis, and follows the background noise in the input signal based on the converted signal on the frequency axis while following the input signal. A method of estimating and forming an estimated noise signal (for example, Japanese Patent No. 3269969) may be used. That is, using the signal spectrum signal Xf, k of noise (background noise) that does not include speech and the estimated noise signal N ′ f, k−1 one analysis frame before, the estimated noise signal N ′ f, k of the current analysis frame is used. A method for calculating the value may be used. In this case, Xf, k and N′f, k−1 are multiplied by a coefficient to give a · N′f, k−1 and (1-a) · Xf, k, respectively, and a · N′f, k− 1 and (1-a) · Xf, k are added to calculate the estimated noise signal N′f, k of the current analysis frame. The voice determination method is a method for determining whether or not a voice is included in the frame for each frame in which a voice signal on which environmental noise is superimposed as an input signal described above is given in chronological order (for example, Japanese Patent No. 3849116). No.) can be used. In other words, the feature quantity that emphasizes the high frequency component of the audio signal feature quantity and the audio frequency feature quantity is calculated, and the determination is made based on the feature quantity that emphasizes the high frequency component in addition to the feature quantity of the audio signal. You can use the method of doing. Further, a method may be used in which the feature amount of the audio signal and the regularity of the resonance period of the audio are calculated and the determination is made based on the regularity of the resonance cycle of the audio in addition to the feature amount of the audio signal.

事前情報保持手段16には、新幹線など、経路やタイムテーブルが定まっている場合に、予め騒音に関する事前情報を取得し、保持する。   The prior information holding means 16 acquires and holds prior information related to noise in advance when a route or a time table is determined for a bullet train or the like.

現存特定騒音情報算出手段18は日時、位置、エンジン回転数の現在情報が与えられた際に、予め定められた特定騒音リストについて、現存するかどうかの判定情報である現存特定騒音情報を算出する。   The existing specific noise information calculation means 18 calculates the existing specific noise information which is determination information as to whether or not the predetermined specific noise list exists when the current information of the date, position, and engine speed is given. .

寄与係数算出手段22−3は、現存特定騒音情報と音声判定結果を用いて、第1の推定騒音周波数特性に対する第2の推定騒音周波数特性の寄与の大きさを表わす寄与係数を算出する。   The contribution coefficient calculation means 22-3 calculates a contribution coefficient representing the magnitude of the contribution of the second estimated noise frequency characteristic to the first estimated noise frequency characteristic using the existing specific noise information and the voice determination result.

第2の騒音推定手段8は、事前情報12と現存特定騒音情報を用いて騒音の周波数特性(第2の推定騒音周波数特性)を算出する。   The second noise estimation means 8 calculates noise frequency characteristics (second estimated noise frequency characteristics) using the prior information 12 and the existing specific noise information.

第3の騒音推定手段20は、第3の実施の形態と同様に、第1の推定騒音周波数特性と第2の推定騒音周波数特性を、寄与係数で重み付け加算することで推定騒音周波数特性を算出する。つまり、騒音推定手段20には第1の推定騒音周波数特性と第2の推定騒音周波数特性とに応じた推定騒音周波数特性出力が得られる。   The third noise estimation means 20 calculates the estimated noise frequency characteristic by weighting and adding the first estimated noise frequency characteristic and the second estimated noise frequency characteristic with the contribution coefficient, as in the third embodiment. To do. That is, the noise estimation means 20 can obtain an estimated noise frequency characteristic output corresponding to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic.

寄与係数算出手段22−3の寄与係数では、第1の騒音推定手段4−3で算出された音声判定結果を用いて、寄与係数を算出する。この場合、音声判定の結果は次の通りである。
coeff2=0 :音声あり判定
coeff2=1 :音声なし判定
この音声判定の結果により、寄与係数は次の通りである。
coeff2:第2の騒音の周波数特性の寄与係数(0〜1の数)
In the contribution coefficient of the contribution coefficient calculation unit 22-3, the contribution coefficient is calculated using the speech determination result calculated by the first noise estimation unit 4-3. In this case, the result of the voice determination is as follows.
coeff2 = 0: audio presence determination
coeff2 = 1: Determination of no voice According to the result of the voice determination, the contribution coefficient is as follows.
coeff2: contribution factor (number of 0 to 1) of the frequency characteristic of the second noise

つぎに、この騒音推定の処理手順について、図16を参照する。図16は騒音推定の処理手順の一例を示している。   Next, FIG. 16 is referred to regarding the processing procedure of this noise estimation. FIG. 16 shows an example of a processing procedure for noise estimation.

図16に示す処理手順は本開示の騒音推定プログラムの一例である。この処理手順では第1の騒音推定手段4−3は入力音から第1の推定騒音周波数特性を算出する(ステップS71)。第1の騒音推定手段4−3は入力音の音声判定を行う(ステップS72)。   The processing procedure illustrated in FIG. 16 is an example of the noise estimation program of the present disclosure. In this processing procedure, the first noise estimation unit 4-3 calculates a first estimated noise frequency characteristic from the input sound (step S71). The first noise estimation unit 4-3 performs voice determination of the input sound (step S72).

現存特定騒音情報算出手段18は現在情報から現存特定騒音情報を算出する(ステップS73)。寄与係数算出手段22−3は、現存特定騒音情報と音声判定結果とから寄与係数を算出する(ステップS74)。   The existing specific noise information calculation means 18 calculates existing specific noise information from the current information (step S73). The contribution coefficient calculation unit 22-3 calculates a contribution coefficient from the existing specific noise information and the voice determination result (step S74).

第2の騒音推定手段8は現存特定騒音情報と事前情報12とから第2の推定騒音周波数特性を算出する(ステップS75)。   The second noise estimation means 8 calculates a second estimated noise frequency characteristic from the existing specific noise information and the prior information 12 (step S75).

第3の騒音推定手段20は第1の推定騒音周波数特性と第2の推定騒音周波数特性と寄与係数とから第3の推定騒音周波数特性を算出し、出力する(ステップS76)。   The third noise estimation means 20 calculates and outputs a third estimated noise frequency characteristic from the first estimated noise frequency characteristic, the second estimated noise frequency characteristic, and the contribution coefficient (step S76).

この騒音推定装置2−5によれば、寄与係数算出手段22−3が第1の騒音推定手段4−1で算出された音声判定結果を用いて、第2の推定騒音周波数特性の寄与度を操作する。このため、音声が含まれる可能性が高い場合には、第2の推定騒音周波数特性の寄与度を大きくすることで、より騒音の推定性能を高めることができる。   According to this noise estimation apparatus 2-5, the contribution coefficient calculation means 22-3 uses the speech determination result calculated by the first noise estimation means 4-1, and calculates the contribution degree of the second estimated noise frequency characteristic. Manipulate. For this reason, when there is a high possibility that speech is included, the noise estimation performance can be further enhanced by increasing the contribution of the second estimated noise frequency characteristic.

〔上記実施の形態のハードウェア〕 [Hardware of the above embodiment]

上記実施の形態に用いられるハードウェアについて、図17を参照する。既述の騒音推定装置2−1、2−2、2−3、2−4、2−5はコンピュータシステム200で構成することができる。図17はコンピュータシステム200の構成例を示している。   FIG. 17 is referred to for hardware used in the above embodiment. The noise estimation apparatuses 2-1, 2-2, 2-3, 2-4, and 2-5 described above can be configured by the computer system 200. FIG. 17 shows a configuration example of the computer system 200.

図17に示すコンピュータシステム200は、CPU(Central Processing Unit :中央処理装置)202、ROM(Read-Only Memory)204、RAM(Random-Access Memory)206、通信インターフェース208、記憶装置210、入出力装置212、読取り装置214、音声入力部215およびバス216によって構成されている。   A computer system 200 shown in FIG. 17 includes a CPU (Central Processing Unit) 202, a ROM (Read-Only Memory) 204, a RAM (Random-Access Memory) 206, a communication interface 208, a storage device 210, and an input / output device. 212, a reader 214, a voice input unit 215, and a bus 216.

ROM204または記憶装置210はプログラム記憶手段の一例である。記憶装置210には、ハードディスク、磁気ディスクなどの記憶装置を用いることができる。この記憶装置210またはROM204には既述のフローチャートで示したプログラムや既述のテーブル21、22、23のデータが格納される。RAM206はワークエリアを構成する。   The ROM 204 or the storage device 210 is an example of a program storage unit. A storage device such as a hard disk or a magnetic disk can be used as the storage device 210. The storage device 210 or the ROM 204 stores the program shown in the above-described flowchart and the data of the tables 21, 22, and 23 described above. The RAM 206 constitutes a work area.

CPU202は、既述のプログラムを実行し、騒音推定、特定騒音情報の算出、事前情報の取得など、既述のフローチャートに記載した処理を実行する。   The CPU 202 executes the above-described program, and executes the processes described in the above-described flowchart, such as noise estimation, calculation of specific noise information, and acquisition of prior information.

音声入力部215は、既述の入力音6、現在情報10、または事前情報12(図1)の入力に用いられ、具体的には、音声入力や騒音入力を電気信号に変換するマイクロフォン217などを備えている。また、音声入力には、ディジタル化された音声を入力するようにしてもよい。   The voice input unit 215 is used to input the input sound 6, the current information 10, or the prior information 12 (FIG. 1) described above, and specifically, a microphone 217 that converts voice input or noise input into an electrical signal, or the like. It has. In addition, digitized voice may be input for voice input.

通信インターフェース208はネットワーク218を介して情報提供手段220と接続され、既述した各種情報の提供を受けることができる。また、この情報は、可搬型記録媒体222に格納された情報を読取り装置214で読み取ってもよい。また、可搬型記録媒体222に既述のプログラムを格納して流通させることも可能である。   The communication interface 208 is connected to the information providing unit 220 via the network 218, and can receive various information described above. As this information, information stored in the portable recording medium 222 may be read by the reading device 214. It is also possible to store the above-described program in the portable recording medium 222 for distribution.

可搬型記憶媒体222はCD−ROM、フレシキブルディスク、光ディスク、光磁気ディスクなど様々な形式の記憶媒体で構成される。したがって、記憶媒体222に格納されたプログラムを読取り装置214で読み取り、既述の騒音推定を行うことができる。   The portable storage medium 222 includes various types of storage media such as a CD-ROM, a flexible disk, an optical disk, and a magneto-optical disk. Therefore, the program stored in the storage medium 222 can be read by the reading device 214 to perform the noise estimation described above.

〔他の実施の形態〕 [Other Embodiments]

(1) 第4の実施の形態の騒音推定装置2−3において、第3の実施の形態で述べた寄与係数算出手段22−1を適用してもよい。   (1) In the noise estimation apparatus 2-3 according to the fourth embodiment, the contribution coefficient calculation unit 22-1 described in the third embodiment may be applied.

(2) 第5の実施の形態の騒音推定装置2−4または第6の実施の形態の騒音推定装置2−5において、第4の実施の形態で述べた事前情報更新手段24を適用してもよい。   (2) In the noise estimation device 2-4 of the fifth embodiment or the noise estimation device 2-5 of the sixth embodiment, the prior information updating means 24 described in the fourth embodiment is applied. Also good.

(3) 上記実施の形態の第3の騒音推定手段20の第1または第2の周波数特性の加算処理において、重み付けはいずれか一方または双方に付与する処理であってもよい。   (3) In the addition process of the first or second frequency characteristic of the third noise estimation unit 20 of the above embodiment, the weighting may be a process given to either one or both.

以上説明したように、騒音推定装置及び騒音推定プログラムの最も好ましい実施の形態等について説明したが、本発明は、上記記載に限定されるものではなく、特許請求の範囲に記載され、又は発明を実施するための形態に開示された発明の要旨に基づき、当業者において様々な変形や変更が可能であることは勿論であり、斯かる変形や変更が、本発明の範囲に含まれることは言うまでもない。
As described above, the most preferred embodiment of the noise estimation device and the noise estimation program has been described. However, the present invention is not limited to the above description, and is described in the scope of claims. It goes without saying that various modifications and changes can be made by those skilled in the art based on the gist of the invention disclosed in the embodiment, and such modifications and changes are included in the scope of the present invention. Yes.

2 騒音推定装置
2−1 騒音推定装置
2−2 騒音推定装置
2−3 騒音推定装置
2−4 騒音推定装置
2−5 騒音推定装置
4 第1の騒音推定手段
4−1 第1の騒音推定手段
4−2 第1の騒音推定手段
6 入力音
8 第2の騒音推定手段
10 現在情報
12 事前情報
14 加算手段
16 事前情報保持手段
18 現存特定騒音情報算出手段
20 第3の騒音推定手段
21 特定騒音情報テーブル
22 特定騒音情報テーブル
22−1 寄与係数算出手段
22−2 寄与係数算出手段
22−3 寄与係数算出手段
23 特定騒音種テーブル
24 事前情報更新手段
26 信頼度算出手段
DESCRIPTION OF SYMBOLS 2 Noise estimation apparatus 2-1 Noise estimation apparatus 2-2 Noise estimation apparatus 2-3 Noise estimation apparatus 2-4 Noise estimation apparatus 2-5 Noise estimation apparatus 4 1st noise estimation means 4-1 1st noise estimation means 4-2 First noise estimation means 6 Input sound 8 Second noise estimation means 10 Current information 12 Prior information 14 Addition means 16 Prior information holding means 18 Existing specific noise information calculation means 20 Third noise estimation means 21 Specific noise Information table 22 Specific noise information table 22-1 Contribution coefficient calculation means 22-2 Contribution coefficient calculation means 22-3 Contribution coefficient calculation means 23 Specific noise type table 24 Prior information update means 26 Reliability calculation means

Claims (7)

入力音に含まれる騒音の騒音周波数特性を示す第1の推定騒音周波数特性を出力する第1の騒音推定手段と、
前記入力音に関する騒音の事前情報と、前記騒音を特定する現在情報とを用いて特定された騒音の第2の推定騒音周波数特性を出力する第2の騒音推定手段と、
前記第1の推定騒音周波数特性と前記第2の推定騒音周波数特性とに応じた推定騒音周波数特性を出力する加算手段と、
を備えることを特徴とする騒音推定装置。
First noise estimation means for outputting a first estimated noise frequency characteristic indicating a noise frequency characteristic of noise included in the input sound;
Second noise estimation means for outputting a second estimated noise frequency characteristic of the noise identified using the prior information of the noise related to the input sound and the current information identifying the noise;
Adding means for outputting an estimated noise frequency characteristic according to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic;
A noise estimation device comprising:
さらに、前記現在情報から騒音情報を算出し、この騒音情報を現存特定騒音情報として出力する現存特定騒音情報算出手段と、
前記事前情報として特定騒音種の周波数特性を保持する事前情報保持手段と、
を備えることを特徴とする請求項1に記載の騒音推定装置。
Further, noise information is calculated from the current information, and the existing specific noise information calculating means for outputting the noise information as existing specific noise information,
Prior information holding means for holding frequency characteristics of a specific noise type as the prior information;
The noise estimation apparatus according to claim 1, further comprising:
さらに、前記現存特定騒音情報を用いて前記第1の推定騒音周波数特性に対する前記第2の推定騒音周波数特性の寄与の大きさを表わす寄与係数を算出する寄与係数算出手段と、
を備え、前記加算手段は、前記寄与係数を反映させた前記第2の推定騒音周波数特性と、前記第1の推定騒音周波数特性とを加算することを特徴とする請求項2に記載の騒音推定装置。
Furthermore, contribution coefficient calculation means for calculating a contribution coefficient representing the magnitude of the contribution of the second estimated noise frequency characteristic to the first estimated noise frequency characteristic using the existing specific noise information;
The noise estimation unit according to claim 2, wherein the adding means adds the second estimated noise frequency characteristic reflecting the contribution coefficient and the first estimated noise frequency characteristic. apparatus.
さらに、前記現存特定騒音情報と、現在の周波数特性と、音声判定結果とを用いて算出した更新情報により前記事前情報を更新する事前情報更新手段と、
を備え、前記第1の騒音推定手段が、前記入力音から前記第1の推定騒音周波数特性と前記現在の周波数特性と、前記音声判定により前記音声判定結果を算出することを特徴とする請求項2又は3に記載の騒音推定装置。
Further, prior information update means for updating the prior information with update information calculated using the existing specific noise information, the current frequency characteristics, and the voice determination result;
The first noise estimation unit calculates the voice determination result from the input sound based on the first estimated noise frequency characteristic, the current frequency characteristic, and the voice determination. The noise estimation apparatus according to 2 or 3.
さらに、前記現存特定騒音情報と、現在の周波数特性と、音声判定結果から前記事前情報の信頼度情報を算出する信頼度算出手段と、
前記現存特定騒音情報と、前記信頼度情報を用いて前記第1の推定騒音周波数特性に対する前記第2の推定騒音周波数特性の寄与の大きさを表わす寄与係数を算出する寄与係数算出手段と、
を備え、前記第1の騒音推定手段は、前記入力音から前記第1の推定騒音周波数特性と前記現在の周波数特性と、前記音声判定結果を算出し、
前記加算手段は、前記寄与係数を反映させた前記第2の推定騒音周波数特性と、前記第1の推定騒音周波数特性に加算する、
ことを特徴とする請求項2ないし4のいずれかに記載の騒音推定装置。
Further, reliability calculation means for calculating reliability information of the prior information from the existing specific noise information, current frequency characteristics, and voice determination results;
Contribution coefficient calculation means for calculating a contribution coefficient representing the magnitude of the contribution of the second estimated noise frequency characteristic to the first estimated noise frequency characteristic using the existing specific noise information and the reliability information;
The first noise estimation means calculates the first estimated noise frequency characteristic, the current frequency characteristic, and the voice determination result from the input sound,
The adding means adds to the second estimated noise frequency characteristic reflecting the contribution coefficient and the first estimated noise frequency characteristic.
The noise estimation apparatus according to any one of claims 2 to 4, wherein
さらに、前記現存特定騒音情報と、音声判定結果から前記第1の推定騒音周波数特性に対する前記第2の推定騒音周波数特性の寄与の大きさを表わす寄与係数を出力する寄与係数算出手段と、
を備え、前記第1の騒音推定手段は前記入力音から前記第1の推定騒音周波数特性と、前記音声判定結果を出力し、
前記加算手段は、前記寄与係数算出手段で算出した前記寄与係数を反映させた前記第2の推定騒音周波数特性と、前記第1の推定騒音周波数特性とを加算する、
ことを特徴とする請求項2ないし5のいずれかに記載の騒音推定装置。
And a contribution coefficient calculation means for outputting a contribution coefficient representing the magnitude of the contribution of the second estimated noise frequency characteristic to the first estimated noise frequency characteristic from the existing specific noise information and the voice determination result;
The first noise estimation means outputs the first estimated noise frequency characteristic and the voice determination result from the input sound,
The adding means adds the second estimated noise frequency characteristic reflecting the contribution coefficient calculated by the contribution coefficient calculating means and the first estimated noise frequency characteristic;
The noise estimation apparatus according to any one of claims 2 to 5, wherein
コンピュータに実行させる騒音推定プログラムであって、
第1の騒音推定手段入力音に含まれる騒音の騒音周波数特性を示す第1の推定騒音周波数特性を出力し、
第2の騒音推定手段が前記入力音に関する騒音の事前情報と、前記騒音を特定する現在情報とを用いて特定された騒音の第2の推定騒音周波数特性を出力し、
加算手段が前記第1の推定騒音周波数特性と前記第2の推定騒音周波数特性とに応じた推定騒音周波数特性を出力する
ことを前記コンピュータが実行する騒音推定プログラム。
A noise estimation program to be executed by a computer,
Outputting a first estimated noise frequency characteristic indicating a noise frequency characteristic of the noise included in the input sound of the first noise estimating means;
A second noise estimation means for outputting a second estimated noise frequency characteristic of the noise identified using the prior information of the noise related to the input sound and the current information identifying the noise;
A noise estimation program executed by the computer, wherein the adding means outputs an estimated noise frequency characteristic corresponding to the first estimated noise frequency characteristic and the second estimated noise frequency characteristic.
JP2011081276A 2011-03-31 2011-03-31 Noise estimation device and noise estimation program Expired - Fee Related JP5768450B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011081276A JP5768450B2 (en) 2011-03-31 2011-03-31 Noise estimation device and noise estimation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011081276A JP5768450B2 (en) 2011-03-31 2011-03-31 Noise estimation device and noise estimation program

Publications (2)

Publication Number Publication Date
JP2012215715A true JP2012215715A (en) 2012-11-08
JP5768450B2 JP5768450B2 (en) 2015-08-26

Family

ID=47268533

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011081276A Expired - Fee Related JP5768450B2 (en) 2011-03-31 2011-03-31 Noise estimation device and noise estimation program

Country Status (1)

Country Link
JP (1) JP5768450B2 (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0573088A (en) * 1991-09-13 1993-03-26 Toshiba Corp Method for generating recognition dictionary, and recognition dictionary generation device and speech recognition device
JP2000321080A (en) * 1999-05-13 2000-11-24 Denso Corp Noise suppressor, voice recognizer and car navigation system
JP2001075595A (en) * 1999-09-02 2001-03-23 Honda Motor Co Ltd On-vehicle voice recognition device
JP2002314637A (en) * 2001-04-09 2002-10-25 Denso Corp Device for reducing noise
JP2007065122A (en) * 2005-08-30 2007-03-15 Aisin Seiki Co Ltd Noise suppressing device of on-vehicle voice recognition device
JP2009031425A (en) * 2007-07-25 2009-02-12 Nec Corp Noise estimation device, method and program
WO2009038013A1 (en) * 2007-09-21 2009-03-26 Nec Corporation Noise removal system, noise removal method, and noise removal program
JP2009231928A (en) * 2008-03-19 2009-10-08 Alpine Electronics Inc Sound improving device and sound improving method
JP2010102204A (en) * 2008-10-24 2010-05-06 Yamaha Corp Noise suppressing device and noise suppressing method
WO2011148860A1 (en) * 2010-05-24 2011-12-01 日本電気株式会社 Signal processing method, information processing device, and signal processing program

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0573088A (en) * 1991-09-13 1993-03-26 Toshiba Corp Method for generating recognition dictionary, and recognition dictionary generation device and speech recognition device
JP2000321080A (en) * 1999-05-13 2000-11-24 Denso Corp Noise suppressor, voice recognizer and car navigation system
JP2001075595A (en) * 1999-09-02 2001-03-23 Honda Motor Co Ltd On-vehicle voice recognition device
JP2002314637A (en) * 2001-04-09 2002-10-25 Denso Corp Device for reducing noise
JP2007065122A (en) * 2005-08-30 2007-03-15 Aisin Seiki Co Ltd Noise suppressing device of on-vehicle voice recognition device
JP2009031425A (en) * 2007-07-25 2009-02-12 Nec Corp Noise estimation device, method and program
WO2009038013A1 (en) * 2007-09-21 2009-03-26 Nec Corporation Noise removal system, noise removal method, and noise removal program
JP2009231928A (en) * 2008-03-19 2009-10-08 Alpine Electronics Inc Sound improving device and sound improving method
JP2010102204A (en) * 2008-10-24 2010-05-06 Yamaha Corp Noise suppressing device and noise suppressing method
WO2011148860A1 (en) * 2010-05-24 2011-12-01 日本電気株式会社 Signal processing method, information processing device, and signal processing program

Also Published As

Publication number Publication date
JP5768450B2 (en) 2015-08-26

Similar Documents

Publication Publication Date Title
US9767790B2 (en) Voice retrieval apparatus, voice retrieval method, and non-transitory recording medium
JP4316583B2 (en) Feature amount correction apparatus, feature amount correction method, and feature amount correction program
JP4202090B2 (en) Speech synthesis system using smoothing filter and method thereof, smoothing filter characteristic control device and method thereof
JP5229234B2 (en) Non-speech segment detection method and non-speech segment detection apparatus
JP6035702B2 (en) Sound processing apparatus and sound processing method
JP2012510666A5 (en)
JP4769673B2 (en) Audio signal interpolation method and audio signal interpolation apparatus
US8630857B2 (en) Speech synthesizing apparatus, method, and program
US20110071823A1 (en) Speech recognition system, speech recognition method, and storage medium storing program for speech recognition
KR102188620B1 (en) Sinusoidal interpolation across missing data
JP4735398B2 (en) Acoustic signal analysis apparatus, acoustic signal analysis method, and acoustic signal analysis program
JP5768450B2 (en) Noise estimation device and noise estimation program
JP5815435B2 (en) Sound source position determination apparatus, sound source position determination method, program
JP6216809B2 (en) Parameter adjustment system, parameter adjustment method, program
JP4486527B2 (en) Acoustic signal analyzing apparatus and method, program, and recording medium
JP5191500B2 (en) Noise suppression filter calculation method, apparatus, and program
JP2004070353A (en) Device and method for inter-signal correlation coefficient determination, and device and method for pitch determination using same
JP3650102B2 (en) Software output generating apparatus and method
JP2009063700A (en) Device, method and program for estimating voice signal section, and storage medium recording the program
KR101756287B1 (en) Apparatus and method for extracting features for speech recognition
JP2016133522A (en) Glottis closing time estimation device, pitch mark time estimation device, pitch waveform connection point estimation device, and method and program thereof
JP6152690B2 (en) Acoustic analyzer
JP2003167599A (en) Speech recognition method and speech recognition device
JP4653673B2 (en) Signal determination apparatus, signal determination method, signal determination program, and recording medium
US20220051657A1 (en) Channel selection apparatus, channel selection method, and program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140204

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20141110

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20141118

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150116

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20150608

R150 Certificate of patent or registration of utility model

Ref document number: 5768450

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees