Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Для нахождения канонического вида B ортогональной матрицы A и ортогональной матрицы Q, нам нужно выполнить процесс диагонализации A. Ортогональная матрица Q будет состоять из собственных векторов A, и B будет иметь диагональную структуру с собственными значениями A на диагонали.
- 1. Найдем собственные значения A, которые будут являться диагональными элементами B. Это можно сделать, решив характеристическое уравнение:
- det(A - λI) = 0, где λ - собственное значение, I - единичная матрица.
- Для матрицы A:
- |2/3 - λ, -1/3, 2/3|
- |2/3, 2/3 - λ, -1/3|
- |-1/3, 2/3, 2/3 - λ|
- Рассчитаем определитель этой матрицы:
- det(A - λI) = (2/3 - λ)[(2/3 - λ)(2/3 - λ) - (-1/3)(2/3)] - (-1/3)[2/3(2/3 - λ) - (-1/3)(-1/3)] + 2/3[2/3(2/3 - λ) - (-1/3)(-1/3)]
- Упростим выражение:
- det(A - λI) = (2/3 - λ)[(4/9 - 4/3λ + λ^2) + 2/9] + (1/9)[(4/3 - 4/3λ) - 1/9] + (2/9)[(4/9 - 4/3λ) - 1/9]
- det(A - λI) = (2/3 - λ)[(2/3 - 4/3λ + λ^2) + 2/9] + (1/9)[(4/3 - 4/3λ) - 1/9] + (2/9)[(4/9 - 4/3λ) - 1/9]
- det(A - λI) = (2/3 - λ)[2/3 - 4/3λ + λ^2 + 2/9] + (1/9)[4/3 - 4/3λ - 1/9] + (2/9)[4/9 - 4/3λ - 1/9]
- det(A - λI) = (2/3 - λ)[10/9 - 4/3λ + λ^2] + (1/9)[26/9 - 4/3λ] + (2/9)[26/9 - 4/3λ]
- Упростим дальше:
- det(A - λI) = (2/3 - λ)[10/9 - 4/3λ + λ^2] + (1/9)[26/9 - 4/3λ] + (2/9)[26/9 - 4/3λ]
- det(A - λI) = (2/3 - λ)[10/9 - 4/3λ + λ^2] + (1/27)[26 - 12λ] + (2/27)[26 - 12λ]
- det(A - λI) = (2/3 - λ)[10/9 - 4/3λ + λ^2] + (1/27)[26 - 12λ + 52 - 12λ]
- det(A - λI) = (2/3 - λ)[10/9 - 4/3λ + λ^2] + (1/27)[78 - 24λ]
- det(A - λI) = (2/3 - λ)[10/9 - 4/3λ + λ^2] + (2/3)[13 - 4λ]
- det(A - λI) = (2/3 - λ)[10/9 - 4/3λ + λ^2 + 26/9 - 8λ]
- det(A - λI) = (2/3 - λ)[36/9 - 12/9λ + 9/9λ^2 + 26/9 - 8λ]
- det(A - λI) = (2/3 - λ)[71/9 - 20/9λ + 9/9λ^2]
- det(A - λI) = (2/3 - λ)[71/9 - 20/9λ + λ^2]
- Теперь решим характеристическое уравнение:
- (2/3 - λ)[71/9 - 20/9λ + λ^2] = 0
- Решения этого уравнения будут собственными значениями A. Найдем их:
- 1. (2/3 - λ) = 0 => λ1 = 2/3
- 2. 71/9 - 20/9λ + λ^2 = 0
- Умножим обе стороны на 9, чтобы избавиться от дробей:
- 71 - 20λ + 9λ^2 = 0
- Это квадратное уравнение. Решим его с помощью дискриминанта:
- D = (-20)^2 - 4 * 9 * 71 = 400 - 2556 = -2156
- Поскольку дискриминант отрицателен, у нас есть два комплексных собственных значения:
- λ2 = (-(-20) + √(-2156)) / (2 * 9) ≈ (20 + 14.7i) / 18 ≈ (10/9 + 7.35i/9)
- λ3 = (-(-20) - √(-2156)) / (2 * 9) ≈ (20 - 14.7i) / 18 ≈ (10/9 - 7.35i/9)
- Теперь у нас есть собственные значения A: λ1 = 2/3, λ2 = (10/9 + 7.35i/9), λ3 = (10/9 - 7.35i/9).
- 2. Теперь нам нужно найти собственные векторы, соответствующие этим собственным значениям. Начнем с λ1 = 2/3:
- (2/3 * I
- - A)v1 = 0, где I - единичная матрица, v1 - собственный вектор, соответствующий λ1.
- (2/3 * I - A)v1 = 0
- (2/3 * [[1, 0, 0], [0, 1, 0], [0, 0, 1]] - [[2/3, -1/3, 2/3], [2/3, 2/3, -1/3], [-1/3, 2/3, 2/3]])v1 = 0
- (2/3 * [[1, 0, 0], [0, 1, 0], [0, 0, 1]] - [[2/3, -1/3, 2/3], [2/3, 2/3, -1/3], [-1/3, 2/3, 2/3]])v1 = 0
- (2/3 - 2/3) * v1 = 0
- 0 * v1 = 0
- Это уравнение сводится к тривиальному уравнению 0 = 0, и его решением будет любой вектор v1. Например, вектор [1, 0, 0] можно выбрать в качестве собственного вектора, соответствующего λ1 = 2/3.
- 3. Теперь найдем собственные векторы, соответствующие комплексным собственным значениям λ2 и λ3. Для этого мы рассмотрим матрицу (2/3 * I - A) и воспользуемся методом Гаусса для решения систем линейных уравнений. Поскольку значения комплексные, я приведу общий подход.
- Рассмотрим матрицу (2/3 * I - A) - λ2 * I:
- [[2/3 - 2/3, -1/3, 2/3], [2/3, 2/3 - 2/3, -1/3], [-1/3, 2/3, 2/3 - 2/3]] - (10/9 + 7.35i/9) * [[1, 0, 0], [0, 1, 0], [0, 0, 1]]
- [[0, -1/3, 2/3], [2/3, 0, -1/3], [-1/3, 2/3, 0]] - (10/9 + 7.35i/9) * [[1, 0, 0], [0, 1, 0], [0, 0, 1]]
- [[0, -1/3, 2/3], [2/3, 0, -1/3], [-1/3, 2/3, 0]] - [[10/9 + 7.35i/9, 0, 0], [0, 10/9 + 7.35i/9, 0], [0, 0, 10/9 + 7.35i/9]]
- Теперь мы решим систему линейных уравнений ((2/3 * I - A) - λ2 * I)v2 = 0, где v2 - собственный вектор, соответствующий λ2. Это можно сделать с помощью метода Гаусса. Рассмотрим расширенную матрицу и приведем ее к ступенчатому виду:
- [[0, -1/3, 2/3, 0], [2/3, 0, -1/3, 0], [-1/3, 2/3, 0, 0]] - [[10/9 + 7.35i/9, 0, 0, 0], [0, 10/9 + 7.35i/9, 0, 0], [0, 0, 10/9 + 7.35i/9, 0]]
- [[0, -1/3, 2/3, 0], [2/3, 0, -1/3, 0], [-1/3, 2/3, 0, 0]] - [[10/9 + 7.35i/9, 0, 0, 0], [0, 10/9 + 7.35i/9, 0, 0], [0, 0, 10/9 + 7.35i/9, 0]]
- [[0, -1/3, 2/3, 0], [2/3, 0, -1/3, 0], [-1/3, 2/3, 0, 0]] - [[10/9 + 7.35i/9, 0, 0, 0], [0, 10/9 + 7.35i/9, 0, 0], [0, 0, 10/9 + 7.35i/9, 0]]
- Теперь используем метод Гаусса для приведения этой матрицы к ступенчатому виду:
- 1. Поменяем местами первую и третью строки:
- [[0, -1/3, 2/3, 0], [2/3, 0, -1/3, 0], [-1/3, 2/3, 0, 0]] - [[10/9 + 7.35i/9, 0, 0, 0], [0, 10/9 + 7.35i/9, 0, 0], [0, 0, 10/9 + 7.35i/9, 0]]
- [[0, -1/3, 2/3, 0], [2/3, 0, -1/3, 0], [-1/3, 2/3, 0, 0]] - [[0, 0, 0, 0], [0, 0, 0, 0], [10/9 + 7.35i/9, 0, 0, 0]]
- 2. Первую строку умножим на (-1/3):
- [[0, 1, -2, 0], [2/3, 0, -1/3, 0], [-1/3, 2/3, 0,
- 0]] - [[0, 0, 0, 0], [0, 0, 0, 0], [10/9 + 7.35i/9, 0, 0, 0]]
- 3. Выразим во второй строке x через y:
- [[0, 1, -2, 0], [1, 0, -1/2, 0], [-1/3, 2/3, 0, 0]] - [[0, 0, 0, 0], [0, 0, 0, 0], [10/9 + 7.35i/9, 0, 0, 0]]
- 4. Вторую строку умножим на (-1) и прибавим к первой строке:
- [[0, 1, -2, 0], [1, 0, -1/2, 0], [0, 2/3, 0, 0]] - [[0, 0, 0, 0], [0, 0, 0, 0], [10/9 + 7.35i/9, 0, 0, 0]]
- 5. Вторую строку умножим на (2/3):
- [[0, 1, -2, 0], [1, 0, -1/2, 0], [0, 2/3, 0, 0]] - [[0, 0, 0, 0], [0, 0, 0, 0], [10/9 + 7.35i/9, 0, 0, 0]]
- 6. Третью строку умножим на (3/2) и прибавим к первой строке:
- [[0, 1, -2, 0], [1, 0, -1/2, 0], [0, 1, 0, 0]] - [[0, 0, 0, 0], [0, 0, 0, 0], [10/9 + 7.35i/9, 0, 0, 0]]
- 7. Выразим в третьей строке y через z:
- [[0, 1, -2, 0], [1, 0, -1/2, 0], [0, 1, 0, 0]] - [[0, 0, 0, 0], [0, 0, 0, 0], [0, 0, 0, 0]]
- Теперь матрица находится в ступенчатом виде, и мы можем найти собственный вектор, соответствующий комплексному собственному значению λ2. Это вектор [0, 0, 1].
- 8. Аналогично, для λ3 = (10/9 - 7.35i/9) решим систему ((2/3 * I - A) - λ3 * I)v3 = 0:
- [[0, -1/3, 2/3, 0], [2/3, 0, -1/3, 0], [-1/3, 2/3, 0, 0]] - [[10/9 - 7.35i/9, 0, 0, 0], [0, 10/9 - 7.35i/9, 0, 0], [0, 0, 10/9 - 7.35i/9, 0]]
- Процесс будет аналогичным, и мы найдем собственный вектор, соответствующий λ3: [0, 0, 1].
- Теперь у нас есть собственные значения и их соответствующие собственные векторы:
- λ1 = 2/3, v1 = [1, 0, 0]
- λ2 = 10/9 + 7.35i/9, v2 = [0, 0, 1]
- λ3 = 10/9 - 7.35i/9, v3 = [0, 0, 1]
- 4. Теперь соберем матрицу Q из собственных векторов:
- Q = [v1, v2, v3]
- Q = [[1, 0, 0], [0, 0, 0], [0, 1, 1]]
- 5. Теперь найдем обратную матрицу Q^(-1). Для этого нам нужно найти определитель матрицы Q и затем вычислить матрицу алгебраических дополнений, транспонировать ее и умножить на 1/det(Q). Определитель матрицы Q равен 1, так как у нее есть нулевые строки.
- Таким образом, Q^(-1) = Q^T, где Q^T - матрица, полученная транспонированием матрицы Q:
- Q^(-1) = [[1, 0, 0], [0, 0, 1], [0, 0, 1]]
- 6. Теперь мы можем найти канонический вид B:
- B = Q^(-1) * A * Q
- B = [[1, 0, 0], [0, 0, 1], [0, 0, 1]] * [[2/3, -1/3, 2/3], [2/3, 2/3, -1/3], [-1/3, 2/3, 2/3]] * [[1, 0, 0], [0, 0, 1], [0, 0, 1]]
- B = [[(1*(2/3) + 0*(2/3) + 0*(-1/3)), (1*(-1/3) + 0*(2/3) + 0*(-1/3)), (1*(2/3) + 0*(-1/3) + 0*(2/3))],
- [(0*(2/3) + 0*(2/3) + 0*(-1/3)), (0*(-1/3) + 0*(2/3) + 1*(-1/3)), (0*(2/3) + 0*(-1/3) + 1*(2/3))],
- [(0*(2/3) + 1*(2/3) + 0*(-1/3)), (0
- *(-1/3) + 0*(2/3) + 0*(-1/3)), (0*(2/3) + 1*(-1/3) + 0*(2/3))]]
- B = [[2/3, 0, 0],
- [0, -1/3, 0],
- [0, 0, 2/3]]
- Таким образом, канонический вид B ортогональной матрицы A равен:
- B = [[2/3, 0, 0],
- [0, -1/3, 0],
- [0, 0, 2/3]]
- Ортогональная матрица Q:
- Q = [[1, 0, 0],
- [0, 0, 1],
- [0, 1, 1]]
- И выше мы нашли, что Q^(-1) = Q^T:
- Q^(-1) = [[1, 0, 0],
- [0, 0, 1],
- [0, 1, 1]]
Advertisement
Add Comment
Please, Sign In to add comment