JP7485061B2

JP7485061B2 - 信号処理装置、信号処理方法、及びプログラム

Info

Publication number: JP7485061B2
Application number: JP2022553244A
Authority: JP
Inventors: 咲子美島; 玲史近藤
Original assignee: NEC Corp
Current assignee: NEC Corp
Priority date: 2020-09-29
Filing date: 2020-09-29
Publication date: 2024-05-16
Anticipated expiration: 2040-09-29
Also published as: JPWO2022070234A1; US20230296409A1; WO2022070234A1

Description

本開示は、信号の変化を検出する技術に関する。

複数種類の信号が混入する時系列信号に対して、信号の変化を検出する技術が存在する。特に、突発的に変化する信号である急変信号の検出に関する技術が、特許文献１及び特許文献２に開示されている。

特許文献１には、入力信号から周波数ごとに位相成分信号の第１の傾きを算出し、算出した第１の傾きの平均値である第２の傾きを算出し、第１の傾きと第２の傾きとの差に基づいて、入力信号における急変の有無を検知する技術が開示されている。

特許文献２には、第１の傾きを振幅に基づいて重み付けし、重み付き第１の傾きの平均を第２の傾きとして算出し、第２の傾きと第１の傾きとの差に基づいて、入力信号における急変の有無を検知する技術が開示されている。

特許６４０６２５８号公報特許６４０６２５７号公報

ここで、信号における急変の有無の検知とは、信号に急変信号が含まれるか否かの検知を示す。入力信号において、急変信号に対して、急変信号でない信号（非急変信号）の割合が大きい場合、すなわち、急変信号対非急変信号比が小さい場合に、急変信号に起因する位相成分信号の傾きが得られる周波数の割合が少なくなる。

特許文献１及び２に記載の技術では、急変の有無を検知したい区間の信号のみから第１の傾きと第２の傾きとを算出し、単に、周波数ごとの第１の傾きと第２の傾きとの差分を算出している。そのため、上記のように、急変信号に起因する位相成分信号の傾きが得られる周波数の割合が少なくなる場合には、急変の有無を検知する精度は悪くなる。

本開示は、上記課題を鑑みてなされたものであり、信号の急変を精度よく検知することができる信号処理装置等を提供することを目的の一つとする。

本開示の一態様にかかる信号処理装置は、入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す切出手段と、前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換する変換手段と、前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する傾き算出手段と、複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出するスコア算出手段と、前記スコアに基づいて、前記判定対象信号の急変の有無を判定する判定手段と、を備える。

本開示の一態様にかかる信号処理方法は、入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出し、前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換し、前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出し、複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出し、前記スコアに基づいて、前記判定対象信号の急変の有無を判定する。

本開示の一態様にかかるコンピュータ読み取り可能な記憶媒体は、入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す処理と前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換する処理と前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する処理と、複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出する処理と、前記スコアに基づいて、前記判定対象信号の急変の有無を判定する処理と、をコンピュータに実行させるプログラムを格納する。

本開示によれば、信号の急変を精度よく検知することができる。

本開示の第１の実施形態の信号処理装置の機能構成の一例を示すブロック図である。本開示の第１の実施形態の信号の切り出し方の一例を示す図である。本開示の第１の実施形態のローカルスコアのイメージの一例を示す図である。本開示の第１の実施形態のスコアを算出するイメージの一例を示す図である。本開示の第１の実施形態の信号処理装置の動作の一例を示すフローチャートである。本開示の第２の実施形態の信号処理装置の機能構成の一例を示すブロック図である。本開示の第２の実施形態の信号処理装置の動作の一例を示すフローチャートである。本開示の第１及び第２の実施形態の信号処理装置を実現するコンピュータ装置のハードウェア構成の一例を示すブロック図である。

以下に、本開示の実施形態について、図面を参照しつつ説明する。

＜第１の実施形態＞
第１の実施形態の信号処理装置について説明する。

図１は、第１の実施形態の信号処理装置１００の機能構成の一例を示すブロック図である。図１に示すように、信号処理装置１００は、切出部１１０、変換部１２０、傾き算出部１３０、スコア算出部１４０、及び判定部１５０を備える。

切出部１１０は、入力信号から、所定区間の判定対象信号と、所定区間と一致しない区間であって、所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す。入力信号は、例えば、信号処理装置１００に入力される時系列信号であり、複数種類の信号が含まれる。切出部１１０は、入力信号のうち、急変信号が含まれるか否か判定を行う所定区間の信号である判定対象信号を切り出す。また、切出部１１０は、入力信号のうち、判定対象信号の区間と一致しない区間であって、判定対象信号の区間の少なくとも一部と重複する区間の信号である重複信号を切り出す。このとき、切出部１１０は、例えば、判定対象信号の先頭時刻、及び、判定対象信号の先頭時刻と重複信号の先頭時刻との時刻差を算出してもよい。切出部１１０は、切り出し手段の一例である。

