Дискуссионный математический форумМатематический форум
Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

Часовой пояс: UTC + 3 часа [ Летнее время ]


Методы вычисления ранга матрицы

Методы вычисления ранга матрицы


Метод окаймляющих миноров


Находить ранг матрицы по определению — вычисляя миноры всех порядков — очень трудоемкая операция. Следующий алгоритм позволяет уменьшить число рассматриваемых миноров.


Пусть дана матрица A размеров m\times n. Будем говорить, что минор M_{i_1i_2\ldots i_ki_{k+1}}^{j_1j_2\ldots j_kj_{k+1}} (k+l)-ro порядка окаймляет (содержит в себе) минор M_{i_1i_2\ldots i_k}^{j_1j_2\ldots j_k} k-го порядка. При описании метода индексы выбранных строк и столбцов, в которых располагается минор, будем указывать, не упорядочивая их по возрастанию. При этом рассматриваемый минор и минор с упорядоченными индексами равны по абсолютной величине и, быть может, отличаются по знаку, но это для метода окаймляющих миноров не имеет никакого значения, поскольку нас интересует только ответ на вопрос: равен минор нулю или нет.


1. Выбираем строку i_1 и столбец j_1 так, чтобы минор 1-го порядка M_{j_1}^{i_1}=a_{i_1j_1} был не равен нулю. Если это возможно, то \operatorname{rg}A\geqslant1, иначе процесс завершается и \operatorname{rg}A=1.


2. Окаймляем минор M_{j_1}^{i_1}\ne0, добавляя к выбранным i_1-ой строке и j_1-му столбцу еще строку i_2\ne i_1 и столбец j_2\ne j_1 так, чтобы минор


M_{j_1j_2}^{i_1i_2}=\begin{vmatrix}a_{i_1j_1}&a_{i_1j_2}\\ a_{i_2j_1}& a_{i_2j_2} \end{vmatrix}\ne0.

Если это возможно, то \operatorname{rg}A\geqslant2, иначе процесс завершается и \operatorname{rg}A=2.


3. Окаймляем минор M_{j_1j_2}^{i_1i_2}\ne0, добавляя к выбранным ранее строкам и столбцам новую строку i_3 и новый столбец j_3 так, чтобы получить минор M_{j_1j_2j_3}^{i_1i_2i_3}\ne0. Если это удалось, то \operatorname{rg}A\geqslant3, иначе процесс завершается и \operatorname{rg}A=3.


Продолжаем процесс окаймления, пока он не завершится. Пусть найден минор r-го порядка M_{j_1\ldots j_r}^{i_1\ldots i_r}\ne0, т.е. \operatorname{rg}A\geqslant r. Однако, все миноры (r+l)-ro порядка, окаймляющие его, равны нулю M_{j_1j_2\ldots j_rj_{r+1}}^{i_1i_2\ldots i_ri_{r+1}}=0 или не существуют (при r=m или r=n). Тогда процесс завершается и \operatorname{rg}A=r.




Пример 3.6. Методом окаймляющих миноров найти ранги матриц


O=\begin{pmatrix}0&0\\0&0\end{pmatrix}\!,\quad A=\begin{pmatrix}3&9\\ 2&4 \end{pmatrix}\!,\quad B=\begin{pmatrix}0&2&3\\ 2&4&6\end{pmatrix}\!,\quad C=\begin{pmatrix} 1&0&2&1&3\\ 2&0&1&1&3\\ 3&0&3&2&5\end{pmatrix}\!.

Решение. Матрица O. 1. В этой матрице нет отличных от нуля миноров первого порядка, так как все ее элементы равны нулю. Поэтому \operatorname{rg}O=o.


Матрица A. 1. Выбираем первую строку (i_1=1) и первый столбец (j_1=1) матрицы A, на пересечении которых стоит ненулевой элемент a_{11}=3\ne0. Получили минор M_1^1=3\ne0. Следовательно, \operatorname{rg}A\geqslant1.


2. Добавляем к выбранным строке и столбцу еще одну строку i_2=2 и еще один столбец j_2=2. Получаем отличный от нуля минор второго порядка


M_{12}^{12}= \det{A}= \begin{vmatrix}3&9\\2&4\end{vmatrix}=-6\ne0 Следовательно, \operatorname{rg}A\geqslant2.

3. Поскольку исчерпаны все строки и все столбцы матрицы A, миноров, окаймляющих M_{12}^{12}\ne0, нет. Следовательно, \operatorname{rg}A=2.


Матрица B. 1. Выбираем первую строку и второй столбец матрицы B, на пересечении которых стоит ненулевой элемент b_{12}=2\ne0. Получили минор M_2^1=2\ne0. Следовательно, \operatorname{rg}B\geqslant1.


