供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、フレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したドライブ波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記フレーム波形生成手段は、制御用PSDの第1の周波数成分に初期位相を与え、隣接する周波数成分には、正規分布に従ったランダムな位相差を順次加算して算出した位相を与えて逆フーリエ変換することにより非ガウス特性のランダム波形を生成すること、を特徴とする振動制御装置。請求項1の振動波形装置において、前記フレーム波形生成手段が、フレーム波形のピーク位置がランダムになるように正規分布の中心値を設定すること、を特徴とする振動制御装置。供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、フレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したドライブ波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記フレーム波形生成手段は、目標非ガウス特性に基づいて、制御用PSDの第1の周波数成分に初期位相を与え、隣接する周波数成分には、正規分布に従ったランダムな位相差を順次加算して算出した位相を与えて逆フーリエ変換することにより非ガウス特性のランダム波形を生成してフレーム波形を生成すること、を特徴とする振動制御装置。供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、フレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したドライブ波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記フレーム波形生成手段は、応答波形から得られた非ガウス特性と目標非ガウス特性とを比較して得た非ガウス特性に基づいて、制御用PSDの第1の周波数成分に初期位相を与え、隣接する周波数成分には、正規分布に従ったランダムな位相差を順次加算して算出した位相を与えて逆フーリエ変換することにより非ガウス特性のランダム波形を生成してフレーム波形を生成すること、を特徴とする振動制御装置。請求項4の振動制御装置において、さらに、前記応答波形から非ガウス特性値を算出して、前記非ガウス特性制御手段に与える非ガウス特性算出手段を備えたこと、を特徴とする振動制御装置。供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、ガウス特性のフレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したランダム波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記ガウス特性のランダム波形を、入力された非ガウス特性に基づいて、非ガウス特性のランダム波形に変換する非ガウス化変換手段、非ガウス特性を前記非ガウス化変換手段に与える非ガウス特性制御手段、を備えたこと、を特徴とする振動制御装置。請求項1~5の何れかの振動制御装置において、前記フレーム波形生成手段が、クルトシスに対応する標準偏差の正規分布に基づいて設定された位相を受けること、を特徴とする振動制御装置。請求項1~5の何れかの振動制御装置において、前記フレーム波形生成手段が、スキューネスに対応する初期位相を受けること、を特徴とする振動制御装置。供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換するとともに、窓関数を掛けることにより、非ガウス特性を有するフレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成した非ガウス特性を有するドライブ波形を出力する重ね合わせ波形生成手段、を備えたこと、を特徴とする振動制御装置。しかしながら、特許文献1の振動制御装置では、目標とするスペクトルを有する振動を供試体4に与えることはできるものの、その振動の確率密度分布は、図14に実線で示すガウス分布正規分布となる。しかし、現実の振動は非ガウス分布図14に点線で示すとなることが多いにもかかわらず、従来の振動制御装置では、非ガウス特性を有する振動を供試体に与えることはできなかった。この発明は、上記の問題を解決して、目標とするスペクトルだけでなく、目標とする非ガウス特性を有する振動を供試体4に与えることを目的とする。さらに詳しくは、目標とする非ガウス特性を有する振動を供試体4に与えるよう制御することができる振動制御装置を提供することを目的とする。1この発明の振動制御装置は、供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、フレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したドライブ波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記フレーム波形生成手段は、制御用PSDの第1の周波数成分に初期位相を与え、隣接する周波数成分には、正規分布に従ったランダムな位相差を順次加算して算出した位相を与えて逆フーリエ変換することにより非ガウス特性のランダム波形を生成すること、を特徴とする。これにより、非ガウス特性を有する振動を供試体に与えて制御することが可能な振動制御装置を提供することができる。2この発明の振動制御装置は、前記フレーム波形生成手段が、フレーム波形のピーク位置がランダムになるように正規分布の中心値を設定すること、これにより、非ガウス特性の劣化を低減することができる。3この発明の振動制御装置は、供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、フレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したドライブ波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記フレーム波形生成手段は、目標非ガウス特性に基づいて、制御用PSDの第1の周波数成分に初期位相を与え、隣接する周波数成分には、正規分布に従ったランダムな位相差を順次加算して算出した位相を与えて逆フーリエ変換することにより非ガウス特性のランダム波形を生成してフレーム波形を生成すること、を特徴とする。これにより、応答波形の非ガウス特性が目標非ガウス特性と一致するようにドライブ波形を制御することができる。4この発明の振動制御装置は、供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、フレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したドライブ波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記フレーム波形生成手段は、応答波形から得られた非ガウス特性と目標非ガウス特性とを比較して得た非ガウス特性に基づいて、制御用PSDの第1の周波数成分に初期位相を与え、隣接する周波数成分には、正規分布に従ったランダムな位相差を順次加算して算出した位相を与えて逆フーリエ変換することにより非ガウス特性のランダム波形を生成してフレーム波形を生成すること、を特徴とする。これにより、応答波形から得られた非ガウス特性と目標非ガウス特性とを比較して得た非ガウス特性に一致するように、ドライブ波形を正確にフィードバック制御することができる。5この発明の振動制御装置は、前記応答波形から非ガウス特性値を算出して、前記非ガウス特性制御手段に与える非ガウス特性算出手段を備えたこと、を特徴とする。これにより、応答波形から算出される非ガウス特性値に基づいて、ドライブ波形を正確にフィードバック制御することができる。6この発明の振動制御装置は、供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換することにより、ガウス特性のフレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成したランダム波形を出力する重ね合わせ波形生成手段、を備えた振動制御装置において、前記ガウス特性のランダム波形を、入力された非ガウス特性に基づいて、非ガウス特性のランダム波形に変換する非ガウス化変換手段、非ガウス特性を前記非ガウス化変換手段に与える非ガウス特性制御手段、を備えたこと、を特徴とする。これにより、応答の周波数成分が変化する場合でも、応答の非ガウス特性が目標値と一致するようにドライブ波形を正確にフィードバック制御することができる。7この発明の振動制御装置は、前記フレーム波形生成手段が、クルトシスに対応する標準偏差の正規分布に基づいて設定された位相を受けること、を特徴とする。これにより、制御しようとするクルトシスを有する非ガウス特性波形を生成することができる。8この発明の振動制御装置は、前記フレーム波形生成手段が、スキューネスに対応する初期位相を受けること、を特徴とする。これにより、制御しようとするスキューネスを有する非ガウス特性波形を生成することができる。9この発明の振動制御装置は、供試体の加速度を測定する加速度センサからの応答波形をフーリエ変換して、応答PSDを算出するPSD算出手段、応答PSDを目標PSDと比較して制御用PSDを求めるPSD比較手段、制御用PSDの各周波数成分に位相を与えて逆フーリエ変換するとともに、窓関数を掛けることにより、非ガウス特性を有するフレーム波形を生成するフレーム波形生成手段、複数のフレーム波形をシフトして重ね合わせることにより生成した非ガウス特性を有するドライブ波形を出力する重ね合わせ波形生成手段、を備えたこと、を特徴とする。これにより、複数のフレーム波形を重ね合わせた非ガウス特性を有する連続したドライブ波形を容易に生成することができる。1.第1の実施形態非ガウス特性のオープンループ図1に本発明の振動制御装置100のブロック図を示す。図1に示すPSD算出手段B01は、振動試験機2から得た応答波形をフーリエ変換して、応答PSDを算出する手段である。図1に示すPSD比較手段B03は、当該応答PSDを予め設定された目標PSDと比較して、制御用PSDを出力する手段である。なお、応答波形が存在しない起動時には、目標PSDをそのまま制御用PSDとして出力する。さらに、図1に示す非ガウス特性波形生成手段B05は、位相設定手段B07から与えられる位相を各周波数成分に与えて逆フーリエ変換し、これに窓関数を掛けることによって、1フレームのランダム波形を生成する手段である。図1に示す中心値設定手段B09は、与える位相の確率密度分布の中心値μを位相設定手段B07に与える。また、非ガウス特性制御手段B13は、予め設定された目標非ガウス特性[K’、S’]から、以下に説明する標準偏差σおよび初期位相φ0を求めて、位相設定手段B07に与える手段である。位相設定手段B07は、まず、最初のライン周波数成分には、非ガウス特性制御手段B13から受けた初期位相φ00~2πを設定し、隣接する各ライン周波数成分には、正規分布[中心値μ、標準偏差σ]に従うランダムな位相差ΔφR0~2πを順次加算して位相を設定する。ここで、重要な点は、各位相φkを一様なランダムで与えたのではないという点である。このように、正規分布に従ったランダムな位相差ΔφRを順に加算して各ラインの位相を算出することで、結果として、非ガウス特性のランダム波形を生成できる。図1に示す重ね合わせ手段B11は、非ガウス特性波形生成手段B05から非ガウス特性ランダム波形を複数受けて、これらをシフトして足し合わせることにより、連続した非ガウス特性ランダム波形を生成し、ドライブ波形として出力する手段である。非ガウス特性制御手段B13は、目標非ガウス特性[K’、S’]を読み出して、目標クルトシスK’に対応する制御用の標準偏差σと、目標スキューネスS’に対応する制御用の初期位相φ0の値を算出する手段である。以上の各手段を備えることで、本発明の振動制御装置100は、目標非ガウス特性を有する連続した非ガウス特性のランダム波形を生成し、振動試験機2に与えることができる。[ハードウェア構成]図1の振動制御装置100のハードウエア構成を、図2に示す。振動制御装置100には、A/D変換器8、D/A変換器10、メモリ12、フラッシュROM14、DSP16が接続されている。振動試験機2によって振動させられた供試体4の振動は、加速度センサ6によって検出され、アナログの応答波形としてA/D変換器8に出力される。A/D変換器8は、アナログの応答波形をデジタルデータに変換する。D/A変換器10は、デジタルのドライブ波形をアナログ信号に変換して、振動試験機2に与える。フラッシュROM14には、非ガウス波形制御のための制御プログラム18が記録されている。制御プログラム18は、目標PSDのデータのほか、制御しようとする目標非ガウス特性[目標クルトシスK’、目標スキューネスS’]の値や、これら目標非ガウス特性[目標クルトシスK’、目標スキューネスS’]に対応する標準偏差σ、初期位相φ0を決定するための計算式テーブルを有する図示せず。図3Aに、目標PSDとして予め設定された、一般的な輸送環境主として道路に対応するPSDデータの例を示すJISZ0232。図3Bに、目標クルトシスK’に対応する標準偏差σを算出するために用いられる、目標クルトシスK’と、位相差Δφの確率密度分布の標準偏差σとの関係を示す。図3Bのグラフを見ると、標準偏差σの値を大きくすれば、確率密度分布の急峻度尖度を示すクルトシスKは3に近づいていくこと、すなわち、ガウス特性のランダム振動が得られることがわかる。一方で、標準偏差σの値を小さくしていくと、クルトシスKの値が大きくなっていくこと、非ガウス特性のランダム振動が得られることがわかる。このようなクルトシスKと標準偏差σの関係に基づき、所望のクルトシスKに対応する標準偏差σの値を設定することが可能となる。分布の非対称性歪度を示すスキューネスSには、初期位相φ0との相関性が存在することが知られている。しかし、クルトシスKと標準偏差σの関係のように一義的に定まるものではなく、クルトシスKとの関係を考慮して探索して決定することが可能である。よって、例えば、予めクルトシスKおよび初期位相φ0との関係で得られる目標スキューネスSの関係を実測してテーブルに記憶しておくことで、目標スキューネスS’に対応する初期位相φ0の値を得ることができる。また、初期位相φ0を設計変数にして、K、Sを目的関数として予めアルゴリズムを作成しておくことで、スキューネスSを導き出すことも可能である。[制御プログラム18のフローチャート]図2に示す制御プログラム18により実行される処理について、図4のフローチャートを用いて以下に説明する。まず、DSP16図2は、振動試験機2からの応答波形を取得するステップS100。これを応答波形フーリエ変換して、各周波数成分についての応答PSDを算出し、メモリ12に記憶するステップS102。さらに、予め設定された目標PSDをフラッシュROM14から読み出しステップS104、応答PSDを目標PSDと比較して、応答PSDが目標PSDと等しくなるように制御用PSDを生成するステップS106。例えば、応答PSDと目標PSDとで、誤差スペクトルを求め制御用PSDを調整する処理を実行する。DSP16は、ステップS106で生成した制御用PSDの各周波数成分に与える位相を算出しステップS108、逆フーリエ変換処理を実行するステップS110。図5に位相設定処理ステップS108の詳細なフローチャートを示す。図5に示すように、DSP16は、フラッシュROM14に記憶された目標非ガウス特性[目標クルトシスK’、目標スキューネスS’]の値を読み出すステップS200。例えば、K’=4、S’=0の値が読み出される。つぎに、目標非ガウス特性[目標クルトシスK’、目標スキューネスS’]に対応する標準偏差σ、初期位相φ0を設定するステップS202。例えば、目標クルトシスK’=4のとき、図3Bに示すグラフから、対応する標準偏差σを220に設定する。目標非ガウス特性Sに対応する初期位相φ0は、目標スキューネスS’および目標クルトシスK’の値に基づいてテーブルを参照して設定することができる。さらに、DSP16は、正規分布の中心値μを設定するステップS204。例えば、中心値μの値として、フレーム毎に0~2πの間で一様にランダムな値を設定することができる。これにより、各フレームのピーク位置をランダムにすることが可能であり、クルトシスの劣化を抑制することができる。DSP16は、ステップS202で設定した初期位相φ0を、フレーム波形生成手段B05に与えるステップS206。さらに、図6に示すような正規分布[標準偏差σ、中心値μ]に従うランダムな位相差ΔφRを各位相ごとに発生させて、これを順次加えた位相φkを算出してフレーム波形生成手段B05に与えるステップS208。このように、最初の周波数成分に隣接する各周波数成分については、k番目の位相φk1を、φk1=φk2ΔφRより算出する。すなわち、第2ラインの周波数成分PSD1に対しては、正規分布に従った0~2πまでのランダムな位相差ΔφR1を、ステップS202で設定した初期位相φ0に加算して得た位相φ1=φ0ΔφR1が与えられる。第3ラインの周波数成分PSD2に対しては、正規分布に従ったランダムな位相差ΔφR2を位相φ1に加算して得た位相φ2=φ1ΔφR2が与えられる。以降のラインについても同様にして、全ての隣接するラインに位相差ΔφRを順次加算して得た位相が与えられる。上記のように与えた位相を設定して、図4のステップS110に示す逆フーリエ変換処理を実行した後、DSP16は、生成した非ガウス特性ランダム波形に窓関数を掛け合わせて窓操作処理を実行する図4のステップS112。なお、この実施形態では、1つの制御用PSDに対して、生成した異なる位相を与えて、複数フレームの非ガウス特性ランダム波形を生成するようにしている。ここで、窓関数を用いるのは、各フレーム間の非ガウス特性ランダム信号を滑らかに繋ぐためである。この実施形態では、窓関数としてハニングウインドウを用いている。窓操作処理を実行した所定フレーム数の非ガウス特性ランダム波形は、シフトした上でステップS114、順に重ね合わせ処理が行われるステップS116。例えば、図7に示すように、所定長だけずらされた状態で、波形W1~Wnが足し合わされる。これにより、図7に示す連続した波形Wが生成される。なお、制御用PSD図4のステップS106で生成は、離散値として算出された線スペクトルであるから、これに基づいて生成される非ガウス特性ランダム波形は、線スペクトルとしての周波数特性を有することになる。しかし、この実施形態では、異なる波形の非ガウス特性ランダム波形を複数生成し、窓操作、重ね合わせを行うことで線スペクトルを連続スペクトル化している。つまり、1つの制御用PSDに対して、生成した異なる位相を与えて、複数フレームの非ガウス特性ランダム波形を生成し、これらに窓操作処理ステップS112を実行し、さらに、シフト処理ステップS114および重ね合わせ処理ステップS116を施すことにより、線スペクトルから連続スペクトルに変えることができる。DSP16は、所定数のフレームを重ね合わせた連続した非ガウス特性ランダム波形をドライブ波形として出力する図4のステップS120。これにより、振動試験機2が、入力されたドライブ波形に基づいて、非ガウス特性の振動を発生させることができる。2.第2の実施形態非ガウス特性のフィードバック上記のようなオープンループ制御の場合、生成した非ガウス特性ランダム波形に比べて、応答波形の非ガウス特性が劣化してしまうことがある例えば、目標のクルトシスKが5であったとき、出力のクルトシスKが4に劣化してしまう。このような場合、目標クルトシスK’よりも大きなクルトシスを与えることが必要となる。そこで、応答波形から得られる非ガウス特性[応答クルトシスKとスキューネスS]をフィードバック制御するようにしてもよい。応答波形から得られる振動加速度の確率密度分布から得られる応答クルトシスK、応答スキューネスSをフィードバック制御した振動制御装置200のブロック図を図8に示す。図8に示すブロック図は、非ガウス特性算出手段B15が追加されている点および非ガウス特性制御手段B13が応答非ガウス特性と、目標非ガウス特性とを比較して、制御用の非ガウス特性を決定する点で、図1に示すブロック図とは異なっている。図8に示す非ガウス特性算出手段B15は、応答波形から非ガウス特性[応答クルトシスKとスキューネスS]を算出して非ガウス特性制御手段13に与える手段である。制御プログラム18が実行する処理は、図4のフローチャートに示す処理と同じであるが、位相設定処理図4のステップS108の詳細が異なる。図9に本実施形態における位相設定処理ステップS108の詳細なフローチャートを示す。図9に示すように、DSP16は、まず、応答波形から応答非ガウス特性[応答クルトシスKとスキューネスS]を算出するステップS300。応答非ガウス特性の算出は、応答波形から求めて行ってもよいし、ヒストグラム分析を行い確率密度分布から求めて行ってもよい。例えば、応答クルトシスK’=4、目標スキューネスS’=0の値が読み出される。つぎに、DSP16は、フラッシュROM14に記憶された目標非ガウス特性[目標クルトシスK’、目標スキューネスS’]の値を読み出すステップS302。例えば、K’=6、S’=0の値が読み出される。さらに、DSP16は、応答非ガウス特性と、目標非ガウス特性とを比較して、制御用の非ガウス特性を決定するステップS304。例えば、応答クルトシスK=4で、目標クルトシスK’=6のとき、制御クルトシスKは、劣化した分を目標クルトシスK’に加算して、K=8に設定する。つぎに、制御非ガウス特性[制御クルトシスK’、制御スキューネスS’]に対応する標準偏差σ、初期位相φ0を設定するステップS306。さらに、DSP16は、正規分布の中心値μを設定するステップS308。例えば、中心値μの値として、フレーム毎に0~2πの間で一様にランダムな値を設定することができる。DSP16は、ステップS202で設定した初期位相φ0を、フレーム波形生成手段B05に与えるステップS310。さらに、図6に示すような正規分布[標準偏差σ、中心値μ]に従うランダムな位相差ΔφRを各位相ごとに発生させて、これを順次加えた位相φkを算出してフレーム波形生成手段B05に与えるステップS312。このように、最初の周波数成分に隣接する各周波数成分については、k番目の位相φk1を、φk1=φk2ΔφRより算出する。上記のように設定した位相を与えて逆フーリエ変換処理図4のステップS110、窓操作処理図4のステップS112を実行して、応答非ガウス特性をフィードバック制御した非ガウス特性波形を生成することができる。3.第3の実施形態ZMNL方式での非ガウス特性のフィードバックなお、上記実施形態では、非ガウス特性波形生成手段B05図1、図8を設けることにより、非ガウス特性ランダム波形を生成したが、図10に示す振動制御装置300のように、代わりにガウス特性波形生成手段B06を設け、ガウス特性ランダム波形を重ね合わせて、連続したガウス特性波形を生成した後で、非ガウス化変換手段B25を適用して、非ガウス波形を生成してもよい。この実施形態では、ガウス特性波形を非ガウス特性波形に変換するために、図11に示すようなZMNL関数による変換を行っている。ZMNL関数は、入力信号の振幅を、当該振幅の大きさによって非線形的に変換して出力するものである。これにより、ガウス特性波形を非ガウス特性波形に変換することができる。ただし、ZMNL変換を行うと波形のPSDが変化してしまう。この点については、図10に示すように、ZMNL変換がPSDを制御するフィードバックループの中に設けられているので、フィードバックループによって、このPSD変化を吸収することができる。つまり、ZMNL変換を挿入したとしても、目標PSDによる振動を達成することができる。なお、上記のようにZMNL変換は非線形変換であるため、特定の周波数領域たとえば高周波においては、フィードバック制御では補えない周波数変化をもたらす可能性もある。ただし、ほとんどの周波数帯域では、フィードバックによる制御が実効性を持ち、ほぼ目標PSDを達成できるといってよい。また、実施形態1、2において示した位相差をランダムに与える手法では、PSDに変化はない制御用PSDを有する非ガウス特性ランダム波形を得ることができる。この点において、実施形態1、2の手法の方が優れている。一方で、ZMNL関数による変換においては、非ガウス特性の制御が容易であるという特徴がある。図10に示すガウス特性波形生成手段B06は、フーリエ変換時に各周波数成分にランダムな位相φkを与える等により、ガウス特性波形を生成する手段であり、生成された波形は、重ね合わせ処理の後、非ガウス特性波形に変換される。図10に非ガウス特性制御手段B15は、応答非ガウス特性と、目標非ガウス特性とを比較して、決定した制御用の非ガウス特性[応答クルトシスK、応答スキューネスS]を非ガウス化変換手段B25に直接与えて、非ガウス特性をフィードバック制御している。非ガウス化変換手段B25としては、例えば、図11に示すZMNL関数を用いることができる。ZMNL関数は、周波数特性が変更してしまう欠点がある反面、クルトシスとスキューネスの値を直接与えることで非ガウス特性の波形を容易に生成できる点で、非ガウス特性の制御がしやすいという利点がある。また、クルトシスとスキューネスの値を算出しなくても、応答波形から得られる確率密度分布をそのままフィードバック制御に利用することができるという利点もある。そこで、この実施形態では、図11に示すようなZMNL関数Zeromemorynonlinersysytem関数を用いて、所望の非ガウス特性を有する信号に変換するようにしているステップS26。図11において、横軸は与えられた信号であり、縦軸は出力信号である。ガウス特性ドライブ波形xを、このZMNL関数によって変換すると、非ガウス特性ドライブ波形yが得られる。ZMNL関数を数式で表すと、数1のとおりである。xがガウス特性信号、yが得られる非ガウス特性信号である。各係数の詳細は、数2に示すとおりである。上式におけるスキューネスS、クルトシスKの値を代入することにより、変換によって所望の非ガウス特性ドライブ波形を生成することのできるZMNL関数を得ることができる。なお、上式では、2.8より小さい目標クルトシスを設定することができない。そこで、数3のような関数を用いるようにしてもよい。係数Cnは、数4に示す誤差変数εを最小にする最適化問題を解くことにより決定することができる。以上のようにして得られた非ガウス特性ドライブ波形を与えることで、非ガウス特性を容易にフィードバックすることができる。4.伝達特性に基づく波形制御第2、第3の実施形態では、生成した非ガウス特性ランダム波形を振動試験機2に与えても応答波形の非ガウス特性が劣化することに鑑み、非ガウス特性のフィードバック制御を行っている。しかし、生成した非ガウス特性ランダム波形と同じ振動を供試体4に与えることができれば、上記の問題は解決する。そこで、この実施形態では、以下のように伝達特性に基づいて波形制御することで、非ガウス特性ランダム波形と同じ振動を与えることとした。図12は、伝達特性に基づいてドライブ波形を修正する振動制御装置400のブロック図である。変換手段B23は、ドライブ波形生成手段32および制御特性更新手段34を備える。図12に示すドライブ波形生成手段32は、振動試験機2および供試体を含む系の伝達特性の逆特性を制御特性として、非ガウス特性ドライブ波形を変形してドライブ波形を生成する。ドライブ波形は、D/A変換器10によってアナログ信号に変換され、振動試験機2に与えられる。図12に示す制御特性更新手段34は、ドライブ波形と応答波形に基づいて、系の伝達特性を算出し、その逆特性である制御特性を更新する。更新された制御特性は、ドライブ波形生成手段32に与えられる。これにより、供試体4に対し、変換手段28の変換特性を設定することにより目標とする非ガウス特性を有し、かつ、目標スペクトルを有する振動を、供試体4に与えることができる。5.その他の実施形態なお、上記実施形態図1、図8では、目標非ガウス特性の値を予め記憶することとしたが、これに限定されるものではない。例えば、目標非ガウス特性として、振動加速度の確率分布グラフを予め記憶しておき、これから目標非ガウス特性の値を算出してもよい。なお、上記実施形態図8では、非ガウス特性算出手段B15が応答波形から算出した応答非ガウス特性の値と、予め記憶された目標非ガウス特性の値とを比較することとしたが、これに限定されるものではない。例えば、図8に示す非ガウス特性算出手段B15を設けずに、非ガウス特性制御手段B13が、応答波形から振動加速度の確率密度分布を生成し、これを目標非ガウス特性として予め記憶した振動加速度の確率密度分布と比較するようにしてもよい。この場合、応答波形から非ガウス特性スキューネス、クルトシスの値を算出する必要がない。なお、上記実施形態では、図1に示す中心値設定手段B09は、正規分布の中心値μを0~2πの間で一様にランダムな値に設定したが図5のステップS204および図9のステップS308、これに限定されるものではない。例えば、正規分布の中心値μを固定値に設定してもよい。なお、上記実施形態では、隣接する周波数成分に加算するランダムな位相差ΔφRを0~2πの間に収まるように算出したが図5のステップS208、図9のステップS312、これに限定されるものではない。例えば、次式ΔφR=2πPσRNDにより、ランダムな位相差ΔφRを算出するようにしてもよい。ここで、Pはフレームにおけるピークの出現位置0~1、σは標準偏差、RNDは正規乱数0~1を示す。なお、図10に示す第3の実施形態では、応答波形の非ガウス特性をフィードバック制御しているが、応答波形の非ガウス特性をフィードバック制御しないようにしてもよい。なお、図12に示す第4の実施形態では、応答波形の非ガウス特性をフィードバック制御していないが、応答波形の非ガウス特性をフィードバック制御するようにしてもよい。なお、上記実施形態では、予め設定した目標非ガウス特性に近づくように非ガウス波形を生成したが、目標非ガウス特性を設定しないて非ガウス波形を生成してもよい。なお、上記実施形態では標準偏差σに対応する目標クルトシスを設定したが、分散σ2に対応する目標クルトシスを設定してもよい。なお、上記実施形態では、DSP16によって処理を実行したが、代わりにCPUを用いて処理を実行してもよい。なお、上記実施形態では、1つの制御用PSDから複数フレームの非ガウス性ランダム波形を生成する際に、フレーム毎に異なった位相を与えるようにしている、すなわち、制御用PSDの各ラインの強度をP1P2P3P4Pnとしたとき、以下のようにφ1φ2φ3φ4φmを与えて、逆フーリエ変換を行っている。1フレーム目P1、P2、P3...Pn2フレーム目P1、P2、P3...Pn3フレーム目P1、P2、P3...Pnmフレーム目P11、P22、P33...Pnn上記のように、フレーム毎に異なった位相を与えることにより、窓関数をかけて重ね合わされた波形の、PSDは線スペクトルではなく連続スペクトルになる。しかし、複数のフレームに同一の位相を与えて、複数フレームの非ガウス性ランダム波形を生成し、窓関数をかけて重ね合わせることにより、PSDとして線スペクトルを有する波形を得ることもできる。すなわち、制御用PSDの各ラインの強度をP1P2P3P4Pnとしたとき、以下のようにφ1φ2φ3φ4φmを与えて、逆フーリエ変換を行っている。1フレーム目P1、P2、P3...Pn2フレーム目P1、P2、P3...Pn3フレーム目P1、P2、P3...Pnmフレーム目P1、P2、P3...Pn換言すると、複数フレームについて同一の非ガウス性ランダム波形を生成し、これに窓関数をかけて重ね合わせた場合には、線スペクトルを有する波形を得ることができる。通常は、連続スペクトルを有する波形の方が好ましいが、試験の目的によっては線スペクトルを用いることもできる。なお、上記実施形態では、異なる位相を与えることにより、1つの制御用PSDから複数フレームの非ガウス性ランダム波形を生成することとしたが、これに限定されるものではない。例えば、図15に示すように、開始位置をランダムにして、波形をずらせることにより、1つの制御用PSDから複数フレームの非ガウス性ランダム波形を生成してもよい。具体的には、1つの制御用PSDから、図15Aに示す1フレームのランダム波形w1を生成する。次に、開始位置をそれぞれランダムな位置T2、T3、T4図15B、C、Dにずらせて、このランダム波形W1を循環的に読み出す。すなわち、1フレームの終わりまで読み出すと、先頭に戻って読み出しを続ける。これにより、図15B、C、Dに示す複数フレームの非ガウス性ランダム波形W2、W3、W4を簡単に生成することができる。