TOP > 国内特許検索 > 回路設計支援方法及びその装置並びに回路設計支援プログラム

回路設計支援方法及びその装置並びに回路設計支援プログラム コモンズ 新技術説明会

国内特許コード P110003639
整理番号 Y04-P037
掲載日 2011年6月27日
出願番号 特願2005-516889
登録番号 特許第4571589号
出願日 平成17年1月7日(2005.1.7)
登録日 平成22年8月20日(2010.8.20)
国際出願番号 JP2005000097
国際公開番号 WO2005066853
国際出願日 平成17年1月7日(2005.1.7)
国際公開日 平成17年7月21日(2005.7.21)
優先権データ
  • 特願2004-003612 (2004.1.9) JP
発明者
  • 関屋 大雄
  • 丹治 裕一
  • 谷萩 隆嗣
  • 呂 建明
出願人
  • 国立研究開発法人科学技術振興機構
発明の名称 回路設計支援方法及びその装置並びに回路設計支援プログラム コモンズ 新技術説明会
発明の概要 【課題】 回路素子値や波形の導出を自動的に計算機で行うことができる様にする。
【解決手段】 設計対象のアナログ電子回路の回路素子値を、ヤコビ行列を含むニュートン法の繰り返し計算式が記述されたプログラムを計算機に実行させ算出する回路設計支援において、前記ヤコビ行列中の要素である偏微分を近似式で置き換えた前記計算式のプログラムを用いる。更に、ヤコビ行列中の前記偏微分の対象となる要素を、前記アナログ電子回路の回路構成の応答として観測される波形から求める。
【選択図】 図1
従来技術、競合技術の概要


近年の情報化社会の発展には著しいものがあり、その中で,計算機のCPUはGHzの動作(クロック)周波数が求められ、無線通信においてはGHz帯の周波数の使用が始まるなど、電子回路の高周波数化が顕著になってきている。これに伴い、CPUなどに電力を供給するための電源回路や無線通信における増幅器等のアナログ回路の開発が極めて重要になってきている。



高周波数化の進んだ電子回路に接続される電源回路や増幅器等のアナログ回路の設計は、
(1)回路構成の決定
(2)回路素子値の決定
(3)回路動作の評価
の3段階に分類される。



これらの設計を計算機上で支援するツールとして「SPICE」などの回路シミュレータが開発されており、これは、上記の(3)に対しての強力な支援ツールとして広く用いられている。しかし、回路動作をシミュレートするためには、回路素子値を与える必要があり、所望の動作を実現するための素子値を導出する上記(2)の作業を欠かすことができない。



アナログ回路の開発者は、回路の素子値を決定するために、初めに手計算で回路解析を行い、その解析結果を用いて回路の素子値を決定している。回路解析には近似を多く含むため、導出した素子値は誤差を含む。従って、導出した素子値をシミュレータに入力し、その結果から、回路の素子値を経験的に調整するという上記(2)と(3)の繰り返し作業を行うことになる。また、素子値の調整は設計者の経験がものをいう分野であるため、設計者の教育も必要となる。



一方、シミュレータを用いてモンテカルロ法で素子値を求めたとしても、満足すべき条件数の増加に伴い、計算時間が指数関数的に増加してしまうという問題もある。つまり、上記(1)と(3)の間の接続作業である上記(2)にあたる作業が、回路開発の時間的ボトルネックとなっている。そのため、経験が浅い設計者でも、回路設計を容易かつ高速に行える設計支援ツールが強く求められている。



そこで、本発明者等は、先に、下記非特許文献1に記載の様に、計算機上で数値的に素子値を導出する素子値算出方法を提案した。以下、この従来の素子値算出方法を説明する。



例えば、次の数1に示される微分方程式で記述される回路について考える。



【数1】




ここで、t∈R,x∈R、及び、λ∈Rは、夫々、時間,n次元状態空間,m次元のシステムパラメータを示す。簡単のため、



【数2】




はC級写像であり、時間領域において周期tの周期性を持つと仮定する。つまり、



【数3】




と表される。更に、上記の数1は、-∞<t<∞において、任意の初期値x∈Rおよびあらゆるシステムパラメータλ∈R:x(0)=φ(0,x,λ)=xに対し、x(t)=φ(t,x,λ)という解を持つとする。



上記の数3で述べた周期性の仮定により、状態空間Rからそれ自身に戻るC級離散写像T:



【数4】




を定義することができる。この写像Tは、ポアンカレ写像と呼ばれる。



数1の微分方程式の解x(t)=φ(t,p,λ)が周期tの周期性を有していた場合、



【数5】




を満たす点p∈Rは写像Tに対する不動点という。



電子回路を設計するには、拘束条件を考えることが多い。もし、条件の数をN(≦m)とすると、時刻t=tckにおける夫々の条件をgと表すことにより、



【数6】




と表すことができる。この場合、設計パラメータとしてN個のパラメータを用意することができる。従って、その他の(m-N)個のパラメータは、設計仕様として与えなくてはならない。



以上より、電子回路の設計は、代数方程式である上記の数5,数6の解法に帰着する。これらの代数方程式(数5,数6)を改めて



【数7】




と書き直す。ここで、T(x,λ),x,λは、夫々、T(x,λ)=〔T(x,λ),T(x,λ),…,T(x,λ)〕,x=x(0)=〔x(0),x(0),…,x(0)〕,λ=〔λ,λ,…,λと表すことができる。更に、λ∈R



