"抽出対象周波数の1周期分に相当する所定時間毎に、その所定時間に入力される入力信号を信号素片として格納する入力バッファと、前記所定時間毎に前記入力バッファから供給される前記信号素片の、前記所定時間内における各時点の信号値をそれぞれ累積加算した累積加算結果を保持するためのレジスタを有し、前記所定時間が経過する毎に、前記レジスタに保持された前記累積加算結果と、新たに前記入力バッファから供給された信号素片との、前記所定時間内における同一時点の信号値をそれぞれ加算して、新たな累積加算結果を求め、その新たな累積加算結果を前記レジスタに保持させる演算部と、を備えることを特徴とする周波数解析装置。前記演算部は、前記累積加算した回数で前記累積加算結果を除算することを特徴とする請求項1に記載の周波数解析装置。前記入力バッファ及び前記演算部は、複数の抽出対象周波数にそれぞれ対応して複数設置される、ことを特徴とす請求項1又は2に記載の周波数解析装置。抽出対象周波数の1周期分に相当する所定時間毎に、その所定時間に入力される入力信号を信号素片として保持する入力処理と、前記所定時間毎に保持される前記信号素片の、前記所定時間内における各時点の信号値をそれぞれ累積加算した累積加算結果を保持し、前記所定時間が経過する毎に、前記保持された前記累積加算結果と、新たに前記入力された信号素片との、前記所定時間内における同一時点の信号値をそれぞれ加算して、新たな累積加算結果を求め、その新たな累積加算結果を前記保持された累積加算結果の代わりに保持する演算処理と、を行うことを特徴する周波数解析方法。前記演算処理では、前記累積加算した回数で前記累積加算結果を除算することを特徴とする請求項4に記載の周波数解析方法。前記入力処理及び前記演算処理は、複数の抽出対象周波数にそれぞれ対応して行われる、ことを特徴とする請求項4又は5に記載の周波数解析方法。コンピュータに抽出対象周波数の1周期分に相当する所定時間毎に、その所定時間に入力される入力信号を信号素片として保持する入力手順と、前記所定時間毎に保持される前記信号素片の、前記所定時間内における各時点の信号値をそれぞれ累積加算した累積加算結果を保持し、前記所定時間が経過する毎に、前記保持された前記累積加算結果と、新たに前記入力された信号素片との、前記所定時間内における同一時点の信号値をそれぞれ加算して、新たな累積加算結果を求め、その新たな累積加算結果を前記保持された累積加算結果の代わりに保持する演算手順と、を実行させることを特徴とするプログラム。前記演算手順では、前記累積加算した回数で前記累積加算結果を除算させることを特徴とする請求項7に記載のプログラム。前記入力手順及び前記演算手順は、複数の抽出対象周波数にそれぞれ対応して実行させる、ことを特徴とする請求項7又は8に記載のプログラム。FFTやDFTを使用した周波数解析方法では、入力信号を窓関数で処理し、さらに、FFT演算或いはDFT演算を行う必要があり、処理量演算量が大きい。このため、処理能力の大きいDSPDigitalSignalProcessorでは実現できるが、処理能力の小さいDSPでは実行困難である。また、バンドパスフィルタを用いた周波数解析方法では、フィルタの通過域を急峻にしても、余計な周波数までも取得してしまうので解析精度が低い。また、レベル計算やフィルタ計算に高度な演算精度と大きな演算量が必要である。このため、処理能力の大きいDSPDigitalSignalProcessorでは実現できるが、処理能力の小さいDSPでは実行困難である。本発明は、以上のような現状に鑑みてなされたものであり、高精度であると共に処理量の小さい周波数解析の手法を提供することを目的とする。上記目的を達成するために、本発明の第1の観点に係る周波数解析装置は、抽出対象周波数の1周期分に相当する所定時間毎に、その所定時間に入力される入力信号を信号素片として格納する入力バッファと、前記所定時間毎に前記入力バッファから供給される前記信号素片の、前記所定時間内における各時点の信号値をそれぞれ累積加算した累積加算結果を保持するためのレジスタを有し、前記所定時間が経過する毎に、前記レジスタに保持された前記累積加算結果と、新たに前記入力バッファから供給された信号素片との、前記所定時間内における同一時点の信号値をそれぞれ加算して、新たな累積加算結果を求め、その新たな累積加算結果を前記レジスタに保持させる演算部と、を備えることを特徴とする。前記演算部は、例えば、前記累積加算した回数で前記累積加算結果を除算する。また、例えば、前記入力バッファ及び前記演算部は、複数の抽出対象周波数にそれぞれ対応して複数設置される。上記目的を達成するために、本発明の第2の観点に係る周波数解析方法は、抽出対象周波数の1周期分に相当する所定時間毎に、その所定時間に入力される入力信号を信号素片として保持する入力処理と、前記所定時間毎に保持される前記信号素片の、前記所定時間内における各時点の信号値をそれぞれ累積加算した累積加算結果を保持し、前記所定時間が経過する毎に、前記保持された前記累積加算結果と、新たに前記入力された信号素片との、前記所定時間内における同一時点の信号値をそれぞれ加算して、新たな累積加算結果を求め、その新たな累積加算結果を前記保持された累積加算結果の代わりに保持する演算処理と、を行うことを特徴する。上記目的を達成するために、本発明の第3の観点に係るコンピュータプログラムは、コンピュータに、抽出対象周波数の1周期分に相当する所定時間毎に、その所定時間に入力される入力信号を信号素片として保持する入力手順と、前記所定時間毎に保持される前記信号素片の、前記所定時間内における各時点の信号値をそれぞれ累積加算した累積加算結果を保持し、前記所定時間が経過する毎に、前記保持された前記累積加算結果と、新たに前記入力された信号素片との、前記所定時間内における同一時点の信号値をそれぞれ加算して、新たな累積加算結果を求め、その新たな累積加算結果を前記保持された累積加算結果の代わりに保持する演算手順と、を実行させることを特徴とする。この発明により、高精度であると共に処理量の小さい周波数解析の手法を提供することができる。以下、本発明の実施の形態に係る周波数解析方法と周波数解析装置を説明する。本実施形態の周波数解析方法は、入力信号を、目的とする周波数の1周期分ずつに分割し、分割信号内の相対的に同一時点の信号強度瞬時値を平均化することで、目的周波数の信号成分を抽出するものである。抽出した信号成分は、波形なのでそこからレベルを求めることができる。より具体的に説明すると、音楽などの自然界の音は一般的に、音が出てからその音の波形が1周期で消える事は無く、同一の波形がある程度繰り返される。本実施形態の周波数解析方法は、この波形が繰り返される特徴を用いて、所望の正弦波を抜き出す方法である。入力信号を構成する複数の正弦波から、例えば、1kHzの波形を抜き出したい場合、波長である1msec毎に入力信号を区切る。1msec毎に区切った1個目の波形をFrame1と呼び、2個目の波形をFrame2、以降Frame3、Frame4と、呼ぶ事にする。1msec毎に区切ったとき、波長が1msecの信号成分は、各Frameで、波形が同一になる。これに対し、波長が1msec以外の場合は、1msecの時間に波長が合わないのでFrame1とFrame2とを比べると波形がずれる。波長が1msecの正弦波は、Frame1、Frame2、Frame3、の各Frameにおいて波形が同一なので、複数のFrameを平均化Frame内の相対的に同一時点の信号成分の平均化しても波形は変化しない。一方、波長が1msec以外の正弦波については、正弦波は周期性があり且つ1波長の平均は0なので、平均化すると0に近い値となる。これにより、音声信号から所望の周波数成分を、加算と割り算という簡単な処理で、抜き出すことが可能となる。この原理を数式で説明する。周波数がf1の正弦波と周波数がf2の正弦波の合成波をサンプリング周波数fsでサンプリングしてデジタル化し、長さf1/fsのFrameに区分し、8つのFrameの相対位置が同一のデータを平均化すると、その値は次のようになる。ここで、f1は、目的抽出対象信号成分の周波数を示している。f2は、入力信号に含まれている他の周波数信号成分の周波数を示している。tは、各Frameにおける時刻を表し、tの添え字fs/f1は、サンプリングされる間隔を表している。nは、8Frame中の何番目のFrameであるかを示している。抽出対象の周波数がf1の場合、Frameはf1/fsの長さとなる。隣接するするFrame間では、次のFrameでは、360度回転するので、前のFrameと同じ波形が出力される。加算平均をとった式の前半のsinの部分は次式のようになる。Frameの長さに合わない周期を持つf2は、正弦波の性質上±1の範囲内なので、加算平均すると、1より小さくなる。よって、後半のsin部分は∞まで加算すると、次式のようになる。このように、ある時間おきに区切って加算平均すると、その時間にあった周波数はいつも波形が同じ形で出る。異なった周波数は、周期がずれて出るので、加算平均することで打ち消しあいゼロに近づく。結果として、角速度が2πf1/fsとなる周期のみの波形が残る形となる。次に、上記周波数解析方法を実現する周波数解析装置の実施形態について具体的に説明する。[第1の実施形態]本実施形態に係る周波数解析装置100は、図1に示すように、入力端子Inに接続された入力バッファ31と、入力バッファ31に接続された演算部32と、制御部33とを備え、入力端子Inから入力されるデジタル入力信号の表す波形から抽出対象の周波数成分の波形を抽出する装置である。入力バッファ31は、入力端子Inを介して逐次入力されるデジタル入力信号の各データ信号素片を、抽出対象の周波数ftの1周期のT=1/ft時間分格納する容量を有する。演算部32は、加算結果レジスタ321と演算回路329とを備える。加算結果レジスタ321は、抽出対象の周波数ftの1周期T=1/ft分のデータを格納する容量を有する。演算回路329は、入力バッファ31と加算結果レジスタ321の相対的に同一時点のデータを加算して、加算結果レジスタ321に格納する。演算回路329は、加算動作をn回繰り返した後、加算結果レジスタ321に保持されている各データをnで除算し、平均値を求める。制御部33は、マイクロプロセッサとメモリから構成され、全体の動作を制御する。演算部32は、マイクロプロセッサ及びメモリで構成してもよく、或いはDSPで構成してもよい。さらに、入力バッファ31及び演算部32の全体を、マイクロプロセッサ及びメモリや、DSPで構成してもよい。次に、図2~図7を参照しつつ、周波数解析装置の動作を説明する。図2は、周波数解析装置100の制御部33の動作を示すフローチャートである。まず、制御部33は、抽出対象の周波数ftを設定するステップS1。次に、制御部33は、デジタル入力信号のサンプリング周波数fsを設定するステップS2。制御部33は、周波数ftの波形の1周期T=1/ftを求めるステップS3。制御部33は、1周期Tにおけるデータ数m=T/{1/fs}=fs/ftを求めるステップS4。制御部33は、演算部32による入力データの加算回数nを定めるステップS5。次に、制御部33は、加算実行回数iを初期値の0に設定するステップS6。次に、制御部33は、加算結果レジスタ321に格納されている全ての加算データを初期値の0に設定するステップS7。一方、入力バッファ31には、入力端子Inを介してデジタル入力信号が供給される。制御部33は、入力バッファ31に、サンプリング周波数fsと同一周波数の転送クロックを供給し、入力バッファ31に、デジタル入力信号をシフトしつつ格納するステップS8。制御部33は、抽出対象の周波数ftの1周期Tが経過したか否か、即ち、入力バッファ31にm個の入力データが格納されたか否かを判別するステップS9。時間Tが経過していない入力データ数がmに達していないと判別した場合ステップS9:NO、ステップS8にリターンし、入力データの格納動作を継続する。一方、T時間が経過した入力データ数がmに達してたと判別した場合ステップS9:YES、制御部33は演算回路329に加算を指示する。この指示に応答して、演算回路329は、入力バッファ31に格納されたm個の入力データをパラレルに取得し、取得したm個のデータそれぞれを、その時点で加算結果レジスタ321に格納されている対応する加算データに加算し、加算結果レジスタ321に再度格納するステップS10。次に、制御部33は、加算実行回数iを1だけインクリメントするステップS11。制御部33は、加算実行回数iが加算回数nに達したか否かを判別しステップS12、加算実行回数iが所定回数nに達していないと判別した場合にはステップS12:NO、処理をステップS8に戻し、入力データを入力バッファ31に新たにm個格納して、加算結果レジスタ321に格納されている加算データと加算する処理を繰り返す。一方、加算実行回数iが加算回数nに達した場合にはステップS12:YES、加算結果レジスタ321には、デジタル入力信号を周期Tで分断し、n周期分のデータについて、各周期T内の相対的に同一時点のデータ同士を加算した加算データが得られたことになる。制御部33は、演算回路329を制御し、加算結果レジスタ321に格納されている各加算データを加算回数nで割るステップS13。これにより、過去nT期間におけるデジタル入力信号の、各時間T内の相対的に同一時点のデータの平均値が加算結果レジスタ321に求められる。演算部32は、加算結果レジスタ321に格納されたデータを出力するステップS14。その後、停止の指示の有無を判別しステップS15、停止の指示があればステップS15:YES、処理を終了する。一方、停止の指示がなければステップS15:NO、加算結果レジスタ321をリセットして全ての加算データを初期値の0に戻しステップS16、加算実行回数iを初期値の0に戻しステップS17、処理をステップS8に戻す。図3は、入力信号の例を示す図である。図4及び図5は、図3に示す入力信号のデータ値を説明する図である。入力データが、図3のように、周波数が600Hzの正弦波と周波数が3000Hzの正弦波と周波数が2500Hzの正弦波との合成信号を表し、図4に示すように、600Hzの信号成分の振幅が0.25、3000Hzの信号成分の振幅が0.4、2500Hzの信号成分の振幅が0.4の割合で加算したものであるとする。また、サンプリング周波数fsを25000Hzとする。この場合、1サンプリング期間Ts=0.00004sec=0.04msec=40μsecである。2500Hzの成分を抽出対象とすると、その周期Ttは0.0004secとなる。この場合、制御部33は、抽出対象周波数ftを2500Hzに設定しステップS1、サンプリング周波数fsを25000Hzに設定しステップS2、抽出対象周波数ftの一周期Tを0.0004secと求めるステップS3。また、1周期Tの間に入力されるデータの数mを、Tt/Ts=10と特定するステップS4。また、加算回数nを設定するステップS5。ここでは、n=8とする。次に、加算実行回数iを0としステップS6、加算結果レジスタ321に格納されている加算データを全て0とするステップS7。続いて、制御部33は、入力バッファ31に周波数fsのシフトクロックを供給し、入力信号を取りこませる。シフトクロックを、期間Tの間に10パルス供給すると、入力バッファ31には、図4の期間T1で示す10個のデータが格納され、ステップS9でYESと判別される。これにより、入力バッファ31の各格納データと加算結果レジスタ321の格納データの対応するもの同士を加算し、加算結果レジスタ321に格納するステップS10。続いて、加算実行回数iが更新されて1となりステップS11、ステップS12で1=i<n=8と判別されるステップS12:NO。制御部33は、入力バッファ31に10クロック分のシフトクロックを供給し、図4の期間T2で示す10個のデータ第11~20のデータが格納され、入力バッファ31の各格納データと加算結果レジスタ321の格納データの対応するもの同士を加算し、加算結果レジスタ321に格納するステップS10。以後、同様の動作が繰り返される、i=8と判別されまで繰り返されるステップS12:NO。図4及び図5の場合、1回目の周期T1で入力バッファ31に格納される入力信号をδ1と表し、2回目の周期T2で入力バッファ31に格納される入力信号をδ2と表し、各周期Tにおけるデータの位置相対的な時点を[]内の数字で示していくと、δ1[0]=0δ1[1]=0.546489δ1[2]=0.853894δ1[3]=0.797907δ2[0]=0.629929δ2[1]=0.846078δ2[2]=0.77058δ2[3]=0.464442となる。演算部32は、ステップS10の加算を8回行った後に、ステップS13で加算結果を所定数n=8で割っている。尚、この割り算は、加算結果レジスタ321に格納されている各加算データを3ビット桁が少なくなる方向にシフトすることでも、実現可能である。以上の加算及び割り算は、入力バッファ31に格納されたデータ群について、周期T内の相対的に同一時点のデータ同士を加算して割るものであり、次式で表される。D1=δ1[0]δ2[0]δ3[0]δ4[0]δ5[0]δ6[0]δ7[0]δ8[0]/8=00.6299290.266447-0.48069-0.44260.2377640.4724540.008907/8=0.086526375D2=δ1[1]δ2[1]δ3[1]δ4[1]δ5[1]δ6[1]δ7[1]δ8[1]/8=0.5464890.8460780.178698-0.407620.0173220.7556040.653113-0.05465/8=0.31687925平均D1,D2,D3の値を図に示すと、図6のようになり、波形で示すと図7のようになり、抽出対象の2500Hzの成分の正弦波の波形が抽出できる。また、抽出対象の2500Hzの成分の振幅も0.4と再現できている。以上説明したように、本実施の形態によれば、簡単な構成及び制御で入力信号から目的周波数ftの信号成分を抽出することができる。なお、以上の実施の形態においては、2500Hzの周波数成分を抽出する例を示したが、抽出対象の周波数は任意である。但し、抽出対象の周波数ftとサンプリング周波数fsとデータ数mとが、1=mft/fsを満たし、mが2以上の自然数となる関係が成立することが望ましい。本実施形態の周波数解析装置100は、次のような利点を有する。1入力信号の格納と、加算処理と、割り算処理或いはビットシフト処理のみで、実現できるので、処理全体が簡単であり、実現しやすい。2演算量の多いFFT解析を行わなくても、周波数解析が可能である。3バンドパスフィルタ等のフィルタを使用しないので、演算精度を気にせず、簡単に計算が可能である。4バンドパスフィルタ等のフィルタを使用するよりも、精度よく所望の周波数成分を抽出できる。変形例上記実施の形態では、ピーク値振幅をレベルとして出力したが、ピーク値/1.4実効値をレベルとして出力してもよい。[第2の実施形態]図8は、本発明の第2の実施形態に係る周波数解析装置100を示す構成図である。この周波数解析装置100は、入力端子Inに共通に接続された複数の例えば4個の入力バッファ31a,31b,31c,31dと、各入力バッファ31a,31b,31c,31dにそれぞれに接続された演算部32a,32b,32c,32dと、装置全体を制御する制御部33とを備え、4つの抽出対象の周波数のレベルを抽出するものである。4つの抽出対象の周波数の1周期をTa,Tb,Tc,Tdとすると、入力バッファ31aは、入力信号のTa時間分を格納する。入力バッファ31bは、入力信号のTb時間分を格納する。入力バッファ31cは、入力信号のTc時間分を格納する。入力バッファ31dは、入力信号のTd時間分を格納する。各入力バッファ31a,31b,31c,31dにそれぞれ接続された演算部32a,32b,32c,32dは、制御部33に制御されて、入力バッファ31a~31dに格納される入力信号に対して、第1の実施形態の演算部32と同様に動作し、4つの抽出対象の周波数のレベルや波形を抽出できる。したがって、本実施形態の周波数解析装置100は、第1の実施形態と同様の利点を奏すると共に、4つの抽出対象の周波数のレベルや波形を抽出できるという利点を有する。[第3の実施形態]図9は、本発明の第3の実施形態に係る周波数解析装置100を示す構成図である。この周波数解析装置100は、入力端子Inに接続された入力バッファ41と、入力バッファ41に接続された4個の演算部42a,42b,42c,42dと、装置全体を制御する制御部33とを備え、4つの抽出対象の周波数のレベルを抽出するものである。入力バッファ41のサイズは、抽出対象の周波数のうちの周波数が一番低い周波数の1周期分を格納できるサイズとし、制御部33からの制御信号等により、入力信号を格納する長さを選択的4段階に変化できるようにしたものである。入力バッファ41に接続された演算部42a,42b,42c,42dは、第2の実施形態の演算部32a,32b,32c,32dと同様のものである。4つの抽出対象の周波数の1周期をTa,Tb,Tc,Tdとすると、入力バッファ41がTa時間の入力信号を格納するように設定されている期間には、入力バッファ41の格納する入力信号に対して演算部42aが、第1の実施形態の演算部32と同様に動作し、1周期がTaとなる抽出対象の周波数のレベル及び波形を抽出する。入力バッファ41がTb時間の入力信号を格納するように設定されている期間には、入力バッファ41の格納する入力信号に対して演算部42bが、第1の実施形態の演算部32と同様に動作し、1周期がTbとなる抽出対象の周波数のレベル及び波形を抽出する。入力バッファ41がTc時間の入力信号を格納するように設定されている期間には、入力バッファ41の格納する入力信号に対して演算部42cが、第1の実施形態の演算部32と同様に動作し、1周期がTcとなる抽出対象の周波数のレベル及び波形を抽出する。入力バッファ41がTd時間の入力信号を格納するように設定されている期間には、入力バッファ41の格納する入力信号に対して演算部42dが、第1の実施形態の演算部32と同様に動作し、1周期がTdとなる抽出対象の周波数のレベル及び波形を抽出する。則ち、本実施形態の周波数解析装置100は、第2の実施形態の周波数解析装置100の4個の入力バッファ31a,31b,31c,31dを、1個の入力バッファ41で置換したものとなっており、構成が第2の実施形態よりも簡素になるという利点を有する。尚、本発明は、上記実施形態に限定されず、種々の変形が可能である。その変形例としては、例えば次のようなものがある。i上記第1~第3の実施形態では、各入力バッファ31,31a~31d,41に格納された入力信号を複数回加算してその平均を求めることにより、抽出対象の周波数成分のレベルまで求められるようにしたが、入力信号に抽出対象の周波数含まれるか否かを解析する場合には、平均まで求める必要はない。この場合、各入力バッファ31,31a~31d,41で格納される入力信号を複数回加算するだけで、抽出対象以外の周波数成分が除去されるので、入力信号に抽出対象の周波数含まれるか否かを解析することができる。ii上記第2の実施形態では、抽出対象の周波数の数と同数の4個の入力バッファ31a~31d及び演算部32a~32dを備えたが、入力バッファ31a~31d及び演算部32a~32dの数を増減することにより、抽出対象の周波数の数を増減できる。iii上記第3の実施形態では、入力バッファ41の入力信号を格納する長さを4段階に変化させる構成としたが、入力バッファ41の入力信号を格納する長さを4段階以外にしてもよい。この場合、演算部42a~42dの数もそれに応じて変化させることにより、抽出対象の周波数の数を増減できる。iv入力信号は、音声を表す信号でもよいし、画像を表す信号でもよいし、各種センサから出力される信号でもよい。v第1~第3の実施形態においては、8Frame毎に、目的周波数の信号成分の有無及びそのレベルを求めたが、Frame単位でその有無レベルを求めるようにしてもよい。以下、この種の処理が可能な周波数解析装置を第4の実施形態として示す。[第4の実施形態]図10は、本発明の第4の実施形態に係る周波数解析装置200を示す構成図である。周波数解析装置200においては、演算部32は、加算結果レジスタ321,演算回路329と、n個のレジスタ341~34nと、装置全体を制御する制御部33とから構成される。レジスタ341~34nは、カスケードに接続されており、制御部33の制御に従って、自己が格納しているm個のデータを後段のレジスタに転送すると共に前段のレジスタから供給されるm個のデータを新たに記憶する。演算回路329は、レジスタ341~34nに格納されているnFrame分のデータのうち、Frame内の相対位置が同一のもの同士を加算し、加算値をnで割った値を加算結果レジスタ321に格納し、これが演算部32の出力となる。入力バッファ31に1Frame分のm個のデータが格納される度に、そのデータがレジスタ341に格納され、さらに、レジスタ341~34n-1に格納されていたデータが後段のレジスタ342~34nに格納される。演算回路329は、Frame単位で、過去nフレーム分の入力信号を平均化し、加算結果レジスタ321に格納する。従って、Frame単位で、目的とする周波数の信号がふくまれているか否か、含まれている場合には、そのレベルを求めることができる。尚、本発明は上記した装置の機能をコンピューターに実現させるためのプログラムを含むものである。これらのプログラムは、記録媒体から読み取られてコンピュータに取り込まれてもよいし、通信ネットワークを介して伝送されてコンピュータに取り込まれてもよい。"