Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ЛР_ОТП_M03434

.pdf
Скачиваний:
4
Добавлен:
07.02.2016
Размер:
744.43 Кб
Скачать

 

11

trace

- сума діагональних елементів матриці

inv

- знаходження оберненої матриці

size

- розмір матриці

reshape

- зміна розміру матриці

length

- довжина вектора

diag

- діагональні елементи матриці

triu

- верхня трикутна частина матриці

tril

- нижня трикутна частина матриці

Виконаємо наступні команди:

 

A=reshape(1:24,4,6)

triu(A)

triu(A,0) triu(A,2) triu(А, -1)

tril(A)

 

 

 

v=1:5 diag(v) diag(v,2) diag(v, -1)

diag(A) diag(A,2) diag(А, -1)

 

A=reshape(1:24,4,6)

rot90(A)

rot90(A,2)

Деякі базові матричні операції

А'

- транспонування матриці

А+В

- додавання (віднімання) матриць

А*В

- добуток матриць

А/В

- праве ділення матриць

А\В

- ліве ділення матриць

А^р

- піднесення матриці до степеня

а^Р

- піднесення скаляра до матричного степеня

А.*В

- поелементний додаток матриць

А./В

- поелементне ділення матриць

А.^р

- поелеменнтне піднесення матриці до степеня

Виконаємо наступні команди:

m=3;n=4; a=rand(m,n) b=rand(m,n) c=rand(1,n)p1=a*b p2=a.*b p3=a*c p4=a.*c

Текстові змінні

Вони беруться в лапки (на букві “є” на латинському регістрі), символ займає 2 байти. Виконаємо наступні команди:

>>t='Аргумент'; f=’Функція’; >>x=[pi:pi/180:pi];y=sin(x); >>disp([t,’ ‘,f]),disp([x’, y’])

PDF created with pdfFactory Pro trial version www.pdffactory.com

12

Формати виводу чисел в пакеті MatLab

format short - з фіксованою точкою та 4 знаками після точки; format long - з фіксованою точкою та 14 знаками після точки; format short e - наукова нотація з 4 десятинними знаками; format long e - наукова нотація з 15 десятинними знаками;

format short g - зручний для виводу цілих. Дійсні числа виводяться з фіксованою точкою та 5 знаками після точки; format long g - зручний для цілих чисел. Дійсні числа виводяться

в науковій нотації з 13 десятинними знаками. Опція format short e дозволяє одержувати рівні стовпці.

1.5 Індивідуальні завдання до лабораторної роботи

Для заданої функції F(a, b, c) і значень параметрів а и b в таблиці 1.1 виконати в системі MatLab:

ввести значення с - вектору аргументів;

зберегти значення a, b, c у mat-файлі;

очистити вікна Сommand Window, Command History, Workspace, завантажити данні з mat-файлу, обчислити функцію F(a,b,c), згідно завдання варіанта;

результуючі значення вектору аргументів с та значень функції F(a, b, c) вивести в вигляді вертикально розташованої таблиці;

повторити вивід цієї таблиці задаючи різні формати даних;

провести запис на диск фрагменту роботи. Створити текстові коментарі до роботи в полученному файлі.

Таблиця 1.1 – Варіанти завдань

 

 

 

 

 

 

 

 

F(a, b, c)

 

a

 

b

 

 

 

 

 

 

 

1

 

 

 

2456

 

0.00078

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

0.02456

 

0.007823

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PDF created with pdfFactory Pro trial version www.pdffactory.com

3

4

5

6

8

9

10

11

12

13

14

15

16

13

0.2456 0.20078

0.024560.007823

0.12456 0.0078 0.024560.007823

0.2456 0.20078

0.024560.007823

0.12456

 

0.078

 

 

 

0.024560.007823

0.2456 0.20078

0.024560.007823

2456 0.00078

0.024560.007823

0.2456 0.20078 0.024560.007823

PDF created with pdfFactory Pro trial version www.pdffactory.com

14

 

 

 

 

 

 

 

17

 

 

 

0.12456

 

0.0078

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18

 

 

 

0.02456

 

0.007823

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19

 

 

 

0.2556

 

0.50078

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

0.02456

 

0.007823

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.6 Завдання на самоперевірку

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

1.>>a=[1 2 3; 4 5 6;7 8 9]

2.>> b1=[1; 2; 3]

3.>>c=[a b1]

4.>> d=zeros(3,4)

5.>> e=[1 2 3; 4 5 ]

6.>>b=[a(1,2:end) a(3,end-1:-1:1)]

7.>> b2=[a(1,2:end); a(3,end-1:-1:1)]

8.>> a[:,2]=[]

Контрольні питання

1.Перелічити тематичні підкаталоги MatLab (toolbox) та їх призначення.

