(57)【特許請求の範囲】
【請求項１】
現時刻からｎ−１時刻前までn個のサンプルデータを入力するサンプルデータ入力手段と、
入力された前記サンプルデータを蓄積するサンプルデータ蓄積手段と、
前記サンプルデータ蓄積手段に蓄積されているデータのうち、最古のサンプルデータ列と、入力された最新のサンプルデータ列に階級番号を付与する階級番号付与手段と、
度数の分布状況、階級番号の分布状況、最も前記サンプルデータの度数が多い最大値及び最も度数が多い階級番号を管理する記憶手段と、
前記サンプルデータが入力された時点において、最古のサンプルデータ列の階級番号と、最新のサンプルデータ列の階級番号を比較して異なる場合に、前記記憶手段の内容を更新し、最大値及び最頻値を取得する最大値・最頻値求解手段と、を有し、
前記階級番号付与手段は、
前記サンプルデータが属する階級に対応する前記階級番号を参照するためのサンプルデータ−階級番号変換値記憶手段を有し、
前記記憶手段は、
現時刻からn―１時刻までのサンプルデータにおける前記階級番号毎の該サンプルデータの度数の分布を示す階級番号別度数分布記憶手段と、
前記サンプルデータの度数別の階級番号の分布を示す度数別階級番号分布記憶手段と、
前記階級番号別度数分布記憶手段において、最も大きいサンプルデータの度数である最大値を管理する最大値記憶手段と、
前記サンプルデータの度数が最も大きい階級番号を管理する最頻値階級番号記憶手段と、
前記度数別階級番号分布記憶手段の更新を管理する更新管理手段と、
前記階級番号別度数分布記憶手段の前記階級番号毎の階級値を管理する階級値記憶手段と、を有し、
前記最大値・最頻値求解手段は、
前記度数別階級番号分布記憶手段が前記度数毎に管理する前記階級番号において、先頭の階級番号が保存される前記度数別階級番号分布記憶手段内の位置を度数ごとに表す度数別先頭位置記憶手段と、
前記度数別階級番号分布記憶手段における前記度数毎に前記階級番号の有効数を管理する度数別階級番号有効数記憶手段と、を有する更新管理手段を有する
ことを特徴とする統計処理装置。
【請求項２】
サンプルデータ入力手段において、現時刻からｎ−１時刻前までn個のサンプルデータを入力し、
サンプルデータ蓄積手段において、前記サンプルデータを蓄積し、記憶手段において、該サンプルデータの度数と前記階級番号を蓄積しておき、
最大値・最頻値求解手段において、新たにサンプリングデータが入力されると、それまで前記記憶手段に蓄積されていた前記サンプルデータの度数と前記階級番号からなるデータのうち、最古のサンプルデータを廃棄し、入力された該サンプルデータを最新のサンプルデータとして追加することにより該記憶手段を更新し、
階級番号付与手段において、前記サンプルデータ蓄積手段に蓄積されている最古のサンプルデータ列と最新のサンプルデータ列に階級番号を付与し、
前記最大値・最頻値求解手段において、更新により最も度数が多い最大値と、最も度数が大きい階級番号である最頻値が変化した場合に、変化した最大値と最頻値を取得し、
前記階級番号付与手段において、前記階級番号を前記サンプルデータ毎に付与する際に、
前記サンプルデータが属する階級に対応する前記階級番号を付与するためのサンプルデータ―階級番号変換テーブルを参照する処理を行い、
前記最大値・最頻値求解手段において、前記最大値・最頻値を取得する際に、
現時刻からn−1時刻までのサンプルデータにおける階級番号別に度数の第１の分布状況を取得し、
前記サンプルデータの度数別の階級番号の第２の分布状況、前記階級番号別に度数の第３の分布を示す状況を取得し、
新たに、サンプルデータが入力された場合に、前記第１、第２、第３の分布状況により前記最大値と前記最頻値が変化したことを確認する処理を行うことを特徴とする統計処理方法。
【請求項３】
前記最大値・最頻値求解手段において、前記最大値・最頻値を取得する際に、
前記度数毎に管理する前記階級番号において、度数毎に、先頭の階級番号が保存される位置と、
前記度数毎に前記階級番号の有効数を参照する請求項２記載の統計処理方法。
【請求項４】
前記サンプルデータ入力手段において、サンプリング対象のサンプルデータを入力し、入力された該サンプルデータに固有の番号である入力点を付与し、
前記サンプルデータ蓄積手段において、前記入力点毎に前記サンプルデータを管理し、前記階級番号付与手段において、最新のサンプルデータ列と最古のサンプルデータ列に階級番号を付与して、最新のサンプルデータの階級番号列と、最古のサンプルデータの階級番号列として出力し、
前記最大値・最頻値求解手段において、前記最新のサンプルデータの階級番号列と、最古のサンプルデータの階級番号列を比較して、前記最新のサンプルデータの階級番号列と、最古のサンプルデータの階級番号列が異なる入力点がある時のみ、前記サンプルデータの度数が最も大きい最大値と、該度数が最も大きい階級値である最頻値を求める請求項２記載の統計処理方法。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、統計処理装置及び統計処理方法に係り、特に、人、車、ベルトコンベア上の物体等の移動物を計数する分野に属し、動画像の各画素毎の輝度、物体の質量、変位等の時々刻々変化するデータを連続サンプリングし、現時刻からｎ−１時刻前まで、ｎ個のサンプルデータ、または、連続的にサンプリングしたｎのデータから連続的に最大値、最頻値を求める統計処理装置及び統計処理方法に関する。
【０００２】
【従来の技術】
図１８は、従来の統計処理装置の構成を示す。
従来の統計処理装置は、現時刻からｎ時刻前までのサンプルデータを順に入力するサンプルデータ入力回路１と、サンプルデータを入力順に蓄積するリングバッファテーブル５を用いてサンプルデータを管理するリングバッファ回路２と、リングバッファテーブル５において現時刻からｎ−１時刻前までの各サンプルデータの属する階級番号毎に度数を管理する階級番号別度数テーブル６より階級番号を取得する階級番号取得部３と階級番号別度数テーブル６から最大値、最頻値を求める最大値・最頻値取得部４を有する。
【０００３】
図１９は、従来の統計処理のフローチャートである。
まず、統計装置のサンプリングデータ入力回路１は、サンプルデータをサンプリングし（ステップ１０）、ｎ回のサンプリングの後に、リングバッファテーブル５と階級番号別度数テーブル６の初期設定を行い（ステップ１１）、リンクバッファ回路２は、ｎ＋１回以降の（ステップ１２）サンプリング毎に得られるサンプルデータにより（ステップ１３）リングバッファテーブル５を更新し（ステップ１４）、階級番号取得部３は、最古と最新のサンプルデータの階級番号求解処理を行い（ステップ１５）、階級番号別度数テーブル６を更新し（ステップ１６）、最大値、最頻値取得部４は、更新後の階級番号別度数テーブル６から最大値、最頻値を求める最大値・最頻値求解処理を繰り返し実施する（ステップ１７）。
【０００４】
【発明が解決しようとする課題】
しかしながら、上記従来の統計装置は、サンプリングごとに行う階級番号求解処理でリングバッファテーブルから得られる最古、最新サンプルデータの属する階級番号を求めるために、最古、最新サンプルデータと入力されたサンプリングデータの階級を逐次比較しなければならず、処理が繁雑である。
【０００５】
また、最大値・最頻値求解処理は、更新した階級番号別度数テーブルから最も大きい度数（最大値）、最も大きい度数を与える階級の階級値（最頻値）をサンプル毎に求めるために、階級番号別度数テーブルにおいて階級番号毎の度数を全て比較する必要があり、１サンプリング毎の処理時間が多くかかり、実時間処理を行うためには、階級の数、サンプリング対象数、サンプリング間隔等に制約が生じるため、高速サンプリング、多種類サンプリングにおける最大値、最頻値求解処理実施が困難であるという問題がある。
【０００６】
本発明は、上記の点に鑑みなされたもので、時々刻々変化する多種類のデータを高速連続サンプリングし、現時刻からｎ−１時刻前までｎ個のサンプルデータから連続的に最大値、最頻値の求解処理が可能な最小構成で実時間処理を可能とする統計処理装置及び統計処理方法を提供することを目的とする。
【０００７】
【課題を解決するための手段】
図１は、本発明の原理構成図である。
本発明の統計処理装置は、現時刻からｎ−１時刻前までn個のサンプルデータを入力するサンプルデータ入力手段３０２と、入力されたサンプルデータを蓄積するサンプルデータ蓄積手段３０３と、サンプルデータ蓄積手段３０３に蓄積されているデータのうち、最古のサンプルデータ列と、入力された最新のサンプルデータ列に階級番号を付与する階級番号付与手段３０４と、度数の分布状況、階級番号の分布状況、最もサンプルデータの度数が多い最大値及び最も度数が多い階級番号を管理する記憶手段３０６と、サンプルデータが入力された時点において、最古のサンプルデータ列の階級番号と、最新のサンプルデータ列の階級番号を比較して異なる場合に、記憶手段３０６の内容を更新し、最大値及び最頻値を取得する最大値・最頻値求解手段３０５と、を有し、
階級番号付与手段３０４は、
サンプルデータが属する階級に対応する階級番号を参照するためのサンプルデータ−階級番号変換値記憶手段を有し、
記憶手段３０６は、
現時刻からn―１時刻までのサンプルデータにおける階級番号毎の該サンプルデータの度数の分布を示す階級番号別度数分布記憶手段と、
サンプルデータの度数別の階級番号の分布を示す度数別階級番号分布記憶手段と、
階級番号別度数分布記憶手段において、最も大きいサンプルデータの度数である最大値を管理する最大値記憶手段と、
サンプルデータの度数が最も大きい階級番号を管理する最頻値階級番号記憶手段と、
度数別階級番号分布記憶手段の更新を管理する更新管理手段と、
階級番号別度数分布記憶手段の階級番号毎の階級値を管理する階級値記憶手段と、を有し、
最大値・最頻値求解手段３０５は、
度数別階級番号分布記憶手段が度数毎に管理する階級番号において、先頭の階級番号が保存される度数別階級番号分布記憶手段内の位置を度数ごとに表す度数別先頭位置記憶手段と、
度数別階級番号分布記憶手段における度数毎に階級番号の有効数を管理する度数別階級番号有効数記憶手段と、を有する更新管理手段を有する。
【００１０】
図２は、本発明の原理を説明するためのフローチャートである。
本発明の統計処理方法は、サンプルデータ入力手段において、現時刻からｎ−１時刻前までn個のサンプルデータを入力し（ステップ１）、サンプルデータ蓄積手段において、サンプルデータを蓄積し、記憶手段において、サンプルデータの度数と階級番号を蓄積しておき（ステップ２）、最大値・最頻値求解手段において、新たにサンプリングデータが入力されると（ステップ４）、それまで記憶手段に蓄積されていたサンプルデータの度数と階級番号からなるデータのうち、最古のサンプルデータを廃棄し、入力された該サンプルデータを最新のサンプルデータとして追加することにより該記憶手段を更新し（ステップ５）、
階級番号付与手段において、サンプルデータ蓄積手段に蓄積されている最古のサンプルデータ列と最新のサンプルデータ列に階級番号を付与し（ステップ６）、最大値・最頻値求解手段において、更新により最も度数が多い最大値と、最も度数が大きい階級番号である最頻値が変化した場合に（ステップ７、Ｙes）、変化した最大値と最頻値を取得し（ステップ８）、
階級番号付与手段において、階級番号をサンプルデータ毎に付与する際に（ステップ６）、
サンプルデータが属する階級に対応する階級番号を付与するためのサンプルデータ―階級番号変換テーブルを参照する処理を行い、
最大値・最頻値求解手段において、最大値・最頻値を取得する際に（ステップ８）、
現時刻からn−1時刻までのサンプルデータにおける階級番号別に度数の第１の分布状況を取得し、
サンプルデータの度数別の階級番号の第２の分布状況、階級番号別に度数の第３の分布を示す状況を取得し、
新たに、サンプルデータが入力された場合に、第１、第２、第３の分布状況により最大値と最頻値が変化したことを確認する処理を行う。
【００１２】
また、最大値・最頻値求解手段において、最大値・最頻値を取得する際に（ステップ８）、
度数毎に管理する階級番号において、度数毎に、先頭の階級番号が保存される位置と、度数毎に階級番号の有効数を参照する。
本発明の統計処理方法は、サンプルデータ入力手段において、サンプリング対象のサンプルデータを入力し、入力された該サンプルデータに固有の番号である入力点を付与し、
サンプルデータ蓄積手段において、入力点毎にサンプルデータを管理し、階級番号付与手段において、最新のサンプルデータ列と最古のサンプルデータ列に階級番号を付与して、最新のサンプルデータの階級番号列と、最古のサンプルデータの階級番号列として出力し、
最大値・最頻値求解手段において、最新のサンプルデータの階級番号列と、最古のサンプルデータの階級番号列を比較して、最新のサンプルデータの階級番号列と、最古のサンプルデータの階級番号列が異なる入力点がある時のみ、サンプルデータの度数が最も大きい最大値と、該度数が最も大きい階級値である最頻値を求める。
【００１３】
【作用】
本発明は、サンプリング毎に得られる最古、最新サンプルデータの属する階級番号を各階級と全て比較せずに、サンプルデータの度数に応じて階級番号に変換することにより、最古サンプルデータ及び、最新サンプルデータが属する階級に対応する階級番号を取得し、この階級番号同士を比較する。これにより、処理時間を短縮することが可能となる。
【００１４】
また、本発明は、階級番号別度数分布表の各階級番号毎の度数を全て比較せずに、度数別に階級番号の分布状況を参照して、最大値と最頻値が変動する条件の時のみ、最大値、最頻値を求めることにより、処理時間を短縮することができる。即ち、本発明は、サンプリング後の最古サンプルデータを破棄することに伴い最大値、最頻値の変動に関する条件と、最新サンプルデータを取得することによる最大値、最頻値の変動に関する条件について監視し、最大値、最頻値の変動条件時のみ、最大値、最頻値の求解処理を実施するものである。
【００１５】
さらに、最大値、最頻値の変動条件に応じた領域の管理構成とすることにより、使用メモリの容量を最小化する。
【００１６】
【実施例】
以下、本発明の実施例を図面と共に説明する。
図３は、本発明の一実施例の統計装置を適用したシステム構成を示す。同図において、対象物の重量を計測する重量センタ１１０とアナログ信号をディジタル信号に変換するＡ／Ｄ変換器１２０を有するデータ入力装置１００、統計処理装置３００、最頻値、最大値を表示する表示装置２００より構成される。
【００１７】
図４は、本発明の一実施例の統計装置の使用例を示す。同図において、ベルトコンベア１０上を移動する生産物２０を重量測定エリア３０において、データ入力装置１００がサンプリングしている状態を示す。サンプリングは、データ入力装置１００の重量センサ１１０において、重量のサンプリングを行う他に、画像処理における動画像の画素毎に輝度、色をサンプリングの対象とすることも可能である。この場合には、輝度、色がサンプルデータとなる。
【００１８】
一般には、統計処理装置３００への入力は、データ入力装置１００で計測されたディジタル信号が逐次入力される。データ入力装置１００は重量センサに限らず、他に、圧力センサ、温湿度センサ、変位センサ等の各種センサや、色、輝度を捕らえるビデオカメラ等の映像入力機器であってもよい。
【００１９】
統計処理装置３００を動作するには、データ入力装置１００から出力されるディタル信号を統計処理装置３００に入力して、各データテーブルを更新、参照し、最大値、最頻値を求め、統計処理結果として表示装置２００に出力する。
図５は、本発明の一実施例の階級番号別度数分布テーブルに基づいて作成したヒストグラムを示す。同図は、最古サンプルデータの属する階級番号と異なる階級番号に属する最新サンプルデータについて、最古、最新サンプルデータにおける処理別に最大値、最頻値が変動する全ての場合を示した状態遷移を表している。同図の記載において、“Ｍｏ”は最頻値、“ＭＡ”は最大値、“ＭＢ”は階級番号を表す。
【００２０】
0003564488.tif 000002
ここで、最頻値は、階級番号別度数分布テーブルにおいて度数が最も大きい階級（複数有る場合はいずれかの階級）の階級値とする。
【００２１】
ヒストグラム（ａ），（ｂ），（ｃ），（ｄ）を用いて、最古サンプルデータに関する処理について以下に説明する。
ヒストグラム（ａ）からヒストグラム（ｃ）への遷移は、ヒストグラム（ａ）における最頻度Ｍｏの属する階級番号と同階級番号に属する最古サンプルデータ２０３１が破棄されることにより、ヒストグラム（ａ）において最大値ＭＡを与える階級番号のうち最頻値Ｍｏの属する階級番号以外の階級番号ＭＢの階級値がヒストグラム（ｃ）において最頻値Ｍｏとなる最頻値変動（以下、最頻値のＡ型変動と呼ぶ）を表す（最頻値Ｍｏ＝階級番号ＭＢ）。
【００２２】
ヒストグラム（ｂ）からヒストグラム（ｄ）への遷移は、ヒストグラム（ｂ）において最頻値Ｍｏの属する階級番号と同階級番号に属する最古サンプルデータ２０４１が破棄されることにより、ヒストグラム（ｂ）において最大値ＭＡが１度数減算された値がヒストグラム（ｄ）において最大値ＭＡとなる最大値変動を表す。
【００２３】
サンプリング前の現時刻からｎ−１時刻前までのサンプルデータにおいて、最頻値の属する階級番号Ａ１と最大値Ｂ、及び階級番号Ａ１以外の最大値Ｂを与える階級番号の数Ｃとその階級番号Ｑｉ  （ｉ＝１，２，…，Ｃ）を管理して、サンプリング後の最古サンプルデータの階級番号Ｄと比較し、最頻値、最大値が変動するときのみ、管理する値から最頻値、最大値を求めることができる。
【００２４】
即ち、ヒストグラム（ａ）からヒストグラム（ｃ）への遷移で示す最頻値のＡ型変動に対しては、最大値を与える階級番号の数が２以上（Ｃ≧１）のとき、かつ、階級番号Ａ１と階級番号Ｄが同階級番号であるときに生じ、最古サンプルデータの破棄処理により、最頻値は階級番号Ｑ１  の階級値となる。また、ヒストグラム（ｂ）からヒストグラム（ｄ）への遷移で示す最大値変動に対しては、最大値を与える階級番号の数が１つ（Ｃ＝０）のとき、かつ、階級番号Ａ１と階級番号Ｃが同階級番号であるときに生じ、最古サンプルデータに関する処理により、最大値は最大値Ｂから１減じた値となる。
【００２５】
次に、ヒストグラム（ｅ），（ｆ），（ｇ），（ｈ）を用いて最新サンプルデータに関する処理について以下に示す。
ヒストグラム（ｅ）からヒストグラム（ｇ）への遷移は、ヒストグラム（ｅ）において最大値を与える階級番号のうち、最頻値Ｍｏの属する階級番号以外の階級番号ＭＣと同じ階級番号に属する最新サンプルデータ２０７１が取得されることにより、階級番号ＭＣの度数が１度数加算され、ヒストグラム（ｅ）の最大値ＭＡに１度数加算された値をヒストグラム（ｇ）において最頻値Ｍｏとする最大値、最頻値変動を表す。
【００２６】
ヒストグラム（ｆ）からヒストグラム（ｈ）への遷移はヒストグラム（ｆ）において最頻値の属する階級番号と同階級番号に属する最新サンプルデータ２０８１が取得されることにより、ヒストグラム（ｆ）において最頻値Ｍｏの属する階級番号の度数に１度数加算された値がヒストグラム（ｈ）において最大値ＭＡとする最大値変動を表す。
【００２７】
サンプリング後の最古サンプルデータに関する処理後の１時刻前からｎ−１時刻前までのサンプルデータにおいて、最頻値を与える階級番号Ｄ１と最大値Ｅ、及び階級番号Ｄ１以外の最大値Ｅを与える階級番号の数Ｆとその階級番号Ｒｉ  （ｉ＝１，２，…，Ｆ）を管理し、サンプリング後の最新サンプルデータの階級番号Ｇと比較し、最頻値、最大値が変動するときのみ、管理する値から最頻値、最大値を求めることができる。
【００２８】
即ち、ヒストグラム（ｅ）からヒストグラム（ｇ）への遷移で示す最大値、最頻値変動に対しては、最大値を与える階級番号の数が２以上（Ｆ≧１）のとき、かつ、階級番号Ｇが階級番号Ｒｉ  （ｉ＝１，２，…，Ｆ）のいずれかと同階級番号であるときに生じ、最新サンプルデータに関する処理により、最頻値は、階級番号Ｇの階級値、最大値は最大値Ｅに１加えた値となる。また、ヒストグラム（ｆ）からヒストグラム（ｈ）への遷移で示す最大値変動に対しては、最大値を与える階級番号の数が１つ（Ｆ＝０）のとき、かつ、階級番号Ｄ１階級番号Ｇが同階級番号であるときに生じ、最新サンプルデータに関する処理により、最大値は最大値Ｅに１加えた値となる。
【００２９】
つまり、最頻値、最大値の求解処理は、サンプリング後の最古のサンプルデータを破棄することに伴う最頻値、最大値の変動に関する条件と最新サンプルデータを取得することによる最頻値、最大値の変動に関する条件について監視し、最頻値、最大値の変動する条件時のみ、最頻値、最大値の求解処理を実施すればよい。
【００３０】
図６は、本発明の一実施例の統計処理装置の構成を示す。
同図に示す統計処理装置３００は、信号入力端子３０１、サンプルデータ入力回路３０２、リングバッファ回路３０３、階級番号取得回路３０４、最大値・最頻値求解回路３０５、記憶回路３０６、出力端子３０７、リングバッファテーブル３０３１、サンプルデータ階級番号変換テーブル３０４１、階級番号別度数分布テーブル３０６１、度数別階級番号分布テーブル３０６２、度数別先頭アドレステーブル３０６３、度数別階級番号有効数テーブル３０６４、最大値メモリ３０６５、最頻値階級番号メモリ３０６６、階級番号−階級値変換テーブル３０６７より構成される。
【００３１】
ここで、リングバッファ回路３０３により参照されるリングバッファテーブル３０３１は、各サンプリング対象毎に現時刻からｎ時刻前までのサンプルデータを入力順に管理するテーブルであり、初期化は全領域に０を設定する。
階級番号取得回路３０４により参照されるサンプルデータ−階級番号変換テーブル３０４１は、サンプルデータの属する階級番号を与えるため、サンプルデータとサンプルデータの属する階級番号を対応づけるテーブルであり、初期化はＡ／Ｄ変換器の分解能と階級間隔に応じて、予め取り得るサンプルデータに対応した階級番号を設定する。例えば、８ビットの分解能を持つＡ／Ｄ変換器では、取り得るサンプルデータは、０〜２５５の２５６通りであり、階級間隔を３とすれば、サンプルデータが０から２のときは、階級番号０、サンプルデータが３から５のときは、階級番号「１」というように、階級番号８５まで８６個の階級について順次設定すればよい。
【００３２】
記憶回路３０６に接続される階級番号別度数分布テーブル３０６１は、各サンプリング対象毎のリングバッファテーブル３０３１で管理する現時刻からｎ−１時刻前までのサンプルデータにおいて、各サンプリング対象毎に、階級番号毎の度数を管理するテーブルであり、初期化は各サンプリング毎の階級番号別度数分布テーブル３０６１において階級番号０の度数をｎ、残りの階級番号の度数を度数０と設定する。
【００３３】
記憶回路３０６に接続される度数別階級番号テーブル３０６２は、各サンプリング対象毎の階級番号別度数分布テーブル３０６１において、度数毎に階級番号を並べ替えたテーブルである。リングバッファテーブル３０３１のサンプリング対象当たりのデータ数をｎとするとき、階級番号別度数分布テーブル３０６１において、同じ度数を持つ複数の階級番号が発生するのは、
１≦（階級番号別度数分布テーブル３０６１の度数）≦ｎ／２
の時である。しかし、階級番号別度数分布テーブル３０６１における階級番号の度数が１の場合は、最大値、最頻値を求める上で管理不要であるため、度数別階級番号テーブル３０６２が管理するのは、
２≦（階級番号別度数分布テーブル３０６１の度数）≦ｎ／２
の場合であり、１サンプリング対象における度数別階級番号テーブル３０６２の容量は、
【００３４】
【数１】
0003564488.tif 000003
【００３５】
で表すことができる。
また、度数別階級番号テーブル３０６２は、最新サンプルデータの属する階級番号と最古サンプルデータの属する階級番号が異なるとき、度数別先頭アドアレステーブル３０６３と度数別階級番号有効数テーブル３０６４を用いて更新され、初期化は必要としない。
【００３６】
記憶回路３０６に接続される度数別先頭アドレステーブル３０６３は、度数別階級番号テーブル３０６２が度数毎に管理する階級番号において、先頭の階級番号が保存されている位置を度数毎に示すアドレステーブルである。なお、位置は、各サンプリング対象毎の度数別階級番号テーブル３０６２の先頭からデータ数で表す。度数別階級番号テーブル３０６２を最小限のメモリ容量で実現するため、度数別階級番号テーブル３０６２が度数毎に管理する階級番号の最大数が度数により異なる（例えば、度数２によって管理する階級番号の数は最大ｎ／２個であり、また、度数（ｎ／２）によって管理する階級番号の数は最大２個となる）ことを利用して度数ｊ、（但し、ｊは、２≦ｊ≦ｎ／２）の先頭アドレスＡＤｊ  は、
ＡＤ２  ＝１                                          （１）
ＡＤｊ  ＝ＡＤ（ｊ−１）  ＋ｉｎｔ（ｎ／（ｊ−１）），
（３≦ｊ≦ｎ／２）  （２）
と初期化される。
【００３７】
記憶回路３０６に接続される度数別階級番号有効数テーブル３０６３は、各サンプリング対象毎の度数別階級番号テーブル３０６２において、度数毎に保存される階級番号の数を度数別に管理するテーブルであり、初期化は全領域に「０」を設定する。
【００３８】
記憶回路３０６に接続される最大値メモリ３０６５は、各サンプリング対象毎の階級番号別度数分布テーブル３０６１における度数の最大値を各サンプリング対象毎に管理するメモリであり、初期化は全領域に度数ｎを設定する。
記憶回路３０６に接続される最頻値階級番号メモリ３０６６は、各サンプリング対象毎の階級番号別度数分布テーブル３０６１において、度数が最大となる階級番号の１つを各サンプリング対象毎に管理するメモリであり、初期化は全領域に階級番号「０」を設定する。
【００３９】
階級番号−階級値変換テーブル３０６７は、階級番号毎の階級値を管理するテーブルである。例えば、階級の数８６個、級間隔３（サンプルデータ−階級番号変換テーブル３０４１の設定条件）とすれば、階級番号０は階級値１とし、順次、階級番号８５まで、階級番号と階級値を対応付ける。
【００４０】
以下に上記の構成の統計処理装置３００の動作を説明する。
図７は、本発明の一実施例の統計処理装置の一連の概要動作のフローチャートである。
ステップ１０１）  データ入力装置１００から信号入力端子３０１を介してサンプルデータ入力回路３０２に、サンプリング対象毎のサンプルデータが入力される。
【００４１】
ステップ１０２）  サンプルデータ入力回路３０２は、サンプリング対象データに対応する空間の座標を与える。即ち、サンプルデータが入力されると、そのときのサンプリング対象の固有番号を付し（以下、サンプリング対象データの固有番号を入力点とよび、サンプリング対象数をｍとする）、サンプルデータ列ｆ（ｘ）（但し、ｘ＝０，１，２，…，ｍ−１）としてリングバッファ回路３０３０に出力する。
【００４２】
ステップ１０３）  リングバッファ回路３０３は、サンプリング毎に各入力点毎のリングバッファテーブル３０３１を更新する。
ステップ１０４）  リングバッファ回路３０３は、各入力点毎のリングバッファテーブル３０３１を用いて入力データを１つ取得すると、ｎ時刻前のサンプルデータ（ｎ時刻までは、各時刻に相当する初期化時のサンプルデータ）を１つ出力させることにより、サンプリング毎に、現時刻のサンプルデータ列ｆ（ｘ）を最新のサンプルデータ列ｆｎｅｗ  （ｘ）として階級番号取得回路３０４に出力する。また、ｎ時刻前のサンプルデータ列ｆ（ｘ）を最古のサンプルデータ列ｆｏｌｄ  （ｘ）として、階級番号取得回路３０４に出力する。但し、上記のｘは、ｘ＝０，１，２，…，ｍ−１とする。
【００４３】
ステップ１０５）  階級番号取得回路３０４は、サンプルデータ−階級番号変換テーブル３０４１を用いて、最新のサンプルデータ列ｆｎｅｗ  （ｘ）の属する階級番号列Ｆｎｅｗ  （ｘ）を最大値・最頻値求解回路３０５に出力する（ただし、ｘ＝０，１，２，…，ｍ−１）。
【００４４】
ステップ１０６）  階級番号取得回路３０４は、サンプルデータ−階級番号変換テーブル３０４１を用いて、最古のサンプルデータ列ｆｏｌｄ  （ｘ）の属する階級番号列Ｆｏｌｄ  （ｘ）を最大値・最頻値求解回路３０５に出力する（ただし、ｘ＝０，１，２，…，ｍ−１）。
【００４５】
ステップ１０７）  最大値・最頻値求解回路３０５は、入力された階級番号列Ｆｎｅｗ  （ｘ）において、各入力点ｘ，（ｘ＝０，１，２，…，ｍ−１）毎に階級番号列Ｆｏｌｄ  （ｘ）と階級番号列Ｆｎｅｗ  （ｘ）の階級番号を比較する。
ステップ１０８）  Ｆｏｌｄ  （ｘ）≠Ｆｎｅｗ  （ｘ）であれば、その入力点ｘについてのみ、最大値、最頻値求解の処理を行い、最大値メモリ３０６５、最頻値階級番号メモリ３０６６を更新する。
【００４６】
ステップ１０９）  最頻値階級番号メモリ３０６６から求めた各入力点毎の最大値を与える階級番号と階級番号−階級値変換テーブル３０６７を用いて、各入力点毎の最頻値Ｍｏ（ｘ）（ｘ＝０，１，２，…，ｍ−１）を求める。
ステップ１１０）  最大値・最頻値求解回路３０５は、最大値メモリ３０６５から各入力点毎の最大値ＭＡ（ｘ）（ｘ＝０，１，２，…，ｍ−１）を求める。
【００４７】
ステップ１１１）  最頻値Ｍｏ（ｘ）（ｘ＝０，１，２，…，ｍ−１）と最大値ＭＡ（ｘ）（ｘ＝０，１，２，…，ｍ−１）を出力端子１０３に出力する。
次に、最大値・最頻値求解の処理を説明する。
図８は、本発明の一実施例の最大値・最頻値求解処理のフローチャートである。同図において、使用されている以下の符号の意味を示す。
Ｆｏｌｄ  ：        入力点ｘの最古サンプルデータの階級番号；
ｄｏｓ１：      入力点ｘの階級番号別度数分布テーブル３０６１から求めた階
級番号Ｆｏｌｄ  （ｘ）の度数；
ｎ：            リングバッファテーブルのデータ数；
Ｆｎｅｗ  （ｘ）：  入力点ｘの最新サンプルデータの階級番号；
ｄｏｓ２：      入力点ｘの階級番号別度数分布テーブル３０６１から求めた階
級番号Ｆｎｅｗ  （ｘ）の度数：
ｄｏｓ３：      入力点ｘの階級番号別度数分布テーブル３０６１から求めた階
級番号Ｆｎｅｗ  （ｘ）の度数：
ｍａ１：        最大値メモリ３０６５から求めた入力点ｘの最大値：
ステップ２０１）  まず、Ｆｏｌｄ  （ｘ）とＦｎｅｗ  （ｘ）が一致するかを調べ、Ｆｏｌｄ  （ｘ）＝Ｆｎｅｗ  （ｘ）であれば、ステップ２０２に移行し、Ｆｏｌｄ  （ｘ）≠Ｆｎｅｗ  （ｘ）であれば、当処理を終了する。
【００４８】
ステップ２０２）  最古サンプルデータの破棄によるテーブル更新処理を行う。ｄｏｓ１が、
２≦ｄｏｓ１≦ｎ／２
の関係にあれば、ステップ２０３に移行し、この関係が成立しない場合には、ステップ２０４に移行する。
【００４９】
ステップ２０３）  図９に示すプロセス１の処理を行う。
ステップ２０４）  図１０に示すプロセス２の処理を行う。
ステップ２０５）  最新サンプルデータの取得によるテーブルの更新処理を行う。入力点ｘの階級番号別度数分布テーブル３０６１から階級番号Ｆｎｅｗ  （ｘ）の度数ｄｏｓ２を求め、２≦（度数ｄｏｓ２）≦（ｎ／２）が成り立つ場合には、ステップ２０６に移行し、この関係が成立しない場合にはステップ２０７に移行する。
【００５０】
ステップ２０６）  図１１に示すプロセス３の処理を行う。
ステップ２０７）  図１２に示すプロセス４の処理を行う。
ステップ２０８）  入力点の階級番号別度数分布テーブル３０６１から求めた階級番号Ｆｎｅｗ  （ｘ）の度数ｄｏｓ３と最大値メモリ３０６５から求めた入力点ｘの最大値ｍａ１についてｍａ１＜ｄｏｓ３が成り立つとき、ステップ２０９に移行し、ｍ１がｄｏｓ３以上（ｍａ１≧ｄｏｓ３）であれば、当該処理を終了する。
【００５１】
ステップ２０９）  最大値メモリ３０６５の入力点ｘの最大値をｄｏｓ３とする。
ステップ２１０）  最頻値階級番号メモリ３０６６の入力点ｘの最頻値を与える階級番号をＦｎｅｗ  （ｘ）とする。
【００５２】
図９は、本発明の一実施例のプロセス１のフローチャートである。
同図に示すフローチャートは、図８のステップ２０２に対応する処理であり、階級番号Ｆｏｌｄ  （ｘ）をｄｏｓ１の管理から削除する処理と、階級番号Ｆｏｌｄ  （ｘ）を（ｄｏｓ１−１）の管理に追加する処理と、最頻値の階級番号を決定する処理、最大値メモリ３０６５における入力点ｘの最大値からデクリメントする処理より構成される。
【００５３】
同図において、使用されている以下の符号の意味を示す。
ｎｏ１：  入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ１で管理す
る階級番号の数；
ｍｏ１：  最頻値階級番号メモリ３０６６から求めた入力点ｘの最頻値の属する
階級番号：
最初に入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ１で管理する階級番号Ｆｏｌｄ  （ｘ）を削除する処理を行う。
【００５４】
ステップ３０１）  この削除処理では、入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ１で管理する階級番号の数ｎｏ１を入力点ｘの度数別階級番号有効数テーブル３０６４から求め、階級番号の数ｎｏ１が「２以上」の場合と「１」の場合とに分けて処理を行う。従って、ｎｏ１が２以上である場合（ｎｏ１≧２）には、ステップ３０２に移行し、ｎｏ１が１の場合（ｎｏ１＜２）にはステップ３０３に移行する。
【００５５】
ステップ３０２）  階級番号の数ｎｏ１が２以上（ｎｏ１≧２）の場合では、入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ１で管理する階級番号において、先頭の階級番号が保存されている位置をアドレスｊｓ１、後尾の階級番号が保存されているアドレスｊｅ１、階級番号Ｆｏｌｄ  （ｘ）が保存されている位置をアドレスｊｒ１とすると、度数別先頭アドレステーブル３０６３から求めたアドレスｊｓ１と階級番号の数ｎｏ１を用いてアドレスｊｅ１は、
ｊｅ１＝ｊｓ１＋ｎｏ１−１
となるため、アドレスｊｒ１は、アドレスｊｓ１からアドレスｊｅ１までに保存される階級番号と階級番号Ｆｏｌｄ  （ｘ）が同階級番号であるかを比較することで求められる。そこで、アドレスｊｒ１にアドレスｊｅ１の階級番号をセットすることで実領域の削除を行う。なお、アドレスは、データ個数を用いてテーブルの先頭からの位置を表す。つまり、アドレスｊｓ１は、度数別階級番号テーブル３０６２の先頭からｊｓ１番目の位置を表している。以下のアドレス表記においても同様とする。
【００５６】
ステップ３０３）  階級番号の数ｎｏ１＝１のとき、入力点ｘの度数別階級番号有効数テーブル３０６４における度数ｄｏｓ１の有効数を１減じることで管理領域の削除を行う。
【００５７】
ステップ３０４）  次に、階級番号Ｆｏｌｄ  （ｘ）を（ｄｏｓ１−１）の管理に追加する処理を行う。度数別階級番号テーブル３０６２が管理する階級番号の度数は２以上であるため、度数（ｄｏｓ１−１≧２）が成り立つかを判断し、成り立つ場合には、ステップ３０５に移行し、成り立たない場合には、ステップ３０７に移行する。
【００５８】
ステップ３０５）  度数（ｄｏｓ１−１）≧２が成り立つとき、入力点ｘの度数別階級番号テーブル３０６２が度数（ｄｏｓ１−１）で管理する階級番号の後尾に階級番号Ｆｏｌｄ  （ｘ）を追加する追加処理を行う。
ステップ３０６）  ステップ３０５の追加処理では、入力点ｘの度数別階級番号テーブル３０６２が度数（ｄｏｓ１−１）で管理する階級番号において、階級番号の数をｎｏ２、先頭の階級番号が保存されている位置をアドレスｊｓ２、後尾の階級番号が保存されている位置アドレスをｊｅ２とすると、度数別先頭アドレステーブル３０６３から求めたアドレスｊｓ２と入力点ｘの度数別階級番号有効数テーブル３０６４から求めたｎｏ２を用いてアドレスｊｅ２は、
ｊｅ２＝ｊｓ２＋ｎｏ２−１
となるため、アドレス（ｊｅ２＋１）に階級番号Ｆｏｌｄ  （ｘ）をセットすることで実領域の追加と、入力点Ｘの度数別階級番号有効数テーブル３０６４における度数（ｄｏｓ１−１）の有効数に１加えることで管理領域の追加を行う。
【００５９】
ステップ３０７）  そして、入力点ｘの階級番号別度数分布テーブル３０６１における階級番号Ｆｏｌｄ  （ｘ）の度数から１減じる。
ステップ３０８）  最頻値階級番号メモリ３０６６から求めた入力点ｘの最頻値の属する階級番号ＭＯ１において、階級番号ＭＯ１＝階級番号Ｆｏｌｄ  （ｘ）かつ、ｎｏ１＞１が成り立つかを判断し、成り立つ時は、ステップ３０９に移行し、成り立たない場合には、ステップ３１０に移行する。
【００６０】
ステップ３０９）  階級番号ＭＯ１＝階級番号Ｆｏｌｄ  （ｘ）かつ、ｎｏ１＞１が成り立つとき、「入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ１で管理する階級番号のうちで先頭に位置する階級番号」つまり、「入力点ｘの度数別階級番号テーブル３０６２において、アドレス（ｊｓ１＋ｎｏ１−１）に保存されている階級番号」を「最頻値階級番号メモリ３０６６で管理する入力点ｘの最頻値の階級番号」とする。
【００６１】
ステップ３１０）  階級番号ＭＯ１＝Ｆｏｌｄ  （ｘ）が成り立つ時には、ステップ３１１に移行し、成立しない場合には当該プロセス１の処理を終了する。
ステップ３１１）  階級番号ＭＯ１＝Ｆｏｌｄ  （ｘ）が成り立つ時には、最大値メモリ３０６５における入力点Ｘの最大値を１減じる。
【００６２】
次に、図８のステップ２０４のプロセス２の処理を説明する。
図１０は、本発明の一実施例のプロセス２のフローチャートである。同図において使用されている以下の符号の意味を示す。
ＭＯ２：  最頻値階級番号メモリ３０６６から求めた入力点ｘの最頻値の属する階級番号；
プロセス２では、度数（ｄｏｓ１）＝（ｎ／２）＋１における追加処理と、入力点ｘの最頻値の属する階級番号ＭＯ２＝階級番号Ｆｏｌｄ  （ｘ）における最大値処理を行う。
【００６３】
ステップ４０１）  入力点ｘの階級番号別度数分布テーブル３０６１における階級番号Ｆｏｌｄ  （ｘ）の度数を１減じる。
ステップ４０２）  度数別階級番号テーブル３０６２が度数毎に管理する階級番号の数は（ｎ／２）以下であるため、度数（ｄｏｓ１−１）＝（ｎ／２）であるかを判断し、度数（ｄｏｓ１−１）＝（ｎ／２）のときは、ステップ４０３に移行し、（ｄｏｓ１−１）≠（ｎ／２）のときは、ステップ４０４に移行する。
【００６４】
ステップ４０３）  入力点ｘの度数別階級番号テーブル３０６２が度数ｎ／２で管理する階級番号の後尾に階級番号Ｆｏｌｄ  （ｘ）を追加する追加処理を行う。ステップ４０４）  ステップ４０３の追加処理では、入力点ｘの度数別階級番号テーブル３０６２が度数（ｄｏｓ−１）で管理する階級番号において、階級番号の数をｎｏ３、先頭の階級番号が保存される位置をアドレスｊｓ３、後尾の階級番号が保存される位置をアドレスｊｅ３とすると、度数別先頭アドレステーブル３０６３から求めたアドレスｊｓ３と入力点ｘの度数別階級番号有効数テーブル３０６４から求めたｎｏ３を用いてアドレスｊｅ３は、
ｊｅ３＝ｊｓ３＋ｎｏ３−１
となるため、アドレス（ｊｅ３＋１）に階級番号Ｆｏｌｄ  （ｘ）をセットすることで、実領域の追加と、入力点ｘの度数別階級番号有効数テーブル３０６４における度数（ｄｏｓ−１）の有効数に１加えることで管理領域の追加を行う。
【００６５】
ステップ４０５）  最大値の処理として、最頻値階級番号メモリ３０６６から求めた入力点ｘの最頻値の属する階級番号ＭＯ２＝階級番号Ｆｏｌｄ  （ｘ）であれば、ステップ４０６に移行し、ＭＯ２≠Ｆｏｌｄ  （ｘ）の場合には当該プロセス２の処理を終了する。
【００６６】
ステップ４０６）  階級番号ＭＯ２＝階級番号Ｆｏｌｄ  （ｘ）のとき、最大値メモリ３０６５における入力点ｘで管理する最大値を１減じる。
次に、図８のステップ２０５〜ステップ２０７における最新サンプルデータの取得によるテーブル更新処理を説明する。
【００６７】
入力点ｘの階級番号別度数分布テーブル３０６１から階級番号Ｆｎｅｗ  （ｘ）の度数ｄｏｓ２を求め、２≦（度数ｄｏｓ２）≦（ｎ／２）が成り立つ場合には、プロセス３を行い、成立しない場合には、プロセス４を実行する。
図１１は、本発明の一実施例のプロセス３のフローチャートである。同図において、使用される符号を説明する。
【００６８】
ｎｏ４：  入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ２で管理する階級番号の数；
最初に入力点ｘの度数階級番号テーブル３０６２が度数ｄｏｓ２で管理する階級番号Ｆｎｅｗ  （ｘ）を削除する削除処理を行う。この削除処理では、入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ２で管理する階級番号の数ｎｏ４を入力点ｘの度数階級番号有効数テーブル３０６４から求め、階級番号ｎｏ４が２以上の場合と１の場合に分けて処理を行う。
【００６９】
ステップ５０１）  階級番号の数ｎｏ４が２以上（ｎｏ４≧２）であるかを判断し、ｎｏ４≧２の場合には、ステップ５０２に移行する。１の場合には、ステップ５０３に移行する。
ステップ５０２）  階級番号の数がｎｏ４≧２の場合には、入力点ｘの度数別階級番号テーブル３０６２が度数ｄｏｓ２で管理する階級番号において、先頭の階級番号が保存されている位置をアドレスｊｓ４、後尾の階級番号が保存されている位置をアドレスｊｅ４、階級番号Ｆｎｅｗ  （ｘ）が保存されている位置をアドレスｊｒ４とすると、度数別先頭アドレステーブル３０６３から求めたアドレスｊｓ４と階級番号の数ｎｏ４を用いてアドレスｊｅ４は、
ｊｅ４＝ｊｓ４＋ｎｏ４−１
となるため、アドレスｊｒ４は、アドレスｊｓ４からアドレスｊｅ４までに保存される階級番号と階級番号Ｆｎｅｗ  （ｘ）が同階級番号であるかを比較することで求められる。そこで、アドレスｊｒ４にアドレスｊｅ４の階級番号をセットすることで実領域の削除と、入力点ｘの度数別階級番号有効数テーブル３０６４における度数ｄｏｓ２の有効数を１減じることで管理領域の削除を行い、ステップ５０３に移行する。
【００７０】
ステップ５０３）  一方、階級番号の数ｎｏ４が１（ｎｏ４＝１）の場合には、入力点ｘの度数別階級番号テーブル３０６２における削除処理を行わず、入力点ｘの度数別階級番号テーブル３０６２における削除処理を行わず、入力点ｘの度数別階級番号有効数テーブル３０６４における度数ｄｏｓ２の有効数を１減じる。
【００７１】
ステップ５０４）  度数別階級番号テーブル３０６２が管理する階級番号の度数は、（ｎ／２）以下であるため、度数（ｄｏｓ２＋１）≦（ｎ／２）が成り立つかを判断し、成り立つ場合には、ステップ５０５に移行し、成立しない場合には、ステップ５０７に移行する。
【００７２】
ステップ５０５）  度数（ｄｏｓ２＋１）≦（ｎ／２）成り立つとき、入力点ｘの度数別階級番号テーブル３０６２が度数（ｄｏｓ２＋１）で管理する階級番号の後尾に階級番号Ｆｎｅｗ  （ｘ）を追加する。
ステップ５０６）  ステップ５０５における追加処理では、入力点ｘの度数別階級番号テーブル３０６２が度数（ｄｏｓ２＋１）で管理する階級番号において、階級番号の数をｎｏ５、先頭の階級番号が保存されている位置アドレスｊｓ５、後尾の階級番号が保存されている位置をアドレスｊｅ５とすると、度数先頭アドレステーブル３０６３から求めたアドレスｊｓ５と入力点ｘの度数別階級番号有効数テーブル３０６４から求めたｎｏ５を用いてｊｅ５は、
ｊｅ５＝ｊｓ５＋ｎｏ５−１
となるため、アドレス（ｊｅ５＋１）に階級番号Ｆｎｅｗ  （ｘ）をセットすることで、実領域の追加と、入力点ｘの度数別階級番号有効数テーブル３０６４における度数（ｄｏｓ２＋１）の有効数を１加えることで、管理領域の追加を行う。
【００７３】
ステップ５０７）  入力点ｘの階級番号別度数分布テーブル３０６１における階級番号Ｆｎｅｗ  （ｘ）の度数を１加える。
次に、プロセス４について説明する。
図１２は、本発明の一実施例のプロセス４のフローチャートである。
【００７４】
プロセス４では、度数（ｄｏｓ１）＝３における追加処理を行う。
ステップ６０１）  入力点ｘの階級番号別度数分布テーブル３０６１における階級番号Ｆｎｅｗ  （ｘ）の度数に１を加える。
ステップ６０２）  度数階級番号テーブル３０６２が度数毎に管理する階級番号の数であるため、度数（ｄｏｓ２＋１）＝２が成り立つかを判断し、成り立つ場合には、ステップ６０３に移行し、成り立たない場合には、当該プロセスの処理を終了する。
【００７５】
ステップ６０３）  度数（ｄｏｓ２＋１）＝２が成り立つとき、入力点ｘの度数別階級番号テーブル３０６２が度数２で管理する階級番号の後尾に階級番号Ｆｎｅｗ  （ｘ）を追加する追加処理を行う。
ステップ６０４）  ステップ６０３の追加処理では、入力点ｘの度数別階級番号テーブル３０６２が度数（ｄｏｓ２＋１）で管理する階級番号において、階級番号の数をｎｏ６、先頭の階級番号が保存されている位置をアドレスｊｓ６、後尾の階級番号が保存されている位置をアドレスｊｅ６とすると、度数別先頭アドレステーブル３０６３から求めたアドレスｊｓ６と、入力点ｘの度数別階級番号有効数テーブル３０６４から求めたｎｏ６を用いてアドレスｊｅ６は、
ｊｅ６＝ｊｓ６＋ｎｏ６−１
となるため、アドレス（ｊｅ６＋１）に階級番号Ｆｎｅｗ  （ｘ）をセットすることで、実領域の追加と、入力点ｘの度数別階級番号有効数テーブル３０６４における（ｄｏｓ２＋１）の有効数に１加えることで管理領域の追加を行う。
【００７６】
記憶回路３０６は、最大値・最頻値求解回路３０５におけるテーブルの更新処理に伴いテーブルからの読み出し、テーブルへの書き込みを管理・実行する。
次に、テーブルの管理方法について説明する。
図１３は、本発明の一実施例の階級番号別度数分布テーブルを説明するための図、図１４は、本発明の一実施例の度数別階級番号分布テーブルを説明するための図、図１５は、本発明の一実施例の度数別階級番号分布テーブルのメモリ構造イメージ図、図１６は、本発明の一実施例の度数別先頭アドレステーブルを説明するための図、図１７は、本発明の一実施例の度数別階級番号有効数テーブルを説明するための図である。
【００７７】
図１３において、縦軸４０１１は、入力点を示し、横軸４０１２は、階級番号を示す。同図において、度数４０１３は、「５」であることを示している。
階級番号別度数分布テーブル３０６１は、階級番号毎の度数を管理するテーブルである。以下、説明の便宜上、例えば、サンプルデータの範囲は、０〜２５５、サンプルデータ数は６０個、入力点の数はｍ個、級間隔は３とし、入力点４０１１の「１」の階級番号別度数分布テーブル３０６１における階級番号別度数は、階級番号４０１２「７１」の度数４０１３「５」、階級番号４０１２「７０」の度数４０１３「１５」、階級番号４０１２「８２」の度数４０１３「２０」、階級番号４０１２「８５」、度数４０１３「２０」、その他の階級番号４０１２の度数４０１３を全て「０」とする。
【００７８】
図１７において、縦軸４０１１は、入力点を示し、横軸４０１３は階級番号を示す。図１７の度数別階級番号有効数テーブル３０６４は、度数別先頭アドレステーブル３０６３における階級番号の有効数を度数毎に管理するテーブルである。入力点１の度数別階級番号有効数テーブル３０６４における階級番号別有効数は、入力点１の階級番号別度数分布テーブル３０６１における階級番号４０「７１」の度数４０１３「５」に対応して、「度数５の有効数が１」、同様に、度数「１５」の有効数が「１」、度数「２０」の有効数が「２」となる。入力点１の階級番号別度数分布テーブル３０６１における度数「０」となる階級番号の有効数は８２個であるが、管理対象の階級番号の有効数は、
２≦（階級番号別度数分布テーブル３０６１の度数）≦ｎ／２
であるから、度数別階級番号有効数テーブル３０６４では、度数「０」における階級番号の有効数を管理しない。
【００７９】
図１６において、度数別先頭アドレス管理テーブル３０６３は、度数別階級番号分布テーブル３０６２が、度数毎に管理する階級番号において先頭の階級番号が保存される位置を表すテーブルである。
度数別階級番号分布テーブル３０６２において、度数「５」を与える階級番号が保存されている先頭アドレスは「６６」である。
【００８０】
ここで、度数ｊ，（２≦ｊ≦ｎ／２）の先頭アドレスＡＤｊ  は、式（１）、（２）の関係式から
ＡＤ２  ＝１
ＡＤｊ  ＝ＡＤ（ｊ−１）  ＋ｉｎｔ（ｎ／（ｊ−１）），（３≦ｊ≦ｎ／２）
である。この関係式は、階級番号別度数分布テーブル３０６１において同じ度数ではあるが、異なる階級番号が発生する頻度に基づいている。
【００８１】
図１４の度数別階級番号分布テーブル３０６２は、階級番号別度数分布テーブル３０６１における階級番号を度数呉に管理するテーブルである。入力点「１」の階級番号別度数分布テーブル３０６１に対応して、入力点「１」の度数別階級番号分布テーブル３０６２は、アドレス６６に階級番号７１、アドレス１３３に階級番号「７０」、アドレス１４９に階級番号「８４」、アドレス１５０に階級番号「８５」を保持する。
【００８２】
まず、度数別階級番号分布テーブル３０６２の管理方法を説明するため、度数別階級番号分布テーブル３０６２における度数別階級番号の参照例を示す。例えば、入力点「１」の階級番号別度数分布テーブル３０６１における階級番号「７１」の度数「５」に対応する、入力点「１」の度数別階級番号分布テーブル３０６２における度数「５」の階級番号は度数別先頭アドレステーブル３０６３における度数「５」の先頭アドレスが「６６」であるから、アドレス「６６」以降に順次保持される。また、入力点「１」の度数別階級番号有効数テーブル３０６４における度数「５」の階級番号有効数が「１」であるから、入力点「１」の度数別階級番号分布テーブル４０４における度数５の階級番号は１個保存する。つまり、入力点「１」の度数別階級番号有効数テーブル３０６４における度数「５」の階級番号は、アドレス「６６」に１個保存し、アドレス「６６」から階級番号「７１」と参照できる。
【００８３】
次に、図１４における度数別階級番号分布テーブル３０６２の管理方法について説明する。
図１５における度数別階級番号分布テーブル３０６２のメモリに示すように、度数別階級番号分布テーブル３０６２の管理方法は、度数別階級番号有効数テーブル３０６４の有効数と度数別先頭アドレステーブル３０６３の先頭アドレスを用いることで、度数毎に必要最小限のメモリ容量で度数別階級番号分布テーブル３０６２を実現することができる。また、度数別階級番号分布テーブル３０６２を直接参照することなく保存する階級番号の数と後尾に位置する階級番号のアドレスが求められるため、処理時間の短縮が図ることができる。
【００８４】
本発明は、上記の実施例に限定されることなく、特許請求の範囲内で種々変更・応用が可能である。
【００８５】
【発明の効果】
上述のように、本発明によれば、常にサンプリング毎に得られる最古・最新サンプルデータが属する階級番号を求めるとき、サンプルデータ−階級番号変換テーブルを用いることと、階級番号別度数分布テーブルにおける度数の最大値、最頻値が変更する条件のときのみ最大値、最頻値を求めることで、処理時間を短縮することが可能である。
【００８６】
また、時々刻々変化するシステムに設定したセンサからのサンプルデータを本発明の統計処理装置に入力すると、最新サンプルデータを含む過去ｎ時間のサンプルデータから連続的に最大値、最頻値の実時間求解が可能であるから、映像に利用すれば、画像処理の特徴量抽出、製造ラインに設置すれば、生産物の実時間の自動検査化に利用できる利点がある。
【００８７】
さらに、最大値、最頻値の変動条件に応じた領域の管理構成をとることにより使用メモリを最小化できる。
【図面の簡単な説明】
【図１】本発明の原理構成図である。
【図２】本発明の原理を説明するためのフローチャートである。
【図３】本発明の一実施例の統計装置を適用したシステム構成図である。
【図４】本発明の一実施例の統計装置の使用例を示す図である。
【図５】本発明の一実施例の階級番号別度数分布テーブルに基づいて作成したヒストグラムを示す図である。
【図６】本発明の一実施例の統計処理装置の構成図である。
【図７】本発明の一実施例の統計処理装置の一連の概要動作のフローチャートである。
【図８】本発明の一実施例の最大値、最頻値求解処理のフローチャートである。
【図９】本発明の一実施例のプロセス１のフローチャートである。
【図１０】本発明の一実施例のプロセス２のフローチャートである。
【図１１】本発明の一実施例のプロセス３のフローチャートである。
【図１２】本発明の一実施例のプロセス４のフローチャートである。
【図１３】本発明の一実施例の階層番号別度数分布テーブルを説明するための図である。
【図１４】本発明の一実施例の度数別階級番号分布テーブルを説明するための図である。
【図１５】本発明の一実施例の度数別階級番号分布テーブルのメモリ構造イメージ図である。
【図１６】本発明の一実施例の度数別先頭アドレステーブルを説明するための図である。
【図１７】本発明の一実施例の度数別階級番号有効数テーブルを説明するための図である。
【図１８】従来の統計処理装置の構成図である。
【図１９】従来の統計処理のフローチャートである。
【符号の説明】
１０  ベルトコンベア
２０  生産物
３０  重量測定エリア
１００  データ入力装置
１１０  重量センサ
１２０  Ａ／Ｄ変換器
２００  表示装置
３００  統計処理装置
３０１  信号入力端子
３０２  サンプルデータ入力回路、サンプルデータ入力手段
３０３  リングバッファ回路、蓄積手段
３０４  階級番号取得回路、階級番号付与手段
３０５  最大値・最頻値求解回路、最大値・最頻値求解手段
３０６  記憶回路
２０３１、２０４１  最古サンプルデータ
２０７１、２０８１  最新サンプルデータ
３０３１  リングバッファテーブル
３０４１  サンプルデータ階層番号変換テーブル
３０５１  更新手段
３０６１  階級番号別度数分布テーブル
３０６２  度数別階級番号分布テーブル
３０６３  度数別先頭アドレステーブル
３０６４  度数別階級番号有効数テーブル
３０６５  最大値メモリ
３０６６  最頻値階級番号メモリ
３０６７  階級番号−階級値変換テーブル
４０１１  入力点
４０１２  階級番号
４０１３  度数
４０２１  有効数
４０３１  先頭アドレス
【図１】
0003564488.tif 000004
【図２】
0003564488.tif 000005
【図３】
0003564488.tif 000006
【図４】
0003564488.tif 000007
【図５】
0003564488.tif 000008
【図６】
0003564488.tif 000009
【図７】
0003564488.tif 000010
【図８】
0003564488.tif 000011
【図９】
0003564488.tif 000012
【図１０】
0003564488.tif 000013
【図１１】
0003564488.tif 000014
【図１２】
0003564488.tif 000015
【図１３】
0003564488.tif 000016
【図１４】
0003564488.tif 000017
【図１５】
0003564488.tif 000018
【図１６】
0003564488.tif 000019
【図１７】
0003564488.tif 000020
【図１８】
0003564488.tif 000021
【図１９】
0003564488.tif 000022