Система Mathcad

.

Вопросы

1. Назначение системы MathCAD. Какие еще пакеты математических программ вы знаете?

2. Интерфейс пользователя в системе MathCAD.

3. Документ в системе MathCAD(заголовок, расширение при сохранении на диск, типы и расположение блоков, точка привязки блока, размеры блоков, сквозная передача данных в документе).

4. Объясните, что такое входной язык системы MathCAD и язык реализации системы MathCAD?

5. Перечислите основные объекты входного языка системы MathCAD. Расскажите об алфавите языка и о встроенных и пользовательских функциях системы MathCAD. Что такое определение функции и обращение к функции?

6. Константы и переменные в системе MathCAD? Как задаются типы данных в MathCAD? Что такое глобальное и локальное присваивание переменных в документе MathCAD? Как вставляется мнимая единица для комплексных чисел? Что такое ранжированная переменная и как она задается?

7. Как задаются массивы в MathCADе? Как можно добавлять строки и столбцы в готовые матрицы? Как удаляются строки и столбцы из матриц?

8. Перечислите операторы входного языка системы MathCAD?

9. Как осуществляется вывод результатов в системе MathCAD? Как можно настроить формат вывода результатов? Как осуществляется управление процессом вычислений в системе MathCAD?

10. Как работать с единицами измерений физических величин в системе MathCAD?

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

Лекция №12

Построение и форматирование графиков

Почти ни один технический расчет не обходится без построения графиков. Графики являются удобнейшим средством представления любой информации. MathCAD обладает обширным арсеналом средств для построения двух- и трехмерных графиков и диаграмм.

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

Графики функций одной переменной (двухмерные графики).

1. Графики в декартовой системе координат.

При построении графика функция представляется в виде набора точек на плоскости. Т.е. для построения графика перебирается определенное количество значений аргумента, и для каждого из них вычисляется значение функции. При упрощенном способе построения, когда диапазон аргумента задается автоматически, по оси абсцисс выбирается диапазон [-10, 10], а количество точек равно 100. Для того, чтобы иметь возможность при построении графика управлять количеством точек, аргумент надо задать как ранжированную переменную. При этом количество точек можно менять либо через шаг ранжированной переменной, либо непосредственно задавать количество точек N.

Изображение нескольких функций на одном графике. Для этого нужно ввести их в поле ввода функции через запятую (при этом каждая из них может зависеть от своего аргумента). Если необходимо, можно также через запятую ввести аргументы каждой функции в поле ввода аргументов. Максимально можно изобразить на одном графике 16 различных функций. Функции, построенные на одном графике, изображаются линиями различного цвета и типа. Цвета и типы линий можно изменить.

Ранжированные переменные дают возможность построить не только график функции, но и графически изобразить одномерный массив или вектор. В этом случае ранжированная переменная играет роль индекса - номер элемента вектора. Поэтому ее шаг обязательно должен равняться единице.

Масштаб.

Черные уголки на графике означают, что границы по осям выбираются автоматически. Границы можно изменить, для этого вместо старого значения надо ввести новое. Если необходимо, чтобы граница снова устанавливалась автоматически, то нужно щелкнуть на этой границе и удалить введенное в нее число, оставив незаполненным поле ввода. Если после этого установить курсор где-нибудь вне графика, то удаленная граница будет выбрана автоматически.

Есть другой более удобный способ изменения масштаба графика: установить курсор на поле графика, вызвать контекстное меню, выбрать команду Zoom, появится диалоговое окно; после чего на поле графика выбрать область, которую надо увеличить, а в диалоговом окне выбрать кнопку Zoom.

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

Трассировка.

Трассировка позволяет очень точно изучить строение графика. Для того, чтобы включить режим трассировки, щелкните в области графика правой кнопкой мыши и выберите в контекстном меню пункт Trace. В результате появится окно трассировки, а в окне графика вы увидите две пересекающиеся пунктирные линии. Перемещая указатель по графику, вы тем самым передвигаете точку пересечения линий трассировки. При этом координаты точки указываются с высокой точностью в окне трассировки в полях X-Value иY-Value. Нажатие кнопки Copy X или Copy Y копирует соответствующее число в буфер обмена. В дальнейшем его можно вставить в любое место документа или в маркер.

2. Графики в полярной системе координат (Polar Plot).