変換部１２０は、判定対象信号と重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換する。変換部１２０は、変換手段の一例である。

傾き算出部１３０は、位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する。傾き算出部１３０は、例えば、位相傾きを、位相成分信号から周波数ごとに算出する。傾き算出部１３０は、例えば、位相傾きを、隣り合う周波数における位相の差を求めることにより算出してもよいし、他の方法で算出してもよい。傾き算出部１３０は、算出手段の一例である。

スコア算出部１４０は、複数の位相傾きに応じて、入力信号の急変性に関するスコアを算出する。スコア算出部１４０は、例えば、切出部１１０によって切り出された信号のそれぞれに対して、位相傾きの出現頻度に関する値を算出する。そして、スコア算出部１４０は、例えば、切り出された信号ごとに算出した位相傾きの出現頻度に関する値を統合してスコアを算出する。より具体的には、スコア算出部１４０は、例えば、算出した位相傾きの出現頻度に関する値のうち、同一時刻に対応する値を加算した結果をスコアとする。このように、位相傾きの出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、判定対象信号の区間と重複信号の区間とが重なる区間における位相傾きの出現頻度に関する値が強調される。なお、スコア算出部１４０は、算出した出現頻度に関する値に、出現頻度に応じた重み付けを行い、重み付けを行った出現頻度に関する値のうち、同一時刻に対応する値を加算した結果をスコアとしてもよい。スコア算出部１４０は、スコア算出手段の一例である。

判定部１５０は、スコアに基づいて、判定対象信号の急変の有無を判定する。例えば、判定部１５０は、スコアの最大値と、スコアが最大となる時刻とに基づいて、判定対象信号に急変信号が含まれるか否かを判定する。このとき、例えば、判定部１５０は、スコアの最大値が所定の閾値以上であり、スコアの最大となる時刻が判定対象信号の区間に含まれる場合に、判定対象信号に急変信号が含まれる、すなわち、判定対象信号に急変が有る、と判定する。判定部１５０は、判定手段の一例である。

このように、第１の実施形態の信号処理装置１００は、判定対象信号の区間と重複する区間の重複信号も用いて、入力信号の急変性に関するスコアを算出する。このときのスコアは、区間ごとに算出された位相傾きの出現頻度に関する値のうち、同一時刻に対応する値を統合することにより算出される。そして、信号処理装置１００は、当該スコアに基づいて判定対象信号の急変の有無を判定する。これにより、急変信号に起因する位相成分信号の傾きに関する情報が強調された上でスコアが算出されるので、第１の実施形態の信号処理装置１００は、信号の急変を精度よく検知することができる、という効果を奏する。

［信号処理装置１００の詳細］
次に、第１の実施形態の信号処理装置１００の各構成の一例について、より具体的に説明する。

切出部１１０は、信号処理装置１００に入力された入力信号を取得する。切出部１１０は、取得した入力信号から複数の信号を切り出す。図２は、信号の切り出し方の一例を示す図である。図２のｘ_ｉｎ（ｔ）は入力信号を示し、ｔは時刻を示す。なお、図２では、図の視認性を高めるために、入力信号の波形は省略している。

切出部１１０は、入力信号ｘ_ｉｎ（ｔ）から、判定対象信号ｘ_０（ｎ_０）を切り出す。このときｘ_０（ｎ_０）は、時刻ｔ_０からＮ_０（＞０）時間経過するまでの区間の信号である。さらに切出部１１０は、判定対象信号ｘ_０（ｎ_０）の少なくとも一部と重複する信号である重複信号を切り出す。図２では、判定対象信号ｘ_０（ｎ_０）と一部が重複する重複信号ｘ₁（ｎ_１）及びｘ_２（ｎ_２）の例を示している。ｘ_１（ｎ_１）は、入力信号に対して、時刻ｔ_０からτ_１（＜０）時間ずれた時刻からＮ₁（＞０）時間を経過するまでの区間の信号である。ｘ_２（ｎ_２）は、入力信号に対して、時刻ｔ_０からτ_２（＞０）時間ずれた時刻からＮ₂（＞０）時間を経過するまでの区間の信号である。本開示において、切出部１１０によって切り出された信号をｘ_ｄ（ｎ_ｄ）（０≦ｄ≦Ｄ）と表す。また、切り出された信号（判定対象信号及び重複信号）の各区間の時間の長さをＮ_ｄとし、切り出された信号の先頭時刻からｔ_０を引いた値、すなわち、切り出された信号と判定対象信号との先頭時刻の時刻差をτ_ｄと表す。このとき、τ_０＝０である。また、切り出された区間における時刻ｎ_ｄを（０≦ｎ_ｄ＜Ｎ_ｄ）と表す。図２では、重複信号の一部が判定対象信号と重複している例を示しているが、重複信号はこの例に限らない。重複信号は、判定対象信号と一致しなければよい。例えば、重複信号は、判定対象信号に包含されてもよいし、判定対象信号を包含していてもよい。なお、本開示において、切出部１１０によって切り出された信号の区間をフレームとも称する。また、図２では、重複信号を二つ切り出しているが、重複信号の数はこの例に限らない。重複信号の数は一以上であればよい。なお、重複信号の数が一である場合に比べ、複数の重複信号が切り出された場合は、第１の実施形態の信号処理装置１００は、急変信号に起因する位相成分信号の傾きに関する情報がより強調されるので、より精度よく信号の急変の有無を検知することができる。

切出部１１０は、例えば、所定の窓関数を用いることにより、入力信号から信号を切り出す。切出部１１０は、例えば、矩形窓を用いて信号を切り出すことができる。この例に限らず、切出部１１０は、ガウス窓、ハニング窓、またはハミング窓等、他の窓関数を用いて切り出してもよい。

変換部１２０は、切出部１１０によって切り出された信号を取得する。変換部１２０は、取得した信号、すなわち、判定対象信号及び重複信号のそれぞれを、周波数ごとに、位相成分信号に変換する。変換部１２０は、例えば、数１に示す離散フーリエ変換を用いて位相成分信号Ｘ_ｄ（ｋ）を算出する。

数１において、ｋ（０≦ｋ＜Ｋ－１：Ｋは自然数）は周波数を示す指標である。ｋの値のそれぞれと、信号に含まれる周波数のそれぞれとが対応している。なお、変換部１２０は、離散フーリエ変換に代えて、フーリエ変換またはウェーブレット変換等、他の変換方法で、位相成分信号を算出してもよい。

傾き算出部１３０は、位相成分信号に基づいて、位相成分の傾きを算出する。位相成分の傾きとは、隣り合う周波数における位相の変化の度合いを示す。本開示において、位相成分の傾きを、位相傾きとも称する。傾き算出部１３０は、位相成分信号Ｘ_ｄ（ｋ）から位相傾きΔθ_ｄ（ｋ）を求める際に、アンラップ位相を用いてもよい。以下の数２、数３、及び数４は、位相傾きを、アンラップ位相を用いて算出する例である。

数２は、各フレームに対して、周波数ごとの位相θ_ｄ（ｋ）を求める式である。位相は、例えば、数２のように位相成分信号を複素数で表現した場合、位相成分信号を示すベクトルにタンジェント（tan）の逆関数を適用することで求められる。

数３は、アンラップ位相を算出する式である。アンラップ位相は、位相の範囲がーπからπの間に限定されない絶対位相である。アンラップ位相は、隣り合う周波数における位相の差に応じて求められる。

数４に示すように、傾き算出部１３０は、隣り合う周波数におけるアンラップ位相の差分に基づいて、各周波数における位相傾きを算出する。

なお、傾き算出部１３０は、以下の数５、数６、数７を用いて位相傾きを算出してもよい。

例えば、傾き算出部１３０は、数５の式を用いて各周波数における位相成分信号を、単位ベクトルに正規化する。そして、傾き算出部１３０は、数６に示すように、隣り合う周波数における単位ベクトルの内積を算出する。これにより、複素平面上において、隣り合う周波数の単位ベクトルが合成されたベクトルが算出される。このとき、算出されたベクトルの位相が、周波数ｋに対する位相傾きに相当する。そのため、傾き算出部１３０は、数６において算出されたベクトルに対して、数７の式のように位相を求めることにより、位相傾きを算出することができる。

次に、スコア算出部１４０は、傾き算出部１３０から取得した位相傾きに応じて、入力信号の急変性に関するスコアを算出する。具体的には、スコア算出部１４０は、以下の数８及び数９に示す式を用いて、入力信号の急変性に関するスコアｙ（ｔ）を算出する。

数８は、あるフレームの時間ごとの、位相傾きの出現頻度を示すローカルスコアＳ_ｄ（ｎ_ｄ）を算出する式である。図３は、あるフレームにおけるローカルスコアＳ_ｄ（ｎ_ｄ）のイメージを表す図である。Ｓ_ｄ（ｎ_ｄ）の数が大きいほど、そのフレーム内時刻において、同じ位相傾きを持つ周波数、すなわち、そのフレームの時刻における位相傾きの出現頻度が多いことを示す。スコア算出部１４０は、図３に示すように、ローカルスコアをフレームごとに集計する。

