 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2179 | 回復: 18 | ||
[求助]
用 Mathematica 分解因式
|
||
|
請教： 如何在實數(shù)（實代數(shù)數(shù)）范圍內(nèi)將整系數(shù)高次多項式分解為一次式與不可再分解的二次式的乘積？比如將 x^4+1 分解為 (x^2-√2 x+1)(x^2+√2 x+1)。以及類似地將有理分式分解為部分分式。 |
回復此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有16人回復
-
304求調(diào)劑
已經(jīng)有6人回復
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復
-
求材料調(diào)劑
已經(jīng)有8人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復
-
一志愿華中科技大學，080502，354分求調(diào)劑
已經(jīng)有4人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有6人回復
-
085410人工智能專碩317求調(diào)劑（0854都可以）
已經(jīng)有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- a^m-b^n因式分解
已經(jīng)有8人回復
- x^3+x^2+3x+4在復數(shù)域分解因式!!!
已經(jīng)有5人回復
- 【軟件】Mathmatica 軟件 7.0版
已經(jīng)有68人回復
- 一元高次方程因式分解
已經(jīng)有14人回復
- 【求助】用matlab能夠?qū)⒓线M行因式分解簡化嗎？
已經(jīng)有5人回復
walk1997
金蟲 (著名寫手)
- 應助: 1 (幼兒園)
- 金幣: 4676.2
- 紅花: 22
- 帖子: 1066
- 在線: 798.1小時
- 蟲號: 416039
- 注冊: 2007-06-29
- 性別: GG
- 專業(yè): 粒子物理學和場論
★
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
- 專業(yè): 粒子物理學和場論
★ ★
jjdg: 金幣+2, 感謝參與 2012-12-16 00:19:27
jjdg: 金幣+2, 感謝參與 2012-12-16 00:19:27
|
我寫了個很笨的代碼求這個問題 可能內(nèi)部有非常簡潔的方法 另外 也沒怎么優(yōu)化 看樣子 能求你說的情況 ---------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 ------------------------------------------------------ 結(jié)果： 解析形式： (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)) 數(shù)值化： (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
- 專業(yè): 粒子物理學和場論
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
- 專業(yè): 粒子物理學和場論
|
下面這個代碼應該可以直接copy過去運行 我用mathematica8.04 第一行的函數(shù)可以改成其他多項式 -------------------- 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
相關版塊跳轉(zhuǎn)
我要訂閱樓主 xxppyy 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 東南大學364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|