Tutorial MATLAB

MATLAB: Conceptos Básicos y Programación Autor: Ing. Luis Rodríguez Ojeda Instituto de Ciencias Matemáticas Escuela Superior Politécnica del Litoral

Objetivo Proporcionar a los interesados los conocimientos básicos para usar el entorno de MATLAB y las facilidades para su programación.

Metodología Mediante explicaciones basadas en los ejemplos incluidos en este manual, el interesado puede adquirir en forma progresiva y autónoma los conocimientos básicos para utilizar MATLAB. Para progresar rápidamente, puede abrir dos ventanas en la pantalla de su computador, una con el programa MATLAB y otra con este manual, entonces puede copiar cada ejemplo del manual y pegarlo en la ventana de MATLAB para probarlo.

Introducción MATLAB (Matrix Laboratory) es un programa interactivo de uso general en Ciencias e Ingeniería La interacción se realiza mediante instrucciones (comandos), y también mediante funciones y programas (scripts). Los objetos básicos con los cuales opera MATLAB son matrices. La especificación y asignación de cada variable la realiza MATLAB en forma dinámica y eficiente, por lo que no son necesarias las declaraciones de variables de algunos lenguajes de computación convencionales. Características de MATLAB • • • • • • • •

Cálculo numérico rápido y con alta precisión Manejo simbólico Graficación y visualización avanzada Programación mediante un lenguaje de alto nivel Programación estructurada y orientada a objetos Soporte básico para diseño de interfaz gráfica Extensa biblioteca de funciones Aplicaciones especializadas para algunas ramas de las ciencias y la ingeniería (toolboxes)

Operación • Simple y eficiente • Interactivo • Sistema de ayuda en línea • Interacción con otros entornos

Una primera mirada a MATLAB

Cargue el programa MATLAB en su computador. El entorno de MATLAB está organizado mediante ventanas Command Window: es la ventana de ejecución de las instrucciones de MATLAB Command History: es el registro de los comandos que han sido ingresados. Workspace: contiene la descripción de las variables usadas en cada sesión. Inicialmente, se sugiere dejar activa únicamente la ventana de comandos, cerrando las otras ventanas con el ícono x. Para restaurarlas use la opción View de la barra de herramientas de MATLAB. Los ejemplos de este manual debe probarlos en la Ventana de Comandos El símbolo >> indica que el programa está listo para recibir sus instrucciones. Escriba cada ejemplo y presione la tecla de ingreso (Enter) para finalizar. MATLAB muestra el resultado inmediatamente, o un mensaje si hubo algún error. Al final de cada ejemplo aparece con letra azul una breve explicación. Por supuesto, usted no tiene que escribir en el computador estas explicaciones. >> exp(2)/3 ans = 2.4630

calcule y muestre inmediatamente el resultado

>> x=exp(2)/3;

calcule y asigne el resultado a la variable x

>> x

escriba la variable para conocer su contenido x= 2.4630

respuesta mostrada por Matlab, ans proviene de answer

respuesta mostrada por MATLAB

>> x=exp(2)/3 x= 2.4630

al omitir ; el resultado se muestra inmediatamente

>> y=2*x+1 y= 5.9260

puede usar el contenido de las variables

>> x=x+1 x= 3.4630

puede modificar el contenido de las variables

respuesta mostrada por MATLAB

respuesta mostrada por MATLAB

respuesta mostrada por MATLAB

Puede reutilizar comandos con las teclas del cursor ↑ ↓ >> x=exp(2)/3; y=2*x+1, z=3*x

. Pruébelo.

Puede escribir y ejecutar varios comandos en una misma línea

y= 5.9260

respuestas mostradas por MATLAB

z= 7.3891 El sistema de ayuda de MATLAB Ofrece una descripción detallada del uso de cada comando. >> help >> help ops >> help exp

despliega temas de ayuda despliega comandos de un tema. Ej. lista de operadores uso de un comando específico. Ej. función exponencial

Adicionalmente, presionando el ícono Help usted puede entrar al sistema de ayuda de MATLAB organizado por contenido, índice, búsqueda y demostraciones. Algunos ejemplos para apreciar el poder de MATLAB

1) Resuelva el sistema: 2x + 3y = 4 5x – 2y = 6 >> a=[2, 3; 5, -2]; >> b=[4; 6]; >> x=inv(a)*b; >> x

ingrese la matriz de coeficientes ingrese el vector columna de constantes obtenga la solución con la función para invertir muestre el vector solución

2) Grafique la función f(x)=sen(x) ex en el intervalo 0 ≤ x ≤ π >> ezplot('sin(x)*exp(x)', 0, pi); escriba la función entre comillas simples >> grid on; muestre cuadrículas en el gráfico

3) Resuelva la ecuación cúbica 5x3 + 2x2 - 3x + 1 = 0; >> a=[5, 2, -3, 1]; >> x=roots(a)

ingrese los coeficientes de la ecuación obtenga y muestre las raíces de la ecuación

4) Obtenga la solución de la ecuación diferencial ordinaria: y´-x-y = 0, y(0)=1 >> y=dsolve('Dy-x-y=0','y(0)=1', 'x'); defina la ecuación, condición y variable >> y muestre la solución analítica obtenida >> ezplot(y, 0, 2); grafique la solución para 0 ≤ x ≤ 2 >> grid on muestre cuadrículas

5) Integre la función f(x) = x sen(x) >> f=int('x*sin(x)') >> s=int('x*sin(x)',0,pi)

obtenga el resultado analítico integre entre 0 y π

6) Manejo simbólico de expresiones >> syms x >> y=x^3-8

defina x con tipo simbólico una expresión con x

>> t=factor(y) >> e=taylor(exp(x), 5);

factorar la espresión asignada a y expandir ex con 5 términos de la serie de Taylor

