1. БАЛАКОВСКИЙ ИНСТИТУТ ТЕХНИКИ, ТЕХНОЛОГИИ И УПРАВЛЕНИЯ (ФИЛИАЛ)

2. ГОУ ВПО «САРАТОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

3. ФАКУЛЬТЕТ ИНЖЕНЕРНО – СТРОИТЕЛЬНЫЙ

4. КАФЕДРА «УПРАВЛЕНИЕ И ИНФОРМАТИКА В ТЕХНИЧЕСКИХ СИСТЕМАХ

1. ЛАБОРАТОРНАЯ РАБОТА №1

по дисциплине

Теория автоматического управления

ИЗУЧЕНИЕ ГРАФИЧЕСКОЙ ОПЕРАЦИОННОЙ СРЕДЫ СИСТЕМЫ MATLAB

Выполнил: ст. гр. УИТ-32

Гусев И.А.

Приняла:

Комлева О.А.

2009

ИЗУЧЕНИЕ ГРАФИЧЕСКОЙ ОПЕРАЦИОННОЙ СРЕДЫ СИСТЕМЫ MATLAB

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

>> A=[1 -2 6 1; 5 0 2 -3; 7 -3 7 2; 3 -7 5 2]

A =

1 -2 6 1

5 0 2 -3

7 -3 7 2

3 -7 5 2

>> B=[88; 88; 181; 99]

B =

88

88

181

99

>> x=A\B

x =

13.3491

0.8876

12.5266

1.2663

>> disp(B')

88 88 181 99

>> disp(B*3)

264

264

543

297

>> disp(B+3)

91

91

184

102

>> disp(B.^3)

681472

681472

5929741

970299

>> disp(A')

1 5 7 3

-2 0 -3 -7

6 2 7 5

1 -3 2 2

>> disp(A^3)

307 -366 610 257

247 -128 288 59

569 -617 1019 452

291 -431 517 288

>> disp(inv(A))

-0.1746 0.0325 0.1509 -0.0148

0.0562 -0.0444 0.1124 -0.2071

0.2367 0.0237 -0.0266 -0.0562

-0.1331 -0.2633 0.2337 -0.0621

>> disp(expm(A))

1.0e+004 *

1.5408 -2.0282 3.0426 1.4939

0.8109 -1.0510 1.5917 0.7730

2.6547 -3.4944 5.2423 2.5738

1.3322 -1.7739 2.6426 1.3079

>> disp(logm(A))

2.3118 -0.8668 0.5097 -0.3732

5.1184 1.2127 -2.4870 -0.2365

0.3278 0.1235 1.8530 0.6578

6.3839 -0.9890 -2.7631 1.1387

>> disp(sqrtm(A))

1.6014 -0.7392 1.3011 -0.1332

3.9134 0.9950 -1.1319 -0.8471

1.4467 -0.2570 2.2913 0.5709

4.1588 -1.9028 -0.7413 1.1684

>> X=[1 -2 6 1; 5 0 2 -3; 7 -3 7 2; 3 -7 5 2]

X =

1 -2 6 1

5 0 2 -3

7 -3 7 2

3 -7 5 2

>> disp(cond(X))

6.7016

>> disp(norm(X,1))

20

>> disp(norm(X,2))

15.0175

>> disp(rcond(X))

0.0955

>> disp(rank(X))

4

>> disp((det(X)))

676

>> disp(trace(X))

10

>>

>> disp(null(X))

>> disp(orth(X))

-0.3834 0.2190 -0.7096 -0.5491

-0.2434 -0.7882 0.2718 -0.4956

-0.6851 -0.2663 -0.2107 0.6445

-0.5696 0.5098 0.6150 -0.1937

>> disp(rref(X))

1 0 0 0

0 1 0 0

0 0 1 0

0 0 0 1

>> p=[1 -1 -29 -71 -140]

p =

1 -1 -29 -71 -140

>> r=(roots(p))

r =

7.0000

-4.0000

-1.0000 + 2.0000i

-1.0000 - 2.0000i

>> disp(poly(r))

1.0000 -1.0000 -29.0000 -71.0000 -140.0000

>> disp(polyval(p,3))

-560

>> disp(polyder(p))

4 -3 -58 -71

>> x=0:0.1:10;

>> y=3*x-14+exp(x)-exp(x);

>> plot(x,y);grid

>> x=0:0.1:1;

>> y=3*x-14+exp(x)-exp(x);

>> [TH,R] = cart2pol(x,y);

>> plot(TH,R),grid

>> P1=[-3 -1 2];

>> P2=[-2 -7 -1];

>> roots(P1);

>> roots(P2);

>> om0=1e-2; omk=1e2;

>> OM=logspace(-2,2,41); p=i*OM;

>> ch=polyval(P1,p);zn=polyval(P2,p);

>> ACH=abs(ch)./abs(zn);

>> FCH=angle(ch./zn)*180/pi;

>> subplot(2,1,1);

>> loglog(OM,ACH);grid

>> figure, semilogx(OM,FCH); grid