次に、スコア算出部１４０は、数９の式を用いて、入力信号の急変性に関するスコアｙ（ｔ）を算出する。なお、本明細書において、入力信号の急変性に関するスコアを、単に「スコア」とも称する。

数９は、スコアｙ（ｔ）を算出する式である。スコアｙ（ｔ）は、ローカルスコアＳ_ｄ（ｎ_ｄ）のそれぞれを、ｔの時間軸に変換し、ｔの時間軸において同一時刻のローカルスコアを統合したものである。図４は、スコアを算出するイメージを表す図である。図４では、判定対象信号ｘ_０（ｎ_０）と、重複信号ｘ_１（ｎ_１）及びｘ_２（ｎ_２）が切り出された場合のスコアの一例を示している。図４では、それぞれの信号に対して、ローカルスコアＳ_０（ｎ_０）、Ｓ_１（ｎ_１）及びＳ_２（ｎ_２）が算出されている。そして、スコア算出部１４０は、ローカルスコアのうち、同一時刻に対応する値を足し合わせることにより、ローカルスコアの統合を行う。

このように、スコア算出部１４０は、判定対象信号及び重複信号のそれぞれにおいて、位相傾きの出現頻度に関する値(例えば、ローカルスコアＳ_ｄ（ｎ_ｄ）)を算出し、算出した当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、スコアを算出する。

判定部１５０は、スコア算出部１４０によって算出されたスコアに基づいて、判定対象信号に急変信号が含まれているか否かを判定する。

数１０は、判定方法を示す式である。数１０では、ｐ＝１のとき、判定部１５０が、判定対象信号に急変信号が含まれると判定したことが示され、ｐ＝０のとき、判定部１５０が、判定対象信号に急変信号が含まれないと判定したことが示されている。判定部１５０は、数１０の式を用いて、スコアの最大値が閾値α以上であり、スコアが最大のときの時刻が判定対象信号のフレームの間の時間である場合に、判定対象信号に急変信号が含まれると判定する。

このように、判定部１５０は、スコアの最大値とスコアが最大となる時刻とに基づいて判定対象信号の急変の有無を判定する。

［信号処理装置１００の動作］
次に信号処理装置１００の動作を、図５を用いて説明する。図５は、信号処理装置１００の動作を説明するフローチャートである。以下に説明する動作においては、信号処理装置１００に、入力信号が入力されたものとする。なお、本明細書において、フローチャートの各ステップを「Ｓ１０１」のように、それぞれのステップに付した番号を用いて表現する。

切出部１１０は、入力信号から、判定対象信号と重複信号とを切り出す（Ｓ１０１）。変換部１２０は、判定対象信号と重複信号とのそれぞれに対して、周波数領域における位相成分信号に変換する（Ｓ１０２）。次に、傾き算出部１３０は、位相成分信号に基づいて、位相傾きを算出する（Ｓ１０３）。そして、スコア算出部１４０は、算出された位相傾きに応じて、入力信号の急変性に関するスコアを算出する（Ｓ１０４）。判定部１５０は、算出されたスコアに基づいて、判定対象信号の急変の有無を判定する（Ｓ１０５）。

このように、第１の実施形態の信号処理装置１００は、入力信号から、所定区間の判定対象信号と、所定区間と一致しない区間であって、所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す。また、信号処理装置１００は、判定対象信号と重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換し、位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する。そして、信号処理装置１００は、複数の位相傾きに応じて、入力信号の急変性に関するスコアを算出し、スコアに基づいて、判定対象信号の急変の有無を判定する。このように、第１の実施形態の信号処理装置１００は、判定対象信号と、判定対象信号の区間と一部重複する重複信号とのそれぞれから得られる位相傾きに応じて、スコアを算出しているので、入力信号において、急変信号対非急変信号比が小さい場合であっても、急変信号に起因する位相成分信号の傾きに関する情報を強調したスコアを算出することができる。そのため、第１の実施形態の信号処理装置１００は、信号の急変を精度よく検知することができる。

［変形例１］
スコア算出部１４０は、位相傾きの出現頻度に関する値Ｓ_ｄ（ｎ_ｄ）に対して、出現頻度に応じた重み付けを行ってもよい。具体的には、スコア算出部１４０は、数１１のように、ｗ_ｄを用いて重み付けを行ってもよい。

ｗ_ｄは、例えば、Ｓ_ｄ（ｎ_ｄ）における分散に基づいて求められる。Ｓ_ｄ（ｎ_ｄ）に対して分散が大きいほど、そのフレームの時間ｎにおいて、急変信号が含まれる可能性が高い。そのため、スコア算出部１４０は、数１２を用いてｗ_ｄを算出してもよい。