Símbolos especiales que usa MATLAB [ ] para definir vectores y matrices ( ) para definir precedencia en expresiones y para subíndices , para separar elementos de un vector use comas o espacios ; para separar filas y para evitar mostrar contenido de variables % para iniciar un comentario (programas y funciones) ... para continuar un comando en la siguiente línea

CÁLCULO NUMÉRICO Formatos de exhibición de números en la pantalla >> format long muestra 14 decimales >> x=exp(2) un ejemplo para visualizar >> format bank 2 decimales >> x >> format rat notación racional (fracciones) >> x >> format short e notación científica >> x >> format long e notación científica con 14 decimales >> x >> format + muestra signos +, , >> x >> format short 4 decimales (MATLAB lo usa por omisión) >> x >> format compact suprime líneas adicionales en la salida >> x >> format loose inserta líneas en blanco en la salida(recomendado) >> x >> format hex formato hexadecimal >> x >> vpa(sqrt(2), 100) variable precision arithmetic (muestra la raíz cuadrada de 2 con 100 dígitos >> format short regrese al formato normal de MATLAB Operadores aritméticos + - * / \ ^ () >> help ops

^ se usa para potenciación / es división a la derecha \ es división a la izquierda liste los operadores y caracteres especiales

Funciones matemáticas exp, log, sqrt, sin, cos, tan, ..... >> help elfun

liste las funciones matemáticas elementales

Practique con expresiones aritméticas >> x=sin(3)+3*exp(2) Operadores relacionales y lógicos < <= > >= == ~= & | ~ los tres últimos corresponden a: == representa al símbolo = ~= representa al símbolo ≠

∧ ∨

Practique con operadores relacionales y lógicos >> t=sin(2) < 0.8 & log(2) > 0.5 el resultado es un valor lógico (0 o 1) Símbolos numéricos especiales >> 2/0 Inf es el símbolo ∞ >> 0/0 NaN significa “Not A Number” (valor indeterminado) >> pi >> eps >> realmin >> realmax >> i

contiene la constante π es la precisión del tipo real en MATLAB el menor número real en MATLAB el mayor número real en MATLAB representa al símbolo complejo √ -1

Practique con números complejos >> x=3+2i asigne un número complejo >> t=2*x + 3 - 5i operación con números complejos >> y=exp(x) el resultado también es complejo >> y=log(-2) el referencial de MATLAB son los complejos Funciones adicionales para números complejos conj, real, imag, abs, angle, complex >> z=3+2i; >> t=conj(z) obtenga el conjugado Variables • No requieren ser declaradas • Su tipo depende del valor asignado • Pueden ser redefinidas • Sensible al tipo de letra (mayúsculas o minúsculas) • ans es la variable por omisión provista por MATLAB >> x=3 >> x='mensaje' >> syms x >> x=[2 7 4] >> x=2+3i >> x >> whos x >> disp(x)

x es de tipo real x es de tipo literal (use comillas simples) x es un símbolo x es un vector x es de tipo complejo muestre el contenido actual de la variable muestre el tipo actual de la variable muestre solamente el contenido



>> x=input('¿dato?'); >> exp(x)/3 >> ans >> y=2*ans

ingrese un valor para una variable desde el teclado la variable ans contiene el último resultado la puede usar

Algunos comandos del Sistema Operativo >> help general lista de comandos >> who lista las variables en uso >> whos lista las variables en uso y su descripción >> clear a b c >> clc >> pwd

clear borra variables. ejemplo despeja la ventana de comandos muestra cual es el directorio actual

>> cd c:\matlab\work >> dir

cd cambia la ruta del directorio actual ejemplo lista el contenido del directorio actual

>> save prueba

save almacena las variables en un archivo .mat ejemplo

>> load prueba

load carga variables y su contenido ejemplo

>> delete prueba.mat >> quit

delete elimina archivo ejemplo para terminar la sesión con MATLAB (no lo digite)

Comandos especiales >> date >> clock >> format rat >> clock >> format short Cadenas de caracteres >> x='Matematica'; >> x(4) >> t=x(2:5); >> n=length(x) >> c=strcat(x, t) >> help strfun Vectores y Matrices >> x=[3, -1, 4, 7, -2] >> x=[3 -1 4 7 -2] >> x(2)=5

fecha fecha hora, vea su uso con help. para visualizar la fecha con mas claridad vuelva al formato normal asignación de una cadena (use comillas simples) manejo de un carácter de la cadena, use un indice En MATLAB los índices se escriben entre paréntesis y son numerados desde 1 manejo de una subcadena, use: (inicio: final) longitud de la cadena concatenación de cadenas liste las funciones para cadenas asignación directa de un vector fila puede separar con comas o con espacios manejo de un componente del vector.

>> y=x(2: 4) >> t=[3; -1; 4; 5] >> t=x'

En MATLAB los índices se escriben entre paréntesis y son numerados desde 1 para asignar parte de un vector use (inicio: final) para asignar un vector columna use ; para obtener la transpuesta de un vector use ' x' es la transpuesta del vector x

>> y=[3, x, -6, 7] >> y=[2:1:10]

puede asignar un vector usando otro vector puede asignar un vector mediante una secuencia En MATLAB las secuencias se escriben: valor inicial : incremento : valor final si el incremento es 1 puede omitirlo

>> y=[2, 5, 4, ... 7, -3] >> x=[3, 5, 2, 0] >> y=2*x >> y=exp(x)

Para continuar en la siguiente línea use ... escriba la continuación de la línea anterior

>> a=[6 3 ; 5 1]

asignación directa de una matriz 2x2 separe elementos con espacios o comas separe filas con punto y coma

>> a(2,1)

manejo de los componentes de una matriz con índices numerados desde 1: (fila, columna)

>> a=[2, -3; 5, 1; 0, 7] >> x=[7, 3] >> a=[x; x] >> b=[5, 6] >> c=[a; b] >> d=[a, b'] >> x=c(1, :) >> x=c(: ,1) >> c(:,2)=[ ]

una matriz 3x2

Matrices especiales >> a=ones(5) >> a=ones(3,5) >> a=zeros(4,5) >> a=eye(5) >> a=magic(4) >> a=hilb(5) >> x=[2, 5, 3, 7]; >> a=vander(x) >> a=[ ]

puede realizar operaciones escalares o crear vectores con funciones

una matriz 2x2 c es una matriz aumentada 3x2 c es una matriz aumentada 2x3 asigne a x la primera fila de c asigne a x la primera columna de c elimine la segunda columna de c matriz 5x5 iniciada con unos matriz 3x5 iniciada con unos matriz 4x5 iniciada con ceros matriz identidad 5x5 cuadrado mágico 4x4 matriz de Hilberth 5x5 matriz de Vandermonde 4x4 usando un vector matriz nula

Generación de números aleatorios >> x=rand genera un número aleatorio entre 0 y 1

>> a=rand(5) >> b=rand(4,5) >> d=fix(rand*10)+1

genera una matriz 5x5 con números aleatorios genera una matriz 4x5 con números aleatorios transfomación para obtener un entero aleatorio entre 1 y 10

Una matriz puede componerse con otras matrices >> a=rand(3); matriz 3x3 con números aleatorios >> b=[5 3 9]; vector de tres componentes >> e=diag(b); matriz 3x3 con b en la diagonal >> c=eye(3); matriz identidad 3x3 >> d=zeros(3); matriz con ceros 3x3 >> t=[a e; c d] matriz compuesta 9x9 Elementos de vectores y matrices pueden manejarse con otro vector o matriz >> x=[ 8 7 9 5 6]; >> p=[2 4 1]; vector para direccionar al vector x >> t=x(p) t contiene los elementos 2, 4 y.1 del vector x >> a=[4 7 3; 5 7 8; 6 0 9]; >> p=[1 3]; vector para direccionar las filas de la matriz a >> q=[2 3]; vector para direccionar las columnas de la matriz a >> t=a(p, q) t contiene las filas 1 y 3, columnas 2 y 3 de a Operaciones con matrices >> a=[3, 2; 1, 4]; >> b=[8, 6; 5, 7]; >> c=a'

transpuesta de a

>> c=2*a >> c=a+b >> c=a*b >> c=a.*b

producto de un escalar por matriz suma de matrices producto de matrices producto elemento por elemento de matrices para operar elemento a elemento use un punto antes del operador

>> c=a^2 >> c=a.^2 >> c=a==b

matriz al cuadrado, equivale a: a*a cada elemento de la matriz a, elevar al cuadrado compare igualdad entre matrices (de igual tamaño) el resultado es una matriz binaria (ceros y unos) compare si dos matrices no son iguales el resultado es una matriz binaria (ceros y unos) compare si cada elemento de a es mayor a 3 el resultado es una matriz binaria (ceros y unos)

>> c=a~=b >> c=a>3

Funciones para matrices >> x=[-2, 0, 6, 5]; un vector para los ejemplos >> a=[1, 2, 3; 4, 5, 6;7, 8, 9]; una matriz para los ejemplos >> n=length(x) longitud del vector x >> [n,m]=size(a) tamaño de la matriz a: el resultado es un vector >> n número de filas: 3 >> m número de columnas: 3

>> exist('c') >> isempty(a) >> any(x) >> any(a) >> t=find(x) >> t=find(x>3) >> [f,c]=find(a)

chequea si la variable existe chequea si un vector o matriz está vacío determina si el vector contiene algún valor no cero igual que arriba, pero por columnas de la matriz obtiene índices de elementos del vector no ceros obtiene los índices de cada elemento > 3 obtiene los índices de filas y columnas de la matriz cuyos elementos son no ceros

>> t=dot(x, x) >> k=rank(a) >> t=trace(a) >> d=det(a) >> b=inv(a) >> h=norm(a, 1) >> h=norm(a, inf) >> h=norm(x, inf) >> c=cond(a) >> t=diag(a) >> t=diag(x) >> t=rot90(a) >> t=fliplr(a) >> t=tril(a) >> t=triu(a) >> b=[5,-1; 3, 4; 2, 7]; >> b=reshape(b, 2, 3)

producto punto entre dos vectores rango de a traza de a determinante de a inversa de a norma de columna de la matriz a norma de fila de la matriz a norma de fila o columna del vector x número de condición de la matriz a vector con la diagonal de la matriz a matriz con x en la diagonal rote a 90 grados (sentido opuesto al reloj) voltee horizontalmente la matriz a obtenga la matriz triangular inferior de a obtenga la matriz triangular superior de a

>> [t,s]=lu(a) >> t >> s >> t*s >> t=cov(a) >> e=eig(a) >> p=poly(a) >> r=roots(ans) >> help matfun

descomposición triangular de a en las matrices triangulares t y s tales que t*s es igual que a

reconfigura la matriz b de 3x2 a 2x3

se obtiene la matriz a matriz de covarianza de a valores propios de a polinomio característico de a valores propios de a liste las funciones para matrices

Ingreso de puntos desde la pantalla con el mouse >> ezplot('sin(x)'); ejemplo para tomar puntos desde un gráfico >> grid on >> [x,y]=ginput(5); ingrese 5 puntos desde la pantalla . Presione el botón del mouse para ingresar cada punto >> x observe las abscisas >> y y las ordenadas ingresadas >> plot(x, y, 'o') grafique los puntos ingresados Funciones adicionales para manejo de datos con vectores y matrices >> x=[2, 5, 4]; un vector >> a=[5,-1; 3, 4; 2, 7]; una matriz

>> t=max(x) >> v=max(a)

el mayor valor del vector x el mayor valor por columnas de la matriz a

>> t=sum(x) >> v=sum(a)

suma de componentes suma de componentes por columnas

>> t=prod(x) >> v=prod(a)

producto escalar producto escalar por columnas

>> t=cumsum(x) >> v=cumsum(a)

suma acumulada suma acumulada por columnas

>> t=cumprod(x) >> v=cumprod(a)

producto acumulado

>> t=mean(x) >> v=mean(a)

media aritmética

>> t=median(x) >> v=median(a)

mediana

>> t=std(x) >> v=std(a)

desviación estándar

>> t=sort(x) >> v=sort(a) >> t=dsort(x)

ordenamiento ascendente

>> bar(x) >> bar(a)

diagrama de barras

>> hist(x) >> stairs(x) >> pie(x) >> pie3(x) >> v=[0,1,0] >> pie3(x,v)

histograma dibuja x mediante escalones gráfico tipo pastel pastel en relieve vector para extraer sectores del pastel gráfico tipo pastel con un sector separado

Polinomios >> a=[2, -3, 0, 5}, >> y=polyval(a,4) >> x=roots(a) >> t=polyval(a, x(1)) >> p=poly(x) >> b=[3, 4, -2]; >> c=conv(a,b) >> [c, r]=deconv(a,b); >> c >> r >> x=[2 3 5 7 8];

ordenamiento descendente

define el polinomio 2x3 – 3x2 + 5 evaluación del polinomio con un valor obtenga un vector con raíces (reales y complejas) verifique una raíz producto de todas las raíces define el polinomio 3x2+ 4x - 2 producto de polinomios división de polinomios cociente residuo abscisas de puntos (x,y)

>> y=[3.2 4.1 5.8 6.4 6.3]; ordenadas de los puntos >> z=3.2; valor para interpolar, z puede ser un vector >> u=interp1(x,y,z,’linear’) resultado de la interpolación lineal >> u=spline(x,y,z) interpolación con un trazador cúbico >> a=polyfit(x, y, 2); polinomio de mínimos cuadrados de grado 2 >> a el vector a contiene los coeficientes Practique: Obtenga y grafique el polinomio de interpolación, la recta de mínimos cuadrados y el trazador cúbico para un conjunto de datos dados >> x=[1 2 4 5 7]; >> y=[5 3 6 7 4]; >> plot(x,y,'o') >> grid on >> hold on

cinco puntos (x, y) para el ejemplo

>> a=polyfit(x,y,4); >> a >> z=[1: 0.1: 7]; >> p=polyval(a,z); >> plot(z,p)

polinomio de interpolación, 5 puntos: grado 4 coeficientes a(1)x^4 + a(2)x^3 + a(3)x^2 + ... puntos para evaluar el polinomio evalúe el polinomio con z obtenga puntos p grafique el polinomio de interpolación

>> b=polyfit(x,y,1); >> b >> t=[1 7]; >> q=polyval(b,t); >> plot(t,q,'r')

recta de mínimos cuadrados (grado 1) coeficientes de la recta: b(1)x + b(2) puntos extremos de la recta (abscisas) obtenga las ordenadas respectivas de la recta grafique la recta en color rojo

>> s=spline(x,y,z); >> plot(z,s,'g') >> hold off

evalúe con z el trazador cúbico y obtenga s grafique el trazador cúbico con verde deshabilite la superposición de gráficos

grafique los datos con círculos poner cuadrículas superponer los siguientes gráficos

MANEJO SIMBÓLICO >> syms x; >> 2*x+3*x >> a=[x 5; 3*x 4]; >> t=inv(a) >> t=solve(2*x^2+3*x-2) >> f=3*x^2+5*x; >> t=factor(f) >> s=expand(t) >> e=taylor(exp(x)) >> syms y; >> f=2*x^3+3*y^2 >> g=diff(f,x) >> u=int(f,x)

definición de variable tipo simbólico suma algebraica matriz con elementos símbolos su inversa también contiene símbolos solución de una ecuación definición simbólica de una función factorar la expresión expandirla expansión con la serie de Taylor una función de dos variables derivada parcial integrar en x

Funciones adicionales con expresiones simbólicas

>> f='2*t+1'; definición de una función en forma literal >> t=3; >> y=eval(f) evaluación de la función >> [a,b] = solve('a^2 + a*b - b = 3','a^2 - 4*b - 5 = 0') resuelve un sistema de dos ecuaciones no lineales >> a=double(a) para expresarlas en forma real >> b=double(b) >> f='exp(x)-pi*x'; >> x=solve(f) >> x=double(x) cambia la solución simbólica a real x= 0.5538 resultados de MATLAB 1.6385 >> x=fzero(f,2) x= 1.6385

solución de una ecuación con un valor inicial

>> x=fzero(f,[1,2]) x= 1.6385

solución usando un rango para la raiz

>> u=diff(f) >> v=int(f) >> r=int(f, 0, 2) >> g='x*exp(-x)'; >> r=int(g, 0, Inf);

diferenciación integración analítica integración entre límites

resultado de MATLAB

resultado de MATLAB

integral impropia

solución analítica y gráfico de una EDO de primer orden con una condición en el inicio >> y=dsolve('Dy=(x-y)/x','y(0)=0','x') >> ezplot(y,0,2); >> grid on EDO de segundo orden con condiciones en el inicio y gráfico de la solución >> y=dsolve('D2y+Dy+2*y-x-3=0','y(0)=0,Dy(0)=1','x') >> ezplot(y,0,2); >> grid on; EDO de segundo orden con condiciones en los bordes >> y=dsolve('D2y-Dy+2*y-5*x-3=0','y(0)=0,y(1)=2','x') >> ezplot(y,0,2); >> grid on >> f=’2*sin(x)*exp(x)’; >> [x,y]=fminbnd(f,-2,2) >> ezplot(f,-2,2)

encuentra un mínimo local de f en 1 ≤ x ≤ 4

Ejercicio: Escriba directamente las instrucciones necesarias para encontrar el valor del radio x de un cilindro de 1000 cc de capacidad, de tal manera que el valor del área sea el mínimo: Primer enfoque: 1) Escriba una función f en términos del radio x 2) Grafique f con ezplot. Localice el intervalo para el mínimo de f(x) 3) Use la función fminbnd para obtener el mínimo x: radio, h: altura >> f='2*pi*x*1000/(pi*x^2)+2*pi*x^2'; >> ezplot(f,0,10) >> grid >> x=fminbnd(f,4,6) x= 5.4193 >> area=eval(f) area = 553.5810 Segundo enfoque 1) Derive f y obtenga la función a minimizar g. 2) Grafique g con ezplot. Localice el intervalo de la raíz de g(x)=0 3) Use la función fzero para obtener la raíz 4) Use la funcion solve para obtener la raíz x: radio, h: altura >> g=diff(f) g= -2000/x^2+4*pi*x >> x=fzero(char(g),[4,6]) x= 5.4193 >> x=solve(g) x= [ 5/pi*4^(1/3)*(pi^2)^(1/3)] [ -5/2/pi*4^(1/3)*(pi^2)^(1/3)+5/2*i*3^(1/2)/pi*4^(1/3)*(pi^2)^(1/3)] [ -5/2/pi*4^(1/3)*(pi^2)^(1/3)-5/2*i*3^(1/2)/pi*4^(1/3)*(pi^2)^(1/3)] >> x=double(x) x= 5.4193 -2.7096 + 4.6932i -2.7096 - 4.6932i Funciones especiales para medir eficiencia de algoritmos >> tic; Inicia cronómetro >> toc; muestra el tiempo transcurrido

