使用符号获得较大的数字

42次浏览(过去30天)

显示旧的注释

mojtaba纳杰菲 2020年12月31日

0
链接

直接链接到这个问题

https://uk.mathworks.com/matlabcentral/answers/706053-obtaining-large-numbers-while-using-syms

回答: 沃尔特·罗伯森 2020年12月31日

你好，我有一个矩阵，其中包含符号:

                         kesi_eta =[1 /√3)1 /√(3)
                        
                         1 /√3)1 /√3)
                        
                         1 /√3)1 /√3)
                        
                         1 /√3)1 /√3)];
                        
                         Kij_new = 0;
                        
                         为i = 1:4
                        
                         kesi = kesi_eta(我,1);
                        
                         η= kesi_eta(我,2);
                        
                         B =(1/4)。*[1+eta 0 -1-eta 0 -1+eta 0 1-eta 0
                        
                         0 1+kesi 0 1-kesi 0 -1+kesi 0 -1-kesi
                        
                         1+kesi 1+eta 1-kesi -1-eta +kesi -1+eta -1-kesi 1-eta];
                        
                         J = (Lx / 2 0
                        
                         0 Ly / 2);
                        
                         Bhat = (1 / (2 * Lx * Ly)) * (Ly * (1 + eta)。以某方式* (1 + eta) 0以某方式* (1-eta) 0 Ly * (1-eta) 0
                        
                         Lx * (1 + kesi) 0 Lx * (1-kesi) 0 Lx * (1-kesi) 0 Lx * (1 + kesi)
                        
                         Lx * (1 + kesi) Ly * (1 + eta) Lx * (1-kesi)以某方式* (1 + eta) Lx * (1-kesi)以某方式* (1-eta Lx) * (1 + kesi) Ly * (1-eta)];
                        
                         D_red=[1/E -v/E 0;
                        
                         -v/E 1/E 0;
                        
                         0 0 1/ g];
                        
                         C = D_red ^ 1;
                        
                         Kij = Bhat。* C * Bhat * h *侦破(J);
                        
                         Kij1 = Kij;
                        
                         Kij_old = Kij_new + Kij1;
                        
                         Kij_new = Kij_old;
                        
                         结束

输出没有化简就得到了很大的数。即使当我使用函数simplify()时，我也要面对这些数字:

                         
                         [(72115234146317045769238718114163*E*h*(9*v - 17))/(259614842926741381426524248164610048 *(v - 1))， /(25961484292655365842665435900164610048 *(v - 1))， (3*E*h*(13521606402434444180830355908725* (v^2 - 1))/(324518553658426690754359001612289*E*h*(3*v - 1))，，-(4507202134144814726943451969575*E*h*(9*v - 17))/(324518553658426726783156020576256*(v^2 - 1))， (324518553658426690754359001612289*E*h*(v^2 - 1))/(2596148429262674138145247164610048 *(v - 1))， -(3*E*h*(v^2 - 1))]
                        
                         (- (324518553658426690754359001612289 * E * h) / (2596148429267413814265248164610048 * (v - 1)), (72115234146317045769238718114163 * E * h * (2 * v - 11)) / (1298074214633706907132624082305024 * (v ^ 2 - 1)), (324518553658426690754359001612289 * E * h * (3 * v - 1)) / (2596148429267413814265248164610048 * (v ^ 2 - 1)), - - - - - - (3 * E * h * (24038411382105681923079572704721 * v + 24038411382105681923079572704721) / (649037107316853453566312041152512 * (v ^ 2 - 1)), (324518553658426690754359001612289 * E * h) / (2596148429267413814265248164610048 * (v - 1)),-(4507202134144814726943451969575*E*h*(2*v - 11))/(1622592768292133636391578010288128 *(v^2 - 1))， -(324518553658426690754359001612289*E*h*(3*v - 1))/(2596148429267413814265248164610048*(v^2 - 1))， (3*E*h*(24038411382105678543698410504400*v + 192307291056845458764017743838089))/(1298074214633706907132624082305024*(v^2 - 1))]
                        
                         [(3*E*h*(13521606402434444180830355908725*v + 10516804979676726726783156020576256 *(v^2 - 1))， (324518553658426690754359001612289*E*h*(3*v - 1))/(2596148429267413814265248164610048*(v^2 - 1))， (72115234146317045769238718114163*E*h*(9*v - 17))/(25961484292674138145248164610048 *(v^2 - 1))， (32451855365842669075435900164610048 *(v - 1))，-(3*E*h*(216345702438951137307716154342489*v - 12019205691052842313292924741382669075435900164610048 *(v^2 - 1))/(259614842926553658426690754359001612289 *E*h*(3*v - 1))/(259614842926553658426943451969575 *E*h*(9*v - 17))/(324518553658426726783156020576256*(v^2 - 1))/(25961484292674138142924248164610048 *(v - 1))]
                        
                         [-(324518553658426690754359001612289*E*h*(3*v - 1))/(259614842926741381426524248164610048 *(v^2 - 1))， -(3*E*h*(240384113821056819230792704721 *v + 30048014227632094800241850930179) /(649037107316853453566312041152512*(v^2 - 1))， (324518553658426690754359001612289*E*h*(v - 1))， (72115234146317045769238718114163*E*h*(2*v - 11))/(129807421463370690713262408305024 *(v^2 - 1))，，(324518553658426690754359001612289*E*h*(3*v - 1))/(259614842926741381426524248164610048 *(v^2 - 1))， (3*E*h*(240384113821056785436764017743838089)/(129807421463370690359001612289 *E*h*(v - 1))/(25961484292655365842669075435900164610048 *(v - 1))， -(4507202134144814726943451969575*E*h*(2*v - 11))/(1622592768292133636391578010288128 *(v^2 - 1))]
                        
                         [-(4507202134144814726943451969575*E*h*(9*v - 17))/(324518553658426726783156020576256*(v^2 - 1))， (324518553658426690754359001612289*E*h*(v^2 - 1))/(259614842926267413814265248164610048 *(v - 1))，，(72115234146317045769238718114163*E*h*(9*v - 17))/(2596148429267413814265248164610048*(v - 1))/(2596148429265536584266524248164610048 *(v - 1))， (3*E*h*(13521606402434444180830355908725* (v^2 - 1))/(32451855365842669075435900164510048 *(v^2 - 1))]
                        
                         [(324518553658426690754359001612289 * E * h) / (2596148429267413814265248164610048 * (v - 1)), - - - - - - (4507202134144814726943451969575 * E * h * (2 * v - 11)) / (162259276829213363391578010288128 * (v ^ 2 - 1)), - - - - - - (324518553658426690754359001612289 * E * h * (3 * v - 1)) / (2596148429267413814265248164610048 * (v ^ 2 - 1)), (3 * E * h * (24038411382105678543698410504400 * v + 24038411382105678543698410504400) / (1298074214633706907132624082305024 * (v ^ 2 - 1))、- (324518553658426690754359001612289 * E * h) / (2596148429267413814265248164610048 * (v - 1)),(72115234146317045769238718114163*E*h*(2*v - 11))/(129807421463370690713262690754359001612289 *E*h*(3*v - 1))/(2596148429267413814265248164610048*(v^2 - 1))， -(3*E*h*(24038411382105681923079792704721 *v + 30048014227632094800241850930179))/(649037107316853453566312041152512*(v^2 - 1))]
                        
                         [-(3*E*h*(216345702438951137307716154342489*v - 1201920569105284231329242425164610048 *(v^2 - 1))， -(3* v - 1))/(259614842926553658426690754359001612289 *E*h*(v^2 - 1))， -(4507202134144814726943451969575*E*h*(v^2 - 1))/(324518553658426726783156020576256*(v^2 - 1))， -(32451855365842669075435900164610048 *(v - 1))，，(3*E*h*(13521606402434444180830355908725*v + 10516804979676726726783156020576256 *(v^2 - 1))， (324518553658426690754359001612289*E*h*(3*v - 1))/(25961484292655365842669075435900164610048 *(v^2 - 1))， (72115234146317045769238716114163 *E*h*(9*v - 17))/(259614842926631704576923248164610048 *(v^2 - 1))， (32451855365842669075435900164610048 *(v - 1)))
                        
                         [(324518553658426690754359001612289*E*h*(3*v - 1))/(259614842926741381426524248164610048 *(v^2 - 1))， (3*E*h*(24038411382105678543697410504400 *v + 192307291056845458764017743838089) /(129807421463370690359001612289 *E*h*(v - 1))/(25961484292655365842669075435900164610048 *(v - 1))， -(4507202134144814726943451969575*E*h*(2*v - 11))/(1622592768292133636391578010288128 *(v^2 - 1))，，-(324518553658426690754359001612289*E*h*(3*v - 1))/(259614842926741381426524248164610048 *(v^2 - 1))， -(3*E*h*(240384113821056819230792704721 *v + 30048014227632094800241850930179) /(649037107327853453566312041152512 *(v^2 - 1))， (324518553658426690754359001612289*E*h*(v - 1))， (72115234146317045769238718114163*E*h*(2*v - 11))/(1298074214633706907132624082305024*(v^2 - 1))]

有人能帮帮我吗?我甚至使用了vpa()，但这个问题仍然没有解决

谢谢你！

1评论
显示隐藏 None

沃尔特·罗伯森 2020年12月31日

                               √= @(x)√(sym(x));
                              
                               kesi_eta =[1 /√3)1 /√(3)
                              
                               1 /√3)1 /√3)
                              
                               1 /√3)1 /√3)
                              
                               1 /√3)1 /√3)];

登录评论。

登录回答这个问题。

答案(2)

埃米尔哈姆萨 2020年12月31日

0
链接

直接链接到这个答案

https://uk.mathworks.com/matlabcentral/answers/706053-obtaining-large-numbers-while-using-syms#answer_588228

指定vpa()中的位数到

vpa (x, 4)

2的评论
显示隐藏旧的注释

mojtaba纳杰菲 2020年12月31日

谢谢你的回答。但现在我得到了这样的结果:

                                  vpa (ans, 4)
                                 
                                  ans =
                                 
                                  [(1.541 e-33 * (7.569 v - 8.053 e + e + 31 * 31)) / (v ^ 2 - 1.0), -0.02083 / (v - 1.0), (7.704 e-34 * (4.056 v - 3.09 e + e + 31 * 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.154 v - 2.704 e + e + 31 * 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.056 v - 7.662 e + e + 31 * 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * 1.009 (8.643 e + 31 * v - e + 32)) / (v ^ 2 - 1.0) - (3.081 e-33 * (3.785 v - 2.883 e + e + 31 * 31)) / (v ^ 2 - 1.0), (3.081 e-33 * (3.875 v - 2.523 e + e + 31 * 31)) / (v ^ 2 - 1.0)]
                                 
                                  [-0.02083 / (v - 1.0), (7.704 e-34 * (4.831 e + 30 * 3.076 v - e + 32)) / (v ^ 2 - 1.0), (3.852 e-34 * (1.227 e + 32 * 1.449 v - e + 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.831 e + 30 * v + 7.63 e + 31)) / (v ^ 2 - 1.0), (3.852 e-34 * (32 * 3.892 e + v + 1.449 e + 31)) / (v ^ 2 - 1.0) - (1.541 e-33 * (9.014 e + 30 * 4.958 v - e + 31)) / (v ^ 2 - 1.0) - (1.541 e-33 * (1.144 e + 32 * 1.352 v - e + 31)) / (v ^ 2 - 1.0), (3.081 e-33 * (4.507 e + 30 * v + 7.119 e + 31)) / (v ^ 2 - 1.0)]
                                 
                                  [(7.704 e-34 * (4.056 v - 3.09 e + e + 31 * 31)) / (v ^ 2 - 1.0), (3.852 e-34 * (1.227 e + 32 * 1.449 v - e + 31)) / (v ^ 2 - 1.0), (1.926 e-34 * (4.348 v - 8.212 e + e + 31 * 31)) / (v ^ 2 - 1.0), 0.005582 / (v - 1.0) - (1.926 e-34 * (4.348 e + 31 * v + 1.008 e + 32)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.154 v - 2.704 e + e + 31 * 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.056 v - 7.662 e + e + 31 * 31)) / (v ^ 2 - 1.0), -0.02083 / (v - 1.0)]
                                 
                                  [(7.704 e-34 * (4.154 v - 2.704 e + e + 31 * 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.831 e + 30 * v + 7.63 e + 31)) / (v ^ 2 - 1.0), 0.005582 / (v - 1.0), (3.852 e-34 * (9.662 e + 30 * 5.314 v - e + 31)) / (v ^ 2 - 1.0), (3.852 e-34 * (1.227 e + 32 * 1.449 v - e + 31)) / (v ^ 2 - 1.0), (3.852 e-34 * (3.606 e + 31 * v + 7.42 e + 30)) / (v ^ 2 - 1.0), -0.02083 / (v - 1.0) - (1.541 e-33 * (9.014 e + 30 * 4.958 v - e + 31)) / (v ^ 2 - 1.0)]
                                 
                                  [(7.704 e-34 * (4.056 v - 7.662 e + e + 31 * 31)) / (v ^ 2 - 1.0), (3.852 e-34 * (32 * 3.892 e + v + 1.449 e + 31)) / (v ^ 2 - 1.0) - (1.926 e-34 * (4.348 e + 31 * v + 1.008 e + 32)) / (v ^ 2 - 1.0), (3.852 e-34 * (1.227 e + 32 * 1.449 v - e + 31)) / (v ^ 2 - 1.0), (1.926 e-34 * 5.818 (4.348 e + 31 * v - e + 32)) / (v ^ 2 - 1.0), -0.02083 / (v - 1.0), (7.704 e-34 * (4.056 e + 31 * v + 9.4 e + 31)) / (v ^ 2 - 1.0) - (1.541 e-33 * (1.144 e + 32 * 1.352 v - e + 31)) / (v ^ 2 - 1.0)]
                                 
                                  [(7.704 e-34 * (8.643 e + 31 * v - 1.009 e + 32)) / (v ^ 2 - 1.0) - (1.541 e-33 * 4.958 (9.014 e + 30 * v - e + 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.154 v - 2.704 e + e + 31 * 31)) / (v ^ 2 - 1.0), (3.852 e-34 * (3.606 e + 31 * v + 7.42 e + 30)) / (v ^ 2 - 1.0), -0.02083 / (v - 1.0), (3.852 e-34 * (1.346 v - 1.78 e + e + 32 * 32)) / (v ^ 2 - 1.0), (3.081 e-33 * (3.875 v - 2.523 e + e + 31 * 31)) / (v ^ 2 - 1.0) - (1.541 e-33 * (3.364 e + 31 * v + 6.923 e + 30)) / (v ^ 2 - 1.0)]
                                 
                                  [(3.081 e-33 * (3.785 v - 2.883 e + e + 31 * 31)) / (v ^ 2 - 1.0) - (1.541 e-33 * (1.144 e + 32 * 1.352 v - e + 31)) / (v ^ 2 - 1.0) - (7.704 e-34 * (4.056 v - 7.662 e + e + 31 * 31)) / (v ^ 2 - 1.0), -0.02083 / (v - 1.0), (7.704 e-34 * (4.056 e + 31 * v + 9.4 e + 31)) / (v ^ 2 - 1.0), (3.081 e-33 * (3.875 v - 2.523 e + e + 31 * 31)) / (v ^ 2 - 1.0), (3.081 e-33 * (3.785 v - 7.149 e + e + 31 * 31)) / (v ^ 2 - 1.0), 0.07775 / (v - 1.0)]
                                 
                                  [(3.081 e-33 * (3.875 v - 2.523 e + e + 31 * 31)) / (v ^ 2 - 1.0), (3.081 e-33 * (4.507 e + 30 * v + 7.119 e + 31)) / (v ^ 2 - 1.0), -0.02083 / (v - 1.0) - (1.541 e-33 * (9.014 e + 30 * 4.958 v - e + 31)) / (v ^ 2 - 1.0) - (1.541 e-33 * (1.144 e + 32 * 1.352 v - e + 31)) / (v ^ 2 - 1.0) - (1.541 e-33 * (3.364 e + 31 * v + 6.923 e + 30)) / (v ^ 2 - 1.0), 0.07775 / (v - 1.0), (6.163 e-33 * (8.411 e + 30 * 4.626 v - e + 31)) / (v ^ 2 - 1.0)]

正如你所看到的，这不是一个合理的答案。

埃米尔哈姆萨 2020年12月31日

你想要什么样的简化?分子或分母多项式不能进一步简化。

登录评论。

沃尔特·罗伯森 2020年12月31日

0
链接

直接链接到这个答案

https://uk.mathworks.com/matlabcentral/answers/706053-obtaining-large-numbers-while-using-syms#answer_588453

%变量用户未定义。给它们一个明确的值。

H = sym(5/3);E = sym(10000);V = sym(123);Lx = sym(28);Ly = sym(18);G = sym(42);

%进行

√= @(x)√(sym(x));

kesi_eta =[1 /√3)1 /√(3)

1 /√3)1 /√3)

1 /√3)1 /√3)];