2.Робочій стіл пакету MatLab (меню, інструменти).

3.Типи та формати даних; системні змінні.

4.Елементарні математичні функції пакету.

5.Основні команди керування командного режиму.

6.Призначення та правила використовування команди help.

PDF created with pdfFactory Pro trial version www.pdffactory.com

15

ЛАБОРАТОРНА РОБОТА № 2 ОФОРМЛЕННЯ РЕЗУЛЬТАТІВ ОБЧИСЛЕНЬ

2.1 Мета роботи

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

2.2 Завдання до лабораторної роботи

Студент повинен:

1.Ознайомитись з основними теоретичними відомостями - розділ 2.3 методичних вказівок та розділи 1, 2, 4 електронного;

2.Ознайомитись з основними теоретичними відомостями методичних вказівок - розділ 2.3 та розділ 9 електронного довідника по системі MatLab.. Виконати наведенні приклади;

3.Виконати в MatLab приклади побудови двомірних графіків - розділ 2.4;

4.Виконати в MatLab індивідуальне завдання - розділ 2.5;

5.Скласти звіт - протокол виконання індивідуального завдання та відповісти на контрольні питання.

2.3 Основні теоретичні відомості 2.3.1 Побудова двовимірних графіків

Основні команди побудови двовимірних графіків:

р1оt(x,y)

- побудова графіка в декартовій системі

координат;

 

loglog(x,y)

- побудова графіка в логарифмічних осях;

semilogx(x,y)

- побудова графіка з логарифмічною віссю

абсцис;

 

semilogy(x,y)

- побудова графіка з логарифмічною віссю

ординат;

 

polar(φ,r)

- побудова графіка в полярній системі ко-

ординат;

 

plotyy(x,y,x,f)

- забезпечує вивід осей як ліворуч, так і

 

праворуч, для кожної функції y та f окре-

 

мо;

fplot(<nameFun>,limits)

- побудова графіка функції в інтервалі.

PDF created with pdfFactory Pro trial version www.pdffactory.com

16

Функція plot

р1оt(x) - побудова графіка залежності значення елемента вектора від його порядкового номера.

р1оt (х,у) - будує графік залежності у=f (х).

р1оt (х1, у1, х2, у2, х3, у3) - на одному графіку виводиться декілька залежностей.

р1оt (x, y, ‘color_style_marker’) – вказівкою рядка color_style_marker (до 3 символів), ідентифікують колір, тип лінії та маркер.

Таблиця 2.1 – Символи, що ідентифікують колір, тип лінії та маркер

 

Колір

 

Тип лінії

Тип марке-

 

 

ра

 

 

 

 

 

 

 

 

 

c

cyan(голубий)

неперервна

+

 

 

 

 

 

k

black(чорний)

− −

штрихова

о

 

 

 

 

 

y

yellow(жовтий)

:

пунктирна

*

 

 

 

 

 

r

red(червоний)

.

штрих пунктир-

х

 

 

 

на

 

g

green(зелений)

none

відсутність лінії

.

 

 

 

 

 

b

blue(синій)

 

 

s -квадрат

 

 

 

 

 

w

white(білий)

 

 

<, >, ^ -

 

 

трикутники

 

 

 

 

 

 

 

 

 

Наприклад:

рlot (х, у, 'у:+' ) - виведе залежність у(х) жовтою пунктирною лінією, на якій кожна пара даних буде відмічена маркером +.

У випадку, якщо визначено тип маркера, але не визначено стиль лінії, на екран будуть виведені тільки точки залежності, задані у вхідних векторах.

Функція fplot

fplot(<ім’я функції> ,limits) - будує графік функції (функцій) в

інтервалі limits=[xmin,xmax].

<ім’я функції> - це ім’я М-файлу в вигляді:

@ MyFun або 'MyFun''

PDF created with pdfFactory Pro trial version www.pdffactory.com

17

або рядок типу:

'sin(x)',

'[sin(x) cos(x) ]',

'[sin(x), myfun1(x), myfun2(x)]'.

limits=[xmin xmax ymin ymax] - встановлює інтервали для для значень аргументу та функції.

Наприклад: fplot(@humps,[0 1]); fplot('humps',[0 1]);

f = inline('abs(exp(-j*x*(0:9))*ones(10,1))'); fplot(f,[0 2*pi]); fplot('[tan(x),sin(x),cos(x)]',2*pi*[-1 1 -1 1]); fplot('sin(1 ./ x)', [0.01 0.1],1e-3).

Декілька допоміжних команд

figure – команда для відкриття нового вікна для виводу графіків. figure(n) – для вибору поточного вікна серед відкритих графіч-

них вікон, де n - номер вікна.

hold on - для виводу нових залежностей у поточному вікні графіка. При цьому MatLab не знищить створений вже графік, а тільки, якщо це необхідно, змінить параметри осей.