>> tic; a=inv(rand(500, 500)); toc tiempo utilizado en invertir una matriz 500x500

GRAFICACIÓN Gráfico de funciones de una variable >> f='exp(x)-3*x'; función para el ejemplo (use comillas simples) >> ezplot(f, 0, 2) función simple para graficar en un dominio dado >> x=[0: 0.1: 2*pi]; puntos para evaluar una función >> y=sin(x); puntos de la función seno >> plot(x,y); función para graficar la función con línea contínua >> plot(x,y,'o') gráfico con círculos. Puede elegir: o . * + x ->> plot(x,y,'r') cambiar a color rojo. Puede elegir r,b,y,m,g,w,k >> plot(x,y,'og') grafique con círculos verdes. >> grid on dibuje las cuadrículas >> title('seno de x') >> gtext('seno de x') >> xlabel('X') >> ylabel('Y')

incluya un título en el gráfico posicione el texto en el gráfico con el mouse rotule el eje horizontal rotule el eje vertical

>> c=[0, 2*pi, -2, 2] >> axis(c) >> hold on >> hold off

defina la región para el gráfico

>> figure(1) >> subplot(2,3,1) >> clf(1) >> clf >> x=[0:0.1:10}; >> y=exp(x); >> semilogx(x,y) >> semilogy(x,y) >> loglog(x,y) >> grid on

