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

03 Графики в MATLAB_2

.pdf
Скачиваний:
73
Добавлен:
12.03.2015
Размер:
4.86 Mб
Скачать

Example 5:

x = cos(t), y = cos(t)2 , z = sin (t), u = sin(t)2 , v = esin(t),

 

w = ecos(t)

>>t = -6.28:0.05:6.28; x = cos(t);

y = cos(t).^2; z = sin(t);

u = sin(t).^2;

v = exp(sin(t)); w = exp(cos(t));

% 6

%

subplot(3,2,1);

plot(t,x); grid on;

title('Graph of function {\it x =cos(t)}'); xlabel('t');

ylabel('x');

>>%

>>subplot(3,2,2)

>>plot(t,y);

>>grid on;

>>title('{\bf y=cos(t)^2}');

>>xlabel('t');

>>ylabel('y');

>>subplot(3,2,3);

>>plot(t,z);

>>grid on;

>>title('{\bf z=sin(t)}');

>>xlabel('t');

>>ylabel('z');

>>subplot(3,2,4);

>>plot(t,u);

>>grid on;

>>title('{\bf u=sin(t)^2}');

>>xlabel('t');

>>ylabel('u');

>>subplot(3,2,5);

>>plot(t,v);

>>grid on;

>>title('{\bf v=exp(sin(t))}');

>>xlabel('t');

>>ylabel('v');

>>subplot(3,2,6);

>>plot(t,w);

>>grid on;

>>title('{\bf w=exp(cos(t))}');

>>xlabel('t');

>>ylabel('w');

25

26

Scatter(x,y)

x

y

 

 

 

 

0.02

12

 

 

 

 

1.2

13

Scatter(x,y)

 

 

2.5

15

 

 

 

 

3.6

19

 

 

 

 

3.8

22

 

 

 

 

4.2

26

 

 

 

 

5.1

33

 

 

 

 

6.2

42

 

 

 

 

6.5

49

 

 

 

 

6.9

55

 

 

 

 

8.9

66

 

 

 

 

11.2

74

 

 

 

 

12.1

82

 

 

 

 

13.0

95

 

 

 

 

13.2

102

 

 

 

 

13.3

109

 

 

 

 

13.6

112

 

 

 

 

15.6

118

 

 

 

 

16.0

119

 

 

 

 

17.2

123

 

 

 

 

17.8

131

 

 

 

 

18.6

133

27

 

 

polyfit ( x, y, n); находить подходяший полином n-того порядка и его Коэффициенты в количестве (n+1), который наилчушим образом проходит через данные точки.

polyval( y, x); находить соответствущие на х, значение ординат, используя предварительно определенный полином.

>>x = [0.02 1.2 2.5 3.6 3.8 4.2 5.1 6.2 6.5 6.9 8.9 11.2 12.1 13.0 13.2 13.3 13.6 15.6 16.0 17.2 17.8 18.6];

>>y = [12 13 15 19 22 26 33 42 49 55 66 74 82 95 102 109 112 118 119 123 131 133];

>>scatter(x,y); %

>>p = polyfit(x,y,3); % 3

%

>>z = polyval(p,x); %

,

>>hold on

>>plot(x,z); %

z

28

MATLAB

!" # $% & &, z = f(x,y), ';

2.z !

3." #

29

" (x,y)

. "

.

 

 

 

 

( )$ & $

 

 

 

 

10

 

 

% .

 

 

 

 

 

5

 

 

x =-3 -2 -1

0

1

2

3

 

 

 

axis-

 

 

 

y =-3 -2 -1

0

1

2

3

 

 

 

0

 

 