subplot (vha) - пакет дозволяє здійснювати поділ графічного вікна на декілька підграфіків зі своїми осями. Такий поділ з одночасним вибором активного вікна здійснюється за допомогою наведеної команди, в якій v та h – число під графіків по вертикалі та горизонталі, а параметр а визначає номер активного вікна.

Наприклад, команда subplot (221), plot (x,y)

створить 4 вікна та виведе графік у перше вікно верхньої половини екрана.

2.3.2 Оформлення графіків

Для оформлення графіків служать команди:

axis ([xmin xmax ymin ymax]) - дає змогу користувачеві само-

му задавати діапазон координат;

axis equal - формує однаковий крок по х- та у-осям;

axis оff - дозволяє не виводити градуіровку координатних осей;

PDF created with pdfFactory Pro trial version www.pdffactory.com

18

grid on/off- дозволяє виводити/не виводити сітку на графік;

хlabel (‘назва’), уlabel (‘назва’) - назви x-,у-осей, відповідно; title (‘назва’) – назва графіка;

text (х, у, ‘текст’) - зробити підпис в будь-якому місці графіка, тут x, y координати початку напису.

2.4 Приклади побудови двовімірних графіків

Приклад 1. Побудови характеристики неробочого ходу генера-

тора.

Маємо результати експерименту:

Ток(Im):

x

Напруга(U):

y

1.5

 

80

 

2.13

 

100

 

2.6

 

120

 

3.09

 

140

 

3.6

 

160

 

4.15

 

180

 

4.52

 

190

 

4.78

 

200

 

5.2

 

212

 

5.59

 

220

 

.6.77

 

240

 

7.74

 

252

 

8.55

 

260

 

9.13

 

264

 

10.58

 

276

 

Функція, що апроксимує ці експериментальні дані має вигляд: y1=239.984*atan (0.222*x) ;

Команди MatLab для побудови графіка характеристики неробочого ходу генератора (експериментальних даних та апроксимуючої функції):

>>y1=239.984*atan(0.222*x)

>>plot (x, y, ‘ko’, x, y1, ‘k: ‘) >>grid

>>axis ([0 11 0 286]) >>ylabel (‘U, [B] ‘) >>xlabel (‘Im, [A] ‘)

>>text (3,60,'o - експериментальні дані')

PDF created with pdfFactory Pro trial version www.pdffactory.com

19

>>text (3,40,'--- - U=239. 984arctg (0.2221Im) ')

>>title (‘Апроксимація характеристики неробочого ходу генератора’)

Результат роботи поданої програми показано на рис.2.1.

Рисунок 2.1- Вигляд графічного вікна після виконання команд програми

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

Приклад 2. Графік функцій, заданих параметрично x(t)=0.5sin(t), y(t)=0.7cos(t) для t € [0, 2π] (еліпс)

будується при допомозі слідуючих команд:

>>t=[0:0.01:2*pi];

>>x=0.5*sin(t);

>>y=0.7*cos(t);

>>plot(x,y)

PDF created with pdfFactory Pro trial version www.pdffactory.com

20

Приклад 3. Графік кусочно-заданої функції

ì

 

 

 

- £ x £ -π ,

ïπ ×sin( x),

y(x) = íπ -

x

,

 

-π < x < π ,

ïπ ×sin 3

x,

π £ x £ 2π.

î

 

 

 

 

будується при допомозі слідуючих команд:

>>x1=[-2*pi:0.01:-pi]; >>y1=pi*sin(x1); >>x2=[-pi:0.01:pi]; >>y2=pi-abs(x2); >>x3=[pi:0.01:2*pi]; >>y3=pi*sin(x3).^3; >>x=[x1 x2 x3]; >>y=[y1 y2 y3]; >>plot(x,y)

Можливо побудувати й по іншому – як три різні функції:

>>plot(x1,y1,x2,y2,x3,y3)

2.5 Індивідуальні завдання до лабораторної роботи

За допомогою двох функцій Matab, plot() та fplot(), побудувати двовимірні графіки згідно індивідуального завдання представленого в табл.2.2 :

1. Графіки функцій y(x) та g(x)(наведені в 1 стовпчику таблиці) вивести:

в одному графічному вікні;

в різних графічних вікнах;

в одному графічному вікні, в різних підграфіках зі своїми

осями.

Побудову кусочно-заданої функції g(x) провести різними вариантами.

2.В 2 стовпчику таблиці наведені параметрично задані функції y(x) та z(x), де параметр x приймає значення від 0 до 2 π. Побудувати графік залежності функції z(x) від y(x ) в новому графічному вікні.

3.Оформити графіки використовуючи:

команди MatLab;

можливості графічного редактора.

PDF created with pdfFactory Pro trial version www.pdffactory.com