puede tener varias figuras abiertas puede dividir una figura en subgráficos borra el gráfico 1 borre todos los gráficos

>> a=0:0.01:2*pi; >> r=sin(3*a); >> polar(a, r);

superponer subsecuentes gráficos deshabilitar opción anterior

grafique en escalas logarítmicas

'rosa' de 3 pétalos grafique en coordenadas polares

Gráfico de funciones de dos variables >> a=[1 3 2; 5 3 7; 4 5 2]; una matriz 3x3 >> mesh(a); grafique sus elementos como puntos sobre el plano.

>> x=-2:0.1:2; >> y=-3:0.1:3;

El siguiente ejemplo es la referencia para graficar funciones de dos variable dominio de la función para el ejemplo: graficar z = x2 – y2

>> [u,v]=meshgrid(x,y); >> z=u.^2 – v.^2; >> mesh(x, y, z) >> contour(x, y, z) >> surfc(x, y, z) >> surf(x, y, z)

u, v: matrices q' contienen cada par ordenado x,y puntos de la función z = x2 – y2 gráfico de malla gráfico de contorno gráfico de superficie y contorno gráfico de superficie

>> xlabel('X') >> title('Silla de montar') >> colormap(copper); >> shading interp;

rotulación de eje x; también con ylabel, zlabel título para el gráfico color del gráfico; también: gray, jet, pink suavizado del gráfico (¡ tiene que verlo !)

