How can i display symbolic expressions using decimals instead of rational fractions?
178 ビュー (過去 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 件のコメント
参考
カテゴリ
Help Center および File Exchange で Symbolic Math Toolbox についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!