How can i display symbolic expressions using decimals instead of rational fractions?

113 次查看(过去 30 天)
Hello,
When I express symbolic expressions using syms with decimal values, I often get very large fractions. Since the result is not constant, I cannot use double() to evaluate the result as a decimal. sym2poly() works in most cases, but one day I will probably need to evaluate a long expression that is not a polynomial. I have also tried sym(expression,'d') but this doesn't work for non-numeric inputs. Is there some way to display symbolic expressions using decimals?
Thanks!
Example:
syms SIG_X SIG_Y SIG_Z
S_bar = [[0.78,-0.35,-0.38];[-0.35,0.78,-0.38];[-0.38,-0.38,0.92]]; %[Pa]
alpha = [-0.018,24.3,24.3]*10^-6; %[/C]
stress = [SIG_X;SIG_Y;SIG_Z];
dT = 25;
strain = S_bar*stress+alpha'*dT';
disp(strain)
Result:
(39*SIG_X)/50 - (7*SIG_Y)/20 - (19*SIG_Z)/50 - 8500259669165361/18889465931478580854784
(39*SIG_Y)/50 - (7*SIG_X)/20 - (19*SIG_Z)/50 + 5603198512389277/9223372036854775808
(23*SIG_Z)/25 - (19*SIG_Y)/50 - (19*SIG_X)/50 + 5603198512389277/9223372036854775808
Desired Result:
0.78*SIG_X - 0.35*SIG_Y - 0.38*SIG_Z - 4.5e-07
0.78*SIG_Y - 0.35*SIG_X - 0.38*SIG_Z + 6.075e-04
0.92*SIG_Z - 0.38*SIG_Y - 0.38*SIG_X + 6.075e-04

采纳的回答

Mischa Kim
Mischa Kim 2015-2-27
Dennis, you could use something like
vpa(strain,5)
  13 个评论
Sambit Supriya Dash
@Walter Roberson first, 'vpa' is simplifying the fraction in decimal form, secondly, 'vpa(strain,a)' is specifying the order of accuracy 'd'.
So, 'vpa(strain)', would also work, if you don't specify any decimal accuracy.
Walter Roberson
Walter Roberson 2021-4-24
When you vpa(strain) then that is the same as vpa(strain, digits()) -- it does the vpa at the number of default digits you have set up, rather than doing a different operation.

请先登录,再进行评论。

更多回答(1 个)

Nitesh khanna
Nitesh khanna 2020-7-27
编辑:Nitesh khanna 2020-7-27
Can we declare variables in matrix as 'real number' . if , y, then how to declare variable of matrix A as 'Real number'
A=[[sym( 'S_x', [N 3]) sym( 'S_y', [N 3]) ]

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by