下式は、状態空間モデルです。

A/B/C/D 行列は、ファイルとして与えます。 4次のrungekuttaで計算しています。呼び出し手順は、

$runge_kutta("File Name",mode,Input_RealArray,Output_RealArray,[X])；

• ”File Name"は、A/B/C/Dの係数行列のテキストファイルを指定します。.
• mode は、initial statementでは、０　それ以外(always) は、１にします。
• 入力vector　u は、Input_RealArrayで与えます.
• 出力結果は、Output_RealArray に入ります。.
• 内部状態Xもオプションで設定可能です。

例として "bessel.vtakprj"　をご参照ください。NRZからベッセルフィルタを通した波形が得られています。.

 

Analog　PLLのシミュレーション

Analog PLLにおいては、下図のループフィルタが安定の鍵になります。

この伝達関数は、下式によります。

Veritakのアナログシミュレータは伝達関数を状態空間表現にする必要があります。手計算でもできますが、MatlabやScilabのtf2ss　コマンドで簡単に変換できます。たとえば、次のようなスクリプトを書いて変換します。
これで、Veritakの行列読み込みファイル"loop_filter.txt"ができました。"pll_test.vtakprj" をLoadしてみます。
PLLがロックする様子が分かります。

TIPS

古典制御理論によるPLLの設計例が、下記にあります。
http://www.sugawara-systems.com/opencores/pll/pll.htm

　途中で内部状態を切り替えるには、適宜、内部状態行列Xをリードライトしてください。ソースが必要な方はサポートまでご連絡ください。