(57)【要約】【課題】通信システム及び通信装置に使用するための、特に規定の（規格化された）ノイズプロファイルに従って上記システム及び装置を試験するための改良された信号を提供することを目的とする。【解決手段】所定の品質を有する該信号を生成するための該方法は、それぞれがスペクトル振幅特性及び位相特性を有する複数の周波数成分を備える第１信号（１０）を表現するステップと、そのスペクトル振幅特性を配列する（１１）ことによって該表現された信号を処理するステップと、該所定の品質基準に従ってそのスペクトル振幅特性を配列する（１１）ことによって、該表現された信号を処理するステップとを含む。【選択図】図２
2008219864.tif 000002
【特許請求の範囲】
【請求項１】
所定のスペクトル振幅の包絡線を有する信号を生成する方法であって、前記方法は、
−それぞれがスペクトル振幅特性及び位相特性を有する、複数の周波数成分からなる第１信号を表現するステップと、
−前記所定のスペクトル振幅の包絡線に従って、前記スペクトル振幅特性を設定すると共に、ランダム位相特性を配列することにより、前記表現された第１信号を処理するステップと、
−振幅歪み関数を前記第１信号に適用し、高い振幅ピークを増幅するステップと、
−歪んだ前記第１信号を周波数整形するステップと、
−前記第１信号の振幅分布の試験波形整形により、前記第１信号が波高因子の要求を満足したことを示すまで、前記振幅歪み関数を適用するステップと周波数整形するステップとを繰り返すステップと、
を含むことを特徴とする方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、概して通信システムに係わり、詳細には、通信システムとともに使用するための信号、このような信号を生成する方法及びこのような信号を生成するためのシステム、このような信号を使用して通信システムの動作を試験する方法、このような方法を作動させるために装置される試験システム及び（電気）通信システムに係わる。
【背景技術】
【０００２】
とりわけ、ｘＤＳＬトランシーバ及びケーブルまたはネットワークなどの通信システム及び通信設備を試験するためには、ケーブルごとに大多数のシステムまたはシステム機器を用いて、実際の配備シナリオにとって典型的となるように、通信システム及び通信装置に応力をかけるための試験信号が必要とされる。
【０００３】
現実的な（騒々しい）試験条件下でシステムまたはシステム機器の伝送性能を測定することによって、人はシステムまたは装置の設計を改良する、及び／またはそれらの性能がＥＴＳＩ、ＩＴＵまたはＡＮＳＩあるいは他の（電気）通信団体によって発行されるような規格に準拠していることを証明することができる。
【０００４】
このような性能試験を実行する方法は、障害として知られている信号を生成することである。さらに具体的には、障害は以下に細かく分けることができる。
【０００５】
（ｉ）スペクトル包絡線及び例えば近隣のｘＤＳＬシステムからのスペクトル振幅分布により特徴付けられるノイズプロファイルを有する漏話雑音と、
（ｉｉ）数多くのディスクリート周波数成分とスペクトル振幅、例えば、無線及びアマチュア放送から生じる変調深度パラメータと変調幅パラメータによって特徴付けられるノイズプロファイルを有する、ｒｆｉ−トーンとも呼ばれる、ディスクリート周波数成分から構成される進入雑音と、
（ｉｉｉ）例えば、切り替え動作及び構成要素により引き起こされる、信号パルスにより特徴付けられるインパルス雑音。
進入雑音のケースでは、周波数は時間が経つと変化する（掃引する）場合がある。
【０００６】
障害を発生させるための装置は、障害ジェネレータして知られており、前記漏話雑音及び進入雑音の少なくとも１つを生成するために、少なくとも１つが特に通信システム内または通信システム上での使用のために装置される。
【０００７】
実際には、通信システム及び通信装置が規格に準拠しているかどうかを試験するために、とりわけ、例えば通信ケーブル内のワイヤ組の長さと数、及び伝送データレートなどのシステムパラメータに従って変化する多様なノイズプロファイルが明確にされてきた。
【０００８】
さらに、ケーブル、銅ケーブルまたは光ファイバまたは他のケーブル種別などの伝送媒体のそれぞれの異なる種類または長さが、異なる雑音信号を必要とする。
【０００９】
ノイズプロファイルを作成するための方法及び装置は、技術で既知である。特にフィルタリング技法及びフィルタが、特定のスペクトル包絡線及びスペクトル振幅分布を有する出力信号を提供する入力信号から雑音を生成するために既知である。
【００１０】
しかしながら、フィルタリング技法及びフィルタを使用することにより、入力信号と出力信号の間に因果関係が確立される。当業者は、このような種別の信号が本物の稼動中の通信システムと通信装置の現実的な模倣により適していないことを認識するだろう。
【００１１】
特許文献１は、振幅の所定ヒストグラムに近づく不規則時間領域信号を生成するための方法及び装置を開示する。第１ステップでは、信号は、白色雑音信号などの雑音信号をフィルタリングし、それにより規定のスペクトル包絡線を有する信号を生成することにより生じる。次のステップでは、振幅の規定ヒストグラムに近づく、必要な時間領域信号を生じさせるために、非線形関数が該フィルタリングされた雑音信号に適用される。追加のステップでは、そのスペクトル包絡線を補正し、必要とされるスペクトル包絡線を有する出力信号を得るために、パルス応答フィルタリングが該時間領域信号に適用される。非線形関数とパルス応答フィルタリング関数の両方とも、提供されるスペクトル包絡線に従って選択される特殊な関数である。
【００１２】
特許文献１は、規定のスペクトル包絡線を有する時間領域信号だけが提供されるにすぎないという意味で制限されている。ＷＯ第００／１６１８１号は、提供される時間領域信号に関して課される他の品質基準、とりわけ位相特性に関して記述していない。
【特許文献１】国際公開ＷＯ００／１６１８１号パンフレット
【発明の開示】
【発明が解決しようとする課題】
【００１３】
通信システム及び通信装置に使用するための、特に規定の（規格化された）ノイズプロファイルに従ってこのようなシステム及び装置を試験するための改良された信号を提供することが、本発明の目的である。
【課題を解決するための手段】
【００１４】
本発明の第１態様では、好ましくは、通信システム内または通信システム上で使用するために、所定の品質基準を有する信号を準備し配列する方法が開示され、該方法は、
−それぞれがスペクトル振幅特性及び位相特性を有する複数の周波数成分を備える第１信号を表現するステップと、
−該所定品質基準またはそれぞれの所定品質基準に従ってスペクトル振幅特性を配列し、無作為位相特性を配列することによって、該表現される第１信号を処理するステップと、
を備える。
【００１５】
スペクトルの包絡線を修正する従来の方法は、デジタルフィルタバンクを使用することである。本発明の目的のために、該表現された第１信号と提供される信号間の因果関係を確立する必要はないため、これは、理想的とは程遠い。本発明に従ったこの内容の理解が、周波数波形整形のアプローチを大幅に簡略化する。
無作為な位相特性を有する第１信号から開始し、該表現された第１信号の周波数波形整形は、例えば該規定品質基準を満たす信号を提供するための適切な動作である可能性がある。本発明による周波数波形整形は、複数の方法で実行できる。
【００１６】
本発明のある実施形態では、第１信号は、各周波数成分のスペクトル振幅及び位相を指定する数の第１集合で表現される。該信号の該無作為位相特性を維持しながら、該周波数領域内で該表現された信号の周波数波形整形を達成するには、各周波数成分の該スペクトル振幅の基準化で十分である。
【００１７】
本発明の追加の実施形態では、該第１信号は、実数部と虚数部分を有する複素数の第２集合により表現され、それらの部分は、組み合わされて各周波数成分のスペクトル振幅及び位相を指定する。周波数波形整形は、該表現された第１信号の基準化後に無作為位相特性を維持するためなどであるが、複素数を適切に基準化することにより達成される。
【００１８】
本発明のまだ追加の実施形態では、該第１信号は、該時間領域内の該第１信号の振幅をそれぞれが指定する、数の第３集合によって表現される。例えば高速フーリエ変換（ＦＦＴ）アルゴリズムを使用して、時間領域から周波数領域にこの数の第３集合を変換することにより、第１信号は、各周波数成分のスペクトル振幅及び位相を指定する数の第４集合によって表現される。この数の第４集合は、数の第１集合と関連して前記に開示されたように、周波数波形整形動作によりさらに処理できる。
【００１９】
ただし、数の該第３集合も、実数部と虚数部分を有する複素数の第５集合により第１信号を表現するために、本発明に従って、時間領域から周波数領域に変換されてよい。前記に開示されたように、周波数波形整形の目的のため、複素数の該第５集合は適切に基準化されなければならない。
【００２０】
無作為ではない位相特性を有する、表現される第１信号のケースでは、無作為位相特性は、該第２集合、該第４集合及び該第５集合を適切に配列することにより近づくことができる。
【００２１】
周波数領域内での基準化は、実数基準化係数または複素数基準化係数を使用して、乗算動作によって呼び出すことができる。ある周波数成分の該スペクトル振幅の乗算のための該基準化係数は、その所望される値をその実際のスペクトル振幅で除算することにより検出される。
【００２２】
該品質基準または所定の品質基準に対するより近い一致を達成するために、本発明による方法の追加実施形態に従って、該処理済みの表現された第１信号の事後処理が提供される。
【００２３】
しかしながら、本発明による通信システム内または通信システム上での使用のため、該周波数領域にこのようにして配列される該表現された第１信号は、例えば、高速フーリエ逆変換（ＩＦＦＴ）アルゴリズムを使用して、該時間領域に変換されなければならない。
【００２４】
さらに、前記に開示された処理ステップは、信号の畳み込みまたは逆重畳または乗算及び付加などの演算も含んでよい。該時間領域では、該所定品質基準またはそれぞれの所定品質基準を満たす該処理済みの表現された第１信号が、とりわけ、該時間領域内の数の第６集合によって表現されてよい。
【００２５】
しかしながら、前記アプローチでは、スペクトル振幅特性及び無作為位相特性などの該周波数領域内の品質基準を満たす、提供された該信号が、まだ、所定の時間領域振幅分布などの該時間領域での品質基準を満たさないことがある。
【００２６】
本発明による該方法のまだ追加の実施形態では、該所定品質基準またはそれぞれの所定品質基準は、所定の時間領域振幅分布及びスペクトル振幅の所定の包絡線を含むグループのどれかを備える。
【００２７】
従って、本発明による該方法の追加実施形態では、該処理済みの表現された第１信号は、所定の時間領域振幅分布に従って配列される。
【００２８】
本発明による該方法の依然として追加の実施形態では、該処理済みの表現された第１信号は、スペクトル振幅の所定の包絡線に従って配列される。
【００２９】
本発明による該方法の依然として追加の実施形態では、該処理済みの表現された第１信号が、スペクトル振幅の所定の包絡線に従って配列される。
【００３０】
本発明に従って、該周波数領域と時間領域の両方で所定の品質基準を正確に満たす信号を提供するために、該時間領域振幅分布及び該スペクトル振幅の包絡線の内の少なくとも１つが、反復プロセスにより取り組まれる。振幅及び周波数波形整形は、両方の形状が妥当な正確度の範囲内で要件を満たすまで必要とされる回数に繰り返されてよい。
【００３１】
本発明のある実施形態では、該反復プロセスは、任意の反復ステップ後の、該処理済みの表現された第１信号の該時間領域振幅分布及びスペクトル振幅の包絡線の何れかと所定の時間領域振幅分布及びスペクトル振幅の所定の包絡線との比較を行う。
【００３２】
該時間領域特性が該要件に十分に近いかどうかを見抜くために、周波数波形整形後に常時ドメイン特性チェックを実行する必要はないことが認められた。該反復を停止するのか、あるいは続行するのかの決定を可能にするには、実際には、波高因子要件の簡略なチェックで十分であることが判明している。該信号の該波高因子は、該信号の該トーンの該平均つまりｒｍｓ値に比較される、該信号の該トーンの該最大、つまりピーク振幅の関係として定義される。
【００３３】
前述されたように本発明による該方法は、特に、とりわけ漏話雑音の生成に適している。
【００３４】
本発明による該方法の第２態様において、進入雑音の該特性を有する信号が生成されなければならない場合、該所定品質基準またはそれぞれの所定の品質基準は、少なくとも１つの変調済み搬送波を備え、該変調済み搬送波またはそれぞれの変調済み搬送波は、搬送周波数、搬送振幅、変調深度、及び変調幅から構成されるグループのどれかを含む。
【００３５】
前記に示された１つまたは複数の品質基準に従って該表現された第１信号を波形整形することにより、ある特定の時間領域振幅分布、及びスペクトル振幅の所定の包絡線を有する、ある特定の種別の進入雑音を表現する信号は、容易に且つ非常に効率的に提供できる。
【００３６】
本発明の該方法に従って、該所定品質基準またはそれぞれの所定の品質基準を満たす該信号は、前記に開示されたように処理される複数の信号を結合することにより提供できる。
【００３７】
例えば、通信網または通信装置の試験で該信号を使用する場合、該処理済みの表現された信号は、とりわけ、例えば、ＦＦＴアルゴリズムを使用して該周波数領域から該時間領域に変換されなければならない。
【００３８】
本発明は、さらに、前記に開示されたように、該第１態様と該第２態様に従って生成される該信号を結合することを提供する。しかしながら、他の信号成分も含まれてよい。
【００３９】
特に、本発明の該方法に従って、該所定品質基準またはそれぞれの所定の品質基準を有する該信号は雑音信号である。
【００４０】
本発明の第３態様では、通信システムの動作を試験する方法が開示され、該方法は、
−前記に開示された本発明の該方法に従って規定の品質基準を有する信号を生成するステップと、
−試験中の該通信システムを通して該信号を転送するステップと、
を備える。
【００４１】
該信号は、コードフォーマットの命令のセットを使用して生成、記憶し、装置上で規定の順序で実行できる。このような命令のセットが、コンピュータでコンパイルされ、該コンピュータまたはコンピュータのネットワーク、フロッピーまたはＣＤ−ＲＯＭに、あるいはインターネットを通じて記憶されるソフトウェアコードである場合がある。該ソフトウェア及び／または生成された信号は、任意波信号発生器（ＡＷＧ）カードに記憶することもでき、該ＡＷＧは、信号を生成するため、あるいはメモリから記憶されている信号を再生するために使用できる。従って、本発明による該方法の実行または使用で使用できるデータ搬送波に記憶されている使用可能な信号のライブラリを有することができる。
【００４２】
該通信システムは、ｘＤＳＬモデムなどの装置、あるいはこのようなモデム内またはモデム用のチップ、あるいはネットワーク内のケーブル、または（電気）通信用のネットワークである場合がある。
【００４３】
本発明の第４態様では、通信システム上または通信システム内で使用するための信号を配列する追加の方法が開示される。好ましくは、該信号は雑音信号である。該信号は、該周波数領域及び該時間領域内の所定の特性を備えた不規則な信号である漏話雑音を備えてよい。該信号は、さらに、ディスクリート周波数スペクトルを有するｒｆｉ−トーンを備えることがある。また、該信号に他の信号成分が含まれることがある。
【００４４】
該方法は、
−時間領域振幅分布を有する時間領域内の第１信号を表現し、該信号が該周波数領域にスペクトル密度を有し、それにより表現された信号を達成するステップと、
−非線形変換に従って該表現された信号を処理し、該非線形変換が少なくとも１つの所定の品質基準を達成するステップと、
−該表現された信号の該時間領域振幅分布が、規定の時間領域振幅分布の少なくとも逆関数で処理されるステップと、
を備える。
【００４５】
該方法は、さらに、該表現された信号の該時間領域振幅分布を、該規定の時間領域振幅領域と比較し、その後、該規定の時間領域振幅分布に近づく時間領域振幅分布を有する処理済みの表現された信号を達成するために、該非線形変換を配列するステップを備えてよい。
【００４６】
本発明の第５態様では、該表現された信号の該時間領域振幅分布を、該規定の時間領域振幅分布と比較し、その後、該規定の時間領域振幅分布に近づく時間領域振幅分布を有する処理済みの表現された信号を達成するために該非線形変換を配列するステップをさらに備える方法が開示される。
【００４７】
本発明の該第５態様に従って、該方法は、時間領域内の第１信号を振幅分布で表現し、該信号が該周波数領域内のスペクトル密度を有し、それにより表現された信号を達成するステップと、該周波数領域内の該信号表現の少なくとも部分を評価するステップ及びその後該周波数領域内で該表現された信号を処理するステップとを含む、該周波数領域内で該表現された信号をフィルタリングするステップとを備えることもある。
【００４８】
本発明の第４態様と第５態様の該方法を結合することができる。本発明の第４態様と第５態様の該方法は、所定の振幅分布を有する、及び／または所定のスペクトル密度を有する、あるいは１つの振幅分布を有する、及び／または所定の品質基準に従って１つのスペクトル密度を有する信号を、さまざまな反復ステップで作成できる。該所定の品質基準が、該信号の波高因子、すなわち該信号の該トーンの平均値つまりｒｍｓ−値に比較される、該信号のトーンの最大つまりピーク値の関係である場合がある。前記に列挙されたような該処理ステップは、高速フーリエ変換（ＦＦＴ）または高速フーリエ逆変換（ＩＦＦＴ）のステップを備えることがある。該処理ステップは、信号の畳み込みまたは逆重畳、あるいは乗算または付加などの演算を含むこともある。
【００４９】
第４態様の該方法では、該表現された信号の該振幅分布は、該規定の振幅分布の逆関数を含むことがある該規定の振幅分布の関数を含めて処理される。
【００５０】
本発明の第４態様と第５態様の列挙されるような該方法は、さらに、該周波数領域内で該第１信号を変換するステップと、該周波数ドメイン内の該第１信号をスペクトル包絡線で乗算し、それにより乗算された信号を達成するステップと、その後時間領域内で該乗算された信号を表現するステップとを備えることがある。
【００５１】
該方法では、該周波数領域内のその表現の中の該第１信号は乱数、好ましくは、複素数、振幅を特徴付ける該複素数の係数、位相を特徴付ける該複素数の該引数の集合として生成することができ、該複素数の原則的にはそれぞれの実数及び／または虚数部分はガウス分布に従って選ぶことができる。該複素数のそれぞれは、該所定のスペクトル密度の該振幅に実質的には等しいことがある。
【００５２】
本発明の第６態様では、少なくとも１つの不規則雑音信号を備える信号が開示され、該不規則信号は、規定の品質基準に従って該時間領域内に振幅分布を、所定の品質基準に従って該周波数領域内にスペクトル密度を有し、該不規則信号は乱数の配列から構成される。該信号は、さらにディスクリート周波数スペクトルを備えることがある。該雑音信号は、コードフォーマットの命令のセットを使用し、規定の順序で実行され、生成できる。このような命令のセットは、コンピュータでコンパイルされ、該コンピュータまたはコンピュータのネットワークまたはフロッピーまたはＣＤ−ＲＯＭに、あるいはインターネットを通して記憶されるソフトウェアコードである場合がある。該ソフトウェアは、任意波信号発生器（ＡＷＧ）カードにも記憶でき、該ＡＷＧは該信号を生成する、あるいは該メモリから記憶されている信号を再生するために使用できる。従って、本発明の第４態様と第５態様の該方法の実行または使用で使用できる、使用可能な信号のライブラリを有することが可能である。
【００５３】
本発明の第７態様では、少なくとも１つの不規則雑音信号を備える信号を生成する方法が開示され、該不規則信号は、規定の品質基準に従って該時間領域内に振幅分布を有し、規定の品質基準に従って該周波数領域内にスペクトル密度を有し、該不規則信号は乱数の配列から構成され、該方法は、コードフォーマットの命令のセットを使用し、規定の順序で実行され、数の無作為な集合を生じさせるステップを備える。該方法は、さらに、ディスクリート周波数スペクトルを生成するステップを備えることがあり、該ディスクリート周波数スペクトルは角度測定機能を使用し、原則的に雑音特性のある該ディスクリート周波数のそれぞれを変調する。該不規則雑音信号及びディスクリート周波数スペクトルは、コードフォーマットの命令のセットを使用し、規定の順序で実行され、結合できる。
【００５４】
本発明の第８態様では、コードフォーマットで、規定の順序で実行可能な命令のセットが開示され、該命令のセットは不規則雑音信号及びディスクリート周波数スペクトルを生成するために配列され、該不規則信号は、所定の品質基準に従って該時間領域に振幅分布を有し、所定の品質基準に従って該周波数領域にスペクトル密度を有する。このような命令のセットは、コンピュータでコンパイルされ、該コンピュータまたはコンピュータのネットワークまたはフロッピーまたはＣＤ−ＲＯＭに、あるいはインターネットを通して記憶されるソフトウェアコードである場合がある。該ソフトウェアは、任意波信号発生器（ＡＷＧ）カードにも記憶でき、該ＡＷＧは、該信号を生成する、または該メモリから記憶されている信号を再生するために使用できる。従って、使用できる、使用可能な信号のライブラリを有することが可能である。該ソフトウェアはＣ−コードであるか、またはＭＡＴＬＡＢ環境でコンパイルできる。
【００５５】
本発明の第９態様では、コードフォーマットを取り、規定の順序で実行可能であり、装置上でコンパイルされる命令のセットを備える通信システムの動作を試験するためのシステムが開示され、該命令のセットは、不規則雑音信号及びディスクリート周波数スペクトルの少なくとも１つを備える雑音信号を生成するために配列され、該不規則信号は、規定の品質基準に従って該時間領域に振幅分布を有し、規定の品質基準に従って該周波数領域にスペクトル密度を有する。本発明の本態様による該試験システムは、該雑音信号を生成するための障害ジェネレータを備えることがある。
【００５６】
該障害ジェネレータを試験される該通信システムに接続する該接続要素（変圧器、能動素子、減衰器等）が、不必要な周波数に依存する応答を有することがある。該不必要な周波数に依存する応答は、例えば、該障害ジェネレータ内で特定の試験信号を生成することによって測定できる。該不必要な周波数に依存する応答は、該接続要素の該不必要な周波数に依存する応答で除算される該信号の該所望されるスペクトル密度を乗算することにより補償できる。
【００５７】
本発明の第１０態様では、ｘＤＳＬモデムなどの通信システムの動作を試験する方法が開示される。該方法は、１つの該モデムによって交信される信号に重畳するステップを備え、信号は不規則雑音信号及びディスクリート周波数スペクトルの少なくとも１つを備え、該不規則信号は、規定の品質基準に従って該時間領域に振幅分布を有し、規定の品質基準に従って該周波数領域内にスペクトル密度を有し、該雑音信号は、さらに乱数の配列から構成される。
【００５８】
本発明の第１１態様では、通信システムの動作の質を試験する方法が開示される。該方法は、１つの該モデムによって交信される信号に重畳するステップを備え、信号は不規則雑音信号及びディスクリート周波数スペクトルの少なくとも１つを備え、該不規則信号は規定の品質基準に従って該時間領域に振幅分布を有し、規定の品質基準に従って該周波数領域にスペクトル密度を有し、該雑音信号は、さらに乱数の配列から構成され、所定の品質基準に従って該交信される信号を評価する。
【００５９】
また本発明の第１２態様では、通信システムの設計及び／または製造を改善する方法が開示され、該方法は、１つの該モデムにより交信される信号上に重畳し、該モデムにより交信される信号上に重畳し、信号は不規則雑音信号及びディスクリート周波数スペクトルの少なくとも１つを備え、該不規則信号が規定の品質基準に従って該時間領域内に振幅分布を有し、規定の品質基準に従って該周波数領域内にスペクトル密度を有し、該雑音信号がさらに乱数の配列から構成されるステップと、規定の品質基準に従って該交信される信号を評価するステップと、該交信される信号を評価するために該品質基準にさらに近く接近するために該モデムの該設計を反復して配列するステップとを備える。
【００６０】
本発明の第１３態様では、不規則雑音信号及びディスクリート周波数スペクトルの少なくとも１つを備える信号を含む電気通信網が開示され、該不規則信号は、規定の品質基準に従って該時間領域内に振幅分布を有し、該規定の品質基準に従って該周波数領域内にスペクトル密度を有し、該雑音信号は、さらに乱数の配列から構成される。
【００６１】
本発明の前述された態様及び実施形態の該特徴は結合できる。
前記に列挙された該信号、該方法及び該命令のセットは、電話線などの媒体または無線媒体上でさらに優れた品質の信号伝送を有することを可能にするだろう。信号のさらに優れた伝送は、通信システムのユーザ向けのさらに多くのサービスのさらに幅広い提供に対処する。
【発明を実施するための最良の形態】
【００６２】
本発明を教示する目的のために、本発明の該信号と方法とシステムの態様及び実施形態が後述される。本発明の他の代替実施形態及び同等な実施形態が、本発明の真の精神から逸脱することなく実践するために考えられ、変形されることは当業者により理解されるだろう。本発明の範囲は、添付請求項によってのみ制限されている。
【００６３】
本発明の実施形態では、ｘＤＳＬトランシーバなどの通信システムの動作を試験するためのシステムが開示される。稼動中のアクセスネットワークにおけるシステムの高い浸透率のシナリオのための試験装置のセットアップが説明される。
【００６４】
通信システム上または通信システム内で使用するための信号を配列する方法が開示される。
【００６５】
伝送性能試験の目的とは、稼動中のアクセスネットワークでのシステムの高い浸透率のシナリオにとって典型的な方法で、ｘＤＳＬトランシーバに応力をかけることである。この高浸透アプローチにより、以下が可能になる。
【００６６】
（ｉ）構成要素及びシステムの設計者が、性能を定量化し、設計を改良し、規格との準拠を証明するためにそれを使用できるようになる。
【００６７】
（ｉｉ）オペレータが、大部分の過剰状態に当てはまる配備規則を明示することができる。
【００６８】
図１は、考えられる試験セットアップ１の機能上の説明を示す。それは、以下を含む。
【００６９】
・本物のケーブルまたはケーブルシミュレータである試験ループ２、
・障害雑音を該試験ループ２に投入するための追加要素３、
・例えば、スペクトラムアナライザまたはｒｍｓ電圧計（図示されていない）などのレベル検出器５と連結される高インピーダンス、及び釣り合いの良い差動電圧プローブ４、及び
・試験中のｘＤＳＬトランシーバ（モデム）６、７
（ＰＯＴＳ信号及びＩＳＤＮ信号について）試験中のｘＤＳＬシステムに外部スプリッタが必要とされるとき、これらのスプリッタを試験中の該モデム６、７の中に含むことができる。
【００７０】
試験装置セットアップ１を通る該信号の流れはポートＴｘからポートＲｘへであり、上流性能及び下流性能を測定するには、トランシーバ位置と試験「ケーブル」端部の交換が必要とされることを意味する。ポートＲｘでの該受信信号レベルは、ポートＲｘだけではなくポートＴｘも試験中のｘＤＳＬトランシーバ（モデム）６、７と成端されるときに、ノードＡ２とＢ２間で測定される水準である。該障害ジェネレータ８は、この測定の間オフに切り替えられる。ポートＴｘでの該伝送信号レベルは、同じ条件下でノードＡ１とＢ１間で測定される水準である。
【００７１】
該障害ジェネレータ８が該試験セットアップ１の中に投入する必要のある該雑音は、周波数に依存している。該障害ジェネレータ８が該試験セットアップ１の中に投入する該雑音は、本物の（スペクトル汚染された）アクセスネットワークの現実的な表現とならねばならず、
（ａ）該試験ループ２の長さに依存し、
（ｂ）下流性能試験と上流性能試験について異なる。
【００７２】
ノードＡ２とＢ２間で測定されるこの障害雑音は、通常、不規則雑音、インパルス性の雑音、及び高調波雑音（ｒｆｉ−トーン）の混合物である。特性の集合が「ノイズプロファイル」として特定される。
【００７３】
信号レベルと雑音レベルは、釣り合いの良い差動電圧プローブ４を用いて厳密に調べられる。
【００７４】
完全に自動化された試験セットアップ１では、該試験ループ２、３、該電圧プローブ４とレベル検出器５、試験中の該モデム６、７、及び該障害ジェネレータ８は、破線により概略して示されるように、中央演算処理装置（ＣＰＵ）９に接続してよい。当業者は、該ＣＰＵ９との該接続が、該ＣＰＵ９による遠隔試験用のデータリンクを含んでよいことを理解するだろう。
【００７５】
該試験装置の使用に関連する定義は、以下のとおりである。
【００７６】
・完全信号バンド上で、このセットアップでｒｍｓ電圧Ｕｒｍｓ［Ｖ］を厳密に調べることは、Ｐ＝１０ｘｌｏｇ１０（Ｕｒｍｓ２／Ｒｖｘ１０００）［ｄＢｍ］に等しいＰ［ｄＢｍ］の電力レベルを意味する：
・（ヘルツ単位の）△ｆという小周波数バンド内で、このセットアップでｒｍｓ−電圧Ｕｒｍｓ［Ｖ］を厳密に調べることは、Ｐ＝１０ｘｌｏｇ１０（Ｕｒｍｓ２／Ｒｖｘ１０００／△ｆ）［ｄＢｍ／Ｈｚ］に等しいそのフィルタリングされたバンド内のＰ［ｄＢｍ／Ｈｚ］という電力スペクトル密度レベルを意味する。
【００７７】
・帯域幅△ｆは、該−３ｄＢ帯域幅ではなく、該フィルタの該雑音帯域幅を特定する。
【００７８】
図２は、本発明による通信システム上または通信システム内で使用するための信号Ｕｉ（ｔ）ｉ＝１，２，３．．．を配列するための該方法の実施形態を概略的に示す。該信号は、漏話雑音、つまり該周波数領域内及び該時間領域内に所定の特性を有する不規則信号を備えてよい。
【００７９】
図２のフロー［１］によって表されるように、「第１信号」と記載されたブロック１０において、本発明における方法は、それぞれがスペクトル振幅特性及び位相特性を有する複数の周波数成分を備える第１集合である第１信号を表現するステップと、「配列」と表示されているブロック１１において、少なくとも一つの所定の品質基準に従ってブロック１０において表現された第１信号にスペクトル振幅特性を配列すると共に、ランダム位相特性を付与するように配列処理を実行し、それにより処理済みの表現された信号を得る工程が示されている。
【００８０】
該第１信号は、各周波数成分のスペクトル振幅及び位相を指定する数の第１集合によって表現されてよい。さらに、該第１信号は、実数部と虚数部分を有する複素数の第２集合によって表現されてよく、部分は組み合わされて各周波数成分のスペクトル振幅及び位相を指定する。すなわち、複素数の係数はスペクトル振幅を特徴付けるが、該複素数の該引数が該周波数成分の該位相を特徴付ける。
【００８１】
本発明に従って、ブロック１０で表現される第１信号は、ブロック１１において、ランダム位相特性を配列するために処理される。ただし、該規定の品質基準に従って該信号の該周波数を波形整形するために、予めランダム位相特性を有する表現される第１信号から配列処理をすると、それは該周波数成分の該スペクトル振幅を波形整形するのに十分である。
【００８２】
本発明におけるブロック１０における方法は、該第１信号を時間領域内において表現するステップを含んでいても良く、当該ステップでは、該第１信号はそれぞれが該時間領域内の該第１信号の振幅を指定する当該第１信号の複数の数からなる第３集合によって表現されるものである。
そして、当該複数の数からなる第３集合を、例えば、ＦＥＴアルゴリズムを使用して、該時間領域から周波数領域に変換することにより、各周波数成分のスペクトル振幅及び位相を特定する複数の数からなる第４集合が得られる。
当該第４集合は、上記したと同じように、少なくとも１つの所定の品質基準に従ってスペクトル振幅特性を配列すると共に、ランダム位相特性を付与するように配列処理を実行し、そのランダム位相特性を配列するだけでなく、該少なくとも１つの所定の品質基準に従って処理されなければならない。
【００８３】
「変換」と表示されたブロック１２において、例えば逆ＦＦＴアルゴリズムを使用して、該処理済みの表現された信号を周波数領域から時間領域に変換することにより、該少なくとも１つの所定の品質基準を有する該信号Ｕ１（ｔ）が最終的に生成される。
【００８４】
該信号Ｕ１（ｔ）は、スペクトル包絡線及び／またはプリエンファシス特性などの該周波数領域内の少なくとも１つの所定の品質基準を満たすが、該時間領域内の所定の品質基準を有する信号を提供することがまだ必要とされる可能性がある。
【００８５】
図２のフロー［２］と［３］に示されているように、「振幅波形整形」を行うブロック１３と「周波数波形整形」を行うブロック１４により、該時間領域内の該品質基準は、所定の振幅分布及び／またはスペクトル振幅の所定の包絡線を備えてもよい。
【００８６】
更に、本発明の方法においては、様々な反復ステップを利用して信号を作成するステップを備えることも出来る。
即ち、図２［４］に示すように、「試験波形整形」と表記されているブロック１５及び帰還ループ１６により、該信号は、規定の時間領域振幅分布及び／又はスペクトル振幅の規定の包絡線及び／又は規定の品質基準に従うスペクトル密度を有することが可能である。
即ち、図２の１のフローは、本発明における雑音信号の基本的な形成方法を示し、図２の２のフローは、図２の１のフローで得られた雑音信号Ｕ１（ｔ）に対して、品質基準に対応するより望ましい雑音信号を得るために振幅波形整形を付加することにより、好ましい振幅分布をもつ雑音信号Ｕ２（ｔ）を形成するものである。
又、図２の３のフローは、図２の２のフローで得られた雑音信号Ｕ２（ｔ）に対して、更に、「周波数波形整形」を付加して、雑音信号Ｕ３（ｔ）を得るものである。
つまり、雑音信号Ｕ１（ｔ）に対して、振幅分布と周波数整形との組合せ処理を行い、品質基準により特定されるスペクトル包絡線をもつ雑音信号Ｕ３（ｔ）を得るものである。
一方、図２の４のフローは、図２の３のフローで得られた雑音信号Ｕ３（ｔ）を当該雑音信号Ｕ３（ｔ）が形成された工程を繰り返すことにより、好ましい振幅分布と好ましいスペクトル包絡線との組合せを有する最終的な雑音信号Ｕ４（ｔ）を得るものである。
【００８７】
該少なくとも１つの所定の品質基準が、該信号の該平均つまりｒｍｓ値に比較される該信号の該トーンの最大つまりピーク値の関係である該信号の波高因子である場合がある。
【００８８】
本発明の上記信号は、コードフォーマットを取り、既知の順序で実行可能であり、装置上でコンパイルされる命令のセットを使用して、生成、記憶できる。このような命令のセットは、コンピュータでコンパイルされ、該コンピュータまたはコンピュータのネットワークまたはフロッピーまたはＣＤ−ＲＯＭに、またはインターネットを通して記憶されるソフトウェアコードである場合がある。該ソフトウェアは、任意波信号発生器（ＡＷＧ）カードにも記憶でき、該ＡＷＧは該信号を生成する、または該メモリから記憶されている信号を再生するために使用できる。従って、本発明の第１態様と第２態様の該方法の実行または使用で使用できる、使用可能な信号のライブラリを有することが可能である。該通信システムが、ｘＤＳＬモデム６、７などの装置、またはこのようなモデム６、７内の、またはモデム６、７用のチップ、または電気通信用のネットワークである場合がある。
【００８９】
図２のブロック１３「振幅波形整形」での該処理は、該時間領域特性に対する変更または制御を達成するために行われる。該信号内の高振幅ピークまたはトーンを増幅する振幅歪み（変換）関数Ｑ（ｘ）が選ばれる。非線形変換関数Ｑ（ｘ）は、該信号の該実際の振幅分布関数及び該規定の振幅分布関数から再構築できる。
【００９０】
０〜Ｔの間の該時間期間ｔ内での雑音信号ｆ（ｔ）の場合、該信号の該振幅分布Ｆ（ａ）は、絶対値の中の該雑音ｆがａより大きい該時間の関数として定義される。つまり、該振幅分布Ｆ（ａ）は、該雑音ｆ（ｔ）の値の絶対値がａ以上である場合の確率を示している。Ｇ（ａ）が（機能拡張ガウスなどの、以下参照）該所定の振幅分布であり、Ｇ−１（ａ）がその逆関数である場合、該雑音信号ｆ（ｔ）から中間信号または最終信号ｇ（ｔ）を作成するための該非線形変換関数Ｑ（ｘ）は、以下のように定義できる。
【００９１】
【数１】
2008219864.tif 000003
【００９２】
【数２】
2008219864.tif 000004
ｘ＜＞０の場合、ｓｉｇｎ（ｘ）＝ｘ／｜ｘ｜。ｘ＝０の場合ｓｉｇｎ（ｘ）＝０
即ち、上記式（２）は、最終信号Ｇ（ｔ）（図２の２のフローの場合であればＧ（ｔ）＝Ｕ２（ｔ））が入力信号値ｆ（ｔ）（図２の２のフローの場合であればｆ（ｔ）＝Ｕ１（ｔ））からどのように計算されるかを示したものであり、上記式（２）は、この演算が非線形変換関数Ｑ（ｘ）を使用して実行されることを示している。
一方、上記式（１）は、当該非線形変換関数Ｑ（ｘ）が振幅分布Ｆと望ましい振幅分布Ｇとから如何に計算されるかを示している。
その結果、ｇ（ｔ）は、該規定の振幅分布Ｇ（ａ）を有するだろう。多くのケースでのＱ（ｘ）は関数である場合があるが、数値的に構築することもできる。機能拡張ガウス関数の例は以下のとおりである。
【００９３】
ガウス型雑音の該振幅分布は以下のとおりである。
【００９４】
【数３】
2008219864.tif 000005
【００９５】
【数４】
2008219864.tif 000006
σは、該信号の該ＲＭＳ値である。
【００９６】
該「機能拡張」ガウス分布は、以下のとおりに定義される。
【００９７】
【数５】
2008219864.tif 000007
ＶＲＭＳが所望された雑音サンプルのＲＭＳ値であり、Ｃｆが所望された波高因子である場合、以下を選択する。
なお、波高因子は、段落００８７で既に定義している公知のファクタである。
【００９８】
【数６】
2008219864.tif 000008
【００９９】
【数７】
2008219864.tif 000009
有効である事が判明したαの典型的な値は、０．００１と０．０１の間の規模であり、これが真のガウス分布から分散される機能拡張ガウスの偏差を表す。
【０１００】
図２のブロック１４、周波数波形整形では、該信号の該周波数領域特性が、該品質基準または各品質基準に対するより近い一致を達成するための提示された処理ステップとして改善される。該補正された周波数曲線は、例えば、該信号Ｕ１（ｔ）に振幅波形整形を施して得られた信号Ｕ２（ｔ）の該測定されたスペクトル密度を使って所定のスペクトル密度を比較する（除算する）ことによって達成できる。この例は、ＦＦＴ関数の畳み込みを用いる続編で説明される最善の態様で示される。
【０１０１】
図２の部分［４］では、ここに前記に詳説された該ステップの反復手順がどのようにして、通信システム内または通信システム上で使用するための該終了信号または最終信号の追加改善策につながるのかが示されている。該反復手順、つまりブロック１５試験波形整形、及び反結合ループ１６による周波数波形整形の試験は、該規定の品質基準（複数の場合がある）が達成されるまで実行される。
【０１０２】
図３は、概略図解で、進入雑音信号の提供のための本発明による該方法の使用を示す。
【０１０３】
前文（段落０００５）に開示されるように、進入雑音は、ディスクリート搬送波周波数ｆｃｉ、ｉ＝１，２，３．．．での複数の周波数成分により特徴付けられてよい。該搬送波周波数ｆｃｉでの該周波数成分は、それぞれ、搬送波振幅Ａｃｉ、ｉ＝１，２，３．．．を有し、当該周波数成分は、該それぞれの搬送波周波数ｆｃｉと関連付けられる該側周波数の該振幅である変調深度を有するだけではなく、該当する場合、変調幅、つまり該関連付けられる搬送波周波数ｆｃｉの左側と右側にある多数のディスクリート周波数も有する。
ここで、変調幅は、上記スペクトル成分が発生する周波数の範囲を意味し、変調深度は、それぞれの周波数成分の振幅を意味する。
【０１０４】
図３は、水平の、つまり周波数軸ｆ、及び垂直の、つまり振幅軸Ａを有する図解で、ほんの一例として、それぞれ搬送波振幅Ａｃ１とＡｃ２を有する２つの搬送波周波数ｆｃ１とｆｃ２から構成される信号を示す。
【０１０５】
そのそれぞれの側面にある該搬送波周波数ｆｃ１の回りには、それぞれが振幅Ａ１を有する３つの側周波数成分が配列され。そのそれぞれの側面上の搬送波周波数ｆｃ２での該周波数成分について、それぞれが振幅Ａ２を有する２つの側周波数成分が配列される。
【０１０６】
本発明に従って、少なくとも１つの所定の品質基準を有する信号を提供するために、該周波数成分の該振幅は、当該所定の品質基準が特定する好ましい特定の周波数に従って、図３の点線１と２によって開示されるように、波形整形されなければならない。
【０１０７】
本発明に従って、無作為位相特性を有する表現される第１信号から開始して、該振幅の該波形整形により、該無作為振幅位相は、該所定の品質基準を有して提供される該信号内に維持される。
【０１０８】
図４は、通信システム内または通信システム上で使用するため、特に、無作為位相特性がすでに規定されている場合に使用するために、信号を配列する方法の追加実施形態を概略して示す。該信号は、該周波数領域内に、及び該時間領域内に規定の特性を備える不規則信号である漏話雑音を備える。該信号は、さらに、既に段落０００５で説明されている公知のディスクリート周波数スペクトルを有するｒｆｉ−トーンを備えることがある。これは、図３に示す信号でもある。また、該信号には他の信号成分を含むことができる。
【０１０９】
該方法は、実質的には、図２で示す信号の配列方法と同じであるが、ブロック１６の雑音生成工程が追加されているものであり、そのほかのブロックには、図２の同一ブロックと同じ符号が付されている。そして、該方法は、さらに、時間領域内の、振幅分布を有する第１信号を表現し、該信号が該周波数領域にスペクトル密度を有し、それにより表現された信号を達成するステップと、非線形変換に従って該表現された信号を処理し、該非線形変換が規定の品質基準を達成するステップとを備えてよい。これは、図４、フロー［２−４］で振幅波形整形として図示される。
【０１１０】
該方法は、さらに、時間領域内の、振幅分布を有する第１信号を表現し、該信号が該周波数領域にスペクトル密度を有し、それにより表現された信号を達成するステップと、規定のスペクトル密度品質基準に従ったスペクトル密度を有する信号が達成されるまで、該表現された信号を処理するステップとを備えてよい。これは、図４、フロー［３−４］で周波数波形整形として示される。該周波数波形整形ステップは、該周波数領域内の該信号表現の少なくとも一部（例えば、フィルタリングが適用される周波数成分等）を評価するステップと、それ以降該周波数領域内の該信号表現を処理するステップとを含む、該周波数領域内の該表現された信号をフィルタリングするステップも含むことがある。
【０１１１】
該方法は、さまざまな反復ステップで信号を作成するステップを備えてもよい。図４、フロー［４］を参照すること。このようにして、該信号は、規定の時間領域振幅分布、及び／または所定のスペクトル密度、あるいは所定の品質基準に従った時間領域振幅分布及び／またはスペクトル密度を有することができる。該所定の品質基準は、該信号の波高因子、つまり該信号の該トーンの平均値に比較される該信号の該トーンの最大つまりピーク値の関係である場合がある。該信号は、コードフォーマットを取り、規定順序で実行可能であり、装置上でコンパイルされる命令のセットを使用して、生成、記憶できる。
【０１１２】
同様に、該命令のセットは、コンピュータでコンパイルされ、該コンピュータ、またはコンピュータのネットワークまたはフロッピーまたＣＤ−ＲＯＭに、あるいはインターネットを通して記憶されるソフトウェアコードである場合がある。該ソフトウェアコードは、任意波信号発生器（ＡＷＧ）カードに記憶することもでき、該ＡＷＧは該信号を生成する、あるいは該メモリから記憶されている信号を再生するために使用できる。従って、本発明の該方法の実行または使用で使用できる、使用可能な信号のライブラリを有することが可能である。該通信システムは、ｘＤＳＬモデムなどの装置、あるいはこのようなモデム内の、またはこのようなモデム用のチップ、あるいは電気通信用ネットワークである場合がある。詳細には、以下の実施形態が図４に示される。
【０１１３】
ソフトウェアを使用して、乱数が生成される、ブロック１６「雑音の生成」。ハードウェアでは、白色雑音を発生できる。該乱数は、規定のスペクトル密度が達成されるまでフィルタリングされる。生成される該乱数は、それぞれ周波数成分を表現する。規定のスペクトル密度を達成するための必要な処理は、該複素数の該振幅を基準化することにより実行され、それ以後、該所望の雑音信号を作成するために、ＩＦＦＴ処理が実行される。該方法を実行する別の方法は、各周波数成分の該位相を表現する乱数を生成することであり、それ以後、該複素数の該振幅は、規定のスペクトル密度に近づくまたは規定のスペクトル密度に等しくなるように配列される。
【０１１４】
ブロック１３、「振幅波形整形」における処理は、該時間領域特性に対する影響または制御を達成するために行われる。図２において説明した非線形変換関数Ｑ（ｘ）の選択の方法が、本具体例でも使用されるものであり、図５は、その一具体例であって、図示されるような緯軸に示される信号の振幅の内、該信号中の高い振幅ピークまたはトーンを増幅する振幅分布（変換）関数Ｑ（ｘ）が選ばれる。該非線形変換関数Ｑ（ｘ）は、該信号の該実際の振幅分布関数、及び等式（１から７）に関連して前記に開示されたような該規定の振幅分布関数から再構築できる。
【０１１５】
図４のブロック１４、「周波数波形整形」では、該信号の該周波数領域特性が改善される。つまり、当該信号のスペクトル包絡線が、好ましいスペクトル包絡線により近くまで接近することを意味する。該補正された周波数曲線は、該（中間）信号の該測定スペクトル密度を通る規定のスペクトル密度を比較する（除算する）ことによって達成される。この例は、ＦＦＴ関数の逆重畳を用いる続編で説明された最善の態様の実施形態に示される。
【０１１６】
図４のフロー［４］では、図２のフロー［４］と同様に、再び、前記に詳説されたステップの反復手順がどのようにして、最終的に通信システム内、又は通信システム上で使用するための良い改良につながるのかが示されている。該反復手順は、該規定の品質基準が達成されるまで実行される。
【０１１７】
前述されたように、本発明による該方法を用いると、漏話雑音及び進入雑音を表現する信号は、漏話雑音信号と進入雑音信号の両方から構成される信号を提供するために装置されてよい障害ジェネレータ８、図１を参照、のような装置を用いて生成できるが、他の信号成分は、必要とされる場合、提供される該出力信号に追加されてよい。
【０１１８】
提供される該信号は、本発明の実施形態において、例えば数の配列などの該時間領域内の数の第６集合として有利に提供できる。
【０１１９】
図６は、フロー型図で、パーソナルコンピュータ２０上で実行している本発明の例の実施形態を示す。該障害雑音は、ブロック２２「白色雑音信号」、ＦＦＴにより提供されるブロック２３「スペクトル波形整形」、ＩＦＦＴによりブロック２３の出力から生じるブロック２４「所望の雑音信号」を備えるＳＰＯＣＳと呼ばれるブロック２１によって生成され、その結果として生じる信号はＡＷＧカード２５上に記憶される。漏話シナリオ、つまりブロック２６では、さらにブロック２３により処理される雑音ＰＳＤが作成される、ブロック２７。
図６のソフトウエアを実行することによって、白色雑音信号２２、スペクトル波形整形信号２３及び所望の雑音信号２４がそれぞれ生成される。
【０１２０】
本発明の該命令のセットの最善の態様の実施形態の一例が参考として以下に開示される。ここに以下に示される該コードは、ＭＡＴＬＡＢ環境でコンパイルされる。該コードの機能性に関連するコメントは、％記号の後に示される。当業者にとって、提供される該コードは自明的である。
【０１２１】
図７から図９は、該最善の態様の実施形態で得られる結果を示す。図８は、該ノイズプロファイルの該ＰＳＤが加えられた、生成済みの雑音サンプルのスペクトルのプロットを示す。図９は、該時間領域内の該生成された雑音サンプルのプロットを示す。図１０は、該生成された雑音サンプルの該分布関数のプロットを示す。図１１は、該生成された雑音サンプルの該累積分布関数のプロットを示す。図７は、グラフィックユーザインタフェース（ＵＩ）及び該ＡＷＧ制御の設定値を示す。
％
関数ＤｅｍｏＩｍｐａｉｒ２；
％
％ＤｅｍｏＩｍｐａｉｒ２
％Ｍａｔｌａｂプログラミング言語でプログラミングされ、障害ジェネレータの基本アルゴリズムを証明するコード。
【０１２２】
％該証明されたアルゴリズムは、
％−連続スペクトルで雑音を生成するときの周波数領域特性と時間領域特性（スペクトル、確率分布）
％−離散スペクトルで雑音を生成するときの搬送波振幅、搬送波周波数、変調深度、及び変調幅
％該所定の品質基準を完全に支配する。
【０１２３】
％両方の種別の雑音は個別に計算され、数を含む配列として該時間領域内で表現される。
【０１２４】
％両方の種別の雑音は、これらの配列を要素の観点で追加することにより同時に使用可能になる。
【０１２５】
％（ｃ）２０００−２００１  ＫＰＮリサーチ
％
％デモ関数
％ＤｅｍｏＸｔａｌｋＮｏｉｓｅ  −連続雑音を作成するプロセスを示す
％ＤｅｍｏＩｎｇｒｅｓｓＮｏｉｓｅ  −離散雑音を作成するプロセスを示す
％
％主関数：雑音は、乱数を含む配列として表現される。
【０１２６】
％ＤｅｆｉｎｅＳｈａｐｅ  −全てのユーザによって定義可能なパラメータを初期化する
％ＣｒｅａｔｅＮｏｉｓｅＣｏｎｔ  −連続雑音を生成する
％ＣｒｅａｔｅＮｏｉｓｅＤｉｓｃｒ＿Ｆａｓｔ  −離散雑音を生成、高速アルゴリズム
％ＣｒｅａｔｅＮｏｉｓｅＤｉｓｃｒ＿Ｓｌｏｗ  −離散雑音を生成、低速アルゴリズム
％ＦｒｅｑｕｅｎｃｙＳｈａｐｅ  −連続雑音のスペクトル密度を修正する
％ＡｍｐｌｉｔｕｄｅＳｈａｐｅ  −連続雑音の振幅分布を修正する
【０１２７】
％
％サポート関数：
％ＣａｌｃＳｐｅｃ  −雑音の該スペクトル密度を計算する
％ＣａｌｃＮＢＳＶ  −雑音の該狭帯域信号電圧を計算する
％ＣａｌｃＣｒｅｓｔ  −雑音の該波高因子を計算する
％ＣａｌｃＤｉｓｔｒｉｂ  −雑音の該確率分布を計算する
％ＣａｌｃＣｕｍＤｉｓｔｒｉｂ  −雑音の該累加分布を計算する
％ＣａｌｃＳｍｏｏｔｈ  −本物のスペクトルアナライザにおいてのようにスペクトルを平滑化する
％ＣａｌｃＥｎｈａｎｃｅｄＧａｕｓｓＤｉｓｔｒｉｂｕｔｉｏｎ  −近ガウス分布のサンプル
％ＣａｌｃＤｅｍｏｄｕｌａｔｉｏｎ  −搬送波で変調された該雑音を計算する
％――――――――――――――――――――――――――――――――
Ｓｈａｐｅ＝ＤｅｆｉｎｅＳｈａｐｅ；
ＤｅｍｏＸｔａｌｋＮｏｉｓｅ（Ｓｈａｐｅ）；
ＤｅｍｏＩｎｇｒｅｓｓＮｏｉｓｅ（Ｓｈａｐｅ）；
％――――――――――――――――――――――――――――――――
関数［Ｕ，ｔ］＝ＤｅｍｏＸｔａｌｋＮｏｉｓｅ（Ｓｈａｐｅ）；
％――――――――――――――――――――――――――――――――
％連続スペクトルを用いた雑音の該生成を証明する
【０１２８】
％例えば、漏話試験用
Ｒ＝Ｓｈａｐｅ．Ｒ；
ＣＦ＿ｍｉｎ＝Ｓｈａｐｅ．Ｘｔａｌｋ．ＣＲ＿ｍｉｎ；
［Ｕ，ｔ］＝ＣｒｅａｔｅＮｏｉｓｅＣｏｎｔ（Ｓｈａｐｅ）；ｐｌｏｔ（ｔ，Ｕ）；    ｔｉｔｌｅ（Ｘ’ｔａｌｋｍｅｔｈｏｄ１’）；ｓｈｇ；ｐａｕｓｅ
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；ｐｌｏｔ（ｆ，Ｘ）；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏ１’）；ｓｈｇ；ｐａｕｓｅ
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；ｐｌｏｔ（ｆ，ｄＢｍ（Ｘ，Ｒ））；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ１’）；ｓｈｇ；ｐａｕｓｅ
［Ｐ，ｕ］＝ＣａｌｃＣｕｍＤｉｓｔｒｉｂ（Ｕ）；ｐｌｏｔ（ｕ，Ｐ）；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ１’）；ｓｈｇ；ｐａｕｓｅ
【０１２９】
％
Ｕ＝ＡｍｐｌｉｔｕｄｅＳｈａｐｅ（Ｕ，Ｓｈａｐｅ）；ｐｌｏｔ（ｔ，Ｕ）；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ２’）；ｓｈｇ；ｐａｕｓｅ
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；ｐｌｏｔ（ｆ，Ｘ）  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ２’）；ｓｈｇ；ｐａｕｓｅ
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；ｐｌｏｔ（ｆ，ｄＢｍ（Ｘ，Ｒ）））  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ２’）；ｓｈｇ；ｐａｕｓｅ
［Ｐ，ｕ］＝ＣａｌｃＣｕｍＤｉｓｔｒｂ（Ｕ）；ｐｌｏｔ（ｕ，Ｐ）；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ２’）；ｓｈｇ；ｐａｕｓｅ
【０１３０】
％
Ｕ＝ＦｒｅｑｕｅｎｃｙＳｈａｐｅ（Ｕ，Ｓｈａｐｅ）；ｐｌｏｔ（ｔ，Ｕ）；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ３’）；ｓｈｇ；ｐａｕｓｅ
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；ｐｌｏｔ（ｆ，Ｘ）  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ３’）；ｓｈｇ；ｐａｕｓｅ
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；ｐｌｏｔ（ｆ，ｄＢｍ（Ｘ，Ｒ）））  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ３’）；ｓｈｇ；ｐａｕｓｅ
［Ｐ，ｕ］＝ＣａｌｃＣｕｍＤｉｓｔｒｉｂ（Ｕ）；ｐｌｏｔ（ｕ，Ｐ）；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ３’）；ｓｈｇ；ｐａｕｓｅ
【０１３１】
％
ｆｏｒ  ｉ＝２：１０
ｉ
Ｕ＝ＡｍｐｌｉｔｕｄｅＳｈａｐｅ（Ｕ，Ｓｈａｐｅ）；
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；
％
Ｕ＝ＦｒｅｑｕｅｎｃｙＳｈａｐｅ（Ｕ，Ｓｈａｐｅ）；
［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；
％
ｉｆＣａｌｃＣｒｅｓｔ（Ｕ）＞ＣＦ＿ｍｉｎ，ｂｒｅａｋ；ｅｎｄ；
ｅｎｄ；
［Ｐ，ｕ］＝ＣａｌｃＣｕｍＤｉｓｔｒｉｂ（Ｕ）；
ｐｌｏｔ（ｔ，Ｕ）；  ｔｉｔｌｅ（’Ｘｔａｌｋｍｅｔｈｏｄ４’）；ｓｈｇ；％ｐａｕｓｅ
ｐｌｏｔ（ｆ，ｄＢｍ（Ｘ，Ｒ））；  （’Ｘｔａｌｋｍｅｔｈｏｄ４’）；ｓｈｇ；％ｐａｕｓｅ
ｐｌｏｔ（ｕ，Ｐ）；  （’Ｘｔａｌｋｍｅｔｈｏｄ４’）；ｓｈｇ；％ｐａｕｓｅ
【０１３２】
％――――――――――――――――――――――――――――――――
関数［Ｕ，ｔ］＝ＤｅｍｏｎＩｎｇｒｅｓｓＮｏｉｓｅ（Ｓｈａｐｅ）；
％――――――――――――――――――――――――――――――――
％例えば、進入試験のために、離散スペクトルを有する雑音の該生成を証明する
Ｒ＝Ｓｈａｐｅ．Ｒ；
［Ｕ，ｔ］＝ＣｒｅａｔｅＮｏｉｓｅＤｉｓｃｒ＿Ｆａｓｔ（Ｓｈａｐｅ）
％［Ｕ，ｔ］＝ＣｒｅａｔｅＮｏｉｓｅＤｉｓｃｒ＿Ｓｌｏｗ（Ｓｈａｐｅ）：％ｇｉｖｅｓｓａｍｅｒｅｓｕｌｔ
［Ｘ，ｆ］＝ＣａｌｃＮＢＳＶ（Ｕ，ｔ）；ｐｌｏｔ（ｆ，ｄＢｍ（Ｘ，Ｒ））；ｔｉｔｌｅ（’Ｉｎｇｒｅｓｓｍｅｔｈｏｄ’）；ｓｈｇ；ｐａｕｓｅｆｏｒＴｏｎｅＮｒ＝［１：２］
［Ｕｓａｃ，Ｕａｃ＿ｒｍｓ］＝ＣａｃｌＤｅｍｏｄｕｌａｔｉｏｎ（Ｕ，ｔ，Ｓｈａｐｅ，ＴｏｎｅＮｒ）：
ｐｌｏｔ（ｔ，Ｕａｃ）；ｔｉｔｌｅ（’ｄｅｍｏｄｕｌａｔｅｄ  ｉｎｇｒｅｓｓ  ｎｏｉｓｅ  ｏｆ  ｏｎｅ  ｃａｒｒｉｅｒ’）；ｓｈｇ；ｐａｕｓｅ
［Ｐ，ｕ］＝ＣａｌｃＤｉｓｔｒｉｂ（Ｕａｃ／Ｕａｃ＿ｒｍｓ）：
ｐｌｏｔ（ｕ，Ｐ）；ｔｉｔｌｅ（’ｄｉｓｔｒｉｂｕｔｉｏｎ  ｏｆ  ｄｅｍｏｄｎｏｉｓｅ’）；ｓｈｇ；ｐａｕｓｅ
ｅｎｄ；
【０１３３】
％――――――――――――――――――――――――――――――――
関数［Ｓｈａｐｅ］＝ＤｅｆｉｎｅＳｈａｐｅ；
％――――――――――――――――――――――――――――――――
％−スペクトル密度（この例では、本質的に方形波）
％−確率分布（この例では、ガウス近く）
％−トーン及び変調
％−平方根あたりのボルト単位のスペクトル（Ｈｚ）
に関して、生成されなければならない雑音のノイズプロファイルを作成する。
【０１３４】
【数８】
2008219864.tif 000010
【０１３５】
％
Ｓｈａｐｅ．Ｎ＝Ｎ；  ％時間サンプルの数
Ｓｈａｐｅ．ｍ＝ｍ；  ％周波数サンプルの数
Ｓｈａｐｅ．ｄＦ＝Ｆｍａｘ．／（ｍ−１）  ％周波数間隔
Ｓｈａｐｅ．ｄＴ＝１／（Ｎ＊Ｓｈａｐｅ．ｄＦ）；  ％時間間隔
Ｓｈａｐｅ．Ｒ＝Ｒ；  ％所望される雑音源のインピーダンス
％漏話雑音ターゲットを定める（スペクトル密度及び振幅分布）
【０１３６】
【数９】
2008219864.tif 000011
【０１３７】
％――――――――――――――――――――――――――――――――
関数［Ｕ，ｔ］＝ＣｒｅａｔｅＮｏｉｓｅＣｏｎｔ（Ｓｈａｐｅ）；
％――――――――――――――――――――――――――――――――
％所定の周波数領域特性（スペクトル）であるが、制御されていない時間領域特性（分布）で雑音電圧Ｕ（ｔ）を作成する
Ｎ＝Ｓｈａｐｅ．Ｎ；  ％生成されるサンプル数
％Ｕ＝ｒａｎｄ（Ｎ，１）；  ％均一分散白色雑音
Ｕ＝ｒａｄｎ（Ｎ，１）；  ％ガウス分散白色雑音
Ｕ＝ＦｒｅｑｕｅｎｃｙＳｈａｐｅ（Ｕ，Ｓｈａｐｅ）；  ％波形整形された雑音
ｔ＝［０：Ｎ−１］’＊Ｓｈａｐｅ．ｄＴ；  ％関連時間軸
【０１３８】
％――――――――――――――――――――――――――――――――
関数［Ｕ，ｔ］＝ＣｒｅａｔｅＮｏｉｓｅＤｉｓｃｒ＿Ｆａｓｔ（Ｓｈａｐｅ）；
％――――――――――――――――――――――――――――――――
％ＡＭ変調された搬送波（ＲＦＩトーン）で電圧Ｕ（ｔ）を作成する。それぞれが、個々の所定の周波数、振幅、変調幅及び変調深度を有する。各搬送波トーンの下方側波帯の無作為な位相は、任意ＱＡＭ変調を、さらに制限されたＡＭ変調（完全ミラーリング）に（ミラーリングなしに）変換するためにミラー化される。
【０１３９】
％このアルゴリズム内でＸは、該（近高調波）進入雑音信号の該フーリエ級数の該成分を指すが、それは（擬似ランダム）漏話雑音信号のケースではスペクトル密度を指す。
【０１４０】
％計算時間は、サンプルの数に伴ないほぼ線形に増加する。
全計算時間の約８０％は、フーリエ逆変換により引き起こされる。
【０１４１】
％――――――――――――――――――――――――――――――――
Ｎ＝Ｓｈａｐｅ．Ｎ；  ％サンプルの数
ｍ＝Ｓｈａｐｅ．ｍ；  ％この数の半分
Ｎｃ＝ｒｏｕｎｄ（Ｓｈａｐｅ．ＩＮｇｒｅｓｓ．ＴｏｎｅＦ／Ｓｈａｐｅ．ｄＦ）＋１；  ％搬送波周波数の指数（正のみ）
Ｎｍ＝ｒｏｕｎｄ（Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．Ｍｏｄｗｉｄｔｈ／Ｓｈａｐｅ．ｄＦ／２）；  ％変調成分の数
Ｘｃ＝０．５＊Ｓｈａｐｅ．ＩＮｇｒｅｓ．Ｔｏｎｅ．Ｕ；％搬送波振幅の振幅
Ｘｍ＝Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＭｏｄＤｅｐｔｈ．＊Ｘｃ．／ｓｑｒｔ（２＊Ｎｍ）；  ％変調バンドの振幅
【０１４２】
Ｘ＝ｚｅｒｏｓ（Ｎ，１）；  ％初期化
Ｘｃ＝Ｘｃ．＊ｅｘｐ（ｊ＊１０００＊ｒａｎｄ（ｓｉｚｅ（Ｘｃ）））；％無作為搬送波位相
Ｘｃｃ＝（Ｘｃ．＊Ｘｃ）．／ａｂｓ（Ｘｃ．＊Ｘｃ）；
ｆｏｒ  ｋ＝１：ｌｅｎｇｔｈ（Ｎｃ）％全変調済み搬送波について、以下を実行する：
Ｎｍｐ＝Ｎｃ（ｋ）＋［１：Ｎｍ（ｋ）］’；  ％上方側波帯周波数の位置を突き止める
Ｎｍｎ＝Ｎｃ（ｋ）−［１：Ｎｍ（ｋ）］’；  ％上方側波帯周波数の位置を突き止める
Ｘｍｐ＝Ｘｍ（ｋ）．＊ｅｘｐ（ｊ＊１０００＊ｒａｎｄ（ｓｉｚｅ（Ｎｍｐ）））；  ％上方側波帯を作成する
【０１４３】
Ｘｍｎ＝ｃｏｎｊ（Ｘｍｐ）＊Ｘｃｃ（ｋ）；  ％下方側波帯をミラー化する
Ｘ（Ｎｍｐ）＝Ｘｍｐ；  ％上方側波帯を挿入する
Ｘ（Ｎｍｎ）＝Ｘｍｎ；  ％下方側波帯を挿入する
ｅｎｄ；
Ｘ（Ｎｃ）＝Ｘｃ；  ％全ての搬送波を挿入する
Ｘ（Ｎ：−１；ｍ＋２）＝ｃｏｎｊ（Ｘ（２：ｃｅｅｉｌ（ｍ）））；％スペクトルを付加する
【０１４４】
％（負の周波数）
％Ｕ＝ｒｅａｌ（ｉｆｆｔ（Ｘ））＊Ｎ；  ％時間領域に変換する
Ｕ＝ｒｅａｌ（ｆｆｔ（Ｘ））；  ％時間領域に変換する
％（１０％高速）
ｔ＝Ｓｈａｐｅ．ｄＴ＊［０：Ｎ−１］；  ％関連時間軸
％――――――――――――――――――――――――――――――――
関数［Ｕ，ｔ］＝ＣｒｅａｔｅＮｏｉｓｅＤｉｓｃｒ＿Ｓｌｏｗ（Ｓｈａｐｅ）；
％――――――――――――――――――――――――――――――――
％所定の周波数でのＲＦＩトーン、振幅と変調帯域幅、及び変調深度で、電圧Ｕ（ｔ）を作成する。
【０１４５】
このアルゴリズムは単純で、非常に効率的であり、証明目的専用である。
それは、ＣｒｅａｔＮｏｉｓｅＤｉｓｃｒ＿Ｆａｓｔが同じ結果を戻すことを証明できる。
％――――――――――――――――――――――――――――――――
Ｎ＝Ｓｈａｐｅ．Ｎ；  ％生成されるサンプルの数
ｍ＝Ｓｈａｐｅ．ｍ；
ｆ＝［０：Ｎ−１］’＊Ｓｈａｐｅ．ｄＦ；
ｔ＝［０：Ｎ−１］’＊Ｓｈａｐｅ．ｄＴ；
Ｆｃ＝Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＴｏｎｅＦ；
Ｆｃ＝Ｓｈａｐｅ．ｄＦ＊ｒｏｕｎｄ（Ｆｃ／Ｓｈａｐｅ．ｄＦ）；  ％期間に整数を強制する
【０１４６】
％
Ｕ＝０；
ｆｏｒｋ＝１：ｌｅｎｇｔｈ（Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＴｏｎｅＦ）；
％］−ｃｒｅａｔｅ  ｎｏｉｓｙ  ｍｏｄｕｌａｔｅ，  ｈａｖｉｎｇＵ＿ａｖｇ＝０  ａｎｄＵ＿ｒｍｓ＝１．
Ｎｍ＝ｒｏｕｎｄ（Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＭｏｄＷｉｄｔｈ（ｋ）／Ｓｈａｐｅ．ｄＦ／２）；
Ｘｎ０＝（［１：Ｎ］＜＝（Ｎｍ＋１）  ％変調雑音振幅を波形整形する
Ｘｎ＝Ｘｎ０．＊ｅｘｐ（ｊ＊１０００＊ｒａｎｄ（Ｎ，１））；  ％変調雑音位相を波形整形する
Ｘｎ（１）＝０；  ％ＤＣ成分を排除する
Ｘｎ（Ｎ：−１ｍ＋２）＝ｃｏｎｊ（Ｘｎ（２：ｃｅｉｌ（ｍ）））：  ％スペクトルを付加する（負の周波数）
【０１４７】
Ｎｏｉｓｅ＝ｒｅａｌ（ｉｆｆｔ（Ｘｎ））；  ％時間領域に変換する
Ｎｏｉｓｅ＝Ｎｏｉｓｅ／ｓｑｒｔ（ｓｕｍ（Ｎｏｉｓｅ．＊Ｎｏｉｓｅ）／Ｎ）  ％ｒｍｓ＝１を強制する
％――変調を実行する
Ｃａｒｒｉｅｒ＝Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＴｏｎｅＵ（ｋ）  ＊ｃｏｓ（２＊ｐｉ＋Ｆｃ（ｋ）＊ｔ＋１０００＊ｒａｎｄ）；
Ｍｏｄｕｌａｔｅ＝Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＭｏｄＤｅｐｔｈ（ｋ）＊Ｎｏｉｓｅ；
Ｕ＝Ｕ＋Ｃａｒｒｉｅｒ．＊（１＋Ｍｏｄｕｌａｔｅ）；
ｅｎｄ；
【０１４８】
％――――――――――――――――――――――――――――――――
関数［Ｕ］＝ＦｒｅｑｕｅｎｃｙＳｈａｐｅ（Ｕ，Ｓｈａｐｅ）
％――――――――――――――――――――――――――――――――
％ターゲット波形整形により指定されるように、サンプルＵのスペクトルを整形し直す
％入力
％Ｕ：  サンプルの連続値
％ｆｓ：  同じ周波数
％スペクトル：  所望のＰＳＤ（Ｖ／ｓｑｒｔ（Ｈｚ）単位）
【０１４９】
％――――――――――――――――――――――――――――――――
Ｎ＝ｌｅｎｇｔｈ（Ｕ）；
ｍ＝ｌｅｎｇｔｈ（Ｓｈａｐｅ．Ｘｔａｌｋ．ｓｐｅｃ）；  ％ｍ＝Ｎ／２
ｔ＝［０：Ｎ−１］’＊Ｓｈａｐｅ．ｄＴ；
％周波数基準化を実行する
Ｓｃａｌｉｎｇ＝Ｓｈａｐｅ．Ｘｔａｌｋ．ｓｐｅｃ．／Ｃａｌｃｐｅｃ（Ｕ，ｔ）；
Ｘ＝ｆｆｔ（Ｕ）；  ％周波数領域に変換する
Ｘ（１）＝０；  ％ＤＣ成分を排除する
Ｘ（２：ｍ＋１）＝Ｘ（２：ｍ＋１）．＊Ｓｃａｌｉｎｇ；  ％スペクトルを基準化する（正の周波数）
Ｘ（Ｎ：−１：ｍ＋２）＝ｃｏｎｊ（Ｘ（２：ｃｅｉｌ（ｍ）））；  ％スペクトルを付加する（負の周波数）
Ｕ＝ｒｅａｌ（ｉｆｆｔ（Ｘ））；  ％時間領域に変換する
％――――――――――――――――――――――――――――――――
関数［Ｕ］＝ＡｍｐｌｉｔｕｄｅＳｈａｐｅ（Ｕ，Ｓｈａｐｅ）
％――――――――――――――――――――――――――――――――  ％
この関数は、関数Ｕの該振幅分布を振幅に依存する（非線形）歪み関数Ｑ（ｘ）で波形整形する。結果はＵ（ｔ）＝Ｑ｛Ｕ（ｔ）｝である。
【０１５０】
％
％ＦＦを該サンプルの実際の累積分布関数とし、ＧＧを所望の累積分布関数とすると、該歪み関数は、以下により示される。
【０１５１】
【数１０】
2008219864.tif 000012
【０１５２】
Ｕ０＝ｓｑｒｔ（ｓｕｍ（Ｕ．＊Ｕ）／ｌｅｎｇｔｈ（Ｕ））；  ％基準化ｆａｒｃｔ（正規化用）
％該歪み関数Ｑを計算する
［ＤｉｓｔＰ１，ＤｉｓｔＵ１］＝ＣａｌｃＣｕｍＤｉｓｔｒｉｂ（Ｕ／Ｕ０）；  ％実際の分布
Ｑ＝ｉｎｔｅｒｐ１（Ｓｈａｐｅ．Ｘｔａｌｋ．ＤｉｓｔＰ，Ｓｈａｐｅ．Ｘｔａｌｋ．ＤｉｓｔＵ，ＤｉｓｔＰ１）；  ％歪み関数
Ｕ＝Ｕ０＊ｉｎｔｅｒｐ１（ＤｉｓｔＵ１，Ｑ，ａｂｓ（Ｕ／Ｕ０））．＊ｓｉｇｎ（Ｕ）  ％歪みを実行する
【０１５３】
ｐｌｏｔ（ＤｉｓｔＵ１，Ｑ）；ｓｈｇ；％ｐａｕｓｅ
％――――――――――――――――――――――――――――――――
％関数［Ｘ，ｆ］＝ＣａｌｃＳｐｅｃ（Ｕ，ｔ）；
％――――――――――――――――――――――――――――――――
％信号の該スペクトル密度を、それが指定解像度帯域幅で「測定される」ときに計算する。
【０１５４】
ＲＢＷ＝１５３  ％ＲＢＷ：Ｕのスペクトルの所望解像度
Ｎ＝ｌｅｎｇｔｈ（Ｕ）；ｍ＝Ｎ／２；
ｄＴ＝ｔ（２）−ｔ（１）；  ％時間間隔
ｄＦ＝１／ｄｔ／Ｎ；  ％周波数間隔
ｆ＝［０：ｍ−１］’＊ｄＦ  ％全ての正の周波数
Ｘ＝ｆｆｔ（Ｕ）＊ｄＴ；  ％周波数領域へ
Ｘ＝ａｂｓ（Ｘ（２：ｍ＋１））；  ％ＤＣなし、負の周波数なし
Ｘ＝ｓｑｒｔ（ＣａｌｃＳｍｏｏｔｈ（Ｘ．＊Ｘ，ｆ，ＲＢＷ））；  ％帯域幅ＲＢＷでそれを平均化する
【０１５５】
％――――――――――――――――――――――――――――――――
関数［Ｘ，ｆ］＝ＣａｌｃＮＢＳＶ（Ｕ，ｔ）；
％――――――――――――――――――――――――――――――――
％信号の狭帯域信号電圧を、それが指定の解像度帯域幅で「測定される」ときに計算する。
【０１５６】
ＲＢＷ＝１Ｅ３  ％ＲＢＷ：Ｕのスペクトルの所望解像度
Ｎ＝ｌｅｎｇｔｈ（Ｕ）；ｍ＝Ｎ／２；
ｄＴ＝ｔ（２）−ｔ（１）  ％時間間隔
ｄＦ＝１／ｄＴ／Ｎ；  ％周波数間隔
ｆ＝［０：ｍ−１］’＊ｄｆ  ％全ての正の周波数
％Ｘ＝ｆｆｔ（Ｕ）＊ｄＴ＊ｓｑｒｔ（ｄＦ）；  ％周波数領域へ
％Ｘｆｆｔ（Ｕ）／Ｎ；  ％周波数領域へ
Ｘ＝ｆｆｔ（Ｕ）／Ｎ＊２；  ％周波数領域へ
Ｘ＝ａｂｓ（Ｘ（２：ｍ＋１））；  ％ＤＣなし、負の周波数なし
％Ｘ＝ｓｑｒｔ（ＣａｌｃＳｍｏｏｔｈ（Ｘ．＊Ｘ，ｆ，ＲＢＷ））；  帯域幅ＲＢＷでそれを平均化する
【０１５７】
％――――――――――――――――――――――――――――――――
関数［ＣＦ］＝ＣａｌｃＣｒｅｓｔ（Ｕ）
％――――――――――――――――――――――――――――――――
％ｒｍｓ−値で除算されるピーク値である、信号（Ｕ（ｔ）の波高因子を計算する。
【０１５８】
【数１１】
2008219864.tif 000013
【０１５９】
％――――――――――――――――――――――――――――――――
関数［Ｆ］＝ＣａｌｃＥｎｈａｎｃｅｄＧａｕｓｓＤｉｓｔｒｉｂｕｔｉｏｎ（ｘ，Ｃｆ）；
％―――――――――――――――――――――――――――――――――
％機能拡張ガウス分布として定義される蓄積ガウス分布関数Ｆ（ｘ）を生成する。
％Ｃｆ＝波高因子
【０１６０】
【数１２】
2008219864.tif 000014
【０１６１】
％――――――――――――――――――――――――――――――――
関数［ＤｉｓｔＰ，ＤｉｓｔＵ，Ｐ］＝ＣａｌｃＤｉｓｔｒｉｂ（Ｕ）
％――――――――――――――――――――――――――――――――
％信号Ｕの振幅分布を計算する。
【０１６２】
Ｎ＝ｌｅｎｇｔｈ（Ｕ）；
Ｎｂｉｎｓ＝１００；
［ｃｕｍｂｉｎ，ｘｘ］＝ｈｉｓｔ（Ｕ，Ｎｂｉｎｓ）；
ｄＸ＝ｘｘ（３）−ｘｘ（２）；
ＤｉｓｔＰ＝ｃｕｍｂｉｎ（：）／ＮｄＸ  ％ｓｕｍ（ＤｉｓｔＰ）＊ｄＸを強制する
【０１６３】
【数１３】
2008219864.tif 000015
【０１６４】
％――――――――――――――――――――――――――――――――
関数［ＤｉｓｔＰ，ＤｉｓｔＵ］＝ＣａｌｃＣｕｍＤｉｓｔｒｉｂ（Ｕ）
％――――――――――――――――――――――――――――――――
％信号Ｕの（逆方向）累積振幅分布を計算する
ｌｅｎ＝ｌｅｎｇｔｈ（Ｕ）；
％――分布関数を評価する
Ｕ＝ａｂｓ（Ｕ／Ｕｅｆｆ）；
【０１６５】
％関数を評価する
Ｎｂｉｎｓ＝ｍｉｎ（［５０，ｆｌｏｏｒ（ｌｅｎ／１０）］）；
［ｃｕｍｂｉｎ，ｘｘ］＝ｈｉｓｔ（Ｕ，Ｎｂｉｎｓ）：
ＢｉｎＷｉｄｔｈ＝ｘｘ（２）−ｘｘ（１）；
ＤｉｓｔＵ＝ｘｘ−ＢｉｎＷｉｄｔｈ／２；％ｓｈｉｆｔ
ｆｏｒ  ｎ＝［Ｎｂｉｎｓ−１：−１：１］；  ｃｕｍｂｉｎ（ｎ）＝ｃｕｍｂｉｎ（ｎ）＋ｃｕｍｂｉｎ（ｎ＋１）；ｅｎｄ
ＤｉｓｔＰ＝ｃｕｍｂｉｎ／ｌｅｎｓ；
％――他のルーチンのために、それらがこの結果を使用するときに、数値安定性を改善する。
【０１６６】
ＤｉｓｔＵ＝［０，ＤｉｓｔＵ（２：ｅｎｄ）］；  ％ｘ＝０で開始する
ＤｉｓｔＰ＝［ＤｉｓｔＰ，１／ｌｅｎ］；
ＤｉｓｔＵ＝［ＤｉｓｔＵ，ｘｘ（Ｎｂｉｎｓ）＋０．９９９＊ＢｉｎＷｉｄｔｈ／２］：  ％最終（単一）ポイントを追加する
ＤｉｓｔＰ＝［ＤｉｓｔＰ，ｌｅ−１００］；
ＤｉｓｔＵ＝［ＤｉｓｔＵ，ｘｘ（Ｎｂｉｎｓ）＋（１．００１）＊ＢｉｎＷｉｄｔｈ／２］；  ％安定性のための係数１．００１
【０１６７】
％――――――――――――――――――――――――――――――――
関数［ＰＳＤ，ｆｒｅｑ］＝ＣａｌｃＳｍｏｏｔｈ（ＰＳＤ，ｆｒｅｑ．ＲＢＷ）
％――――――――――――――――――――――――――――――――
％有限解像度帯域幅、及びガウス波形整形済み帯域フィルタで、本物のスペクトルアナライザを模倣する。
【０１６８】
％ＰＳＤ＝「スペクトル密度」のｄｅ  ｓｑｕａｒｅである「電力スペクトル密度」。平方ヘルツあたりボルト単位。
【０１６９】
％
Ｎ＝ｌｅｎｇｔｈ（ＰＳＤ）；
ｄｆ＝ｆｒｅｑ（２）−ｆｒｅｑ（１）；
ｂｒ＝３＊ｆｌｏｏｒ（ＲＢＷ／ｄｆ）；
ｆａｃｔｏｒ＝２＊ｂｒ＋１；
ｉｆ（ｆａｃｔｏｒ＞１）
．．ｆｆ＝ｄｆ＊（−ｂｒ：ｂｒ）；  ％平滑化間隔
【０１７０】
【数１４】
2008219864.tif 000016
【０１７１】
．．ｍａｓｋ＝ｍａｓｋ／ｓｕｍ（ｍａｓｋ）；  ％解像度帯域フィルタのガウスマスク
．．ｘｈｅｌｐ＝［ＰＳＤ；ｚｅｒｏｓ（２＊ｂｒ，１）］；
１ｐ＝ｆｉｌｔｅｒ（ｍａｓｋ，１，ｘｈｕｌｐ）；  ％スマートな畳み込み
ＰＳＤ＝ｙｈｅｌｐ（ｂｒ＋１：ｅｎｄ−ｂｒ）；
ｅｎｄ
％――――――――――――――――――――――――――――――――
関数［Ｕａｃ，Ｕａｃ＿ｒｍｓ］＝ＣａｌｃＤｅｍｏｄｕｌａｔｉｏｎ（ｕ，ｔ，Ｓｈａｐｅ，ＴｏｎｅｒＮｒ）；
％――――――――――――――――――――――――――――――――
％離散雑音の搬送波で変調された雑音を復調する。それは、離散雑音が該所定のパラメータを満たすことを証明するためだけの証明目的用である。
【０１７２】
％
％復調器は、同相にロックされない同期検出を使用する
％結果は、完全復調バンド上での未知の減衰である
％これは、ＤＣレベルを測定することによって後に補正される。そして、このＤＣレベルが１ボルトに正規化されるまで復調された信号を増幅する。
【０１７３】
％
％証拠（ｐｓｉは未知である）
ｌｅｔ  Ｕｒｆ＝ｃｏｓ（ｗ＊ｔ＋ｐｓｉ）＊（１＋Ｕａｃ）；  ％＝「１＋Ｕａｃ」で変調された搬送波
Ｕｃ＝ｃｏｓ（ｗ＊ｔ）；  ％＝搬送波
Ｕｄ＝Ｕｒｆ＊Ｕｃ；  ％＝同期検出信号
％ｔｈｅｎＵｄ＝１／２＊（ｃｏｓ（ｐｓｉ＋２＊ｗ＊ｔ）＋ｃｏｓ（ｐｓｉ））＊（１＋Ｕａｃ）；
Ｕｌｆ＝ｃｏｓ（ｐｓｉ）／２；（１＋Ｕａｃ）；  ％低域濾波後
Ｕｄｃ＝ｃｏｓ（ｐｓｉ）／２；  ％＝Ｕｌｆを平均化することによる
％Ｕａｃ＝（Ｕｌｆ／Ｕｄｃ）−１；
【０１７４】
％――――――――――――――――――――――――――――――――
Ｎ＝Ｓｈａｐｅ．Ｎ；
Ｆｃ＝Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＴｏｎｅＦ（ＴｏｎｅＮｒ）；  ％搬送波周波数を選択する
Ｆｃ＝Ｓｈａｐｅ．ｄＦ＊ｒｏｕｎｄ（Ｆｃ／Ｓｈａｐｅ．ｄＦ）；  ％期間の整数を強制する
ＭｏｄＷｉｄｔｈ＝Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＭｏｄＷｉｄｔｈ（ＴｏｎｅＮｒ）；
ＭｏｄＤｅｐｔｈ＝Ｓｈａｐｅ．Ｉｎｇｒｅｓｓ．ＭｏｄＤｅｐｔｈ（ＴｏｎｅＮｒ）；
【０１７５】
％
Ｕｄ＝Ｕ．＊ｃｏｓ（２＊ｐｉ＊Ｆｃ＊ｔ）；  ％変調済み搬送波の同期検出
Ｎｍ＝ｒｏｕｎｄ（１．１＊ＭｏｄＷｉｄｔｈ／Ｓｈａｐｅ．ｄＦ／２）；  ％フィルタ周波数を計算する
ｍａｓｋ＝ｚｅｒｏｓ（Ｎ，１）；ｍａｓｋ（［１：Ｎｍ，Ｎ−Ｎｍ：Ｎ］）＝１  ％フィルタを作成する
Ｕ１ｆ＝ｒｅａｌ（ｉｆｆｔ（ｆｆｔ（Ｕｄ）．＊ｍａｓｋ））；  ％低域濾波を実行する
Ｕｄｃ＝ｓｕｍ（Ｕｌｆ）／Ｎ；  ％正規化していないＤＣレベルを見つける
Ｕａｃ＝Ｕｌｆ／Ｕｄｃ−１；  ％全レベルを正規化し、ＤＣを削除する
【０１７６】
％
Ｕａｃ＿ｒｍｓ＝ｓｑｒｔ（ｓｕｍ（Ｕａｃ．＊Ｕａｃ／Ｎ））；  ％Ｕｄｃ＝１であるため、ＭｏｄＤｅｐｔｈに等しくなくてはならない
Ｓｃａｌｅ＝Ｕａｃ＿ｒｍｓ／ＭｏｄＤｅｐｔｈ；  ％「１」でなければならない
【図面の簡単な説明】
【０１７７】
【図１】図１はブロック図で、通信システムにおける、本発明の該方法に従って動作する障害ジェネレータを使用する性能試験のためのセットアップを示す。
【図２】図２は本発明による該方法の実施形態の流れ図種類を示す。
【図３】図３は進入雑音を生成するための本発明による該方法のある実施形態を図解で示す。
【図４】図４は本発明による方法の追加実施形態の流れ図を示す。
【図５】図５は信号内の高振幅ピークまたはトーンを増幅する振幅歪み（非線形変換）関数Ｑ（ｘ）を示す。
【図６】図６は保温発明の例の実施形態を流れ図で示す。
【図７】図７は本発明の実施形態に従って得られる結果を示す。
【図８】図８は本発明の実施形態に従って得られる結果を示す。
【図９】図９は本発明の実施形態に従って得られる結果を示す。
【図１０】図１０は本発明の実施形態に従って得られる結果を示す。
【図１１】図１１は本発明の実施形態に従って得られる結果を示す。
【符号の説明】
【０１７８】
１  試験セットアップ
２  試験ループ
３  追加要素
４  差動電圧プローブ
５  レベル検出器
６、７  試験中のｘＤＳＬトランシーバ（モデム）
８  障害ジェネレータ
９  中央演算処理装置（ＣＰＵ）
【図１】
2008219864.tif 000017
【図２】
2008219864.tif 000018
【図３】
2008219864.tif 000019
【図４】
2008219864.tif 000020
【図５】
2008219864.tif 000021
【図６】
2008219864.tif 000022
【図７】
2008219864.tif 000023
【図８】
2008219864.tif 000024
【図９】
2008219864.tif 000025
【図１０】
2008219864.tif 000026
【図１１】
2008219864.tif 000027