Gráfico de funciones implícitas >> x=0: 0.1: 4; intervalos para el siguiente ejemplo: >> y=0: 0.1: 4; grafique la circunferencia (x - 2)2 + (y - 2)2 = 4 >> [u,v]=meshgrid(x,y); u, v: matrices q' contienen cada par ordenado x,y >> f=(u-2).^2+ (v-2).^2 - 4; puntos de la función f = (x - 2)2 + (y - 2)2 - 4 >> contour(x, y, f, [0,0]) gráfico del primer nivel del contorno de f (circunferencia en el plano X -Y) >> axis square gráfico con ejes con escala proporcional >> title('Círcunferencia') título para el gráfico

PROGRAMACIÓN CON MATLAB Archivos de comandos Para crear un archivo de comandos (programa, o script), seleccione en la barra de herramientas de MATLAB la opción: File → New → M-file o presione en el ícono respectivo. Se abre una ventana de edición En la ventana de edición escriba los comandos y almacénelos en un archivo con algún nombre. Puede incluir comentarios con el símbolo % Ej. Escriba en la ventana de edición las instrucciones para graficar sen(x) entre 0 y 2π x=0:0.1:2*pi; y=sin(x); plot(x,y); grid on Guarde el archivo con algún nombre, ejemplo: prueba Para ejecutar el programa escriba en la ventana de comandos >> prueba Para editar un archivo de comandos, seleccione en la barra de herramientas de MATLAB la opción: File → Open o seleccione el ícono respectivo.