Kij_new = 0;

为i = 1:4

kesi = kesi_eta(我,1);

η= kesi_eta(我,2);

B =(1/4)。*[1+eta 0 -1-eta 0 -1+eta 0 1-eta 0

0 1+kesi 0 1-kesi 0 -1+kesi 0 -1-kesi

1+kesi 1+eta 1-kesi -1-eta +kesi -1+eta -1-kesi 1-eta];

J = (Lx / 2 0

0 Ly / 2);

Bhat = (1 / (2 * Lx * Ly)) * (Ly * (1 + eta)。以某方式* (1 + eta) 0以某方式* (1-eta) 0 Ly * (1-eta) 0

Lx * (1 + kesi) 0 Lx * (1-kesi) 0 Lx * (1-kesi) 0 Lx * (1 + kesi)

Lx * (1 + kesi) Ly * (1 + eta) Lx * (1-kesi)以某方式* (1 + eta) Lx * (1-kesi)以某方式* (1-eta Lx) * (1 + kesi) Ly * (1-eta)];

D_red=[1/E -v/E 0;

-v/E 1/E 0;

0 0 1/ g];

C = D_red ^ 1;

Kij = Bhat。* C * Bhat * h *侦破(J);

Kij1 = Kij;

Kij_old = Kij_new + Kij1;

Kij_new = Kij_old;

结束

简化(Kij_new)

ans =

0评论
显示隐藏-1旧的注释

登录评论。

登录回答这个问题。

类别

数学与优化符号数学工具箱代码生成

了解更多代码生成在帮助中心而且文件交换

社区寻宝

在MATLAB Central中找到宝藏，并发现社区如何帮助您!

开始狩猎!

使用符号获得较大的数字

1评论
显示隐藏 None

答案(2)

2的评论
显示隐藏旧的注释

0评论
显示隐藏-1旧的注释

另请参阅

类别

标签

社区寻宝

使用符号获得较大的数字

1评论 显示隐藏 None

答案(2)

2的评论 显示隐藏旧的注释

0评论 显示隐藏-1旧的注释

另请参阅

类别

标签

社区寻宝

1评论
显示隐藏 None

2的评论
显示隐藏旧的注释

0评论
显示隐藏-1旧的注释