| 查看: 2239 | 回复: 18 | |||
[求助]
用 Mathematica 分解因式
|
|||
|
请教： 如何在实数（实代数数）范围内将整系数高次多项式分解为一次式与不可再分解的二次式的乘积？比如将 x^4+1 分解为 (x^2-√2 x+1)(x^2+√2 x+1)。以及类似地将有理分式分解为部分分式。 |
回复此楼» 猜你喜欢
2本，初试303，0860求调剂
已经有4人回复
-
286求调剂
已经有26人回复
-
335求调剂
已经有11人回复
-
293求调剂
已经有9人回复
-
290调剂生物0860
已经有25人回复
-
295分求调剂
已经有11人回复
-
一志愿华中农微生物，288分，三年实验经历
已经有11人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有4人回复
-
求助调剂，跨调
已经有12人回复
-
308求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
