找回密码
 欢迎注册
楼主: chyanog

[提问] 如何求一个一元八次方程的符号解

[复制链接]
发表于 2013-1-20 20:35:17 | 显示全部楼层
60# mathe 漂亮!! 根据这个提示我构造了一个
  1. MinimalPolynomial[AlgebraicNumber[Root[#^8 - 2 &, 1], {0, 1, 1, 1, 1, 1}], x]
复制代码
$x^8-24 x^6-96 x^5-132 x^4-208 x^3-248x^2-112 x-62=0$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-1-20 21:04:57 | 显示全部楼层
这个恐怕是最简洁的非平凡的有根式解的8次方程了。 $2 - 16 x - 40 x^2 - 32 x^3 - 4 x^4 + x^8=0$ 分解因式就是 $(x - a - a^2) (x + a - a^2) (x^2 + a^5x + a^2 + a^4 - a^7) (x^2 - a^5x + a^2 + a^4 + a^7) (x^2 + 2 a^2x+ a^2 + a^4) $ a表示2的开8次方,1.090507733
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-1-21 10:35:18 | 显示全部楼层
谁来试验一下这个东东: http://www.gap-system.org/Releases/index.html
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-1-21 11:30:38 | 显示全部楼层
4# wayne Mathematica9装在办公室的电脑上,周末在家,没得用,只好用本本上的Maxima。我倒是想到过20#的分解方法,可惜在Maxima用得不如Mathematica9熟练,不知道能不能实现。 将1#那个方程略作处理,可以得到很漂亮的解。 由于$x^8$的系数不是1,所以方程的根不是代数整数,但是常数项为1,所以1/x是代数整数,故我们作代换y=1/x,将方程化为
  1. 256 + 1536 y - 4608 y^2 + 768 y^3 + 1040 y^4 - 96 y^5 - 72 y^6 + y^8==0
复制代码
然后来分解因式:
  1. In[7]:= Factor[ 256 + 1536 y - 4608 y^2 + 768 y^3 + 1040 y^4 - 96 y^5 - 72 y^6 + y^8, Extension -> {Sqrt[2], Sqrt[3], Sqrt[7]}]
复制代码
捕获.JPG

评分

参与人数 1鲜花 +12 收起 理由
wayne + 12 太神奇了,这都能发现!!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-1-21 12:11:23 | 显示全部楼层
64# hujunhua ,Mathematica自动化简的没这么简洁
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-1-21 12:23:32 | 显示全部楼层
我突然想起,在Mathematica中,根式分母有理化有了一个很好的方法!(Maple有自带的函数rationalize)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-1-21 12:36:22 | 显示全部楼层
  1. rad = -(1/(Sqrt[2] - Sqrt[3] + Sqrt[5] + Sqrt[6]));
  2. Select[x /. Solve[Factor[MinimalPolynomial[rad, x],
  3. Extension -> {Sqrt[2], Sqrt[3], Sqrt[5]}] == 0, x],
  4. Chop[#1 - rad] == 0 & ]//First
复制代码
这是一种根式分母有理化的方法,效率不算高,应该能优化的
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2013-1-21 13:25:29 | 显示全部楼层
mathematica分解这个很慢呀
  1. Factor[1 - 216 x^2 - 192 x^3 + 16140 x^4 + 18816 x^5 - 547528 x^6 - 687168 x^7 + 8960886 x^8 + 12394752 x^9 - 67518888 x^10 - 108989760 x^11 + 178031596 x^12 + 374357376 x^13 + 61149384 x^14 - 96214464 x^15 + 3999249 x^16, Extension -> {Sqrt[2], Sqrt[3], Sqrt[5], Sqrt[6], Sqrt[11]}]
复制代码
Maple的比较快
  1. factor(1-216*x^2-192*x^3+16140*x^4+18816*x^5-547528*x^6-687168*x^7+8960886*x^8+12394752*x^9-67518888*x^10-108989760*x^11+178031596*x^12+374357376*x^13+61149384*x^14-96214464*x^15+3999249*x^16, [sqrt(2), sqrt(3), sqrt(5), sqrt(6), sqrt(11)])
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-1-21 15:13:57 | 显示全部楼层
其中Sqrt[6]不需要。 如果比较慢的话,可以先Extension -> {Sqrt[2], Sqrt[3]},再将分解出来的每个因式分解,Extension -> {Sqrt[5], Sqrt[11]}
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-1-21 15:58:42 | 显示全部楼层

点评

nyy
你居然也知道这个,而且比我早了十年这样  发表于 2024-2-19 14:35
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

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

GMT+8, 2024-11-24 07:10 , Processed in 0.033816 second(s), 18 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

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