Система MatLab имеет функцию roots (P), возвращающую вектор, коэффициенты которого являются корнями заданного многочлена Р.

Многочлен задается в виде вектора коэффициентов при переменной, начиная со старшего – например: S(x)=x5+8x4-31x3+80x2+94x-20 в MatLab будет записан:

>> S=[1 8 -31 80 94 -20]

S =

1 8 -31 80 94 -20

>> roots(S) – решение заданного многочлена:

ans =

-11.3026

2.0475 + 2.3585i

2.0475 - 2.3585i

-0.9779

0.1855 ,

где ans—название результирующей переменной по умолчанию.

Операции с матрицами можно выполнять по правилам вычисления матриц, а так же выполнять действия по элементно – сравните:

>> A=[1 2 3;4 5 6;1 1 1];

>> B=[1 2 1; 3 3 3; 3 2 1];

>> C=A*B

C =

16 14 10

37 35 25

7 7 5

>> D=A.*B

D =

1 4 3

12 15 18

3 2 1.

Работа с графическими средствами.

Основные операторы графики:

plot – построение графика в линейном масштабе,

polar – построение графика в полярной системе координат,

mesh – построение графика трехмерной поверхности,

plot3 – построение графика трехмерной поверхности,

surf– построение графика трехмерной поверхности,

text – вывод надписи в заданное место графика,

title – задание титульной надписи,

xlabel – задание надписи по х,

ylabel – задание надписи по у,

grid – задание пунктирной масштабной сетки,

clg – очистка графического экрана,

hold – сохранение предшествующих построений,

shg – показ графического экрана,

pause – задание паузы в вычислениях.

 

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

Тип линий Тип точек Цвета
Черта - Плюс + Красный r
Пунктир -- Звёздочка * Зелёный G
Двоеточие : Кружок о Синий B
Штрих-пунктир -. Крестик х Белый W

 

Примеры.

1. Построить графики затухающих колебаний: x(t)=e-0.2tsin(t), y(t)=e-0.2tcos(t), где t изменяется от 0 до 10 с шагом 0,1. Выполнение этого задания должно выглядеть следующим образом:

t=0:.1:10;

x=exp(-.2*t).*sin(t);

y= exp(-.2*t).*cos(t);