Guiadematlab
Guiadematlab
Guiadematlab
DE MATLAB
ÍNDICE
INTRODUCCIÓN.........................................................................................................4
CARACTERÍSTICAS BÁSICAS................................................................................4
EL ESPACIO DE TRABAJO DE MATLAB........................................................................................4
MATEMÁTICA SENCILLA.................................................................................................................5
ALMACENAR Y RECUPERAR DATOS.............................................................................................7
FORMATOS DE VISUALIZACIÓN DE NÚMEROS..........................................................................7
ACERCA DE LAS VARIABLES..........................................................................................................8
OTRAS CARACTERÍSTICAS BÁSICAS.............................................................................................9
AYUDA EN LÍNEA.......................................................................................................9
VECTORES Y MATRICES.......................................................................................12
CÓMO DEFINIRLOS.......................................................................................................................12
DIRECCIONAMIENTO DE ELEMETOS DE VECTORES Y MATRICES......................................13
CONSTRUCCIÓN ABREVIADA DE ALGUNOS VECTORES........................................................15
CONSTRUCCIÓN DE ALGUNAS MATRICES................................................................................16
OPERACIONES BÁSICAS CON MATRICES..................................................................................17
FUNCIONES PARA OPERAR CON VECTORES............................................................................18
FUNCIONES PARA EL ANÁLISIS DE MATRICES........................................................................19
OTRAS OPERACIONES CON MATRICES......................................................................................20
TEXTO.........................................................................................................................21
HIPERMATRICES.....................................................................................................22
CÓMO DEFINIRLAS.......................................................................................................................22
OPERACIONES CON HIPERMATRICES.......................................................................................23
ESTRUCTURAS..........................................................................................................24
CÓMO DEFINIRLAS.......................................................................................................................24
OPERAR CON ESTRUCTURAS......................................................................................................25
2
Manual de MATLAB Servicios Informáticos U.C.M.
GRÁFICAS...................................................................................................................30
2-D....................................................................................................................................................30
3-D....................................................................................................................................................33
PROGRAMACIÓN DE MATLAB............................................................................41
SENTENCIA FOR.............................................................................................................................41
SENTENCIA WHILE........................................................................................................................41
SENTENCIA IF.................................................................................................................................42
SENTENCIA BREAK........................................................................................................................43
SENTENCIA CONTINUE.................................................................................................................43
FUNCIONES EN M-ARCHIVOS..............................................................................44
ANÁLISIS DE DATOS...............................................................................................45
POLINOMIOS.............................................................................................................47
RAÍCES.............................................................................................................................................47
OTRAS CARACTERÍSTICAS............................................................................................................48
ANÁLISIS NUMÉRICO.............................................................................................49
REPRESENTACIÓN GRÁFICA.......................................................................................................49
OTRAS CARACTERÍSTICAS............................................................................................................50
EJERCICIOS PROPUESTOS...................................................................................54
3
INTRODUCCIÓN
MATLAB es el nombre abreviado de “MATriz LABoratory”. Es un programa para realizar cálculos
numéricos con vectores y matrices, y por tanto se puede trabajar también con números escalares
(tanto reales como complejos), con cadenas de caracteres y con otras estructuras de información más
complejas.
Matlab es un lenguaje de alto rendimiento para cálculos técnicos, es al mismo tiempo un entorno y un
lenguaje de programación. Uno de sus puntos fuertes es que permite construir nuestras propias
herramientas reutilizables. Podemos crear fácilmente nuestras propias funciones y programas
especiales (conocidos como M-archivos) en código Matlab, los podemos agrupar en Toolbox
(también llamadas librerías): colección especializada de M-archivos para trabajar en clases
particulares de problemas.
Matlab, a parte del cálculo matricial y álgebra lineal, también puede manejar polinomios, funciones,
ecuaciones diferenciales ordinarias, gráficos …
CARACTERÍSTICAS BÁSICAS
Nada más abrir Matlab (podemos hacerlo pinchando en el icono que aparece en el escritorio o en su
defecto en Inicio->Todos los programas) aparecerá una pantalla como la siguiente:
Todas las sentencias que vamos a utilizar las escribiremos en la ventana Command Window (ventana
de comandos). Es la ventana de mayor tamaño.
Si queremos información acerca de las variables que estamos utilizando en Matlab podemos verlas en
la ventana Workspace (espacio de trabajo) o usar:
Para ver esta ventana tenemos que pinchar en la pestaña que tienen este nombre. Está en la parte
superior izquierda:
Si lo que queremos es conocer el valor que tiene una variable lo hacemos escribiendo el nombre de la
variable y pulsando Intro.
Para recordar órdenes previas usamos las flechas del teclado ↑ y ↓. También podemos verlas en la
ventana Command History, ventana situada en la parte inferior izquierda:
MATEMÁTICA SENCILLA
Si queremos que Matlab evalúe la línea pero que no escriba la respuesta, basta escribir punto y coma
(;) al final de la sentencia.
Si la sentencia es demasiado larga para que quepa en una sola línea podemos poner tres puntos (…)
seguido de la tecla Intro para indicar que continúa en la línea siguiente.
Ejemplos:
>> a / b
ans =
1.7500
>> a ^ b
ans =
2401
>> 5 * a
ans =
35
>> whos % da una lista de las variables usadas más completa que la anterior
Matlab permite guardar y cargar datos de los archivos del computador. En el menú File, la opción
Save Workspace as… guarda todas las variables actuales y Import Data… carga variables de un
espacio de trabajo guardado previamente.
Otra forma sería guardar el estado de una sesión de trabajo con el comando save antes de salir:
>> save
>> load
Matlab no cambia la representación interna de un número cuando se escogen distintos formatos, sólo
se modifica la forma de visualizarlo.
Las variables son sensibles a las mayúsculas, deben comenzar siempre con una letra, no pueden
contener espacios en blanco y pueden nombrarse hasta con 63 caracteres (en versiones anteriores no
permitía tantos caracteres). Si se nombra una variable con más de 63 caracteres truncará el nombre de
dicha variable.
Tecleando clear podemos borrar todas las variables del espacio de trabajo, pero no borra lo de las
demás ventanas, es decir, no desaparece lo que hay escrito en la ventana de comandos.
Tecleando clc borramos lo que hay en la ventana de comandos pero no borra las variables de la
memoria del espacio de trabajo.
Algunos comandos de Matlab nos facilitan información sobre la fecha, como clock, date o calendar.
>> clock % año mes día hora minutos y segundos, en este orden
ans =
1.0e+003 *
2.0060 0.0110 0.0140 0.0120 0.0190 0.0437
Los comentarios se escriben después del símbolo de tanto por ciento (%), de este modo todo lo que
se escriba a continuación en la misma línea no será leído por Matlab. Podemos colocar varias órdenes
en una línea si se separan correctamente, puede ser:
Para cerrar Matlab podemos hacerlo tecleando quit, cerrando con el aspa típico de Windows,
entrando en File->Exit Matlab o con las teclas Ctrl+Q.
AYUDA EN LÍNEA
Matlab proporciona asistencia de varios modos.
Por último con la orden lookfor <palabra>, busca en todas las primeras líneas de las ayudas de los
temas de Matlab y devuelve aquellos que contienen la palabra clave que hemos escrito. No es
necesario que la palabra clave sea una orden de Matlab.
APROXIMACIONES
(con x escalar, vector o matriz, pero redondearía en cada caso los elemento individualmente)
Ejemplo:
Ejemplos:
Ejemplos:
VECTORES Y MATRICES
CÓMO DEFINIRLOS
Para crear un vector introducimos los valores deseados separados por espacios (o comas) todo ello
entre corchetes []. Si lo que queremos es crear una matriz lo hacemos de forma análoga pero
separando las filas con puntos y comas (;).
Generalmente usaremos letras mayúsculas cuando nombremos a las matrices y minúsculas para
vectores y escalares. Esto no es imprescindible y Matlab no lo exige, pero resulta útil.
Ejemplos:
>> y = [2,1,3,7] % es otro vector, los elementos los separamos con comas
y=
2 1 3 7
>> z = [0 1 2,3 4,5] % es otro vector, da igual separar los elementos por comas o espacios
z=
0 1 2 3 4 5
Para acceder a los elementos individuales de un vector lo haremos utilizando subíndices, así x(n)
sería el n-ésimo elemento del vector x. Si queremos acceder al último podemos indicarlo usando end
como subíndice.
>> x = [5 7 -2 4 -6];
>> x (2) % segundo elemento del vector x
ans =
7
Para acceder a un bloque de elementos a la vez, se usa la notación de dos puntos (:), así x (m:n) nos
da todos los elementos desde el m-ésimo hasta el n-ésimo del vector x.
Si introducimos un número entre el primero y el segundo también separado por dos puntos (:) se
mostrarán los elementos del primero al último indicado, incrementados según el número que aparece
en el centro (o decrementados si el número es negativo).
Otra forma de obtener un conjunto concreto de elementos del vector es indicando entre corchetes []
las posiciones de los elementos que queremos obtener poniendo paréntesis fuera de los corchetes.
Para acceder a los elementos de una matriz necesitamos dar dos valores, el primero indica la fila y el
segundo la columna.
>> A = [1 2 3; 4 5 6];
>> A (2,1) % elemento de la matriz que está en la fila 2 y en la columna 1
ans =
4
Si queremos que escriba toda una fila usaremos los dos puntos para indicar que queremos todos los
elementos.
Al igual que con los vectores podemos indicar que escriba una serie de filas o columnas, la manera de
hacerlo sería muy parecido.
Matlab tiene además otra forma de identificar cada elemento de una matriz, de modo que podemos
acceder a un elemento de una matriz indicando sólo un valor y no dos, pero debemos saber que el
orden elegido por Matlab es por columnas así los elementos de la matriz A serían denominados:
A(1)A(3)A(5)
A(2)A(4)A(6)
Ejemplo:
>> A (5) % accede al elemento 5 de la matriz, es decir, igual que si escribiéramos A (1,3)
ans =
3
Pero es preferible para evitar confusiones trabajar con los elementos de las matrices indicando la fila
y la columna correspondiente.
CONSTRUCCIÓN ABREVIADA DE ALGUNOS VECTORES
A parte de definir un vector introduciendo cada uno de sus elementos, también podemos crearlo
haciendo uso de las siguientes sentencias:
linspace (a,b,c) genera un vector linealmente espaciado entre los valores a y b con c elementos.
linspace (a,b) genera un vector linealmente espaciado entre los valores a y b con 100 elementos.
logspace (a,b,c) genera un vector logarítmicamente espaciado entre los valores 10^a y 10^b con c
elementos.
logspace (a,b) genera un vector logarítmicamente espaciado entre los valores 10^a y 10^b con 50
elementos.
Ejemplos:
>> logspace (0,2,4) % genera un vector logarítmicamente espaciado entre 10^0 y 10^2 con 4
elementos
ans =
1.0000 4.6416 21.5443 100.0000
CONSTRUCCIÓN DE ALGUNAS MATRICES
Al igual que pasa con los vectores, existen unas sentencias que nos ayudan a crear más rápidamente
algunas matrices que Matlab ya tiene predefinidas (m y n deben tomar valores naturales):
rand (n) crea una matriz cuadrada n x n de números aleatorios con distribución uniforme (0,1).
rand (m,n) crea una matriz m x n de números aleatorios con distribución uniforme (0,1).
randn (n) crea una matriz cuadrada n x n de números aleatorios con distribución normal (0,1).
randn (m,n) crea una matriz m x n de números aleatorios con distribución normal (0,1).
eye (n) crea una matriz cuadrada n x n de unos en la diagonal y ceros el resto.
magic (n) crea una matriz cuadrada n x n de enteros de modo que sumen lo mismo las filas y las
columnas.
hilb (n) crea una matriz cuadrada n x n de Hilbert, es decir, los elementos (i,j) responden a la
expresión (1/(i+j-1)).
invhilb (n) crea una matriz cuadrada n x n que es la inversa de la matriz de Hilbert.
Ejemplos:
A= B= C=
1 2 1 1 1.0000 + 1.0000i 2.0000 + 2.0000i
3 4 0 1 3.0000 + 1.0000i 4.0000 + 7.0000i
Ejemplos:
Ejemplos:
>> v = [1 2 3];
>> diag (v) % crea una matriz diagonal a partir del vector v
ans =
1 0 0
0 2 0
0 0 3
>> A = [1 2 3 4; 7 8 9 2; 2 4 6 8]
A=
1 2 3 4
7 8 9 2
2 4 6 8
>> size (A) % devuelve las dimensiones de la matriz como un vector fila
ans =
3 4
19
>> length (A) % devuelve la mayor de las dos dimensiones de la matriz
ans =
4
TEXTO
Una cadena de caracteres es texto rodeado por comillas simples (') y se manejan como vectores filas.
Se direccionan y manipulan igual que los vectores. Son posibles las operaciones matemáticas sobre
cadenas. Una vez hecha una operación matemática sobre una cadena, ésta se ve como un vector de
números en ASCII.
Para ver la representación ASCII de una cadena, podemos utilizar las funciones abs, double o
sumamos cero. Para restaurarla y verla de nuevo como cadena de caracteres, usamos la función
setstr. Si queremos cambiar a minúsculas añadiremos la diferencia entre 'a' y 'A'.
>> a = 'casa'; b = 'gato'; % a y b son cadenas de caracteres (se manejarán como vectores)
>> a + b
ans =
202 194 231 208
>> double (a) % otra tercera forma de ver la representación ASCII de la cadena
ans =
99 97 115 97
>> abs ('a') – abs ('A') % calculamos la diferencia entre minúsculas y mayúsculas
ans =
32
>> disp ('escribe esto') % escribe el texto que vaya entre las comillas
escribe esto
HIPERMATRICES
CÓMO DEFINIRLAS
Matlab permite trabajar con matrices de más de dos dimensiones. Los elementos de una hipermatriz
pueden ser números, caracteres, estructuras y vectores o matrices de celdas. Las funciones que operan
con matrices de más de dos dimensiones son análogas a las funciones vistas anteriormente aunque
con algunas diferencias, por ejemplo, a la hora de definirlas:
>> HM(:,:,1) = [1 2 3; 4 5 6]; % definimos la primera capa
>> HM(:,:,2) = [7 8 9; 10 11 12] % definimos la segunda capa
HM(:,:,1) =
1 2 3
4 5 6
HM(:,:,2) =
7 8 9
10 11 12
Algunas funciones para generar matrices admiten más de dos subíndices y pueden ser utilizadas para
generar hipermatrices como rand, randn, zeros y ones, también se pueden emplear con
hipermatrices las funciones size y reshape entre otras. La función cat permite concatenar matrices
según las distintas “dimensiones”.
Ejemplos:
>> A = zeros (2,3); B = ones (2,3); % definimos dos matrices de las mismas dimensiones
>> cat (1,A,B) % las concatena una debajo de la otra
ans =
0 0 0
0 0 0
1 1 1
1 1 1
>> cat (3,A,B) % las concatena como distintas capas de una hipermatriz
ans(:,:,1) =
0 0 0
0 0 0
a ns(:,:,2) =
1 1 1
1 1 1
1. Las funciones que operan sobre escalares, como sin, cos, etc., se aplican sobre hipermatrices
elemento a elemento (igual que ocurre al aplicarlas sobre vectores y matrices).
2. Las funciones que operan sobre vectores, como sum, max, etc., se aplican a matrices e
hipermatrices según la primera dimensión, resultando un array de una dimensión inferior.
3. Las funciones matriciales propias del álgebra lineal, como det, inv, etc., no se pueden aplicar
a hipermatrices, para aplicarlas habría que extraer las matrices correspondientes.
ESTRUCTURAS
CÓMO DEFINIRLAS
Es una agrupación de datos de tipo diferente bajo un mismo nombre. A los datos les llamamos
campos. No hace falta definir previamente el modelo de la estructura, podemos ir creando los
distintos campos uno a uno o bien con el comando struct, donde los nombres de los campos se
escriben entre apóstrofos (') seguidos del valor que se les quiere asignar.
Ejemplos:
nombre: 'Pablo'
apellido1: 'Fernández'
apellido2: 'García'
edad: 15
nombre: 'Fermín'
apellido1: 'Martínez'
apellido2: 'Gil'
edad: 16
>> alumno (1) % nos devuelve los datos del primer elemento del vector de la estructura
ans =
nombre: 'Pablo'
apellido1: 'fernández'
apellido2: 'García'
edad: 15
>> alumno (2) % nos devuelve los datos del segundo elemento del vector de la estructura
ans =
nombre: 'Fermín'
apellido1: 'Martínez'
apellido2: 'Gil'
edad: 16
>> alumno.nombre % escribe los datos de todos los campo nombre de la estructura en orden
ans =
Pablo
ans =
Fermín
Ejemplos:
>> isfield (alumno,'nombre') % devuelve 1 por ser cierto que nombre es un campo de alumno
ans =
1
>> isstruct (alumno) % devuelve 1 porque es cierto que alumno es una estructura
ans =
1
CÓMO DEFINIRLOS
Un vector de celdas es un vector cuyos elementos son cada uno de ellos una variable de cualquier
tipo. En todo vector sus elementos pueden ser números o cadenas de caracteres, pero en un vector de
celdas el primer elemento puede ser un número, el segundo una matriz, el tercero una estructura, etc.
>> celda (1) = {[0 1 2]}; % creamos un vector de celdas definiendo celda a celda
>> celda (2) = {'cadena de caracteres'};
>> celda (3) = {eye(2)};
>> celda (4) = {-7};
>> celda
celda =
[1x3 double] [1x20 char] [2x2 double] [-7]
>> cel {1} = [0 1 2]; % creamos otro vector de celdas definiendo celda a celda de forma distinta
>> cel {2} = 'cadena de caracteres';
>> cel {3} = eye (2);
>> cel {4} = -7;
>> cel
cel =
[1x3 double] [1x20 char] [2x2 double] [-7]
Si queremos crear una matriz o una hipermatriz de celdas se haría de forma similar.
Ejemplos:
OPERADORES LÓGICOS
Además de los operadores relacionales y lógicos básicos anteriores, Matlab proporciona una serie de
funciones relacionales y lógicas adicionales que incluyen:
Podemos ver muchos más casos pero todos serían similares: ischar, isempty, isequal, isfloat,
isinteger, islogical, isnumeric, isprime, isreal, isscalar, isspace, …
Ejemplos:
>> r1 = a<6 % pregunta si a es menor que 6, devuelve 1 cuando es verdadero y 0 cuando es falso
r1 =
1 1 1 1 1 0 0 0 0
>> r4 = (a>b)&(b>-3)% pregunta si a>b y b>-3, devuelve 1 cuando es verdadero y 0 cuando es falso
r4 =
0 0 1 1 1 1 1 0 0
>> isnan (c) % pregunta cuando c es NaN, devuelve 1 cuando es verdadero y 0 cuando es falso
ans =
0 0 0 0 1 0
>> isinf (c) % pregunta cuando c es Inf, devuelve 1 cuando es verdadero y 0 cuando es falso
ans =
1 0 0 0 0 0
>> isfinite (c) % pregunta cuando c es finito, devuelve 1 cuando es verdadero y 0 cuando es falso
ans =
0 1 1 1 0 1
GRÁFICAS
2-D
La orden plot genera una gráfica. Los argumentos deben ser vectores de la misma longitud.
Ejemplo:
Para obtener una información más detallada se recomienda utilizar la ayuda de Matlab:
Una ventana gráfica se puede dividir en m particiones horizontales y en n verticales, de modo que
cada subventana tiene sus propios ejes, y para hacer esto vamos a usar subplot (m,n,p) donde p
indica la subdivisión que se convierte en activa.
>> x = 1:360; y1 = sind (x); y2 = cosd (x); y3 = exp (x); y4 = exp (-x);
>> subplot (2,2,1), plot (x,y1), title ('seno')
>> subplot (2,2,2), plot (x,y2), title ('coseno')
>> subplot (2,2,3), plot (x,y3), title ('exponencial')
>> subplot (2,2,4), plot (x,y4), title ('-exponencial')
Para volver al modo por defecto basta escribir: subplot (1,1,1).
Para dibujar polígonos podemos usar la función plot pero teniendo en cuenta que el último punto de
ambos vectores deben coincidir para que la gráfica quede cerrada. Pero si lo que queremos es que
quede coloreado todo el interior del polígono debemos usar mejor la función fill, tiene tres
argumentos, los dos vectores que forman los puntos y un tercer argumento para indicar el color.
3-D
Gráficos de línea:
También podemos crear gráficas en 3 dimensiones, se trata de extender la orden de plot (2-D) a plot3
(3-D) donde el formato será igual pero los datos estarán en tripletes:
>> x = -720:720; y = sind (x); z = cosd (x);
Podemos hacer girar la gráfica usando de la barra de herramientas el botón o hacerla más grande o
más pequeña con . Al igual que ocurría con las gráficas en dos dimensiones podemos nombrar
los ejes o hacer modificaciones entrando en opciones con el botón .
Superficie de malla:
La orden [X,Y]=meshgrid(x,y) crea una matriz X cuyas filas son copias del vector x y una matriz Y
cuyas columnas son copias del vector y. Para generar la gráfica de malla se usa la orden
mesh(X,Y,Z), mesh acepta un argumento opcional para controlar los colores. También puede tomar
una matriz simple como argumento: mesh(Z).
Ejemplo:
Gráfica de superficie:
Es similar a la gráfica de malla, pero aquí se rellenan los espacios entre líneas. La orden que usamos
es surf con los mismos argumentos que para mesh.
Ejemplo:
Ejemplo:
Ejemplo:
Manipulación de gráficos:
Función Colores
Jet
HSV
Hot
Cool
Spring
Summer
Autumn
Winter
Gray
Bone
Copper
Pink
Lines
Ejemplos:
SENTENCIA FOR
Un bloque for en cada iteración asigna a la variable la columna i-ésima de la expresión y ejecuta las
órdenes. En la práctica las expresiones suelen ser del tipo escalar:escalar en cuyo caso las columnas
son escalares.
Ejemplo:
SENTENCIA WHILE
Un bloque while ejecuta las órdenes mientras todos los elementos de la expresión sean verdaderos.
while <expresión>
<orden>
<orden>
…
<orden>
end
Ejemplo:
>> a=3;
>> while a < 5
disp ('a es menor que 5 ya que vale')
disp (a)
a = a + 1;
end
SENTENCIA IF
Un bloque if puede escribirse de varias maneras distintas. Lo que hace es evaluar una expresión
lógica y si es cierta ejecuta las órdenes que encuentre antes del end.
if <expresión>
<órdenes evaluadas si la expresión es verdadera>
end
Puede que nos interese que en caso de no ejecutar dicha orden ejecute otra distinta. Esto se lo
indicaremos usando else dentro del bloque.
if <expresión>
<órdenes evaluadas si la expresión es verdadera>
else
<órdenes evaluadas si la expresión es falsa>
end
Si queremos dar una estructura mucho más completa, usaremos la más general donde sólo se evalúan
las órdenes asociadas con la primera expresión verdadera de todas. En cuanto la evalúe deja de leer el
resto y se dirige directamente al end.
if <expresión1>
<órdenes evaluadas si la expresión1 es verdadera>
elseif <expresión2>
<órdenes evaluadas si la expresión2 es verdadera>
elseif <expresión3>
<órdenes evaluadas si la expresión3 es verdadera>
elseif
…
…
else
<órdenes evaluadas si ninguna otra expresión es verdadera>
end
Ejemplo:
>> b = 2;
>> if b == 0 % ponemos == porque no es una asignación sino una expresión lógica
disp ('b vale 0')
elseif b == 1
disp ('b vale 1')
elseif b == 2
disp ('b vale 2')
elseif b == 3
disp ('b vale 3')
else
disp ('b no vale ni 0 ni 1 ni 2 ni 3')
end
b vale 2 % es lo que devuelve por pantalla
SENTENCIA BREAK
Si queremos que en un momento dado termine la ejecución de un bucle for o un bucle while
usaremos break.
SENTENCIA CONTINUE
La sentencia continue hace que se pase inmediatamente a la siguiente iteración del bucle for o del
bucle while saltando todas las órdenes que hay entre el continue y el fin del bucle en la iteración
actual.
Ejemplo:
Podemos mezclar en un programa varias sentencias de este estilo. Aquí podemos ver un programa
que escribe por pantalla los primos del 1 al 100 usando las sentencias if, while y for.
Las funciones permiten definir funciones análogas a las de Matlab, con su nombre, argumentos y
valores de salida. La primera línea que no sea comentario debe empezar por la palabra function,
seguida por los valores de salida (entre corchetes [ ] y separados por comas si hay más de uno), el
signo igual (=) y el nombre de la función seguido de los argumentos (entre paréntesis ( ) y separados
por comas):
function [a,b,c] = nombre_función (x,y,z)
En las líneas siguientes escribimos los argumentos de salida a partir de los de entrada. El nombre de
la función y el nombre del archivo deben ser idénticos y no empezar por cifra sino por letra.
Todas las variables dentro de una función se aíslan del espacio de trabajo de Matlab. Las únicas
conexiones entre las variables dentro de una función y el espacio de trabajo de Matlab son las
variables de entrada y salida.
El número de variables de entrada pasadas a una función está disponible dentro de la función en la
variable nargin y el número de variables de salida solicitadas cuando una función es llamada, está
disponible dentro de la función en la variable nargout.
Debemos tener siempre en cuenta que los argumentos pueden ser vectores, luego si queremos que las
operaciones se hagan elemento a elemento y no vectorialmente debemos usar el punto.
Ejemplo:
ANÁLISIS DE DATOS
Matlab ejecuta análisis estadístico sobre conjuntos de datos. Estos conjuntos de datos se almacenan en
matrices orientadas por columnas. Matlab incluye, entre otras, las siguientes funciones estadísticas:
>> X = [5 7 9 2 9; 3 1 7 5 1; 3 9 2 7 5; 1 5 5 1 8]
X=
5 7 9 2 9
3 1 7 5 1
3 9 2 7 5
1 5 5 1 8
>> Y = [5 7 9 2 9 3 1 7 5 1 3 9 2 7 5 1 5 5 1 8];
>> tabulate (Y) % tabla de frecuencias generada a partir de una serie de valores
Value Count Percent
1 4 20.00%
2 2 10.00%
3 2 10.00%
4 0 0.00%
5 5 25.00%
6 0 0.00%
7 3 15.00%
8 1 5.00%
9 3 15.00%
>> range (X) % rango de cada columna (diferencia entre el máximo y el mínimo)
ans =
4 8 7 6 8
POLINOMIOS
RAÍCES
Un polinomio se representa por un vector fila con sus coeficientes en orden descendiente, no
debemos olvidar colocar los términos con coeficiente nulo.
Así por ejemplo si queremos indicar el polinomio 5x4 + 2x2 – x + 7 escribiríamos [5 0 2 -1 7].
Matlab trabaja con los polinomios como vectores fila y con las raíces como vectores columnas.
Ejemplos:
>> poly (ans) % devuelve el polinomio generado por esas cuatro raíces
ans =
1.0000 -9.0000 13.0000 9.0000 -14.0000
OTRAS CARACTERÍSTICAS
Ejemplos:
>> p = [1 2 7];
>> q = [1 3 6]; % polinomios
>> polyval (p, [0 1 2; -1 -2 -3; 4 0 7] ) % igual pero toma los valores de una matriz
ans =
7 10 15
6 7 10
31 7 70
ANÁLISIS NUMÉRICO
REPRESENTACIÓN GRÁFICA
Existe la función fplot que evalúa la función que se desea representar en la gráfica de salida. Como
entrada, necesita conocer el nombre de la función como una cadena de caracteres y el rango de
representación como un vector de dos elementos: fplot ('nombre', [ valor min, valor max] ).
Ejemplo:
Matlab incorpora una serie de funciones para resolver ecuaciones diferenciales ordinarias. Si se trata
de un problema rígido deberíamos usar: ode15s, ode23s, ode23t u ode23tb, si por el contrario se trata
de un problema sin rigidez: ode113, ode 23 y ode45. Para saber más de estas funciones consultar la
ayuda de Matlab.
Ejemplos:
>> feval ('cos', [0 pi/3 pi] ) % para evaluar en varios puntos debemos darlo como un vector
ans =
1.0000 0.5000 -1.0000
>> feval (@cos, [0 pi/3 pi] ) % es lo mismo que lo anterior, da igual comillas que el @
ans =
1.0000 0.5000 -1.0000
>> fminbnd (@sind,0,360) % valor del dominio donde la función toma el mínimo
ans =
270.0000
>> fzero ('sind',100) % el valor más próximo a 100 donde la función seno vale cero
ans =
180
function nombre (el nombre de la función debe coincidir con el nombre del fichero)
Ahora debemos situarnos en el directorio donde tengamos el fichero que queremos transformar
usando el comando cd, por ejemplo:
Lo que debemos escribir a continuación es el comando mcc seguido de –m y el nombre del fichero:
Con esto nos aparecerá en el mismo directorio donde estamos un ejecutable con el mismo nombre.
También aparecerán una carpeta y varios archivos.
Ejemplo:
Creamos un fichero que va a ser una función que a su vez va a llamar a otras dos funciones que
también hemos creado nosotros:
Fichero algebra.m:
% algebra
function algebra
x = input ('Escribe un número: ');
y = input ('Escribe otro número: ');
disp ('La suma es...')
suma (x,y)
disp ('La resta es...')
resta (x,y)
pause %para que no se cierre la ventana automáticamente al ejecutarse
Fichero suma.m:
% suma
function m = suma (tt,xx)
m = tt + xx;
Fichero resta.m:
% resta
function m = resta (tt,xx)
m = tt - xx;
Para ello usaremos un botón que se encuentra en la ventana workspace. Vemos que al situar el
ratón sobre él aparece un letrero diciendo para lo que sirve (import data):
Al pinchar en él se abre una ventana. Debemos localizar el fichero que queremos importar y pinchar
en el botón donde pone Abrir. Aparecerá una nueva ventana similar a ésta:
(Aquí estamos importando un fichero de datos .xls de Excel con el nombre Libro1 pero al importarlo
lo renombra como Hoja1 ya que el fichero en cuestión tenía 3 hojas, pero sólo la Hoja1 tenía datos)
Pinchamos en el botón Finish y ya tendremos el fichero de datos convertido en una matriz en Matlab.
Si queremos visualizarla sólo tenemos que llamarla ya que se almacenará con el mismo nombre.
Ejemplo:
Queremos importar un fichero de datos de Excel con dos hojas (la Hoja3 está vacía):
Al importar el fichero nos aparece la ventana siguiente (sólo aparecen dos matrices porque la Hoja3
está vacía):
Pinchamos en Finish y aceptamos. Si queremos ver cómo ha guardado los datos basta llamar a las
matrices con el nombre que hayan sido almacenadas. (Los espacios en blanco los ha guardado como
NaN).
>> Hoja1
Hoja1 =
10 6 -1
5 15 0
10 2 6
-8 8 4
>> Hoja2
Hoja2 =
-2 4 8
6 NaN 4
1 3 NaN
Para exportar una matriz podemos convertirla en texto haciendo lo siguiente:
Primero escribimos:
A=
0 1
2 3
4 5
6 7
8 9
Escribimos:
Ahora para abrirlo buscamos el documento llamado my_data.out que se encuentra en:
C:\Archivos de programa\MATLAB\R2006b\work
Haz que aparezcan por pantalla las variables que estás utilizando:
>> who
Your variables are:
ans x y z
>> whos
Name Size Bytes Class Attributes
ans 1x1 8 double
x 1x1 8 double
y 1x1 8 double
z 1x1 8 double
Borra la variable z:
>> clear z
Vuelve a hacer que aparezcan por pantalla las variables que estás utilizando:
>> who
Your variables are:
ans x y
Crea un vector v3 donde sus elementos vayan desde el 2 al 20 y que en total tenga 10 elementos:
>> v3=linspace(2,20,10)
v3 =
2 4 6 8 10 12 14 16 18 20
1 2 3 4
Crea la matriz M= 5 6 7 8 :
9 10
11 12
>> M=[1 2 3 4;5 6 7 8;9 10 11 12]
M=
1 2 3 4
5 6 7 8
9 10 11 12
Calcula el rango de M:
>> rank(M)
ans =
2
Calcula la traza de la matriz M:
>> trace(M)
ans =
19
Crea una matriz llamada M2 que tenga por diagonal el vector v y el resto sean todo ceros:
>> M2=diag(v)
M2 =
1 0 0 0
0 2 0 0
0 0 3 0
0 0 0 4
Calcula el determinante de M:
>> det(M)
ans =
0
Crea un polinomio q que tenga por raíces los elementos del vector v:
>> q=poly(v)
q=
1 -10 35 -50 24
Crea un pequeño programa que escriba por pantalla BUENOS DÍAS y a continuación los 15
primeros números pares (usando un for):
disp('BUENOS DÍAS')
for i=1:15
disp(2*i)
end