24小时热门版块排行榜    

查看: 991  |  回复: 3

wangzgenyu

新虫 (小有名气)

[求助] mathematica不能化简得问题 已有1人参与

程序编了有一段时间了,只是一个简单的二元二次方程组求解问题,但就是不能很好的化简,现在代码拷贝在下面,希望有人能指点一些。另外将代码的规则写出以后,还带有TextCell[""],不知道是为什么。
Clear["Global`*"]
Di = 1.98;
tc = 0.01;
Ri = 1.99;
hs = 0.01;
ts = 0.01;
p = 1.6 10^6;
nu = 0.3;
r = 0.5*(Di + tc);
R = Ri + 0.5*hs;
nfy = ArcSin[r/R];
kc = Power[3 (1 - nu^2), (4)^-1]/Sqrt[r*tc];
nlam = Power[3 (1 - nu^2), (4)^-1] Sqrt[R/hs];
k21 = 1 - (1 - 2 nu)/(2 Power[3 (1 - nu^2), (4)^-1]) Sqrt[hs/R]
     Cot[nfy];
k22 = 1 - (1 + 2 nu)/(2 Power[3 (1 - nu^2), (4)^-1]) Sqrt[hs/R]
     Cot[nfy];
k2d = 1 -
   nu/Power[3 (1 - nu^2), (4)^-1] Sqrt[hs/R] Cot[nfy] + (1 - 2 nu)^2/(
    8 Sqrt[3 (1 - nu^2)]) hs/R (Cot[nfy])^2;
k23 = 1 + (2 - nu)/(2 nu Power[3 (1 - nu^2), (4)^-1]) hs/R Cot[nfy];

a11 = Sqrt[r/tc]*(1 + Sqrt[Sin[nfy]]*k2d);
a12 = (1 - tc^2/(k21*hs^2)) ;
a21 = (1 - tc^2/(k21*hs^2));
a22 = 2 ((3 (1 - nu^2))^0.75)*Sqrt[tc/
   r] (1 + 1/(k21 Sqrt[Sin[nfy]]) (tc/hs)^2.5);
a10 = 1/4 ((1 - nu)/Sin[nfy] tc/hs -
     2 Power[3 (1 - nu^2), (4)^-1] Sqrt[R/hs] tc/hs k2d*
      Cos[nfy] - (2 - nu));
a20 = (p*r)/(2.0*k21)*Cot[nfy];
Solve[{a11 x + a12 y == a10, a21 x + a22 y == a20}, {x, y}]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangzgenyu

新虫 (小有名气)

这个方法我把具体表达式写出来,跳过solve命令:
azz = (a20 a11 - a10 a21)/(a11 a22 - a12 a21),试着求解,结果还是不化简,是这样的;
2.37196*10^7/(174.572 - 0.000380934 TextCell[""])
2楼2016-05-18 16:18:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangzgenyu

新虫 (小有名气)

终于解决了,自己粗心引入了一个不知道哪里来的不知道是什么的一个符号。而且还是透明的。唉!
3楼2016-05-18 16:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kill_speed

铁虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
这代码写的挺流畅,能帮我用这软件解决两个问题么?

http://muchong.com/bbs/viewthread.php?tid=10194430

http://muchong.com/bbs/viewthread.php?tid=10209670
天,有多高?不上去怎会知道。海,有多深?不下去怎会知晓。
4楼2016-05-19 12:36:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangzgenyu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-02 00:29 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 00:08 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 23/1150 2026-11-01 23:49 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 23:45 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 21:52 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-01 12:51 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
信息提示
请填处理意见