Матричные вычисления.

MathCAD поддерживает два вида массивов – одномерные (векторы) и двумерные (матрицы). Элементами массива могут быть числа, строки, математические выражения и даже другие массивы. Основные операции для работы с векторами и матрицами собраны на панели математических инструментов Matrix. Учтите, что элементы матрицы по умолчанию нумеруются с 0, если хотите, чтобы элементы матрицы нумеровались с 1, нужно в начале документа ввести ORIGIN: =1.

1. Способы задания матрицы.

Матрицу можно целиком ввести с клавиатуры, либо с помощью функциональной зависимости элемента массива от его индексов.

1.1 Введение элементов матрицы с клавиатуры

Для того чтобы ввести элементы матрицы с клавиатуры, выберите команду меню Insert → Matrix или на панели инструментов Matrix щелкните на кнопке Matrix or Vector. В открывшемся диалоговом окне введите количество строк (Rows) и столбцов (Columns). При нажатии клавиши OK появится шаблон матрицы, в который можно вводить ее элементы.

Это же диалоговое окно (Insert Matrix или Вставить Матрицу) позволяет добавлять и удалять несколько строк и столбцов в уже имеющейся матрице.

Для того, чтобы добавить строки и столбцы в матрицу, установите курсор на элемент матрицы, справа от которого вы хотите вставить столбцы и ниже которого вы хотите вставить строки. Введите количество вставляемых строк и столбцов и щелкните на кнопке Insert (Вставить).

Например:

Для того, чтобы удалить строки и столбцы из матрицы, установите курсор на на элемент матрицы, справа от которого вы хотите удалить столбцы и ниже которого вы хотите удалить строки. Введите количество удаляемых строк и столбцов и щелкните на кнопке Delete (Удалить). ПРИ ЭТОМ СТРОКА И СТОЛБЕЦ, НА ПЕРЕСЕЧЕНИИ КОТОРЫХ СТОЯЛ УСТАНОВЛЕННЫЙ КУРСОР ТОЖЕ БУДУТ УДАЛЕНЫ.

Для доступа к элементу матрицы нужно указать номер строки и столбца нужного элемента в виде индексов.

Если ввести элемент матрицы, которого не существует, то матрица автоматически будет увеличена до размера, вмещающего введенный элемент. Например:

Кроме доступа к отдельным элементам матрицы MathCAD дает возможность выводить и изменять отдельный столбец или строку матрицы. Для того, чтобы обратиться к столбцу матрицы введите ее имя, щелкните на кнопке с изображением M^{< >} на панели инструментов Matrix и в появившемся поле введите номер столбца.

Для выделения аналогичным образом строки, матрицу нужно предварительно транспонировать. Например:

1.2 Задание матрицы, элементы которой являются функциями индексов.

Например, для того, чтобы задать матрицу A размером 4х 5, каждый элемент которой равен сумме номера строки и удвоенного номера столбца, то нужно

задать переменные диапазона i: =1..4, j: =1..5 и задать формулу A_{i, j}: =i+2k.

Можно задать такую матрицу без использования переменных диапазона, с помощью специальной функции matrix(4, 5, f), описав заранее функцию f. Например:

1. Основные операторы и функции для работы с массивами.

К основным матричным операциям можно отнести:

- арифметические операции: поэлементное сложение (+), вычитание (-), матричное умножение (*), которые вводятся с помощью клавиш

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

Некоторые матричные операции заданы в MathCAD в виде функций (для их вставки следует выбрать команду меню Insert → Function→ Vector and Matrix). Рассмотрим некоторые из них:

• identity(n) – возвращает единичную матрицу размера nxn;
• diag (v) – возвращает диагональную матрицу, у которой на диагонали расположены элементы вектора v;
• rank (M) – возвращает ранг матрицы M;
• tr (M) – возвращает след (сумму диагональных элементов) матрицы M;
• norme (M) – возвращает евклидову норму матрицы M (корень из суммы квадратов всех элементов).

1. Оператор векторизации.

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

1. Объединение матриц и выделение подматрицы.

В MathCAD можно присоединять матрицы СПРАВА с помощью функции augment (аргументами может быть любое количество матриц с одинаковым количеством строк) и СНИЗУ с помощью функции stack (аргументами должны быть матрицы с одинаковым количеством столбцов). Например:

Выделение подматрицы осуществляется с помощью функции submatrix(M, imin, imax, jmin, jmax), где M – исходная матрица, imin, imax – номера первой и последней строк исходной матрицы, входящих в выделяемый блок, jmin, jmax – номера первого и последнего столбцов исходной матрицы, входящих в выделяемый блок.

1. Собственные вектора и собственные числа.

Для поиска собственных векторов и собственных чисел в MathCAD предусмотрены следующие функции:

• eigenvals(A) – собственные числа матрицы A;
• eigenvecs(A) – собственные векторы матрицы A;
• eigenvec (A, l) – cобственный вектор, соответствующий собственному значению l.

Например:

1. Символьные операции с матрицами.

Для того, чтобы производить преобразования матриц в символьном виде используется оператор символьного вычисления →. Например: