РЕШЕНИЕ МАТРИЧНЫХ УРАВНЕНИЙ

Решение матричных уравнений

Матричные уравнения бывают трех типов.

  • 1. [latex]A \cdot X=B[/latex]
  • 2. [latex]X \cdot A=B[/latex]
  • 3. [latex]C \cdot X \cdot A=B[/latex]
  • Пример 1. Чтобы решить уравнение первого типа нужно обе части уравнения умножить на обратную к матрице [latex]A[/latex] слева.
    [latex]\begin{pmatrix} 1 & 2 \\ 3 & 4 \\ \end{pmatrix} \cdot X=[/latex] [latex]\begin{pmatrix} 3 & 5 \\ 5 & 9 \\ \end{pmatrix}[/latex], [latex]\det \begin{pmatrix} 1 & 2 \\ 3 & 4 \\ \end{pmatrix}=-2[/latex]
    [latex]A_{11}=(-1)^{1+1} \cdot 4=4[/latex]
    [latex]A_{12}=(-1)^{1+2} \cdot 3=-3[/latex]
    [latex]A_{21}=(-1)^{2+1} \cdot 2=-2[/latex]
    [latex]A_{22}=(-1)^{2+2} \cdot 1=1[/latex]
    [latex]\begin{pmatrix} 4 & -3 \\ -2 & 1 \\ \end{pmatrix}[/latex], полученную матрицу транспонируем и умножим на [latex]\det^{-1} \begin{pmatrix} 1 & 2 \\ 3 & 4 \\ \end{pmatrix}=-1/2[/latex]. Обратная матрица к [latex]\begin{pmatrix} 1 & 2 \\ 3 & 4 \\ \end{pmatrix}[/latex] равна [latex]\begin{pmatrix} -2 & 1 \\ 3/2 & -1/2 \\ \end{pmatrix}[/latex].
    [latex]X=\begin{pmatrix} -2 & 1 \\ 3/2 & -1/2 \\ \end{pmatrix} \cdot [/latex] [latex] \begin{pmatrix} 3 & 5 \\ 5 & 9 \\ \end{pmatrix}[/latex], [latex]X= \begin{pmatrix} -1 & -1 \\ 2 & 3 \\ \end{pmatrix}[/latex]. Сделаем проверку [latex]\begin{pmatrix} 1 & 2 \\ 3 & 4 \\ \end{pmatrix} \cdot \begin{pmatrix} -1 & -1 \\ 2 & 3 \\ \end{pmatrix}=[/latex][latex]\begin{pmatrix} 3 & 5 \\ 5 & 9 \\ \end{pmatrix}[/latex]. Уравнение решили правильно.
    Пример 2. Чтобы решить уравнение второго типа нужно обе части уравнения умножить на обратную к матрице [latex]A[/latex] справа.
    [latex]X \cdot \begin{pmatrix} 3 & -2 \\ 5 & -4 \\ \end{pmatrix}=[/latex] [latex]\begin{pmatrix} -1 & 2 \\ -5 & 6 \\ \end{pmatrix}.[/latex] Матрица обратная к [latex]\begin{pmatrix} 3 & -2 \\ 5 & -4 \\ \end{pmatrix}[/latex] равна [latex]\begin{pmatrix} 2 & -1 \\ 5/2 & -3/2 \\ \end{pmatrix}.[/latex] [latex]X=\begin{pmatrix} -1 & 2 \\ -5 & 6 \\ \end{pmatrix} \cdot \begin{pmatrix} 2 & -1 \\ 5/2 & -3/2 \\ \end{pmatrix},[/latex] [latex]X=\begin{pmatrix} 3 & -2 \\ 5 & -4 \\ \end{pmatrix}[/latex].
    Пример 3. Чтобы решить уравнение третьего типа нужно обе части уравнения умножить на обратную к матрице [latex]A[/latex] справа и на обратную матрице [latex]C[/latex] слева.
    [latex]\begin{pmatrix} 3 & -1 \\ 5 & -2 \\ \end{pmatrix} \cdot X \cdot \begin{pmatrix} 5 & 6 \\ 7 & 8 \\ \end{pmatrix}=[/latex] [latex]\begin{pmatrix} 14 & 16 \\ 9 & 10 \\ \end{pmatrix}[/latex]. Обратная матрица к [latex]\begin{pmatrix} 3 & -1 \\ 5 & -2 \\ \end{pmatrix}[/latex] равна [latex]\begin{pmatrix} 2 & -1 \\ 5 & -3 \\ \end{pmatrix},[/latex] обратная матрица к [latex]\begin{pmatrix} 5 & 6 \\ 7 & 8 \\ \end{pmatrix}[/latex] равна [latex]\begin{pmatrix} -4 & 3 \\ 7/2 & -5/2 \\ \end{pmatrix}[/latex]. [latex]X=\begin{pmatrix} 2 & -1 \\ 5 & -3 \\ \end{pmatrix} \cdot \begin{pmatrix} 14 & 16 \\ 9 & 10 \\ \end{pmatrix} \cdot [/latex][latex] \begin{pmatrix} -4 & 3 \\ 7/2 & -5/2 \\ \end{pmatrix}=\begin{pmatrix} 1 & 2 \\ 3 & 4 \\ \end{pmatrix}[/latex].
    Проверка [latex]\begin{pmatrix} 3 & -1 \\ 5 & -2 \\ \end{pmatrix} \cdot \begin{pmatrix} 1 & 2 \\ 3 & 4 \\ \end{pmatrix} \cdot \begin{pmatrix} 5 & 6 \\ 7 & 8 \\ \end{pmatrix}=[/latex] [latex]\begin{pmatrix} 14 & 16 \\ 9 & 10 \\ \end{pmatrix}[/latex].
    Пример 4. Случай когда обратная матрица не существует.
    [latex]X \cdot \begin{pmatrix} 3 & 6 \\ 4 & 8 \\ \end{pmatrix}=[/latex] [latex]\begin{pmatrix} 2 & 4 \\ 9 & 18 \\ \end{pmatrix}[/latex].
    Матрицу [latex]X[/latex] запишем как [latex]\begin{pmatrix} x_{1} & x_{2} \\ x_{3} & x_{4} \\ \end{pmatrix}[/latex], [latex]\begin{pmatrix} 3 \cdot x_{1}+4 \cdot x_{2} & 6 \cdot x_{1}+8 \cdot x_{2} \\ 3 \cdot x_{3}+4 \cdot x_{4} & 6 \cdot x_{3}+8 \cdot x_{4} \\ \end{pmatrix}=[/latex][latex]\begin{pmatrix} 2 & 4 \\ 9 & 18 \\ \end{pmatrix}[/latex].

    \begin{cases}
    3 \cdot x_{1}+4 \cdot x_{2} = 2\\
    6 \cdot x_{1}+8 \cdot x_{2} = 4\\
    3 \cdot x_{3}+4 \cdot x_{4} = 9\\
    6 \cdot x_{3}+8 \cdot x_{4}=18
    \end{cases}
    Эта система эквивалентна
    \begin{cases}
    3 \cdot x_{1}+4 \cdot x_{2} = 2\\
    3 \cdot x_{3}+4 \cdot x_{4} = 9
    \end{cases}
    Решив данную систему получим общей вид решения [latex]X=\begin{pmatrix} x_{1} & (2-3x_{1})/4 \\ x_{3} & (9-4x_{1})/3 \\ \end{pmatrix}[/latex]
    Литература

  • 1. Белозёров Г. С. Конспект по алгебре и геометрии
  • 2. Линейная алгебра. Воеводин. В. В. М.: Наука. Главная редакция физико-математической литературы, 1980 год, стр. 211-213.
  • Сборник задач по линейной алгебре. Проскуряков. И. В. М. 1961 год, стр. 118-119.
  • Решение матричных уравнений

    Обращение матриц. Решение матричных уравнений

    Таблица лучших: Решение матричных уравнений

    максимум из 2 баллов
    Место Имя Записано Баллы Результат
    Таблица загружается
    Нет данных

    M1498. Решение одной системы n-уравнений второй степени


    Условие

    Решите при каждом [latex]n>1[/latex] систему уравнений
    [latex]\left\{\begin{matrix}\alpha =\frac{\pi (2m+1)}{2(n+1)} x_{1}x_{n}=2, \\x_{2}(x_{n}-x_{1})=1, \\ …, \\x_{n-1}(x_{n}-x_{n-2})=1, \\x_{n}(x_{n}-x_{n-1})=1 \end{matrix}\right.[/latex]

    Решение

    При нескольких первых значениях [latex]n(n=2,3,4,5)[/latex] систему удается решить «в лоб»: положить [latex]x_{n}=z[/latex], можно вырвзить через [latex]z[/latex] последовательно [latex]x_{1},x_{2},…,[/latex], и наконец из последнего уравнения системы получить уравнение вида [latex]P_{n}(z)=0[/latex], где [latex]P_{n}[/latex] — многочлен. Например, при [latex]n=2[/latex] получим [latex]z=\pm \sqrt{3}[/latex], при [latex]n=3[/latex] — [latex]z=\pm \sqrt{2\pm \sqrt{2}}[/latex], при [latex]n=4[/latex] в ответе появляется корень из [latex]5[/latex]. Это может привести на мысль сделать тригонометрическую заменну переменной (и даже — какую именно). Положим [latex]x_{n}=2\cos \alpha[/latex]. Тогда [latex] x_{1}=\frac{1}{\cos \alpha }, x_{2}=\frac{1}{2\cos \alpha-\frac{1}{\cos \alpha } }=\frac{\cos \alpha }{\cos 2\alpha }:[/latex] и далее по индукции — предположив, что [latex] x_{k}=\frac{\cos (k-1)\alpha }{\cos k\alpha },[/latex] найдем [latex] x_{k+1}=\frac{1}{2\cos \alpha -\frac{\cos (k-1)\alpha }{\cos k\alpha }}=\frac{\cos k\alpha }{\cos (k+1)\alpha },[/latex] поскольку [latex]2\cos \alpha \cos \beta =\cos (\beta +\alpha )+\cos (\beta -\alpha )[/latex]. Последнее уравнение системы даст: [latex] x_{n}=\frac{\cos (n-1)\alpha }{\cos n\alpha }=2\cos \alpha [/latex] и преобразуется к виду [latex]\cos (n+1)\alpha =0[/latex]. Откуда [latex] \alpha =\frac{\pi (2m+1)}{2(n+1)};[/latex] при этом

    [latex] x_{k}=\frac{\cos (k-1)\alpha }{\cos k\alpha }(k=1,2, \ldots,n).(*)[/latex]

    Разные значения [latex]\cos \alpha[/latex] получаются при [latex]0< \frac{\pi (2m+1)}{2(n+1)} < \pi[/latex], т.е. при [latex]m=0,1, \ldots,n [/latex]. Однако не все они годятся: чтобы ни одно из чисел [latex]\cos k\alpha (k=1, \ldots,n)[/latex] не обращалось в [latex]0[/latex], необходимо и достаточно, чтобы [latex]2m+1[/latex] и [latex]n+1[/latex] не имели общего делителя, большего [latex]1[/latex] (если [latex]2m+1=dp[/latex], [latex]n+1=dp[/latex], [latex]d> 1[/latex], то [latex]p[/latex] — нечетно и [latex]\cos q\alpha =\cos \frac{\pi dpq}{2dq}=\cos \frac{p\pi}{2}=0[/latex]; легко доказать и обратное).

    Итак, к строчке [latex](*)[/latex], дающей ответ надо добавить условие: НОД [latex](2m+1,n+1)=1[/latex], [latex]0\leq m\leq n[/latex].

    Нужно еще показать, что найдены все решения. Из сказанного выше следует, что нет других решений, для которых [latex]\left | x_{n} \right |\leq 2[/latex]. Вот один из способов доказать, что решения с [latex]\left | x_{n} \right |> 2[/latex] быть не может.

    Обозначим [latex]\cosh\alpha =\frac{e^{\alpha }+e^{-\alpha }}{2}[/latex], где [latex]e[/latex] — основание натуральных логарифмов — что, впрочем, здесь не важно: нам понадобиться лишь, что [latex]e> 0[/latex] и что, как и для [latex]\cos \alpha [/latex], [latex]2\cosh\alpha \cosh\beta = \cosh(\alpha +\beta )+\cosh(\alpha -\beta ) [/latex]

    (Тем, кто знаком с комплексными числами, напомним, что [latex]\cos \alpha =\frac{e^{i\alpha }+e^{-i\alpha }}{2} [/latex], так что «гиперболический косинус» [latex]\cosh\alpha[/latex] — это просто [latex]\cos (i\alpha )[/latex].) Рассуждая так же, как и выше, — положив [latex]x_{n}=\pm 2\cosh\alpha [/latex], — найдем, что [latex]\cosh(n+1)\alpha =0[/latex]. Но функция [latex]\cosh[/latex] вообще не обращается в [latex]0[/latex] ([latex]\cosh\alpha \geq 1[/latex] при любом [latex]\alpha[/latex]), так что решений с [latex]\left | x_{n} \right |> 2[/latex] нет.

    К задаче M1498 Рассказ об этой задаче был бы неполон без объяснения, откуда возникла такая странная на первый взгляд система уравнений. Ее источник — геометрия. Построим равнобедренный треугольник [latex]ABC[/latex] с боковыми сторонами [latex]AB=BC=1[/latex] и углами при основании [latex]\alpha= \frac{\pi }{2(n+1)}[/latex]. Пусть [latex]K[/latex] — середина основания. Отметим на отрезке [latex]KC[/latex] точки [latex]M_{1},…,M_{n-1}[/latex] такие, что [latex]\angle M_{k-1}BM_{k}=\alpha [/latex] (здесь и ниже [latex]k=1,2,…,n[/latex]; [latex] M_{0}=K [/latex], [/latex] M_{k}=C[/latex], см. рисунок).

    Треугольники [latex]ABM_{k}[/latex] и [latex]CM_{k-1}B[/latex] подобны (их углы: [latex]\alpha [/latex], [latex](k+n)\alpha [/latex], [latex](n+1-k)\alpha [/latex]), так что [latex]AM_{k}\cdot M_{k-1}C =AB\cdot BC [/latex]. Положим [latex]x_{k}=AM_{k}[/latex], в частности, [latex]x_{n}=AC[/latex] тогда [latex]M_{k-1}C=x_{n}-x_{k-1}[/latex], поэтому [latex]x_{k}(x_{n}-x_{k-1})=1[/latex] и (поскольку [latex]AM_{0}=x_{0}/2[/latex]) [latex]x_{1}x_{n}=2[/latex]. Легко видеть, что (см. рисунок) [latex]AM_{k}=\cos (k-1)\alpha /\cos k\alpha [/latex], в частности, [latex]AM_{1}=1/\cos \alpha [/latex], [latex]AC=2\cos \alpha [/latex]. Таким образом, мы получим иллюстрацию «основного» решения системы с [latex]m=1[/latex].

    Заметим, что наш рисунок — фрагмент правильного [latex]2(n+1)[/latex]-угольника со стороной [latex]1[/latex]; [latex]x_{k}[/latex] — это кусочки, высекаемые на одной диагонали [latex]AC[/latex] диагоналями, выходящими из вершины [latex]B[/latex]. Решения системы, отвечающие значемиям [latex]m> 1[/latex], можно интерпретировать аналогичным образом как кусочки диагоналей ( или их продолжений ) правильной [latex]2(n+1)[/latex]-угольной звезды.

    Эта геометрическая интерпретация позволяет выяснить, при каких [latex]n[/latex] решения системы выражаются в квадратных радикалах ( через рациональные числа ): при тех, для которых можно построить правильный [latex](n+1)[/latex]-угольник ( а значит, и [latex]2(n+1)[/latex]-угольник ) циркулем и линейкой. Это — в точности те [latex]n[/latex], для которых число решений системы — степень двойки. Вот несколько первых значений [latex]n:2,3,4,5,7,9,11,14,15,16,19,23,…[/latex] ( см. статью А.Кириллова «О правильных многоугольниках, функции Эйлера и числах Ферма», «Квант» №6 за 1994 год).

    И.Васильев