1

我的问题与这个问题相似,但我认为它更笼统。
我使用 Matlab 的符号数学工具箱来求解一个方程:

MAZ = 0.5;
MAU = 1.0;
XI = 1.0;
ALPHA = 2.0;

DRG = 0.5;
SRG = 1.0;

PHI = 1 / (2 * MAU);

syms L;

f = 1 - DRG - sqrt(1 + (ALPHA * XI - L / (2 * XI * PHI) ) ^ 2 ) / ...
    sqrt(1 + (ALPHA * XI) ^ 2) + L / (4 * PHI * SRG * sqrt(1 + (ALPHA * XI)^2));

a = solve(f,L,'Real',true);

答案是:

a =
5^(1/2)/3 + (10*((4*5^(1/2))/25 + 6/25)^(1/2))/3 + 8/3
5^(1/2)/3 - (10*((4*5^(1/2))/25 + 6/25)^(1/2))/3 + 8/3  

如何自动将这些a不包含任何符号表达式的表达式转换为浮点数,以便稍后在我的代码中访问它们?

4

2 回答 2

3

Matlab 的符号数学工具箱包含一个名为double.
此函数将函数的结果solve- 即 a sym- 转换为双精度数:

double(a)

ans =

      5.98921078320145
     0.834834535131742
于 2013-04-24T19:58:42.997 回答
1

虽然double是正确的,但它仅限于 64 位浮点的精度。vpa如果需要更高的精度,请使用:

>> vpa(a)

ans =

  5.9892107832014511063435699584181
 0.83483453513174202459587915406938

>> digits(99)

>> vpa(a)

ans =

  5.98921078320145110634356995841813862213621375395128614574627036653958858547362556353272837962692249
 0.834834535131742024595879154069378868157531819123064337100994463734092031618244369410214559292265698
于 2015-08-03T16:50:36.260 回答