Большинство разработчиков с трудом представляет, как его синтаксис, так и возможности. Всё дело в том, что язык напрямую связан с популярным программным продуктом, стоимость которого может достигать потрясающих воображение значений. Итак, главный вопрос: так ли хорош непосредственно язык Matlab? И может ли он быть полезен именно вам.
Начнём не со стандартного экскурса в историю и обсуждения плюсов и минусов языка, а с программной среды MATLAB/Simulink - единственного места, где герой этого текста может быть полезен. Просто представьте себе графический редактор, в котором вы сможете реализовать любую свою задумку, не имея за плечами нескольких лет опыта и соответствующего образования. И создав один раз схему взаимодействия инструментов, получить качественный скрипт для многократного использования.
MATLAB - именно такой редактор в мире данных. Область его применения безгранично широка: IoT, финансы, медицина, космос, автоматика, робототехника, беспроводные системы и многое-многое другое. В общем почти неограниченные возможности по сбору и визуализации данных, а также прогнозированию, но только если есть возможность купить соответствующий пакет.
Что касается цены, то в верхней границы почти нет, а вот нижняя находится в район 99$. Чтобы урвать столь мощный продукт за относительно небольшие деньги, вам необходимо быть студентом ВУЗа. И конечно же вы получите довольно ограниченный продукт.
Язык MATLAB - инструмент, обеспечивающий взаимодействие оператора (часто даже не программиста) со всеми доступными возможностями анализа, сбора и представления данных. У него есть очевидные плюсы и минусы, свойственные языку живущему в замкнутой экосистеме.
Недостатки:
Медленный и перегруженный операторами, командами, функциями язык, основной целью которого является улучшение визуального восприятия.
Узконаправленный. Нет никакой больше программной платформы, где бы MATLAB был полезен.
Дороговизна ПО. Если вы не студент - либо готовьтесь опустошить карманы или перейти границу закона. И даже если студент - цена приличная.
Невысокий спрос. Несмотря на большой интерес к MATLAB практически во всех сферах, фактически и легально его используют лишь немногие.
Достоинства:
Язык легок для изучения, обладает простым и понятным синтаксисом.
Огромные возможности. Но это скорее преимущество всего продукта в целом.
Частые обновления, как правило заметные положительные преобразования происходят не реже пары раз в год.
Программная среда позволяет преобразовывать его в “быстрый” код на С, С++.
Разумеется, MATLAB нужен далеко не всем. Несмотря на широчайшую область применения, трудно представить, что рядовому разработчику приложений может понадобиться знание этого языка. MATLAB крайне полезен в областях, требующих особой надёжности при обработке данных, например, в системах автопилота в автомобилях или бортовых электронных системах самолёта.
То есть если вы не очень программист, но так или иначе ваша профессия связана с необходимостью программной обработки данных, то продукт MATLAB/Simulink с соответствующим языком способны сильно упростить ваши каждодневные задачи.
Завершаем обзор языка как всегда списком учебной литературы. Само-собой среди них вы не отыщите книг исключительно по языку, но от этого восприятие языка будет только проще:
А у вас есть опыт работы с MATLAB? И какой?
Для тех, кто хочет стать программистом - .
Те, кто имеет дело с высшей математикой, прекрасно знают, с какими математическими «чудовищами» иногда приходится сталкиваться. Например, на вычисление какого-нибудь гигантского тройного интеграла можно потратить настоящую уйму времени, душевных сил и не восстанавливающихся нервных клеток. Конечно, это очень интересно, бросить вызов интегралу, и взять его. Но, что делать, если вместо этого интеграл грозиться взять Вас? Или, что еще хуже, кубический трехчлен вышел из-под контроля и разбушевался? Такого и врагу не пожелаешь.
Раньше вариантов было всего два: плюнуть на все и пойти гулять или вступить в многочасовую схватку с интегралом. Ну, кому многочасовую, кому многоминутную – кто как учился. Но суть не в этом. Двадцатый век и неумолимо движущийся прогресс предлагают нам третий способ, а именно позволяют взять самый сложный интеграл «по-быстрому». То же самое касается решения всевозможных уравнений, построения графиков функций в виде кубических гиперболоидов и т.д.
Для таких неординарных, но периодически случающихся среди студентов ситуаций существует мощное математическое оружие. Встречайте, кто еще не знает – пакет программ MATLAB.
Матлаб и решит уравнение, и аппроксимирует, и построит график функции. Понимаете, что это значит, друзья?
Это значит, что – один из мощнейших на сегодняшний день пакетов обработки данных. Название расшифровывается как Matrix Laboratory. Матричная Лаборатория, если по-русски. Возможности программы покрывают практически все области математики. Так, пользуясь матлабом, Вы сможете:
Помимо всего этого возможности MATLAB позволяют визуализировать данные вплоть до построения трехмерных графиков и создания анимированных роликов.
Наше описание матлаб, конечно, далеко не полное. Помимо предусмотренных производителем возможностей и функций существует огромное количество инструментов матлаб, написанных просто энтузиастами или другими компаниями.
А еще – это язык программирования, используемый непосредственно при работе с программой. Не будем вдаваться в подробности, скажем только, что программы, написанные на языке MATLAB, бывают двух видов: функции и скрипты.
Основной рабочий файл программы – М-файл. Это бесконечный текстовый файл, и именно в нем происходит непосредственно программирование вычислений. Кстати, пусть Вас не пугает это слово – для того, чтобы работать в MATLAB, вовсе не нужно быть профессиональным программистом.
М-файлы делятся на
Для того чтобы наглядно показать, как происходит работа в MATLAB, приведем ниже пример создания функции в матлабе. Данная функция будет вычислять среднее значение вектора.
f
unction y = average (x)
% AVERAGE Среднее значение элементов вектора.
% AVERAGE(X), где X - вектор. Вычисляет среднее значение элементов вектора.
% Если входной аргумент не является вектором, генерируется ошибка.
= size(x);
if (~((m == 1) | (n == 1)) | (m == 1 & n == 1))
error("Входной массив должен быть вектором’)
end
y =sum(x)/length(x); % Собственно вычисление
Строка определения функции сообщает системе MATLAB, что файл является М-функцией, а также определяет список входных аргументов. Так, строка определения функции average имеет вид:
function y = average(x)
Где:
Итак, чтобы написать функцию в матлабе, необходимо помнить, что каждая функция в системе MATLAB содержит строку определения функции, подобную приведенной.
Безусловно, такой мощный пакет нужен не только для того, чтобы облегчить жизнь студентам. В настоящее время MATLAB, с одной стороны, очень популярен среди специалистов многих научных и инженерных отраслей. С другой стороны, возможность работы с большими матрицами делает MATLAB незаменимым инструментом финансовых аналитиков, позволяющим решить намного больше задач, чем, к примеру, известный всем Excel. Подробнее о том, вы можете прочитать в обзорной статье.
Какие есть трудности в работе с MATLAB? Трудность, пожалуй, всего одна. Но фундаментальная. Чтобы полностью раскрыть возможности MATLAB и с легкостью решать встающие перед Вами задачи, придется попотеть и сначала разобраться с самим матлабом (как создать файл, как создать функцию и др.). А это не так просто, ибо мощность и широкие возможности требуют жертв.
При всем желании нельзя сказать, что MATLAB – простая программа. Тем не менее, надеемся, все вышеперечисленное будет достаточным аргументом для того, чтобы взяться за ее освоение.
И напоследок. Если Вы не знаете, почему все в Вашей жизни пошло так, а не иначе, спросите об этом у матлаба. Просто наберите в командной строке “why” (почему). Он ответит. Попробуйте!
Теперь вы знаете возможности Матлаб. В области образования MATLAB часто используется в преподавании численных методов и линейной алгебры. Многим студентам не обойтись без него при обработке результатов эксперимента, проведенного в ходе лабораторной работы. Для быстрого и качественного освоения основ работы с MATLAB Вы всегда можете обратиться к , в любой момент готовым ответить на любой Ваш вопрос.
Инструкция
В среде MATLAB есть несколько режимов работы. Самый простой – это ввод команд непосредственно в окно команд (Command Window
).
Если оно не видно в интерфейсе программы, значит нужно его открыть. Найти окно команд можно через меню Desktop
-> Command Window
.
Давайте для примера введём в это окно последовательно друг за другом команды "x = ; y = sqrt(x); plot(y);", и нажмём клавишу "Ввод" (Enter
). Программа моментально создаст переменные X, создаст переменную Y и посчитает её значения по заданной функции, а затем построит её график.
Стрелками клавиатуры "Вверх" и "Вниз" в окне команд мы можем переключаться между всеми введёнными командами, тут же изменять их при необходимости, и по нажатию Enter снова отправлять среде MATLAB на исполнение.
Удобно? Безусловно. И главное - очень быстро. Все эти действия занимают несколько секунд.
Но что если нужна более сложная организация команд? Если нужно циклическое исполнение каких-то команд? Вводить команды вручную по одной, а потом долго искать их в истории может быть довольно утомительным делом.
Чтобы упростить жизнь учёному, инженеру или студенту, служит окно редактора (Editor
). Давайте откроем окно редактора через меню Desktop
-> Editor
.
Здесь можно создавать новые переменные, строить графики, писать программы (скрипты), создавать компоненты для обмена с другими средами, создавать приложения с пользовательским интерфейсом (GUI), а также редактировать уже имеющиеся. Но нас в данный момент интересует написание программы, содержащей функции для повторного использования в будущем. Поэтому идём в меню File
и выбираем New
-> M-File
.
В поле редактора напишем простую программу, но чуть-чуть усложним её:
function draw_plot(x)
y = log(x); % Задаём первую функцию
subplot(1, 2, 1), plot(x, y); % Строим первый график
y = sqrt(x); % Задаём вторую функцию
subplot(1, 2, 2), plot(x, y); % Строим второй график
Мы добавили вторую функцию и будем выводить сразу два графика рядом друг с другом. Знаком процента обозначаются в среде MATLAB комментарии.
Не забудем сохранить программу. Стандартное расширение файла с программой Матлаб – *.m
.
Теперь закройте редактор и окно с графиком, который мы построили ранее.
Переходим обратно в окно команд.
Можно очистить историю команд, чтобы лишняя информация нас не отвлекала. Для этого кликните правой кнопкой мыши на поле ввода команд и в открывшемся контекстном меню выберите пункт Clear Command Window
.
Переменная X у нас осталась после предыдущего эксперимента, мы её не изменяли и не удаляли. Поэтому в окно команд можно сразу ввести:
draw_plot(x);
Вы увидите, что MATLAB прочитает нашу функцию из файла и выполнит её, нарисовав график.
В среде MATLAB есть несколько режимов работы. Самый простой - это ввод команд непосредственно в окно команд (Command Window ).
Окно команд MATLAB
Если оно не видно в интерфейсе программы, откроем его. Найти окно команд можно через меню Desktop Command Window .
Давайте для примера введём в это окно последовательно друг за другом команды
X = ; y = sqrt(x); plot(y);
и нажмём клавишу «Ввод» (Enter ). Программа моментально создаст переменную X, создаст переменную Y и посчитает её значения по заданной функции, а затем построит её график.
Стрелками клавиатуры вверх и вниз в окне команд мы можем переключаться между введёнными командами, тут же изменять их, а по нажатию Enter отправлять среде MATLAB на исполнение. Стрелками влево и вправо можно перемещаться по введённой команде и редактировать её. Если в конце команды стоит точка с запятой, то результат будет посчитан, но не будет выведен в окно команд; в противном случае результат выполнения команды будет отображён тут же. По любой функции в среде MATLAB есть подробная встроенная справка. Например, чтобы получить справку по команде plot , выделите эту команду, нажмите на неё правой кнопкой мыши, и в открывшемся контекстном меню выберите пункт Help on Selection или нажмите клавишу F1 .
Удобно? Безусловно. И главное - очень быстро. Все эти действия занимают несколько секунд.
Но что если нужна более сложная организация команд? Если нужно циклическое исполнение каких-то команд? Вводить команды вручную по одной, а потом долго искать их в истории может быть довольно утомительным делом.
Чтобы упростить жизнь учёному, инженеру или студенту, служит окно редактора (Editor ). Давайте откроем окно редактора через меню Desktop Editor .
В окне редактора можно создавать новые переменные, строить графики, писать программы (скрипты), создавать компоненты для обмена с другими средами, создавать приложения с пользовательским интерфейсом (GUI), а также редактировать имеющиеся.
Нас в данный момент интересует написание программы, содержащей функции для повторного использования в будущем. Поэтому идём в меню File редактора и выбираем New M-File .
М-файлами в среде МАТЛАБ называются файлы, содержащие текст программ (скрипты) или определённые пользователем функции.
Давайте напишем в редакторе простую функцию draw_plot :
function draw_plot(x) % Задаём первую функцию: y = log(x); % Строим первый график: subplot(1, 2, 1), plot(x, y); % Задаём вторую функцию: y = sqrt(x); % Строим второй график: subplot(1, 2, 2), plot(x, y);Переходим обратно в окно команд.
Можно очистить историю команд, чтобы лишняя информация нас не отвлекала. Для этого кликните правой кнопкой мыши на поле ввода команд и в открывшемся контекстном меню выберите пункт Clear Command Window .
Переменная X у нас осталась после предыдущего эксперимента, мы её не изменяли и не удаляли. Поэтому в окно команд можно сразу ввести:
Draw_plot(x);
Вы увидите, что MATLAB прочитает нашу функцию из файла и выполнит её, нарисовав график.
Если MATLAB при выполнении программы выдаст сообщение, Undefined function or method "draw_plot" for input arguments of type "double". (т.е. вызывается неизвестная функция), нажмите в окне редактора зелёную кнопку со стрелкой (Run ) или через меню редактора: Debug Run draw_plot.m . MATLAB сообщит, что директория, в которой находится файл с нашей программой (draw_plot.m), не является рабочей директорией. Нажмите в диалоговом окне кнопку Add to Path , чтобы пакет MATLAB добавил директорию к рабочему пути и мог использовать наш M-файл. После этого программа должна запуститься нормально.
Среда MATLAB включает интерпретатор команд на языке высокого уровня, графическую систему, пакеты расширений и реализована на языке C. Вся работа организуется через командное окно (Command Window), которое появляется при запуске программы matlab.exe. В процессе работы данные располагаются в памяти (Workspace), для изображения кривых, поверхностей и других графиков создаются графические окна.
В командном окне в режиме диалога проводятся вычисления. Пользователь вводит команды или запускает на выполнение файлы с текстами на языке MATLAB. Интерпретатор обрабатывает введенное и выдает результаты: числовые и строковые данные, предупреждения и сообщения об ошибках. Строка ввода помечена знаком >>. В командном окне показываются вводимые с клавиатуры числа, переменные, а также результаты вычислений. Имена переменных должны начинаться с буквы. Знак = соответствует операции присваивания. Нажатие клавиши Enter заставляет систему вычислить выражение и показать результат. Наберите с клавиатуры в строке ввода:
Нажмите клавишу Enter, на экране в зоне просмотра появится результат вычисления:
Все значения переменных, вычисленные в течение текущего сеанса работы, сохраняются в специально зарезервированной области памяти компьютера, называемой рабочим пространством системы MATLAB (Workspace). Командой clc можно стереть содержимое командного окна, однако это не затронет содержимого рабочего пространства. Когда исчезает необходимость в хранении ряда переменных в текущем сеансе работы, их можно стереть из памяти компьютера командой clear или clear(имя1, имя2, ...). Первая команда удаляет из памяти все переменные, а вторая - переменные с именами имя1 и имя2. Командой who можно вывести список всех переменных, входящих в данный момент в рабочее пространство системы. Для просмотра значения любой переменной из текущего рабочего пространства системы достаточно набрать ее имя и нажать клавишу Enter.
После окончания сеанса работы с системой MATLAB все ранее вычисленные переменные теряются. Чтобы сохранить в файле на диске компьютера содержимое рабочего пространства системы MATLAB, нужно выполнить команду меню File / Save Workspace As ... .По умолчанию расширение имени файла mat, поэтому такие файлы принято называть МАТ-файлами. Для загрузки в память компьютера ранее сохраненного на диске рабочего пространства нужно выполнить команду меню: File / Load Workspace ... .
Вещественные числа и тип данных double
Система MATLAB представляет на машинном уровне все действительные числа заданные мантиссой и показателем степени, например, 2.85093Е+11, где буквой Е обозначается основание степени равное 10. Этот основной тип данных носит название double. MATLAB по умолчанию использует формат short для вывода вещественных чисел, при котором показываются только четыре десятичных цифры после запятой.
Введите с клавиатуры пример:
» res=5.345*2.868/3.14-99.455+1.274
Получите результат вычисления:
Если требуется полное представление вещественного числа res, введите с клавиатуры команду:
нажмите клавишу Enter и получите более подробную информацию:
res = -93.29900636942675
Теперь все результаты вычислений будут показываться с такой высокой точностью в течение данного сеанса работы в среде системы MATLAB. Если требуется до прекращения текущего сеанса работы вернуться к старой точности визуального представления вещественных чисел в командном окне, нужно ввести и исполнить (нажав клавишу Enter) команду:
Целые числа показываются системой в командном окне в виде целых чисел.
Над вещественными числами и переменными типа double производятся арифметические операции: сложения +, вычитания -, умножения *, деления / и возведения в степень ^ . Приоритет в выполнении арифметических операций обычный. Операции одинакового приоритета выполняются в порядке слева направо, но круглые скобки могут изменить этот порядок.
Если нет необходимости видеть в командном окне результат вычисления некоторого выражения, то в конце введенного выражения следует поставить точку с запятой и только после этого нажать Enter.
В системе MATLAB присутствуют все основные элементарные функции для вычислений с вещественными числами. Любая функция характеризуется своим именем, списком входных аргументов (перечисляются через запятую и стоят внутри круглых скобок, следующих за именем функции) и вычисляемым (возвращаемым) значением. Список всех имеющихся в системе элементарных математических функций может быть получен по команде help elfun. В Приложении 1 перечислены стандартные функции вещественного аргумента.
Вычислите выражение, включающее вычисление функции арксинус:
Убедитесь, что получился следующий результат:
соответствующее числу «пи». В системе MATLAB для вычисления числа «пи» есть специальное обозначение: pi. (Список системных переменных MATLAB находится в Приложении 2).
MATLAB имеет также логические функции, функции, связанные с целочисленной арифметикой (округления до ближайшего целого: round, усечение дробной части числа: fix). Есть еще функция mod - остаток от деления с учетом знака, sign - знак числа, lcm - наименьшее общее кратное, perms - вычисление числа перестановок и nchoosek - числа сочетаний и много других. Многие из функций имеют область определения, отличную от множества всех действительных чисел.
Помимо арифметических операций над операндами типа double выполняются еще операции отношения и логические операции. Операции отношения сравнивают между собой два операнда по величине. Эти операции записываются следующими знаками или комбинациями знаков (Таблица 1):
Таблица 1
В случае истинности операции отношения ее величина равна 1, а в случае ложности - 0. Операции отношения имеют более низкий приоритет, чем арифметические операции.
Наберите с клавиатуры выражение с операциями отношения и вычислите
» a=1; b=2; c=3;
» res=(a
Вы получите следующий результат:
Логические операции над вещественными числами обозначаются знаками, перечисленными в таблице 2:
Таблица 2
& | | | ~ |
И | ИЛИ | НЕ |
Первые две из этих операций являются бинарными (двухоперандными), а последняя - унарной (однооперандной). Логические операции трактуют свои операнды как «истинные» (не равные нулю) или «ложные» (равные нулю). Если оба операнда операции «И» истинны (не равны нулю), то результат этой операции равен 1 («истина»); во всех остальных случаях операция «И» вырабатывает значение 0 («ложь»). Операция «ИЛИ» вырабатывает 0 («ложь») только в случае, когда являются ложными (равными нулю) оба операнда. Операция «НЕ» инвертирует «ложь» на «истину». Логические операции имеют самый низкий приоритет.
Комплексные числа и комплексные функции
Комплексные переменные, как и вещественные автоматически имеют тип double и не требуют никакого предварительного описания. Для записи мнимой единицы зарезервированы буквы i или j. В случае, когда коэффициентом перед мнимой единицей является не число, а переменная, между ними следует обязательно использовать знак умножения. Итак, комплексные числа можно записывать следующим образом:
» 2+3i; -6.789+0.834e-2*i; 4-2j; x+y*i;
Почти все элементарные функции допускают вычисления с комплексными аргументами. Вычислите выражение:
» res=sin(2+3i)*atan(4i)/(1 -6i)
Получится результат:
1.8009 - 1.91901
Специально для работы с комплексными числами предназначены следующие функции: abs (абсолютное значение комплексного числа), conj (комплексно сопряженное число), imag (мнимая часть комплексного числа), real (действительная часть комплексного числа), angle (аргумент комплексного числа), isreal («истина», если число действительное). Функции комплексного переменного перечислены в Приложении 1.
В отношении арифметических операций ничего нового для комплексных чисел (по сравнению с вещественными) сказать невозможно. То же самое относится и к операциям отношения «равно» и «не равно». Остальные операции отношения вырабатывают результат исходя только из действительных частей этих операндов.
Введите выражение, получите результат и объясните его:
» c=2+3i; d=2i; » c>d
Логические операции трактуют операнды как ложные, если они равны нулю. Если же у комплексного операнда не равна нулю хотя бы одна его часть (вещественная или мнимая), то такой операнд трактуется как истинный.
Числовые массивы
Для создания одномерного массива можно использовать операцию конкатенации, которая обозначается с помощью квадратных скобок . Элементы массива помещаются между скобками и отделяются друг от друга пробелом или запятой:
» al=; d=;
Для доступа к индивидуальному элементу массива нужно применить операцию индексации, для чего после имени элемента указать в круглых скобках индекс элемента.
Можно изменять элементы уже сформированного массива путем применения операций индексации и присваивания. Например, введя:
мы изменим третий элемент массива. Или, после введения:
» al(2)=(al(1)+al(3))/2;
второй элемент массива станет равным среднему арифметическому первого и третьего элементов. Запись несуществующего элемента вполне допустима - она означает добавление нового элемента к уже существующему массиву:
Применяя после выполнения этой операции к массиву а1 функцию length, находим, что количество элементов в массиве возросло до четырех:
Тоже самое действие - «удлинение массива а1» - можно выполнить и с помощью операции конкатенации:
Можно задать массив, прописывая все его элементы по отдельности:
» a3(1)=67; a3(2)=7.8; a3(3)=0.017;
Однако этот способ создания не является эффективным. Еще один способ создания одномерного массива основан на применении специальной функции, обозначаемой двоеточием (операция формирования диапазона числовых значений). Через двоеточие следует набрать первое число диапазона, шаг (приращение) и конечное число диапазона. Например:
» diap=3.7:0.3:8.974;
Если не нужно выводить на экран весь получившийся массив, то в конце набора (после конечного числа диапазона) следует набрать точку с запятой. Чтобы узнать, сколько элементов в массиве, следует вызвать функцию length (имя массива).
Для создания двумерного массива (матрицы) также можно использовать операцию конкатенацию. Элементы массива набираются один за другим согласно их расположению в строках, в качестве разделителя строк используется точка с запятой.
Введите с клавиатуры:
» a=
Нажмите ENTER, получим:
Полученную матрицу а размером 3x2 (первым указывается число строк, вторым - число столбцов) можно сформировать также вертикальной конкатенацией вектор-строк:
» a=[;;];
или горизонтальной конкатенацией вектор-столбцов:
» a=[,];
Структуру созданных массивов можно узнать с помощью команды whos(имя массива), размерность массива - функцией ndims, а размер массива - size.
Двумерные массивы можно задать также с помощью операции индексации, прописывая по отдельности его элементы. Номер строки и столбца, на пересечении которых находится задаваемый элемент массива, указываются через запятую в круглых скобках. Например:
» a(1,1)=1; a(1,2)=2; a(2,1)=3; » a(2,2)=4; a(3,1)=5; a(3,2)=6;
Однако будет намного эффективнее, если до начала прописывания элементов массива, создать массив нужного размера функциями ones (m,n) или zeros(m,n), заполненный единицами или нулями (m - число строк, n - число столбцов). При вызове этих функций предварительно выделяется память под заданный размер массива, после этого постепенное прописывание элементов нужными значениями не требует перестройки структуры памяти, отведенной под массив. Использование этих функций возможно и при задании массивов других размерностей.
Если после формирования массива Х потребуется, не изменяя элементов массива, изменить его размеры, можно воспользоваться функцией reshape (Х, М, N), где M и N - новые размеры массива Х
Объяснить работу этой функции можно, только исходя из способа, каким система MATLAB хранит элементы массивов в памяти компьютера. Она хранит их в непрерывной области памяти упорядоченно по столбцам: сначала располагаются элементы первого столбца, вслед за ними расположены элементы второго столбца и т.д. Помимо собственно данных (элементов массива) в памяти компьютера хранится также управляющая информация: тип массива (например, double), размерность и размер массива, другая служебная информация. Этой информации достаточно для определения границ столбцов. Отсюда следует, что для переформирования матрицы функцией reshape достаточно изменить только служебную информацию и не трогать собственные данные.
Поменять местами строки матрицы с ее столбцам можно операцией транспортирования, которая обозначается знаком." (точка и апостроф). Например,
» A=;
Операция " (апостроф) выполняет транспонирование для вещественных матриц и транспонирование с одновременным комплексным сопряжением для комплексных матриц.
Объекты, с которыми работает MATLAB, являются массивами. Даже одно заданное число во внутреннем представлении MATLAB является массивом, состоящим из одного элемента. MATLAB позволяет делать вычисления с огромными массивами чисел также легко как и с одиночными числами, и это является одним из самых заметных и важных преимуществ системы MATLAB над другими программными пакетами, ориентированными на вычисления и программирование. Помимо памяти, необходимой для хранения числовых элементов (по 8 байт на каждый в случае вещественных чисел и по 16 байт в случае комплексных чисел), MATLAB автоматически при создании массивов выделяет еще и память для управляющей информации.
Вычисления с массивами
В традиционных языках программирования вычисления с массивами осуществляются поэлементно в том смысле, что нужно запрограммировать каждую отдельную операцию над отдельным элементом массива. В М-языке системы MATLAB допускаются мощные групповые операции над всем массивом сразу. Именно групповые операции системы MATLAB позволяют чрезвычайно компактно задавать выражения, при вычислении которых реально выполняется гигантский объем работы.
Операции сложения и вычитания матриц обозначаются стандартными знаками + и -.
Задайте матрицы А и В и выполните операцию сложения матриц:
» A=; B=;
Если используются операнды разных размеров, выдается сообщение об ошибке, за исключением случая, когда один из операндов является скаляром. При выполнении операции А + скаляр (А - матрица) система расширит скаляр до массива размера А, который и складывается далее поэлементно с А.
Для поэлементного перемножения и поэлементного деления массивов одинаковых размеров, а также поэлементного возведения в степень массивов, применяются операции, обозначаемые комбинациями двух символов: .* , ./, и.^. Использование комбинаций символов объясняется тем, что символами * и / обозначены специальные операции линейной алгебры над векторами и матрицами.
Кроме операции./, называемой операцией правого поэлементного деления, есть еще операция левого поэлементного деления.\. Разница между этими операциями: выражение А./ В приводит к матрице с элементами А (k, m) /В (k, m), а выражение А.\В приводит к матрице с элементами В (k, m) /А (k, m).
Знак * закреплен за перемножением матриц и векторов в смысле линейной алгебры.
Знак \ закреплен в системе MATLAB за решением довольно сложной задачи линейной алгебры - нахождением корней системы линейных уравнений. Например, если требуется решить систему линейных уравнений Ay = b, где А - заданная квадратная матрица размера N´N, b - заданный вектор- столбец длины N, то для нахождения неизвестного вектор-столбца у достаточно вычислить выражение А\b (это равносильно операции: A -1 B).
Типичные задачи аналитической геометрии в пространстве, связанные с нахождением длин векторов и углов между ними, с вычислением скалярного и векторного произведений, легко решаются разнообразными средствами системы MATLAB. Например, для нахождения векторного произведения векторов предназначена специальная функция cross, например:
» u=; v=;
Скалярное произведение векторов можно вычислить с помощью функции общего назначения sum, вычисляющей сумму всех элементов векторов (для матриц эта функция вычисляет суммы для всех столбцов). Скалярное произведение, как известно, равно сумме произведений соответствующих координат (элементов) векторов. Таким образом, выражение: » sum(u.*v)
вычисляет скалярное произведение двух векторов u и v. Скалярное произведение можно также вычислить как: u*v".
Длина вектора вычисляется с помощью скалярного произведения и функции извлечения квадратного корня, например:
» sqrt(sum(u.*u))
Ранее рассмотренные для скаляров операции отношения и логические операции выполняются в случае массивов поэлементно. Оба операнда должны быть одинаковых размеров, при этом операция возвращает результат такого же размера. В случае, когда один из операндов скаляр, производится его предварительное расширение, смысл которого уже был пояснен на примере арифметических операций.
Среди функций, генерирующих матрицы с заданными свойствами, часто испольльзуется функция eye , производящую единичные квадратные матрицы, а также широко применяемую на практике функцию rand, генерирующую массив со случайными элементами, равномерно распределенными на интервале от 0 до 1. Например, выражение
порождает массив случайных чисел размером 3х3 с элементами, равномерно распределенными на интервале от 0 до 1.
Если вызвать эту функцию с двумя аргументами, например R=rand(2,3), то получится матрица R случайных элементов размером 2x3. При вызове функции rand с тремя и более скалярными аргументами производятся многомерные массивы случайных чисел.
Определитель квадратной матрицы вычисляется с помощью функции det. Среди функций, производящих простейшие вычисления над массивами, помимо рассмотренной выше функции sum, используется еще и функция prod, которая во всем аналогична функции sum, только вычисляет она не сумму элементов, а их произведение. Функции max и min ищут соответственно максимальный и минимальный элементы массивов. Для векторов они возвращают единственное числовое значение, а для матриц они порождают набор экстремальных элементов, вычисленных для каждого столбца. Функция sort сортирует в возрастающем порядке элементы одномерных массивов, а для матриц она производит такую сортировку для каждого столбца отдельно.
В MATLAB есть уникальная возможность производить групповые вычисления над массивами, используя обычные математические функции, которые в традиционных языках программирования работают только со скалярными аргументами. В результате с помощью крайне компактных записей, удобных для ввода с клавиатуры в интерактивном режиме работы с командным окном системы MATLAB, удается произвести большой объем вычислений. Например, всего два коротких выражения
» x=0:0.01:pi/2; y=sin(x);
вычисляют значения функции sin сразу в 158 точках, формируя два вектора x и у со 158 элементами каждый.
Построение графиков функций
Графические возможности системы MATLAB являются мощными и разнообразными. Изучим наиболее простые в использовании возможности (высокоуровневую графику).
Сформируйте два вектора х и у:
» x=0:0.01:2; y=sin(x);
Вызовите функцию:
и вы получите на экране график функции (рис. 1).
Рис. 1. График функции y=sin(x)
MATLAB показывает графические объекты в специальных графических окнах, имеющих в заголовке слово Figure. Не убирая с экрана дисплея первое графическое окно, введите с клавиатуры выражения
и получите новый график функции в том же самом графическом окне (при этом старые оси координат и график пропадают - этого также можно добиться командой clf, командой cla удаляют только график с приведением осей координат к их стандартным диапазонам от 0 до 1).
Если нужно второй график провести «поверх первого графика», то перед вторичным вызовом графической функции plot нужно выполнить команду hold on, которая предназначена для удержания текущего графического окна:
» x=0:0.01:2; y=sin(x);
Практически тоже самое получится (рис. 2), если набрать:
» x=0:0.01:2; y=sin(x); z=cos(x);
» plot(x,y,x,z)
Рис. 2. Графики функций y=sin(x), z=cos(x), построенные в одном графическом окне
Если нужно одновременно визуализировать несколько графиков так, чтобы они не мешали друг другу, то это можно сделать двумя способами. Первым решением является построение их в разных графических окнах. Для этого перед вторичным вызовом функции plot следует набрать команду figure, которая создает новое графическое окно и заставляет все последующие за ней функции построения графиков выводить их туда.
Вторым решением показа нескольких графиков без конфликта диапазонов осей координат является использование функции subplot. Эта функция позволяет разбить область вывода графической информации на несколько подобластей, в каждую из которых можно вывести графики различных функций.
Например, для ранее выполненных вычислений с функциями sin и cos постройте графики этих двух функций в первой подобласти, а график функции exp(x) - во второй подобласти одного и того же графического окна (рис. 3):
» subplot(1,2,1); plot(x,y,x,z)
» subplot(1,2,2); plot(x,w)
Рис. 3. Графики функций y=sin(x), z=cos(x) и w=exp(x), построенные в двух подобластях одного графического окна
Диапазоны изменения переменных на осях координат этих подобластей независимы друг от друга. Функция subplot принимает три числовых аргумента, первый из которых равен числу рядов подобластей, второй равен числу колонок подобластей, а третий аргумент - номеру подобласти (номер отсчитывается вдоль рядов с переходом на новый ряд по исчерпании). Снять действие функции subplot можно командой:
» subplot(1,1,1)
Если для одиночного графика диапазоны изменения переменных вдоль одной или обеих осей координат слишком велик, то можно воспользоваться функциями построения графиков в логарифмических масштабах. Для этого предназначены функции semilogx, semilogy и loglog.
Построить график функции в полярных координатах (рис. 4) можно с помощью графической функции polar.
» phi=0:0.01:2*pi; r=sin(3*phi);
Рис. 4. График функции r=sin(3*phi) в полярных координатах
Рассмотрим дополнительные возможности, связанные с управлением внешним видом графиков - задание цвета и стиля линий, а также размещение различных надписей в пределах графического окна. Например, команды
» x=0:0.1:3; y=sin(x);
» plot(x,y,"r-",x,y, "ko")
позволяют придать графику вид красной сплошной линии (рис. 5), на которой в дискретных вычисляемых точках проставляют черные окружности. Здесь функция plot дважды строит график одной и той же функции, но в двух разных стилях. Первый из этих стилей отмечен как "r-", что означает проведение линии красным цветом (буква r), а штрих означает проведение сплошной линии. Второй стиль, помечен как "ko", означает проведение черным цветом (буква k) окружностей (буква o) на месте вычисляемых точек.
Рис. 5. Построение графика функции y=sin(x) в двух разных стилях
В общем случае функция plot (x1, y1, s1, x2, y2, s2, ...) позволяет объединить в одном графическом окне несколько графиков функций y1(x1), y2(x2), ... проведя их со стилями s1, s2, ... и т. д.
Стили s1, s2,... задаются в виде набора трех символьных маркеров, заключенных в одиночные кавычки (апострофы). Один из этих маркеров задает тип линии (Таблица 3). Другой маркер задает цвет (Таблица 4). Последний маркер задает тип проставляемых «точек» (Таблица 5). Можно указывать не все три маркера. Тогда используются маркеры, установленные по умолчанию. Порядок, в котором указывают маркеры, не является существенным, то есть "r+-" и "-+r" приводит к одинаковому результату.
Таблица 3. Маркеры, задающие тип линии
Таблица 4 Маркеры, задающие цвет линии
Таблица 5 Маркеры, задающие тип точки
Если в строке стиля поставить маркер на тип точки, но не проставить маркер на тип линии, то тогда отображаются только вычисляемые точки, а непрерывной линией они не соединяются.
Система MATLAB устанавливает пределы на горизонтальной оси равными тем значениям, что указаны пользователем для независимой переменной. Для зависимой переменной по вертикальной оси MATLAB самостоятельно вычисляет диапазон изменения значений функции. Если надо отказаться от этой особенности масштабирования при построении графиков в системе MATLAB, то нужно явным образом навязать свои пределы изменения переменных по осям координат. Это делается с помощью функции axis().
Для проставления различных надписей на полученном рисунке применяют функции xlabel, ylabel, title и text. Функция xlabel создает подпись у горизонтальной оси, функция ylabel - тоже для вертикальной оси (причем эти надписи ориентированы вдоль осей координат). Если требуется разместить надпись в произвольном месте рисунка, применяют функцию text. Общий заголовок для графика создается функцией title. Кроме того, используя команду grid on, можно нанести измерительную сетку на всю область построения графика. Например (рис. 6):
» x=0:0.1:3; y=sin(x);
» plot(x,y,"r-",x,y,"ko")
»title("Function sin(x) graph");
» xlabel("xcoordinate"); ylabel("sin(x)");
» text(2.1, 0.9, "\leftarrowsin(x)"); grid on
Надпись функцией text помещается начиная от точки с координатами, указанными первыми двумя аргументами. По умолчанию координаты задаются в тех же единицах измерения, что и координаты, указанные на горизонтальной и вертикальной осях. Специальные управляющие символы вводятся внутри текста после символа \ (обратная косая черта).
Трехмерная графика
Каждая точка в пространстве характеризуется тремя координатами. Набор точек, принадлежащих некоторой линии в пространстве, нужно задать в виде трех векторов, первый из которых содержит первые координаты этих точек, второй вектор - вторые их координаты, третий вектор - третьи координаты. После чего эти три вектора можно подать на вход функции plot3, которая и осуществит проектирование соответствующей трехмерной линии на плоскость и построит результирующее изображение (рис. 7). Введите с клавиатуры:
» t=0:pi/50:10*pi; x=sin(t);
» y=cos(t); plot3(x,y,t); grid on
Рис. 7. График винтовой линии, построенный с помощью функции plot3
Эту же функцию plot3 можно применить и для изображения поверхностей в пространстве, если, конечно, провести не одну линию, а много. Наберите с клавиатуры:
» u=-2:0.1:2; v=-1:0.1:1;
» =meshgrid(u,v);
» z=exp(-X.^2-Y.^2);
Получите трехмерное изображение графика функции (рис. 8).
Функция plot3 строит график в виде набора линий в пространстве, каждая из которых является сечением трехмерной поверхности плоскостями, параллельными плоскости yOz. Помимо этой простейшей функции система MATLAB располагает еще рядом функций, позволяющих добиваться большей реалистичности в изображении трехмерных графиков.
Рис. 8. График поверхности в пространстве, построенный с помощью функции plot3
Сценарии и m-файлы.
Для простых операций удобен интерактивный режим, но если вычисления нужно многократно выполнять или необходимо реализовывать сложные алгоритмы, то следует использовать m-файлы MATLAB (расширение файла состоит из одной буквы m). script-m-файл (или сценарий) - текстовый файл, содержащий инструкции на языке MATLAB, подлежащими исполнению в автоматическом пакетном режиме. Создать такой файл удобнее с помощью редактора системы MATLAB. Он вызывается из командного окна системы MATLAB командой меню File/New/M-file (или самой левой кнопкой на полосе инструментов, на которой изображен чистый белый лист бумаги). Записанные в script-файлы команды будут выполнены, если в командной строке ввести имя script-файла (без расширения). Переменные, определяемые в командном окне и переменные, определяемые в сценариях, составляют единое рабочее пространство системы MATLAB, причем переменные, определяемые в сценариях, являются глобальными, их значения заместят значения таких же переменных, которые были использованы до вызова данного script- файла.
После создания текста сценария его надо сохранить на диске. Путь к каталогу с файлом обязательно должен быть известен системе MATLAB. Командой File/Set Path вызывается диалоговое окно просмотрщика путей доступа к каталогам. Для добавления нового каталога в список путей доступа необходимо выполнить далее команду меню Path/Add to path.