(57)【要約】【課題】ゼロシフト成分からなる偽信号成分を含む原信号から、偽信号成分を効果的に除去する。【解決手段】原信号を時間で積分して被積分原信号を導出し、ウェーブレット変換により被積分原信号を１つのＡｐｐｒｏｘｉｍａｔｉｏｎと１つ以上のＤｅｔａｉｌとの和に分解した場合に、当該１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの原信号のエネルギーに対する比率が所定の閾値と所定の関係になる閾レベル数に達するまでレベル数を深めながらウェーブレット変換を実行し、閾レベル数におけるＡｐｐｒｏｘｉｍａｔｉｏｎをゼロに近似させ、Ａｐｐｒｏｘｉｍａｔｉｏｎ及び閾レベル数までのＤｅｔａｉｌをウェーブレット逆変換により再構成して時間で微分することによって偽信号が除去された主信号を導出する。【選択図】　　　　　図１
2004093455.tif 000002
【特許請求の範囲】
【請求項１】
ゼロシフト成分を有しない主信号とゼロシフト成分からなる偽信号とを含む原信号から当該偽信号を除去する偽信号除去方法であって、
前記原信号を時間で積分することによって被積分原信号を導出するステップと、
ウェーブレット変換により前記被積分原信号を１つのＡｐｐｒｏｘｉｍａｔｉｏｎと１つ以上のＤｅｔａｉｌとの和に分解した場合に、当該１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの前記原信号のエネルギーに対する比率が所定の閾値と所定の関係になる閾レベル数に達するまでレベル数を深めながらウェーブレット変換を実行するステップと、
前記閾レベル数における前記Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分として閾値処理によりゼロに近似させるステップと、
前記閾値処理されたＡｐｐｒｏｘｉｍａｔｉｏｎ及び前記閾レベル数までのＤｅｔａｉｌをウェーブレット逆変換により再構成することによって被積分主信号を導出するステップと、
前記被積分主信号を時間で微分することによって前記主信号を導出するステップと、を有することを特徴とする偽信号除去方法。
【請求項２】
前記閾レベル数は、前記１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの前記原信号のエネルギーに対する比率が所定の閾値より小さくなるようなレベル数の内で最小のものであることを特徴とする請求項１に記載の偽信号除去方法。
【請求項３】
前記閾値処理は、前記閾レベル数における前記Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分としてゼロにすることを特徴とする請求項１又は２に記載の偽信号除去方法。
【請求項４】
前記原信号は、衝撃試験でセンサによって測定される加速度信号であり、及び
前記偽信号は、前記衝撃試験の衝撃によって前記センサに誘導されるゼロシフト成分であることを特徴とする請求項１から３のいずれか１項に記載の偽信号除去方法。
【請求項５】
ゼロシフト成分を有しない主信号とゼロシフト成分からなる偽信号とを含む原信号から当該偽信号を除去するため、コンピュータに、
前記原信号を時間で積分することによって被積分原信号を導出するステップと、
ウェーブレット変換により前記被積分原信号を１つのＡｐｐｒｏｘｉｍａｔｉｏｎと１つ以上のＤｅｔａｉｌとの和に分解した場合に、当該１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの前記原信号のエネルギーに対する比率が所定の閾値と所定の関係になる閾レベル数に達するまでレベル数を深めながらウェーブレット変換を実行するステップと、
前記閾レベル数における前記Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分として閾値処理によりゼロに近似させるステップと、
前記閾値処理されたＡｐｐｒｏｘｉｍａｔｉｏｎ及び前記閾レベル数までのＤｅｔａｉｌをウェーブレット逆変換により再構成することによって被積分主信号を導出するステップと、
前記被積分主信号を時間で微分することによって前記主信号を導出するステップと、を実行させるための偽信号除去プログラム。
【請求項６】
前記閾レベル数は、前記１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの前記原信号のエネルギーに対する比率が所定の閾値より小さくなるようなレベル数の内で最小のものであることを特徴とする請求項５に記載の偽信号除去プログラム。
【請求項７】
前記閾値処理は、前記Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分としてゼロにすることを特徴とする請求項５又は６に記載の偽信号除去プログラム。
【請求項８】
前記原信号は、衝撃試験でセンサによって測定される加速度信号であり、及び
前記偽信号は、前記衝撃試験の衝撃によって前記センサに誘導されるゼロシフト成分であることを特徴とする請求項５から８のいずれか１項に記載の偽信号除去プログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、信号処理方法に関し、より詳しくは、信号に含まれるエネルギーを多く含む偽信号の除去方法に関する。
【０００２】
【従来の技術】
従来のフーリエ変換による信号処理方法では、周波数領域で信号を把握、処理することに優れており、多くの分野で使われている技術である。しかし、フーリエ変換は、近年各分野で研究開発が進むウェーブレット変換による信号処理方法とは異なり、時間領域で信号を把握、処理することには適していない。人工衛星等の衝撃試験に於いて、計測された加速度信号に偽信号が含まれる場合、従来の技術では偽信号を除去するための有効な方法が無く、従って、偽信号発生の程度の大きいものに関しては試験データを評価対象外とせざるを得なかった。
【０００３】
従来のフーリエ変換による信号処理方法は、信号の周期性を仮定していること、及び周波数領域での信号処理方法であるため、試験データとして有効な信号成分（以下、「主信号成分」）と偽信号成分とを分離して処理することが出来ない。また、ウェーブレット変換による信号処理方法をそのまま用いた場合、主信号成分に対して、偽信号成分が顕著でない場合、主信号成分まで除去してしまうという欠点があった。
【０００４】
【発明が解決しようとする課題】
本発明は、上記の課題に鑑みてなされたものであり、ゼロシフト成分からなる偽信号成分を含む原信号から、偽信号成分を除去することを目的とする。
【０００５】
【課題を解決するための手段】
上記の課題は、以下のような特徴を有する本発明によって解決される。請求項１に記載の発明は、ゼロシフト成分を有しない主信号とゼロシフト成分からなる偽信号とを含む原信号から当該偽信号を除去する偽信号除去方法であって、当該原信号を時間で積分することによって被積分原信号を導出するステップと、ウェーブレット変換により当該被積分原信号を１つのＡｐｐｒｏｘｉｍａｔｉｏｎと１つ以上のＤｅｔａｉｌとの和に分解した場合に、当該１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの当該原信号のエネルギーに対する比率が所定の閾値と所定の関係になる閾レベル数に達するまでレベル数を深めながらウェーブレット変換を実行するステップと、当該閾レベル数における当該Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分として閾値処理によりゼロに近似させるステップと、当該閾値処理されたＡｐｐｒｏｘｉｍａｔｉｏｎ及び当該閾レベル数までのＤｅｔａｉｌをウェーブレット逆変換により再構成することによって被積分主信号を導出するステップと、当該被積分主信号を時間で微分することによって当該主信号を導出するステップと、を有することを特徴とする。
【０００６】
請求項２に記載の発明は、請求項１に記載の発明の特徴に加えて、当該閾レベル数は、当該１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの当該原信号のエネルギーに対する比率が所定の閾値より小さくなるようなレベル数の内で最小のものであることを特徴とする。
【０００７】
請求項３に記載の発明は、請求項１又は２に記載の発明の特徴に加えて、当該閾値処理は、当該閾レベル数における当該Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分としてゼロにすることを特徴とする。
【０００８】
請求項４に記載の発明は、請求項１から３のいずれか１項に記載の発明の特徴に加えて、当該原信号は、衝撃試験でセンサによって測定される加速度信号であり、及び当該偽信号は、当該衝撃試験の衝撃によって当該センサに誘導されるゼロシフト成分であることを特徴とする。
【０００９】
請求項５に記載の発明は、ゼロシフト成分を有しない主信号とゼロシフト成分からなる偽信号とを含む原信号から当該偽信号を除去するため、コンピュータに、当該原信号を時間で積分することによって被積分原信号を導出するステップと、ウェーブレット変換により当該被積分原信号を１つのＡｐｐｒｏｘｉｍａｔｉｏｎと１つ以上のＤｅｔａｉｌとの和に分解した場合に、当該１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの当該原信号のエネルギーに対する比率が所定の閾値と所定の関係になる閾レベル数に達するまでレベル数を深めながらウェーブレット変換を実行するステップと、当該閾レベル数における当該Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分として閾値処理によりゼロに近似させるステップと、当該閾値処理されたＡｐｐｒｏｘｉｍａｔｉｏｎ及び当該閾レベル数までのＤｅｔａｉｌをウェーブレット逆変換により再構成することによって被積分主信号を導出するステップと、当該被積分主信号を時間で微分することによって当該主信号を導出するステップと、を実行させることを特徴とする。
【００１０】
請求項６に記載の発明は、請求項５に記載の発明の特徴に加えて、当該閾レベル数は、当該１つのＡｐｐｒｏｘｉｍａｔｉｏｎのエネルギーの当該原信号のエネルギーに対する比率が所定の閾値より小さくなるようなレベル数の内で最小のものであることを特徴とする。
【００１１】
請求項７に記載の発明は、請求項５又は６に記載の発明の特徴に加えて、当該閾値処理は、当該閾レベル数における当該Ａｐｐｒｏｘｉｍａｔｉｏｎを偽信号成分としてゼロにすることを特徴とする。
【００１２】
請求項８に記載の発明は、請求項５から７のいずれか１項に記載の発明の特徴に加えて、当該原信号は、衝撃試験でセンサによって測定される加速度信号であり、及び当該偽信号は、当該衝撃試験の衝撃によって当該センサに誘導されるゼロシフト成分であることを特徴とする。
【００１３】
【発明の実施の形態】
まず、ウェーブレットの一般論について説明する。ウェーブレット分解とは、マザー・ウェーブレットと称される波形関数のスケーリング（拡大）とシフト（移動）で信号をそれぞれが部分的に近似するウェーブレット係数の数列に分解し、それで信号を表現する手法である。通常これらのマザー・ウェーブレットは、ローカル的な波形（サポート・コンパクト）である。これによって、変動の激しい信号の波形は小さなスケールで、変動のゆっくりした信号の波形は大きなスケールで表現することができる。一般的にゼロシフトのような偽信号成分は、変動のゆっくりした波形に偏っているという特徴があるので、それを利用してゼロシフト成分を弁別し、除去することができる。一般的なウェーブレット解析の処理手順を以下に示す。まず、信号ｓ（ｔ）からウェーブレット変換により、以下の式で表わされるウェーブレット係数Ｃ（ｊ，ｋ）を求める。ｓ（ｔ）は、一般に時系列の波形サンプル値の数列として表わされる。
【００１４】
【数１】
2004093455.tif 000003
ここで、離散化スケールはａ＝２ｊ（ａ＝１，２，３，…，Ｌ（Ｌはウェーブレットのレベル（階層数）））、シフトはｂ＝ｋ２ｊ（ｋ＝１，２，３，…，Ｌ）、及びΨはマザー・ウェーブレットである。本発明では、Ｄａｕｂｅｃｈｉｅｓウェーブレットを使用し、ウェーブレット係数Ｃ（ｊ，ｋ）によって信号を以下の式に示すように階層で分離する。
【００１５】
【数２】
2004093455.tif 000004
ここで、
【数３】
2004093455.tif 000005
である。ＡＪ（ｔ）は、変動がゆっくりした、スケールの大きな信号成分を表現する。そして、信号ＡＪ（ｔ）と原信号ｓ（ｔ）とのエネルギー比
【数４】
2004093455.tif 000006
を計算して、エネルギー比ＥＪが所定の閾値と所定の関係になるレベル数（本願では「閾レベル数」と称する。）に達するまでレベル数を深めながらウェーブレット変換を実行する。すなわち、エネルギー比ＥＪを必要なレベル数を見極めるための評価基準とする。好適には、そのような閾レベル数は、エネルギー比ＥＪが所定の閾値より小さくなるようなレベル数の内で最小のものである。従って、エネルギー比ＥＪが所定の閾値より小さくなるようなレベルの内で最小のものをレベルＪとして決定する。閾値は、例えば９９．７％などとすることができる。これは、変動のゆっくりしたＡ（ｔ）のエネルギーは、原信号ｓ（ｔ）のエネルギーの大部分を占めるためである。このレベルＪを閾レベル数と呼ぶことにする。ＡＪは、偽信号成分とみなすことができる。
【００１６】
以下に、信号分離のための具体的な数値計算方法の一連の手順を説明する。
１）ウェーブレットの種類を選択
本発明では、Ｄｂ３（Ｄａｕｂｅｃｈｉｅｓウェーブレット、次数Ｎ＝３）を選択した。以下の表に、様々なＮの値に対応する離散マザー・ウェーブレット関数点の値を示す。
【００１７】
【表１】
2004093455.tif 000007
【００１８】
２）分解、再構成のローパス・フィルタ、ハイパス・フィルタの計算
Ｄｂ３の離散点ｗ（ｋ），ｋ＝１，２，．．．，６を用いてフィルタを計算する。
【数５】
2004093455.tif 000008
2004093455.tif 000009
３）２進ダウンサンプリング
以下の式のように、信号ｓ（ｔ）とフィルタの畳み込み積を計算し、２進ダウンサンプリングを実施する。
【００１９】
【数６】
2004093455.tif 000010
２進ダウンサンプリングは以下の２進ダウンサンプリング関数により実施する。
【００２０】
【数７】
2004093455.tif 000011
2004093455.tif 000012
この操作は、データ数を半分に減らすことでもある。レベルｊから下位レベルｊ＋１への操作を以下のシステム流れ図で表す。
【数８】
2004093455.tif 000013
【００２１】
４）構成信号のレベルＪの決定
エネルギー比
【数９】
2004093455.tif 000014
より、構成信号のレベルＪを決定する。そのためには、エネルギー比ＥＪが所定の閾値と所定の関係になる閾レベル数に達するまでレベル数を深めながらウェーブレット変換を実行する。すなわち、エネルギー比ＥＪを必要なレベル数を見極めるための評価基準とする。好適には、そのような閾レベル数は、エネルギー比ＥＪが所定の閾値より小さくなるようなレベル数の内で最小のものである。それは、このような閾レベル数のウェーブレット分解によって、偽信号成分のほとんどを表わすＡｐｐｒｏｘｉｍａｔｉｏｎを得られるからである。従って、エネルギー比ＥＪが所定の閾値より小さくなるようなレベルの内で最小のものをレベルＪとして決定する。本発明では、積分した信号のエネルギーは、偽信号成分によって支配されるため、例えばエネルギー比の閾値を９９．７％にする。このようにすることによって、偽信号成分を弁別することができる。すなわち、レベルＪのＡｐｐｒｏｘｉｍａｔｉｏｎを除去することによって、偽信号成分を除去する。
【００２２】
５）アップサンプリング
データの間（偶数間）にゼロを挿入するアップサンプリングにより、データ数を２倍にする。
【００２３】
６）信号の再構成
再構成ローパス・フィルタ、ハイパス・フィルタとの畳込み積により、２つの信号を再構成する。
【００２４】
７）信号の構築
上記の６）で作成した２つの信号を加え合わせて信号を構築する。レベルｊ＋１から上位レベルｊへの操作を以下のシステム流れ図で表す。
【数１０】
2004093455.tif 000015
以上が、一般的なウェーブレット解析の手順である。
【００２５】
次に、上記の手順を具体的に応用した、偽信号除去方法について説明する。図１は、本発明に係る偽信号除去方法のフロー図である。この方法では、加速度信号を積分することによって偽信号成分が強調されて、元の信号に対する偽信号成分のエネルギー比が大きくなる性質を利用する。本発明に係る方法は、コンピュータなどを用いて実施することができる。すなわち、コンピュータに原信号を入力し、コンピュータの演算手段を使用して、その原信号に、積分操作、ウェーブレット分解、閾値処理、再構成、及び微分操作を行って主信号を導出することによって、実施することができる。
【００２６】
偽信号を含む原信号をａ（ｔ）とする。原信号は、図２（ａ）に示すような衝撃試験で圧電型加速度センサ１３によって測定される加速度信号である。まず、原信号ａ（ｔ）を積分し、速度信号ｖ（ｔ）を得る（ステップＳ１０１）。
【００２７】
次に速度信号ｖ（ｔ）に対してウェーブレット分解を施して、ウェーブレット係数を得る。ウェーブレット分解のパラメータは、例えばマザーウェーブレットにＤａｕｂｅｃｈｉｅｓ　Ｎ＝３とする。なお、このウェーブレット係数から不要な成分を除いた後に再構築を行って元の信号の構成信号を得る手続きを「ウェーブレット処理」と呼ぶ事とする。速度信号ｖ（ｔ）を所定のウェーブレット分解によって偽信号成分（ウェーブレット変換でいうところのＡｐｐｒｏｘｉｍａｔｉｏｎ）Ａ（ｔ）と主信号成分（ウェーブレット変換でいうところのＤｅｔａｉｌ）Ｄ（ｔ）に分解する（ステップＳ１０３）。次に、所定のウェーブレット処理を実施する。ここで所定のウェーブレット処理とは、以下に示す条件を満たすような処理である。まず、Ａ（ｔ）の元の信号ｖ（ｔ）に対するエネルギー比Ｅ＝｜｜Ａ（ｔ）｜｜２／｜｜ｖ（ｔ）｜｜２を計算する。そして、エネルギー比Ｅがある閾値より小さくなるようなＡ（ｔ）を求め、そのようなレベル数の内で最小のものである閾レベル数を求める。これにより、その閾レベル数でのウェーブレット分解により、その閾レベル数個のＤ（ｔ）と１つのＡ（ｔ）が得られる。ここでこの閾値は、偽信号成分が元の信号に対して占めるエネルギー比の閾値を表し、例えば９９．７％などとすることができる。このようにできる理由は、本発明に係る方法では、加速度信号である原信号を積分して求めた速度信号においては、主信号はゼロシフト成分を有しないために積分してもほぼゼロであるために、偽信号成分が大きなエネルギーを占めるようになるためである。即ちこのようなＡ（ｔ）を偽信号成分とみなすことができる。上記の所定のウェーブレット処理とは、この閾レベル数のＡ（ｔ）を後述の閾値処理でゼロに近似した後に、Ａ（ｔ）とＤ（ｔ）とで信号を再構成することである。なお、ウェーブレット分解のレベルは、上記の閾レベル数が求まるまでの深さまでで充分である。
【００２８】
これらのＤ（ｔ）とＡ（ｔ）に対してウェーブレット処理を実行する（ステップＳ１０５）。ウェーブレット処理は、偽信号成分の除去と再構成である。偽信号成分の除去については、まず、閾レベル数のウェーブレット分解におけるＡ（ｔ）を偽信号成分として除去する（ゼロにする）方法がある。すなわち、残りのＤ（ｔ）が主信号を表わすこととなる。ここで、偽信号成分Ａ（ｔ）を完全には除去せずに、絶対値が特定の閾値以内の偽信号成分Ａ（ｔ）をゼロにする、いわゆる閾値処理を行ってもよい。その際の閾値処理は、ハードスレッショルド処理（絶対値が特定の閾値内のところはゼロとし、絶対値が特定の閾値外のところは変更しない処理。絶対値が特定の閾値のところでＡ（ｔ）が不連続となる。）でも、ソフトスレッショルド処理（絶対値が特定の閾値内のところはゼロとし、絶対値が特定の閾値外のところは徐々にゼロに近づける処理。絶対値が特定の閾値のところでＡ（ｔ）が連続となる。）でもよい。偽信号成分が除去された信号は、ウェーブレット再構成されて、主信号成分に対応する速度信号ｖｃ（ｔ）として出力される。ｖｅｒｒ（ｔ）は、Ａ（ｔ）を偽信号成分として除去した場合、Ｄ（ｔ）を再構成したものとなる。一方、偽信号成分は、ウェーブレット再構成されて、偽信号に対応する速度信号ｖｅｒｒ（ｔ）として出力される。ｖｅｒｒ（ｔ）は、Ａ（ｔ）を偽信号成分として除去した場合、Ａ（ｔ）を再構成したものとなる。
【００２９】
元の速度信号ｖ（ｔ）から偽信号成分Ａ（ｔ）が除去された主信号成分に対応するｖｃ（ｔ）を時間微分して、補正された（偽信号成分が除去された）加速度信号ａｃ（ｔ）を得る（ステップＳ１０７）。このようにして、ゼロシフト成分からなる偽信号が除去された加速度信号ａｃ（ｔ）を得ることができる。なお偽信号が表わす速度信号ｖｅｒｒ（ｔ）を時間微分すると、偽信号Ａｅｒｒ（ｔ）が得られる（ステップＳ１０９）。
【００３０】
ウェーブレット変換の性質上、マザーウェーブレット及びそのレベルに関しては、対象となる信号の性質に応じて、他のマザーウェーブレット及びレベルを選択することも可能である。
【００３１】
【実施例】
図２（ａ）に示すように、圧電型加速度センサー１３によって計測した、偽信号を含む衝撃応答加速度信号ａ（ｔ）を本発明に係る信号処理の対象とした。なお比較例として、ゼロシフトによる偽信号が原理的に生じない特公平６−５２２７０号公報に開示された衝撃加速度計の動的応答特性測定方法を用いて、図２（ｂ）に示されるように歪みゲージによる衝撃応答加速度信号を測定し、それを基準として本発明による偽信号除去の効果を定量化した。
【００３２】
図３に、これから偽信号を除去しようとする衝撃応答加速度信号ａ（ｔ）を表わすグラフを示す。図３に示す衝撃応答加速度信号ａ（ｔ）を積分し、速度信号ｖ（ｔ）（図４）を得る。ａ（ｔ）、ｖ（ｔ）は、丸棒端面での反射波の影響でこのような形状となる。ｖ（ｔ）をウェーブレット分解し（この例ではウェーブレットとしてＤａｕｂｅｃｈｉｅｓ　Ｎ＝３を用いた）、主信号成分Ｄ（ｔ）（図５）と偽信号成分Ａ（ｔ）（図６）に分解した。この際、偽信号成分Ａ（ｔ）の元の信号ｖ（ｔ）に対するエネルギー比Ｅを９９．５７％として、閾値処理を行った。すなわち、エネルギー比Ｅが９９．５７％より小さくなるようなレベルの内で最小のものである閾レベル数におけるウェーブレット分解で得られたＡ（ｔ）を偽信号成分として除去した。偽信号成分Ａ（ｔ）が除去された主信号成分Ｄ（ｔ）は、補正された速度信号ｖｃ（ｔ）に対応するが、これを時間微分して、補正された加速度信号ａｃ（ｔ）（図７）を得た。なお、偽信号成分Ａ（ｔ）を時間微分すると、偽信号成分の加速度信号ａｅｒｒ（ｔ）（図８）を得ることができる。ａ（ｔ），ａｃ（ｔ），ａｅｒｒ（ｔ）及び正解値となる歪みゲージによって計測した比較例による衝撃応答加速度信号ａｔｒｕｅ（ｔ）のＳｈｏｃｋ　Ｒｅｓｐｏｎｓｅ　Ｓｐｅｃｔｒｕｍ（ＳＲＳ）解析結果を図９に示す。図９中で、「○」印がａ（ｔ）、「＋」印がａｃ（ｔ）、「□」印がａｔｒｕｅ（ｔ）、「▲」印がａｅｒｒ（ｔ）のＳＲＳである。この図９より、ａｃ（ｔ）とａｔｒｕｅ（ｔ）とがよく一致しており、本発明に係る方法によって偽信号成分が除去されたことが分かる。
【００３３】
本発明に係る方法によって、従来の手法に対してゼロシフト除去の精度があがったことは以下のことから確認できる。図１０は、元のゼロシフトした加速度信号（上側の図）及びそのＳＲＳ（下側の図）を示す。加速度信号を見ると、それは初期の大振幅の振動から一定値に収束しているが、そのいずれの時点においても、中心がゼロより大きく上側に偏っており、ゼロシフトしていることがわかる。また、ＳＲＳを見ると、「○」で示されたポジティブ成分と「△」で示されたネガティブ成分とが大きく離れており、加速度信号がポジティブ側に偏っていることが一見して確認できる。図１１は、従来の手法に基づき、加速度信号にそのまま（積分を行わずに）ウェーブレット処理を施して解析した結果の加速度信号（上側の図）及びそのＳＲＳ（下側の図）を示す。加速度信号を見ると、それは初期の大振幅の振動から一定値に収束しているが、特に初期においては、中心がゼロより若干上側に偏っており、ゼロシフトしていることがわかる。また、ＳＲＳを見ると、「○」で示されたポジティブ成分と「△」で示されたネガティブ成分とが一致してはおらず、加速度信号がポジティブ側に若干偏っていることがやはり一見して確認できる。図１２は、本発明に係る手法によって、加速度信号を積分して速度信号に変換しゼロシフト成分を強調した上でウェーブレット処理を施した結果である。加速度信号を見ると、それは初期の大振幅の振動から一定値に収束しているが、そのいずれの時点においてもは、中心はゼロ近傍に位置しており、ゼロシフトがほとんど見られないことがわかる。また、ＳＲＳを見ると、「○」で示されたポジティブ成分と「△」で示されたネガティブ成分とが全域においてほぼ一致しており、加速度信号にゼロシフトによる偏りがなく、偽信号成分が効果的に除去されていることが一見して確認できる。これらが示すように、本手法によって従来の手法より精度が大きく向上している。
【００３４】
【発明の効果】
上記のように、本発明によれば、ゼロシフト成分からなる偽信号成分を含む原信号から、偽信号成分を効果的に除去することができるという効果を得ることができる。
【図面の簡単な説明】
【図１】本発明に係る偽信号除去方法のフロー図である。
【図２】（ａ）は本発明の実施例に用いた衝撃加速度計の動的応答特性測定方法に基づく実験装置の概略図であり、及び（ｂ）は本発明の比較例に用いた特公平６−５２２７０号公報に開示された衝撃加速度計の動的応答特性測定方法に基づく実験装置の概略図である。
【図３】偽信号を含む衝撃応答加速度信号ａ（ｔ）を表わすグラフである。
【図４】偽信号を含む衝撃応答速度信号ｖ（ｔ）を表わすグラフである。
【図５】ウェーブレットによって分解された主信号成分Ａ（ｔ）を表わすグラフである。
【図６】ウェーブレットによって分解された偽信号成分Ｄ（ｔ）を表わすグラフである。
【図７】偽信号成分が除去された衝撃応答加速度波形ａｃ（ｔ）を表わすグラフである。
【図８】偽信号成分の衝撃応答加速度波形ａｅｒｒ（ｔ）を表わすグラフである。
【図９】Ｓｈｏｃｋ　Ｒｅｓｐｏｎｓｅ　Ｓｐｅｃｔｒｕｍ（ＳＲＳ）による偽信号除去の評価結果を表わすグラフである。
【図１０】ゼロシフトを含む加速度信号とそのＳＲＳを表わすグラフである。
【図１１】従来の方法によって偽信号除去を行った加速度信号とそのＳＲＳを表わすグラフである。
【図１２】本方法によって偽信号除去を行った加速度信号とそのＳＲＳを表わすグラフである。
【符号の説明】
１１，２１　衝撃発生用の弾丸
１２，２２　丸棒
１３，２３　圧電型加速度センサ
２４　　　　ひずみゲージ
２５　　　　距離Ｌ
ａ（ｔ）　　偽信号を含む衝撃応答加速度信号。
ｖ（ｔ）　　偽信号を含む衝撃応答速度信号。
Ｄ（ｔ）　　ウェーブレットによって分解された主信号成分。
Ａ（ｔ）　　ウェーブレットによって分解された偽信号成分。
ａｃ（ｔ）　　偽信号成分が除去された衝撃応答加速度信号。
ａｅｒｒ（ｔ）　偽信号成分の衝撃応答加速度信号。
【図１】
2004093455.tif 000016
【図２】
2004093455.tif 000017
【図３】
2004093455.tif 000018
【図４】
2004093455.tif 000019
【図５】
2004093455.tif 000020
【図６】
2004093455.tif 000021
【図７】
2004093455.tif 000022
【図８】
2004093455.tif 000023
【図９】
2004093455.tif 000024
【図１０】
2004093455.tif 000025
【図１１】
2004093455.tif 000026
【図１２】
2004093455.tif 000027