| 查看: 2261 | 回复: 18 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用 Mathematica 分解因式
|
|||
|
请教： 如何在实数（实代数数）范围内将整系数高次多项式分解为一次式与不可再分解的二次式的乘积？比如将 x^4+1 分解为 (x^2-√2 x+1)(x^2+√2 x+1)。以及类似地将有理分式分解为部分分式。 |
回复此楼» 猜你喜欢
二苯甲酮酸类衍生物
已经有4人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有11人回复
-
294求调剂
已经有6人回复
-
335求调剂
已经有23人回复
-
279求调剂
已经有13人回复
-
304求调剂
已经有6人回复
-
322求调剂
已经有3人回复
-
271求调剂
已经有39人回复
-
收到复试调剂但是去不了
已经有6人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有5人回复
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