2. Добавляем к уже выбранным вторую строку и третий столбец. Получаем минор второго порядка M_{23}^{12}=\begin{vmatrix}2&3\\4&6\end{vmatrix}=0. Выбор оказался неудачным, так как получили нулевой минор. Вместо третьего столбца возьмем первый. Тогда получим отличный от нуля минор второго порядка M_{21}^{12}= \begin{vmatrix}2&0\\4&2\end{vmatrix}=4\ne0. Следовательно, \operatorname{rg}B\geqslant2.


3. Все строки матрицы B исчерпаны. Миноров третьего порядка нет. Поэтому \operatorname{rg}B=2.


Матрица C. 1. Выбираем первую строку (i_1=1) и первый столбец (j_1=1) матрицы C, на пересечении которых стоит ненулевой элемент a_{11}=1\ne0. Получили минор M_{1}^{1}=1\ne0. Следовательно, \operatorname{rg}C\geqslant1.


2. Добавляем к выбранным строке и столбцу еще одну строку i_2=2 и еще один столбец j_2=2. Получили минор второго порядка M_{12}^{12}= \begin{vmatrix}1&0\\2&0\end{vmatrix}. Выбор второго столбца оказался неудачным, так как получили минор, равный нулю. Возьмем вместо второго третий столбец (j_2=3). Получим минор M_{13}^{12}= \begin{vmatrix}1&2\\2&1\end{vmatrix}=-3\ne0. Следовательно, \operatorname{rg}C\geqslant2.


3. Окаймляем минор M_{13}^{12}\ne0. Имеется три окаймляющих минора


M_{134}^{123}= \begin{pmatrix}1&2&1\\ 2&1&1\\ 3&3&2\end{pmatrix}=0,\quad M_{135}^{123}= \begin{pmatrix}1&2&3\\ 2&1&2\\ 3&3&5\end{pmatrix}=0,\quad M_{132}^{123}= \begin{pmatrix}1&2&0\\ 2&1&0\\ 3&3&0\end{pmatrix}=0.

Три определителя равны нулю, так как третья строка равна сумме первых двух строк. Следовательно, нельзя найти отличный от нуля окаймляющий минор 3-го порядка, т.е. ранг матрицы C равен 2.


Замечание 3.4. Метод окаймляющих миноров позволяет уменьшить по сравнению с определением количество рассматриваемых миноров. Если в матрице размеров m\times n выбран минор r-го порядка (r<\min\{m,n\}), то количество окаймляющих его миноров (r+l)-ro порядка равно (m-r)(n-r), а общее количество миноров (r+1)-го порядка гораздо больше.




Метод Гаусса нахождения ранга матрицы


Пусть дана матрица A размеров m\times n. Для нахождения ее ранга нужно выполнить следующие действия.


1. Привести матрицу к ступенчатому виду (см. метод Гаусса).


2. В полученной матрице вычислить количество r ненулевых строк. Это число равно рангу матрицы A.


Замечания 3.5.


1. Обоснованием этого метода служит следствие 2 теоремы 3.4. Базисным минором в матрице ступенчатого вида (см. рис. 1.4) является минор


M=\begin{vmatrix}1&\ast&\cdots&\ast\\ 0&1&\cdots&\ast\\ \vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&1 \end{vmatrix},

составленный из столбцов, содержащих единичные элементы (в начале каждой "ступеньки"). Этот определитель треугольного вида отличен от нуля (равен 1), а любой его окаймляющий минор (если такой найдется) равен нулю, так как содержит нулевую строку.

2. Метод Гаусса для нахождения ранга произвольной матрицы наиболее экономичен, так как требует меньше вычислений, чем другие методы. Конечно, для матриц какого-либо специального вида (блочных, разреженных и т.п.) можно предложить более эффективные методы.




Пример 3.7. Методом Гаусса найти ранги матриц


O=\begin{pmatrix}0&0\\0&0\end{pmatrix}\!,\quad A=\begin{pmatrix}3&9\\2&4 \end{pmatrix}\!,\quad B=\begin{pmatrix}0&2&3\\2&4&6\end{pmatrix}\!,\quad C=\begin{pmatrix} 1&0&2&1&3\\ 2&0&1&1&2\\ 3&0&3&2&5\end{pmatrix}\!,\quad D=\begin{pmatrix}1&2&3\\ 0&0&0\\ 2&1&3\\ 1&1&2\\ 3&2&5 \end{pmatrix}\!.

Решение. Матрица O. 1. Нулевая матрица уже имеет ступенчатый вид (см. п.2 замечаний 1.8).


2. Количество ненулевых строк равно нулю. Следовательно, \operatorname{rg}O=0.


Матрица A. 1. Приводим матрицу A к ступенчатому виду (см. пример 1.29):


A\sim \begin{pmatrix}1&3\\0&1\end{pmatrix}\!.

2. В этой матрице две ненулевые строки. Следовательно, \operatorname{rg}A=2.


Матрица B. 1. Приводим матрицу B к ступенчатому виду (см. пример 1.29):


