(57)【要約】【課題】  実測値に基づく伝送線路のモデル化に際し、伝送線路特性を計測する特定の計測器が有する固有の誤差、ノイズ等の影響を抑制し、モデルを実際の伝送線路特性にフィッティングさせるために、従来のような試行錯誤によるモデル化作業の長時間化を解消した伝送線路特性のモデル化方法およびそのモデル化装置を提供することにある。【解決手段】  実測値と伝送線路特性のシミュレーション結果とを比較しながら、伝送線路特性を実測値に近づけるフィッティング手法を用いて伝送線路特性の決定を行う。【選択図】  図１
2008116415.tif 000002
【特許請求の範囲】
【請求項１】
伝送線路の複数の伝送特性をそれぞれ測定し、伝送線路モデルの計算結果と、前記伝送線路の複数の測定結果とを同時に比較しながらフィッティングを行う伝送線路のモデル化方法において、多目的最適化手法を用いて複数のパレート最適解を求めることを特徴とする伝送線路のモデル化方法。
【請求項２】
前記複数の伝送特性として、周波数特性と時間軸応答特性を測定し、周波数特性のフィッティングと、時間軸応答特性のフィッティングを同時に行うことを特徴とする請求項1記載の伝送線路のモデル化方法。
【請求項３】
前記周波数特性のフィッティングにおいて、伝送線路の特性を実測した特性と伝送線路モデルから計算した伝送特性の結果とを比較し、その振幅特性と位相特性それぞれの誤差が最小となるように伝送線路モデルを調整することを特徴とする請求項２記載の伝送線路のモデル化方法。
【請求項４】
前記時間軸応答特性のフィッティングにおいて、インパルス波形を伝送線路に通して実測した受信波形と、伝送線路モデルから算出したインパルス応答波形の誤差を求め、その誤差が最小となるように伝送線路モデルを調整することを特徴とする請求項２記載の伝送線路のモデル化方法。
【請求項５】
前記伝送線路モデルとして、表皮効果の周波数特性を考慮した伝搬定数γを用いることを特徴とする請求項１記載の伝送線路のモデル化方法。
【請求項６】
前記多目的最適化手法として、遺伝的アルゴリズムを用いることを特徴とする請求項1乃至５のいずれか１項記載の伝送線路のモデル化方法。
【請求項７】
請求項１乃至６のいずれか１項記載のモデル化方法を動作可能にプログラムとして組み込んだことを特徴とするモデル化装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、伝送線路の測定結果から、高精度の伝送シミュレーションを実現するための伝送線路モデルを作成する方法および装置に関し、複数の特性間の重み付けに伴う試行錯誤を不要とし、複数の特性のフィッティングを自動化する方法および装置に関する。
【背景技術】
【０００２】
伝送線路の特性をモデル化する方法として、特許文献１にあるように、伝送線路の構造、材質等の設計情報を基に計算で導出する方法が提案されている。また、伝送線路の特性を実測値とのパラメータフィッティングによってモデル化する手法が、特許文献２において提案されている。
【特許文献１】特開２００２−２５９４８３号公報
【特許文献２】特開２００４−３２８０６５号公報
【非特許文献１】飯島洋祐, 河西  勇二, 高橋栄一, 樋口哲也,「表皮効果および誘電損を考慮した伝送線路モデルの提案と検証」、電子情報通信学会論文誌Ｃ編、Vol.J89-C、No.5、pp.292-300、2006.
【発明の開示】
【発明が解決しようとする課題】
【０００３】
信頼性の高い、高速伝送回路の設計、開発を行うためには、伝送路特性の実測値に基づいた高精度のシミュレーションを行えることが必要不可欠である。しかし、従来の伝送線路の構造、材質等の設計情報を基に計算で導出する方法では、伝送線路の製造上の誤差や材質のバラツキ等の影響を考慮できない。
一方、従来の実測に基づいた伝送線路特性のモデル化手法では、正確な特性を得ることが難しい。なぜならば、従来の手法では、周波数特性もしくは時間軸応答特性の一方のみの測定結果から導出しているため、測定に固有の誤差、ノイズ等の影響がそのまま伝送線路特性に反映されてしまう。
【０００４】
本発明の目的は、伝送線路のモデル化に際し、伝送線路の特性を測定する際に含まれる、測定器が有する固有の誤差、ノイズ等の影響を抑制し、正確な伝送線路の特性にフィッティングさせるために、従来のような試行錯誤による解析の長時間化を解消した伝送線路特性のモデル化方法およびそのモデル化装置を提供することにある。
【課題を解決するための手段】
【０００５】
本発明の解決手段は、概略、実測した周波数特性ならびに時間軸応答特性の二つを同時に用いたモデル化方法にある。二つの特性をモデルにフィッティングするために、多目的最適化方法を用いることを特徴とする。これにより、フィッティングを自動化することが可能となる。また、測定系の誤差に影響をうけにくくなるので、信頼性の高い高精度のシミュレーションを行うことができる。
【０００６】
具体的には、以下のようになる。
（１）伝送線路のモデル化方法は、伝送線路の複数の伝送特性をそれぞれ測定し、伝送線路モデルの計算結果と、前記伝送線路の複数の測定結果とを同時に比較しながらフィッティングを行う伝送線路のモデル化方法において、多目的最適化手法を用いて複数のパレート最適解を求めることを特徴とする。
（２）上記（１）記載の伝送線路のモデル化方法は、前記複数の伝送特性として、周波数特性と時間軸応答特性を測定し、周波数特性のフィッティングと、時間軸応答特性のフィッティングを同時に行うことを特徴とする。
（３）上記（２）記載の伝送線路のモデル化方法は、前記周波数特性のフィッティングにおいて、伝送線路の特性を実測した特性と伝送線路モデルから計算した伝送特性の結果とを比較し、その振幅特性と位相特性それぞれの誤差が最小となるように伝送線路モデルを調整することを特徴とする。
（４）上記（２）記載の伝送線路のモデル化方法は、前記時間軸応答特性のフィッティングにおいて、インパルス波形を伝送線路に通して実測した受信波形と、伝送線路モデルから算出したインパルス応答波形の誤差を求め、その誤差が最小となるように伝送線路モデルを調整することを特徴とする。
（５）上記（１）記載の伝送線路のモデル化方法は、前記伝送線路モデルとして、表皮効果の周波数特性を考慮した伝搬定数γを用いることを特徴とする。
（６）上記（１）乃至（５）のいずれか１項記載の伝送線路のモデル化方法は、前記多目的最適化手法として、遺伝的アルゴリズムを用いることを特徴とする。
（７）モデル化装置は、上記（１）乃至（６）のいずれか１項記載のモデル化方法を動作可能にプログラムとして組み込んだことを特徴とする。
【０００７】
上記解決手段の基礎となる事項は、以下のように説明できる。
１．基本的な考え
近年、半導体素子の動作周波数が急速に向上する中、信号伝送に対する高速化の要求は日々高まっている。その理由は、半導体素子の高周波数化に比べて、伝送速度の高速化が遅れており、この遅れが伝送システム全体の高速化のボトルネックになっているからである。高速化の遅れの要因には次の２つがある。
・周波数が高くなるに従って伝送線路の振る舞いが複雑化することで、設計、開発が困難になること。
・マイクロ波以上の周波数領域で表皮効果、誘電損の影響によって、伝送線路の扱い方がより一層難しくなること、
これら２つの要因から、現在の高速伝送回路の設計、開発においてシミュレーションの重要性は増し、高精度の伝送線路シミュレーション性能が要求されている。
【０００８】
しかし、理論的な伝送線路シミュレーションでは、実際の複雑な伝送線路の振る舞いが表現できず、もはや精度面で限界に達している。そのため、実測に基づいたシミュレーションは必要不可欠となっている。現在、伝送線路シミュレーションにおいては、伝送線路の等価回路モデルが広く用いられているが、等価回路モデルを用いて実測結果に基づいた信頼性の高い、高精度のシミュレーションを実現するためには、図１に示す次の３つの手順が必要である：（１）伝送線路の測定（２）伝送線路特性の決定（３）等価回路モデルの作成、筆者らは文献［１］で表皮効果、誘電損を考慮した等価回路モデルを提案し、（２）から（３）の導出、つまり伝送線路特性から等価回路モデルを作成することを可能にした。そこで、本発明では（１）から（２）の導出、すなわち実測値に基づいて伝送線路特性を決定する手法を提案する。これにより、（１）から（３）の手順が統合されれば、ＳＰＩＣＥ等の回路シミュレータで実測に基づいた高精度なシミュレーションが可能になる。高速伝送回路の設計、開発のためのシミュレーションにおいては、伝送線路特性が高精度で表現され、かつ周波数特性、時間軸応答特性が解析可能であることが重要である。しかし、伝送線路特性の従来の決定手法では、これらを同時に満たすことができない。
【０００９】
そこで、本発明では実測値と伝送線路特性のシミュレーション結果とを比較しながら、伝送線路特性を実測値に近づけるフィッティング手法を用いて伝送線路特性の決定を行う。本発明では次の２つの特徴を持つフィッティングを用いて伝送線路特性を決定する手法を提案する。（１）周波数特性、時間軸応答特性の実測した２つの特性を同時に用いたフィッティング（２）多目的最適化手法を用いたフィッティング。このような特長をもつ提案手法の有効性を検証するために、本発明では実際に作製したマイクロストリップ線路を用いた実験を行う。
【００１０】
２．伝送線路特性の決定における課題
伝送線路特性として、一般的に伝搬定数γが用いられる。発明者は非特許文献［１］で表皮効果を考慮した伝播定数γを示し、γから等価回路モデルへの展開を可能にした。そこで、ここでは、２．１で非特許文献［１］で示した伝播定数γと、γから等価回路モデルへの展開方法について説明する。次に２．２で、この展開を実現する上で問題となる、伝送線路特性の決定における課題について述べる。
【００１１】
２．１  伝搬定数γから等価回路モデルの作成
一般的に伝送線路特性は伝搬定数γで表される。発明者らは非特許文献［１］で、表皮効果の影響を考慮した伝播定数γを次式数１のように表現した。
【数１】
2008116415.tif 000003
上記数１のαは振幅定数、βは位相定数、ωは角周波数を示す。LとCは、それぞれ伝送線路の単位長さ当たりのインダクタンス成分とキャパシタンス成分を示す。ここで、ωｏにおける表皮効果の大きさをＲｓｅとして、表皮効果の影響は下記数２で表される。
【数２】
2008116415.tif 000004
コンダクタンス成分Ｇは誘電損失をｔａｎδとして、下記数３で表される。
【数３】
2008116415.tif 000005
で与えられる。非特許文献［１］において数１で表現する伝搬定数γを表現することが可能な等価回路モデルを提案し、その回路パラメータの決定方法を示した。提案した等価回路モデルを図２に示す。図２の単位長さ当たりの各回路パラメータは次のように決定される。
【００１２】
【数４】
2008116415.tif 000006
ここで、Ｇｘはある角周波数ωｘでのωＣｔａｎδの値である。
以上から、等価回路モデルへの展開を行うためには、伝搬定数γが決定していることが必要であり、等価回路モデルによって高精度のシミュレーションを行うためには伝搬定数γの決定方法が重要になる。
【００１３】
２．２伝送線路特柱の決定における課題
高速伝送回路の設計、開発に使用するための伝送線路特性の決定には次の２つの課題がある。
（１）実測値に基づいた高精度の伝送線路特性が決定できること
この理由は、信頼性の高い高速伝送回路の設計、開発を行うためには、実測値に基づいた高精度なシミュレーションが行えることが必要不可欠だからである。伝送線路の構造、材質等の設計情報を基に導出する従来の伝送線路特性では、伝送線路の製造上の誤差や材質のバラつき等の影響を考慮できないため、正確な伝送線路の特性を表現しにくい。
（２）周波数特性、時間軸応答特性の両方の特性が表現可能な伝送線路特性が決定できること
【００１４】
この理由は、フィルタ等の主に周波数領域の特性を扱う回路の設計、開発には周波数特性が必要であり、イコライザ等の時間領域の特性を扱う回路の設計、開発には時間軸応答特性が必要だからである。しかし、実際の測定した結果には、測定に固有の誤差、ノイズ等が含まれるため、従来の実測に基づいた伝送線路特性の決定手法では正確な伝送線路特性を得ることが難しい。なぜなら、従来の手法では周波数特性か時間軸応答特性の一方のみの測定結果から伝送線路特性を直接導出するため、測定に固有の誤差、ノイズ等の影響がそのまま伝送線路特性に反映されるからである。このため、従来手法を用いて決定した伝送線路特性で高精度な周波数特性と時間軸応答特性の両方を表現することは困難である。
【００１５】
３．伝送線踏特性の決定手法の提案
ここでは、３．１で本発明で提案する伝送線路特性の決定手法の２つの特徴について説明し、その後の３．２、３．３でそれぞれの特徴について説明する。３．１提案手法の特徴２．２で述べた課題を解決する手段として、本発明では以下の２つの特徴を持つ、フィッティングを用いた伝送線路特性の決定手法を提案する。ここで、フィッティングとは実測値と伝送線路特性のシミュレーション結果とを比較し、この誤差が最小になるように伝送線路特性を調整し、決定する作業である。フィッティングを用いることで、従来手法のように伝送線路特性にノイズ等の影響が直接反映されてしまうことがなくなり、高精度な伝送線路特性を決定することが可能になる。
（１）周波数特性、時間軸応答特性の実測した２つの特性を用いたフィッティング
周波数特性と時間軸応答特性の実測した２つの測定結果を同時に用いたフィッティングにより、伝送線路特性を決定する。この2つの測定では、別々の測定系を使用しているため、測定に固有の誤差、ノイズ等が異なる。このため、2つの測定結果の間で重み付けを調整してフィッティングを行うことで、フィッティングに対するそれらの影響を低減できる。これによって、従来手法のように誤差、ノイズ等の影響が直接フィッティング結果に反映されることを回避でき、周波数特性、時間軸応答特性の両方で高精度なシミュレーションが実現できる。詳細は３．２で述べる。
（２）多目的最適化手法を用いたフィッティング
（１）のフィッティングを行うことで、評価関数が周波数特性と時間軸応答特性の複数になる。しかしながら、これによって、人手による複数の特性間の重み付け作業が必要になる。このため、筆者らは多目的最適化手法を用いたフィッティング手法を提案する。これによって、人手による重み付け作業が不要となり、フィッティングが自動化される。なお、多目的最適化の結果には複数の解候補が得られ、その中からユーザが誤差、ノイズ等の影響を低減する解を選択できる。詳細は３．３で述べる。
【００１６】
３．２周波数特性、時聞軸応答特性の２つを同時に用いたフィッティング
図３に提案する周波数特性（Ｆｒｅｑｕｅｎｃｙ  ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ））と時間軸応答特性（Ｔｉｍｅ−ｂａｓｅ  ｒｅｓｐｏｎｓｅ  ｗａｖｅｆｏｒｍ）の２つを同時に用いたフィッティング手法の概略を示す。図３のように伝送線路の周波数特性と時間軸応答特性を別々に測定し、伝送線路特性と、これら２つの測定結果とを同時に比較しながらフィッティングを行う。本発明では、伝送線路特性として数１の式で与えられる伝搬定数γを用いる。フィッティングにおける周波数特性と時間軸応答特性それぞれの評価関数を以下に示す。周波数特性のフィッテッィングでは、図３のようにネットワークアナライザを用いて実測したＳ２１特性と伝搬定数γとを比較し、数５の式で求められる△ｆｒｅｑを評価関数とする。具体的には、図４のように周波数特性の振幅特性、位相特性のシミュレーション値(Simulation result)を実測値(measurement result)とフィッティングする。ここで、S21特性の測定の前に、測定ケーブルを含めた状態でネットワークアナライザの校正を行い、機器と測定対象の線路を接続するケーブルの影響を除去する。α［ｎｅｐｅｒ／ｍ］は振幅特性、β［ｒａｄ／ｍ］は位相特性を表す。周波数特性を比較する際に、測定したＳ２１と伝搬定数γでは単位が異なるため、単位長さ当たりに換算したＳ２１特性を次のように変換して伝搬定数γとの誤差△ｆｒｅｑを下記数５の式で計算する。
【００１７】
【数５】
2008116415.tif 000007
ここで、ΔAmp、ΔPhaseはそれぞれ振幅特性、位相特性の誤差を示す。α、βはそれぞれ伝搬定数γの振幅定数、位相定数である。
時間軸応答特性のフィッティングでは、図３のように信号発生器で発生させたインパルス波形を伝送線路に通してオシロスコープで実測した受信波形と、伝搬定数γから算出したインパルス応答を比較し、数８の式で求められる２乗誤差和△ｗａｖｅを評価関数とする。具体的には、図５のように、Ｔｒａｎｓｍｉｔｔｅｄ  ｗａｖｅ（送信波）とＲｅｃｅｉｖｅｄ  ｗａｖｅ（受信波）に基づき、受信波(Simulation result)の時間軸応答特性の振幅特性、遅延特性を実測値(measurement result)とフィッティングする。インパルス応答は伝搬定数γから逆フーリエ変換（ＩＦＦＴ）した結果と、実測した入力波形との畳み込み演算を行うことで得られる。このように、時間軸応答特性のフィッティングでは、実測した入力波形から算出したインパルス応答と、実測した受信波形を比較する。これによって、伝送線路のみの特性を比較することができ、機器と測定対象の線路を接続するケーブルの影響を取り除くことができる。インパルス応答ＷＩＦＦＴ（＝ｙ0，ｙ1，・・・，ｙｋ）は、サンプリング時間をｔＤ、サンプリング点数をＮ、入力波形をｚ０，ｘ１，・・・，ｘｍとしたとき、次のように算出できる。
【００１８】
【数６】
2008116415.tif 000008
【数７】
2008116415.tif 000009
ここで、Hkは周波数fkでの伝送線路特性、hkはそのインパルス応答を示す。jは虚数単位である。この算出した結果ＷＩＦＦＴと、実測値Ｗｍｅａｓ（=a0、a1、・・・ak）から２乗誤差和△ｗａｖｅを次の数８の式のように計算する。
【数８】
2008116415.tif 000010
【００１９】
以上のように求められた周波数特性と時間軸応答特性の両方の評価関数△ｆｒｅｑと△ｗａｖｅを最小化することで、周波数領域、時間領域で実測値に基づいた高精度な伝送線路特性が得られる。さらに、この2つの特性間で重み付けを行うことで、フィッティングに対する誤差、ノイズ等の影響を低減することができる。具体的には、フィッティング結果に対して、誤差、ノイズ等の影響が最も小さくなるような特性間の重み付け作業を行う。しかし、この重み付け作業を人の手で行うためには、試行錯誤が必要であり、手間と時間が必要になる。このため、本提案では、フィッティングに多目的最適化手法を適用し、この重み付け作業を自動化した。
【００２０】
３．３  多目的最適化手法を用いたフィッティング
ここでは、３．３．１で多目的最適化手法の概念とその結果得られるパレード最適解について説明したのち、３．３．２では３．２で定義した２つの評価関数に対して多目的最適化を適用するフィッティング手法について提案する。
【００２１】
３．３．１  多目的最適化問題およびパレート最適解
多目的最適化問題とは、複数個の評価関数を同時に最適化する問題である。直感的にも分かるように、一般にはすべての評価関数を同時に最小にすることはできず、ほとんどの場合に単独解が存在しない。しかし、評価関数間に競合的な関係があることが、多目的最適化の本質である。このため、この多目的最適化問題の合理的な解として、パレート最適解と呼ばれるものを考える。パレート解(Pareto-optional Solutions)とは、いずれかの評価関数については必ず他の解よりも優れているような解の集合である。
パレート最適解の例として、２目的最適化問題、
【数９】
2008116415.tif 000011
を考える。図６は横軸にｆ１、縦軸にｆ２の値をプロットしたものである。この問題のパレート最適解に対応する（ｆ１，ｆ２）を、図６中の太線部分として示す。ｘ＝０とｘ＝２の間にある解すべては、２つの評価関数間でトレードオフが成立しており、パレート最適解である。たとえば、ｘ＝０という解は、ｆ１に関しては最適であるが、ｆ２に関しては最適ではない。このようなパレート最適解の集合を求めることを、多目的最適化手法と呼ぶ。
【００２２】
３．３．２  多目的最適化を用いたフィッティング手法の提案
３．２で述べた周波数特性と時間軸応答特性の２つの評価関数を用いたフィッティングを行うために、本発明では多目的最適化手法を用いた伝送線路特性のフィッティング手法を提案する。多目的最適化を用いることで、複数の評価関数間での重み付け作業を必要としないフィッティングが可能になり、さらに最終的に得られるパレート最適解の中から、シミュレーションの目的に合わせた最終的な解をユーザが選択できる利点も得られる。
本提案のフィッティングでは、多目的最適化の評価関数として次の２つを用いる。
（１）実測値との周波数特性の誤差△ｆｒｅｑ
（２）実測値との時間軸応答特性の誤差△ｗａｖｅ
多目的最適化を用いたフィッティングによって、これらの２つの誤差を最小化するように、以下に示す伝搬定数γの回路パラメータを調整し、パレート最適解を求める。
本発明で行う伝送線路特性のフィッティングでは、数１の式の伝搬定数γの回路パラメータＲｓｅ、ωｏ、Ｌ、Ｃ、ｔａｎδを調整する。そこで、これらのパラメータを多目的最適化の６つの調整パラメータｐａｒａｍ［０］〜［５］で以下のように表現する。
【００２３】
【数１０】
2008116415.tif 000012
ここで、Ｌ０、Ｃ０は、実測した特性インピーダンスＺ０と位相特性βから次のように求める。
【数１１】
2008116415.tif 000013
【数１２】
2008116415.tif 000014
【００２４】
以上の多目的最適化を用いたフィッティング手法によって、実測した周波数特性、時間軸応答特性から伝搬定数γを決定する。多目的最適化では結果として複数の解候補、つまりパレート最適解が得られるので、この中から最終的に１つの解をユーザがシミュレーションの目的に応じて、選択すればよい。
【００２５】
（モデル化装置）
上記伝送線路特性を決定するモデル化方法の手順をプログラムとして格納および実行するために、装置として記憶装置とＩ／Ｏ装置を備えたマイクロプロセッサを有する演算処理装置、例えば、パーソナルコンピュータ等を用いる。この演算処理装置は、持ち運び可能であるため、実測現場での適切なモデル化処理が可能となる。
【発明の効果】
【００２６】
本発明では、（１）多目的最適化手法を用いたフィッティングによる伝送線路特性の決定手法の提案、（２）実際のマイクロストリップ線路を用いた提案手法の検証実験の２つを行った。
前者の提案手法では、実測した周波数特性と時間軸応答特性の２つの結果を同時に用いて多目的最適化によるフィッティングを行うことで、それぞれの特性間の重み付け作業を必要とせずに、フィッティングを自動化できた。さらに、得られたパレート解の中からシミュレーションの目的に合わせて最適な伝送線路特性を決定することができた。
後者の実験結果から、最適化の結果として得られたパレート最適解から周波数特性、時間軸応答特性で高精度のシミュレーションが可能なパラメータを決定し、周波数軸、時間軸上で高精度かつ、実測値を用いた高精度のシミュレーションが可能になることを示すことができた。
本発明のモデル化装置は、上記モデル化方法を好ましい形で動作可能にプログラムとして組み込むことができたので、モデル化方法の効果をそのまま達成することができる。
【発明を実施するための最良の形態】
【００２７】
本発明の実施の形態を図に基づいて詳細に説明する。
【実施例】
【００２８】
４．マイクロストリップ線路を用いた提案手法の検証
ここでは、３．で提案した伝送線路特性の決定手法の有効性を検証するために、実際に作製したマイクロストリップ線路を用いて実験を行なった。４．１では実験の条件を述べ、４．２で実験の結果を述べ、その中で得られたパレート最適化の特性について評価した。
【００２９】
４．１  実験条件
図７に作製したマイクロストリップ線路の寸法を示す。作製したマイクロストリップ線路は、特性インピーダンスを５０［Ω］、長さ１００［ｃｍ］である。線路の誘電体にはＦＲ−４（ｔａｎδ＝０．０２）を用い、導体(Conductor)はＣｕを使用した。
まず、作製したマイクロストリップ線路の周波数特性と時間軸応答特性を測定する。周波数特性は、ネットワークアナライザを用いてＳ２１特性を測定した。時間軸応答特性は、時間幅０．２４［ｎｓ］のインパルス波形を入力し、マイクロストリップ線路１００［ｃｍ］の伝送波形をサンプリングした。この結果を用いて、３．で述べた多目的最適化を用いたフィッティング手法によって伝搬定数γを決定した。本発明では、多目的最適化手法として遺伝的アルゴリズム（Ｇｅｎｅｔｉｃ  Ａｌｇｏｒｉｔｈｍ；ＧＡ）を用いた代表的な多目的最適化手法であるＮＳＧＡ−２を用いた。
【００３０】
４．２  実験結果
図８にフィッティングの結果得られたパレート最適解の集合を示す。図８の横軸は周波数特性の実測値とシミュレーション結果との誤差△ｆｒｅｑ、縦軸は時間軸応答特性についての誤差△ｗａｖｅを示す。
ここで、図8に示すようにパレート最適解を次の3つの領域（A、B、C）にわける。Aは時間軸応答特性の誤差が小さい領域、Bは周波数特性の誤差が小さい領域、CはAとBの中間にあたる領域である。そして、これら3つの領域の中から、代表的な解（パラメータ）a、b、cを１つずつ選択し、それぞれのパラメータについてシミュレーション特性を評価した。このシミュレーションした結果を図9、10に示す。
図９は実測値と、図８のａ、ｂ、ｃの３点のパラメータを用いてシミュレーションした周波数特性を示す。図９からわかるように、ｂ点のパラメータを用いたシミュレーション結果(simulation)が実測値(measured)と比較して最も誤差が小さく、ｒｍｓ誤差は０．１２０［ｎｅｐｅｒ／ｍ］である。これに対して、ａ点のパラメータを用いたシミュレーション結果は実測値との誤差が最大であり、ｒｍｓ誤差は０．３０９［ｎｅｐｅｒ／ｍ］である。ｃ点のパラメータを用いた結果は、実測値とのｒｍｓ誤差が０．１４４［ｎｅｐｅｒ／ｍ］である。このｃ点のシミュレーションでの誤差はノイズレベルと比較して十分に小さく、周波数軸上でシミュレーションする上で問題とならない範囲である。
【００３１】
図１０は実測値と、図８のａ、ｂ、ｃの３点のパラメータを用いてシミュレーションした時間軸応答特性を示す。Ｔｒａｎｓｍｉｔｔｅｄ  ｗａｖｅ（送信波）とＲｅｃｅｉｖｅｄ  ｗａｖｅ（受信波）のｔｉｍｅｓ（時間）対Ａｍｐ（振幅）特性を示す。図１０からわかるように、ａ点のパラメータを用いたシミュレーション結果が実測値と比較して最も誤差が小さく、ｍｓ誤差は０．２０１［ｍＶ］であり、実測値にほぼ一致している。これに対して、ｂ点のパラメータを用いたシミュレーション結果は実測値と比較し、遅廷特性のずれが０．３［ｎｓ］であり、ｒｍｓ誤差は０．８１１［ｍＶ］である。ｃ点のパラメータを用いた結果は、実測値と比較してｒｍｓ誤差は０．３７１［ｍＶ］である。図１０のシミュレーション結果を見てわかるように、この０．３７１［ｍＶ］のｒｍｓ誤差は時間軸上のシミュレーションを行う上で十分に小さい。
【００３２】
５．まとめ
４で行ったシミュレーション結果から、次の3つのようなシミュレーションの目的に応じて、図8のA、B、Cの領域を使いわける。これによって、時間軸応答特性と周波数特性の高精度なシミュレーションが実現できる。
・時間領域の特性を扱うシミュレーション
主に時間領域の伝送線路特性を扱うイコライザ等のシミュレーションの場合、時間軸応答特性の誤差が最も小さい図8に示すAの領域のパラメータを使用する。これによって、時間軸上で高精度なシミュレーションが行える。
・周波数領域の特性を扱うシミュレーション
主に周波数領域の伝送線路特性を扱う増幅器等のシミュレーションの場合、周波数特性の誤差が最も小さい図8に示すBの領域のパラメータを使用する。これによって、高精度な周波数特性のシミュレーションが行える。
・時間領域、周波数領域の両方の特性を扱うシミュレーション
この場合には、図8に示すCの領域のパラメータを使用する。Cのパラメータは、時間軸応答特性、周波数特性ともに誤差が十分に小さいため、幅広い高速伝送シミュレーションで使用することが可能である。
【００３３】
加えて、図8に示すa、b、cのような最終的なパラメータを選択する際に、測定条件等を考慮して最終的な解をユーザが選択することで、より信頼性の高い、高精度なシミュレーションが可能になる。例えば、時間軸応答特性の測定結果の方が信頼性が高い場合は、よりＡに近い値を選択することで信頼性の高いパラメータが得られる。提案手法に対して、一般的な単目的最適化手法を用いたフィッティングを行った場合には次のような欠点がある。
・重み付けの値によって図８のうちの１点の解しか得られないため、シミュレーションの目的にあった解を最終的に得るためには、重み付けを何度か変更して試行錯誤する必要がある。
・適用するアプリケーションが変わった場合には、また試行錯誤を繰り返してフィッティングをやり直す必要がある。
さらに、本発明では、回路シミュレータでの計算結果を評価関数に追加することで、図１に示した（１）から（３）までの全ての手順が自動化でき、より高精度な伝送線路解析を実現できる。その際のフィッティングにも、多目的最適化を用いた本提案手法は、評価関数の数を選ばず、重み付け作業が不要である点で有効である。
【産業上の利用可能性】
【００３４】
本発明は、ＥＤＡ（Electronic Design Automation）ツールとしてソフトウェア化が可能である。具体的には、高速なデータ伝送を行うプリント基板の設計や、高速なデータ伝送のための増幅器、イコライザ等のＬＳＩ設計用に特化したＣＡＤ（computer-aided design）ソフトウェアとして、製品化が可能である。現在、素子間のデータ伝送速度がボトルネックとなり、半導体素子の高周波数化と比較して、伝送システム全体の高速化が遅れている。このボトルネックを解消するために本発明は有効であり、ＣＡＤや計測器への適用が可能である。
【図面の簡単な説明】
【００３５】
【図１】伝送線路シミュレーションフローである。
【図２】伝送線路の等価回路モデルである。
【図３】周波数特性と時間軸応答特性の２つを同時に用いたフィッティング手法の概略図である。
【図４】周波数特性の振幅特性、位相特性のシミュレーション値を実測値とフィッティングする図である。
【図５】受信波の時間軸応答特性の振幅特性、遅延特性を実測値とフィッティングさせる図である。
【図６】２目的最適化問題のパラメータｆ１を横軸に、ｆ２を縦軸にプロットした図である。
【図７】マイクロストリップ線路の構成図である。
【図８】フィッティングの結果得られたパレート最適解の集合を示す図である。
【図９】図８に示す３つ点（ａ、ｂ、ｃ）での解（パラメータ）を選択し、実測値と、図８のａ、ｂ、ｃの３点のパラメータを用いてシミュレーションした周波数特性を示す。
【図１０】実測値と、図８のａ、ｂ、ｃの３点のパラメータを用いてシミュレーションした時間軸応答特性を示す。
【符号の説明】
【００３６】
ｄｘＬｎ    インダクタンス成分
ｄｘＣｎ    容量成分
ｄｘＧｎ    コンダクタンス成分
【図１】
2008116415.tif 000015
【図２】
2008116415.tif 000016
【図３】
2008116415.tif 000017
【図４】
2008116415.tif 000018
【図５】
2008116415.tif 000019
【図６】
2008116415.tif 000020
【図７】
2008116415.tif 000021
【図８】
2008116415.tif 000022
【図９】
2008116415.tif 000023
【図１０】
2008116415.tif 000024