Mathcad №2

Системы линейных алгебраических уравнений

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

В задачах линейной алгебры практически всегда возникает не­обходимость выполнять различные операции с матрицами. Предва­рительно матрицу нужно определить и ввести в рабочий документ Mathcad. Для того чтобы определить матрицу, необходимо ввести с клавиатуры имя матрицы и знак присваивания. Затем щелкнуть по кнопке в панели математических инструментов, чтобы открыть панель операций с ма­трицами и векторами (рис.1).

Рис.1

Далее нужно щелчком по кнопке открыть окно диалога Вставить матрицу, определить число строк и число столбцов, а затем щелкнуть по кнопке ОК.

В рабочем документе справа от знака присваивания откроется поле ввода матрицы с помеченными позициями для ввода элементов.

Функции, предназначенные для решения задач линейной алгебры, можно разделить на три группы: функции определения матриц и операции с блоками матриц, функции вычисления различных числовых характеристик матриц и функции, реализующие численные алгоритмы решения задач линейной алгебры. Опишем наиболее часто используемые функции.

Функции определения матриц и операции с блоками ма­триц:

diag(v) - создает диагональную матрицу, элементы главной диагонали которой хранятся в векторе v;

identity(n) — создает единичную матрицу порядка n;

augment(A, В) — формирует матрицу, в первых столбцах которой содержится матрица А, а в последних — матрица В (матрицы А и B должны иметь одинаковое число строк);

stack(A, В) — формирует матрицу, в первых строках которой содержится матрица А, а в последних — матрица В (матрицы А и В должны иметь одинаковое число столбцов);

submatrix(A, ir, jr, ic, jc) — формирует матрицу, которая является блоком матрицы А, расположенным в строках с ir no jr и в столбцах с ic по jc, ir < jr, ic < jc.

Номер первой строки (столбца) матрицы или первой компоненты вектора хранится в Mathcad в переменной ORIGIN. По умолчанию в Mathcad координаты векторов, столбцы и строки матрицы нумеру­ются начиная с 0 (ORIGIN := 0). Поскольку в математической записи чаще используется нумерация с 1, то в дальнейшем перед началом работы с матрицами будем определять значение переменной ORIGIN равным 1, т.е. будем прежде всего выполнять команду ORIGIN := 1.

Фрагмент рабочего документа Mathcad, содержащий примеры исполнения перечисленных выше функций, приведен ниже.

Функции вычисления различных числовых характеристик матриц:

last(v) — вычисление номера последней компоненты вектора v;

length (v) — вычисление количества компонент вектора v;

rows(A) — вычисление числа строк в матрице А;

cols(A) — вычисление числа столбцов в матрице А;

mах(А) — вычисление наибольшего элемента в матрице А;

rank(A) — вычисление ранга матрицы А;

norm1(A), norm2(A), norme(A), normi(A) — вычисление норм квадратной матрицы А.

Примеры использования названных функций для матриц, определенных в приведенном ранее рабочем документе.

last(v) = 4 length(v) = 4

rows(D)=3 cols(D)=7 max(D)=5 rank(A) =2

norm1(B)=4

Функции, реализующие численные алгоритмы решении задач линейной алгебры:

rref(A) - приведение матрицы к ступенчатому виду с единичным базисным минором (выполняются элементарные операции со строками матрицы);

lsolve(A, b) — решение системы линейных алгебраических уравнений Ax=b.

Матричная форма записи линейных систем.

Рассмотрим систему n линейных алгебраических уравнений относительно n неизвестных x₁, x₂,…, x_n:

(1)

В соответствии с правилами умножения матриц рассмотренная система линейных уравнений может быть записана в матричном виде Ax = b, где

, ,

Матрица А называется матрицей системы; матрица-столбец b, элементами которой являются правые части уравнений системы, называется матрицей правой части. Матрица-столбец x, элементы которой – искомые неизвестные, называется решением системы.

Метод Крамера

Справедливо следующее утверждение. Если определитель  = det A матрицы системы Ax = b отличен от нуля, то система имеет единственное решение x₁, x₂,…, x_n, определяемое формулами Крамера , где _i – определитель матрицы n-го порядка, полученной из матрицы системы заменой i-го столбца столбцом правых частей.