В полярной системе координат каждая точка задается двумя координатами: радиус-вектором и полярным углом. При изображении графика в полярной системе координат функция обычно задает зависимость радиус-вектора от угла. Если аргумент не был описан ранее как ранжированная переменная, то выбирается стандартный диапазон изменения угла - 2. Для полярного графика можно изменять границы только радиус-вектора. Поля для изменения этих границ расположены справа от графика. Здесь проявляется особенность полярных графиков в MatCAD - возможность установки для нижней границы радиус-вектора значений, отличных от нуля, т.е. можно получить график, для которого в начале координат радиус будет равняться, например, пяти.

3. График параметрически заданной функции.

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

Выше приведенные команды Zoom и Trace также справедливы для полярных графиков и для параметрически заданных.

Форматирование двухмерных графиков.

Вызов диалогового окна для форматирования графика осуществляется двойным щелчком мыши по полю графика.

Можно форматировать оси графика (X-Y Axes), линии на графике (Traces) - это форматирование понятно из контекста.

Для нанесения текстовых подписей используется вкладка Labels. Текстовые подписи форматируются стилем Variable, который используется для всех буквенных выражений в вычисляемых блоках. Для того чтобы этот стиль правильно отображал русские буквы, установите курсор на любом буквенном выражении в блоке формул или на графике. При этом в поле, отображающем текущий стиль, должно быть написано - Variable. Теперь выберите из раскрывающегося списка шрифтов шрифт, поддерживающий кириллицу, Times New Roman Cyr. После этого все введенные на русском языке подписи к графикам отобразятся правильно. Этот способ вставки подписей имеет несколько существенных недостатков: нельзя изменить их положение относительно осей, нельзя изменить форматирование шрифта конкретной подписи, нельзя создать подписи непосредственно в области осей координат для большей наглядности графика. Поэтому есть также другой, лишенный перечисленных недостатков, способ - ввести все необходимые подписи в виде текстовых блоков, отформатировать их нужным образом, а затем «перетащить» их в нужное место графика. После надо вызвать контекстное меню и выбрать команду «Bring to Front» (поместить перед графиком), иначе подпись будет невидимой, она закрыта полем графика.

Графики функций двух переменных (трехмерные графики).

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

- Surface Plot (График поверхности) - график в виде трехмерной поверхности;

- Contour Plot - (Контурный график или линии уровня) - также график поверхности, но спроектированный на одну из координатных плоскостей, а изменение функции по оставшейся оси изображается в виде линии уровня;

- Scatter Plot - отображение набора точек в трехмерном пространстве;

- Vector Field Plot - изображение двухмерного векторного поля;

- Bar Plot - трехмерная гистограмма;

- Patch Plot - также гистограмма, но отображающая только верхние грани столбиков.

Для каждого из этих графиков, кроме Patch Plot, есть отдельная кнопка на панели инструментов Graph. Также можно выбрать тип графика уже после создания в окне форматирования.

После выбора соответствующего типа графика появится шаблон с тремя координатными осями и местозаполнителем. В этот местозаполнитель следует ввести либо имя z функции двух переменных z(x,y) для быстрого построения трехмерного графика, либо имя матричной переменной Z, которая задает распределение данных ZX,Y на плоскости XY. Для графиков, основой которых служат матрицы, шкалу плоскости XY приходится задавать вручную. MathCAD просто рисует поверхность, точки в пространстве или линии уровня, основываясь на двумерной структуре этой матрицы. При быстром же построении графиков имеется возможность строить их в различном диапазоне аргументов, подобно двумерным графикам.

Вид трехмерного графика зависит от того, как он развернут относительно точки просмотра. Можно: 1) задавать угол обзора из окна форматирования мышью; 2) при нажатой клавише Ctrl можно мышью удалять или приближать объект к наблюдателю; 3) при нажатой клавише Shift можно инициировать анимацию («живую» картинку), для останова вращения надо щелкнуть кнопкой мыши.

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

График линий уровня. Линией уровня функции двух переменных X, Y называется геометрическое место точек в плоскости xOy, в которых функция принимает одно и то же значение. Линии уровня функции z = f(x,y) определяются уравнением f(x,y) = c, где с = const.

Показ цифр на графике можно включить вкладки окна форматирования трехмерного графика.

Векторное поле.

График векторного поля несколько отличается от остальных типов двухмерных графиков. Его смысл заключается в построении некоторого вектора в каждой точке плоскости XY. Чтобы задать вектор на плоскости, требуется два скалярных числа. Поэтому в MathCAD принято, что векторное поле задает комплексная матрица. Действительные части каждого ее элемента определяют проекцию вектора на ось X, а мнимые - на ось У.

Можно матрицу набрать вручную. А можно ее представить и как векторную функцию.

График параметрически заданной поверхности и график параметрически заданной кривой.

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

Кроме параметрической поверхности, можно также изобразить на трехмерном графике параметрическую кривую. Для этой цели подойдет трехмерный график типа Scatter Plot (точечный график). Параметрическая кривая в пространстве задается тремя функциями одной переменной. Для того чтобы получить приведенный выше график, на вкладке окна форматирования графика в группе Range был установлен диапазон изменения параметрической переменной от 0 (поле ввода Start) до 20 (поле End) и количество точек - 50.

Остальные типы графиков (точечные, гистограмма) пояснений не требуют.

Форматирование трехмерных графиков.

MathCAD представляет очень богатые возможности для форматирования трехмерных

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

Вкладка позволяет изменить тип графика (поле Plot), ориентацию графика в пространстве (Rotation - вращение, Tilt - наклон, Twist - поворот, Zoom - масштаб), менять стиль осей (Perimeter - периметр, Corner - углом, None - оси отсутствуют), можно задать изображение куба (флажок Show Box).

Вкладка содержит три вложенных вкладки , в которых задаются параметры для каждой из трех координатных осей.

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

На вкладке надо включить флажок . После этого станут активными все параметры на вкладке. Здесь можно включить и настроить до восьми различных источников света, для каждого из них есть отдельная вкладка и т.д. Но для того чтобы быстро создать освещение для графика, лучше воспользоваться заготовленными схемами освещения. Их можно выбрать из раскрывающегося списка (схема освещения).

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

Во вкладке задается заголовок графика, который можно расположить как сверху графика, так и снизу.

Во вкладке (плоскость заднего плана) задается показ проекций координатной сетки на три скрытые плоскости трехмерного графика.

Во вкладке можно изменить диапазон изменения координат X и Y. На этой вкладке есть группы элементов (диапазон 1) и (диапазон 2). В полях ввода Start и End этих групп можно указать диапазоны изменения обеих переменных. Также здесь можно указать тип координатной системы : - Декартова, - Сферическая, - Цилиндрическая.

Пример:

На этом рисунке представлен один и тот же график до и после форматирования. Для форматирования были применены следующие действия:

1. во вкладке установить переключатель и ;

2. во вкладке в группе установлен переключатель (это приведет к тому, что положение осей перестанет быть жестким, теперь они будут располагаться по периметру куба, ограничивающего график, и всегда на его передних гранях);

3. во вкладке на каждой из вкладок в группе установить флажок (отобразить линии), это приведет к отображению линий сетки, активизируя поле можно установить необходимый цвет линий сетки;

4. во вкладке включить флажок , из раскрывающегося списка выбрать схему освещения ;

5. на не отформатированном графике линии на графике изображены в виде так называемого каркаса (вкладка ), это очень загромождает график. Лучше линии отобразить в виде линий уровня по оси z. Для этого во вкладке в группе (параметры линий уровня) выбрать из раскрывающегося списка - линии уровня по z, далее в этой же группе следует установить флажок (отобразить линии); если плотность линий уровня оказалась слишком маленткой или слишком большой, то надо сбросить флажок и в поле ввода ввести необходимое число, в нашем случае это число 19.

Анимация

Во многих случаях самый зрелищный способ представления результатов математических расчетов - это анимация. MathCAD позволяет создавать анимационные ролики и сохранять их в видеофайлах.

Основной принцип анимации - покадровая анимация - это просто последовательность кадров, представляющих собой некоторый участок документа, который выделяется пользователем. Расчеты производятся обособленно для каждого кадра, причем формулы и графики, которые в нем содержатся, должны быть функцией от номера кадра. Номер кадра задается системной переменной , которая может принимать лишь натуральные значения. По умолчанию, если не включен режим подготовки анимации, =0.

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

Введите в документ необходимые выражения и графики, в которых участвует переменная номера кадра . Подготовьте часть документа, которую вы желаете сделать анимацией, таким образом, чтобы она находилась в поле вашего зрения на экране. В нашем примере подготовка сводится к определению функции f(x,t) = sin(x-t) и создании ее декартова графика f(x,FRAME).

Выполните команду Tools|Animation|Record.

В диалоговом окне Animate задайте номер первого кадра в поле From(От), номер последнего кадра в поле To (До) и скорость анимации в поле At (Скорость) в кадрах в секунду.

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

В диалоговом окне Animate нажмите кнопку Animate. После этого в окошке диалогового окна Animate будут появляться результаты расчетов выделенной области, сопровождающиеся выводом текущего значения переменной . По окончании этого процесса на экране появится окно проигрывателя анимации.

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

В случае если вид анимации вас устраивает, сохраните ее в виде видеофайла, нажав кнопку Save As в диалоговом окне Animate.

Закройте диалог Animate.

Сохраненный видеофайл можно использовать за пределами MathCAD. Если в Проводнике Windows дважды щелкнуть на имени этого файла, он будет загружен в проигрыватель видеофайлов Windows, и вы увидите его на экране компьютера. Таким образом, запуская видеофайлы обычным образом, можно устроить красочную презентацию результатов работы как на своем, так и на другом компьютере.

Ввод/вывод во внешние файлы.

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

READPRN (“file”) - чтение данных в матрицу из текстового файла;

WRITEPRN(“file”) - запись данных в текстовый файл;

APPENDPRN(“file”) - дозапись данных в существующий текстовый файл,

где file - путь к файлу.

Вопросы

1. Каким образом представляется функция при построении графика? График в декартовой системе координат: создание, задание аргумента по умолчанию, в виде ранжированной переменной, изображение нескольких функций на одном графике, работа с командами «масштаб» и «трассировка».

2. Построение графика в полярной системе координат и графика параметрически заданной функции.

3. Форматирование 2-х мерных графиков. Размещение надписей на поле графика.

4. Перечислите графики функций двух переменных, которые может строить MathCAD. Объясните, как строятся график поверхности? Каким образом на трехмерном графике можно изменить масштаб изображения, угол поворота? И как можно задать «живую» картинку?

5. Расскажите о контурных и векторных графиках в MathCAD. Объясните их.

6. Как строятся в MathCAD графики параметрически заданной поверхности и параметрически заданной кривой?

7. Расскажите о форматировании 3-хмерных графиков.

8. Что такое анимация и как она создается в системе MathCAD?

9. Как осуществляется ввод/вывод во внешние файлы в системе MathCAD?

Лекция №13

Задачи линейной алгебры

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

Операторы , вам уже знакомы. Напомним только, что оператор вычисляет только детерминант матрицы, а модуль вектора, который равен квадратному корню из суммы квадратов его элементов, вычисляется с помощью оператора , который расположен на панели Calculator. К сожалению, по внешнему виду они не отличаются.

При попытке вычислить модуль вектора с панели Matrix будет ошибочное состояние. Точно также будет ошибочное состояние при попытке вычислить детерминант матрицы с панели Calculator.

Рассмотрим неизвестные вам до сих пор операторы панели Matrix.

- 1

скалярное и векторное произведение векторов. Скалярное произведение векторов определяется как скаляр, равный сумме попарных произведений соответствующих элементов (идентичен обычному оператору умножения). Векторы должны иметь одинаковый размер. Для обозначения скалярного произведения используется символ «точка». Векторное произведение двух векторов u и v с углом между ними равно вектору с модулем , направленным перпендикулярно плоскости векторов u и v. Векторное произведение векторов применимо только для трехкомпонентных векторов. Обозначают векторное произведение символом х, который можно ввести нажатием кнопки на панели Matrix/

1

- 1

сумма элементов вектора.

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

Для решения задач линейной алгебры в MathCAD встроены матричные функции. Их можно разделить на три основные группы:

- функции определения (генерации) матриц и операции с блоками матриц;

- функции вычисления различных числовых характеристик матриц;

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

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

Первая группа:

1. matrix(m, n, f) - создает и заполняет матрицу размерности m x n, элемент которой, расположенный в i-ой строке и j-ом столбце равен значению f(i, j) функции f(x, y);

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

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

4. augment(A, B) -объединяет матрицы A и B; матрица B располагается справа от матрицы A, при этом матрицы должны иметь одинаковое число строк;

5. stack(A, B) - объединяет матрицы A и B, матрица В располагается внизу под матрицей А, при этом матрицы должны иметь одинаковое число столбцов;

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

Вторая группа:

1. last(v) - вычисляет номер последнего элемента вектора V;

2. length(v) - вычисляет количество элементов вектора V;

3. min(v), max(v) - вычисляет минимальное и максимальное значения вектора V;

4. Re(v) - создает вектор из реальных частей комплексных элементов вектора V;

5. Im(v) - создает вектор из мнимых частей комплексных элементов вектора V;

6. sort(V) - сортировка элементов вектора V по возрастанию;

7. reverse (sort(v)) - сортировка элементов вектора V по убыванию;

8. csort (A,n) - сортировка элементов n - го столбца матрицы А по возрастанию (перестановкой строк);

9. rsort (A,n) - сортировка элементов n - ой строки матрица А по возрастанию (перестановкой столбцов);

