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

Графические возможности системы MATLAB являются мощными и разно­образными. Изучим наиболее простые в использовании возможности (высоко­уровневую графику).

Сформируйте два вектора х и y:

» x=0:0.01:2; y=sin(x);

Вызовите функцию:

» plot(x,y)

и вы получите на экране график функции (рис. 1).

Рис. 1. График функции y=sin(x)

MATLAB показывает графические объекты в специальных графических окнах, имеющих в заголовке слово Figure.

Не убирая с экрана дисплея первое графическое окно, введите с клавиату­ры выражения

» z=cos(x);

» plot(x,z)

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

Если нужно второй график провести «поверх первого графика», то перед вторичным вызовом графической функции plot нужно выполнить команду hold on, которая предназначена для удержания текущего графического окна:

» x=0:0.01:2; y=sin(x);

» plot(x,y)

» z=cos(x);

» hold on

» plot(x,z)

Практически тоже самое получится (рис. 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(х) - во второй подобласти одного и того же графического окна (рис. 3):

» w=exp(x);

» 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);

» polar(phi,r)

Рис. 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

 

Маркер	Цвет линии	Маркер	Цвет линии
c	Голубой	g	Зеленый
m	Фиолетовый	b	Синий
y	Желтый	w	Белый
r	Красный	k	Черный

 

Маркеры, задающие тип точки

Таблица 5

 

Маркер		+	*	о	X
Тип точки	Точка	Плюс	Звездочка	Кружок	Крестик

Если в строке стиля поставить маркер на тип точки, но не проставить мар­кер на тип линии, то тогда отображаются только вычисляемые точки, а непре­рывной линией они не соединяются.

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

axis ([ xmin, xmax, ymin, ymax ])

Для проставления различных надписей на полученном рисунке применяют функции xlabеl, ylabеl, title и text. Функция xlabеl создает подпись у горизон­тальной оси, функция ylabеl - тоже для вертикальной оси (причем эти надписи ориентированы вдоль осей координат). Если требуется разместить надпись в произвольном месте рисунка, применяем функцию text. Общий заголовок для графика создается функцией title. Кроме того, используя команду grid on, мож­но нанести измерительную сетку на всю область построения графика. Напри­мер (рис. 6):

Рис. 6. График функции y=sin(x), построенный в двух стилях, с надпи­сями на координатных осях и на рисунке

» 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 помещается начиная от точки с координатами, указанными первыми двумя аргументами. По умолчанию координаты задаются в тех же единицах измерения, что и координаты, указанные на горизонтальной и вертикальной осях. Специальные управляющие символы вводятся внутри текста после символа \ (обратная косая черта).

Трехмерная графика

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

Каждая точка в пространстве характеризуется тремя координатами. Набор точек, принадлежащих некоторой линии в пространстве, нужно задать в виде трех векторов, первый из которых содержит первые координаты этих точек, второй вектор - вторые их координаты, ну а третий вектор - третьи координа­ты. После чего эти три вектора можно подать на вход функции 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;

» [X,Y]=meshgrid(u,v);

» z=exp(-X.^Ù2-Y.^Ù2);

» plot3(X,Y,z)

Получите трехмерное изображение графика функции (рис. 8).

Рис. 8. График поверхности в пространстве, построенный с помо­щью функции plot3

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

Сценарии и 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.

⇐ Предыдущая123

Поделиться:

Воспользуйтесь поиском по сайту: