double

シンボリック値を MATLAB の倍精度に変換

構文

d = double(s)

説明

d = double(s) は、シンボリック値 s を倍精度に変換します。シンボリック値の倍精度への変換は、MATLAB^® 関数でシンボリック値を受け取らない場合に役立ちます。シンボリック数と倍精度数の違いについては、数値演算またはシンボリック演算の選択を参照してください。

例

シンボリック数を倍精度に変換

double を使用してシンボリック数を倍精度に変換します。シンボリック数は厳密ですが、倍精度数には丸め誤差があります。

$\pi$ と $\frac{1}{3}$ をシンボリック型から倍精度に変換します。

symN = sym([pi 1/3])

symN = 
$$\left(\begin{array}{cc}\pi & \frac{1}{3}\end{array}\right)$$

doubleN = double(symN)

doubleN = 1×2

    3.1416    0.3333

丸め誤差については、丸め誤差の認識と回避を参照してください。

可変精度の倍精度への変換

vpaで作成される可変精度の数値はシンボリック値です。MATLAB 関数でシンボリック値を受け取らない場合、double を使用して可変精度を倍精度に変換します。

$\pi$ と $\frac{1}{3}$ を可変精度から倍精度に変換します。

vpaN = vpa([pi 1/3])

vpaN = $$\left(\begin{array}{cc}3.1415926535897932384626433832795& 0.33333333333333333333333333333333\end{array}\right)$$

doubleN = double(vpaN)

doubleN = 1×2

    3.1416    0.3333

シンボリック行列の倍精度行列への変換

行列 symM のシンボリック数を、double を使用して倍精度数に変換します。

a = sym(sqrt(2));
b = sym(2/3);
symM = [a b; a*b b/a]

symM = 
$$\left(\begin{array}{cc}\sqrt{2}& \frac{2}{3}\\ \frac{2\hspace{0.17em}\sqrt{2}}{3}& \frac{\sqrt{2}}{3}\end{array}\right)$$

doubleM = double(symM)

doubleM = 2×2

    1.4142    0.6667
    0.9428    0.4714

高精度な変換

内部桁落ちや丸め誤差の影響を受けるシンボリック式を変換する場合、数値変換の前に、digitsを使用して作業精度を引き上げます。

数値的に不安定な式 Y を double を使用して変換します。次に、digits を使用して精度を 100 桁に引き上げ、再度 Y を変換します。低い精度の変換が正確ではないのに対し、こうした高い精度の変換は正確です。

Y = ((exp(sym(200)) + 1)/(exp(sym(200)) - 1)) - 1;
lowPrecisionY = double(Y)

lowPrecisionY = 
0

digitsOld = digits(100);
highPrecisionY = double(Y) 

highPrecisionY = 
2.7678e-87

digits を使用して古い精度に戻し、計算を続けます。

digits(digitsOld)

シンボリックな結果から倍精度への変換

solve を使用して三角関数方程式 $\sin \left(2\mathit{x}\right)+\cos \left(\mathit{x}\right)=0$ を解きます。完全解、解で使用されているパラメーター、およびそれらのパラメーターの条件を返すように、ReturnConditions オプションを true に設定します。

syms x
eqn = sin(2*x) + cos(x) == 0;
[solx,params,conds] = solve(eqn,x,ReturnConditions=true)

solx = 
$$\left(\begin{array}{c}\frac{\pi }{2}+\pi \hspace{0.17em}k\\ 2\hspace{0.17em}\pi \hspace{0.17em}k-\frac{\pi }{6}\\ \frac{7\hspace{0.17em}\pi }{6}+2\hspace{0.17em}\pi \hspace{0.17em}k\end{array}\right)$$

params = $$k$$

conds = 
$$\left(\begin{array}{c}k\in \mathbb{Z}\\ k\in \mathbb{Z}\\ k\in \mathbb{Z}\end{array}\right)$$

ソルバーは、このパラメーターの変数 $\mathit{k}$ を MATLAB® ワークスペースに作成しません。この変数を作成します。subs を使用して $\mathit{k}=2$ の解を求めます。

syms k
sols_k2 = subs(solx,k,2)

sols_k2 = 
$$\left(\begin{array}{c}\frac{5\hspace{0.17em}\pi }{2}\\ \frac{23\hspace{0.17em}\pi }{6}\\ \frac{31\hspace{0.17em}\pi }{6}\end{array}\right)$$

正確なシンボリック数が解として得られます。これらの数値を倍精度の数値に変換します。

doublesols_k2 = double(sols_k2)

doublesols_k2 = 3×1

    7.8540
   12.0428
   16.2316

数値のシンボリック行列変数から倍精度への変換

${\mathit{A}}^2-2\mathit{A}+{\mathit{I}}_2$ を表すシンボリック式 S を作成します。ここで、$\mathit{A}$ は 2 行 2 列のシンボリック行列変数です。

syms A 2 matrix
S = A*A - 2*A + eye(2)

S = $${\mathrm{I}}_2-2\hspace{0.17em}A+A^2$$

$\mathit{A}$ にシンボリック数 $\left[\begin{array}{cc}\cos \left(\frac{\pi }{5}\right)& \sin \left(\frac{\pi }{4}\right)\\ -1& 0\end{array}\right]$ を代入します。

Aval = [cos(sym(pi)/5) sin(pi/4); -1 0];
symS = subs(S,A,Aval)

symS = 
$$\begin{array}{l}-2\hspace{0.17em}{\Sigma }_1+{{\Sigma }_1}^2+{\mathrm{I}}_2\\ \\ \mathrm{where}\\ \\ \mathrm{ }{\Sigma }_1=\left(\begin{array}{cc}\frac{\sqrt{5}}{4}+\frac{1}{4}& \frac{\sqrt{2}}{2}\\ -1& 0\end{array}\right)\end{array}$$

結果を倍精度の行列に変換します。

doubleS = double(symS)

doubleS = 2×2

   -0.6706   -0.8422
    1.1910    0.2929

入力引数

s — シンボリック入力

シンボリック数 | シンボリック数の配列 | 数値のシンボリック行列変数

シンボリック入力。シンボリック数、シンボリック数の配列、または数値のシンボリック行列変数として指定します。入力に数値のみが含まれる場合 (シンボリック数が含まれない場合)、MATLAB 関数 double が呼び出されます。

データ型: sym | symmatrix

バージョン履歴

R2006a より前に導入

R2022b: 数値のシンボリック行列変数から倍精度への変換

数値のシンボリック行列変数 (symmatrix 型) を double 型の倍精度行列に変換できます。例については、数値のシンボリック行列変数から倍精度への変換を参照してください。

参考

sym | vpa | subs

トピック

• 数値演算またはシンボリック演算の選択
• 数値計算精度の引き上げ
• 丸め誤差の認識と回避
• 精度の引き下げによる速度の向上
• Change Output Format of Symbolic and Variable-Precision Arithmetic