10. rows(A) - вычисляет число строк в матрице А;

11. cols(A) - вычисляет число столбцов в матрице А;

12. max(A), min(A) - определяет максимальное и минимальное значения матрицы А;

13. tr(A) - вычисляет след квадратной матрицы А (след матрицы равен сумме ее диагональных элементов по главной диагонали);

14. mean(A) - среднее значение элементов матрица А.

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

Третья группа:

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

2. rank(A) - вычисляет ранг матрицы А (количество линейно-независимых строк или это число ненулевых строк ступенчатой матрицы rref(A));

3. eigenvals(A) - вычисление собственных значений квадратной матрицы А;

4. eigenvecs (A) - вычисление собственных векторов квадратной матрицы А, значением функции является матрица, столбцы которой есть собственные векторы матрицы А, причем порядок следования векторов отвечает порядку следования собственных значений, вычисленных с помощью функции eigenvals(A);

5. eigenvec(A,e) - вычисление собственного вектора матрицы А, отвечающего собственному значению e;

6. normi(A) - max - норма, или - норма (infinity norm). в линейной алгебре используются различные матричные нормы, которые ставят в соответствие матрице некоторую скалярную числовую характеристику;

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

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

В самом простом виде задача на собственные значения матрицы формулируется следующим образом: требуется найти такие значения , чтобы матричное уравнение имело решение. В таком случае число называют собственным числом матрицы А, а n- компонентный вектор Х, приводящий уравнение с заданным в тождество - собственным вектором. В вышеприведенном примере собственные вектора матрицы А получены в матрице MS. Проверка проведена для первого столбца матрицы MS и соответствующего ему собственного числа 0=5.439.

Решение систем линейных алгебраических уравнений. Этот вопрос является центральным в вычислительной линейной алгебре.

В математике рассматриваются системы линейных уравнений двух видов - однородные и неоднородные.

Неоднородная система уравнений в матричном виде записывается следующим образом: . Здесь А - матрица коэффициентов системы, В - вектор свободных членов, Х - вектор неизвестных системы.

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

- метод обратной матрицы, он вам уже известен;

- метод исключений Гаусса;

- метод Крамера.

Неоднородная система линейных уравнений в случае равенства ее определителя нулю имеет множество решений, если ранг матрицы системы равен рангу расширенной матрицы системы, либо не имеет решения, если это условие не выполняется. Решить такие системы в MatCADe можно методом Гаусса.

В выше приведенном примере получили систему из трех уравнений с пятью неизвестными, поэтому решение системы будет иметь два свободных параметра (x4, x5).

Однородная система линейных алгебраических уравнений может быть представлена в виде , т.е. правая часть уравнения представляет вектор из нулевых элементов. Как известно, для того чтобы однородная система линейных уравнений имела решение, определитель соответствующей матрицы должен равняться нулю. Это означает, что количество независимых уравнений в системе (т.е. ранг матрицы) меньше, чем количество неизвестных (т.е. порядок матрицы): rank(A) < n. Но вначале нужно выделить в системе эти самые независимые уравнения. Это делается с помощью функции rref, которая с помощью метода исключений Гаусса приводит матрицу к ступенчатому виду.

Приближенное решение систем линейных уравнений. Точные методы решения линейных систем применяют для решения линейных систем относительно небольшой размерности. Для решения систем большой размерности используют итерационные методы. Итерационные методы хороши для систем с разреженными матрицами. Простейший итерационный метод - метод простых итераций.

Метод состоит в том, что система уравнений преобразуется к виду и ее решение вычисляется как предел последовательности

В качестве условия окончания итерационного процесса можно взять условие

,

где - заданная погрешность приближенного решения.

Для сходимости метода простых итераций достаточно, чтобы выполнялось условие

normi(A) < 1.

Рассмотрим следующий пример:

Вопросы

1. Поясните работу команд панели Matrix - скалярное и векторное произведение, детерминант матрицы, сумма элементов вектора, операция векторизации.

2. Перечислите три основные группы матричных функций. Расскажите о матричных функциях, возвращающих числовые характеристики. Приведите примеры.

3. Матричные функции, реализующие генерацию матриц и операции работы с блоками матриц.

4. Перечислите матричные функции, реализующие численные алгоритмы решения задач линейной алгебры. Объясните, как работают функции rref и rank.

5. Какие функции вычисляют собственные вектора и собственные числа квадратной матрицы?

6. Решение в системе MathCAD неоднородных систем линейных уравнений, когда определитель матрицы не равен нулю. Три способа.

Страницы: 1, 2, 3, 4