(57)【要約】要約書なし。【特許請求の範囲】
【請求項１】
１／ｆノイズの乱数の少なくとも１つの数列を生成するための方法において、
望ましいスペクトル値βを決定する工程と、
１／ｆノイズの生成される乱数の数を決定する工程と、
強度定数ｃｏｎｓｔを決定する工程と、
経過変数ｎのための初期値を決定する工程であって、
長さｎのベクトルｙの、望ましい数の要素ｙ（ｎ）が、１／ｆ分散される乱数から演算されるまでの間ずっとループ状に繰り返される以下の（１）〜（８）のステップを含む工程とを有する方法。
（１）経過変数ｎの実際の数値を１だけ増加させるステップ、
（２）シミュレーション時間工程［ｔｎ−１；ｔｎ］を決定するステップ、
（３）以下の規則に基づいて、次元（ｎ×ｎ）の分散行列Ｃの要素Ｃｉｊを決定するステップ、
Ｃｉｊ：＝ｃｏｎｓｔ・（−｜ｔｊ−ｔｉ｜β＋１＋｜ｔｊ−１−ｔｉ｜β＋１
＋｜ｔｊ−ｔｉ−１｜β＋１−｜ｔｊ−１−ｔｉ−１｜β＋１）
ｉ，ｊ＝１，．．．，ｎ
（４）分散行列Ｃを反転することにより、行列Ｃ−１を決定するステップ、
（５）以下の規則に基づいて、変数σを決定する工程（ｓｑｒｔは、関数「平方根」を表し、ｅ（ｎ，ｎ）は、（ｎ，ｎ）によって表示される、反転された分散行列Ｃ−１の要素を表す）、
σ＝ｓｑｒｔ（１／ｅ（ｎ，ｎ））
（６）長さｎのベクトルｘのｎ番目の構成要素を形成している、（０、１）正規分布乱数を決定するステップ、
（７）反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素と、先行している（ｎ−１）個のシミュレーション時間工程のために演算されている、ベクトルｙの（ｎ−１）個の要素とから、特に、以下の規則に基づいて変数μを形成するステップ（ｙ（ｎ−１）が、ベクトルｙの最初の（ｎ−１）個の構成要素を表しており、Ｃ−１・，ｎが、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素を表しており、Ｃ−１ｎ，ｎが、（ｎ，ｎ）によって表示される、反転された分散行列Ｃ−１の要素を表す）、
【数１】
2004524607.tif 000029
（８）以下の規則に基づいて、１／ｆ分散された乱数から長さｎのベクトルｙの要素ｙ（ｎ）を決定するステップ
Ｙ（ｎ）＝ｘ（ｎ）＊σ＋μ
【請求項２】
１／ｆノイズの乱数のｑの数列を同時に演算する方法であって、請求項１に記載のループ状に繰り返される（１）〜（８）のステップのうち（６）〜（８）のステップの代わりに、
（９）長さｎ、ｋ＝１，・・・，ｑのベクトルｘｋのそれぞれ最後の構成要素を形成する、（０，１）正規分散された量ｑの乱数ｘｋ，ｎを決定するステップ、
（１０）以下の規則に基づいて、ｑの変数μｋを形成するステップ（ｙ（ｎ−１），ｋが、先行しているシミュレーション時間工程のために計算されている、ベクトルｙｋの最初の（ｎ−１）個の構成要素を表しており、Ｃ−１・，ｎが、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素を表しており、Ｃ−１ｎ，ｎが、反転された分散行列Ｃ−１の、（ｎ，ｎ）により表示される要素を表している。これは、ｋ＝１，．．．，ｑに関して実施される）、
【数２】
2004524607.tif 000030
（１１）１／ｆ分散された乱数から、特に、以下の規則に基づいて、長さｎのベクトルｙｋの各ｎ番目の構成要素を形成しているｑ個の要素ｙｋ，ｎを演算するステップを含むことを特徴とする請求項１に記載の方法。
ｙｋ，ｎ＝ｘｋ，ｎ＊σ＋μｋ（ただしｋ＝１，．．．，ｑ）
【請求項３】
１／ｆノイズの対象となる技術的なシステムをシミュレーションするための方法であって、上記システムの入力部チャンネルに存在する変数を形成および／または決定する際に、請求項１または２に記載の方法に基づいて生成された乱数を使用する方法。
【請求項４】
請求項１〜３のいずれかに記載の方法を実施できるように設計されている、１／ｆノイズの乱数の数列を決定するためのコンピュータプログラム。
【請求項５】
請求項４に記載のコンピュータプログラムを有するデータ記録媒体。
【請求項６】
インターネット等のような電子データ網から、データ網に接続されているコンピュータにダウンロードされる請求項４に記載のコンピュータプログラム。
【請求項７】
請求項１〜３のいずれかに記載の１／ｆノイズの乱数の数列を決定するための方法を実施できるコンピュータシステム。
【発明の詳細な説明】
本発明は、１／ｆノイズの乱数（Ｚｕｆａｌｌｓｚａｈｌｅｎ）の数列（列、連続、Ｆｏｌｇｅｎ）を生成する（発生させる）ための方法に関するものである。
【０００１】
１／ｆノイズの乱数は、例えば、ノイズ影響を考慮した短期回路シミュレーションの場合に使用できる。１／ｆノイズを特定の周波数スペクトルを有する確率過程と解釈する。この周波数スペクトルは、下記数式
【０００２】
【数３】
2004524607.tif 000002
【０００３】
によって表せる。
【０００４】
１／ｆノイズ源は、多数の技術的および物理的システムにおけるノイズ影響のモデル化（Ｍｏｄｅｌｌｉｅｒｕｎｇ、形成化）、および、金融市場における事件についての判断と予測とのためのシステムに適している。特に、多数の電子部品、例えば、ｐｎダイオードおよびＭＯＳ電解効果トランジスタなどは、１／ｆノイズ源を備えている。
【０００５】
１／ｆノイズ源を、周波数スペクトルとして各１つのローレンツスペクトルを有する多数のノイズ源を合計することにより概算できる。このようなノイズ源は、例えば、線形の時間的に変化するシステム（ＬＴＩシステムとも称される）のシステム応答によって、モデル化される。そのシステム入力部には、白色雑音が印加される。この方法では、数により計算される微分方程式システムが、きわめて増大されることが短所である。このことにより、コンピュータシステムの演算時間が長くなり、メモリー需要が高くなる。このコンピュータは、１／ｆノイズの影響を受けるシステムをシミュレーションするために使用される。
【０００６】
本発明の目的は、高速で、僅かな演算経費により実施できる、１／ｆノイズの乱数の数列（乱数列）を生成するための方法を提供することである。本発明のほかの目的は、１／ｆノイズの対象となる（影響を受ける）技術的なシステムをシミュレーションするための改良された方法を提供することである。最後に、高速に実施でき、コンピュータシステムの資源をあまり必要としない、１／ｆノイズの乱数の数列を決定するためのコンピュータプログラムを有するコンピュータシステムも提供する必要がある。
【０００７】
この目的は、独立請求項の要旨により達成される。改良は、各従属請求項に記載されている。
【０００８】
本発明では、シミュレーションされるシステムをモデル化する場合のノイズシミュレーションの問題が、乱数数列の生成の問題に移行される。本発明に基づき、乱数の相関が決定される、このことは、相当する乱数数列を簡単かつ正確に生成するために使用される。
【０００９】
この場合、１／ｆノイズの乱数の少なくともひとつの数列を生成するための、本発明に基づく方法は、
スペクトル値βを決定（特定）する工程と、
強度定数ｃｏｎｓｔを決定する工程とをまず備えている。
これにより、シミュレーションされる１／ｆノイズの特徴が決定される。
【００１０】
次に、生成される、１／ｆノイズの乱数の数、およびシミュレーションに使用される経過変数（Ｌａｕｆｖａｒｉａｂｌｅ）ｎのための初期値を決定する。
【００１１】
長さｎの１つあるいはそれ以上のベクトルｙの、望ましい数の要素ｙ（ｎ）が、１／ｆ分散（分配、拡散）される（ｖｅｒｔｅｉｌｔｅｎ）乱数から演算されるまでの間ずっと、環状に繰り返される次のようなステップを本発明は備えている。次のようなステップとは、すなわち、
（１）経過変数ｎの実際の数値（実数）を１だけ増加（上昇）させるステップと、
（２）シミュレーション時間工程［ｔｎ−１；ｔｎ］を決定するステップと、
（３）以下の規則
Ｃｉｊ：＝ｃｏｎｓｔ・（−｜ｔｊ−ｔｉ｜β＋１＋｜ｔｊ−１−ｔｉ｜β＋１
＋｜ｔｊ−ｔｉ−１｜β＋１−｜ｔｊ−１−ｔｉ−１｜β＋１），
ｉ，ｊ＝１，．．．，ｎ
に基づいて、次元（寸法）（ｎ×ｎ）の分散行列（分散マトリックス）Ｃの要素Ｃｉｊを決定（特定）するステップと、
（４）分散行列Ｃを反転することにより、行列Ｃ−１を決定するステップと、
（５）ｓｑｒｔは、関数「平方根」を表し、ｅ（ｎ，ｎ）は、（ｎ，ｎ）によって表示される、反転された分散行列Ｃ−１の要素を表している以下の規則
σ＝ｓｑｒｔ（１／ｅ（ｎ，ｎ））
に基づき、変数（Ｇｒｏｅｓｓｅ）σを決定するステップと、
（６）長さｎのベクトルｘのｎ番目の構成要素を形成している、（０、１）正規分布乱数を決定するステップと、
（７）反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素と、先行している（ｎ−１）個のシミュレーション時間工程のために演算されている、ベクトルｙの（ｎ−１）個の要素とから、つまり、ｙ（ｎ−１）が、ベクトルｙの最初の（ｎ−１）個の構成要素を表しており、Ｃ−１・，ｎが、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素を表しており、Ｃ−１ｎ，ｎが、（ｎ，ｎ）によって表示される、反転された分散行列Ｃ−１の要素を表している以下の規則
【００１２】
【数４】
2004524607.tif 000003
【００１３】
に基づいて、変数μを形成するステップと、
（８）以下の規則
Ｙ（ｎ）＝ｘ（ｎ）＊σ＋μ
に基づいて、１／ｆ分散された乱数から長さｎのベクトルｙの要素ｙ（ｎ）を決定するステップとを含む方法。
【００１４】
なお、本明細書および図面におけるＣは、下記数式
【００１５】
【数５】
2004524607.tif 000004
【００１６】
を表し、Ｃ−１は下記数式
【００１７】
【数６】
2004524607.tif 000005
【００１８】
を表し、Ｃ−１・，ｎは下記数式
【００１９】
【数７】
2004524607.tif 000006
【００２０】
を表し、Ｃ−１ｎ，ｎは下記数式
【００２１】
【数８】
2004524607.tif 000007
【００２２】
を表し、Ｃｉｊは下記数式
【００２３】
【数９】
2004524607.tif 000008
【００２４】
を表す。
【００２５】
本発明に基づく方法を用い、技術的なシステムのシミュレーションを任意に延長できる。このため、既に生成された１／ｆ乱数が存在する場合、１／ｆ分散されている乱数を簡単な方法で生成できる。さらに、前もってシミュレーションされている時間間隔の結果をシミュレーションの基礎にできる。いわゆる、再起動能力は、シミュレーション実践のために非常に重要な特性である。このことは、ちょうど１／ｆノイズ源のためには、達成するのが難しい。なぜなら、ある時間間隔のための１／ｆノイズ源をシミュレーションする乱数は、既に数として（数により）決定されている、先行する時間間隔のための乱数に依存しているからである。本発明により、適応工程サイズ制御（ａｄａｐｔｉｖｅｎ Ｓｃｈｒｉｔｔｗｅｉｔｅｎｓｔｅｕｅｒｕｎｇ）を使用でき、このことにより、技術的なシステムをシミュレーションするための演算時間が著しく延長されることはない。このような適応工程サイズ制御は、シミュレーションされる技術システムのダイナミックを数によって決定する場合に、正確さ、および、演算時間効率を著しく上昇させる。
【００２６】
本発明に基づく方法では、シミュレーションする時間間隔を前もって設定しておく必要がもうない。変数の工程サイズを備えるだけで、実際のシステムダイナミックへ適応できる。このことは、シミュレーションの正確さを高める。
【００２７】
本発明は、１／ｆ分散された乱数の数列を、連続して、つまり、要素ごとに生成するための１方法を示している。この場合、この方法は、それぞれの新しく生成された乱数が、正しい方法で、確率論的な意味で、前もって生成されている１／ｆ分散された乱数に依存していることを保証する。このことにより、回路の数によるシミュレーションの間に、それぞれ必要な乱数を生成できる。
【００２８】
本発明は、１／ｆ分散された乱数を生成するために、条件付確率密度を使用する。この１／ｆ分散された乱数は、この乱数の確率論的な相関を、既に生成されており、前もって行われているシミュレーション工程のために必要な乱数によって保証する。
【００２９】
本発明の特に有効な形態では、１／ｆノイズの乱数のｑ数列が同時に演算される。この場合、環状に繰り返される以下の工程、すなわち、
（６）長さｎのベクトルｘのｎ番目の構成要素を形成している、（０、１）正規分布乱数を決定するステップと、
（７）反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素と、先行している（ｎ−１）個のシミュレーション時間工程のために演算されている、ベクトルｙの（ｎ−１）個の要素とから、つまり、ｙ（ｎ−１）が、ベクトルｙの最初の（ｎ−１）個の構成要素を表しており、Ｃ−１・，ｎが、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素を表しており、Ｃ−１ｎ，ｎが、（ｎ，ｎ）によって表示される、反転された分散行列Ｃ−１の要素を表している以下の規則
【００３０】
【数１０】
2004524607.tif 000009
【００３１】
に基づいて変数μを形成するステップと、
（８）以下の規則
Ｙ（ｎ）＝ｘ（ｎ）＊σ＋μ
に基づいて、１／ｆ分散された乱数から長さｎのベクトルｙの要素要素ｙ（ｎ）を決定するステップとの代わりに、
（９）長さｎ、ｋ＝１，・・・，ｑのベクトルｘｋのそれぞれ最後の構成要素を形成する、（０，１）正規分散された量ｑの乱数ｘｋ，ｎを決定するステップ（ただし、本工程では、ベクトルｘｋのそれぞれ最初の（ｎ−１）個の構成要素が、既に前もって演算されている）と、
（１０）ｙ（ｎ−１），ｋが、先行しているシミュレーション時間工程のために計算されている、ベクトルｙｋの最初の（ｎ−１）個の構成要素を表しており、Ｃ−１・，ｎが、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素を表しており、Ｃ−１ｎ，ｎが、反転された分散行列Ｃ−１の、（ｎ，ｎ）により表示される要素を表している、以下の規則
【００３２】
【数１１】
2004524607.tif 000010
【００３３】
に基づいて、ｋ＝１，．．．，ｑに関して実施される、ｑ個の変数μｋを形成するステップと、
（１１）１／ｆ分散された乱数から、つまり、以下の規則、
ｙｋ，ｎ＝ｘｋ，ｎ＊σ＋μｋ（ただしｋ＝１，．．．，ｑ）
に基づいて、長さｎのベクトルｙｋの各ｎ番目の構成要素を形成しているｑの要素ｙｋ，ｎを演算するステップとが備えられている。
【００３４】
１／ｆ分散された乱数のうち長さｎのｑのベクトルｙｋ（ｋ＝１，．．．，ｑ）は、行列ＮＯＩＳＥに特に有効に配置されている。この行列ＮＯＩＳＥは、シミュレーションにおいて、シミュレーションするシステムの１／ｆノイズ影響を示している。
【００３５】
本発明に基づき以下のような筋道の考え方（Ｇｅｄａｎｋｅｎｇａｎｇ）は、１／ｆノイズのシミュレーションのためのコンセプトを基礎としている。確率論的な影響を受けるシステムのダイナミックは、確率過程によって十分にモデル化される。このようなシステムダイナミックをシミュレーションするために、一般的に、基礎となっている確率過程の個々のランダム実現（いわゆる経路（Ｐｆａｄｅ））において、数によって演算されている。１／ｆノイズ源を有するシステムをシミュレーションするために、
【００３６】
【数１２】
2004524607.tif 000011
【００３７】
の形状の確率積分の経路を数によって演算することが有効である。この場合、ｓ（積算数）およびｔ（集積上限）は、時間、η１／ｆ（ｓ）ｄｓは、１／ｆノイズ源、Ｙ（ｓ）は、変数（例えば、回路シミュレーションにおける電圧）の時間的なダイナミックを表す確率過程を示している。
【００３８】
微分（数学的には：超関数の意味での微分）が、１／ｆノイズプロセスη１／ｆ（ｓ）を生じる確率過程を、ＢＦＢＭ（ｓ）によって表す場合、演算される確率積分は、
【００３９】
【数１３】
2004524607.tif 000012
【００４０】
と記載できる。
【００４１】
右側の積分は、積分器（Ｉｎｔｅｇｒａｔｏｒ、積分要素）としてプロセスＢＦＢＭ（ｓ）を含む確率過程Ｙ（ｓ）のリーマン―スチルテス積分と解釈できる。この積分器は、０≡ｔ０＜ｔ１＜．．．＜ｔ≡ｔに基づく積分間隔［０，ｔ］を、ｎ個に分散された部分間隔［ｔｉ，ｔｉ−１］，ｉ＝１，．．．ｎ，に分解することにより、合計によって概算される。
【００４２】
【数１４】
2004524607.tif 000013
【００４３】
この合計は、ランダム変数である。確率的実験の結果ωへの依存は、常に省略されている。
【００４４】
プロセスＢＦＢＭ（ｓ）の普遍化される微分は、１／ｆスペクトルを備え、このプロセスＢＦＢＭ（ｓ）は、文献では、「分数ブラウン運動（Ｆｒａｃｔｉｏｎａｌ Ｂｒｏｗｎｉａｎ Ｍｏｔｉｏｎ）」の名称で知られている。ＢＦＢＭ（ｓ）は、ガウス確率過程であり、このような過程として、その予測値
【００４５】
【数１５】
2004524607.tif 000014
【００４６】
および、その共分散関数
【００４７】
【数１６】
2004524607.tif 000015
【００４８】
により、完全に特徴付けられている。
【００４９】
適切な乱数を必要に応じて生成するための本発明に基づく方法では、１／ｆノイズ影響のシミュレーションが、乱数［ＢＦＢＭ（ｔｉ）−ＢＦＢＭ（ｔｉ−１）］の実現を発生させること（Ｅｒｚｅｕｇｕｎｇ ｖｏｎ Ｒｅａｌｉｓｉｅｒｕｎｇｅｎ）、つまり分数ブラウン運動の増進を発生させるること（ｇｅｎｅｒａｔｉｏｎ ｏｆ ｉｎｃｒｅｍｅｎｔｓ）に基本的に由来している。
【００５０】
本発明により、ランダム変数ΔＢＦＢＭ（ｉ）を、オンラインで必要に応じて実現でき、つまり、システム方程式の連続した積分の経過の間に生成できる。このことから、結果として以下の２つの方法が必要である。
【００５１】
（ａ）乱数の数列の長さｎ
｛ΔＢＦＢＭ（１），．．．，ΔＢＦＢＭ（ｎ）｝は、シミュレーション経過の間、変数のままである必要がある。特に、いつの時点でもシミュレーションを延長できる必要がある（再起動性能）。このことは、方法の性能を意味しており、この性能は、このために必要とされる付加的な乱数を、正しい方法で既に生成されている部分数列（部分列、サブ連続）と相関するように生成する。
【００５２】
（ｂ）ｔｉを、シミュレーション経過の間に、実際に達成される時間とする。そうすると、時間間隔［ｔｉ，ｔｉ＋１］、つまり、次の積分工程の工程サイズを、目下のシステムダイナミックから、つまり、適応させて、決定できる必要がある。
【００５３】
本発明は、｛ΔＢＦＢＭ（１），．．．，ΔＢＦＢＭ（ｎ）｝の生成、つまり、乱数の数列を、連続してつまり要素ごとに、どのようにして生成できるかという規則を示すことにより、２つの必要性に適合する。この場合、工程サイズΔｔｉ：＝ｔｉ−ｔｉ−１を各新しい乱数のために自由に選択できる。
【００５４】
まず、第１工程では、いわゆる、「条件付密度」のための手がかりが調査される。
【００５５】
また、第１工程では、ランダム変数ベクトル（ΔＢＦＢＭ（１），．．．，ΔＢＦＢＭ（ｎ））の分布（Ｖｅｒｔｅｉｌｕｎｇ）が考慮される。
【００５６】
個々のランダム変数ΔＢＦＢＭ（ｉ）は、ガウス確率過程の増加であるので、ランダム変数ベクトル（ΔＢＦＢＭ（１），．．．，ΔＢＦＢＭ（ｎ））は、ｎ次元のガウス分散されたランダム変数であり、すなわち、その（ｎ次元の）予測値Ｅおよびその分散行列Ｃにより完全に決定されている。２つの変数は、定式（１．３）および（１．４）から、
Ｅ（ΔＢＦＢＭ（ｉ））＝０，ｉ＝１，．．．，ｎ      （３．５）
【００５７】
【数１７】
2004524607.tif 000016
【００５８】
であると演算される。
この場合、本発明に基づくオンライン方法を、完全な帰納法のやり方で示す必要がある。
【００５９】
帰納法の始まり、つまり、方法の出発点は、予想値０および変数
【００６０】
【数１８】
2004524607.tif 000017
【００６１】
を有する実数値のガウス分散を実現することである。
【００６２】
帰納終結（Ｉｎｄｕｋｔｉｏｎｓｓｃｈｌｕｓｓｅｓ）の目的で、完全な（ΔＢＦＢＭ（１），．．．，ΔＢＦＢＭ（ｎ））の実現（Ｒｅａｌｉｓｉｅｒｕｎｇ）の生成（ｅｒｇｉｂｔ）のために、ΔＢＦＢＭ（ｎ）の実現によって、どのようにして（ΔＢＦＢＭ（１），．．．，ΔＢＦＢＭ（ｎ−１））の実現を拡張（ｅｒｗｅｉｔｅｒｎ）するかを示す必要がある。書き方を簡単にするために、乱数の既に「２乗されている」部分数列を、（ｙ１，．．．，ｙｎ−１）＝：ｙ（ｎ−１）Ｔで示し、ΔＢＦＢＭ（ｎ）のまだ２乗されていない実現をｙｎによって示す。
【００６３】
この場合、問題を以下のように定式化できる。
ｎ次元を、分散行列Ｃを有し、平均値の定まっていない（ｍｉｔｔｅｌｗｅｒｔｆｒｅｉ）ガウスランダム変数Ｚとする。Ｚの実現の最初のｎ−１個の要素は、乱数ベクトルｙ（ｎ−１）の形状に既に２乗されており、既知である。このとき見出されるのは、分散であり、この分散から、Ｚの実現化ｙ＝（ｙ（ｎ−１），ｙｎ）のためのｙ（ｎ−１）を完了するため、ｎ番目の要素ｙｎを引き出す必要がある。
【００６４】
本発明の解決法は、（ｙ（ｎ−１）が既に決定されているという条件下で）ｙｎのための条件付確率密度ｆ（ｙｎ｜ｙ（ｎ−１））を考える場合に見出せる。この変数は、ガウス正規分布が上記の場合には、
【００６５】
【数１９】
2004524607.tif 000018
【００６６】
と演算される。
この場合、変数Ｃ−１ｎ，ｎが、以下のような反転された分散行列Ｃ−１：
【００６７】
【数２０】
2004524607.tif 000019
【００６８】
の書き方から生じる。ただし、
【００６９】
【数２１】
2004524607.tif 000020
【００７０】
【数２２】
2004524607.tif 000021
【００７１】
【数２３】
2004524607.tif 000022
【００７２】
である。
【００７３】
変数μは、
【００７４】
【数２４】
2004524607.tif 000023
【００７５】
を表している。
条件付密度ｆ（ｙｎ｜ｙ（ｎ−１））は、平均値μおよび変数１／Ｃ−１ｎ，ｎを有する、ガウス正規分布の確率密度でもある。
【００７６】
上記のような変数が存在するように、Ｃ−１ｎ，ｎ≠０が有効である必要がある。このことは、以下の論証に基づき保証される。
Ｃ−１およびＣ−１は、同一固有方向（ｓｅｌｂｅｎ Ｅｉｇｅｎｒｉｃｈｔｕｎｇｅｎ）と反転固有値（ｉｎｖｅｒｓｅ Ｅｉｇｅｎｗｅｒｔｅ）とを有している。従って、その結果、行列Ｃ−１の固有値０は、ランダム変数ベクトル（ΔＢＦＢＭ（１），．．．，ΔＢＦＢＭ（ｎ））の異なる変数を有する可能性もある。それゆえ、Ｃ−１の全ての固有値が、ゼロに等しくないという条件を付けられる。Ｃ−１の固有値は、あらゆる場合に負ではないので、従って：行列Ｃ−１は、対称であり正に定義されている、ということが有効である。座標軸の名前の付け替えにより、この行列を、形式（３．８）から以下の形式にできる。
【００７７】
【数２５】
2004524607.tif 000024
【００７８】
この行列は、その構造の意味するように、同じく対称であり、正に定義されている。対称であり、正に定義されている行列のためのジルベスタ条件（Ｓｙｌｖｅｓｔｅｒ−Ｋｒｉｔｅｒｉｕｍｓ）に基づき、
【００７９】
【数２６】
2004524607.tif 000025
【００８０】
が推定され、定理が示される。本発明に基づく方法によって、１／ｆノイズ源のシミュレーションは、ガウス分散されている乱数の生成に由来している。
【００８１】
要求される方法により、既に生成されている数列ｙ（ｎ−１）と相関する乱数ｙｎを生成するために、反転された分散行列Ｃ−１（ｎ×ｎ行列）が必要である。厳密に言えば、この行列のｎ番目の行を認識すること、つまり、（（Ｃ−１・，ｎ）Ｔ，Ｃ−１ｎ，ｎ）を認識することだけが必要である。形式（３．６）から読み取れるように、分散行列Ｃは、分散された部分間隔（工程サイズ）［ｔｉ−１，ｔｉ］におけるシミュレーション間隔［０，ｔｎ］の分解に応じている。特に、Ｃの最後の列は、（最後の行と同一のＣとの対称が原因で）ｔｎ、すなわち、実際の工程サイズΔｔｎ−ｔｎ−ｔｎ−１に応じている。
【００８２】
ｎＸｎ分散行列Ｃの左辺上側（ｎ−１）×（ｎ−１）部分行列
【００８３】
【数２７】
2004524607.tif 000026
【００８４】
は、正確に長さｎ−１の乱数数列のための分散行列である。この分散行列は、ｙ（ｎ−１）の演算（もしくは最後の要素ｙ（ｎ−１）の演算）のために既に決定されており、反転されている必要がある。従って、方法をより迅速にするため、例えば、シューア補集合（Ｓｃｈｕｒ−Ｋｏｍｐｌｅｍｅｎｔｓ）を用いて、行列反転に増加の方法を再び利用できる。
【００８５】
本発明は、１／ｆノイズの影響を受ける技術的なシステムをシミュレーションするための方法においても実現されている。この場合、モデル化する場合、および／または、システムの入力部チャンネルに存在する変数を決定する場合に、乱数が使用される。この乱数は、本発明の方法に基づいて決定されている。
【００８６】
同じように、コンピュータシステムおよび／またはコンピュータプログラムが、１／ｆノイズの乱数の数列を決定するため、または、本発明に基づくほかの方法を実施するために備えられている。本発明は、このようなコンピュータプログラムを有するデータ記録媒体（データキャリア）においても実現されている。さらに、本発明は、本発明に基づくコンピュータプログラムを、例えば、データ網と接続されているコンピュータにおけるインターネットからのように、電子デー網からダウンロードするという方法において実現されている。
【００８７】
本発明を実施例を参考に、図で説明する。
【００８８】
図１は、シミュレーションされる技術的なシステムの概略図である。図２は、１／ｆノイズの乱数の数列を決定するための組織図である。図３は、図３ａ〜３ｆに基づく、第１シミュレーション時間工程のための演算例を示す図である。図４は、図４ａ〜４ｆに基づく、第２シミュレーション時間工程のための演算例を示す図である。図５は、図５ａ〜５ｆに基づく、第３シミュレーション時間工程のための演算例を示す図である。
【００８９】
図１は、シミュレーションされる必要があり、ノイズを含むシステムの概略図を示す。
【００９０】
このシステムは、システム特性を表す、四角で示されるシステムモデル１によって示されている。このシステム特性は、ベクトルＩＮＰＵＴとしても示される入力部チャンネル２、および、ＯＵＴＰＵＴとしても示される出力部チャンネル３から生じている。さらに、システムに条件付けられているノイズが含まれている。このノイズは、ノイズ入力部チャンネル４に存在し、ベクトル、もしくは、行列ＮＯＩＳＥとしても表される。従って、この行列ＮＯＩＳＥは、ノイズを複数のチャンネルについて考慮する場合に存在し、このとき、行列ＮＯＩＳＥの各列は、ノイズ入力部チャンネルに存在しているノイズ値のベクトルを含んでいる。
【００９１】
ノイズ入力部チャンネル４におけるノイズを、ノイズに条件付けられている、システムモデル１の変化と解釈することが好ましい。
【００９２】
入力部チャンネル２と出力部チャンネル３との特性を、微分方程式のシステム、または、代数微分方程式のシステムによって説明できる。その結果、システム特性の信頼性のある予測が可能である。
【００９３】
図１に示すシステムのシミュレーションの各時間工程のときに、入力部チャンネル２に存在するベクトルＩＮＰＵＴと、ノイズ入力部チャンネル４に存在するベクトルＮＯＩＳＥとのために、出力部チャンネル３のベクトルＯＵＴＰＵＴが演算される。
【００９４】
合理的には、より長い期間にわたるシミュレーションのために、ベクトルＩＮＰＵＴ，ＯＵＴＰＵＴ，ＮＯＩＳＥは、行列として示され、この場合、当該行列の各一つの列ｋは、当該ＩＮＰＵＴ，ＯＵＴＰＵＴ，ＮＯＩＳＥの対応する時間数列（Ｚｅｉｔｒｅｉｈｅ）の値を含んでいる。
【００９５】
図２は、システムモデル１のノイズ入力部チャンネル４のための行列ＮＯＩＳＥの列ｋを形成している各ベクトルｙｋを、どのようにして求めるかということを具体的に示している。各ベクトルｙｋを、ノイズ源のシミュレーションに使用する。
【００９６】
第１工程では、望ましいスペクトル値βおよび強度定数ｃｏｎｓｔを決定する。さらに、実際のシミュレーション時間間隔の計数器ｎを０に設定する。
【００９７】
さて、各シミュレーション時間工程のために、演算工程を以下の順序で連続して実施する。
【００９８】
まず、実際のシミュレーション時間工程を決定する。これと同じく、実際のシミュレーション時間工程の終了も決定する。このことにより、次の考察時点（Ｂｅｔｒａｃｈｔｕｎｇｓｚｅｉｔｐｕｎｋｔ）が生じる。
【００９９】
その後、実際のシミュレーション時間工程の計数器の計数が１だけ増える。
【０１００】
続いて、次元（ｎ×ｎ）の分散行列Ｃを、等式（３．６）に基づいて決定する。
【０１０１】
その後、例えば、コレスキー分解を用いる行列Ｃの反転の工程が続く。このとき、例えば、シュール補足技術（Ｓｃｈｕｒｋｏｍｐｌｅｍｅｎｔ−Ｔｅｃｈｎｉｋｅｎ）を使用する場合、効率を上げるために先行する工程の反転行列もアクセスされる。
【０１０２】
次に、変数σを、定式
σ＝ｓｑｒｔ（１／ｅ（ｎ，ｎ））
により演算する。ただし、ｓｑｒｔは、関数「平方根」を表し、ｅ（ｎ，ｎ）は、（ｎ，ｎ）によって表示される、反転された分散行列Ｃ−１の要素を表している。
【０１０３】
さらに、（０，１）正規分散されたランダム変数Ｘｋの値が引き出され、これにより、正規分散された乱数のベクトルｘｋを補足する。引き出された乱数は、予測値０および変数１を備えている。この工程を、シミュレーションされる各ノイズ源のために実施する。
【０１０４】
さらに、変数μｋを形成する。この変数は、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素と、先行する（ｎ−１）個のシミュレーション時間工程のために演算されており、１／ｆ分散されている、（ｎ−１）個の乱数の数列とから形成される。このため、定式（３．９）に基づく処理を行う。この工程を、シミュレーションされる各ノイズ源ｋのために実施する。
【０１０５】
最後に、列インデックスｋがシミュレーションされるノイズ源を示しており、時間インデックスがｎに等しい、行列ＮＯＩＳＥの要素を演算する。これにより、実際のシミュレーション時間工程が表示される。実際に演算された、行列ＮＯＩＳＥの要素ｒ（ｋ，ｎ）は、乱数である。この乱数は、ＮＯＩＳＥの同じ列ｋの要素であって、上にある（ｎ−１）個の要素と共に、長さｎのベクトルｙｋを、１／ｆ分散された乱数から形成している。このベクトルｙｋは、最初のｎ個のシミュレーション時間工程のためのノイズ源のうちの１つをシミュレーションするのに使用される。
【０１０６】
次に、ＮＯＩＳＥのｎ番目の行の各要素ｙｋを、等式（３．７）〜（３．９）に基づき、ベクトルｘの最後の乱数ｘｋ、ならびに、変数μｋおよびσから決定する。つまり、ｙｋを以下の規則：
ｙｋ＝ｘｋ＊σ＋μｋ
に基づいて決定する。
【０１０７】
図３〜図５に、具体的な演算結果を出す実施例を示す。
【０１０８】
この場合、スペクトル値βの値は、常に０，５とする。強度ｃｏｎｓｔの値は、任意に１，０とする。３つのノイズ源のシミュレーションに応じて、それぞれ３つの乱数が同時に処理される。これら３つのノイズ源は、別々のチャンネルにおいて、シミュレーションされるシステムに同時に影響を与え、それぞれ１つのベクトルｙｋに配置されている。ただし、ｋは、１〜３の整数値である。
【０１０９】
図３は、そのサブ図３ａ〜３ｆに基づく、第１シミュレーション時間工程［ｔ０、ｔ１］＝［０、０，５］のための演算例を示す。
【０１１０】
図３ａは、シミュレーション工程サイズ（Ｓｉｍｕｌａｔｉｏｎｓ−Ｓｃｈｒｉｔｔｗｅｉｔｅ）の場合に、１／ｆ分散されている乱数を生成するための、次元１×１の分散行列Ｃを示す。このとき、Ｃは、値０，７０を有する単なる（ただ１つの）スカラーである。なぜなら、Ｃ（１，１）―つまり、ただしｉ＝ｊ＝１―は、等式（３．６）を使用して、
１，０・（−｜ｔ１−ｔ１｜０，５＋１＋｜ｔ１−１−ｔ１｜０，５＋１
＋｜ｔ１−ｔ１−１｜０，５＋１＋｜ｔ１−１−ｔ１−１｜０，５＋１）＝
０＋０，５１，５＋０，５１，５−０＝０，７０７１０６．．．
となる。
【０１１１】
図３ｂは、図３ａの分散行列Ｃの反転を示す。このことは、ここでは詳しく説明していないコレスキー分解により行われたものである。（Ｃ  Ｃ−１）＝（０，７０７１０６．．．０，７０７１０６．．．−１）の検証により、正しい値１が生じる。このことは、Ｃ（１，１）のための値の正しさを証明している。
【０１１２】
図３ｃは、第１シミュレーション工程ｎ＝１のための変数σを示している。この変数は、等式
σ＝ｓｑｒｔ（１／０，７０７１０６．．．）
により算出される。ただし、ｓｑｒｔは、平方根を、ｅ（１，１）は、（１，１）によって表示される、反転された分散行列Ｃ−１の要素０，７０７１０６．．．を表している。
【０１１３】
図３ｄは、各１つのシミュレーションされるノイズ源のための、（０，１）正規分布されたランダム変数Ｘｋの３つの値ｘ１、ｘ２、ｘ３を示している。これらの値は、正規分散された乱数のベクトルｘｋ毎に第１要素を形成する。引き出された乱数は、予想値０および変数１を備えている。
【０１１４】
図３ｅは、３つのシミュレーションされるノイズ源のぞれぞれのための３つの変数μｋを示している。この変数μｋは、定式（３．９）に基づき、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素と、先行する（ｎ−１）個のシミュレーション時間工程のために演算されており、１／ｆ分散されている、（ｎ−１）片の乱数の数列とから生じる。第１シミュレーション工程では、これら２つのベクトルが、長さ０をそれぞれ有している。従って、第１シミュレーション工程における全ての変数μｋのために、：μｋ＝０が生じる。
【０１１５】
図３ｆは、１／ｆ分散されている乱数を有する長さ１の３つのベクトルｙｋを示している。これら３つのベクトルｙｋは、１／ｆ分散されている３つのノイズ源の特性を、第１シミュレーション時間工程［０，ｔ１］＝［０，０．５］のためにシミュレーションする。行列ＮＯＩＳＥが、これら３つのベクトルｙｋから生じる。ただし、値ｋは、１〜３の整数値である。ＮＯＩＳＥの第１行目の各要素ｙｋは、等式（３．７）〜（３．９）により、関係するベクトルｘの最後の乱数ｘｋ、ならびに、変数μｋおよびσ１から、以下の規則に基づいて決定される。例として、以下に、第１ベクトルｙ１の第１要素ｙ１（ｎ＝１）を演算する：
ｙ１（ｎ＝１）＝ｘ１（ｎ＝１）＊σ＋μ１＝
＝−０，３５．．．＊０，８４．．．＋０，００．．．＝
＝−０，３０．．．；
ｙ２（ｎ＝１）およびｙ３（ｎ＝１）は、これと同じように演算される。
【０１１６】
図４は、そのサブ図４ａ〜４ｆに基づく、第２シミュレーション時間工程［ｔ１，ｔ２］＝［０，５、０，７５］のための演算例を示す。第２シミュレーション時間工程のためのｎの値は、常に２に等しい。
【０１１７】
図４ａは、ノイズ源毎に各１つのほかの乱数を生成するために必要な次元（ｎ×ｎ）＝２ｘ２の分散行列Ｃを示す。このように新しく生成される乱数は、図３ｆに基づく結果と共に、長さ２のベクトルｙｋを、１／ｆ分散された乱数から形成する。この場合、各ノイズ源は、ベクトルｙｋを生成する。この場合、分散行列Ｃは、等式（３．６）に基づいて決定される。
【０１１８】
例えば、要素Ｃ（２，１）―ただしｉ＝２、ｊ＝１―においてこのことが行われる。等式（３．６）を使用して、Ｃ（２，１）が、
１，０・（−｜ｔ１−ｔ２｜０，５＋１＋｜ｔ１−１−ｔ２｜０，５＋１＋｜ｔ１−ｔ２−１｜０，５＋１−｜ｔ１−１−ｔ２−１｜０，５＋１）＝
＝（−｜０，５−０，７５｜０，５＋１＋｜０−０，７５｜０，５＋１＋｜０，５−０，５｜０，５＋１−｜０−０，５｜０，５＋１）＝  −０，１２５＋０，６４９５．．＋０−０，３５３５．．．＝０，１７０９．．．；
と算出される。
【０１１９】
図４ｂは、図４ａの分散行列Ｃの反転を示している。条件（Ｃ  Ｃ−１）の検証（ここには示さず）は、次元２×２の行列を生じる。この場合、（１，１）および（２，２）により表示される要素は、１に等しく、他の要素は、値０を有している。
【０１２０】
図４ｃは、変数σを示している。この変数σは、工程４ｂの反転された分散行列Ｃ−１から演算される。変数σは、
σ＝ｓｑｒｔ（１／ｅ（ｎ，ｎ）＝ｓｑｒｔ（１／ｅ（２，２））＝
＝ｓｑｒｔ（１／４，７９．．．）＝０，４５．．．
として生じる。ただし、ｓｑｒｔは、平方根を表し、ｅ（２，２）は、（２，２）によって表示される、図４ｂの反転された分散行列Ｃ−１の要素を表している。
【０１２１】
図４ｄは、独立している（０，１）正規分布乱数の３つのベクトルｘｋを示している。この際、ベクトルｘｋは、長さ２をそれぞれ有している。シミュレーションされるノイズ源毎に、（０，１）正規分布されている１つの乱数ｘｋが引き出される。この引き出された乱数は、予測値０と変数１とをそれぞれ備えている。これにより、図３ｄの正規分散されている乱数のベクトルｘｋが補足され、その結果、図４ｄの正規分布されている乱数のベクトルｘｋが生じる。
【０１２２】
図４ｅは、３つの変数μｋを示している。このμｋは、工程４ｂに基づく反転された分散行列Ｃ−１および工程３ｆに基づく３つの乱数から演算されている。各シミュレーションされるノイズ源のために、変数μｋが、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素と、先行している（ｎ−１）個のシミュレーション時間工程のために演算されており、１／ｆ分散されている、（ｎ−１）片の乱数の数列とから演算される。第２シミュレーション工程では、変数μｋが、Ｃ−１の第２行目の第１構成要素と、ベクトルｙｋの第１構成要素とから演算される。例として、これを、値μ１に基づき実施する。
【０１２３】
【数２８】
2004524607.tif 000027
【０１２４】
図４ｆは、１／ｆ分散されている乱数を有する長さ２の３つのベクトルｙｋを示している。これら３つのベクトルｙｋは、１／ｆ分散されている３つのノイズ源の特性を、第２シミュレーション時間工程［ｔ１，ｔ２］＝［０，５、０，７５］のためにシミュレーションする。行列ＮＯＩＳＥが、これら３つのベクトルｙｋから生じる。ただし、値ｋは、１〜３の整数値である。ＮＯＩＳＥの第２行目の各要素ｙｋは、等式（３．７）〜（３．９）により、関係するベクトルｘの最後の乱数ｘｋ、ならびに、変数μｋおよびσ１から、以下の規則に基づいて決定される。例として、以下に、第１ベクトルｙ１の第２要素ｙ１（ｎ＝２）を演算する。
【０１２５】
ｙ１（ｎ＝２）＝ｘ１（ｎ＝２）＊σ＋μ１＝
＝−０，３９．．．＊０，４５．．．−０，０７．．．＝
＝−０，１０．．．；
図５は、その図５ａ〜５ｆに基づく、第３シミュレーション時間工程［ｔ２、ｔ３］＝［０，７５、１，２５］のための演算例を示す。第３シミュレーション時間工程の間のｎの値は、常に３に等しい。
【０１２６】
図５ａは、ノイズ源毎に各１つのほかの乱数を生成するために必要な次元（ｎ×ｎ）＝３×３の分散行列Ｃを示す。このように新しく生成される乱数は、図４ｆに基づく結果と共に、長さ３のベクトルｙｋを、１／ｆ分散された乱数から形成する。この場合、各ノイズ源は、ベクトルｙｋを生成する。この場合、分散行列Ｃは、等式（３．６）に基づいて決定される。
【０１２７】
例として、要素Ｃ（３，１）―ただしｉ＝３、ｊ＝１―についてこのことが行われる。等式（３．６）を使用して、Ｃ（３，１）が、
１，０・（−｜ｔ１−ｔ３｜０，５＋１＋｜ｔ１−１−ｔ３｜０，５＋１＋｜ｔ１−ｔ３−１｜０，５＋１−｜ｔ１−１−ｔ３−１｜０，５＋１）＝
＝（−｜０，５−１，２５｜０，５＋１＋｜０−１，２５｜０，５＋１＋｜０，５−０，７５｜０，５＋１−｜０−０，７５｜０，５＋１）＝
−０．６４９５．．．＋１，３９５７．．＋０，１２５−０，６４９５．．．＝０，２
２．．．
と算出される。
【０１２８】
図５ｂは、図５ａの分散行列Ｃの反転Ｃ−１を示している。条件（Ｃ  Ｃ−１）の検証（ここには示さず）は、次元３×３の行列を生じる。この場合（１，１）、（２，２）および（３，３）により表示される要素は、１に等しく、他の要素は、値０を有している。
【０１２９】
図５ｃは、変数σを示している。この変数σは、工程５ｂの反転された分散行列Ｃ−１から演算される。変数σは、
σ＝ｓｑｒｔ（１／ｅ（ｎ，ｎ））＝ｓｑｒｔ（１／ｅ（３，３））＝
＝ｓｑｒｔ（１／１，７５．．．）＝０，７５．．．
として生じる。ただし、ｓｑｒｔは、平方根を表し、ｅ（３，３）は、（３，３）によって表示される、図５ｂの反転された分散行列Ｃ−１の要素を表している。
【０１３０】
図５ｄは、独立している（０，１）正規分布乱数の３つのベクトルｘｋを示している。この際、ベクトルｘｋは、長さ３をそれぞれ有している。シミュレーションされるノイズ源毎に、（０，１）正規分散されている１つの乱数ｘｋが引き出される。この引き出された乱数は、予測値０と変数１とをそれぞれ備えている。これにより、図４ｄの正規分散されている乱数のベクトルｘｋが補足され、その結果、図５ｄの正規分散されている乱数のベクトルｘｋが生じる。
【０１３１】
図５ｅは、３つの変数μｋを示している。このμｋは、工程５ｂに基づく反転された分散行列Ｃ−１および工程４ｆに基づく３つの乱数から演算されている。各シミュレーションされるノイズ源のために、変数μｋが、反転された分散行列Ｃ−１のｎ番目の行の最初の（ｎ−１）個の構成要素と、先行している（ｎ−１）個のシミュレーション時間工程のために定式（３．９）に基づいて演算されており、１／ｆ分散されている、（ｎ−１）片の乱数の数列とから演算される。第２シミュレーション工程では、変数μｋが、Ｃ−１の第３行目の最初の２つの構成要素と、ベクトルｙｋの最初の２つの構成要素とから演算される。例として、これを、値μ１に基づき実施する。
【０１３２】
【数２９】
2004524607.tif 000028
【０１３３】
図５ｆは、１／ｆ分散乱数を有する長さ３の３つのベクトルｙｋを示している。これら３つのベクトルｙｋは、１／ｆ分散されている３つのノイズ源の特性を、第３シミュレーション時間工程［ｔ２，ｔ３］＝［０，７５、１，２５］のためにシミュレーションする。行列ＮＯＩＳＥが、これら３つのベクトルｙｋから生じる。ただし、値ｋは、１〜３の整数値である。ＮＯＩＳＥの第３行目の各要素ｙｋ（ｎ＝３）は、等式（３．７）−（３．９）により、関係するベクトルｘの最後の乱数ｘｋ（ｎ＝３）、ならびに、変数μｋおよびσ１から、以下の規則に基づいて決定される。例として、以下に、第１ベクトルｙ１の第３要素ｙ１（ｎ＝３）を演算する。
【０１３４】
ｙ１（ｎ＝３）＝ｘ１（ｎ＝３）＊σ＋μ１＝
＝−０，９０．．．＊０，７５．．．＋０，００．．．＝
＝−０，６７．．．；
上記演算例を具体的に実施するため、さらに以下の条件が考えられる。
【０１３５】
図３，図４および図５に示す数値は、図２を参照して説明された第１、第２および第３シミュレーション間隔のための計算工程の中間結果および最終結果を再現している。このとき、正確な数の演算に基づく全ての値は、これら結果をよりよく再現できるように、小数点第２位以下が切り捨てられている。従って、実施例を正確に理解するには、上記ｘベクトルから始めて、上記ｙベクトルにまで達成するように、図に示されている中間値によってではなく、正確な中間値によってさらに計算する必要がある。
【０１３６】
図３ｃ、図４ｃおよび図５ｃに、（０，１）正規分散ランダム変数の各ベクトルを示す。この場合、ランダム変数は、それぞれノイズ源である。簡易化するため、予測値０と変数１とを有するこのような乱数を、どのようにして得るかをここには示していない。このことは、当業者に周知である。
【図面の簡単な説明】
【図１】
シミュレーションされる技術的なシステムの概略図である。
【図２】
１／ｆノイズの乱数の数列を決定するための組織図である。
【図３】
図３ａ〜３ｆに基づく、第１シミュレーション時間工程のための演算例を示す図である。
【図４】
図４ａ〜４ｆに基づく、第２シミュレーション時間工程のための演算例を示す図である。
【図５】
図５ａ〜５ｆに基づく、第３シミュレーション時間工程のための演算例を示す図である。
【図１】
2004524607.tif 000031
【図２】
2004524607.tif 000032
【図３】
2004524607.tif 000033
【図４】
2004524607.tif 000034
【図５】
2004524607.tif 000035
【国際公開パンフレット】
2004524607.tif 000036
2004524607.tif 000037
2004524607.tif 000038
2004524607.tif 000039
2004524607.tif 000040
2004524607.tif 000041
2004524607.tif 000042
2004524607.tif 000043
2004524607.tif 000044
2004524607.tif 000045
2004524607.tif 000046
2004524607.tif 000047
2004524607.tif 000048
2004524607.tif 000049
2004524607.tif 000050
2004524607.tif 000051
2004524607.tif 000052
2004524607.tif 000053
2004524607.tif 000054
2004524607.tif 000055
2004524607.tif 000056
2004524607.tif 000057
2004524607.tif 000058
2004524607.tif 000059
2004524607.tif 000060
2004524607.tif 000061
2004524607.tif 000062
2004524607.tif 000063
2004524607.tif 000064
2004524607.tif 000065
2004524607.tif 000066
2004524607.tif 000067
2004524607.tif 000068
2004524607.tif 000069
2004524607.tif 000070
2004524607.tif 000071
2004524607.tif 000072
2004524607.tif 000073
2004524607.tif 000074
2004524607.tif 000075
2004524607.tif 000076
【国際公開パンフレット（コレクトバージョン）】
2004524607.tif 000077
2004524607.tif 000078
2004524607.tif 000079
2004524607.tif 000080
2004524607.tif 000081