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
0 个评论
采纳的回答
Mischa Kim
2015-2-27
Dennis, you could use something like
vpa(strain,5)
13 个评论
Sambit Supriya Dash
2021-4-23
@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.
Thanks @Mischa Kim
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
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]) ]
2 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!