En el archivo abierto en la ventana de edición haga los cambios y guárdelo nuevamente con el mismo nombre.

Estructuras de control de flujo en MATLAB Instrucciones de Entrada y Salida Ingreso de un dato desde el teclado: Ej. x = input('ingrese un dato ');

variable=input(‘mensaje’);

Salida de un resultado a la pantalla: Ej. x=exp(2); disp(x);

disp(valor)

Salida de mas de un resultado a pantalla: Ej. x=2^7; y=sqrt(pi); disk([x, y]);

disp([valor, valor, …]);

Salida de resultados formateados a pantalla: fprintf(‘formatos’,variables) Ej. x=2^7; y=sqrt(pi); fprintf('%d %f’',x,y); Puede especificar cantidad de columnas y decimales: Ej. x=2^7; y=sqrt(pi); fprintf('%5d %8.3f’',x,y); Otras especificaciones de formato puede verlas con help fprintf

Decisiones (instrucción if) if condición instrucciones end

if condición instrucciones else instrucciones end

Ej. Escriba y almacene el siguiente programa para mostrar el mayor entre dos datos: a=input('ingrese el primer dato '); b=input('ingrese el segundo dato ');

if a>b m=a; else m=b; end disp(m); Guárdelo con el nombre prueba y úselo desde la ventana de comandos: >> prueba ingrese el primer dato 5 interacción con MATLAB ingrese el segundo dato 8 8

Decisiones múltiples (instrucción switch) Ej. Escriba el siguiente programa para instrumentar la definición: x =1  3,  y = 2x + 1, x = 2,3,4  ex, otro x  x=input('dato '); switch x case 1, y=3; case {2,3,4}, y=2*x+1; otherwise y=exp(x); end disp(y); Almacénelo con algún nombre. Ejemplo prueba2 Para usarlo escriba en la ventana de comandos >> prueba2

Repetición condicionada al inicio (instrucción while) while condición instrucciones end Ej. Sume los n primeros términos de la serie armónica: n=input('cantidad de terminos '); s=0; while n>0 s=s+1/n; n=n-1;

end disp(s); Almacénelo con algún nombre. Ejemplo prueba3 Para usarlo escriba en la ventana de comandos >> prueba3

Repetición condicionada a una secuencia (instrucción for) for variable=inicio: incremento: final instrucciones end Ej. Sume los n primeros términos de la serie armónica: n=input('cantidad de terminos '); s=0; for i=1:n s=s+1/i; end disp(s); Almacénelo con algún nombre. Ejemplo prueba4 Para usarlo escriba en la ventana de comandos >> prueba4

Puede interrumpir una repetición (instrucción break) Ej. Lea n datos. Calcule y muestre la raíz cuadrada. Pero si entra un valor negativo, muestre un mensaje y termine n=input(‘cantidad de datos ‘); for i=1:n x=input(‘ingrese siguiente dato ‘); if x<0 disp(‘Error’); break; else r=sqrt(x); disp([x,r]); end end Para ver la descripción de las estructuras del lenguaje de MATLAB, escriba >> help lang

Funciones en MATLAB

En general una función en los lenguajes de programación es un conjunto de instrucciones que se escriben separadamente del programa y que realizan alguna tarea especificada. Los usuarios pueden definir funciones y agregarlas a las funciones propias de MATLAB. El mecanismo usual para trasmitir datos a las funciones es mediante una lista de variables que se denominan parámetros. Sin embargo, a diferencia de los programas, las variables que se usan dentro de una función, no están disponibles fuera de ella, a menos que se use una declaración explícita y que se verá mas adelante. Declaración de una función en MATLAB function variable = nombre (parámetros) instrucciones variable contendrá el resultado que entrega la función parámetros son variable que reciben los datos que entran a la función nombre identifica a la función instrucciones se incluyen en la función según la tarea especificada Las funciones se escriben en la ventana de edición de Matlab y se las almacena en alguna carpeta. Es conveniente que el nombre asignado sea igual al nombre usado en la declaración de la función. El uso de una función es similar al uso de las funciones comunes de MATLAB. El nombre debe coincidir con el nombre asignado, aunque los parámetros pueden tener nombres diferentes, pero su uso debe ser coherente. Ej. Escriba una función para elegir el mayor entre dos números Abra un documento nuevo en la ventana de edición y escriba: function m = mayor(a, b) if a>b m = a; else m = b; end m mayor a, b

es la variable que entrega el resultado es el nombre de la función son los parámetros que ingresan los datos a la función

Almacene esta función en el disco con el nombre mayor Suponer que quiere escoger el mayor entre eπ y πe. Escriba en la ventana de comandos:

>> a = exp(pi); >> b = pi^exp(1); >> m = mayor(a, b) 23.1407

(respuesta que muestra MATLAB)

Los nombres de las variables pueden ser diferentes: >> x = exp(pi); >> y = pi^exp(1); >> t = mayor(x, y)

23.1407

(respuesta que muestra MATLAB)

Ej. Escriba una función que reciba un número y determine si es un número primo. El resultado que entrega la función será 1 o 0 según corresponda; function p = primo( x ) c = 0; for d = 1: x if mod(x, d) == 0 c = c + 1; end end if c > 2 p = 0; else p = 1; end Guarde la función en el disco con el nombre primo Pruebe la función desde la ventana de comandos >> x = 25; >> p = primo(x)

0

(resultado que muestra MATLAB)

>> x = 43; >> p = primo(x)

1

(resultado que muestra MATLAB)

