ログインしてさらにmixiを楽しもう

コメントを投稿して情報交換!
更新通知を受け取って、最新情報をゲット!

独学ノート(土筆の子)コミュのもかいの問題の清書

  • mixiチェック
  • このエントリーをはてなブックマークに追加
http://userdisk.webry.biglobe.ne.jp/020/691/47/N000/000/001/132625293685913118455_index_gr_1_20120111123536.gif

x^3+3x^2-1=0 ①
の一つの解をαとし、他の2つの解をβ、γとする。

α^3+3α^2-1=0より、
(1/α)=α^2+3α ②

①より、
αβγ=1 したがって、βγ=(1/α)=α^2+3α ③
α+β+γ=-3 したがって、β+γ=-3-α     ④

③、④の関係から、求めるβ、γは次の2次方程式の解である。

X^2+(3+α)X+(α^2+3α)=0

これは、2次方程式の解の公式から、

β、γ=(1/2){(-3-α)±√{α+3)^2-4(α^2+3α)}

ルートの中は、
α^2+6α+9-4α^2-12α=-3α^2-6α+9  

とここまでで、前半。

ルートの中が2乗になるように係数を求める。
-3α^2-6α+9=(A*α^2 + B*α + C)^2

面倒なので、Mathmaticaを使いながら計算することにする。

(A*α^2 + B*α + C)^2
=C^2+2BCα+B^2α^2 + 2 A Cα^2 + 2 A Bα^3 + A^2α^4
=C^2+2BCα+B^2α^2 +2 ACα^2 +2 AB - 6 ABα^2 + A^2α-3 A^2(1 - 3α^2)
= -3 A^2 + 2 AB + C^2 + A^2α+2BCα+9 A^2α^2-6ABα^2 + 2 ACα^2 +B^2 α^2}
= -3 A^2 + 2 AB + C^2 + (A^2 + 2 BC)α+ (9 A^2 - 6 AB + 2 AC + B^2)α^2

Solve[-3*A^2 + 2*A*B + C^2 == 9 && A^2 + 2 B*C == -6 &&
9*A^2 - 6 A*B + 2*A*C + B^2 == -3, {A, B, C}]

{{A -> -2, B -> -5, C -> 1}, {A -> 2, B -> 5,
C -> -1}, {A ->
1/3 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3)),
B -> (1/2176)(6016/
3 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3)) -
832/3 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3))^3 +
176/9 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3))^5 -
1/3 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3))^7),
C -> (1/2176)(-(14336/
3) (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3)) +
1184/3 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3))^3 -
104/9 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3))^5 +
1/9 (2 + 2/(1/2 (-1 + I Sqrt[3]))^(1/3) +
2^(2/3) (-1 + I Sqrt[3])^(1/3))^7)},
以下省略。

コメント(5)

結局以下の2つが当てはまる。プラスの方をとる。

{A -> -2, B -> -5, C -> 1}, {A -> 2, B -> 5, C -> -1}

