(57)【要約】  機器から発せられる異常音のうち、振幅強度（音圧）の変化が生じない、もしくは僅少な場合でかつ、位相変化が生じている異常音を検知することを目的とする。  診断対象の機器から採取した音データを入力信号として前記機器の異常診断を行う異常音検出システムにおいて、前記入力信号の１周期分の音データの零点の大きさを計算する零点大きさ計算部と、前記零点の大きさから零点密度分布を計算する零点密度分布計算部と、前記入力信号の零点密度分布と、正常な音データの零点密度分布とを比較して前記診断対象の機器が正常か異常かを判定する正常／異常診断部とを有する。
2014016914.tif 000002
【特許請求の範囲】
【請求項１】
診断対象の機器から採取した音データを入力信号として前記機器の異常診断を行う異常音検出システムにおいて、
前記入力信号の１周期分の音データの零点の大きさを計算する零点大きさ計算部と、
前記零点の大きさから零点密度分布を計算する零点密度分布計算部と、
前記入力信号の零点密度分布と、正常な音データの零点密度分布とを比較して前記診断対象の機器が正常か異常かを判定する正常／異常診断部とを有することを特徴とする異常音検出システム。
【請求項２】
前記正常／異常診断部は、前記入力信号の零点密度分布と前記正常な音データの零点密度分布との相関係数を計算し、前記相関係数と正常異常判定閾値とを比較することにより前記診断対象の機器が正常か異常かを判定することを特徴とする請求項１に記載の異常音検出システム。
【請求項３】
前記診断対象の機器が正常か異常かの判定結果を表示させる表示制御部を有することを特徴とする請求項１または２に記載の異常音検出システム。
【請求項４】
前記表示制御部は、前記入力信号の零点密度分布と、前記正常な音データの零点密度分布とを表示させることを特徴とする請求項３に記載の異常音検出システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、異常音検出システムに関する。
【背景技術】
【０００２】
本技術分野の背景技術として、機器の稼働率維持のため、機器の診断が重要となっている。機器の診断時、診断対象部位にセンサの取り付けが困難な場合の対応として、診断対象部位に接触が不要である音響による診断が注目されている。そのため、従来より、音響を用いた診断が用いられている。【０００３】
特許文献１には、「異常が発生した場合に所定の周期で異音を生じる音源の診断を支援する音響診断支援装置」であって、「異常を容易に把握することができるようにする」と記載されている。そのために、「集音器は、回転機から発生した音を所定周波数で回転周期よりも長い所定時間サンプリングした測定データを記録する。診断支援装置は、集音器から測定データを取得する。診断支援装置は、１番目の測定データにおける先頭から１周期分の基準データ系列を抽出し、ｎ番目の測定データにおける抽出位置をずらしながら、１周期分の比較データ系列を抽出し、基準データ系列と比較データ系列の相関度を算出していき、相関度の最も大きい抽出位置をシフト量として、２番目の測定データをシフトした上で、１番目の測定データに足し合わせる」と記載されている。【０００４】
また、特許文献２には、「診断対象から発生する音響信号に基づいて、前記診断対象で発生する異常を診断する方法であって、前記診断対象が正常な状態にあるときに、前記音響信号を周波数分析して得た周波数または周波数帯域と音圧との関係に基づいて、各前記周波数または各前記周波数帯域の音圧とその前記周波数または周波数帯域の周辺の他の１以上の前記周波数または周波数帯域に対応する音圧との相対的な差を相対的な音圧差として求める異常診断方法」と記載されている。【０００５】
さらに、特許文献３には、「蓄積された信号系列の中から所定の信号、またはその一部に類似した信号を探索して検出する信号検出方法に関するものであり、例えば音響信号検出に適用可能である。従来、信号検出方法に関しては、蓄積信号中で目的信号に類似した箇所を検出することを目的とした信号検索方法が知られている。しかし、局所的な枝刈りのみを用いていたため、膨大な蓄積信号を対象とする場合には、検索に長時間を要するという欠点があった。Ｌ１距離に基づいて、大局的なグループ化および局所的なグループ化を行い、探索空間を効率的に絞り込むことによって、探索制度を保ちつつ、高速に効果的な部分信号検出ができるという利点がある」と記載されている。【０００６】
このように、機器から発せられる異常音検出方法および特徴量、類似情報の検索については、従前より取り組まれている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１２−９３０９４号公報
【特許文献２】特開２００５−２５７４６０号公報
【特許文献３】国際公開第２００６／００９０３５号パンフレット
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、従来手法では、異常音の変化が振幅強度（音圧）に現れることが前提となっている。そのため、振幅強度（音圧）に異常が生じなく、位相変化として出現する場合、たとえば、空気吸入口つまりといった振幅変化が小さい異常に関しては、検知を行うことができない。【０００９】
特許文献１は、異常が生じた場合、音響データの周期データを複数取得し、同期をとったうえでそれらを重ね合わせることで周期的に発生している音源の異常を把握する発明である。このため、周期的に発生していない音源の異常、例えば位相変化が生じる異常音については把握することができない。そのため、振幅変化の小さい、空気吸入口つまりといった異常、初期の異常については検出することができない。【００１０】
特許文献２は、音響信号を周波数分析して得ることができる周波数または周波数帯域と音圧との関係に基づいて、異常判定用レベル差上限閾値と異常判定用レベル差下限閾値とを定め、診断時においては、周波数成分もしくはその音圧を用いて診断を行う発明である。そのため、異常音が周波数帯域ならびに音圧値（大きさ）に出現しない異常については検出することができない。【００１１】
特許文献３では、Ｌ１距離を用いて高速な探索を行うために特徴量をヒストグラム化して検索を行う発明である。また、映像の色を特徴量としている。Ｌ１距離とは、距離の差が１乗に基づく距離と定義されている。【００１２】
このように、周波数帯域、振幅強度（音圧）に変化がなく、位相変化が生じている異常音を検知する手法については、先行技術文献には記載されていない。【００１３】
本発明は、機器から発せられる異常音のうち、振幅強度（音圧）の変化が生じない、もしくは僅少な場合でかつ、位相変化が生じている異常音を検知することを目的とする。
【課題を解決するための手段】
【００１４】
そこで、本発明では、異常診断に用いる特徴量として、零点に着目した診断を行う。零点とは、時間領域の波形での零交差であり、周波数領域上においてはエネルギーが零となる点のことである。零点の大きさの密度分布は、位相変化に対応した分布となる。これは、時間領域での零交差の間隔が零点の大きさとして、複素数平面上に出現し、複素数平面上においては、大きさ毎に零点の数を数えることが可能となる。すなわち、位相変化が生じない場合においては、分布が一意に定まり、位相変化が生じる場合においては、零点の密度分布に変化が生じる。【００１５】
以上のことから、零点密度分布を用いることにより、従来手法で検知することが不可能である、振幅強度（音圧）変化が生じない異常音のうち、位相変化を伴う異常音の検知を行う。【００１６】
本発明は、例えば、診断対象の機器から採取した音データを入力信号として前記機器の異常診断を行う異常音検出システムにおいて、前記入力信号の１周期分の音データの零点の大きさを計算する零点大きさ計算部と、前記零点の大きさから零点密度分布を計算する零点密度分布計算部と、前記入力信号の零点密度分布と、正常な音データの零点密度分布とを比較して前記診断対象の機器が正常か異常かを判定する正常／異常診断部とを有することを特徴とする。
【発明の効果】
【００１７】
本発明によれば、振幅強度の変化が小さく位相変化が生じている異常音の検知が可能となる。【００１８】
上記した以外の課題、構成および効果は、以下の実施形態の説明により明らかにされる。
【図面の簡単な説明】
【００１９】
【図１】異常音検出システムの構成図の例である。
【図２】ハードウェア構成の例である。
【図３】システムフローチャートの例である。
【図４】１周期計算部の処理の例である。
【図５】零点密度分布データ構成の例である。
【図６】零点密度分布の例である。
【図７】零点密度分布表示の例である。
【図８】モーター音の診断例である。
【発明を実施するための形態】
【００２０】
本発明の実施例を、図面を参照しながら説明する。尚、各図において、同一又は類似の構成要素には同じ符号を付し、説明を省略する。
【００２１】
図１は、異常音検出システムの構成図の例である。図２は、ハードウェア構成の例である。図１に示すように、異常音検出システム１は、１周期計算部１０１、零点計算部１０２、零点大きさ計算部１０３、零点密度分布計算部１０４、蓄積信号照合部１０５、正常／異常診断部１０６、零点密度データベース１０７を有している。図１に示す異常音検出システム１は、図２に示すデータ処理部Ｈ０２に対応しており、例えばコンピュータやマイコンなどにより処理が実現される。また、図１に示す密度分布などの結果を表示するディスプレイ１０８は、図２における出力結果表示部Ｈ０３に対応する。
【００２２】
異常音検出システム１は、入力信号Ｉ１を１周期計算部１０１への入力とし、１周期計算部１０１は、入力された入力信号Ｉ１のうち１周期分を切り出し、入力信号１周期分の信号として零点計算部１０２へと出力する。入力信号Ｉ１は、診断対象の機器から採取された音データであり、入力信号Ｉ１の入力手段については、図２に示すマイクＨ０１からデータ処理部Ｈ０２へ直接入力する場合と、マイクＨ０１からデータロガー等を介してデータ処理部Ｈ０２へ入力されるケースが考えられる。
【００２３】
零点計算部１０２は、１周期計算部１０１により計算された入力信号１周期分の信号から零点の算出を行い、算出された零点を零点大きさ計算部１０３へと出力する。零点の算出方法については、ｎ次多項式近似を行い、ｎ次の近似式を解くことにより算出される。多項式近似方法としては、べき級数展開による方法、ラグランジュ補間による手法等、数値計算（参考文献：数値計算法、長嶋秀世著）の手法が公知となっている。
【００２４】
零点大きさ計算部１０３では、零点計算部１０２で算出された零点から、零点の大きさを計算する。零点は、複素数で定義されている。そのため、複素数の絶対値をとることで大きさを算出する。零点大きさ計算部１０３で計算された零点の大きさは、零点密度分布計算部１０４へと入力する。
【００２５】
零点密度分布計算部１０４は、入力された零点の大きさのデータから、零点密度分布を計算する。具体的には、零点の全個数における特定の零点の大きさを持つ零点の個数の割合を計算する。零点は、近似した次数の数の分だけ算出される。例えば、６次関数に近似した場合、零点は最大６個算出される。また、６周期分の波形が存在しており、６次関数に近似した場合は、零点の合計個数は、３６個である。零点密度分布計算部１０４で計算された入力信号Ｉ１の零点密度分布は、蓄積信号照合部１０５に出力される。
【００２６】
また、異常音検出システム１は、蓄積信号Ｉ２でも、１周期計算部１０１、零点計算部１０２、零点大きさ計算部１０３、零点密度分布計算部１０４において同様の計算を行う。ここで、蓄積信号Ｉ２は、機器が正常な状態で採取された正常音である。蓄積信号Ｉ２によって計算された蓄積信号Ｉ２の零点密度分布は、零点密度データベース１０７に記録される。
【００２７】
蓄積信号照合部１０５は、零点密度データベース１０７に構築されている蓄積信号Ｉ２の零点密度分布のデータベースの中から、入力信号Ｉ１により計算された入力信号Ｉ１の零点密度分布と比較するための蓄積信号Ｉ２の零点密度分布を読み出し、入力信号Ｉ１の零点密度分布と蓄積信号Ｉ２の零点密度分布とを、正常／異常診断部１０６へと出力する。
【００２８】
零点密度データベース１０７には、正常データとして蓄積信号Ｉ２の零点密度分布を蓄積している。
【００２９】
正常／異常診断部１０６は、入力信号Ｉ１の零点密度分布と蓄積信号Ｉ２の零点密度分布とを比較し、正常または異常の診断を行い、診断結果を出力する。入力信号Ｉ１の零点密度分布と蓄積信号Ｉ２の零点密度分布との比較方法については、入力信号Ｉ１の零点密度分布と蓄積信号Ｉ２の零点密度分布との相関係数を算出することにより比較する方法等が考えられる。
【００３０】
相関係数を用いて比較する場合、相関係数が１．０である場合は相関が強いことを示しており、蓄積信号Ｉ２と入力信号Ｉ１は同じであることを示しており、相関係数が０．０である場合は相関がないことを示しており、蓄積信号Ｉ２と入力信号Ｉ１は異なっていることを示す。相関係数については、確率変数間の類似度を統計的に評価する手法として既に公知となっている。なお、相関係数を用いて正常と異常とを識別するための閾値については、診断対象の機器、検知を行う機器、故障の進展度合い等により定める必要がある。例えば、正常異常判定閾値を０．８とし、相関係数が正常異常判定閾値以上であれば正常、正常異常判定閾値未満であれば異常と判定し、判定結果をディスプレイ１０８に出力する。
【００３１】
ディスプレイ１０８は、正常／異常診断部１０６で判定された正常または異常の判定結果を表示する。異常音検出システム１は、図示しない表示制御部を有しており、ディスプレイ１０８への表示を制御する。尚、ディスプレイ１０８は、例えば後述する図７で説明するように、その他の情報も一緒に表示しても良い。
【００３２】
図３は、システムフローチャートの例である。１周期計算ステップＦ０１では、１周期計算部１０１に、入力信号Ｉ１もしくは蓄積信号Ｉ２と、診断対象の機器の回転数Ｉ１０、診断対象の機器の機器名Ｉ５０が入力される。回転数Ｉ１０、機器名Ｉ５０の入力方法は、キーボード、オンライン等が考えられ、入力方法は問わない。１周期計算部１０１は、入力信号１周期分（あるいは蓄積信号１周期分）の信号を計算して出力する。
【００３３】
１周期計算ステップＦ０１にて計算された入力信号１周期分の信号（あるいは蓄積信号１周期分）に基づいて、零点計算ステップＦ０２にて、零点計算部１０２で零点が計算されて出力される。零点大きさ計算ステップＦ０３では、零点大きさ計算部１０３で零点の大きさを計算する。
【００３４】
零点密度分布計算ステップＦ０４では、零点密度分布計算部１０４によって零点密度分布の計算がなされる。零点密度分布計算部１０４では、密度分布の間隔Ｉ６０も入力される。間隔Ｉ６０の入力方法は、キーボードからの直接入力、オンライン等を問わない。間隔Ｉ６０とは、密度分布算出時において、零点の大きさ毎の個数を算出するための間隔である。例えば、間隔Ｉ６０が１０であるということは、零点の大きさを１０刻みで計算していくことである。
【００３５】
蓄積信号照合ステップＦ０５では、蓄積信号照合部１０５に、零点密度分布計算ステップＦ０４にて計算された入力信号Ｉ１の零点密度分布、入力信号Ｉ１に対応する診断対象の機器の回転数Ｉ１０、入力信号Ｉ１に対応する診断対象の機器名Ｉ５０が入力され、入力信号Ｉ１に対応する診断対象の機器の回転数Ｉ１０、入力信号Ｉ１に対応する診断対象の機器名Ｉ５０に対応する蓄積信号Ｉ２の零点密度分布が零点密度データベース１０７から読み出される。そして、正常／異常診断ステップＦ０６では、正常／異常診断部１０６により、正常／異常の判定が行われる。
【００３６】
図４は、１周期計算部の処理の例である。１周期計算部１０１においては、入力信号Ｉ１の回転数Ｉ１０、収録時間Ｉ２０、サンプリング周波数Ｉ３０から、入力信号Ｉ１の１周期を計算し、１周期分を切り出して入力信号１周期分Ｉ４０を出力する。例えば、回転数６０（Ｈｚ）、収録時間２０（秒）、サンプリング周波数５００００（Ｈｚ）の場合、回転数Ｉ１０には６０（Ｈｚ）、収録時間Ｉ２０には２０（秒）、サンプリング周波数Ｉ３０には５００００（Ｈｚ）が入力される。入力手段については、キーボード、オンライン等のデータを問わない。尚、図４では入力信号Ｉ１の場合の処理を例に示したが、蓄積信号Ｉ２の場合の処理も同様である。
【００３７】
図５は、零点密度分布データ構成の例である。図５に示す零点密度データベース１０７のデータ構成Ｋ０１は、機器名Ｉ５０、回転数Ｉ１０、零点密度分布Ｉ７０、間隔Ｉ６０により構成される。回転数Ｉ１０については、１分あたりの回転数であるｒｐｍ、もしくは、１秒あたりの回転数であるＨｚのどちらかが考えられる。１分あたりの回転数であるｒｐｍの場合、別途、１秒あたりの回転数であるＨｚへ変換する必要がある。例えば、６００ｒｐｍが入力された場合、１０Ｈｚへと変換する必要がある。図５の例では、機器Ｘ、回転数１０（Ｈｚ）、零点密度分布は零点の大きさの間隔が１０毎にそれぞれ１０％、２０％、４０％、２０％、１０％であることを示している。
【００３８】
図６は、零点密度分布の例である。図６に示している零点密度分布Ｍ０１は、横軸に零点の大きさ、縦軸に零点密度（％）を示している。零点密度（％）は、全零点の個数における、特定の零点の大きさを持つ零点の個数の割合のことである。例えば、図６において、全零点の個数が６００個の場合で、零点の大きさが０以上１０未満の個数が６０個、同様に大きさが１０以上２０未満の個数が１２０個、２０以上３０未満の個数が２４０個、３０以上４０未満の個数が１２０個、４０以上５０未満の個数が６０個の場合、零点の大きさが０以上１０未満の密度は１０（％）、１０以上２０未満の密度が２０（％）、２０以上３０未満の密度が４０（％）、３０以上４０未満の密度が２０（％）、４０以上５０未満の密度が１０（％）であることを示している。
【００３９】
図７は、零点密度分布表示の例である。図７に示している零点密度分布表示Ｄ０１は、ディスプレイ１０８に表示された画面の一例である。入力信号Ｉ１の零点密度分布Ｍ１を点線で、入力信号Ｉ１に対応する蓄積信号Ｉ２の零点密度分布Ｍ２を実線で表示しており、横軸は零点の大きさ、縦軸は零点密度（％）である。また、図７では、正常／異常診断部１０６で計算された指標も表示した例を示している。例えば、相関係数を用いて正常音と異常音とを判定する場合で、計算された相関係数が０．７、正常異常判定閾値が０．８の場合、図７において、数値評価指標は「相関係数０．７」、正常異常判定閾値は「相関係数０．８」、判定結果は「異常」と表示した例を示している。尚、これらの表示の制御は、異常音検出システム１の図示しない表示制御部によって行われる。
【００４０】
図８は、モーター音の診断例である。図８では、診断対象の機器をモーターＪ０１とした場合を示している。モーターＪ０１には、モーター駆動のための電源Ｊ０２が接続されている。
【００４１】
そして、モーターＪ０１付近にマイクＨ０１を設置する。マイクＨ０１から収録された音声（音データ）は、データロガーＪ０３へと入力される。データロガーＪ０３に入力された収録音は、データ処理部Ｈ０２へ入力信号Ｉ１として入力され、診断処理が行われ、診断結果が出力結果表示部Ｈ０３に表示される。
【００４２】
尚、蓄積信号Ｉ２として正常音を予め収録し、データ処理部Ｈ０２で蓄積信号Ｉ２の零点密度分布を計算し、零点密度データベース１０７に零点密度分布データとして収録しておく。
【００４３】
蓄積信号Ｉ２の採取方法としては、例えば６０秒毎に複数回に分けて収録を行う方法や、機器を連続稼動し信号を蓄積していく方法が考えられる。入力信号Ｉ１に関しても同様に、複数回に分けて入力する方法や、連続稼動しながら信号を入力し診断を行う方法が考えられる。
【００４４】
以上のように、異常診断の際に、零点密度分布を用いることにより、従来手法で検知することが不可能である、振幅強度（音圧）変化が生じない異常音のうち、位相変化を伴う異常音の検知を行うことができる。
【００４５】
以上、本発明の実施例を説明してきたが、これまでの実施例で説明した構成はあくまで一例であり、本発明は、技術思想を逸脱しない範囲内で適宜変更が可能である。
【符号の説明】
【００４６】
１        異常音検出システム
１０１    １周期計算部
１０２    零点計算部
１０３    零点大きさ計算部
１０４    零点密度分布計算部
１０５    蓄積信号照合部
１０６    正常／異常診断部
１０７    零点密度データベース
１０８    ディスプレイ
Ｈ０１    マイク
Ｈ０２    データ処理部
Ｈ０３    出力結果表示部
Ｆ０１    １周期計算ステップ
Ｆ０２    零点計算ステップ
Ｆ０３    零点大きさ計算ステップ
Ｆ０４    零点密度分布計算ステップ
Ｆ０５    蓄積信号照合ステップ
Ｆ０６    正常／異常診断ステップ
Ｉ１      入力信号
Ｉ２      蓄積信号
Ｉ１０    回転数
Ｉ２０    収録時間
Ｉ３０    サンプリング周波数
Ｉ４０    入力信号１周期分
Ｉ５０    機器名
Ｉ６０    間隔
Ｉ７０    零点密度分布
Ｋ０１    零点密度分布データ構成
Ｍ０１    零点密度分布
Ｍ１      入力信号の零点密度分布
Ｍ２      蓄積信号の零点密度分布
Ｄ０１    零点密度分布表示
Ｊ０１    モーター
Ｊ０２    電源
Ｊ０３    データロガー
【図１】
2014016914.tif 000003
【図２】
2014016914.tif 000004
【図３】
2014016914.tif 000005
【図４】
2014016914.tif 000006
【図５】
2014016914.tif 000007
【図６】
2014016914.tif 000008
【図７】
2014016914.tif 000009
【図８】
2014016914.tif 000010
【手続補正書】
【提出日】平成27年12月3日(2015.12.3)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
診断対象の機器から採取した音データを入力信号として前記機器の異常診断を行う異常音検出システムにおいて、
前記診断対象の機器の回転数と前記入力信号とが入力され前記入力信号の１周期分の信号を出力する１周期計算部と、
前記１周期分の信号を周波数領域上でｎ次多項式近似を行い、ｎ次の近似式を解くことによりｎ個の零点を複素数で計算する零点計算部と、
前記零点の複素数の絶対値をとることで零点の大きさを計算する零点大きさ計算部と、
前記入力信号の複数周期分の信号から得られた前記零点の大きさから前記零点の大きさの密度分布である零点密度分布を計算する零点密度分布計算部と、
前記入力信号の零点密度分布と、正常な音データの零点密度分布とを比較して前記診断対象の機器が正常か異常かを判定する正常／異常診断部とを有することを特徴とする異常音検出システム。
【請求項２】
前記正常／異常診断部は、前記入力信号の零点密度分布と前記正常な音データの零点密度分布との相関係数を計算し、前記相関係数と正常異常判定閾値とを比較することにより前記診断対象の機器が正常か異常かを判定することを特徴とする請求項１に記載の異常音検出システム。
【請求項３】
前記診断対象の機器が正常か異常かの判定結果を表示させる表示制御部を有することを特徴とする請求項１または２に記載の異常音検出システム。
【請求項４】
前記表示制御部は、前記入力信号の零点密度分布と、前記正常な音データの零点密度分布とを表示させることを特徴とする請求項３に記載の異常音検出システム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１４
【補正方法】変更
【補正の内容】
【００１４】
そこで、本発明では、異常診断に用いる特徴量として、零点に着目した診断を行う。零点とは、時間領域の波形での零交差であり、周波数領域上においてはエネルギーが零となる点のことである。周波数領域上での零点の大きさの密度分布は、位相変化に対応した分布となる。すなわち、位相変化が生じない場合においては、分布が一意に定まり、位相変化が生じる場合においては、零点の密度分布に変化が生じる。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
本発明は、例えば、診断対象の機器から採取した音データを入力信号として前記機器の異常診断を行う異常音検出システムにおいて、前記診断対象の機器の回転数と前記入力信号とが入力され前記入力信号の１周期分の信号を出力する１周期計算部と、前記１周期分の信号を周波数領域上でｎ次多項式近似を行い、ｎ次の近似式を解くことによりｎ個の零点を複素数で計算する零点計算部と、前記零点の複素数の絶対値をとることで零点の大きさを計算する零点大きさ計算部と、前記入力信号の複数周期分の信号から得られた前記零点の大きさから前記零点の大きさの密度分布である零点密度分布を計算する零点密度分布計算部と、前記入力信号の零点密度分布と、正常な音データの零点密度分布とを比較して前記診断対象の機器が正常か異常かを判定する正常／異常診断部とを有することを特徴とする。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００２３
【補正方法】変更
【補正の内容】
【００２３】
零点計算部１０２は、１周期計算部１０１により計算された入力信号１周期分の信号から零点の算出を行い、算出された零点を零点大きさ計算部１０３へと出力する。零点の算出方法については、入力信号１周期分の信号を周波数領域上でｎ次多項式近似を行い、ｎ次の近似式を解くことにより算出される。多項式近似方法としては、べき級数展開による方法、ラグランジュ補間による手法等、数値計算（参考文献：数値計算法、長嶋秀世著）の手法が公知となっている。
【国際調査報告】
2014016914.tif 000011
2014016914.tif 000012
2014016914.tif 000013
2014016914.tif 000014