Escriba en una nueva ventana de edición un programa que use la función primo para encontrar todos los números primos menores a 20: for x = 1: 20 if primo(x) == 1 disp(x); end end Almacene su programa en el disco con el nombre prueba En la ventana de comandos pruebe su programa:

>> prueba 1 2 3 5 7 11 13 17 19 Ejercicios.-

(resultados mostrados por MATLAB)

1) Escriba una función que reciba 3 números y entregue como resultado el mayor con el siguiente encabezado: function m=mayor(a,b,c) Escriba un programa que lea las tres dimensiones de un bloque, calcule el valor de las diagonales de sus caras y use la función anterior para encontrar el valor de la mayor diagonal 2) Escriba una función que muestre cual es el valor de la suma de los cubos de los primeros n números naturales con el siguiente encabezado: function s=suma(n) Escriba n programa que use la función suma y la función primo anterior para determinar si las suma de los cubos de los 7 primeros números naturales es un número primo Una función puede entregar más de un resultado Las variables que entregan los resultados deben definirse entre [ ] Ej. Escriba una función que entregue el área y el volumen de un cilindro dados su radio (r) y su altura (h) function [area, vol] = cilindro(r, h) area = 2*pi*r*h + 2*pi*r^2; vol = pi*r^2*h; Escriba y almacene la función con el nombre cilindro. Use la función para calcular el área y el volumen de una lata de cilíndrica que tiene un diámetro de 10cm y una altura de 12cm Escriba en la ventana de comandos: >> r = 5; >> h = 12;

>> [a, v] = cilindro(r,h); >> a >> v MATLAB mostrará los resultados almacenados en a y en v Las variables definidas dentro de una función son locales, es decir que a diferencia de los programas, no son visibles fuera de la función Ej. Escriba la función: function x=fn(a, b) c = a + b; x = 2*c; Almacene con el nombre fn y úsela desde la ventana de comandos: >> a = 3; >> b = 5; >> t = fn(a, b) t = 16 >> c

(resultado que muestra MATLAB) (intentamos conocer el valor de c en la función)

??? Undefined function or variable 'c'. (mensaje de error de MATLAB) Compare con lo que ocurre si escribe un programa en vez de la función; a = input(‘ingrese dato ‘); b = input(‘ingrese dato ‘); c = a + b; x = 2*c; disp(x); Almacene con el nombre prueba y active el programa: >> prueba ingrese dato 3 ingrese dato 5 16 >> c c=8

(interacción para ingreso de datos) (resultado que muestra MATLAB) (la variable c puede ser utilizada)

Es posible hacer que las variables de una función sean visibles fuera de su ámbito, mediante la declaración global Ej. Modifique la función fn para que la variable c sea visible: function x=fn(a, b) global c; c = a + b;

x = 2*c; Almacene con el nombre fn y use la función: >> a = 3; >> b = 5; >> t = fn(a, b) t = 16 >> c c=8

(resultado que muestra MATLAB) (intentamos conocer el valor de c en la función) (la variable c está disponible ahora)

Una función puede no necesitar parámetros Ej. Escriba una función que lea y valide un entero entre 1 y 5 function n=entero x=0; while x==0 n=input(‘ingrese un entero entre 1 y 5 ‘); if n>0 & n<6 x=1; end end Una función puede no entregar resultados ni usar parámetros Ej. Escriba una función que imprima un menú function menú disp(‘1) ingresar’); disp(‘2) borrar’); disp(‘3) salir’); para usa esta función escriba >> menu Una función puede recibir como parámetros vectores o matrices. Ej. Escriba una función que reciba un vector y entregue el promedio del valor de sus elementos. function p=prom(x) n=length(x); s=0; for i=1:n s=s+x(i); end p=s/n;

Para usar esta función debe definir el vector antes de llamar a la función. La función determina la longitud del vector con la función length >> x=[2 7 3 5 4 7 6]; >> t=prom(x) t = 4.8571

(es el resultado que muestra MATLAB)

Una función puede entregar como resultado un vector o una matriz Ej. Escriba una función que entregue un vector de longitud n conteniendo números aleatorios enteros con valor entre 1 y 6: function d=dados(n) for i=1:n d(i)=fix(rand*6+1); end Para usar esta función debe enviar un valor para el parámetro n: >> t=dados(5) t=6 3 4 3 2

(es el vector resultante que entrega MATLAB)

Una función puede recibir y entregar vectores o matrices Ej. Escriba una función que reciba dos vectores A, B y entregue un tercer vector que contenga los elementos que están en ambos vectores: function C=interseccion(A,B) n=length(A); m=length(B); k=1; for i=1:n for j=1:m if A(i) == B(j) C(k) = A(i); k = k + 1; end end end Para usar esta función debe definir los vectores que entran. Recuerde que pueden tener nombres diferentes a los que usa la función: >> A=[2 7 5 4 3 8]; >> B=[7 1 3 9 0]; >> C=interseccion(A,B)

C =

7

3

(Es el vector resultante que entrega MATLAB)

Si la salida de una función es antes del final, puede usar el comando return Ej. Escriba una función para determinar si los elementos de un vector están en orden creciente: function t=orden(x) t=0; n=length(x); for i=1:n-1 if x(i) > x(i+1) return; end end t=1; Un programa puede llamar a funciones Ej. Escriba una función para eliminar espacios intermedios de una frase: function x=compactar(f) n=length(f); x=''; for i=1:n if f(i) ~= ' ' x = strcat(x, f(i)); end end Ahora escriba un programa que lea una frase, use la función compactar para eliminar los espacios intermedios, y luego muestre un mensaje en caso de que sea simétrica: sus caracteres opuestos son iguales f=input('ingrese una frase '); f=compactar(f); n=length(f); sim=1; for i=1:n/2 if f(i) ~= f(n-i+1) sim=0; end end if sim == 1 disp('la frase es simetrica'); else disp('la frase no es simetrica'); end

