Strange output from symbolic calculation

Question

David Winthrop 2022-11-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1840988-strange-output-from-symbolic-calculation

评论： Steven Lord 2022-11-1

Can someone explain this output (in particular, the numbers)

syms f0 T t0 Omegan t n
subs(2*f0*pi^2/(Omegan*T*(- Omegan^2*t0^2 + 4*pi^2)),Omegan,2*n*pi/T)
ans =
-(2778046668940015*f0)/(281474976710656*n*pi*((4*n^2*t0^2*pi^2)/T^2 - 2778046668940015/70368744177664))

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

David Winthrop 2022-11-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1840988-strange-output-from-symbolic-calculation#answer_1088893

在 MATLAB Online 中打开

I figured out what was going on. It seems the sym interpreter was giving a rational approximation for pi only some of the time. Very strange. For future reference, this is how I got what is expected:

syms f0 T t0 Omegan t n
p = sym(pi); % Key step here.
simplify(subs(2*f0*p^2/(Omegan*T*(- Omegan^2*t0^2 + 4*p^2)),Omegan,2*n*p/T))
ans =
-(T^2*f0)/(4*n*pi*(- T^2 + n^2*t0^2))

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Gautam Chettiar 2022-11-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1840988-strange-output-from-symbolic-calculation#answer_1088873

编辑：Gautam Chettiar 2022-11-1

在 MATLAB Online 中打开

% The numbers come because of your substitution of Omegan to 2*n*pi/T, with

% some of the variables cancelling out, and most of the numbers coming from

% 2 and pi

syms f0 T t0 Omegan t n

ans = subs(2*f0*pi^2/(Omegan*T*(- Omegan^2*t0^2 + 4*pi^2)),Omegan,2*n*pi/T)

ans =

vpa(simplify(ans), 3)

ans =

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 3

James Tursa 2022-11-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1840988-strange-output-from-symbolic-calculation#answer_1088908

编辑：James Tursa 2022-11-1

在 MATLAB Online 中打开

This isn't strange at all. You are asking the Symbolic Engine to examine all floating point expressions involving pi, which isn't even represented exactly in IEEE floating point, and extract the pi part. There are an infinite number of such expressions, of course, so you can't expect the Symbolic Engine to get them all. As such, it only gets the simple ones. E.g.,

sym(pi)

ans =

π

sym(4*pi)

ans =

sym(pi/3)

ans =

sym(pi^2)

ans =

sym(pi)^2

ans =

sym(exp(1))

ans =

exp(sym(1))

ans =

e

sym(exp(1)/3)

ans =

exp(sym(1))/3

ans =

sym(pi*exp(1))

ans =

sym(pi)*exp(sym(1))

ans =

The solution, as you have already discovered, is to insert symbolic constants up front rather than passing results of floating point expressions to the Symbolic Engine and expecting it to divine the source in all cases. Even the simple expressions have a limit for when the Symbolic Engine stops looking. E.g.,

sym(pi/100000)

ans =

sym(pi/100001)

ans =

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Steven Lord 2022-11-1

See the description of the 'r' value for the flag input argument on the sym function's documentation page for a list of the expressions sym "recognizes". Most of the examples you posted that convert pi to π fall into the p*pi/q (for modest-sized integers p and q) case.

请先登录，再进行评论。