【数8】




と定義する。



設計方程式(数7)を解くには、代数方程式を解くための最も一般的アルゴリズムであるニュートン法を用いる。数7中の未知数をu∈Rn+N:u=〔x,λと表すと、繰り返し計算



【数9】




を終了条件∥uk+1-u∥<δを満足するまで行う。ここで、F’∈R(n+N)×(n+N)はFのヤコビ行列つまり、



【数10】




を表す。更に、kは繰り返しの回数であり、δ≪1である。ここでは、δ=10-9とする。このとき、uk+1が数7の解となる。



数9を計算するためには数10中のT(u)を求める必要があるが、このT(u)は、数1にルンゲクッタ法を適用することで求めることができる。更に、数10のヤコビ行列F’(u)の要素は、一次変分方程式



【数11】




を解くことにより導出する。数11は、∂φ(t,u)/∂x,∂φ(t,u)/∂λの微分方程式と見なせる。つまり、数11を0からtまたはtckまでルンゲクッタ法を用いて解けば、∂φ(t,u)/∂x,∂φ(tck,u)/∂x,∂φ(t,u)∂λ,∂φ(tck,u)/∂λを導出することができる。Tはφの関数であり、gは一般的にφの関数として表すことができるため、これらの値から、ヤコビ行列の要素であるT(u)/∂x,T(u)/∂λ,g(tck,u)/∂x,g(tck,u)/∂λを数値的に導出することができる。



従来は、以上の計算により、未知係数uを求め、その結果、電子回路の設計値であるλを決定していた。



図2は、上述した従来の計算手順を用いて回路設計を行う処理手順を示すフローチャートである。先ず、電子回路の構成を決定する回路構成決定ステップ1と、決定された回路構成から回路設計者が手作業で回路方程式(数1)を決定する回路方程式構築ステップ2と、この回路方程式から変分方程式(数11)を回路設計者が手作業で導出する変分方程式導出ステップ3と、制約条件(数7)を導出する制約条件導出ステップ4と、回路方程式(数1)と変分方程式(数11)を用いて条件式(数7)を解くことによって回路素子値と波形とを導出するステップ5とからなり、ステップ5で満足のいく解が得られないときはステップ1に戻る。



【非特許文献1】
IEEE TRANSACTION ON CIRCUITS AND SYSTEM-1: FUNDAMENTAL THEORY AND APPLICATIONS.VOL.49.NO.7, JULY 2002

産業上の利用分野


本発明は回路設計支援方法及びその装置並びに回路設計支援プログラムに係り、特に、微分方程式でなる回路方程式を与えることなく回路設計を設計者と計算機との間で対話的に自動で行うことが可能な回路設計支援方法及びその装置並びに回路設計支援プログラムに関する。

特許請求の範囲 【請求項1】
設計対象のアナログ電子回路の設計に関わる拘束条件を、回路素子値を未知変数とする代数方程式で表し、該代数方程式を、計算手段がニュートン法の繰り返し計算で解く回路シミュレータの回路設計支援方法において、前記ニュートン法の繰り返し計算で用いるヤコビ行列の要素である偏導関数の値を計算するための要素を、前記回路シミュレータで観測される前記アナログ電子回路の回路応答としての出力波形から求めると共に、前記計算手段は、前記ヤコビ行列の要素である前記偏導関数を、該偏導関数の定義に基づいた微小区間における平均変化値を求める近似式で代用して計算を行う回路設計支援方法。

【請求項2】
請求項1において、前記拘束条件として、前記出力波形の観測値の最大値,最小値,平均値を用いる回路設計支援方法。

【請求項3】
設計対象のアナログ電子回路の設計に関わる拘束条件を、回路素子値を未知変数とする代数方程式で表し、該代数方程式を、ニュートン法の繰り返し計算で解く回路設計支援装置において、前記ニュートン法の繰り返し計算で用いるヤコビ行列の要素である偏導関数の値を計算するための要素を、該回路設計支援装置で観測される前記アナログ電子回路の回路応答としての出力波形から求めると共に、該ヤコビ行列の要素である前記偏導関数を、該偏導関数の定義に基づいた微小区間における平均変化値を求める近似式で代用する計算手段を備える回路設計支援装置。

【請求項4】
請求項3において、前記拘束条件として、前記出力波形の観測値の最大値,最小値,平均値を用いる回路設計支援装置。

【請求項5】
設計対象のアナログ電子回路の設計に関わる拘束条件を、回路素子値を未知変数とする代数方程式で表し、該代数方程式を、計算手段がニュートン法の繰り返し計算で解く回路シミュレータの回路設計支援プログラムにおいて、前記ニュートン法の繰り返し計算で用いるヤコビ行列の要素である偏導関数の値を計算するための要素を、前記回路シミュレータで観測される前記アナログ電子回路の回路応答としての出力波形から求めると共に、前記ヤコビ行列の要素である前記偏導関数を、該偏導関数の定義に基づいた微小区間における平均変化値を求める近似式で代用して前記計算手段が計算を行うステップを備える回路設計支援プログラム。

【請求項6】
請求項5において、前記拘束条件として、前記出力波形の観測値の最大値,最小値,平均値を用いる回路設計支援プログラム。
国際特許分類(IPC)
Fターム
出願権利状態 登録
ライセンスをご希望の方、特許の内容に興味を持たれた方は、問合せボタンを押してください。


PAGE TOP

close
close
close
close
close
close
close