数１２のＳ_ａｖｅは、Ｓ_ｄ（ｎ_ｄ）の平均である。このとき、Ｓ_ｄ（ｎ_ｄ）に対する分散が大きいほど、ｗ_ｄの値は大きくなる。

また、ｗ_ｄは、ｎ_ｄに対する分散に基づいて求められてもよい。ｎ_ｄに対して分散が小さいほど、一部の時間に、Ｓ_ｄ（ｎ_ｄ）の値が偏って算出される。すなわち、その時間に急変信号が含まれる可能性が高い。そのため、スコア算出部１４０は、数１３を用いて、ｗ_ｄを算出してもよい。

数１３のｎ_ａｖｅは、ｎ_ｄの平均である。このとき、ｎ_ｄに対する分散が小さいほど、ｗ_ｄの値は大きくなる。

このように、信号処理装置１００は、位相傾きの出現頻度に関する値に対して、出現頻度に応じた重み付けを行い、重み付けを行った出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、スコアを算出してもよい。この構成により、変形例１の信号処理装置１００は、急変信号が含まれる可能性の高い時刻における、出現頻度に関する値を強調させることができるので、より精度よく信号の急変の有無を検知することができる。

［変形例２］
判定部１５０は、スコアの推定分布値を用いて、判定対象信号に、急変信号が含まれるか否かを判定してもよい。具体的には、スコアｙ（ｔ）に対して、正規分布に従う確率密度関数をフィッティングさせ、そのときの確率密度関数のパラメータから、判定対象信号に急変信号が含まれるか否かを判定する。

数１４は、正規分布に従う確率密度関数ｇ（ｔ）を示す式である。数１５は、スコアｙ（ｔ）と確率密度関数ｇ（ｔ）との差を示す関数の式である。判定部１５０は、確率密度関数ｇ（ｔ）を数９あるいは数１１に示すスコアｙ（ｔ）にフィッティングさせる。すなわち、判定部１５０は、関数ｆ（μ，σ）が最小となるときの平均μ及び標準偏差σの値を算出する。そして、判定部１５０は、算出されたμ及びσを用いて、以下の数１６の式に基づき、判定を行う。

判定部１５０は、例えば、ｇ（μ）×１／σの値が閾値α以上であり、μの値が、判定対象信号のフレームの間の値のとき、判定対象信号に急変信号が含まれると判定する。

ここで、ｇ（μ）は、スコアｙ（ｔ）の最大値に対応する値である。また、ｇ（ｔ）に関して、σの値が大きいほど、ｇ（ｔ）における分散は大きい。すなわち、σの値が大きいほど、ｇ（ｔ）はなだらかな変化を行う関数となる。ｇ（ｔ）がなだらかな変化を行う場合、ｇ（μ）が閾値αを超えていても、急変信号が判定対象信号に含まれる可能性は低くなる。そのため、判定部１５０は、σの値に応じて、閾値αと比較する値を変換する。

このように、変形例２の信号処理装置１００は、スコアに対してフィッティングした確率密度関数における最大値を、確率密度関数の分散に応じて変換した値と、確率密度関数の値が最大となる時刻と、に基づいて、判定対象信号の急変の有無を判定する。この構成により、変形例２の信号処理装置１００は、判定対象信号の急変の有無を判定する際に、スコアの最大値を用いるときよりも、精度よく判定することができる。

＜第２の実施形態＞
次に、第２の実施形態の信号処理装置について説明する。

図６は、第２の実施形態の信号処理装置１０１の機能構成の一例を示すブロック図である。図６に示す通り、信号処理装置１０１は、第１の実施形態の信号処理装置１００の構成に加え、信頼度算出部１６０を備える。すなわち、信号処理装置１０１は、切出部１１０と、変換部１２０と、傾き算出部１３０と、スコア算出部１４０と、判定部１５０と、信頼度算出部１６０と、を備える。なお、図６に示す信号処理装置１０１の構成及び動作が、第１の実施形態の説明と重複する内容については、説明を省略する。

［信号処理装置１０１の詳細］
判定部１５０によって判定対象信号に急変があると判定された場合、信頼度算出部１６０は、スコアに応じて、判定対象信号の各周波数に対して、急変を示す信号（すなわち、急変信号）の周波数である可能性を示す信頼度を算出する。例えば、信頼度算出部１６０は、急変信号の位相傾きとして推定される値を、判定対象信号の周波数ごとの位相傾きと比較する。そして、信頼度算出部１６０は、判定対象信号の周波数ごとに、急変信号の周波数である可能性を示す値を、信頼度として算出する。

数１７は、急変信号の位相傾きの推定値を算出する式である。信頼度算出部１６０は、例えば、数１７の式を用いて、急変信号の位相傾きの推定値を算出する。このとき、信頼度算出部１６０は、スコアｙ（ｔ）が最大値となるときの時刻における、入力信号の位相傾きを、急変信号の位相傾きの推定値として算出する。

