Как я могу найти критерии того, что многочлен имеет и имеет только двойные корни?
Aug 19 2020
В справочном документе функции SolveAlways
я видел условие, что многочлен имеет тройной корень.
(*Find a condition for a cubic polynomial to have a triple root:*)
f[x_] := x^3 + a x^2 + b x + c;
SolveAlways[Implies[f[x] == 0 && f[y] == 0, x == y], {x, y}]
Теперь я хочу воспроизвести этот пример, чтобы найти условие, при котором многочлен имеет только двойной корень (не существует нескольких корней выше квадратичного).
A = {{1, 2, -3}, {-1, 4, -3}, {1, a, 5}};
f[λ_] := CharacteristicPolynomial[A, λ]
Reduce[Exists[{x, y}, Implies[f[x] == 0 && f[y] == 0, x == y]], {a}]
Solve[(f[x] /. a -> -2) == 0, x]
SolveAlways[Implies[(f[x] == f'[x] == 0), f''[x] != 0], x]
Но вывод приведенного выше кода не является условием оценки. Что я могу сделать, чтобы решить эту проблему?
Ответы
7 DanielLichtblau Aug 19 2020 at 22:46
Вы хотите, чтобы дискриминант обращался в нуль, а дискриминант производной не обращался в нуль.
mat = {{1, 2, -3}, {-1, 4, -3}, {1, a, 5}};
cpoly = CharacteristicPolynomial[mat, x];
disc1 = Discriminant[cpoly, x]
disc2 = Discriminant[D[cpoly, x], x]
(* Out[616]= -288 - 720 a - 504 a^2 - 108 a^3
Out[617]= -4 (2 + 9 a) *)
Теперь разберитесь в этом.
Reduce[disc1 == 0 && disc2 != 0]
(* Out[618]= a == -2 || a == -(2/3) *)
7 user64494 Aug 19 2020 at 23:47
Используя определение корня кратности 2, получаем
A = {{1, 2, -3},{-1, 4, -3},{1, a, 5}};f[\[Lambda]_]:= CharacteristicPolynomial[A,\[Lambda]]
Reduce[Exists[x, f[x] == 0 && f'[x] == 0 && f''[x] != 0], a, Reals]
(*a == -2 || a == -(2/3)*)
Этот подход работает не только для многочленов.
6 J.M.'sennui Aug 19 2020 at 14:02
Вот для чего можно использовать дискриминант :
Reduce[Discriminant[CharacteristicPolynomial[{{1, 2, -3},
{-1, 4, -3},
{1, a, 5}}, x], x] == 0, a]
a == -2 || a == -2/3
Проверьте:
With[{a = -2/3}, Eigenvalues[{{1, 2, -3}, {-1, 4, -3}, {1, a, 5}}]]
{4, 4, 2}
With[{a = -2}, Eigenvalues[{{1, 2, -3}, {-1, 4, -3}, {1, a, 5}}]]
{6, 2, 2}