Графическое построение функций двух переменных, заданных на прямоугольной области.

 

Рассмотрим некоторые возможности МАТЛАБ на примере функции

,

 

Примеры отдельных текстов М-файлов

 

Построение сеточной области

xy=input('[xa ya;xb yb]=');

dl=xy(2,:)-xy(1,:);

n=input('[n1 n2]=');

 

h=dl./n;

x1=xy(1,1):h(1):xy(2,1);

x2=xy(1,2):h(2):xy(2,2);

[x y] = meshgrid(x1,x2);

 

U=-log(x.^2+y.^2+1);

 

Ввод исходных данных в командном окне

[xa ya;xb yb]=[-3 -4;3 4]

[n1 n2]=[10 20]

 

Использование функции surf

 

figure()

surf(x,y,U)

title('SERF: function U')

colormap hsv

colorbar;

axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

 

figure()

surf(x,y,U)

title('SERF: function U')

colormap hsv

colorbar;

%axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

 

figure()

surf(x,y,U)

title('SERF: function U')

colormap copper

colorbar;

axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

shading interp

 

figure()

surf(x,y,U)

title('SERF: function U')

colormap copper

colorbar;

axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

%shading interp

 

Использование функции plot3

 

figure()

plot3(x,y,U,'LineWidth',8)

grid on

title('PLOT3: function U')

colormap hsv

axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

figure()

%plot3(x,y,U,'LineWidth',8)

plot3(x,y,U)

grid on

title('PLOT3: function U')

colormap hsv

axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

 

Использование функции contour3

figure()

[C h]=contour3(x,y,U,10); set(h,'LineWidth',5)

grid on

title('COUNTOUR3: function U')

colormap hsv

colorbar;

axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

figure()

%[C h]=contour3(x,y,U,10); set(h,'LineWidth',5)

contour3(x,y,U,10);

grid on

title('COUNTOUR3: function U')

colormap hsv

colorbar;

axis image;

xlabel('X'),ylabel('Y'),zlabel('z=U(x,y)')

 

Использование функции contourf

 

figure()

contourf(x,y,U,8)

title('CONTOURF: function U')

colormap hsv

colorbar;

axis image;

xlabel('X'),ylabel('Y')

 

figure()

contourf(x,y,U,8)

title('CONTOURF: function U')

%colormap hsv

colorbar;

axis image;

xlabel('X'),ylabel('Y')

 

Использование функции contour

 

figure()

[C h]=contour(x,y,U,8); set(h,'LineWidth',2),clabel(C,h);

title('CONTOUR: function U')

colormap gray

colorbar

axis image;

xlabel('X'),ylabel('Y')

figure()

%[C h]=contour(x,y,U,8); set(h,'LineWidth',2),clabel(C,h);

contour(x,y,U,8);

title('CONTOUR: function U')

%colormap gray

colorbar

%axis image;

xlabel('X'),ylabel('Y')

5.6. Некоторые возможности построения графиков функций трех
переменных с использованием стандартной функции slice

Пример: , ,

Пример М-файла

function PLOT_SLICE

clear all

close all

 

n=input('[n1 n2 n3]=');

dl=input('[dl1 dl2 dl3]=');

 

h=dl./n;

x1=0:h(1):dl(1);

x2=0:h(2):dl(2);

x3=0:h(3):dl(3);

[x y z] = meshgrid(x1,x2,x3);

 

U=sin(pi*x/dl(1))+10*cos(pi*y/dl(2)-pi/2)+(z/dl(3)).^2;

 

figure()

slice(x,y,z,U,[0 dl(1)],[0 dl(2)],[0 dl(3)])

title('SLICE: [0 dl(1)],[0 dl(2)],[0 dl(3)]')

colormap hsv

axis image

colorbar

 

figure()

slice(x,y,z,U,[dl(1)/2 dl(1)],[0 dl(2)],[0 dl(3)])

title('SLICE: [dl(1)/2 dl(1)],[0 dl(2)],[0 dl(3)]')

colormap hsv

axis image

colorbar

 

figure()

slice(x,y,z,U,[0 dl(1)],[dl(2)/2 dl(2)],[0 dl(3)])

title('SLICE: [0 dl(1)],[dl(2)/2 dl(2)],[0 dl(3)]')

colormap hsv

axis image

colorbar

 

figure()

slice(x,y,z,U,[0 dl(1)],[0 dl(2)],[0 dl(3)/2])

title('SLICE: [0 dl(1)],[0 dl(2)],[0 dl(3)/2]')

colormap hsv

axis image

colorbar

 

Ввод исходных данных

[n1 n2 n3]=[10 10 20]

[dl1 dl2 dl3]=[10 10 5]

>>