z=f(x,y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

-5

 

 

) & )$ ( $ .

 

 

 

 

 

 

 

 

 

 

 

-10

 

 

xx =-3

-2

-1

 

0

1

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-3

-2

-1

 

0

1

2

3

 

 

2

 

 

-3

-2

-1

 

0

1

2

3

 

 

 

 

2

 

 

 

0

 

-3

-2

-1

 

0

1

2

3

 

 

0

 

 

 

 

 

 

-3

-2

-1

 

0

1

2

3

 

 

-2

-2

 

 

 

 

 

-3

-2

-1

 

0

1

2

3

 

 

y-axis

 

 

 

 

x-axis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-3

-2

-1

 

0

1

2

3

 

 

 

 

 

 

 

 

 

 

 

yy =-3

-3

-3

 

-3

-3

-3

-3

 

 

 

 

 

-2

-2

-2

 

-2

-2

-2

-2

 

 

 

 

 

-1

-1

-1

 

-1

-1

-1

-1

 

 

, xx ,

 

0

0

0

 

0

0

0

0

 

 

 

1

1

1

 

1

1

1

1

 

 

#

 

 

 

 

 

 

 

2

2

2

 

2

2

2

2

 

 

 

 

 

 

 

 

 

 

 

3

3

3

 

3

3

3

3

30

 

 

 

 

z=f(xx,yy)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

$# meshgrid()

% , #

# meshgrid

>>x=-3.5:3.5;

>>y=-2.5:2.5;

>>[xx,yy]=meshgrid(x,y)

xx =

 

 

 

 

 

 

 

-3.5

-2.5

-1.5

-0.5

0.5

1.5

2.5

3.5

-3.5

-2.5

-1.5

-0.5

0.5

1.5

2.5

3.5

-3.5

-2.5

-1.5

-0.5

0.5

1.5

2.5

3.5

-3.5

-2.5

-1.5

-0.5

0.5

1.5

2.5

3.5

-3.5

-2.5

-1.5

-0.5

0.5

1.5

2.5

3.5

-3.5

-2.5

-1.5

-0.5

0.5

1.5

2.5

3.5

yy =

 

 

 

 

 

 

 

-2.5

-2.5

-2.5

-2.5

-2.5

-2.5

-2.5

-2.5

-1.5

-1.5

-1.5

-1.5

-1.5

-1.5

-1.5

-1.5

-0.5

-0.5

-0.5

-0.5

-0.5

-0.5 -0.5 -0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

0.5

1.5

1.5

1.5

1.5

1.5

1.5

1.5

1.5

2.5

2.5

2.5

2.5

2.5

2.5

2.5

2.5

(i, j), !

# ,

yy ! :

xx(3,1) = -3.5 yy(3,1) = -0.5

& , (x,y) = (-3.5, -0.5)

31

' z=f(x,y)

>>x=-3:3;

>>y=-3:3;

>> [xx,yy]=meshgrid(x,y)

 

 

 

20

 

xx =

 

 

 

 

 

 

 

 

3

-2

-1

0

1

2

3

15

 

-3

-2

-1

0

1

2

3

 

 

 

-3

-2

-1

0

1

2

3

 

 

-3

-2

-1

0

1

2

3

10

 

-3

-2

-1

0

1

2

3

 

 

-3

-2

-1

0

1

2

3

5

 

-3

-2

-1

0

1

2

3

 

 

 

yy =

 

 

 

 

 

 

 

 

-3

-3

-3

-3

-3

-3

-3

0

 

-2

-2

-2

-2

-2

-2

-2

4

 

 

 

-1

-1

-1

-1

-1

-1

-1

2

4

 

 

0

0

0

0

0

0

0

0

2

0

1

1

1

1

1

1

1

-2

-2

2

2

2

2

2

2

2

 

-4

-4

3

3

3

3

3

3

3

 

 

>> zz=xx.^2 + yy.^2;

 

 

 

, * ) *

 

>> mesh(xx,yy,zz)

 

 

 

 

 

% * (+ % '$ * ,

 

 

 

 

 

 

 

& % z,

( ) ,(

32

 

 

 

>>x=-3:0.1:3;

>>y=-3:0.1:3;

>>[xx,yy]=meshgrid(x,y);

>>zz=xx.^2 + yy.^2;

>>mesh(xx,yy,zz)

33

>>x=-3:0.1:3;

>>y=-3:0.1:3;

>>[xx,yy]=meshgrid(x,y);

>>zz=xx.^2 + yy.^2;

>>mesh(xx,yy,zz)

axis tight box on view(120,30)

34