| 查看: 2195 | 回复: 18 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用 Mathematica 分解因式
|
|||
|
请教： 如何在实数（实代数数）范围内将整系数高次多项式分解为一次式与不可再分解的二次式的乘积？比如将 x^4+1 分解为 (x^2-√2 x+1)(x^2+√2 x+1)。以及类似地将有理分式分解为部分分式。 |
回复此楼» 猜你喜欢
化学求调剂
已经有5人回复
-
一志愿中国石油大学（华东） 本科齐鲁工业大学
已经有3人回复
-
332求调剂
已经有4人回复
-
265求调剂
已经有9人回复
-
328求调剂，英语六级551，有科研经历
已经有8人回复
-
求调剂
已经有3人回复
-
278求调剂
已经有6人回复
-
307求调剂
已经有10人回复
-
279分求调剂 一志愿211
已经有12人回复
-
297求调剂
已经有9人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
5.2版本时 上面代码中有些命令不支持 下面这个可以在5.2上运行 你试试4.0 代码写的很笨 可能有很大空间可以优化 也许有极简单的内部命令完成 就是不清楚 哪个选项 只是计算简单 不讲究效率也可以.... -------------------------------- 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]; t2 = Cases[t2, x_ /; Element[x, Reals]]; eq1 := {ComplexExpand[Im[t1[[j]]^2 - a*t1[[j]] + b]] == 0}; r1 = Table[Solve[eq1, a], {j, 1, Length[t1]}] // Flatten; r1t = N[r1, 10]; r1 = Table[If[! MemberQ[Take[r1t, j - 1], r1t[[j]]], r1[[j]], xxxx], {j, 1, Length[r1t]}]; r1 = DeleteCases[r1, xxxx]; (*r1 = DeleteDuplicates[r1]*) eq2 := {ComplexExpand[Re[t1[[j]]^2 - a*t1[[j]] + b]] == 0}; r2 = Table[Solve[eq2, b], {j, 1, Length[t1]}] // Flatten; (*r2 = DeleteDuplicates[r2]*) r2t = N[r1, 10]; r2 = Table[If[! MemberQ[Take[ r2t, j - 1], r2t[[j]]], r2[[j]], xxxx], {j, 1, Length[r2t]}]; r2 = DeleteCases[r2, xxxx]; r3 = Table[{r1[[j]], r2[[j]] /. r1[[j]]}, {j, 1, Length[r1]}]; ff1 = Table[temp1 = x^2 - a*x + b /. r3[[j]] // Expand; Print[temp1]; temp1, {j, 1, Length[r3]}]; ff2 = Table[temp1 = x - t2[[j]]; Print[temp1]; temp1, {j, 1, Length[t2]}]; ff = Join[ff1, ff2]; N[Apply[Times, ff], 20] |
17楼2012-12-16 13:27:41
2楼2012-12-11 11:25:31
3楼2012-12-11 17:55:38
4楼2012-12-11 17:57:38
