RootFinding[Parametric] - パラメトリック多項式システムの計算用パッケージ
|
使い方
|
|
RootFinding[Parametric][command](arguments)
command(arguments)
|
|
説明
|
|
•
|
RootFinding[Parametric] パッケージは多項式および不等式のパラメトリックシステム向けに設計されています。このパッケージはパラメトリック多項式システム
|
|
を構造解析し、実数解を持つパラメータ値やシステムにおける解の数を得るためのプログラムが含まれています。式で、f および g は不定元 (変数) およびパラメータと呼ばれる 2 種類の変数の多項式を表します。
|
•
|
一方で補集合 は一般にセルと呼ばれる関連した開集合の有限和で、入力システムの実数解の数が有限かつ同一セルではパラメータ値を変更してもその数は一定であるという重要な特性を持っています。これにより、システムの実数解の数を基にしたパラメータ空間の完全な分類ができます。
|
•
|
このパッケージに含まれているコマンドは、円筒代数分解 (Cylindrical Algebraic Decomposition; CAD) の各結果セルにおいてシステムの解の数が一定であるという主要特性を保つようにして、CAD でセル をさらに分割します。CAD 結果の最大次元の各オープンセル (open cell) に対しては、完全にセル内部に位置するサンプル点が計算されます。
|
•
|
入力システムは以下の要件を満たすものとします。
|
–
|
ほとんどの複素パラメータ値に対して最低でも 1 つ、最大でも有限数の複素解が存在する (すなわち、システムは一般的に可解で、一般的に 0 次元である)。
|
–
|
ほとんどの複素パラメータ値に対して解の重複度は 1 を上回らない (すなわち、システムは一般的に根基である)。具体的には、入力方程式に平方因子は含まれていない。
|
|
|
RootFinding[Parametric] パッケージのコマンドへのアクセス
|
|
•
|
RootFinding[Parametric] パッケージの各コマンドは long form (長形式) または short form (短形式) のいずれかのフォーマットでコマンド名を指定することによりアクセスできます。RootFinding[Parametric] パッケージはモジュールとして実装されているため、RootFinding[Parametric]:-command の形式を使用してもパッケージのコマンドにアクセスできます。
|
|
|
RootFinding[Parametric] パッケージのコマンド一覧
|
|
|
|
解レコードのデータ構造
|
|
•
|
RootFinding[Parametric] パッケージのコマンドの多くはパラメトリック多項式システムの構造情報を含む解レコードを出力するか、解レコードを入力として受け付けます。これは以下のフィールドを持つ record (レコード) です。
|
–
|
Equations: システムの方程式を表す多項式 のリスト
|
–
|
Inequalities: 形式のシステムの不等式を表す多項式 のリスト
|
–
|
Parameters: 名前のリスト; パラメータ
|
–
|
DiscriminantVariety: パラメータの多項式のリストリスト; discriminant variety は各多項式の内部リストの判別式の和になります。
|
–
|
ProjectionPolynomials: CAD (円筒代数分解) の意味での射影多項式を表すパラメータの多項式のリストリスト
|
–
|
SamplePoints: parameter=rational number の形式の方程式のリストリスト。リストリストは CAD の各オープンセルに対して 1 つ存在し、各リストは完全にセルの内部に位置するサンプル点を表します。
|
•
|
解レコード のフィールドは、たとえば m:-Equations, m:-Inequalities のように、選択演算子 (:-) を使用することでアクセスできます。
|
•
|
解レコードのフィールドは変更可能ですが、これは推奨されません。
|
|
|
例
|
|
この例では 1 変数のパラメトリック多項式 について、パラメータ a および b の値による実数解の数を調べます。この場合は既知の判別式に対応する discriminant variery を計算します。
>
|
|
| (6.1) |
>
|
|
| (6.2) |
>
|
|
| (6.3) |
>
|
|
| (6.4) |
2 次元のパラメータ空間を最大次元のオープンセルに分解し、各セルの内部に位置するサンプル点を含めてプロットします。
>
|
|
| (6.5) |
>
|
|
| (6.6) |
>
|
|
合計 4 つのセルが存在することにご注目ください。各セルに対して、システムの実数解の数を計算します。
>
|
|
| (6.7) |
>
|
|
| (6.8) |
パラメータ値を黄色またはピンクのセルから選択した場合、システムの実数解数は 1 で、青または緑のセルから選択した場合、システムの実数解数は 3 になります。該当するサンプル点でシステムの解を求めることでセル 1 および 2 に対する結論が正しいことを検証します。
>
|
|
| (6.9) |
>
|
|
| (6.10) |
2 つ目のセルの幾何学的記述を求めます。
>
|
|
| (6.11) |
結果は次のように解釈します;パラメータ空間における点 は
の場合に限って 2 つ目のセルに位置します。
点 はこの要件を満たしています。確かに 2 つ目のセルに位置していることを確認します。
>
|
|
| (6.12) |
|
|
参考文献
|
|
Lazard, D., and Rouillier, F. "Solving parametric polynomial systems." Journal of Symbolic Computation. (2007).
|
|
Download Help Document
Was this information helpful?