数学研发论坛

 找回密码
 欢迎注册
12
返回列表 发新帖
楼主: mathematica

[提问] 如何用mathematica得到方程组的实数根?

[复制链接]
 楼主| 发表于 2011-2-26 21:52:40 | 显示全部楼层
真不知道这个符号解能干什么!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2011-2-26 22:01:25 | 显示全部楼层
(*清除所有变量*)
Clear["Global`*"];
(*求解方程,并得到10位有效数字*)
eq1=NSolve[{x^2+c^2==4^2,y^2+c^2==5^2,
            2/x==b/c,2/y==a/c,
            a+b==c},{a,b,c,x,y},10]
(*只保留数据,去掉变量符号*)
eq2={a,b,c,x,y}/.eq1
(*提取结果都是实数的结果*)
eq3=Cases[eq2,{_Real,_Real,_Real,_Real,_Real}]
(*提取都大于零的结果*)
eq4=Select[eq3,(#[[1]]>0&&#[[2]]>0&&#[[3]]>0&&#[[3]]>0&&#[[4]]>0&&#[[5]]>0)&]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-2-26 23:17:30 | 显示全部楼层
4# mathematica
  1. a=List@ToRules@Reduce[{x^2 + y^3 == 1, x^4 + y^4 == 2}, {x, y}, Reals]
复制代码
  1. N[{x, y} //. a, 40]
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-2-26 23:20:52 | 显示全部楼层
8# mathematica
  1. aa=List@ToRules@Reduce[{x^2 + c^2 == 4^2, y^2 + c^2 == 5^2, 2/x == b/c, 2/y == a/c, a + b == c}, {a, b, c, x, y}, Reals]
复制代码
  1. N[{a, b, c, x, y} //. aa, 20]
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-2-26 23:21:47 | 显示全部楼层
如果是Mathematica 8的话,还简单些,直接用Solve
  1. NSolve[{x^2 + c^2 == 4^2, y^2 + c^2 == 5^2, 2/x == b/c, 2/y == a/c, a + b == c}, {a, b, c, x, y}, Reals, WorkingPrecision -> 60]
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2011-2-27 12:03:50 | 显示全部楼层
如果是Mathematica 8的话,还简单些,直接用SolveNSolve[{x^2 + c^2 == 4^2, y^2 + c^2 == 5^2, 2/x == b/c, 2/y == a/c, a + b == c}, {a, b, c, x, y}, Reals, WorkingPrecision -> 60]
wayne 发表于 2011-2-26 23:21


我觉得mathematica有不如maple的地方,maple这个功能早就有了,
mathematica到了8才有这个功能,当然maple在群论这方面似乎也
比mathematica好。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2011-2-27 12:58:33 | 显示全部楼层
既然把我屏蔽了,我只好把我以前的代码重新发一下了:
(*清除所有变量*)
Clear["Global`*"];
(*求解方程,并得到10位有效数字*)
eq1=NSolve[{x^2+c^2==4^2,y^2+c^2==5^2,
            2/x==b/c,2/y==a/c,
            a+b==c},{a,b,c,x,y},10]
(*只保留数据,去掉变量符号*)
eq2={a,b,c,x,y}/.eq1
(*根根据虚部都等于零的结果*)
eq3=Select[eq2,Im[#]=={0,0,0,0,0}&]
(*提取结果都大于零的结果*)
eq4=Select[eq3,(#[[1]]>0&&#[[2]]>0&&#[[3]]>0&&#[[3]]>0&&#[[4]]>0&&#[[5]]>0)&]
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2020-10-22 10:59 , Processed in 0.246738 second(s), 14 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表