そして、信頼度算出部１６０は、数１８の式を用いて、急変信号の位相傾きの推定値と、判定対象信号の位相傾きとの差を、判定対象信号の周波数ごとに求めることにより、信頼度Ｙ（ｋ）を算出する。このとき、推定値と判定対象信号の位相傾きとの差が小さいほど、比較された判定対象信号の位相傾きをもつ周波数が、急変信号の周波数である可能性が高くなる。一方で、推定値と判定対象信号の位相傾きとの差が大きいほど、比較された判定対象信号の位相傾きをもつ周波数は、急変信号の周波数でない可能性が高くなる、すなわち、急変信号でない信号の周波数である可能性が高くなる。

このように、信頼度算出部１６０は、判定対象信号に急変があると判定された場合、判定対象信号の周波数ごとの位相傾きと、スコアに応じて算出された、急変を示す信号の位相傾きの推定値と、に基づいて、判定対象信号の各周波数に対して、急変を示す信号の周波数である可能性を示す信頼度を算出する。信頼度算出部１６０は、信頼度算出手段の一例である。

［信号処理装置１０１の動作］
次に、信号処理装置１０１の動作を、図７を用いて説明する。図７は、信号処理装置１０１の動作の一例を説明するフローチャートである。なお、図７において、Ｓ１０１乃至Ｓ１０５の処理は、図５におけるＳ１０１乃至Ｓ１０５の処理と同様であるため、その説明を省略する。

信号処理装置１０１は、Ｓ１０５の処理において、判定対象信号に急変無し、と判定された場合（Ｓ１０６の「Ｎｏ」）、フローを終了する。判定対象信号に急変ありと判定された場合（Ｓ１０６の「Ｙｅｓ」）、信頼度算出部１６０は、スコアに応じて算出された急変信号の位相傾きの推定値と、判定対象信号の周波数ごとの位相傾きと、に基づいて、判定対象信号の周波数ごとに信頼度を算出する（Ｓ１０７）。

以上のように、第２の実施形態の信号処理装置１０１は、判定対象信号に急変があると判定された場合、判定対象信号の周波数ごとの位相傾きと、スコアに応じて算出された、急変を示す信号の位相傾きの推定値と、に基づいて、判定対象信号の各周波数に対して、急変を示す信号の周波数である可能性を示す信頼度を算出する。この構成により、第２の実施形態の信号処理装置１０１は、判定対象信号において、どの周波数の信号が、急変信号に起因する周波数の信号である可能性が高いか、または、急変信号でない信号に起因する周波数の信号である可能性が高いか判定することができる。この判定結果により、例えば、ユーザは、判定対象信号から非急変信号を取り除いたり、反対に、判定対象信号から急変信号を取り除いたりすることができる。

＜信号処理装置のハードウェアの構成例＞
上述した第１、及び第２の実施形態の信号処理装置を構成するハードウェアについて説明する。図８は、各実施形態における信号処理装置を実現するコンピュータ装置のハードウェア構成の一例を示すブロック図である。図８が示す各ブロックは、各実施形態における信号処理装置及び信号処理方法を実現するコンピュータ装置１０と、ソフトウェアとの組み合わせにより実現できる。

図８に示すように、コンピュータ装置１０は、プロセッサ１１、ＲＡＭ（ＲａｎｄｏｍＡｃｃｅｓｓＭｅｍｏｒｙ）１２、ＲＯＭ（ＲｅａｄＯｎｌｙＭｅｍｏｒｙ）１３、記憶装置１４、入出力インタフェース１５、バス１６、及びドライブ装置１７を備える。なお、信号処理装置は、複数の電気回路によって実現されてもよい。

記憶装置１４は、プログラム（コンピュータプログラム）１８を格納する。プロセッサ１１は、ＲＡＭ１２を用いて本信号処理装置のプログラム１８を実行する。具体的には、例えば、プログラム１８は、図５及び図６で説明した信号処理装置の処理をコンピュータに実行させるプログラムを含む。プロセッサ１１が、プログラム１８を実行することに応じて、本信号処理装置の各構成要素（上述した、切出部１１０、変換部１２０、傾き算出部１３０、スコア算出部１４０、判定部１５０、信頼度算出部１６０）の機能が実現される。プログラム１８は、ＲＯＭ１３に記憶されていてもよい。また、プログラム１８は、記憶媒体２０に記録され、ドライブ装置１７を用いて読み出されてもよいし、図示しない外部装置から図示しないネットワークを介してコンピュータ装置１０に送信されてもよい。

