| 24小時(shí)熱門版塊排行榜 |
| 查看: 2208 | 回復(fù): 18 | |||
xxppyy木蟲(chóng) (小有名氣)
|
[求助]
用 Mathematica 分解因式
|
||
|
請(qǐng)教: 如何在實(shí)數(shù)(實(shí)代數(shù)數(shù))范圍內(nèi)將整系數(shù)高次多項(xiàng)式分解為一次式與不可再分解的二次式的乘積?比如將 x^4+1 分解為 (x^2-√2 x+1)(x^2+√2 x+1)。以及類似地將有理分式分解為部分分式。 |
金蟲(chóng) (著名寫(xiě)手)
|
上面代碼顯示好像不對(duì) 下面這個(gè) 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 |
新蟲(chóng) (初入文壇)

新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
木蟲(chóng) (小有名氣)
金蟲(chóng) (著名寫(xiě)手)
|
我寫(xiě)了個(gè)很笨的代碼求這個(gè)問(wèn)題 可能內(nèi)部有非常簡(jiǎn)潔的方法 另外 也沒(méi)怎么優(yōu)化 看樣子 能求你說(shuō)的情況 ---------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) |
金蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (小有名氣)
金蟲(chóng) (著名寫(xiě)手)
|
下面這個(gè)代碼應(yīng)該可以直接copy過(guò)去運(yùn)行 我用mathematica8.04 第一行的函數(shù)可以改成其他多項(xiàng)式 -------------------- 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] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +11 | AZMK 2026-04-01 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 安全工程 285 求調(diào)劑 +3 | Xinyu56 2026-04-01 | 4/200 |
|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +10 | GENJIOW 2026-03-31 | 11/550 |
|
|
[考研] 307分求調(diào)劑 +14 | (o~o) 2026-03-31 | 15/750 |
|
|
[考研] 324求調(diào)劑 +10 | hanamiko 2026-03-26 | 12/600 |
|
|
[考研] 326求調(diào)劑 +4 | 崽崽仔 2026-03-31 | 4/200 |
|
|
[考研] 0710生物學(xué)求調(diào)劑! +6 | 敘述文 2026-03-31 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +10 | @taotao 2026-03-29 | 10/500 |
|
|
[考研] 311(085601)求調(diào)劑 +12 | liziyeyeye 2026-03-28 | 13/650 |
|
|
[考研] 08開(kāi)頭看過(guò)來(lái)!! +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
[考研] 本科211安全工程,初試290分,求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚(yú) 2026-03-27 | 12/600 |
|
|
[考研]
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|