Probamos este programa suponiendo que lo hemos almacenado con el nombre prueba: >> prueba ingrese una frase 'anita lava la tina'; la frase es simetrica

(dato que ingresamos) (resultado de MATLAB)

Una función puede llamarse a si misma Estas funciones de denominan recursivas Ej. Use la siguiente definición recursiva para calcular el máximo común divisor entre dos números enteros:  mcd(a − b, b), a > b  mcd(a, b) =  mcd(a, b - a), b > a a, a=b  Escriba y almacene una función para instrumentar esta definición: function c=mcd(a, b) if a>b c=mcd(a-b, b); else if b>a c=mcd(a, b-a); else c=a; end end Use la función: >> x=mcd(36, 48) Para desplegar un mensaje de error y terminar la ejecución use error Ej. if d<0 error('valor incorrecto'); end

---------------------------------------------------------------------------------Practique: Funciones del Módulo de Estadística normpdf: distribución de probabilidad normal >> x=4:0.1:10; >> p=normpdf(x,7,1); >> plot(x,p)

normcdf: Distribución de probabilidad normal acumulada Calcule la probabilidad que una variable aleatoria distribuida normalmente, con media 10 y desviación estándar 2, tenga un valor entre 11.4 y 13.5 >> p=normcdf(13.5, 10, 2) - normcdf(11.4, 10, 2); p= 0.2019 norminv: Inversa de la distribución normal acumulada Encuentre un intervalo de 99% de probabilidad para una variable con distribución normal con media 2 y desviación estándar 1 >> p=[0.005 0.995]; >> x=norminv(p,2,0.5) x= 0.7121 3.2879 Para chequear >> p=normcdf(3.2879,2,0.5)-normcdf(0.7121,2,0.5) p= 0.9900

Prueba de Hipótesis [h,p,c]=ztest(x,u,sn); >> [h,p,c]=ztest(2,2.05,0.04) h= 0 p= 0.2113 c= 1.9216 2.0784 Practique: Una función que entrega un resultado analítico Escriba y almacene una función para obtener el polinomio de interpolación en forma analítica con la fórmula de Lagrange: function p=Lagrange(x,y) p=0; syms t; n=length(x); for i=1:n

L=1; for j=1:n if i~=j L=L*(t-x(j))/(x(i)-x(j)); end end p=p+L*y(i); p=simplify(p); end Use la función desde la línea de comandos >> x=[3, 5, 6, 8, 9]; >> y=[3.5, 4.7, 4.9, 5.2, 6.8]; >> p=Lagrange(x,y) >> plot(x,y,'o'), grid, hold >> ezplot(p,x(1),x(5))

Manejo de archivos en MATLAB El siguiente ejemplo completo es una referencia para conocer las instrucciones de MATLAB para manejo de archivos Una aplicación de manejo de un archivo con datos de estudiantes. Los datos que se almacenan en el archivo incluyen: código del estudiante (entero de 5 dígitos) código de la materia (entero de 5 dígitos) calificación obtenida por el estudiante en la materia (entero de 3 dígitos) El programa debe pedir un nombre para el archivo, abrirlo para agregar datos y ofrecer las siguientes opciones mediante un menú: 1) Ingreso de un nuevo dato 2) Dado el código de un estudiante, mostrar sus materias y calificaciones 3) Dado el código de una materia, mostrar los estudiantes y su calificación 4) Salir Variables que son usadas en el programa e codigo de estudiante m codigo de materia c calificacion arch nombre del archivo dado por el usuario

function control global arch; clc; arch=input('ingrese el nombre del archivo '); op=0; while op~=4 op=menu; switch op case 1, ingresar; case 2, consultae; case 3, consultam;

end end function op=menu clc; disp('1) ingreso'); disp('2) consulta por estudiante'); disp('3) consulta por materia'); disp('4) salir'); op=input('ingrese una opcion '); function ingresar global arch; f=fopen(arch,'a+'); e=input('codigo de estudiante '); m=input('codigo de materia '); c=input('calificacion '); fprintf(f,'%5d%5d%3d',e,m,c); fclose(f); function consultae global arch; f=fopen(arch,'r'); if f<0 disp('archivo no existe '); pause; return; end x=input('codigo de estudiante '); while ~feof(f) e=fscanf(f,'%d',1); m=fscanf(f,'%d',1); c=fscanf(f,'%d',1); if e==x fprintf('\n codigo de materia %5d',m); fprintf('\n calificacion %3d',c); end end pause; fclose(f); function consultam global arch; f=fopen(arch,'r'); if f<0 disp('archivo no existe '); pause; return; end x=input('codigo de materia ');

while ~feof(f) e=fscanf(f,'%d',1); m=fscanf(f,'%d',1); c=fscanf(f,'%d',1); if m==x fprintf('\n codigo de estudiante %5d',e); fprintf('\n calificacion %3d',c); end end pause; fclose(f);

INTERACCIÓN CON OTROS ENTORNOS Interacción con EXCEL 1) Importar una tabla de datos desde Excel a una matriz en MATLAB

a) En Excel cree la tabla y almacénela con formato tipo texto delimitado con tabulaciones. Elija algún nombre. Ejemplo T.txt

b) En MATLAB cargue la tabla T y úsela como una matriz: >> load T.txt; >> A=T 2) Exportar una matriz de datos desde MATLAB a una tabla en Excel

a) En MATLAB cree una matriz y almacénela con el comando save con el siguiente formato. Elija los nombres. Ejemplo A: nombre de la matriz en MATLAB T: nombre para la tabla almacenada >> save T A -ascii

b) En Excel abra el archivo T y úselo como una tabla de datos -------------------------------------------------------------------------------------------------Este manual es para uso y distribución libre. Comentarios y sugerencias envíelos a: [email protected] --------------------------------------------------------------------------------------------------