入出力インタフェース１５は、周辺機器（キーボード、マウス、表示装置など）１９とデータをやり取りする。入出力インタフェース１５は、データを取得または出力する手段として機能する。バス１６は、各構成要素を接続する。

なお、信号処理装置の実現方法には様々な変形例がある。例えば、信号処理装置は、専用の装置として実現することができる。また、信号処理装置は、複数の装置の組み合わせに基づいて実現することができる。

各実施形態の機能における各構成要素を実現するためのプログラムを記憶媒体に記録させ、該記憶媒体に記録されたプログラムをコードとして読み出し、コンピュータにおいて実行する処理方法も各実施形態の範疇に含まれる。すなわち、コンピュータ読取可能な記憶媒体も各実施形態の範囲に含まれる。また、上述のプログラムが記録された記憶媒体、及びそのプログラム自体も各実施形態に含まれる。

該記憶媒体は、例えばフロッピー（登録商標）ディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ（ＣｏｍｐａｃｔＤｉｓｃ）－ＲＯＭ、磁気テープ、不揮発性メモリカード、またはＲＯＭであるが、この例に限らない。また該記憶媒体に記録されたプログラムは、単体で処理を実行しているプログラムに限らず、他のソフトウェア、拡張ボードの機能と共同して、ＯＳ（ＯｐｅｒａｔｉｎｇＳｙｓｔｅｍ）上で動作して処理を実行するプログラムも各実施形態の範疇に含まれる。

以上、実施形態を参照して本願発明を説明したが、本願発明は上記実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で当業者が理解しうる様々な変更をすることができる。

上記の実施形態の一部または全部は、以下の付記のようにも記載されうるが、以下には限られない。

＜付記＞
［付記１］
入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す切出手段と、
前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換する変換手段と、
前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する傾き算出手段と、
複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出するスコア算出手段と、
前記スコアに基づいて、前記判定対象信号の急変の有無を判定する判定手段と、を備える、
信号処理装置。

［付記２］
前記スコア算出手段は、前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
付記１に記載の信号処理装置。

［付記３］
前記スコア算出手段は、前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値に対して、前記出現頻度に応じた重み付けを行い、重み付けを行った当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
付記１に記載の信号処理装置。

［付記４］
前記判定手段は、前記スコアの最大値と前記スコアが最大となる時刻とに基づいて前記判定対象信号の急変の有無を判定する、
付記１乃至３のいずれかに記載の信号処理装置。

［付記５］
前記判定手段は、前記スコアに対してフィッティングした確率密度関数における最大値を、前記確率密度関数の分散に応じて変換した値と、前記確率密度関数の値が最大となる時刻と、に基づいて、前記判定対象信号の急変の有無を判定する、
付記１乃至３のいずれかに記載の信号処理装置。

［付記６］
前記切出手段は、複数の前記重複信号を切り出す、
付記１乃至５のいずれかに記載の信号処理装置。

［付記７］
前記判定対象信号に急変があると判定された場合、前記判定対象信号の周波数ごとの位相傾きと、前記スコアに応じて算出された、前記急変を示す信号の位相傾きの推定値と、に基づいて、前記判定対象信号の各周波数に対して、前記急変を示す信号の周波数である可能性を示す信頼度を算出する信頼度算出手段を備える、
付記１乃至６のいずれかに記載の信号処理装置。

［付記８］
入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出し、
前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換し、
前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出し、
複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出し、
前記スコアに基づいて、前記判定対象信号の急変の有無を判定する、
信号処理方法。

［付記９］
前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
付記８に記載の信号処理方法。

［付記１０］
前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値に対して、前記出現頻度に応じた重み付けを行い、重み付けを行った当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
付記８に記載の信号処理方法。

［付記１１］
前記スコアの最大値と前記スコアが最大となる時刻とに基づいて前記判定対象信号の急変の有無を判定する、
付記８乃至１０のいずれかに記載の信号処理方法。

［付記１２］
前記スコアに対してフィッティングした確率密度関数における最大値を、前記確率密度関数の分散に応じて変換した値と、前記確率密度関数の値が最大となる時刻と、に基づいて、前記判定対象信号の急変の有無を判定する、
付記８乃至１０のいずれかに記載の信号処理方法。

［付記１３］
複数の前記重複信号を切り出す、
付記８乃至１２のいずれかに記載の信号処理方法。

［付記１４］
前記判定対象信号に急変があると判定された場合、前記判定対象信号の周波数ごとの位相傾きと、前記スコアに応じて算出された、前記急変を示す信号の位相傾きの推定値と、に基づいて、前記判定対象信号の各周波数に対して、前記急変を示す信号の周波数である可能性を示す信頼度を算出する、
付記８乃至１３のいずれかに記載の信号処理方法。

