2-graficas con matlab

Centro de Estudios EnergéticosDepartamento de Ingeniería Mecánica

Dirección de Investigación y PostgradoUNEXPO Vicerrectorado de Puerto Ordaz

Introducción al MatLab y Simulink


Elaborado por: Prof. Edgar Gutiérrez, Ing. MSc.http://edgar-gutierrez.blogspot.com/



Contenidos:• Introducción al Curso.• Fundamentos del MatLab.• Gráficas con MatLab.• Algebra Lineal.• Integración y Diferenciación Numérica.• Raíces de Ecuaciones no Lineales.• Resolución de Ecuaciones Diferenciales.• Simulink de MatLab.




Gráficas con MatLab.


Resumen de Comandos para Graficar.





Funciones Gráficas 2D.





Funciones Gráficas 3D.





Grafica lineal simple, y la estructura es la siguiente: plot(x,y) grafica y vs x.plot(y) grafica la columna de y vs su subíndice.plot(x1,y1,S1,x2,y2,S2,x3,y3,S3,...) Grafica en una sola gráfica las curvas x1 vs y1, x2 vs y2, x3 vs y3 … y las identifica con S1, S2, S3, … El identificador es una cadena de tres caracteres como máximo. Los caracteres del identificador son:

Función Plot.





Ejemplo:Función Plot.

x=linspace(0,2*pi,20);y=sin(x);z=cos(x);plot(x,y,’b:’,x,z,’rv’);





Ejemplo:Característica de la Ventana Gráfica

x = [0:.2:20];y = sin(x)./sqrt(x+1);y(2,:) = sin(x/2)./sqrt(x+1);y(3,:) = sin(x/3)./sqrt(x+1);plot(x,y)

Barra de herramientas

Ventana de gráfica de Matlab

Anclar/desanclar del desktop de Matlab

Eje coordenado Gráficas

Para acceder al editor de las gráficas se puede hacer desde la gráfica mediante el menú, o usando el comando:

plottools





Modificación de Características por comandosSe pueden modificar:

– Tipo de letra.– Tamaño de letra.– Anchura de las líneas.– Letras griegas.

Ejemplo:t=0:.05:2*pi;w=1;x=cos(w*t);plot(t,x,’r’,’LineWidth’,4);title(’Grafica x_1 con \omega=1’);text(2,0.2,’COS’,’fontname’,’arial’,’fontsize’,14);





xlabel, ylabel, zlabel: Anota una cadena de carácter en los ejes respectivos. La estructura es la siguiente:

Anotaciones sobre las Gráfica

xlabel(‘Eje x’); ylabel(‘Eje y’);zlabel(‘Eje z’);

title: Anota un a cadena de carácter como titulo de un gráfico. La estructura es la siguiente:

title(‘Titulo’);

h=xlabel(‘Eje x’); La variable h almacena un valor que si es -1 indica que no se ejecutó el comando

h=title(‘Titulo’); La variable h almacena un valor que si es -1 indica que no se ejecutó el comando

text: Anota un texto dentro de la gráfica. La estructura es:text(x,y,‘Texto’); %x,y coordenada en el planotext(x,y,‘Texto’); %x,y,z coordenada en el espacio





texlabel: Anota una cadena de carácter en formato TeX. La estructura es:


texlabel(f) texlabel(f,'literal')

Ejemplo:text(.5,.5,texlabel('lambda12^(3/2)/pi - pi*delta^(2/3)'));text(.25,.25,texlabel('lambda12^(3/2)/pi - pi*delta^(2/3)','literal'));

gtext: Anota una cadena de carácter con el mouse. La estructura es:

gtext(‘Texto')





Anotaciones sobre las Gráficalegend: Muestra la leyenda de las curvas graficadas. La estructura es:legend('string1','string2',...)legend('string1','string2',...,'Location',location)h = legend('string1','string2',...)La h es una cadena que se usa para etiquetar el gráfico.

Especificador Localización en el cuadroNorth Interior superiorSouth Interior inferiorEast Interior derechoWest Interior izquierdoNorthEast Interior superior derechoNorthWest Interior superior izquierdoSouthEast Interior inferior derechoSouthWest Interior inferior izquierdoNorthOutside Exterior superiorSouthOutside Exterior inferiorEastOutside Exterior derechoWestOutside Exterior izquierdoNorthEastOutside Exterior superior derechoNorthWestOutside Exterior superior izquierdoSouthEastOutside Exterior inferior derechoSouthWestOutside Exterior inferior izquierdoBest La mejor ubicación interiorBestOutside La mejor ubicación exterior

grid: Activa do desactiva la rejilla de la figura.grid on: Activa la rejilla mayor de los ejes.grid off: Desactiva la rejilla mayor (principal) y menor (secundaria) de los ejes.grid minor: Cambia a la rejilla menor.grid: Cambia a la rejilla mayor.






Ejemplo:x = -pi:pi/20:pi;plot(x,cos(x),'-ro',x,sin(x),'-.b');h = legend('cos_x','sin_x',2);set(h,'Interpreter','none')





Otros Parámetros de Modificación de GráficaAxis: Controla la escala de los ejes coordenados y apariencia. Las estructura es la siguiente:

axis([XMIN XMAX YMIN YMAX]): escala los ejes x, y del plot actual.axis([XMIN XMAX YMIN YMAX ZMIN ZMAX]) : escala los ejes x, y, z de plot 3D actual.axis auto: Retorna los ejes a la escala por defecto.axis equal: Iguala los ejes coordenados.axis off: Desactiva las etiquetas de los ejes, marcas de graduación y fondo.axis on: Reactiva las etiquetas de los ejes, marcas de graduación y fondo.

Hold: Mantiene el gráfica actual:hold on: Mantiene la características actual y las propiedades de los ejes para la siguiente operaciones de gráficas.hold off: Regresa a las propiedades por defectos.hold all: Mantiene las características, color y estilos para las siguientes operaciones.





Otros Parámetros de Modificación de GráficaFigures: Crea una ventana gráfica. Las estructura es la siguiente:

figure(H): H es el un número entero que representa la ventana gráfica activa.

Gcf: Devuelve el identificador de la gráfica actual activa.H=gcf: H almacena el valor de la gráfica activa.

Get: Obtiene las propiedades del objeto.v=get(H,'PropertyName'): Devuelve en v las propiedades del gráfico H.v=get(H,'Default'): Devuelve en v las propiedades por defecto del gráfico H.

Set: Establece las propiedades del objeto.set(H,'PropertyName',Value): Establece el valor de la propiedad en la gráfica H. Para varias gráfica H es un vector.set(H,'PropertyName1',Value1,'PropertyName2',Value2,...): Establece el valor de varias propiedad en la gráfica H.





Comando SubplotCrea varias gráficas en una misma figura. Las estructura es la siguiente:

H=supplot(m,n,p): Crea varias figuras en una ventana gráfica posicionada en forma matricial mxn y p representa la ubicación de la figuraEjemplo:income = [3.2 4.1 5.0 5.6];outgo = [2.5 4.0 3.35 4.9];subplot(2,1,1); plot(income);title('Income');subplot(2,1,2); plot(outgo);title('Outgo')





Comando Subplot





Comando SubplotEjemplo:

Figurefor i=1:12

subplot(12,1,i);plot (sin(1:100)*10^(i-1));set(gca,'xtick',[],'ytick',[]);

End% Resetea los eje y se crea% los datos del Subplot inferior.set(gca,'xtickMode', 'auto')





Comando AxesEl comando axes crea ejes arbitrarios en una figura. La estructura es la siguiente:axes: Función simple para creación de ejes.axes('PropertyName',value,...):Crea los ejes con las propiedades especificadas.h = axes(...): Devuelve en h el identificador de los ejes creados.Ejemplo:axes('position',[.1 .1 .8 .6]);mesh(peaks(20));axes('position',[.1 .7 .8 .2]);pcolor([1:10;1:10]);





Permite graficar una función definiendo los límites del intervalo: fplot(FUN,LIMS)

Donde FUN es la función a graficar y LIMS = [XMIN XMAX] es un vector que define el dominio a graficar. También se puede limitar el rango definiendo LIMS = [XMIN XMAX YMIN YMAX].

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

Función fplot.





Permite graficar con dos escalas en y, una en el eje izquierdo y otra en el derecho. La sintaxis es la siguiente:

plotyy(X1,Y1,X2,Y2): Y1 con datos X1 en el eje izquierdo y Y2 con datos X2 en eje derecho.plotyy(X1,Y1,X2,Y2,FUN): Utiliza la función FUN para especificar el tipo de gráfico. FUN puede ser una función o un cadena que especifica el tipo de gráfica (semilogx, semilogy, loglog, stem, etc,) o cualquier función que acepta la sintaxis H = FUN (X, Y). Por ejemplo:

plotyy(x1,y1,x2,y2,@loglog) % Control por función plotyy(x1,y1,x2,y2,'loglog') % Control se especifica en una cadena

plotyy(X1,Y1,X2,Y2,'FUN1','FUN2'): Utiliza FUN1(X1, Y1) para representar los datos en el eje izquierdo y FUN2(X2, Y2) para representar los datos en el eje derecho.[AX,H1,H2] = plotyy(...): AX es un vector que devuelve los controles de los ejes creados, donde AX(1) es el eje izquierdo y AX(2) es el eje derecho. H1 y H2 devuelve los controles los objetos gráficos.

Función plotyy.





Ejemplo: x = 0:0.01:20;y1 = 200*exp(-0.05*x).*sin(x);y2 = 0.8*exp(-0.5*x).*sin(10*x);[AX,H1,H2] = plotyy(x,y1,x,y2,'plot');set(get(AX(1),'Ylabel'),'String','Slow Decay'); set(get(AX(2),'Ylabel'),'String','Fast Decay');xlabel('Time (\musec)');title('Multiple Decay Rates');set(H1,'LineStyle','--');set(H2,'LineStyle',':');

Función plotyy.

0 2 4 6 8 10 12 14 16 18 20-200

-150

-100

-50

0

50

100

150

200

Slo

w D

ecay

Time (sec)

Multiple Decay Rates

0 2 4 6 8 10 12 14 16 18 20-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

Fast

Dec

ay





Gráfica log-log. La sintaxis es:loglog(Y)loglog(X1,Y1,...)loglog(X1,Y1,LineSpec,...)loglog(...,'PropertyName',Value,...)h = loglog(...)hlines = loglog('v6',...)

Ejemplo:x = logspace(-1,2);loglog(x,exp(x),'-s');grid on

Función loglog.

10-1 100 101 102100

1010

1020

1030

1040

1050





Semilogx genera grafica semilogarítmica en x, y semilogy lo realiza en y. La sintaxis es:semilogx(Y) ó semilogy(Y)semilogx(X1,Y1,...) semilogx(X1,Y1,LineSpec,...) semilogx(...,'PropertyName',Value,...) h = semilogx(...)

Ejemplo:x = 0:.1:10;semilogy(x,10.^x)

Función semilogx, semilogy).





Crea líneas. La sintaxis es:line(X,Y)line(X,Y,Z)line(X,Y,Z,'PropertyName',value,...)line('XData',x,'YData',y,'ZData',z,...)h = line(...)

Ejemplo:t = 0:pi/20:2*pi;hl1=plot(t,sin(t),'k');hl2= line(t+.06,sin(t),'LineWidth',4,'Color',[.8 .8 .8]);

Función line.

0 1 2 3 4 5 6 7-1

-0.8

-0.6

-0.4

-0.2

0

0.2

0.4

0.6

0.8

1





Gráfica en Coordenadas polares. La sintaxis es:polar(theta,rho)polar(theta,rho,LineSpec)h = polar(...)

Ejemplo:t = 0:.01:2*pi;polar(t,sin(2*t).*cos(2*t),'--r');

Función polar.





Gráfica en Coordenadas polares de forma directa. La sintaxis es:

ezpolar(fun)ezpolar(fun,[a,b])h = ezpolar(...)

Ejemplo1:ezpolar('1+cos(t)');

Ejemplo 2:fh = @(t) t.^2.*cos(t);ezpolar(fh);

Ejemplo 3:Si se define la función:function s = myfun(t,k1,k2)

s = sin(k1*t).*cos(k2*t);

Función ezpolar.

La gráfica se realiza por:ezpolar(@(t)myfun(t,2,3))





Gráficas Estadísticas en 2D.Otras instrucciones de dibujos en dos dimensiones:

– pie Representa gráficos de torta.– bar Gráficos de barras.– bar3 Gráficos de barras en 3 D.– stairs Gráficos de escalera.– hist Histogramas.– pareto Diagrama de pareto.– Complejos:

• compass• Feather

Nota: Ver detalles en el help de Matlab.





Gráfica lineal en 3D. La sintaxis es:

plot3(X1,Y1,Z1,...)plot3(X1,Y1,Z1,LineSpec,...)plot3(...,'PropertyName',Value,...)h = plot3(...)

Ejemplo:t = 0:pi/50:10*pi;plot3(sin(t),cos(t),t)grid onaxis square

Función plot3.

Una función equivalente pero que se usa en forma similar ezpolar, es:

ezplot3('s/2','2*s','s^2') ófh1 = @(s) s./2; fh2 = @(s) 2.*s; fh3 = @(s) s.^2;ezplot3(fh1,fh2,fh3)





Para graficas las superficies se debe siguiente procedimiento:Gráficas de Superficies.

• Evalúa la función Z=f(x,y).• La fila i-ésima de la matriz Z será f(x,y(i)), donde se mantiene

constante la variable x durante toda la evaluación de la fila.• Del mismo modo, la columna j-ésima corresponde a f(x(j),y),

donde se mantiene constante la variable y.Es decir, para graficar Z=f(x,y), usando los puntos dados por la coordenadas:

x = [-2 –1 0]; y=[1 2 3]

Z debe evaluarse en: Z = f(-2,1) f(-2,2) f(-2,3)f(-1,1) f(-1,2) f(-1,3)f(0,1) f(0,2) f(0,3)

Ctte. la variable x

Ctte. la variable y





Para evaluar los puntos que generan la matriz con los pares coordenados se unas la función meshgrid(x,y). La sintaxis es:

Gráficas de Superficies.

[X,Y] = meshgrid(x,y): Transforma el dominio dado por los vectores x,yen matrices rectangulares X,Y.[X,Y] = meshgrid(x): Similar que [X,Y] = meshgrid(x,x).[X,Y,Z] = meshgrid(x,y,z): Similar a la primera opción pero genera una matriz tridimensional.

Ejemplo:x=[-2 –1 0];y=[1 2 3;[X,Y]=meshgrid(x,y);

La salida es:

X =

-2 -1 0

-2 -1 0

-2 -1 0

Y =

1 1 1

2 2 2

3 3 3





Las función mesh permite grafica una función z=f(x,y) en un sistema coordenado 3D utilizado un modelo de alambre. Meshc, es similar, pero en el plano xy gráfica las líneas de contorno. Meshz proyecta el borde del dominio hacia el plano xy. La sintaxis es la siguiente:mesh(X,Y,Z) mesh(Z): La gráfica se realiza para X = 1:n, Y = 1:m, donde [m, n]=size(Z).mesh(...,C) mesh(...,'PropertyName',PropertyValue,...) mesh(axes_handles,...) meshc(...) meshz(...) h = mesh(...)hsurface = mesh('v6',...)hsurface = meshc('v6',...),

Funciones mesh, meshc, meshz y ezmesh.



-10-5

05

10

-10-5

0

510

-0.5

0

0.5

1

Sombrero Mexicano



Ejemplo: Grafique la función: en (-7.5, 7.5)

x=-7.5:.5:7.5;y=x;[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;mesh(X,Y,Z);title('Sombrero Mexicano');

Funciones mesh, meshc, meshz y ezmesh.





Funciones mesh, meshc, meshz y ezmesh.x=-7.5:.5:7.5;y=x;[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;meshc(X,Y,Z);title('Sombrero Mexicano');

-10-5

05

10

-10-5

0

510

-0.5

0

0.5

1

Sombrero Mexicano





Funciones mesh, meshc, meshz y ezmesh.x=-7.5:.5:7.5;y=x;[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;meshz(X,Y,Z);title('Sombrero Mexicano');

-10-5

05

10

-10-5

0

510

-0.5

0

0.5

1

Sombrero Mexicano





Funciones mesh, meshc, meshz y ezmesh.ezmesh es similar a ezplot y ezplot3.

Ejemplo: El grid es de 40x40

fh = @(x,y) x.*exp(-x.^2-y.^2);ezmesh(fh,40)colormap([0 0 1])

-20

2

-2-1

01

2

-0.4

-0.2

0

0.2

0.4

x

x exp(-x2-y2)

y



-10-5

05

10

-10-5

0

510

-0.5

0

0.5

1

Sombrero Mexicano



Función waterfall.

Ejemplo:x=-7.5:.5:7.5;y=x;[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;meshz(X,Y,Z);title('Sombrero Mexicano');

Similar a meshz pero las líneas de alambre se realizan sólo en una sola dirección, paralela a x o y. La sintaxis es:

waterfall(Z) waterfall(X,Y,Z) waterfall(...,C) waterfall(axes_handles,...) h = waterfall(...)



-10-5

05

10

-10-5

0

510

-0.5

0

0.5

1

Sombrero Mexicano



Funciones surf y surfc.

Ejemplo:x=-7.5:.5:7.5;y=x;[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;surf(X,Y,Z);title('Sombrero Mexicano');

Similar a mesh y meshc pero el se usa contornos de colores. La sintaxis es:surf(Z); surf(Z,C) surf(X,Y,Z); surf(X,Y,Z,C) surf(...,'PropertyName',Value) surf(axes_handles,...)surfc(...) h = surf(...)hsurface = surf('v6',...)hsurface = surfc('v6',...)



-10-5

05

10

-10-5

0

510

-0.5

0

0.5

1

Sombrero Mexicano



Funciones surf y surfc.Ejemplo:x=-7.5:.5:7.5;y=x;[X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;surfs(X,Y,Z);title('Sombrero Mexicano');

Nota: Hay otras funciones similares, como son:

surfl, ezsurf y ezsurfc





Función contour.Grafica las líneas de contornos (isocurvas) de una función. La sintaxis es:

contour(Z)contour(Z,n): n es el número de contornoscontour(Z,v): v es un vector que almacena los valores de contornos.contour(X,Y,Z)contour(X,Y,Z,n)contour(X,Y,Z,v)contour(...,LineSpec)contour(axes_handle,...)[C,h] = contour(...)[C,h] = contour('v6',...)





Función contour.Ejemplo:[C,h] = contour(peaks(20),10);

2 4 6 8 10 12 14 16 18 20

2

4

6

8

10

12

14

16

18

20



-5.0954

-3.8302

-3.8302

-2.5

651

-2.5

651

-2.5

651

-1.2999

-1.2999

-1.2999

-1.2999

-1.2

999

-0.034787

-0.034787

-0.0

3478

7

-0.034787

-0.03

4787

-0.0347871.2304

1.2304

1.2304

1.2304

1.2304

1.2304

1.2304

2.4955

2.4955

2.49

55

2.4955

2.4955

3.7606

3.7606

5.02585.0258

6.2909

2 4 6 8 10 12 14 16 18 20

2

4

6

8

10

12

14

16

18

20



Función contour.Ejemplo:[C,h] = contour(peaks(20),10);text_handle = clabel(C,h);





Función contour.Ejemplo:Z = peaks;[C,h] = contour(interp2(Z,4));text_handle = clabel(C,h);set(text_handle,'BackgroundColor',…[1 1 .6], 'Edgecolor',[.7 .7 .7])

-6

-4

-4

-2

-2

-2

-2 -2

0

0 0

00

0

0

2

2

2

2

2

2

2

4

4 4

6

68

100 200 300 400 500 600 700

100

200

300

400

500

600

700





Función contourf.Similar a contour pero hora la isocurvas son contornos de colores.Ejemplo:[C,h] = contourf(peaks(20),10);

2 4 6 8 10 12 14 16 18 20

2

4

6

8

10

12

14

16

18

20



-0.0

3478

7

-0.034787

-0.034787-0.034787

-0.034787

1.2304

1.2304

1.2304

1.2304

1.2304

2.4955

2.49

55

2.4955

-1.2999

-1.2999

-1.2

999

3.7606

3.7606

-2.5

651

-2.5

651

-1.2999

-1.2999

1.2304

1.2304

5.02585.0258

-3.8302

-3.8302

2.4955

6.2909

2.4955

-5.0954

-2.5

651

2 4 6 8 10 12 14 16 18 20

2

4

6

8

10

12

14

16

18

20



Función contourf.Ejemplo:[C,h] = contourf(peaks(20),10);text_handle = clabel(C,h);



-0.0

3478

7

-0.034787

-0.034787-0.034787

-0.034787

1.2304

1.2304

1.2304

1.23041.2304

2.4955

2.49

55

2.4955

-1.2999

-1.2999

-1.2

999

3.7606

3.7606

-2.5

651

-2.5

651

-1.2999

-1.2999

1.2304

1.2304

5.02585.0258

-3.8302

-3.8302

2.4955

6.2909

2.4955

-5.0954

-2.5

651

2 4 6 8 10 12 14 16 18 20

2

4

6

8

10

12

14

16

18

20

-6

-4

-2

0

2

4

6



Función contourf.Ejemplo:[C,h] = contourf(peaks(20),10);text_handle = clabel(C,h);colorbar

colorbar, se utiliza para activar la barra de colores de contorno que permite visualizar los valores de cada contorno de color.



-2-1

01

2

-2-1

0

12

-0.4

-0.2

0

0.2

0.4



Función contour3.Similar a contour pero en 3D: La sintaxis es:

contour3(Z)contour3(Z,n)contour3(Z,v)contour3(X,Y,Z)contour3(X,Y,Z,n)contour3(X,Y,Z,v)contour3(...,LineSpec)contour3(axes_handle,...)[C,h] = contour3(...)Ejemplo:

[X,Y] = meshgrid([-2:.25:2]);Z = X.*exp(-X.^2-Y.^2);contour3(X,Y,Z,30)



-2 -1.5 -1 -0.5 0 0.5 1 1.5 2-2

-1.5

-1

-0.5

0

0.5

1

1.5

2



Función quiver.Se usa para grafica campo de velocidad en forma vectorial: La sintaxis es:

quiver(x,y,u,v)quiver(u,v)quiver(...,scale)quiver(...,LineSpec)quiver(...,LineSpec,'filled')quiver(axes_handle,...)h = quiver(...)hlines = quiver('v6',...)

Ejemplo:[X,Y] = meshgrid(-2:.2:2);Z = X.*exp(-X.^2 - Y.^2);[DX,DY] = gradient(Z,.2,.2);contour(X,Y,Z)hold onquiver(X,Y,DX,DY)





Función quiver.Ejemplo:

[X,Y] = meshgrid(-2:.2:2);Z = X.*exp(-X.^2 - Y.^2);[DX,DY] = gradient(Z,.2,.2);contour(X,Y,Z)hold onquiver(X,Y,DX,DY);colormap hsvhold off

-2 -1.5 -1 -0.5 0 0.5 1 1.5 2-2

-1.5

-1

-0.5

0

0.5

1

1.5

2





Función quiver3.Similar a quiver pero en 3D: La sintaxis es:

quiver3(x,y,z,u,v,w)quiver3(z,u,v,w)quiver3(...,scale)quiver3(...,LineSpec)quiver3(...,LineSpec,'filled')quiver3(axes_handle,...)h = quiver3(...)





Función quiver3.Ejemplo:

[X,Y] = meshgrid(-2:0.25:2,-1:0.2:1);Z = X.* exp(-X.^2 - Y.^2);[U,V,W] = surfnorm(X,Y,Z);quiver3(X,Y,Z,U,V,W,0.5);hold onsurf(X,Y,Z);colormap hsvview(-35,45)axis ([-2 2 -1 1 -.6 .6])hold off

-2-1

01

2

-1

-0.5

0

0.5

1

-0.4

-0.2

0

0.2

0.4

0.6





Función streamlice.Dibuja las líneas de corriente en el plano: La sintaxis es:

streamslice(X,Y,Z,U,V,W,startx,starty,startz) streamslice(U,V,W,startx,starty,startz) streamslice(X,Y,U,V) streamslice(U,V) streamslice(...,density) streamslice(...,'arrowsmode') streamslice(...,'method') streamslice(axes_handle,...) h = streamslice(...) [vertices arrowvertices] = streamslice(...)



80 90 100 110 120 130

20

25

30

35

40

45

50

55



Función streamlice.Ejemplo: Grafique las líneas de corriente de la data del viento existente.

load winddaspect([1 1 1])streamslice(x,y,z,u,v,w,[],[],[5])axis tight





Función streamline.Grafica líneas de corriente desde una data en 2D o 3D. La sintaxis es:

streamline(X,Y,Z,U,V,W,startx,starty,startz) streamline(U,V,W,startx,starty,startz) streamline(XYZ) streamline(X,Y,U,V,startx,starty) streamline(U,V,startx,starty) streamline(XY) streamline(...,options) streamline(axes_handle,...) h = streamline(...)



6080

100120

140

0

20

40

60-5

0

5

10

15

20



Función streamline.Ejemplo:

load wind[sx,sy,sz] = meshgrid(80,20:10:50,0:5:15);h = streamline(x,y,z,u,v,w,sx,sy,sz);set(h,'Color','red')view(3)

2-graficas con matlab

Documents