| 查看: 2436 | 回复: 18 | ||
[求助]
用 Mathematica 分解因式
|
||
|
请教： 如何在实数（实代数数）范围内将整系数高次多项式分解为一次式与不可再分解的二次式的乘积？比如将 x^4+1 分解为 (x^2-√2 x+1)(x^2+√2 x+1)。以及类似地将有理分式分解为部分分式。 |
回复此楼» 猜你喜欢
filecode，4个jtjc了
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有117人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
-
今天放榜吗？
已经有13人回复
-
今天基金会出结果吗？20260819
已经有8人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
-
时间戳变了，能看出什么问题？
已经有14人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有30人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- a^m-b^n因式分解
已经有8人回复
- x^3+x^2+3x+4在复数域分解因式!!!
已经有5人回复
- 【软件】Mathmatica 软件 7.0版
已经有68人回复
- 一元高次方程因式分解
已经有14人回复
- 【求助】用matlab能够将集合进行因式分解简化吗？
已经有5人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★
jjdg: 金币+1, 感谢参与 2012-12-16 00:19:35
jjdg: 金币+1, 感谢参与 2012-12-16 00:19:35
|
上面代码显示好像不对 下面这个 f[x_] := x^7 + 2 t1 = x /. Solve[f[x] == 0, x]; t1 = DeleteCases[t1, x_ /; Element[x, Reals]]; t2 = x /. Solve[f[x] == 0, x, Reals]; r1 = Table[ Solve[ComplexExpand[Im[t1[]^2 - a*t1[] + b]] == 0, a], {i, 1, Length[t1]}] // Flatten; r1 = DeleteDuplicates[r1]; r2 = Table[ Solve[ComplexExpand[Re[t1[]^2 - a*t1[] + b]] == 0, b], {i, 1, Length[t1]}] // Flatten; r2 = DeleteDuplicates[r2]; r3 = Table[{r1[], r2[] /. r1[]}, {i, 1, Length[r1]}]; ff1 = Table[temp1 = x^2 - a*x + b /. r3[] // Expand; Print[temp1]; temp1, {i, 1, Length[r3]}]; ff2 = Table[temp1 = x - t2[]; Print[temp1]; temp1, {i, 1, Length[t2]}]; ff = Join[ff1, ff2]; Apply[Times, ff] Apply[Times, ff] // N // Chop |
7楼2012-12-15 13:52:59
2楼2012-12-11 11:25:31
3楼2012-12-11 17:55:38
4楼2012-12-11 17:57:38
5楼2012-12-13 15:02:59
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★ ★
jjdg: 金币+2, 感谢参与 2012-12-16 00:19:27
jjdg: 金币+2, 感谢参与 2012-12-16 00:19:27
|
我写了个很笨的代码求这个问题 可能内部有非常简洁的方法 另外 也没怎么优化 看样子 能求你说的情况 ---------code------------- Clear["Global`*"]; f[x_] := x^7 + 2 t1 = x /. Solve[f[x] == 0, x]; t1 = DeleteCases[t1, x_ /; Element[x, Reals]]; t2 = x /. Solve[f[x] == 0, x, Reals]; r1 = Table[ Solve[ComplexExpand[Im[t1[]^2 - a*t1[] + b]] == 0, a], {i, 1, Length[t1]}] // Flatten; r1 = DeleteDuplicates[r1]; r2 = Table[ Solve[ComplexExpand[Re[t1[]^2 - a*t1[] + b]] == 0, b], {i, 1, Length[t1]}] // Flatten; r2 = DeleteDuplicates[r2]; r3 = Table[{r1[], r2[] /. r1[]}, {i, 1, Length[r1]}]; ff1 = Table[temp1 = x^2 - a*x + b /. r3[] // Expand; Print[temp1]; temp1, {i, 1, Length[r3]}]; ff2 = Table[temp1 = x - t2[]; Print[temp1]; temp1, {i, 1, Length[t2]}]; ff = Join[ff1, ff2]; Apply[Times, ff] Apply[Times, ff] // N // Chop ------------------------------------------------------ 结果： 解析形式： (x+Power[2, (7)^-1]) (x^2-Power[2, (7)^-1] x cos((3 \[Pi])/14) csc(\[Pi]/7)-2^(2/7) sin((3 \[Pi])/14)+2^(2/7) cos((3 \[Pi])/14) cot(\[Pi]/7)) (x^2-Power[2, (7)^-1] x sin(\[Pi]/7) sec(\[Pi]/14)+2^(2/7) cos(\[Pi]/7)+2^(2/7) sin(\[Pi]/7) tan(\[Pi]/14)) (x^2+Power[2, (7)^-1] x cos(\[Pi]/14) sec((3 \[Pi])/14)+2^(2/7) sin(\[Pi]/14)+2^(2/7) cos(\[Pi]/14) tan((3 \[Pi])/14)) 数值化： (x+1.10409) (x^2-1.9895 x+1.21901) (x^2-0.491366 x+1.21901) (x^2+1.37678 x+1.21901) |
6楼2012-12-15 13:50:02
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
8楼2012-12-15 13:55:08
9楼2012-12-15 18:32:54
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
下面这个代码应该可以直接copy过去运行 我用mathematica8.04 第一行的函数可以改成其他多项式 -------------------- f[x_] := x^7 + 2 t1 = x /. Solve[f[x] == 0, x]; t1 = DeleteCases[t1, x_ /; Element[x, Reals]]; t2 = x /. Solve[f[x] == 0, x, Reals]; eq1 := {ComplexExpand[Im[t1[]^2 - a*t1[] + b]] == 0}; r1 = Table[Solve[eq1, a], {i, 1, Length[t1]}] // Flatten; r1 = DeleteDuplicates[r1]; eq2 := {ComplexExpand[Re[t1[]^2 - a*t1[] + b]] == 0}; r2 = Table[Solve[eq2, b], {i, 1, Length[t1]}] // Flatten; r2 = DeleteDuplicates[r2]; r3 = Table[{r1[], r2[] /. r1[]}, {i, 1, Length[r1]}]; ff1 = Table[temp1 = x^2 - a*x + b /. r3[] // Expand; Print[temp1]; temp1, {i, 1, Length[r3]}]; ff2 = Table[temp1 = x - t2[]; Print[temp1]; temp1, {i, 1, Length[t2]}]; ff = Join[ff1, ff2]; N[Apply[Times, ff], 20] |
10楼2012-12-16 10:52:12
