(57)【要約】【課題】アナログ回路やアナログ回路とディジタル回路が混在する回路において、パルス信号や急激な波形変化を伴う信号波形の過渡解析を数値演算によってシミュレーションし、解析時間を短縮する。【解決手段】数値演算を行うための回路ブロックのモデル化の第１の形態は、回路ブロック内の電圧および電流を変数とし、この変数を微分した微分変数を含む状態方程式を用いて記述して数値演算のための演算ブロックを形成し、モデル化の第２の形態は、電気回路による重ねの理を適用し、回路ブロックの伝達関数と電圧源又は電流源の積の和を求め、この積の和を重ねの理によってアナログ能動要素の入力信号として記述し、入力信号の伝達関数とアナログ能動要素の伝達関数との積によって回路ブロックの総合的な伝達関数を記述して数値演算のための演算ブロックを形成する。【選択図】図１
2010092434.tif 000002
【特許請求の範囲】
【請求項１】
アナログ回路の動作を数値演算によってシミュレートする方法であって、
アナログ回路を複数の回路ブロックに分割し、
当該分割した回路ブロックにおいて、
出力を求めるために時間変化を解析する過渡解析を要する回路ブロックについて、
回路ブロック内の電圧および電流を変数とし、前記変数を微分した微分変数を含む状態方程式を用いて記述して数値演算のための演算ブロックを形成し、当該演算ブロックによって回路ブロックをモデル化し、微分変数を数値演算で解くことによって当該回路ブロックの出力の時間変化を求め、
前記過渡解析を要さない回路ブロックについて、
当該回路ブロックのアナログ能動要素に入力する電圧源および又は電流源について、前記電圧源又は電流源が単独で存在するときの伝達関数を、当該電圧源又は電流源からアナログ能動要素までをＳ関数で記述し、当該伝達関数とＳ関数で記述した前記電圧源の電圧又は電流源の電流との積の和を求め、この積の和を重ねの理によってアナログ能動要素の入力信号として記述し、
前記アナログ能動要素の伝達関数をＳ関数で記述し、
前記入力信号の伝達関数と前記アナログ能動要素の伝達関数との積によって回路ブロックの総合的な伝達関数を記述して数値演算のための演算ブロックを形成し、当該演算ブロックによって回路ブロックをモデル化し、前記総合的な伝達関数をラプラス逆変換することによって回路ブロックの出力の時間変化を求めることを特徴とする、アナログ回路の数値演算によるシミュレーション方法。
【請求項２】
入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路と
前記電圧パルスを平滑化して負荷に供給する出力平滑回路と、
前記出力平滑回路の出力電圧と参照電圧の差分電圧を入力信号とし、当該入力信号を電圧増幅し、低域通過させた電圧を出力し、リミタ電流を出力端から入力端に帰還する誤差増幅器と、
前記出力平滑回路に流れる電流を検出し、当該電流を電圧変換して出力する前記電流検出回路と、
誤差増幅器出力電圧と電流検出回路の出力電圧とを比較する比較器を備え、
前記スイッチング回路のスイッチング期間を制御するディジタル制御回路とを備えるスイッチング電源回路の動作を数値演算によってシミュレートする方法であって、
前記出力平滑回路は、前記スイッチング回路を構成する出力トランジスタの出力電圧と負荷電流とを入力信号とし、
前記出力平滑回路の出力電圧と当該出力平滑回路が備えるインダクタンス素子に流れる電流とを出力信号とし、
当該出力平滑回路が備える容量素子の両端の電圧と前記インダクタンス素子に流れる電流を変数とし、前記変数を微分した微分変数を含む状態方程式を用いて前記出力平滑回路を記述して数値演算を行う演算ブロックを形成し、当該演算ブロックを用いてスイッチング回路および出力平滑回路をモデル化し、
前記インダクタンス素子を流れる電流および前記出力平滑回路の出力直流電圧の時間変化を求め、
前記誤差増幅器は、
入力信号電圧源から前記誤差増幅器の入力端までの伝達関数と、
出力信号電圧源から前記誤差増幅器の入力端までの伝達関数と、
参照電圧源から前記誤差増幅器の入力端までの伝達関数と、
入力端と出力端との間に接続されたダイオードリミタに流れる電流を表すリミタ電流源から前記誤差増幅器の入力端までの伝達関数とについて、それぞれが単独で存在するときの伝達関数をＳ関数で記述し、
当該各伝達関数とＳ関数で記述した前記入力信号電圧源、出力信号電圧源、参照電圧源、リミタ電流源の電圧又は電流源の電流との積の和を求め、この積の和を重ねの理によって誤差増幅器の入力信号として記述し、
前記誤差増幅器の伝達関数をＳ関数で記述し、
前記入力信号のＳ関数と前記誤差増幅器の伝達関数との積によって誤差増幅器の総合的な伝達関数を記述して数値演算を行う演算ブロックを形成し、当該演算ブロックを用いて誤差増幅器をモデル化し、
前記総合的な伝達関数をラプラス逆変換することによって誤差増幅器の出力の時間変化を求め、
前記電流検出回路は、
前記スイッチングトランジスタに流れる電流に対して所定比率の電流を流す検出用トランジスタを前記スイッチングトランジスタに並列接続し、
前記検出用トランジスタに流れる電流と抵抗との積によって記述して数値演算を行う演算ブロックを形成し、当該演算ブロックを用いて電流検出回路をモデル化し、
前記所定比率と前記抵抗との積を一定値の電流−電圧変換係数として数値演算を行うことによって電流を検出することによって、
前記スイッチング回路と前記出力平滑回路と前記誤差増幅器と前記電流検出回路の各アナログ回路の動作を数値演算によってシミュレートすることを特徴とする、スイッチング電源回路の数値演算によるシミュレーション方法。
【請求項３】
入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路と、
前記電圧パルスを平滑化して負荷に供給する出力平滑回路と、
前記出力平滑回路の出力電圧と参照電圧の差分電圧を入力信号とし、当該入力信号を電圧増幅し、低域通過させた電圧を出力し、リミタ電流を出力端から入力端に帰還する誤差増幅器と、
前記出力平滑回路に流れる電流を検出し、当該電流を電圧変換して前記電流検出回路と、
前記スイッチング回路のスイッチング期間を制御するディジタル制御回路とを備えるスイッチング電源回路において、
入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路と、前記電圧パルスを平滑化して負荷に供給する出力平滑回路、前記誤差増幅器、前記電流検出回路、および前記ディジタル制御回路を含む開ループの利得および位相の周波数特性を数値演算によってシミュレートする方法であって、
前記スイッチング回路を構成する出力トランジスタの出力電圧と負荷電流とを入力信号とし、
前記出力平滑回路の出力電圧と当該出力平滑回路が備えるインダクタンス素子に流れる電流とを出力信号とし、
当該出力平滑回路が備える容量素子の両端の電圧と前記インダクタンス素子に流れる電流を変数とし、前記変数を微分した微分変数を含む状態方程式を用いて前記出力平滑回路を記述して数値演算を行う演算ブロックを形成し、当該演算ブロックを用いてスイッチング回路および出力平滑回路をモデル化し、
前記インダクタンス素子を流れる電流および前記出力平滑回路の出力直流電圧の時間変化を求めることを特徴とする、出力平滑回路の数値演算によるシミュレーション方法。
【請求項４】
入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路と、
前記電圧パルスを平滑化して負荷に供給する出力平滑回路と、
前記出力平滑回路の出力電圧と参照電圧の差分電圧を入力信号とし、当該入力信号を電圧増幅し、低域通過させた電圧を出力し、リミタ電流を出力端から入力端に帰還する誤差増幅器と、
前記出力平滑回路に流れる電流を検出し、当該電流を電圧変換して前記電流検出回路と、
前記スイッチング回路のスイッチング期間を制御するディジタル制御回路とを備えるスイッチング電源回路において、
入力信号電圧と参照電圧の差分電圧を入力信号とし、当該入力信号を電圧増幅し、低域通過させた電圧を出力し、リミタ電流を出力端から入力端に帰還する誤差増幅器の周波数特性を数値演算によってシミュレートする方法であって、
入力信号電圧源から前記誤差増幅器の入力端までの伝達関数と、
出力信号電圧源から前記誤差増幅器の入力端までの伝達関数と、
参照電圧源から前記誤差増幅器の入力端までの伝達関数と、
入力端と出力端との間に接続されたダイオードリミタに流れる電流を表すリミタ電流源から前記誤差増幅器の入力端までの伝達関数とについて、それぞれが単独で存在するときの伝達関数をＳ関数で記述し、
当該各伝達関数とＳ関数で記述した前記入力信号電圧源、出力信号電圧源、参照電圧源、リミタ電流源の電圧又は電流源の電流との積の和を求め、この積の和を重ねの理によって誤差増幅器の入力信号として記述し、
前記誤差増幅器の伝達関数をＳ関数で記述し、
前記入力信号のＳ関数と前記誤差増幅器の伝達関数との積によって誤差増幅器の総合的な伝達関数を記述して数値演算を行う演算ブロックを形成し、当該演算ブロックを用いて誤差増幅器をモデル化し、
前記総合的な伝達関数をラプラス逆変換することによって誤差増幅器の出力の時間変化を求めることを特徴とする、誤差増幅器の数値演算によるシミュレーション方法。
【請求項５】
入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路と、
前記電圧パルスを平滑化して負荷に供給する出力平滑回路と、
前記出力平滑回路の出力電圧と参照電圧の差分電圧を入力信号とし、当該入力信号を電圧増幅し、低域通過させた電圧を出力し、リミタ電流を出力端から入力端に帰還する誤差増幅器と、
前記出力平滑回路に流れる電流を検出し、当該電流を電圧変換して前記電流検出回路と、
前記スイッチング回路のスイッチング期間を制御するディジタル制御回路とを備えるスイッチング電源回路において、
入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路のスイッチング期間を制御するディジタル制御回路の制御動作を数値演算によってシミュレートする方法であって、
電流検出回路の出力電圧と、誤差増幅器の出力電圧との差分を零と比較して論理値を出力し、
当該比較結果の論理値によって前記スイッチング期間を制御するディジタル制御回路の制御動作をブール代数によって記述して数値演算を行う演算ブロックを形成し、当該演算ブロックを用いてディジタル制御回路をモデル化し、
前記ブール代数の数値演算によってスイッチング期間を制御することを特徴とする、ディジタル制御回路の数値演算によるシミュレーション方法。
【請求項６】
入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路と、
前記電圧パルスを平滑化して負荷に供給する出力平滑回路と、
前記出力平滑回路の出力電圧と参照電圧の差分電圧を入力信号とし、当該入力信号を電圧増幅し、低域通過させた電圧を出力し、リミタ電流を出力端から入力端に帰還する誤差増幅器と、
前記出力平滑回路に流れる電流を検出し、当該電流を電圧変換して前記電流検出回路と、
前記スイッチング回路のスイッチング期間を制御するディジタル制御回路とを備えるスイッチング電源回路において、
スイッチング回路を構成するスイッチングトランジスタのスイッチングによって、入力直流電圧をスイッチングして電圧パルスを出力し、当該電圧パルスを平滑化して得られる電流を検出し、当該電流を電圧変換して出力する電流検出回路の動作を数値演算によってシミュレートする方法であって、
前記スイッチングトランジスタに流れる電流に対して所定比率の電流を流す検出用トランジスタを前記スイッチングトランジスタに並列接続し、
前記検出用トランジスタに流れる電流と抵抗との積によって記述して数値演算を行う演算ブロックを形成し、当該演算ブロックを用いて電流検出回路をモデル化し、
前記所定比率と前記抵抗との積を一定値の電流−電圧変換係数として数値演算を行うことによって電流を検出することを特徴とする、電流検出回路の数値演算によるシミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、数値演算によるシミュレーション方法に関し、アナログ回路、およびアナログ回路とディジタル回路が混在するＤＣ−ＤＣコンバータに代表されるスイッチング電源回路において周波数特性を数値演算によってシミュレートする方法に関する。
【背景技術】
【０００２】
アナログ電子回路のシミュレーション方法として、回路方程式をガウス消去法、ＬＵ分解法を用いて解く他、非線形回路の解析ではNewton-Rahpson法を用いる方法が知られている。また、回路の動特性をシミュレートするには、回路中の素子を時間刻み毎に変化する電流源、電圧源に置き換え、この時間刻み単位で非線形回路計算を繰り返すことが行われている。また、ディジタル回路のシミュレーション方法としてイベントドリブン法と呼ばれる方法が知られている（特許文献１）。
【０００３】
上記した特許文献１には、アナログディジタル混在回路において、アナログ部のシミュレーションは通常高精度を目標として作られており時間刻みの間隔を制御する必要があるため計算時間が長くなるという課題を指摘し、この課題を解決する方法として、Ｓ関数で記述されたアナログ回路部の伝達関数をシミュレーション精度に応じた周期を持つサンプリング周期によるＺ関数に変換し、変換したＺ関数に基づいてディジタル回路網を形成し、論理シミュレーションを行うことが開示されている。
【０００４】
また、アナログ回路をシミュレートするプログラムとしてSPICE（登録商標）と呼ばれる回路シミュレーションプログラムが知られ、数値演算プログラムとしてMATLAB（登録商標）と呼ばれるプログラムが知られている。
【特許文献１】特許第３０１９３６８号（従来技術の項、発明が解決しようとする課題の項、課題を解決するための手段の項）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
前記した、アナログ回路やアナログディジタル混在回路のシミュレーション方法では、アナログ回路部の動作を論理シミュレーションするために、アナログ回路をＳ関数による伝達関数で表現し、Ｓ関数をＺ関数に変換することによってディジタル回路を形成する手法を用いている。
【０００６】
しかしながら、パルス波形や急激な電圧変化や電流変化を持つ信号波形を扱う回路では、その回路構成を伝達関数で表現することが難しい場合があり、また、その過渡解析を行うには急激に変化する部分を細かく解析するために解析ステップ時間を細かく設定する必要があるためシミュレーション時間が増大するという問題がある。また、複数の信号源が接続されるような回路構成では、伝達関数で表現することが難しい場合がある。
【０００７】
特に、アナログディジタル混在回路のように、ディジタル回路の解析のために数値演算によってシミュレートを行う必要がある場合には、アナログ回路の動作解析を数値演算によるシミュレートに適した構成とすることが求められる。この際、アナログ回路部分において、パルス波形や急激な電圧変化や電流変化を持つ信号波形のシミュレーションに要する時間を短縮することが求められ、また、アナログ回路の構成を数値演算のシミュレーションに適したモデル化が成されることが求められる。
【０００８】
以下、この課題について、ＤＣ−ＤＣコンバータを例にして説明する。一般に、ＤＣ−ＤＣコンバータは、入力直流電圧Ｖinをスイッチングして電圧パルス波形に変換し、この電圧パルス波形をインダクタＬとコンデンサＣとにより平滑することで入力電圧とは異なる所望の直流出力電圧Ｖoutを得る構成となっている。
【０００９】
上記電圧パルス波形のパルス時間幅ｔBとパルス波形の繰り返し周期Ｔとの比ｔB／Ｔをデューティ比Ｄで表すと、ＤＣ−ＤＣコンバータの直流出力電圧Ｖoutは、降圧形（Ｖin＞Ｖout）の場合、
【数１】
2010092434.tif 000003
のようになり、昇圧形（Ｖin＜Ｖout）の場合には
【数２】
2010092434.tif 000004
のようになる。
【００１０】
式（１）あるいは式（２）の動作において、入力直流電圧Ｖinをスイッチして電圧パルス波形として平滑回路に送り出す役目は、スイッチング回路を構成する出力トランジスタが担っている。大きな電流値を供給しなければならないような場合であっても、スイッチがオンとなっている状態の出力トランジスタの両端の電圧は小さいので、出力トランジスタにおける電力消費は少ない。この理由によりＤＣ−ＤＣコンバータは電力変換効率の高い変換器として多用されている。
【００１１】
しかしながらＤＣ−ＤＣコンバータに要求される電流供給能力は負荷の状況により常に変動する。負荷がスタンバイモードになった場合には供給する電流値は小さくなり、負荷がアクティブな場合にはその電流値は大となる。電流値が変動すると、直流出力電圧Ｖoutも変動する。そこで負帰還をかけてＶoutを安定にする必要がある。
【００１２】
図１２は、ＬＳＩ化降圧形ＤＣ−ＤＣコンバータ（Ｖin＞Ｖout）の構成図を示している。図１２に示すＤＣ−ＤＣコンバータは、入力直流電圧Ｖinをスイッチングするスイッチング回路６０、スイッチング回路６０の出力を平滑化して直流出力電圧Ｖoutを形成する出力平滑回路２０、直流出力電圧Ｖoutを定数（Ｋ）倍した電圧ＫＶoutと参照電圧Ｖrefとの電圧差（Δ）を増幅する誤差増幅器部４０、出力平滑回路２０のインダクタ２１に流れる電流ｉLを形成する電流検出器５０、誤差増幅器部４０の出力電圧ｖcと電流検出器５０で検出した電流にスロープ補償電流を加えて電圧変換した電圧値ｖcfbとを比する比較器３０、比較器３０の出力ｖfを帰還してスイッチング回路６０のスイッチング動作を制御するディジタル制御回路１０によって構成されている。ディジタル制御回路１０は発振器（ＯＳＣ）１１の周期信号で駆動される。
【００１３】
出力電圧Ｖoutを定数（Ｋ）倍した電圧ＫＶoutと参照電圧Ｖrefの電圧差Ｅ（Δ）を誤差増幅器部４０で増幅して、その出力電圧ｖcを比較器３０に送る。出力電圧Ｖoutを定数（Ｋ）倍した電圧ＫＶoutが参照電圧Ｖrefの電圧より大である場合には、誤差増幅器部４０の出力電圧ｖcは小となる。
【００１４】
一方、出力平滑回路２０のインダクタ２１（Ｌ）に流れる電流ｉLは、スイッチング回路６０の出力トランジスタＭｐあるいは出力トランジスタＭｎに流れる電流と一致している。電流ｉLの電流値の変化は誤差増幅器部４０では十分検出できないため、この電流ｉLの電流値の変化は電流検出器５０で検出し更にスロープ補償５１を施して比較器３０に入力される。
【００１５】
以上のように、直流出力電圧の電圧情報（Ｖout）およびインダクタの電流情報（ｉL）は比較器３０を通してディジタル制御回路１０へと帰還され、該ディジタル制御回路１０にてスイッチング回路６０の出力トランジスタＭｐあるいはＭｎのオン期間を制御する。これにより直流出力電圧（Ｖout）が変化した場合、あるいは出力電流（ｉL）が変化して直流出力電圧（Ｖout）が変化した場合でも、直流出力電圧（Ｖout）は常に一定に保たれる。
【００１６】
図１２の回路において、誤差増幅器部４０のアナログ回路は主に変化の緩やかな波形を扱うが、出力平滑回路２０や電流検出器５０等の他のアナログ回路ではパルス波形や急激な変化のある波形を扱っている。
【００１７】
ここで、アナログ回路の回路シミュレーションツールを用いて、急激な電圧変化あるいは電流変化を持つ信号波形を過渡解析すると、回路シミュレーションツールは急激な変化のある部分で細かく精密な解析を行おうとしてその解析ステップ時間を非常に細かく設定する。このため過渡解析に要するシミュレーション時間が急激に増大する。
【００１８】
さらにアナログ回路シミュレーションツールはトランジスタをベースとしたシミュレーションツールであるから、図１２の各回路構成に含まれる全ての素子を含んだ形で解析を行うため、回路規模が大となる。
【００１９】
上述したように、パルス的な動作を解析するために、解析の時間ステップを非常に細かく設定する必要があること、回路を構成する各素子を基本とした解析をする必要があるため回路規模が大となることの２つの要因によって、図１２に示したような回路の過渡解析を行うと解析時間に長時間を要するという課題がある。
【００２０】
また図１２ではパルス的な信号波形を扱い、多くのスイッチやディジタル回路機能が含まれるので、アナログシミュレータで用いられる小信号モデルに基づく伝達関数を使ったＡＣ解析を行うことができない。開ループ特性を求める場合には、直流出力電圧Ｖoutと誤差増幅器部４０のマイナス入力端子間に５０Ω程の低抵抗を接続し、この低抵抗の両端に誤差増幅器部４０のマイナス入力端子側が正となるような信号を与えて過渡解析を行う。このとき直流出力電圧Ｖoutに重畳して入力信号成分が現れるので、過渡解析結果に含まれる開ループを通過してきた入力信号成分を抽出する。
【００２１】
ただし上記の操作では一つの周波数の信号に対する開ループの応答しか得られない。開ループの周波数特性を求めるには、開ループに様々な周波数の信号を入力してそれぞれの周波数に対する応答のデータを集める必要がある。一つの過渡解析に要する解析時間が大である上、さらに、周波数特性を求める場合には、この動作を繰り返して行う必要があり、膨大な時間を費やさないと求まらないという課題がある。
【００２２】
そこで、本発明は上述した課題を解決して、アナログ回路やアナログ回路とディジタル回路が混在する回路において、パルス信号や急激な波形変化を伴う信号波形の過渡解析を数値演算によってシミュレーションすることを目的とし、解析時間を短縮することを目的とする。
【課題を解決するための手段】
【００２３】
本発明は、アナログ回路やアナログ回路とディジタル回路が混在する回路の動作を数値演算によってシミュレートする第１の態様、アナログ回路とディジタル回路が混在するスイッチング電源回路の動作を数値演算によってシミュレートする第２の態様とを含む。
【００２４】
本発明の第１の態様は、アナログ回路やアナログ回路とディジタル回路が混在する回路の動作を数値演算によってシミュレートする方法であり、アナログ回路を複数の回路ブロックに分割し、この分割した回路ブロックにおいて２つの形態によって数値演算を行う回路ブロックをモデル化し、モデル化した回路ブロックを数値演算で解くことによって回路ブロックの出力の時間変化を求める。これによって、アナログ回路やアナログ回路とディジタル回路が混在する回路において、パルス信号や急激な波形変化を伴う信号波形の過渡解析の数値演算によるシミュレーションを可能とし、また、解析時間を短縮することができる。
【００２５】
数値演算を行うための回路ブロックのモデル化の第１の形態は、出力に時間変化を解析する過渡解析を要する回路ブロックに適用するものである。この形態のモデル化では、回路ブロック内の電圧および電流を変数とし、この変数を微分した微分変数を含む状態方程式を用いて記述して数値演算のための演算ブロックを形成し、形成した演算ブロックによって回路ブロックをモデル化する。
【００２６】
モデル化した回路ブロックにおいて、微分変数を数値演算で解くことによって出力の時間変化を求める。
【００２７】
状態方程式は微分変数を含んで記述されるため、時間変化を数値演算で求めることができる。また、数値演算であるため、膨大な解析ステップを繰り返すことなくパルス信号や急激な波形変化を伴う信号波形の過渡解析を数値演算でシミュレートすることができる。
【００２８】
数値演算を行うための回路ブロックのモデル化の第２の形態は、前記過渡解析を要さない回路ブロックに適用するものである。この形態のモデル化は、電気回路による重ねの理を適用するものであり、回路ブロックのアナログ能動要素に入力する電圧源および又は電流源について、電圧源又は電流源が単独で存在するときの伝達関数を、電圧源又は電流源からアナログ能動要素までをＳ関数で記述し、この回路ブロックの伝達関数とＳ関数で記述した電圧源の電圧又は電流源の電流との積の和を求め、この積の和を重ねの理によってアナログ能動要素の入力信号として記述する。さらに、アナログ能動要素の伝達関数をＳ関数で記述し、入力信号の伝達関数とアナログ能動要素の伝達関数との積によって回路ブロックの総合的な伝達関数を記述して数値演算のための演算ブロックを形成し、この演算ブロックによって回路ブロックをモデル化し、この総合的な伝達関数をラプラス逆変換することによって回路ブロックの出力の時間変化を求める。
【００２９】
本発明の第２の態様は、アナログ回路とディジタル回路が混在するスイッチング電源回路の動作を数値演算によってシミュレートする方法である。
【００３０】
本発明のシミュレーション方法を適用するスイッチング電源回路は、入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路と、電圧パルスを平滑化して負荷に供給する出力平滑回路と、出力平滑回路の出力電圧の定数倍の電圧と参照電圧の差分電圧を入力信号とし、この入力信号を電圧増幅し、低域通過させた電圧を出力し、リミタ電流を出力端から入力端に帰還する誤差増幅器と、出力平滑回路に流れる電流を検出し、この電流を電圧変換して比較器に帰還する電流検出回路と、スイッチング回路のスイッチング期間を制御するディジタル制御回路とを備える。
【００３１】
上記したスイッチング電源回路において、出力平滑回路、誤差増幅器、電流検出回路、およびディジタル制御回路をモデル化して数値演算によって信号の時間変化を求める。
【００３２】
出力平滑回路は、スイッチング回路を構成する出力トランジスタの出力電圧と負荷電流とを入力信号とし、出力平滑回路の出力電圧と出力平滑回路が備えるインダクタンス素子に流れる電流とを出力信号とし、出力平滑回路が備える容量素子の両端の電圧とインダクタンス素子に流れる電流を変数とし、この変数を微分した微分変数を含む状態方程式を用いて出力平滑回路を記述して数値演算を行う演算ブロックを形成し、形成した演算ブロックを用いてスイッチング回路および出力平滑回路をモデル化し、このモデル化においてインダクタンス素子を流れる電流および出力平滑回路の出力直流電圧の時間変化を求める。
【００３３】
このモデル化によれば、状態方程式において、微分した変数を積分することは変数の時間変化を求めることに対応するため、微分変数を含む状態方程式を数値演算で解くことによって、膨大な解析ステップを繰り返すことなくパルス信号や急激な波形変化を伴う信号波形の過渡解析を数値演算でシミュレートすることができる。
【００３４】
誤差増幅器は、入力信号電圧源から前記誤差増幅器の入力端までの伝達関数と、出力信号電圧源から前記誤差増幅器の入力端までの伝達関数と、参照電圧源から前記誤差増幅器の入力端までの伝達関数と、入力端と出力端との間に接続されたダイオードリミタに流れる電流を表すリミタ電流源から前記誤差増幅器の入力端までの伝達関数とについて、それぞれが単独で存在するときの伝達関数をＳ関数で記述し、各伝達関数とＳ関数で記述した入力信号電圧源、出力信号電圧源、参照電圧源、リミタ電流源の電圧又は電流源の電流との積の和を求め、この積の和を重ねの理によって誤差増幅器の入力信号として記述する。
【００３５】
さらに、誤差増幅器の伝達関数をＳ関数で記述し、入力信号の伝達関数と誤差増幅器の伝達関数との積によって誤差増幅器の総合的な伝達関数を記述して数値演算を行う演算ブロックを形成して誤差増幅器をモデル化し、このモデル化において総合的な伝達関数をラプラス逆変換することによって誤差増幅器の出力の時間変化を求める。
【００３６】
ここで、誤差増幅器の信号波形は急激な波形変化を伴わないため、伝達関数によってモデル化することができる。本発明では、複数の外部要素の伝達関数を組み合わせることによってモデル化する。より詳細には、重ねの理を用い、それぞれが単独で存在するときの伝達関数と入力信号の伝達関数との積の和によって求める。
【００３７】
電流検出回路は、スイッチング回路のスイッチングトランジスタに流れる電流に対して所定比率の電流を流す検出用トランジスタをスイッチングトランジスタに並列接続し、この検出用トランジスタに流れる電流と抵抗との積によって記述して数値演算を行う演算ブロックを形成し、形成した演算ブロックを用いて電流検出回路をモデル化する。
【００３８】
このモデル化において、所定比率と抵抗との積を一定値の電流−電圧変換係数として数値演算を行うことによって電流を検出する。
【００３９】
上記した各アナログ回路のモデル化によって、スイッチング回路と、出力平滑回路と誤差増幅器と電流検出回路の各アナログ回路の動作を数値演算によってシミュレートすることができる。
【００４０】
また、スイッチング電源回路において、入力直流電圧をスイッチングして電圧パルスを出力するスイッチング回路のスイッチング期間を制御するディジタル制御回路の制御動作を数値演算によってシミュレートする構成では、電流検出回路の出力電圧と、誤差増幅器の出力電圧との差分を零と比較して論理値を出力し、この比較結果の論理値によって前記スイッチング期間を制御するディジタル制御回路の制御動作をブール代数によって記述して数値演算を行う演算ブロックを形成し、この演算ブロックを用いてディジタル制御回路をモデル化する。このモデル化において、ブール代数の数値演算によってスイッチング期間を制御する。
【発明の効果】
【００４１】
本発明によれば、アナログ回路やアナログ回路とディジタル回路が混在する回路において、パルス信号や急激な波形変化を伴う信号波形の過渡解析を数値演算によってシミュレーションすることができ、解析時間を短縮することができる。その結果、短時間で周波数特性を求めることができる。
【発明を実施するための最良の形態】
【００４２】
以下、本発明の実施の形態を実施例に基づき詳細に説明する。図１は、本発明において、スイッチング電源回路をモデル化して形成されたブロック構成を説明するための図である。
【００４３】
図１において、スイッチング電源回路は、前記図１２に示すＤＣ−ＤＣコンバータのディジタル制御回路１０、スイッチング回路６０、出力平滑回路２０、誤差増幅器部４０、電流検出器５０、比較器３０の各回路をモデル化した回路ブロックによって構成される。この回路ブロックは、ディジタル制御回路ブロック１、スイッチング回路ブロック６、出力平滑回路ブロック２、誤差増幅器ブロック４、電流検出回路ブロック５によって構成される。なお、図１では、ディジタル制御回路ブロック１は、電流検出回路ブロック５の出力電圧ｖcfbと誤差増幅器ブロック４の出力電圧ｖcとを引き算によって差分電圧ｖfを求める比較ブロック３を含む構成としている。
【００４４】
以下、スイッチング回路ブロック６、出力平滑回路ブロック２、誤差増幅器ブロック４、ディジタル制御回路ブロック１、および電流検出回路ブロック５の各部分のブロック化について説明する。
【００４５】
はじめに、スイッチング回路ブロック６および出力平滑回路ブロック２について説明する。図２は、本発明を適用するスイッチング電源回路のスイッチング回路６０の出力トランジスタ（Ｍｐ、Ｍｎ）、出力平滑回路２０、および負荷の部分を示した図である。
【００４６】
図１２に示される出力平滑回路２０において、インダクタＬは実際には抵抗成分ＲLを持ち、コンデンサＣにも抵抗成分Ｒcが存在する。また実際は直流出力電圧Ｖoutを抵抗で分圧する回路を誤差増幅器部４０の入力に含むため、この抵抗としてＲfが接続される。ここで、負荷は電流源Ｉoutで示している。
【００４７】
図２において、キルヒホフの法則を適用して式をたてると、
【数３】
2010092434.tif 000005
がなりたつため、
【数４】
2010092434.tif 000006
という状態方程式が得られる。
【００４８】
式（４）は、
【数５】
2010092434.tif 000007
の形で表すことができ、変数ｘの時間変化つまりｖC（ｔ），ｉL（ｔ）の時間変化を計算し、Ｖout（ｔ），ｉL（ｔ）の時間変化を求めている。ただしｘ'は時間ｔによるｘの微分を示している。ここで、ｘ、ｙ、ｕは以下で表される。
【００４９】
【数６】
2010092434.tif 000008
【００５０】
したがって変数ｙを求めることはｉL（ｔ）およびＶout（ｔ）の過渡解析に相当する。式（５）の係数にあたるＡ、Ｂ、ＣおよびＤは式（４）では回路定数から簡単に計算することができる定数であり、またｖc（ｔ），ｉL（ｔ）の時間変化波形も理想的な変化波形として考えられているので、数値演算の高速化が期待される。
【００５１】
一方、アナログ回路のシミュレーションソフト（例えば、ＳＰＩＣＥ（登録商標）と呼ばれるプログラム）を用いた場合には、シミュレーションにおける時間ステップ幅を急峻な変化時には細かく設定する必要があり、処理時間が長時間化する。
【００５２】
図３は式（４）を数値演算するためのＭＡＴＬＡＢ（登録商標）と呼ばれるプログラムの入力例である。
【００５３】
通常、プログラム入力にはこのようなシンボル表記が用いられている。式（５）の形態の状態変数方程式を計算するには外部入力Ｕと出力ｙを指定する必要がある。この場合、外部入力ＵとしてＶLX（ｔ），Ｉout（ｔ）が指定され、出力ｙにはＶout（ｔ），ｉL（ｔ）が指定されている。状態変数ｘや係数Ａ，Ｂ，Ｃ，Ｄは別途タイプ入力で指定する。
【００５４】
次に、誤差増幅器ブロック４について説明する。図４は、本発明を適用するスイッチング電源回路の誤差増幅器部４０を示した図であり、図５は誤差増幅器部４０をモデル化した回路ブロックを示した図である。ここで、Ｖoutは出力直流電圧、ｖcは誤差アンプ４６の出力であり、ｖiは誤差アンプ４６の入力である。図４の誤差増幅器回路はアナログ回路であるから、モデル化に際しては伝達関数で表現するのが適当である。
【００５５】
ここで、誤差増幅器部４０の伝達関数をｖc／Ｖoutで表すと、誤差アンプ４６の周波数特性をもこの伝達関数内に含めて記述する必要がある。誤差アンプ４６を独立して定義し、その出力にＬＰＦ（低域通過フィルタ）を接続する形で記述することによってシミュレーションの自由度を増すことができる。したがって、ここでは誤差アンプ４６の入力ｖiを重ねの理を用いて、その後に誤差アンプの電圧増幅部分（図５中の電圧増幅度（−ｋ１）で表している）と低域通過フィルタ部分（ＬＰＦ）を接続する構成によって表している。
【００５６】
入力ｖiは、直流出力電圧Ｖout、誤差アンプ４６の出力ｖc、参照電圧源４５の参照電圧Ｖref、リミタ４７を流れる電流Ｉlimに重ねの理を適用して、
【数７】
2010092434.tif 000009
という形で記述される。
【００５７】
式（７）において、伝達関数Ｈout（ｓ）は、直流出力電圧Ｖout（ｓ）以外の電圧ｖc（ｓ）、参照電圧Ｖref（ｓ）、リミタを流れる電流Ｉlim（ｓ）などの電圧源および電流源の値をゼロとして、直流出力電圧Ｖout（ｓ）のみの入力ｖi（ｓ）への影響を示す伝達関数である。ここで、各伝達関数、電圧、および電流はＳ関数で記述している。誤差アンプ４６の出力部からの伝達関数Ｈc（ｓ）、リミタ４７の伝達関数Ｈlim（ｓ）も伝達関数Ｈout（ｓ）と同様に、該当外の電圧源や電流源をゼロとして一入力のみによって求めている。このようにして求めた伝達関数Ｈout（ｓ）、伝達関数Ｈc（ｓ）、および伝達関数Ｈlim（ｓ）は、それぞれ以下の式で表される。
【００５８】
【数８】
2010092434.tif 000010
【００５９】
ただし、ＺCR＝Ｒeff＋１／ｓＣeffである。またＭlim４７は誤差アンプ４６の出力端子の電圧ｖcが大にならないようにするためのリミタである。また、抵抗４１（Ｒf1)，抵抗４２（Ｒf2)は直流出力電圧Ｖoutを分圧する抵抗であり、抵抗４８（Ｒf3)は誤差アンプ４６に入力端に接続する入力抵抗であり、電圧源４５（Ｖref）は参照電圧源である。容量４３（Ｃerr），抵抗４４（Ｒerr)は誤差アンプ４６の帰還路を構成している。
【００６０】
ここで、リミタ４７が動作した場合には、
【数９】
2010092434.tif 000011
の電流が流れる。ただしβlimはリミタ４７のトランジスタの定数（ゲート係数と呼ばれる）である。
【００６１】
図５は、式（７）を用いて誤差増幅器部４０をＭＡＴＬＡＢ（登録商標）上でモデル化した図である。誤差アンプ４６は電圧増幅部分（電圧増幅度−Ｋ１）と低域通過フィルタ部分（ＬＰＦ）の従属接続でモデル化することができる。それ以外の部分は、式（７）に従って構成した例を示している。
【００６２】
次に、ディジタル制御回路ブロック１について説明する。図６は、本発明を適用するスイッチング電源回路のディジタル制御回路１０をブロック化した回路ブロックを示した図である。図６は、図１２におけるディジタル制御回路１０をブール代数を基本として表現している。ディジタル制御回路ブロック１は、ＮＯＴ（１ｂ）、ＡＮＤ（１ｃ）、×（積）（１ｅ）、Ｓ−Ｒフリップフロップ（１ｄ）、比較結果を判定する＞＝０機能（１ａ）により構成され、その出力はＶLXで示している。これらは、実回路においてもディジタル回路の基本回路である。
【００６３】
なおＯＳＣ（発振器）はトリガ機能を担うだけであるため、ＯＳＣに限らず任意の繰り返しパルス波形を用いることができる。
【００６４】
また、図１２中の比較器３０の機能については、引き算によってモデル化することができる。なお、この引き算によるモデル化の構成は図６には示していない。図６中の比較動作１ａは、比較器３０の機能の出力ｖfを零と比較する構成部である。図７は、図６で示す各部分の信号を示している。図７（ｃ）に示す比較結果は、比較器３０の出力ｖfと零との比較結果であり、この論理値に応じてＲ−Ｓフリップフロップの出力Ｑのオンとなる期間が制御されて電圧Ｖinが出力される期間が制御され、スイッチング回路のスイッチングトランジスタの出力電圧ＶLXが制御される。
【００６５】
本構成によれば、単純な論理演算だけで出力電圧ＶLXが求まる。アナログ回路のシミュレーションに用いるＳＰＩＣＥ（登録商標）回路シミュレーションプログラムを使って計算する場合には、クロック波形の時間軸に沿って解析が行われるため、波形が急峻に変化する度に細かく時間を区切って計算をする必要があり、そのため膨大な計算時間となる。
【００６６】
これに対して、本構成によれば実際の波形と違って理想的な波形を考えて処理することができ、ＶLXに求められるのは出力トランジスタのオン・オフ制御であって簡単な機能であるため誤差も影響しにくいという効果も奏することができる。
【００６７】
次に、電流検出回路ブロック５について説明する。図８は、本発明を適用するスイッチング電源回路の電流検出回路５０を示した図である。図８には図１２のインダクタ２１に流れる電流をコピーして電流帰還ループの制御電圧として変換する回路を示している。
【００６８】
図８において、ＰＭＯＳトランジスタＭpc（５１３）は、出カトランジスタＭｐと素子サイズが異なるのみで他の条件は同一のトランジスタである。アンプＡ１（５１４）とＭpcfb（５１５）の働きでｃf1とｃf2の端子電圧が等しくなるので、電流ｉLが流れている状態ではＭpに流れる電流が素子サイズの比でＭpc（５１３）にコピーされる。Ｍpc（５１３）に流れる電流は抵抗Ｒcfb（５１６）に流れて、制御電圧ｖcfbに変換される。
【００６９】
図８ではアンプは電圧増幅度Ａ１を持つ素子として扱い、周波数特性は考慮していない。しかしながら、トランジスタＭpcは線形領域で動作するため、その出力コンダクタンスｇdsmpcは変換比に影響する。またトランジスタＭpcfbの伝達コンダクタンスｇmcfbも変換比に影響する。これらは回路構成より考えられる影響と言える。それらの影響を考慮すれば、
【数１０】
2010092434.tif 000012
で表される。
【００７０】
実際にはスイッチＳＷ１（５１１）、ＳＷ２（５１２）の影響を考慮し、また出カトランジスタＭpがオンの場合はスイッチＳＷ１（５１１）がオン、ＳＷ２（５１２）がオフとなり（図９の１の場合）、その他の場合はＳＷ１がオフ、ＳＷ２がオンする（図９の２の場合）２つの状態が考えられるので、モデル図は図９のようになる。
【００７１】
図９は、電流検出回路をモデル化した回路ブロックを示し、電流検出回路ブロック５は、トランジスタＭpの伝達コンダクタンスｇdsmp（５ａ）、スイッチＳＷの伝達コンダクタンスｇdsssw（５ｂ）、抵抗Ｒcfb（５ｄ，５ｅ）、倍率（５ｃ）等で記述することができる。
【００７２】
なお、ｇdsswはスイッチＳＷ１およびＳＷ２のコンダクタンス、ｉrefは一定の電流値、×１０は１０倍する機能を表している。また、Ｋをあらかじめ求めておけば、図９中の１の場合および２の場合について、ｖcfbを容易に計算することができる。
【００７３】
次に、図１０，図１１を用いてシミュレーションの実施例を示し、本発明による効果を示す。
【００７４】
図１０には、図１２に示す回路構成を、通常のＳＰＩＣＥ（登録商標）回路シミュレーションプログラムにより構成した回路（以下ＳＰＩＣＥ（登録商標）の回路と言う）と、本発明による数値演算用のプログラムとしてマトラボ（ＭＡＴＬＡＢ（登録商標））を用いて構成した回路（以下ＭＡＴＬＡＢ（登録商標）の回路と言う）の２種類のプログラムを用いて、図１２の回路の電源を０Ｖから３Ｖまで立ち上げた場合の直流出力電圧Ｖoutの変化をシミュレートした結果を示している。
【００７５】
図１０では時刻ｔ＝０〜２０μＳでＶinを立ち上げ、過渡解析の時間間隔を１ｎＳとして直流出力電圧Ｖoutの電圧変化を１００μＳまで観測したものを示している。なお、シミュレーション条件は、ＯＳＣ、参照電圧、および参照電流は理想状態とし、注入信号を印加しないものとし、Ｖin＝３Ｖ、Ｖout＝２Ｖ、Ｖref＝０．５Ｖ、出力電流は１００ｍＡとしている。
【００７６】
図１０では、ＳＰＩＣＥ（登録商標）の回路による結果とＭＡＴＬＡＢ（登録商標）の回路による結果は、ほぼ同じであることを示している。
【００７７】
本過渡解析において、ＳＰＩＣＥ（登録商標）の回路では約３，６００秒のＣＰＵ時間を要したが、ＭＡＴＬＡＢ（登録商標）の回路では約２５秒ほどであった。これは本発明の手法が、１００倍以上の高速化を達成できることを意味している。
【００７８】
図１１は図１２の回路の周波数特性を示す図である。周波数特性は、図１２において直流出力電圧Ｖoutと誤差増幅器部４０のマイナス（−）端子間に５０Ω程の低抵抗を接続し、この低抵抗の両端に誤差増幅器部４０のマイナス入力端子側が正となるような信号を加えて直流出力電圧Ｖoutに現れる信号成分を観測することで得られる。すなわち、等価的にループを切り離して開ループ特性を評価している。図１１（ａ）は振幅特性を示し、図１１（ｂ）は位相特性を示している。
【００７９】
図１１において、濃い実線がＭＡＴＬＡＢ（本発明の手法でプログラム）を使用したシミュレーション値を表し、薄い実線は図１２の回路を０．３５μｍＣＭＯＳプロセスを用いて実際に試作したＩＣチップより得られた実測値を表している。また、×印はＳＰＩＣＥ（登録商標）によるシミュレーション値を表している。図１２の回路を５ＭＨｚのクロックで動作させた時の特性を示している。
【００８０】
図１１によれば、振幅特性および位相特性のシミュレーション値が実測値に良好に合致していることが観察される。ただし。クロック信号の周期が５ＭＨｚであるため、特性は２．５ＭＨｚまでしか有効ではない。また、ＳＰＩＣＥ（登録商標）での解析は少なくとも１時間／１点の時間を要するため、図１１では４点しか表していない。
【００８１】
なお、図１１において、全体の周波数特性を計算する時間は、ＭＡＴＬＡＢ（登録商標）上でのシミュレーションでは約２５分を要した。
【００８２】
以上のように、図１１での比較によれば、本発明の手法を適用してＭＡＴＬＡＢ（登録商標）の数値演算プログラムを用いれば、図１２の回路の周波数特性を従来の手法と同程度の精度で高速にシミュレーションすることができる。
【００８３】
本発明の数値演算によるシミュレーション方法は、ＤＣ−ＤＣコンバータに代表されるスイッチング電源の設計に利用することができ、従来の手法に比べて同様の精度で、短時間で回路シミュレーションを行うことができるため、設計の効率を大幅に向上させることができる。
【産業上の利用可能性】
【００８４】
本発明の数値演算によるシミュレーション方法は、ＤＣ−ＤＣコンバータに代表されるスイッチング電源の設計に利用する他、パルス波形や急激な電圧変化や電流変化を持つ信号波形を処理する回路に適用することができ、アナログ回路とディジタル回路が混在する回路のシミュレーションに適用することができる。
【図面の簡単な説明】
【００８５】
【図１】本発明のスイッチング電源回路をモデル化して形成されたブロック構成を説明するための図である。
【図２】本発明を適用するスイッチング電源回路のスイッチング回路の出力トランジスタ、出力平滑回路、および負荷の部分を示す図である。
【図３】数値演算するためのプログラムの入力例である。
【図４】本発明を適用するスイッチング電源回路の誤差増幅器部を示した図である。
【図５】本発明により誤差増幅器部をモデル化した回路ブロックを示した図である。
【図６】本発明によりディジタル制御回路をモデル化した回路ブロックを示した図である。
【図７】ディジタル制御回路の回路ブロックの各部分の信号を示す図である。
【図８】本発明を適用するスイッチング電源回路の電流検出回路を示した図である。
【図９】本発明により、電流検出回路をモデル化した回路ブロックを示した図である。
【図１０】本発明による手法をＭＡＴＬＡＢ（登録商標）上に適用したプログラムと、従来のＳＰＩＣＥ（登録商標）プログラムとで同一回路構成の過渡解析を行い、実行時間を比較した図である。
【図１１】回路の周波数特性を、本発明による手法をＭＡＴＬＡＢ（登録商標）上に適用したプログラムと、従来のＳＰＩＣＥ（登録商標）プログラムと、試作した実測値とを比較した図である。
【図１２】ＬＳＩ化降圧形ＤＣ−ＤＣコンバータ（Ｖin＞Ｖout）の構成図を示す図である。
【符号の説明】
【００８６】
１    ディジタル制御回路ブロック
２    出力平滑回路ブロック
３    比較ブロック
４    誤差増幅器ブロック
５    電流検出回路ブロック
６    スイッチング回路ブロック
１０    ディジタル制御回路
２０    出力平滑回路
２１    インダクタ
３０    比較器
４０    誤差増幅器部
４５    参照電圧源
４６    誤差アンプ
４７    リミタ
５０    電流検出器
５１    スロープ補償
６０    スイッチング回路
【図１】
2010092434.tif 000013
【図２】
2010092434.tif 000014
【図３】
2010092434.tif 000015
【図４】
2010092434.tif 000016
【図５】
2010092434.tif 000017
【図６】
2010092434.tif 000018
【図７】
2010092434.tif 000019
【図８】
2010092434.tif 000020
【図９】
2010092434.tif 000021
【図１０】
2010092434.tif 000022
【図１１】
2010092434.tif 000023
【図１２】
2010092434.tif 000024