B\sim \begin{pmatrix}1&2&3\\0&1&1,\!5\end{pmatrix}\!.

В этой матрице две ненулевые строки. Следовательно, \operatorname{rg}B=2.


Матрица C. 1. Приводим матрицу C к ступенчатому виду. Взяв в качестве ведущего элемента a_{11}=1, делаем равными нулю остальные элементы первого столбца: ко второй строке прибавляем первую, умноженную на (-2), к третьей строке — первую, умноженную на (-3). Получаем матрицу


C=\begin{pmatrix}1&0&2&1&3\\ 2&0&1&1&2\\ 3&0&3&2&5\end{pmatrix}\sim \begin{pmatrix}1&0&2&1&3\\ 0&0&-3&-1&-4\\ 0&0&-3&-1&-4\end{pmatrix}\!,

У которой имеются две равные строки. По следствию 1 теоремы 3.3 одну из равных строк вычеркиваем:


C\sim \begin{pmatrix}1&0&2&1&3\\ 0&0&-3&-1&-4\end{pmatrix}\!.

Получили матрицу ступенчатого вида (см. п. 1 замечаний 1.8).


2. В этой матрице две ненулевые строки. Следовательно, \operatorname{rg}C=2.


Матрица D. 1. Приводим матрицу D к ступенчатому виду. Вычеркнув предварительно нулевую строку, берем в качестве ведущего элемента a_{11}=1, и делаем равными нулю остальные элементы первого столбца:


D\sim \begin{pmatrix}1&2&3\\ 2&1&3\\ 1&1&2\\ 3&2&5 \end{pmatrix}\sim \begin{pmatrix}1&2&3\\ 0&-3&-3\\ 0&-1&-1\\ 0&-4&-4\end{pmatrix}\!.

Последние три строки матрицы пропорциональны. По следствию 1 теоремы 3.3 две из них можно вычеркнуть:


D\sim \begin{pmatrix}1&2&3\\ 0&-1&-1\end{pmatrix}\!.

Получили матрицу ступенчатого вида (см. п. 1 замечаний 1.8).


2. В этой матрице две ненулевые строки. Следовательно, \operatorname{rg}D=2.


Заметим, что \operatorname{rg}C=\operatorname{rg}D, так как D=C^T (см. следствие 1 теоремы 3.4).




Пример 3.8. Даны матрицы


A=\begin{pmatrix}1&0\\0&0\end{pmatrix}\!,\quad B=\begin{pmatrix}0&0\\0&1 \end{pmatrix}\!,\quad C=\begin{pmatrix}0&1\\1&1\end{pmatrix}\!.

Найти ранги матриц: A+B;~A+C;~AB;~AC.


Решение. По определению имеем \operatorname{rg}A=1,~ \operatorname{rg}B=1,~ \operatorname{rg}C=2. Находим суммы и произведения данных матриц, а также их ранги:


A+B= \begin{pmatrix}1&0\\0&0\end{pmatrix}+ \begin{pmatrix}0&0\\0&1\end{pmatrix}= \begin{pmatrix}1&0\\0&1\end{pmatrix}= \operatorname{rg}(A+B)=2, то есть \operatorname{rg}(A+B)= \operatorname{rg}A+ \operatorname{rg}B;

A+C= \begin{pmatrix}1&0\\0&0\end{pmatrix}+ \begin{pmatrix}0&1\\1&1\end{pmatrix}= \begin{pmatrix}1&1\\1&1\end{pmatrix}= \operatorname{rg}(A+C)=1, то есть \operatorname{rg}(A+C)< \operatorname{rg}A+ \operatorname{rg}C;

A\cdot B= \begin{pmatrix}1&0\\0&0\end{pmatrix}\!\cdot \!\begin{pmatrix}0&0\\0&1\end{pmatrix}= \begin{pmatrix}0&0\\0&0\end{pmatrix}= \operatorname{rg}(A\cdot B)=0, то есть \operatorname{rg}(A\cdot B)< \min\{\operatorname{rg}A, \operatorname{rg}B\};

A\cdot C= \begin{pmatrix}1&0\\0&0\end{pmatrix}\!\cdot \!\begin{pmatrix}0&1\\1&1\end{pmatrix}= \begin{pmatrix}0&1\\0&0\end{pmatrix}= \operatorname{rg}(A\cdot C)=1, то есть \operatorname{rg}(A\cdot C)= \min\{\operatorname{rg}A, \operatorname{rg}C\};

Полученные результаты иллюстрируют справедливость теорем 3.5, 3.6.


См. также Ранг системы столбцов (строк)

Математический форум (помощь с решением задач, обсуждение вопросов по математике).
Кнопка "Поделиться"
Если заметили ошибку, опечатку или есть предложения, напишите в комментариях.

Часовой пояс: UTC + 3 часа [ Летнее время ]


Яндекс.Метрика

Copyright © 2010-2023 MathHelpPlanet.com. All rights reserved