［付記１５］
入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す処理と
前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換する処理と
前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する処理と、
複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出する処理と、
前記スコアに基づいて、前記判定対象信号の急変の有無を判定する処理と、をコンピュータに実行させるプログラムを格納する、
コンピュータ読み取り可能な記憶媒体。

［付記１６］
前記スコアを算出する処理において、前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
付記１５に記載のコンピュータ読み取り可能な記憶媒体。

［付記１７］
前記スコアを算出する処理において、前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値に対して、前記出現頻度に応じた重み付けを行い、重み付けを行った当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
付記１５に記載のコンピュータ読み取り可能な記憶媒体。

［付記１８］
前記判定する処理において、前記スコアの最大値と前記スコアが最大となる時刻とに基づいて前記判定対象信号の急変の有無を判定する、
付記１５乃至１７のいずれかに記載のコンピュータ読み取り可能な記憶媒体。

［付記１９］
前記判定する処理において、前記スコアに対してフィッティングした確率密度関数における最大値を、前記スコアの確率密度関数に応じて変換した値と、前記確率密度関数の値が最大となる時刻と、に基づいて、前記判定対象信号の急変の有無を判定する、
付記１５乃至１７のいずれかに記載のコンピュータ読み取り可能な記憶媒体。

［付記２０］
前記切り出す処理において、複数の前記重複信号を切り出す、
付記１５乃至１９のいずれかに記載のコンピュータ読み取り可能な記憶媒体。

［付記２１］
前記判定対象信号に急変があると判定された場合、前記判定対象信号の周波数ごとの位相傾きと、前記スコアに応じて算出された、前記急変を示す信号の位相傾きの推定値と、に基づいて、前記判定対象信号の各周波数に対して、前記急変を示す信号の周波数である可能性を示す信頼度を算出する処理をさらに実行させるプログラムを格納する、
付記１５乃至２０のいずれかに記載のコンピュータ読み取り可能な記憶媒体。

１００、１０１信号処理装置
１１０切出部
１２０変換部
１３０傾き算出部
１４０スコア算出部
１５０判定部
１６０信頼度算出部

Claims

入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す切出手段と、
前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換する変換手段と、
前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する傾き算出手段と、
複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出するスコア算出手段と、
前記スコアに基づいて、前記判定対象信号の急変の有無を判定する判定手段と、を備える、
信号処理装置。
前記スコア算出手段は、前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
請求項１に記載の信号処理装置。
前記スコア算出手段は、前記判定対象信号及び前記重複信号のそれぞれにおいて、前記位相傾きの出現頻度に関する値を算出し、算出した当該出現頻度に関する値に対して、前記出現頻度に応じた重み付けを行い、重み付けを行った当該出現頻度に関する値のうち、同一時刻に対応する値を統合することにより、前記スコアを算出する、
請求項１に記載の信号処理装置。
前記判定手段は、前記スコアの最大値と前記スコアが最大となる時刻とに基づいて前記判定対象信号の急変の有無を判定する、
請求項１乃至３のいずれかに記載の信号処理装置。
前記判定手段は、前記スコアに対してフィッティングした確率密度関数における最大値を、前記確率密度関数の分散に応じて変換した値と、前記確率密度関数の値が最大となる時刻と、に基づいて、前記判定対象信号の急変の有無を判定する、
請求項１乃至３のいずれかに記載の信号処理装置。
前記切出手段は、複数の前記重複信号を切り出す、
請求項１乃至５のいずれかに記載の信号処理装置。
前記判定対象信号に急変があると判定された場合、前記判定対象信号の周波数ごとの位相傾きと、前記スコアに応じて算出された、前記急変を示す信号の位相傾きの推定値と、に基づいて、前記判定対象信号の各周波数に対して、前記急変を示す信号の周波数である可能性を示す信頼度を算出する信頼度算出手段を備える、
請求項１乃至６のいずれかに記載の信号処理装置。
入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出し、
前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換し、
前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出し、
複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出し、
前記スコアに基づいて、前記判定対象信号の急変の有無を判定する、
信号処理方法。
入力信号から、所定区間の判定対象信号と、前記所定区間と一致しない区間であって、前記所定区間の少なくとも一部と重複する区間の重複信号と、を切り出す処理と
前記判定対象信号と前記重複信号とのそれぞれを、周波数ごとに、周波数領域における位相成分信号に変換する処理と
前記位相成分信号に基づいて、周波数における位相成分の傾きである位相傾きを、各周波数において算出する処理と、
複数の前記位相傾きに応じて、前記入力信号の急変性に関するスコアを算出する処理と、
前記スコアに基づいて、前記判定対象信号の急変の有無を判定する処理と、をコンピュータに実行させるプログラム。