(57)【特許請求の範囲】
【請求項１】
複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により得られる第１の面積と、前記平面において前記複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して前記第１の要素値列から補間計算される第２の要素値列により得られる第２の面積との差が、許容範囲内であるか判断する判断処理と、
前記判断処理を１又は複数の間隔について実施させることで、前記複数の散乱パラメータの全要素について前記差が前記許容範囲内となる間隔を特定する特定処理と、
を、コンピュータに実行させるためのプログラム。
【請求項２】
前記複数の第１の周波数における間隔を抽出する処理
をさらに前記コンピュータに実行させ、
前記特定処理において、抽出された前記間隔のうち最も長い間隔から順番に前記判断処理を実施させる
請求項１記載のプログラム。
【請求項３】
前記複数の第１の周波数における間隔が均等である場合には、第１の周波数の数が２の冪乗個となるように前記第１の周波数の最大値を決定する処理
をさらに前記コンピュータに実行させる請求項２記載のプログラム。
【請求項４】
前記複数の第１の周波数の最大値又は設定された周波数の最大値まで、特定された前記間隔でサンプリングした場合のサンプリング数が２の冪乗個となっていない場合には、サンプリング数が２の冪乗個となるように、特定された前記間隔を短縮する処理
をさらに前記コンピュータに実行させるための請求項１乃至３のいずれか１つ記載のプログラム。
【請求項５】
前記第１の面積は、前記第１の要素値列のうち連続する２つの第１の要素値を直線で結んだ辺と、当該２つの第１の要素値から周波数の軸に下ろした垂線による辺と、当該２つの第１の要素値に対応する２つの第１の周波数の間隔に対応する辺とを有する台形を用いて算出され、
前記第２の面積は、前記第２の要素値列のうち連続する２つの第２の要素値を直線で結んだ辺と、当該２つの第２の要素値から周波数の軸に下ろした垂線による辺と、当該２つの第２の要素値に対応する２つの第２の周波数の間隔に対応する辺とを有する台形を用いて算出される
請求項１乃至４のいずれか１つ記載のプログラム。
【請求項６】
複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により得られる第１の面積と、前記平面において前記複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して前記第１の要素値列から補間計算される第２の要素値列により得られる第２の面積との差が、許容範囲内であるか判断する判断処理と、
前記判断処理を１又は複数の間隔について実施させることで、前記複数の散乱パラメータの全要素について前記差が前記許容範囲内となる間隔を特定する特定処理と、
を含み、コンピュータにより実行される情報処理方法。
【請求項７】
複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により得られる第１の面積と、前記平面において前記複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して前記第１の要素値列から補間計算される第２の要素値列により得られる第２の面積との差が、許容範囲内であるか判定する判定部と、
前記判断処理を１又は複数の間隔について実施させることで、前記複数の散乱パラメータの全要素について前記差が前記許容範囲内となる間隔を特定する特定部と、
を有する情報処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は、回路解析システムのための制御パラメータの設定技術に関する。
【背景技術】
【０００２】
Ｓパラメータとは、高周波電子回路や高周波電子部品の特性を表わすために用いられる回路網パラメータの一つであり、散乱パラメータ（Scattering parameter）とも呼ばれる。ｎ個のポート（Ｓパラメータで特性化される系では、一般的に端子のことをポートと呼ぶ）を有し、Ｓパラメータで特性化されている電子回路系（一個の高周波電子部品も１つの小さな系と考えられる）を考える。１つのポートに正弦波を入力したとき、その正弦波が他のポートへ透過する割合、その正弦波が正弦波を入力したポートへ反射する割合、及びそれら透過波及び反射波の入力正弦波に対する位相遅れを表わしたものが、その系のＳパラメータである。
【０００３】
Ｓパラメータは複素数であり、複素数の絶対値で透過割合（反射割合）を、複素数の偏角で位相遅れを表現している。透過割合と透過波の位相遅れをまとめて透過係数、反射割合と反射波の位相遅れをまとめて反射係数と呼ぶ。
【０００４】
系に複数のポートがあるとき、それらのポートの各々について正弦波を入力し、そのときの透過係数と反射係数を得る。即ち、１つの入力ポートに対して、反射係数が１個、透過係数が（ｎ−１）個の合計ｎ個のパラメータが得られる。これがｎ入力ポート分あるので、ｎ×ｎ個のパラメータがある。これらのＳパラメータが記述されたファイルが、Ｓパラメータファイルである。
【０００５】
入力正弦波の周波数が変わると特性が変わるため、Ｓパラメータファイルには、入力正弦波の周波数毎にｎ×ｎ個のパラメータが行列形式で記述されている。
【０００６】
プリント基板設計において、ＳＩ（Signal Integrity）検証のための伝送線路解析や、ＰＩ（Power Integrity）検証のための電源基板解析では回路シミュレータを用いる。解析対象モデルは、抵抗Ｒ、インダクタＬ、キャパシタＣなどの受動素子の他に、フィルタ、ヴィア、コネクタ、伝送線路系などをＳパラメータで特性を表現した素子（以下、Ｓ素子）も用いられる。
【０００７】
回路シミュレータは、これら受動素子、Ｓ素子、ダイオードやＭＯＳ（Complementary Metal Oxide Semiconductor）トランジスタ等の能動素子、各種電源が電気的に接続した回路網を認識する。回路網で電気的に等電位を示す「節点」には、素子が接続している。そして、回路シミュレータは、回路網理論の基本法則であるキルヒホッフの法則に従い、各節点の電圧や各素子に流入する電流の時間変化を計算する。
【０００８】
回路シミュレータは、Ｓ素子に流入する電流値を算出するために、Ｓパラメータを逆フーリエ変換してインパルス応答関数を求め、時間領域においてＳ素子の端子電圧とインパルス応答関数の畳み込積分を計算することが一般的である。更に、計算機を使った逆フーリエ変換のアルゴリズムとして高速逆フーリエ変換（Inverse Fast Fourier Transform）アルゴリズムを用いることが一般的である。ＩＦＦＴ関数に渡すデータをＳパラメータから等間隔周波数でＮ個抽出するが、ＩＦＦＴのアルゴリズムの条件から、Ｎは２の累乗ということになっている。
【０００９】
Ｓパラメータを得るには、スペクトルアナライザによる実測を行うか、Ｓパラメータ抽出シミュレータを使ってＳパラメータを算出するか、又は回路部品ベンダーから提供されているＳパラメータライブラリを入手するなどの方法が考えられる。ところで、これらの方法で入手したＳパラメータは周波数が等間隔で記述されているとは限らず、周波数の個数（データ項目数又はサンプリング数とも呼ぶ）も２の累乗であるとは限らない。
【００１０】
Ｓ素子を取り扱うことができる回路シミュレータには、通常、ＩＦＦＴのサンプリングに関連する様々な制御パラメータが用意されている。具体的には、Ｓパラメータのサンプリング周波数間隔（FBASEと呼ぶ）、Ｓパラメータをどの周波数までサンプリングするのかを示す最大周波数（FMAXと呼ぶ）、データ補間が発生する場合、補間方法（INTERPOLATION）、補外方法（LOWPASS,HIGHPASS）などである。
【００１１】
FBASE, FMAX, INTERPOLATION, LOWPASS, HIGHPASSの値を変更すると、ＩＦＦＴによって計算される結果（すなわちインパルス応答関数）は異なる。そうすると、インパルス応答関数と信号の畳み込積分結果も異なって、結果として回路シミュレーション結果が異なる。しかしながら、従来このような制御パラメータを適切に設定することはなかった。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２０００−３４６８９１号公報
【特許文献２】特開平１１−２９５３６５号公報
【特許文献３】特開平８−２６２０８０号公報
【特許文献４】特開平４−１６８９０７号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
本技術の目的は、一側面において、回路解析システムのための制御パラメータを適切に設定するための技術を提供することである。
【課題を解決するための手段】
【００１４】
本技術の情報処理方法は、（Ａ）複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により得られる第１の面積と、上記平面において複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して第１の要素値列から補間計算される第２の要素値列により得られる第２の面積との差が、許容範囲内であるか判断する判断処理と、（Ｂ）判断処理を１又は複数の間隔について実施させることで、複数の散乱パラメータの全要素について上記差が許容範囲内となる間隔を特定する特定処理とを含む。
【発明の効果】
【００１５】
回路解析システムのための制御パラメータを適切に設定できるようになる。
【図面の簡単な説明】
【００１６】
【図１】図１は、実施の形態に係るシステムの概要を示す図である。
【図２】図２は、実施の形態に係る処理フローを示す図である。
【図３】図３は、第１のテーブルの一例を示す図である。
【図４】図４は、第２のテーブルの一例を示す図である。
【図５】図５は、実施の形態に係る処理フローを示す図である。
【図６】図６は、台形面積を説明するための図である。
【図７】図７は、台形面積を説明するための図である。
【図８】図８は、周波数間隔を均一にして線形補間を行った後の台形面積を説明するための図である。
【図９】図９は、周波数間隔を均一にして線形補間を行った後の台形面積を説明するための図である。
【図１０】図１０は、本実施の形態に係る処理フローを示す図である。
【図１１】図１１は、効果を説明するための波形図である。
【図１２】図１２は、効果を説明するための波形図である。
【図１３】図１３は、コンピュータの機能ブロック図である。
【発明を実施するための形態】
【００１７】
本技術の実施の形態に係るシステムの構成例を図１に示す。本実施の形態に係るシステムは、情報処理装置１００と、当該情報処理装置１００とネットワークなどで接続される回路シミュレータを含む回路解析システム２００とを有する。なお、情報処理装置１００と回路解析システム２００とは一体の場合もある。
【００１８】
回路解析システム２００には、情報処理装置１００の出力である修正後Ｓパラメータファイル３００が入力されるが、回路解析システム２００については従来と変わらないので、これ以上説明しない。
【００１９】
情報処理装置１００は、入力部１１０と、第１データ格納部１２０と、制御パラメータ決定部１３０と、第２データ格納部１４０と、出力部１５０とを有する。制御パラメータ決定部１３０は、最大値決定部１３１と、間隔決定部１３２とを有する。間隔決定部１３２は、誤差判定部１３２１を有する。
【００２０】
本実施の形態では、入力部１１０は、ユーザが何らかの手段により入手したＳパラメータファイル及び各種パラメータを受け付け、第１データ格納部１２０に格納する。制御パラメータ決定部１３０の最大値決定部１３１は、第１データ格納部１２０に格納されているデータを用いて、制御パラメータである周波数最大値FMAXを決定する。また、間隔決定部１３２は、第１データ格納部１２０に格納されているデータを用いて、制御パラメータである周波数間隔FBASEを決定する。誤差判定部１３２１は、周波数間隔FBASEを決定する際に後に述べる誤差が許容範囲内であるか否かを判断する。制御パラメータ決定部１３０は、第１データ格納部１２０に格納されているＳパラメータファイルに、決定された制御パラメータを書き込んで、第２データ格納部１４０に格納する。出力部１５０は、第２データ格納部１４０に格納されている修正後Ｓパラメータファイル３００を表示装置や他のコンピュータ（回路解析システム２００を含む）に出力する。
【００２１】
次に、図２乃至図１０を用いて本情報処理装置１００の処理内容について説明する。まず、入力部１１０は、ユーザからのＳパラメータファイルの入力を受け付けるか、他のコンピュータに格納されているＳパラメータファイルを読み込んで、第１データ格納部１２０に格納すると共に、当該Ｓパラメータファイルから第１のテーブルを生成し、第１データ格納部１２０に格納する（図２：ステップＳ１）。例えば、図３に示すような第１のテーブルが生成される。図３は、２ポートの場合のＳパラメータの例を示しており、周波数（サンプリング点とも呼ぶ）毎に、ＳパラメータであるＳ１１、Ｓ１２、Ｓ２１及びＳ２２の値が登録されるようになっている。なお、Ｓ１１乃至Ｓ２２の各々は複素数であり、実部と虚部との値がペアになっている。Ｓパラメータファイルには、このようなデータが登録されているので、このようなデータを読み出してテーブルに登録する。
【００２２】
また、入力部１１０は、以下で算出される面積差の許容誤差（許容範囲とも呼ぶ）εの入力を受け付け、第１データ格納部１２０に格納する（ステップＳ３）。さらに、入力部１１０は、解析のターゲット周波数Ｆ_target及び高調波（ターゲット周波数のＭ倍）の倍数Ｍの入力を受け付け、第１データ格納部１２０に格納する（ステップＳ５）。なお、Ｍについては例えば５や１０といったデフォルト値をそのまま用いるようにしても良い。
【００２３】
そして、制御パラメータ決定部１３０の間隔決定部１３２は、第１データ格納部１２０に格納されている第１のテーブルから、周波数間隔を抽出して、第１データ格納部１２０における第２のテーブルに登録する（ステップＳ７）。第２のテーブルに格納されるデータの一例を図４に示す。図４の例では、抽出された間隔の種類が列挙される。なお、複数の種類の間隔が抽出される場合もあれば、１つの間隔のみが抽出される場合もある。図３の例では、複数の種類の間隔が抽出される。
【００２４】
また、間隔決定部１３２は、周波数間隔が等間隔か否かを、第２のテーブルから特定し、等間隔でない場合には周波数間隔の最大値及び最小値を特定する（ステップＳ９）。例えば、周波数間隔が等間隔である場合には、変数Equal_delta_fに周波数間隔を代入し、等間隔ではない場合には、変数Equal_delta_fに「−１」を代入する。また、周波数間隔の最大値を変数f_step_maxに代入し、周波数間隔の最小値を変数f_step_minに代入する。
【００２５】
そして、制御パラメータ決定部１３０の最大値決定部１３１は、第１のテーブルにおける周波数の最大値が、ターゲット周波数Ｆ_target＊Ｍ以上であるか判断する（ステップＳ１１）。周波数の最大値がターゲット周波数Ｆ_target＊Ｍ以上である場合には、最大値決定部１３１は、仮に、周波数最大値FMAX＝Ｆ_target＊Ｍと設定する（ステップＳ１３）。そして処理は端子Ａを介して図５の処理に移行する。一方、周波数の最大値がターゲット周波数Ｆ_target＊Ｍ未満である場合には、Ｆ_target＊Ｍまではデータがないので、最大値決定部１３１は、周波数最大値FMAXに第１のテーブルにおける周波数最大値を設定する（ステップＳ１５）。そして処理は端子Ａを介して図５の処理に移行する。
【００２６】
図５の処理の説明に移行して、間隔決定部１３２は、周波数間隔が均一であるか判断する（ステップＳ１７）。すなわち、変数Equal_delta_fに「−１」ではない値が設定されているか判断する。均一であれば、間隔決定部１３２は、この周波数間隔を周波数間隔FBASEに設定する（ステップＳ１９）。このように周波数が等間隔であれば問題は無いので、最大値決定部１３１の処理に移行する。すなわち、最大値決定部１３１は、第１のテーブルにおいてサンプリング数（第１のテーブルにおけるレコード数）が２の冪乗個であるか判断し、２の冪乗個でない場合には周波数最大値FMAX以下で最大の２の冪乗番目の周波数を特定して周波数最大値FMAXに設定する（ステップＳ２１）。そして処理は端子Ｃを介して図１０の処理に移行する。
【００２７】
一方、周波数間隔が均一ではない場合には、間隔決定部１３２の誤差判定部１３２１は、第１のテーブルにおけるＳパラメータの各要素（Ｓ１１乃至Ｓ２１の各実部要素Ｓmn）について、当該要素（無単位）と周波数とで張られる平面において各周波数について当該周波数とその要素値（パラメータ値とも呼ぶ）と次の周波数とその要素値とで形成される台形（場合によっては三角形）の面積を算出すると共にその総和Area_Smnを算出する（ステップＳ２３）。
【００２８】
Ｓパラメータファイルにおける最低周波数FMINからFMAXまでの各サンプリング点についてＳパラメータのＳ１１の実部の値を、周波数とＳ１１の実部とで張られる平面においてプロットすると、例えば図６のようになる。このような場合、最低周波数FMINから順番に各周波数について、処理に係る周波数Ｆiと、当該処理に係る周波数ＦiにおけるＳ１１の実部の値Ｓmn（Ｆi）と、当該処理に係る周波数の次の周波数Ｆ(i+1)と、当該次の周波数Ｆ(i+1)におけるＳ１１の実部Ｓmn(Ｆ(i+1))とで囲まれる台形の面積（図６のハッチング部分）を算出する。そしてその総和Area_Smnを算出する。これが基準値となる。
【００２９】
また、FMAXからFMINまでの各サンプリング点についてＳパラメータＳ２１の実部の値を、周波数とＳ２１の実部とで張られる平面においてプロットすると、例えば図７のようになる。この場合にも同様に台形面積を各周波数について算出して、その総和Area_Smnを算出する。
【００３０】
また、間隔決定部１３２の誤差判定部１３２１は、周波数間隔FBASEに周波数間隔の最大値f_step_maxを設定する（ステップＳ２５）。
【００３１】
その後、間隔決定部１３２の誤差判定部１３２１は、第１のテーブルに格納されているデータから、FBASE間隔のＳパラメータ値を線形補間で算出し、例えば第１データ格納部１２０に格納する（ステップＳ２７）。そして、間隔決定部１３２の誤差判定部１３２１は、Ｓパラメータの各要素について、上記の平面においてFBASE間隔の各周波数について当該周波数とその補間要素値と次の周波数とその補間要素値とで形成される台形（場合によっては三角形）の面積を算出すると共にその総和Area_Smn_fbaseを算出する（ステップＳ２９）。
【００３２】
FMINからFMAXまでFBASE間隔で線形補間されたＳ１１の値（実部）を、周波数とＳ１１の実部とで張られる平面においてプロットすると、例えば図８に示すようになる。図８は一部分を拡大して示しているが、周波数間隔はFBASEで同一となっており、FBASE間隔で設定された新サンプリング点Ｆnと、当該新サンプリング点ＦnにおけるＳ１１の実部Ｓ’mn(Ｆn)と、新サンプリング点の次の新サンプリング点Ｆ(n+1)と、当該次の新サンプリング点Ｆ(n+1)におけるＳ１１の実部Ｓ’mn(Ｆ(n+1)）とで囲まれる台形の面積（図８のハッチング部分）を算出する。そしてその総和Area_Smn_fbaseを算出する。
【００３３】
また、FMAXからFMINまでFBASE間隔で線形補間されたＳ２１の値（実部）を、周波数とＳ２１の実部とで張られる平面においてプロットすると、例えば図９のようになる。この場合にも同様に台形面積を各新サンプリング点について算出して、その総和Area_Smn_fbaseを算出する。
【００３４】
そして、間隔決定部１３２の誤差判定部１３２１は、Ｓパラメータの各要素について｜Area_Smn−Area_Smn_fbase｜＜εが成り立つか判断する（ステップＳ３１）。全ての要素についてこの条件が満たされる場合には、すなわちこのような面積についての誤差が許容範囲内である場合、端子Ｄを介して図１０の処理に移行する。
【００３５】
一方、いずれかの要素についてこの条件を満たさない場合には、間隔決定部１３２の誤差判定部１３２１は、第２のテーブルにおいて次に小さな周波数間隔をFBASEに設定する（ステップＳ３３）。但し、次に小さな周波数間隔が存在しない場合には０を設定する。そして、間隔決定部１３２の誤差判定部１３２１は、FBASEが最小周波数間隔f_step_min以上であるか判断する（ステップＳ３５）。FBASEが最小周波数間隔f_step_min以上である場合には処理はステップＳ２７に戻る。一方、FBASEが最小周波数間隔f_step_min未満となっている場合には、端子Ｅを介して図１０の処理に移行する。
【００３６】
図１０の処理の説明に移行して、端子Ｅの後に、間隔決定部１３２は、FBASEに最小周波数間隔f_step_minを設定し、エラーメッセージと共にFBASEを、出力部１５０に表示装置などに対して出力させる（ステップＳ３７）。例えば、「これ以上サンプリング間隔を細かくしても精度は改善しません。Ｓパラメータのデータ間隔を細かくして測定し直してください」といったようなメッセージを出力する。そして処理はステップＳ４１に移行する。
【００３７】
一方、端子Ｄの後、すなわち面積ベースで誤差が許容範囲内である場合、間隔決定部１３２は、FBASEでのサンプリング数が２の冪乗個になるように、再度FBASEを調整する（ステップＳ３９）。例えば、FBASE間隔でFMAXまでサンプリングした場合のサンプリング数が２の冪乗個となっていればFBASEは調整しない。一方、FBASE間隔でFMAXまでサンプリングした場合のサンプリング数が２の冪乗個となっていない場合には、Ｎx≧FMAX／FBASEを満たす最小の整数Ｎxを計算し、２k≧Ｎxを満たす最小のｋを算出し、FBASE＝FMAX／２kを調整後のFBASEとして設定する。このようにすれば、FMAXを変更することなく、精度誤差が許容範囲内と判断された周波数間隔以下でサンプリング数が２の冪乗個という要件を満たす最も大きい周波数間隔を決定できるようになる。そして処理はステップＳ４１に移行する。
【００３８】
ステップＳ４１では、制御パラメータ決定部１３０は、これまでに確定したFBASE及びFMAXを第１データ格納部１２０に格納されているＳパラメータファイルのコメントとして例えばヘッダ部分に書き出して、第２データ格納部１４０に格納する（ステップＳ４１）。
【００３９】
その後、ユーザからの指示などに従って、出力部１５０は、修正後Ｓパラメータファイル３００を、第２データ格納部１４０から読み出して、回路解析システム２００に出力する（ステップＳ４３）。なお、本ステップは、エラーが発生していたり、ユーザが別の理由で指示しない場合などには行われないので、点線ブロックとして示されている。
【００４０】
以上のような処理を実施することで、Ｓパラメータファイルに規定される周波数における周波数間隔のうち最も大きい周波数間隔で精度低下が許容範囲内であるものを基準に周波数間隔を決定できるようになる。すなわち、できるだけ大きい周波数間隔を設定できるようになるので、回路解析システム２００における処理負荷の増加を抑えることができるようになる。
【００４１】
なお、FBASE及びFMAXといった制御パラメータを従来のようにデフォルト値をそのまま用いるような場合には、シミュレーション結果としては例えば図１１に示すようなオーバーシュートやアンダーシュートも大きな波形の信号が生成されてしまう。一方、上で述べたような処理によってFBASE及びFMAXを設定すれば、図１２に示すように、同じ信号でもオーバーシュートもアンダーシュートも小さくなった波形が生成される。すなわち、より適切な回路シミュレーションが実施されるようになる。
【００４２】
以上本技術の実施の形態について説明したが、本技術はこれに限定されるものではない。例えば上で述べた例では、当初Ｓパラメータファイルに規定されている周波数の周波数間隔に基づいて図５のステップＳ２７乃至Ｓ３５のループが回るようになっている。しかしながら、周波数間隔について固定値ずつ減少させてこのループを回すようにしても良い。また、サンプリング数を２の冪乗個に設定する処理については、様々なバリエーションが可能である。ステップＳ２１で述べた処理をステップＳ３９でも実施するようにしても良い。さらに、多少回路解析システム２００における処理負荷が上がることを許容できれば、ステップＳ３９でもｋをさらに小さな値に設定する場合もある。ステップＳ２１においても、同様に最大周波数をさらに小さな値を採用するようにしても良い。
【００４３】
また、図１に示したシステム構成図は一例であり、必ずしも実際のプログラムモジュール構成とは一致しない場合もある。また、処理フローについても処理結果が変わらない限り、処理を並列に実施したり、処理順番を入れ替えたりしても良い。
【００４４】
さらに、線形補間については、例えばＦｎの両側に規定されているＦi及びＦ(i+1)についてのＳmn(Fi)及びＳmn(F(i+1))から線形補間するだけではなく、より多くのサンプリング点の要素値から補間するようにしても良い。
【００４５】
また、上で述べた例ではＳパラメータファイルにFBASE及びFMAXを書き込む例を示したが、別途回路解析システム２００に入力する場合もある。
【００４６】
なお、上で述べた情報処理装置１００は、コンピュータ装置であって、図１３に示すように、メモリ２５０１とＣＰＵ（Central Processing Unit）２５０３とハードディスク・ドライブ（ＨＤＤ：Hard Disk Drive）２５０５と表示装置２５０９に接続される表示制御部２５０７とリムーバブル・ディスク２５１１用のドライブ装置２５１３と入力装置２５１５とネットワークに接続するための通信制御部２５１７とがバス２５１９で接続されている。オペレーティング・システム（ＯＳ：Operating System）及び本実施例における処理を実施するためのアプリケーション・プログラムは、ＨＤＤ２５０５に格納されており、ＣＰＵ２５０３により実行される際にはＨＤＤ２５０５からメモリ２５０１に読み出される。ＣＰＵ２５０３は、アプリケーション・プログラムの処理内容に応じて表示制御部２５０７、通信制御部２５１７、ドライブ装置２５１３を制御して、所定の動作を行わせる。また、処理途中のデータについては、主としてメモリ２５０１に格納されるが、ＨＤＤ２５０５に格納されるようにしてもよい。本技術の実施例では、上で述べた処理を実施するためのアプリケーション・プログラムはコンピュータ読み取り可能なリムーバブル・ディスク２５１１に格納されて頒布され、ドライブ装置２５１３からＨＤＤ２５０５にインストールされる。インターネットなどのネットワーク及び通信制御部２５１７を経由して、ＨＤＤ２５０５にインストールされる場合もある。このようなコンピュータ装置は、上で述べたＣＰＵ２５０３、メモリ２５０１などのハードウエアとＯＳ及びアプリケーション・プログラムなどのプログラムとが有機的に協働することにより、上で述べたような各種機能を実現する。
【００４７】
以上述べた本実施の形態をまとめると、以下のようになる。
【００４８】
本実施の形態に係る情報処理方法は、（Ａ）複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により囲まれる第１の面積と、上記平面において複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して第１の要素値列から補間計算される第２の要素値列により囲まれる第２の面積との差が、許容範囲内であるか判断する判断処理と、（Ｂ）判断処理を１又は複数の間隔について実施させることで、複数の散乱パラメータの全要素について上記差が許容範囲内となる間隔を特定する特定処理とを含む。
【００４９】
このようにすれば精度低下を許容範囲に抑えることができるような周波数間隔を特定できるようになる。これによって適切な回路シミュレーションなどの回路解析を実施することができるようになる。
【００５０】
さらに、上で述べた情報処理方法は、複数の第１の周波数における間隔を抽出する処理をさらに含むようにしても良い。この場合、上で述べた特定処理において、抽出された間隔のうち最も長い間隔から順番に判断処理を実施させるようにしても良い。このようにすれば、できるだけ長い間隔を特定できるようになるので、回路シミュレーション等の回路解析の処理負荷の増加を抑えることができるようになる。
【００５１】
また、本情報処理方法は、複数の第１の周波数における間隔が均等である場合には、第１の周波数の数が２の冪乗個となるように第１の周波数の最大値を決定する処理をさらに含むようにしても良い。周波数の数についても２の冪乗個となるように最大周波数を設定できるようになる。
【００５２】
また、本情報処理方法は、複数の第１の周波数の最大値又は設定された周波数の最大値まで、特定された間隔でサンプリングした場合のサンプリング数が２の冪乗個となっていない場合には、サンプリング数が２の冪乗個となるように、特定された間隔を短縮する処理をさらに含むようにしても良い。このように周波数の数についても２の冪乗個となるように間隔を調整するようにしても良い。
【００５３】
また、上で述べた第１の面積は、第１の要素値列のうち連続する２つの第１の要素値を直線で結んだ辺と、当該２つの第１の要素値から周波数の軸に下ろした垂線による辺と、当該２つの第１の要素値に対応する２つの第１の周波数の間隔に対応する辺とを有する台形を用いて算出される場合もある。また、第２の面積は、第２の要素値列のうち連続する２つの第２の要素値を直線で結んだ辺と、当該２つの第２の要素値から周波数の軸に下ろした垂線による辺と、当該２つの第２の要素値に対応する２つの第２の周波数の間隔に対応する辺とを有する台形を用いて算出されるようにしても良い。台形面積の総和を用いることで比較的演算量を抑えることができるようになる。
【００５４】
なお、上で述べたような処理をコンピュータに実施させるためのプログラムを作成することができ、当該プログラムは、例えばフレキシブル・ディスク、ＣＤ−ＲＯＭなどの光ディスク、光磁気ディスク、半導体メモリ（例えばＲＯＭ）、ハードディスク等のコンピュータ読み取り可能な記憶媒体又は記憶装置に格納される。
【００５５】
以上の実施例を含む実施形態に関し、さらに以下の付記を開示する。
【００５６】
（付記１）
複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により囲まれる第１の面積と、前記平面において前記複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して前記第１の要素値列から補間計算される第２の要素値列により囲まれる第２の面積との差が、許容範囲内であるか判断する判断処理と、
前記判断処理を１又は複数の間隔について実施させることで、前記複数の散乱パラメータの全要素について前記差が前記許容範囲内となる間隔を特定する特定処理と、
を、コンピュータに実行させるためのプログラム。
【００５７】
（付記２）
前記複数の第１の周波数における間隔を抽出する処理
をさらに前記コンピュータに実行させ、
前記特定処理において、抽出された前記間隔のうち最も長い間隔から順番に前記判断処理を実施させる
付記１記載のプログラム。
【００５８】
（付記３）
前記複数の第１の周波数における間隔が均等である場合には、第１の周波数の数が２の冪乗個となるように前記第１の周波数の最大値を決定する処理
をさらに前記コンピュータに実行させる付記２記載のプログラム。
【００５９】
（付記４）  
前記複数の第１の周波数の最大値又は設定された周波数の最大値まで、特定された前記間隔でサンプリングした場合のサンプリング数が２の冪乗個となっていない場合には、サンプリング数が２の冪乗個となるように、特定された前記間隔を短縮する処理
をさらに前記コンピュータに実行させるための付記１乃至３のいずれか１つ記載のプログラム。
【００６０】
（付記５）
前記第１の面積は、前記第１の要素値列のうち連続する２つの第１の要素値を直線で結んだ辺と、当該２つの第１の要素値から周波数の軸に下ろした垂線による辺と、当該２つの第１の要素値に対応する２つの第１の周波数の間隔に対応する辺とを有する台形を用いて算出され、
前記第２の面積は、前記第２の要素値列のうち連続する２つの第２の要素値を直線で結んだ辺と、当該２つの第２の要素値から周波数の軸に下ろした垂線による辺と、当該２つの第２の要素値に対応する２つの第２の周波数の間隔に対応する辺とを有する台形を用いて算出される
付記１乃至４のいずれか１つ記載のプログラム。
【００６１】
（付記６）
複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により囲まれる第１の面積と、前記平面において前記複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して前記第１の要素値列から補間計算される第２の要素値列により囲まれる第２の面積との差が、許容範囲内であるか判断する判断処理と、
前記判断処理を１又は複数の間隔について実施させることで、前記複数の散乱パラメータの全要素について前記差が前記許容範囲内となる間隔を特定する特定処理と、
を含み、コンピュータにより実行される情報処理方法。
【００６２】
（付記７）
複数の散乱パラメータの各要素について、周波数と当該要素とで張られる平面において複数の第１の周波数について予め規定されている当該要素の第１の要素値列により囲まれる第１の面積と、前記平面において前記複数の第１の周波数とは異なり且つ均一な間隔を有する複数の第２の周波数に対して前記第１の要素値列から補間計算される第２の要素値列により囲まれる第２の面積との差が、許容範囲内であるか判定する判定部と、
前記判断処理を１又は複数の間隔について実施させることで、前記複数の散乱パラメータの全要素について前記差が前記許容範囲内となる間隔を特定する特定部と、
を有する情報処理装置。
【符号の説明】
【００６３】
１１０  入力部
１２０  第１データ格納部
１３０  制御パラメータ決定部
１４０  第２データ格納部
１５０  出力部
【図１】
0005920039.tif 000002
【図２】
0005920039.tif 000003
【図３】
0005920039.tif 000004
【図４】
0005920039.tif 000005
【図５】
0005920039.tif 000006
【図６】
0005920039.tif 000007
【図７】
0005920039.tif 000008
【図８】
0005920039.tif 000009
【図９】
0005920039.tif 000010
【図１０】
0005920039.tif 000011
【図１１】
0005920039.tif 000012
【図１２】
0005920039.tif 000013
【図１３】
0005920039.tif 000014