β、γ=(1/2){(-3-α)±(2α^2+5α-1)

= (1/2){-4+4α+2α^2},(1/2){-2-6α-2α^2}

=-2+2α+α^2, -1-3α-α^2


検算してみる。

確かに、α+β+γ=-3

α*β*γはアルファの文字が分かりにくくなるのでTにしてみる。

In[28]:= Expand[(-2 + 2 T + T^2)*(-1 - 3 T - T^2)*T]

Out[28]= 2 T + 4 T^2 - 5 T^3 - 5 T^4 - T^5

2T+4T^2-5(1-3T^2)-5T*(1-3T^2)-T^2(1-3T^2)

=-5 - 3 T + 18 T^2 + 15 T^3 + 3 T^4
=-5 - 3 T + 18 T^2 + 15(1-3T^2)+3T*(1-3T^2)

In[38]:= Expand[-5 - 3 T + 18 T^2 + 15 (1 - 3 T^2) + 3 T (1 - 3 T^2)]
Out[38]= 10 - 27 T^2 - 9 T^3

=10-27T^2-9(1-3T^2)
=10-9
=1

はじめて、最後まで、たどり着きました。 感動!!

それにしても、もっとうまいMathmaticaの使い方があるはず。

解と係数の関係で、不思議なのは、補間法。
(x,y)が(0,-1),(1,3),(-1,1),(2,19)を通れば、
3次式の係数は決まって、元の式が構成される。

In[14]:= (* Y=x^3+3x^2-1の近似式を求めるためのラグランジの補間法*)

(-1)*(x - 1) (x + 1) (x - 2)/{(0 - 1) (0 + 1) (0 - 2)} + (3)*(x -
0) (x + 1) (x - 2)/{(1 - 0) (1 + 1) (1 - 2)} + (1)*(x - 0) (x -
1) (x - 2)/{(-1 - 0) (-1 - 1) (-1 - 2)} + (19)*(x - 0) (x -
1) (x + 1)/{(2 - 0) (2 - 1) (2 + 1)}

Out[14]= {1/6 (2 - x) (-1 + x) x - 1/2 (-2 + x) (-1 + x) (1 + x) +
3/2 (2 - x) x (1 + x) + 19/6 (-1 + x) x (1 + x)}

In[15]:= Expand[{1/6 (2 - x) (-1 + x) x -
1/2 (-2 + x) (-1 + x) (1 + x) + 3/2 (2 - x) x (1 + x) +
19/6 (-1 + x) x (1 + x)}]

Out[15]= {-1 + 3 x^2 + x^3}
In[16]:= (y1)*(x - x2) (x -
x3) (x - x4)/{(x1 - x2) (x1 - x3) (x1 - x4)} + (y2)*(x - x1) (x -
x3) (x - x4)/{(x2 - x1) (x2 - x3) (x2 - x4)} + (y3)*(x - x1) (x -
x2) (x - x4)/{(x3 - x1) (x3 - x2) (x3 - x4)} + (y4)*(x - x1) (x -
x2) (x - x3)/{(x4 - x1) (x4 - x2) (x4 - x3)}

Out[16]= {((x - x2) (x - x3) (x - x4) y1)/((x1 - x2) (x1 - x3) (x1 -
x4)) + ((x - x1) (x - x3) (x - x4) y2)/((-x1 + x2) (x2 -
x3) (x2 - x4)) + ((x - x1) (x - x2) (x - x4) y3)/((-x1 +
x3) (-x2 + x3) (x3 - x4)) + ((x - x1) (x - x2) (x -
x3) y4)/((-x1 + x4) (-x2 + x4) (-x3 + x4))}

{(x^3 (-x2^2 x3 y1 + x2 x3^2 y1 + x2^2 x4 y1 - x3^2 x4 y1 -
x2 x4^2 y1 + x3 x4^2 y1 + x1^2 x3 y2 - x1 x3^2 y2 -
x1^2 x4 y2 + x3^2 x4 y2 + x1 x4^2 y2 - x3 x4^2 y2 -
x1^2 x2 y3 + x1 x2^2 y3 + x1^2 x4 y3 - x2^2 x4 y3 -
x1 x4^2 y3 + x2 x4^2 y3 + x1^2 x2 y4 - x1 x2^2 y4 -
x1^2 x3 y4 + x2^2 x3 y4 + x1 x3^2 y4 - x2 x3^2 y4))/
  ((-x1 +x2) (x1 - x3) (x2 - x3) (x1 - x4) (x2 - x4) (x3 -x4))

+ (x^2 (x2^3 x3 y1 - x2 x3^3 y1 - x2^3 x4 y1 + x3^3 x4 y1 + x2 x4^3 y1 - x3 x4^3 y1 - x1^3 x3 y2 +
x1 x3^3 y2 + x1^3 x4 y2 - x3^3 x4 y2 - x1 x4^3 y2 +
x3 x4^3 y2 + x1^3 x2 y3 - x1 x2^3 y3 - x1^3 x4 y3 +
x2^3 x4 y3 + x1 x4^3 y3 - x2 x4^3 y3 - x1^3 x2 y4 +
x1 x2^3 y4 + x1^3 x3 y4 - x2^3 x3 y4 - x1 x3^3 y4 +
x2 x3^3 y4))/((-x1 + x2) (x1 - x3) (x2 - x3) (x1 - x4) (x2 - x4) (x3 - x4)) +

(x (-x2^3 x3^2 y1 + x2^2 x3^3 y1 +
x2^3 x4^2 y1 - x3^3 x4^2 y1 - x2^2 x4^3 y1 + x3^2 x4^3 y1 +
x1^3 x3^2 y2 - x1^2 x3^3 y2 - x1^3 x4^2 y2 + x3^3 x4^2 y2 +
x1^2 x4^3 y2 - x3^2 x4^3 y2 - x1^3 x2^2 y3 + x1^2 x2^3 y3 +
x1^3 x4^2 y3 - x2^3 x4^2 y3 - x1^2 x4^3 y3 + x2^2 x4^3 y3 +
x1^3 x2^2 y4 - x1^2 x2^3 y4 - x1^3 x3^2 y4 + x2^3 x3^2 y4 +
x1^2 x3^3 y4 - x2^2 x3^3 y4))/((-x1 + x2) (x1 - x3) (x2 -
x3) (x1 - x4) (x2 - x4) (x3 - x4)) +

(x2^3 x3^2 x4 y1 -
x2^2 x3^3 x4 y1 - x2^3 x3 x4^2 y1 + x2 x3^3 x4^2 y1 +
x2^2 x3 x4^3 y1 - x2 x3^2 x4^3 y1 - x1^3 x3^2 x4 y2 +
x1^2 x3^3 x4 y2 + x1^3 x3 x4^2 y2 - x1 x3^3 x4^2 y2 -
x1^2 x3 x4^3 y2 + x1 x3^2 x4^3 y2 + x1^3 x2^2 x4 y3 -
x1^2 x2^3 x4 y3 - x1^3 x2 x4^2 y3 + x1 x2^3 x4^2 y3 +
x1^2 x2 x4^3 y3 - x1 x2^2 x4^3 y3 - x1^3 x2^2 x3 y4 +
x1^2 x2^3 x3 y4 + x1^3 x2 x3^2 y4 - x1 x2^3 x3^2 y4 -
x1^2 x2 x3^3 y4 +
x1 x2^2 x3^3 y4)/((-x1 + x2) (x1 - x3) (x2 - x3) (x1 - x4) (x2 -
x4) (x3 - x4))}

ログインすると、みんなのコメントがもっと見れるよ

mixiユーザー
ログインしてコメントしよう!

独学ノート(土筆の子) 更新情報

独学ノート(土筆の子)のメンバーはこんなコミュニティにも参加しています

星印の数は、共通して参加しているメンバーが多いほど増えます。

人気コミュニティランキング