(57)【特許請求の範囲】
【請求項１】
指定振幅確率分布と指定交差分布から計算された互いに異なる推移確率をそれぞれ格納する複数のメモリと、
これら各メモリにそれぞれ格納された推移確率を用いて前記指定振幅確率分布と前記指定交差率分布に従う信号をそれぞれ生成する複数の任意分布乱数生成部と、
これら任意分布乱数生成部がそれぞれ生成した前記信号の一つを選択して出力するスイッチ回路と、
このスイッチ回路から出力された前記信号に所定の遅延を与えて前記スイッチの選択制御をする遅延回路と、
前記スイッチにより選択されて出力された信号をアナログ値に変換するＤ／Ａコンバータと
を備え、
前記指定振幅確率分布と前記指定交差率分布に従う擬似雑音を生成する擬似雑音発生装置。
【請求項２】
指定振幅確率分布と指定交差分布から推移確率を示す複数のディジタル情報を計算して出力する計算機と、
前記複数のディジタル情報をそれぞれ格納する複数のメモリと、
該複数のメモリにそれぞれ格納された前記複数のディジタル情報を用いて、該複数のディジタル情報にそれぞれ対応して前記指定振幅確率分布と前記指定交差率分布に従う複数のディジタル出力信号を生成する複数の任意分布乱数生成部と、
これら任意分布乱数生成部がそれぞれ生成した前記ディジタル出力信号の一つを選択して出力するスイッチ回路と、
このスイッチ回路から出力された前記ディジタル出力信号に所定の遅延を与えて前記スイッチの選択制御をする遅延回路と、
前記スイッチ回路から出力された前記ディジタル出力信号をアナログ値に変換するＤ／Ａコンバータと
を備え、
前記指定振幅確率分布と前記指定交差率分布に従う擬似雑音を生成する擬似雑音発生装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、デジタル通信システムの通信品質の評価、電磁妨害波に対する電子機器の耐妨害性（Immunity）の評価を行うため等に用いる擬似雑音を発生する擬似雑音発生装置に関するものである。
【０００２】
【従来の技術】
多くのノイズが加算された場合や熱雑音、都市雑音の場合などでは、ノイズの振幅は正規分布となる。これらをシミュレートするためのノイズ・ダイオードを利用しホワイトノイズを発生させるものがある。
狭帯域デジタル通信では、通信システムのビット誤り率（ＢＥＲ：Bit Error Rate）と電磁妨害波の振幅確率分布（ＡＰＤ：Amplitude Probability Distribution）との間に相関関係があり、電磁妨害波の振幅確率分布（ＡＰＤ）から通信システムのＢＥＲを推定することができるという報告がされている。この任意分布乱数発生器を用いて、指定された振幅確率分布（ＡＰＤ）に従うノイズを発生させるものなどがある。
【０００３】
【発明が解決しようとする課題】
ノイズの性質を表す際の重要な統計パラメータとして振幅確率分布（ＡＰＤ），交差率分布（ＣＲＤ：Crosing Rate Distribution ）などがある。
振幅確率分布は、観測した電磁妨害波等の振幅が観測した振幅レベルを超えている時間率で定義される。交差率分布は、単位時間内で電磁妨害波が観測振幅レベルを正交差する回数で定義される。
ホワイトノイズを発生させる擬似雑音発生器は、ノイズの分散と平均を指定することはできるが、ノイズの分布は正規分布に限定されており、振幅確率分布を指定することはできない。
また、任意分布乱数発生器を用いたものは、ノイズの振幅確率分布を指定し、任意の振幅確率分布を持つノイズを発生することはできるが、時間相関のない独立事象を対象としている。
【０００４】
しかし、本発明装置により擬似しようとする雑音対象に含まれる電子レンジや一般の電子機器からのノイズは、電源電圧の周期や動作クロックの周期に依存した非独立事象であり、従来の任意分布乱数発生器を使用した擬似雑音発生器から発生されるような、交差率分布，パルス幅分布，パルス間隔分布などが時間相関を有する非独立事象のノイズとは異なっている。
【０００５】
本発明の目的は、ノイズの統計パラメータのうち、振幅確率分布と交差率分布を指定した擬似雑音を発生することができる擬似雑音発生装置を提供することである。
【０００６】
【課題を解決するための手段】
この課題を解決するため、本発明による擬似雑音発生装置は、指定振幅確率分布と指定交差分布から計算された推移確率をそれぞれ格納する複数のメモリと、これら各メモリにそれぞれ格納された推移確率を用いて前記指定振幅確率分布と前記指定交差率分布に従う信号をそれぞれ生成する複数の任意分布乱数生成部と、これら任意分布乱数生成部がそれぞれ生成した前記信号の一つを選択して出力するスイッチ回路と、このスイッチ回路から出力された前記信号に所定の遅延を与えて前記スイッチの選択制御をする遅延回路と、前記スイッチにより選択されて出力された信号をアナログ値に変換するＤ／Ａコンバータとを備え、前記指定振幅確率分布と前記指定交差率分布に従う擬似雑音を生成するように構成されている。
【０００７】
また、本発明による擬似雑音発生装置は、指定振幅確率分布と指定交差分布から推移確率を示す複数のディジタル情報を計算して出力する計算機と、前記複数のディジタル情報をそれぞれ格納する複数のメモリと、該複数のメモリにそれぞれ格納された前記複数のディジタル情報を用いて、該複数のディジタル情報にそれぞれ対応して前記指定振幅確率分布と前記指定交差率分布に従う複数のディジタル出力信号を生成する複数の任意分布乱数生成部と、これら任意分布乱数生成部がそれぞれ生成した前記ディジタル出力信号の一つを選択して出力するスイッチ回路と、このスイッチ回路から出力された前記ディジタル出力信号に所定の遅延を与えて前記スイッチの選択制御をする遅延回路と、前記スイッチ回路から出力された前記ディジタル出力信号をアナログ値に変換するＤ／Ａコンバータとを備え、前記指定振幅確率分布と前記指定交差率分布に従う擬似雑音を生成するように構成することができる。
【０００８】
【発明の原理】
図１のような波形の振幅レベルｋにおける交差率分布ｃk （０≦ｋ≦ｎ−２）、すなわち時刻ｔ−１に振幅レベルｉ（０≦ｉ≦ｋ）の範囲に存在する信号が時刻ｔに振幅レベルｊ（ｋ＜ｊ≦ｎ−１）の範囲に移動する確率は、式（１）で表される。
【０００９】
【数１】
0003978580.tif 000002
【００１０】
ここで、ｔは離散形時間（ｔ＝０，１，…）、ｐi は振幅レベルｉにおける信号振幅の確率密度である。また、ｗijは信号が振幅レベルｉから振幅レベルｊに移動する確率（推移確率）であり、逆に振幅レベルｊから振幅レベルｉに移動する確率はｗjiで表される。
【００１１】
図２で示すと、推移確率ｗijは図示のようなｎ×ｎ行列の対角線から右上の部分、推移確率ｗjiは対角線の左下の部分を表す。
ここで、推移確率ｗij，ｗjiと確率密度ｐi ，ｐj の間に式（２）が成立すれば、推移確率ｗxy（０≦ｘ，ｙ≦ｎ−１）の極限推移確率はｐy （０≦ｙ≦ｎ−１）と一致する。すわなち、この推移確率ｗxyを用いて擬似雑音を生成すると、擬似雑音の確率密度はｐy となる。また、推移確率ｗxyの対角線より右上の成分ｗijに式（１）が成立すれば、推移確率ｗxyを用いて生成した擬似雑音の交差率分布は振幅レベルｋにおける交差率分布ｃk と一致する。従って、推移確率ｗxyが式（１）と（２）を同時に満足すれば、推移確率ｗxyから確率密度関数ｐy と交差率分布ｃk に従う擬似雑音を生成可能である。
【００１２】
【数２】
0003978580.tif 000003
【００１３】
推移確率ｗxyについて更に検討する。図３において、点線で囲む各行のグループは、時刻ｔ−１における振幅レベルが、上段から順次０，１，ｘ，ｎ−１のときの確率密度に対応する推移確率分布ｗxyである。図３に示すごとく、点線で囲まれた推移確率ｗxyのｘ行成分ｗxy（０≦ｙ≦ｎ−１）は、時刻ｔ−１における信号の振幅レベルがｘであるときに、時刻ｔに発生する信号の振幅レベルｙ（０≦ｙ≦ｎ−１）の確率密度を表す。この確率密度をｐｘy （０≦ｙ≦ｎ−１）とすると、推移確率ｗxy（０≦ｙ≦ｎ−１）のｘ行成分は式（３）で表される。
【００１４】
【数３】
0003978580.tif 000004
【００１５】
また、時刻ｔ−１における信号の振幅レベルがｘであるときの確率密度ｐｘy （０≦ｙ≦ｎ−１）から、時刻ｔ−１における信号の振幅レベルがｘであるときの振幅確率分布ａｘy （０≦ｙ≦ｎ−１）は、式（４）を用いて計算することができる。振幅確率分布ａｘy に従う擬似雑音は、従来から存在する任意分布乱数発生器で生成可能である。
【００１６】
【数４】
0003978580.tif 000005
【００１７】
従って、振幅レベルがｎ個に分割されている場合（時刻ｔ−１において取り得る信号振幅レベルｘが０≦ｘ≦ｎ−１のｎ個の場合）、図４に示すごとくｘ行の振幅確率分布ａｘy を計算することができる。すなわち、図４は推移確率ｗxyから振幅確率分布ａｘy への変換を説明するための図であり、推移確率ｗxy（０≦ｙ≦ｎ−１）について、上段は、推移確率ｗ0yを確率密度ｐ０y とし振幅確率分布ａ０y を計算し、中段は、推移確率ｗxyを確率密度ｐｘy とし振幅確率分布ａｘy を計算し、下段は、推移確率ｗn-1,y を確率密度ｐ（ｎ−１）y とし振幅確率分布ａ（ｎ−１）y を計算する。
【００１８】
また、図５に示すごとく、ｎ個の任意分布乱数発生器ｘ（０≦ｘ≦ｎ−１）を用意しておき、任意分布乱数発生器Ｇx において時刻ｔ−１における信号の振幅レベルがｘであったときの振幅確率分布ａｘy に従う信号ｙx を発生させる。ここで、Ｇ0 はａ０y （０≦ｙ≦ｎ−１）に従う任意分布乱数発生器、Ｇx はａｘy （０≦ｙ≦ｎ−１）に従う任意分布乱数発生器、Ｇn-1 はａ（ｎ−１）y （０≦ｘ≦ｎ−１）である。さらに、時刻ｔ−１における信号の振幅レベルｘに対応して任意分布乱数発生器Ｇx の出力ｙ0 ，…，ｙx ，…ｙn-1 を選択することにより、ｐy およびｃk に従う擬似雑音ｙを生成することが可能である。
【００１９】
【実施例】
図６に、以上の原理に基づく、本発明による擬似雑音発生装置の構成例を示す。この擬似雑音発生装置は、計算機１と任意分布乱数発生部２、Ｄ／Ａ変換部３から構成されている。計算機１は、例えば、パーソナルコンピュータである。任意分布乱数発生部２は、パーソナルコンピュータ１からのディジタル情報ｄ0 ，ｄ1 ，…，ｄx ，…，ｄn-1 をそれぞれ保持するメモリＭ0 ，Ｍ1 ，…，Ｍx ，…，Ｍn-1 と、各メモリＭ0 ，Ｍ1 ，…，Ｍx ，…，Ｍn-1 に対応する任意乱数発生器Ｇ0 ，Ｇ1 ，…，Ｇx ，…，Ｇn-1 より構成されている。Ｄ／Ａ変換部３は、任意乱数発生器Ｇ0 ，Ｇ1 ，…，Ｇx ，…，Ｇn-1 の各出力ｙ0 ，ｙ1 ，…，ｙx ，…，ｙn-1 を選択するスイッチ３−１と、スイッチ３−１の出力ｙに所要の遅延を与えてその出力ｘによりスイッチ３−１の選択制御をする遅延回路と、スイッチ３−１の出力ｙをＤＡ変換するＤ／Ａ変換器３−３とを備えている。
【００２０】
この擬似雑音発生装置の動作について図７を参照して説明する。
最初に、計算機１上で確率密度ｐy （０≦ｙ≦ｎ−１）と交差率分布ｃk （０≦ｋ≦ｎ−２）から推移確率ｗxy（０≦ｘ，ｙ≦ｎ−１）を求め、推移確率ｗxyのｘ行の成分ｐｘy （０≦ｙ≦ｎ−１）から任意分布乱数発生器Ｇx に設定するデータｄｘy （０≦ｙ≦ｎ−１）を作成する。
確率密度関数ｐy と交差率分布ｃk から状態推移確率ｗxyを決定するために、最適化を用いる。すなわち、適当にｗxyの初期値を定め、式（１），（２）を満たすように少しずつｗxyを修正する。以下にｗxyの決定方法を示す。
【００２１】
まず、ｗxyの初期値を式（５）のように決定する（Ｓ１）。ここで、ｐy は擬似雑音の振幅の確率密度である。
【００２２】
【数５】
0003978580.tif 000006
【００２３】
次に、式（６）を用い、現在の推移確率ｗijを用いたレベルｋでの交差率分布ｓk （０≦ｋ≦ｎ−２）を求める（Ｓ２）。さらに、式（７）より、このｓk と指定交差率分布ｃk の誤差ｅk を求めてその比ｅk ／ｃk が許容値th以内かをテストし（Ｓ３）、そのテストがＹｅｓのときに、式（８）のようにｅk を用いｗijを修正する（Ｓ４）。以上（６）から（８）をｋ＝０，１，２，…，ｎ−２において計算する（Ｓ５）。
【００２４】
【数６】
0003978580.tif 000007
【００２５】
次に、式（２）を満たすようにｗjiを修正する。ｗjiの修正には、式（２）を変形した式（９）を用いる（Ｓ６）。
【００２６】
【数７】
0003978580.tif 000008
【００２７】
式（１０）を用いて修正後の推移確率ｗxyの各行の和 sumx を求めると、 sumx ≠１となるが、本来推移確率ｗxyの各行の和は１である。従って、式（１０）で和を計算し、推移確率ｗxyの和が１となるように、式（１１）で推移確率ｗxyを調整する（Ｓ７）。以上の計算処理を、式（１０），（１１）を全ての行（０≦ｘ≦ｎ−１）について実行する（Ｓ８）。
【００２８】
【数８】
0003978580.tif 000009
【００２９】
再度、全てのｋ（０≦ｋ≦ｎ−２）において、修正後のｅk を計算し、式（１２）を満足すれば、ｐy とｃk を満たすｗxyが求まったとして処理を終了する。ここで、thは誤差の許容値である。それ以外は、再び式（６）に戻りｗxyの修正を行う（Ｓ１０）。
【００３０】
【数９】
0003978580.tif 000010
【００３１】
推移確率ｗxyが求まった後、引続き計算機１上で任意分布乱数発生器Ｇx （０≦ｘ≦ｎ−１）に与えるデータｄｘy （０≦ｙ≦ｎ−１）を計算する。この推移確率ｗxy（０≦ｘ，ｙ≦ｎ−１）のｘ行成分を確率密度ｐｘy （０≦ｙ≦ｎ−１）とし、確率密度関数ｐｘy から振幅確率分布ａｘy （０≦ｙ≦ｎ−１）を計算し、ｄｘy に変換する（Ｓ１１）。
ｎ個の任意分布乱数発生器を用意しておき、計算機から任意分布乱数発生器ｘにデータｄｘy を設定する。
【００３２】
ｄｘy にもとづき任意分布乱数発生器Ｇx において振幅確率分布ａｘy に従う信号ｙx （０≦ｙx ≦ｎ−１）を発生させ、時刻ｔ−１における信号の振幅レベルｘ（（０≦ｘ≦ｎ−１）に応じて任意分布乱数発生器Ｇx の出力を選択することにより、ｐy およびｃk に従う擬似雑音ｙを生成する。さらに、Ｄ／Ａ変換し出力する（Ｓ１２）。
【００３３】
本発明による処理のシミュレーション結果を、図８乃至図１１に示す。これらの図では、実線が設定値、記号がシミュレーション結果を表す。図８，９は、従来からある振幅確率分布（確率密度の累積分布）のみを指定することができる任意分布乱数発生器から出力された擬似雑音のシミュレーション結果である。図８に示すごとく、確率密度関数は設定値と一致するが、図９に示すように、交差率分布は設定値をシミュレーション結果が異なっている。一方、本発明による処理を用いて発生させた擬似雑音のシミュレーション結果では、図１０に示す確率密度関数も、図１１に示す交差率分布も、設定値とシミュレーション結果が一致している。
【００３４】
【発明の効果】
以上の説明から明らかなように、本発明によれば、指定した振幅確率分布，交差率分布に従う擬似雑音を生成することができる。
【図面の簡単な説明】
【図１】本発明に用いる交差率を説明するための図である。
【図２】本発明に用いる推移確率ｗxyの全体図である。
【図３】本発明に用いる推移確率ｗxyのｘ行と時刻ｔ−１の振幅レベルｘに依存した確率密度ｐｘy の対応図である。
【図４】本発明に用いる推移確率ｗxyから振幅確率分布ａｘy への変換を示した図である。
【図５】本発明に用いる振幅確率分布ａｘy を用いた擬似雑音発生の概念図である。
【図６】本発明による擬似雑音発生装置構成図である。
【図７】本発明に用いる推移確率ｗxyおよび擬似雑音発生装置用設定データｄｘy を作成するフローチャートである。
【図８】本発明に用いる振幅確率分布のみを指定した擬似雑音の確率密度関数を示す図である。
【図９】本発明に用いる振幅確率分布のみを指定した擬似雑音の交差率分布を示す図である。
【図１０】本発明に用いる振幅確率分布と交差率分布を指定した擬似雑音の確率密度関数を示す図である。
【図１１】本発明に用いる振幅確率分布と交差率分布を指定した擬似雑音の交差率分布を示す図である。
【符号の説明】
１  計算機
２  任意分布乱数発生部
３  Ｄ／Ａ変換部
３−１  スイッチ
３−２  遅延回路
３−３  Ｄ／Ａ変換器
Ｇ0 ，Ｇ1 ，…Ｇx …Ｇn-1   任意分布乱数発生器
Ｍ0 ，Ｍ1 ，…Ｍx …Ｍn-1   メモリ
【図１】
0003978580.tif 000011
【図２】
0003978580.tif 000012
【図３】
0003978580.tif 000013
【図４】
0003978580.tif 000014
【図５】
0003978580.tif 000015
【図６】
0003978580.tif 000016
【図７】
0003978580.tif 000017
【図８】
0003978580.tif 000018
【図９】
0003978580.tif 000019
【図１０】
0003978580.tif 000020
【図１１】
0003978580.tif 000021