Post on 14-Feb-2018
L. Maria Pimentel Herrera
UNI - CEPS
M A T L A B
Contenido Introduccin 5 Entorno de trabajo de matlab 6 Path de matlab 7 Workspace browser y array editor 7 Formatos de salida 8 Guardar variables y estados de una sesin 9 Guardar sesin y copiar salidas 9 Lneas de comentarios 10 Help 10 Medida de tiempos y de esfuerzo de clculo 15 Operaciones con nmeros reales 17 Operaciones aritmticas 17 Constantes matemticas 19 Referencia de funciones 20 Operaciones con Matrices y Vectores 21 Definicin de variables 21 Definicin de matrices desde teclado 21 Tipos de matrices 22 Matrices especiales 25 Operaciones de Matrices 29 Funciones matemticas elementales 33 Operaciones elemento a elemento 35 Generador de vectores 36 Comandos relacionados con tamao de datos 36 Cambiar elementos en una matriz 37 Crear submatrices de una matriz 37 Concatenacin de vectores y matrices 40 Nmeros y matrices asociados a A 43 Operadores relacionales 47 Operadores lgicos 47 Variables lgicas 48 Ejercicios 49 Polinomios 52 Ejercicios 54 lgebra Lineal 55 Matriz de logaritmos y exponenciales 55 Anlisis de matriz 56 Valores propios y valores singulares 57
http://www.mat.ucm.es/~infante/matlab/notas.htm#lgicas#lgicas
Anlisis de datos 59 Operaciones bsicas 59 Estadstica Descriptiva 60 Derivadas e integrales 61 Ejercicio 62 Programacin 63 Ordenes de gestin de archivos 64 Bifurcaciones y bucles 64 Sentencia if 65 Sentencia switch 69 Sentencia for 70 Sentencia break 73 Sentencia continue 73 Aplicacin 73 Funciones 74 Funciones de usuario propias 77 Funcin que devuelve una sola variable 77 Funcin que devuelve mltiples variables 77 Funcin que utiliza otra funcin 78 Grficos 79 Grficos 2D 79 Grficos Estadsticos 88 Grficos 3D 96 Superficies de revolucin 102 Grficos de funciones complejas 102 Grficas en movimiento 102
Introduccin
MatLab significa MATrix LABoratory. Es un programa para hacer computacin numrica. Fue diseado para manipular matrices y ploteo de datos. Ahora incluye funciones para: analizar datos, procesar seales, optimizar funciones. Contiene funciones para los grficos 2-Dy 3-D con su respectiva animacin. Matlab permite leer y escribir archivos .MAT, .TXT, etc. Tiene interfaces con otros lenguajes. Permite la computacin simblica con el Maple. reas de aplicacin
Ingeniera elctrica y mecnica
Informtica
Matemticas
Ingeniera aeroespacial y automotriz
Ingeniera qumica y biomdica
Economa y finanzas Matlab
Trabaja nmeros escalares (reales y complejos), con caracteres y otras estructuras de datos.
Tiene un lenguaje de programacin propio. Permite un rpido prototipeo de aplicaciones cientficas. Pero puede ser ms lento que C/C++ o Fortran. Dispone de cdigo bsico y toolboxes.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
6 | L. Maria Pimentel Herrera
Entorno de trabajo de matlab El entorno de trabajo de MATLAB es muy grfico e i ntuitivo, similar al de ot ras aplicaciones profesionales de Windows.
Las componentes ms importantes del entorno de trabajo de MATLAB son las siguientes:
1. El Escritorio de Matlab (Matlab Desktop), que es la ventana o contenedor de mximo nivel en la que se pueden situar (to dock) las dems componentes.
2. Las componentes individuales, orientadas a tareas concretas, entre las que se puede citar:
a. La ventana de comandos (Command Window),
b. La ventana histrica de comandos (Command History),
c. El espacio de trabajo (Workspace),
d. La plataforma de lanzamiento (Launch Pad),
e. El directorio actual (Current Directory),
f. La ventana de ayuda (Help)
g. El editor de ficheros y depurador de errores (Editor&Debugger),
h. El editor de vectores y matrices (Array Editor).
i. La ventana que permite estudiar cmo se emplea el tiempo de ejecucin (Profiler).
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 7
Path de matlab
establecer el camino de bsqueda (search path) MATLAB puede l lamar a una g ran variedad de f unciones, tanto propias como programadas por los usuarios. Puede incluso haber funciones distintas con el mismo nombre. Interesa saber cules son las reglas que determinan qu funcin o qu fichero *.m es el que se va a ejecutar cuando su nombre aparezca en una lnea de comandos del programa. Esto queda determinado por el camino de bsqueda (search path) que el programa utiliza cuando encuentra el nombre de una funcin.
El search path de MATLAB es una lista de directorios que se puede ver y modificar a partir de la lnea de comandos, o utilizando el cuadro de dilogo Set Path, del men File. El comando path hace que se escriba el search path de MATLAB (el resultado depende de en qu directorio est instalado MATLAB):
>> path
Workspace browser y array editor El espacio de t rabajo de MATLAB (Workspace) es el conjunto de v ariables y de funciones de usuario que en un d eterminado momento estn definidas en la memoria del programa o de la funcin que se est ejecutando. Para obtener informacin sobre el Workspace desde la lnea de co mandos se pueden utilizar los comandos who y whos. El segundo proporciona una informacin ms detallada que el primero.
>>x=9
x=
9
>>y=15
y=
15
>>whos
Name Size Bytes Class x 1x1 8 double array y 1x1 8 double array Grand total is 2 elements using 16 bytes
La ventana Workspace constituye un entorno grfico para ver las variables definidas en el espacio de trabajo. Se activa con el comando View/Workspace. La ventana Workspace cuando se abre desde un determinado programa. Haciendo doble clic por ejemplo sobre la matriz BARS aparece una nueva ventana (o pestaa, si la ventana ya exista) del Array Editor, en l a que se muestran y pueden se r modificados los elementos de dicha matriz.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
8 | L. Maria Pimentel Herrera
Formatos de salida Respecto a l os formatos numricos con que MATLAB muestra los resultados (recurdese que siempre calcula con doble precisin, es decir con unas 16 cifras decimales equivalentes), y son las siguientes:
short coma fija con 4 decimales (defecto)
long coma fija con 15 decimales
hex cifras hexadecimales
bank nmeros con dos cifras decimales
short e notacin cientfica con 4 decimales
short g notacin cientfica o decimal, dependiendo del valor
long e notacin cientfica con 15 decimales
long g notacin cientfica o decimal, dependiendo del valor
rational expresa los nmeros racionales como cocientes de enteros
Estos formatos se pueden cambiar tambin desde la lnea de comandos anteponiendo la palabra format. Por ejemplo, para ver las matrices en formato long habr que ejecutar el comando:
>>format long
>>1/3
ans =
0.33333333333333
>>format % Vuelve al formato estndar que es el de 4 cifras decimales
>> m=17/3;
>> c=9/1974;
Comando Representacin de m Representacin de c
format short 5.6667
format short e 5.6667e+000
format long 5.666666666666667
format long e 5.666666666666667e+000
format hex 4016aaaaaaaaaaab
format bank 5.67
format rat 17/3
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 9
Guardar variables y estados de una sesin
Comandos save y load En muchas ocasiones puede resultar interesante interrumpir el trabajo con MATLAB y poderlo recuperar ms tarde en el mismo punto en el que se dej (con las mismas variables definidas, con los mismos resultados intermedios, etc.). Hay que tener en cuenta que al salir del programa todo el contenido de la memoria se borra automticamente.
Para guardar el estado de una sesin de trabajo existe el comando save. Si se teclea:
>> save
antes de abandonar el programa, se crea en el directorio actual un fichero binario llamado matlab.mat (o matlab) con el estado de la sesin (excepto los grficos, que por ocupar mucha memoria hay que guardar aparte). Dicho estado puede recuperarse la siguiente vez que se arranque el programa con el comando:
>> load
Esta es la forma ms bsica de los comandos save y load. Se pueden guardar tambin matrices y vectores de forma selectiva y en f icheros con nombre especificado por el usuario. Por ejemplo, el comando (sin comas entre los nombres de variables):
>> save filename A x y
guarda las variables A, x e y en un f ichero binario llamado filename.mat (o filename). Para recuperarlas
en otra sesin basta teclear:
>> load filename
Si no se indica ninguna variable, se guardan todas las variables creadas en esa sesin.
Guardar sesin y copiar salidas
Comando diary Los comandos save y load crean ficheros binarios o ASCII con el estado de la sesin. Existe otra forma ms sencilla de almacenar en un fichero un texto que describa lo que el programa va haciendo (la entrada y salida de los comandos utilizados). Esto se hace con el comando diary en la forma siguiente:
>> diary filename.txt
...
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
10 | L. Maria Pimentel Herrera
>> diary off
...
>> diary on
...
El comando diary off suspende la ejecucin de diary y diary on la reanuda. El simple comando diary pasa de on a off y viceversa. Para poder acceder al fichero filename.txt con Notepad o Word es necesario que diary est en off. Si en el comando diary no se incluye el nombre del fichero se utiliza por defecto un fichero llamado diary (sin extensin).
Lneas de comentarios El carcter tanto por ciento (%) indica comienzo de comentario. Cuando aparece en una lnea de comandos, el programa supone que todo lo que va desde ese carcter hasta el fin de la lnea es un comentario. Otra forma de comentar bloques de sentencias (similar a la utilizada en C/C++ con /* y */) es encerrar las lneas que se desea inutilizar entre los caracteres %{ y %}. Los bloques comentados pueden incluirse dentro de otros bloques comentados ms amplios (bloques anidados).
Help! Este comando nos permite solicitar ayuda sobre cualquier comando o funcin que se encuentre instalada en Matlab.
Escribiendo help en la lnea de comando, el programa devuelve un listado de todas las libreras instaladas. Entonces:
>>help
HELP topics:
matlab\general - General purpose commands.
matlab\ops - Operators and special characters.
matlab\lang - Programming language constructs.
matlab\elmat - Elementary matrices and matrix manipulation.
matlab\randfun - Random matrices and random streams.
matlab\elfun - Elementary math functions.
matlab\specfun - Specialized math functions.
matlab\matfun - Matrix functions - numerical linear algebra.
matlab\datafun - Data analysis and Fourier transforms.
matlab\polyfun - Interpolation and polynomials.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 11
matlab\funfun - Function functions and ODE solvers.
matlab\sparfun - Sparse matrices.
matlab\scribe - Annotation and Plot Editing.
matlab\graph2d - Two dimensional graphs.
matlab\graph3d - Three dimensional graphs.
..
clic en matlab\general te mostrar
General purpose commands.
MATLAB Version 7.8 (R2009a) 15-Jan-2009
General information.
syntax - Help on MATLAB command syntax.
demo - Run demonstrations.
ver - MATLAB, Simulink and toolbox version information.
version - MATLAB version information.
verLessThan - Compare version of toolbox to specified version string.
Managing the workspace.
who - List current variables.
whos - List current variables, long form.
clear - Clear variables and functions from memory.
onCleanup - Specify cleanup work to be done on function completion.
pack - Consolidate workspace memory.
load - Load workspace variables from disk.
save - Save workspace variables to disk.
saveas - Save Figure or model to desired output format.
memory - Help for memory limitations.
recycle - Set option to move deleted files to recycle folder.
quit - Quit MATLAB session.
exit - Exit from MATLAB.
Managing commands and functions.
what - List MATLAB-specific files in directory.
type - List M-file.
open - Open files by extension.
which - Locate functions and files.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
12 | L. Maria Pimentel Herrera
pcode - Create pre-parsed pseudo-code file (P-file).
mex - Compile MEX-function.
inmem - List functions in memory.
namelengthmax - Maximum length of MATLAB function or variable name.
Managing the search path.
path - Get/set search path.
addpath - Add directory to search path.
rmpath - Remove directory from search path.
rehash - Refresh function and file system caches.
import - Import packages into the current scope.
finfo - Identify file type against standard file handlers on path.
genpath - Generate recursive toolbox path.
savepath - Save the current MATLAB path in the pathdef.m file.
Managing the java search path.
javaaddpath - Add directories to the dynamic java path.
javaclasspath - Get and set java path.
javarmpath - Remove directory from dynamic java path.
Controlling the command window.
echo - Echo commands in M-files.
more - Control paged output in command window.
diary - Save text of MATLAB session.
format - Set output format.
beep - Produce beep sound.
desktop - Start and query the MATLAB Desktop.
preferences - Bring up MATLAB user settable preferences dialog.
Operating system commands.
cd - Change current working directory.
copyfile - Copy file or directory.
movefile - Move file or directory.
delete - Delete file or graphics object.
pwd - Show (print) current working directory.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 13
dir - List directory.
ls - List directory.
fileattrib - Set or get attributes of files and directories.
isdir - True if argument is a directory.
mkdir - Make new directory.
rmdir - Remove directory.
getenv - Get environment variable.
! - Execute operating system command (see PUNCT).
dos - Execute DOS command and return result.
unix - Execute UNIX command and return result.
system - Execute system command and return result.
perl - Execute Perl command and return the result.
computer - Computer type.
isunix - True for the UNIX version of MATLAB.
ispc - True for the PC (Windows) version of MATLAB.
Debugging.
debug - List debugging commands.
Tools to locate dependent functions of an M-file.
depfun - Locate dependent functions of an M-file or P-file.
depdir - Locate dependent directories of an M-file or P-file.
Loading and calling shared libraries.
calllib - Call a function in an external library.
libpointer - Creates a pointer object for use with external libraries.
libstruct - Creates a structure pointer for use with external libraries.
libisloaded - True if the specified shared library is loaded.
loadlibrary - Load a shared library into MATLAB.
libfunctions - Return information on functions in an external library.
libfunctionsview - View the functions in an external library.
unloadlibrary - Unload a shared library loaded with LOADLIBRARY.
java - Using Java from within MATLAB.
usejava - True if the specified Java feature is supported in MATLAB.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
14 | L. Maria Pimentel Herrera
Controlling multithreading setting.
maxNumCompThreads - Controls the maximum number of computational threads.
See also lang, datatypes, iofun, graphics, ops, strfun, timefun,
matfun, demos, graphics, datafun, uitools, doc, punct, arith.
Para pedir mas detalles sobre las funciones que pertenecen a una l ibrera dada, ingresamos help seguido del nombre de la librera. Por ejemplo:
>>help stats
La librera stats agrupa diferentes rutinas tiles en probabilidad y estadstica.
Al final de l a ayuda nos remite a al gunos temas relacionados para que podam os continuar la bsqueda, si es que no terminamos de encontrar lo que buscbamos.
Si quisiramos ver con mas detalle de algn algn tem de la lista, basta con escribir help .
>>help fliplr
Ante cualquier duda sobre el help escribe:
>>help help
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 15
Hay otro modo de ayuda, un poco mas cmodo, que se puede acceder desde el men desplegable Help. El contenido es el mismo que el de la lnea de comandos, solo que disponemos de un pequeo navegador.
>>demos
Medida de tiempos y de esfuerzo de clculo MATLAB dispone de funciones que permiten calcular el tiempo empleado en l as operaciones matemticas realizadas. Algunas de estas funciones son las siguientes:
cputime devuelve el tiempo de C PU (con precisin de centsimas de segundo) desdeque el programa arranc. Llamando antes y despus de realizar una oper acin y restando los valores devueltos, se puede saber el tiempo de CPU empleado en esa operacin. Este tiempo sigue corriendo aunque MATLAB est inactivo.
etime(t2, t1) tiempo transcurrido entre los vectores t1 y t2 (atencin al orden!), obtenidos como respuesta al comando clock.
tic ops toc imprime el tiempo en segundos requerido por ops. El comando tic pone el reloj a cero y toc obtiene el tiempo transcurrido.
A modo de ejemplo, el siguiente cdigo mide de varias formas el tiempo necesario para resolver un sistema de 1000 ecuaciones con 1000 incgnitas. Tngase en cuenta que los tiempos pequeos (del orden de las dcimas o centsimas de segundo), no se pueden medir con gran precisin.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
16 | L. Maria Pimentel Herrera
>> n=1000; A=rand(n); b=rand(n,1); x=zeros(n,1);
>> tiempoIni=clock; x=A\b; tiempo=etime(clock, tiempoIni)
>> time=cputime; x=A\b; time=cputime-time
>> tic; x=A\b; toc
donde se han puest o varias sentencias en la misma lnea para que se ejecuten todas sin tiempos muertos al pulsar intro. Esto es especialmente importante en la lnea de comandos en la que se quiere medir los tiempos. Todas las sentencias de clculos matriciales van seguidas de punto y coma (;) con objeto de ev itar la impresin de resultados. Conviene ejecutar dos o tres veces cada sentencia para obtener tiempos ptimos, ya que la primera vez que se ejecutan se emplea un cierto tiempo en cargar las funciones a memoria.
Nota: Un punto y coma al final de una sentencia hace que no se vea el resultado de la operacin. Esto es muy importante, porque Matlab tiene la costumbre de ir mostrando todos los resultados obtenidos. El punto y coma nos permite filtrar los resultados parciales y en todo caso solo mostrar aquellos que sean de nuestro inters.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 17
Operaciones con nmeros reales Los clculos que no se asignan a una variable en concreto se asignan a la variable de respuesta por defecto que es ans (del ingls, answer):
>>2+3
ans =
5
Sin embargo, si el clculo se asigna a una variable, el resultado queda guardado en ella:
>>x=2+3
x =
5
Para conocer el valor de una variable, basta teclear su nombre:
>>x
x =
5
Si se aade un punto y coma (;) al final de la instruccin, la mquina no muestra la respuesta...
>>y=5*4;
... pero no por ello deja de realizarse el clculo.
>>y
y =
20
Operaciones aritmticas Las operaciones se evalan por orden de prioridad: primero las potencias, despus las multiplicaciones y divisiones y, finalmente, las sumas y restas. Las operaciones de igual prioridad se evalan de izquierda a derecha.
Prioridad Operador Significado Ejemplo Resultado
1 ^ Exponente 78^2
2 / Divisin derecha 29/23
2 \ Divisin izquierda 29\23
3 * Multiplicacin 17*28
4 + Adicin 15.3+19.5
4 - Sustraccin 15.3-19.5
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
18 | L. Maria Pimentel Herrera
>>2/4*3
ans =
1.5000
>>2/(4*3)
ans =
0.1667
Se pueden utilizar las funciones matemticas habituales. As, por ejemplo, la funcin coseno,
>>cos(pi) % pi es una variable con valor predeterminado 3.14159...
ans =
-1
o la funcin exponencial
>>exp(1) % Funcin exponencial evaluada en 1, es decir, el nmero e
ans =
2.7183
Adems de la variable pi , MATLAB tiene otras variables con valor predeterminado; ste se pierde si se les asigna otro valor distinto. Por ejemplo:
>>eps % psilon de la mquina. Obsrvese que MATLAB trabaja en doble precisin
ans =
2.2204e-016
pero...
>>eps=7
eps =
7
Otro ejemplo de funcin matemtica: la raz cuadrada; como puede verse, trabajar con complejos no da ni ngn tipo de pr oblema. La uni dad imaginaria se representa en MATLAB como i o j, variables con dicho valor como predeterminado:
>>sqrt(-4)
ans =
0+ 2.0000i
El usuario puede controlar el nmero de deci males con que aparece en pant alla el valor de las variables, sin olvidar que ello no est relacionado con la precisin con la que se hacen los clculos, sino con el aspecto con que stos se muestran:
>>1/3
ans =
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 19
0.3333
Para conocer las variables que se han usado hasta el momento:
>>who
Your variables are: ans eps x y
o, si se quiere ms informacin (obsrvese que todas las variables son arrays):
>>whos
Name Size Bytes Class ans 1x1 8 double array eps 1x1 8 double array x 1x1 8 double array y 1x1 8 double array Grand total is 4 elements using 32 bytes
Para deshacerse de una variable
>>clear y
>>who
Your variables are: ans eps x
Constantes matemticas
Variable Valor
eps Nmero ms pequeo tal que, cuando se le suma 1, crea un nmero en coma flotante en el computador mayor que 1.
pi Relacin entre la circunferencia del crculo a su dimetro
intmax Mayor valor de tipo entero especificado.
intmin Menor valor de tipo entero especificado.
inf Infinito.
NaN Magnitud no numrica.
i y j i=j=-1, unidad imaginaria.
realmin El nmero real positivo ms pequeo que es utilizable.
relamas El nmero real positivo ms grande que es utilizable.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
20 | L. Maria Pimentel Herrera
Referencia de funciones
Escritorios y el Desarrollo para el Medio Ambiente
De inicio, ventana de comandos, ayuda, edicin y depuracin, tuning, otras funciones generales
Importacin y exportacin de datos
General y de bajo nivel / S de archivos, adems de formatos de archivo especficos, como de audio, hoja de clculo, HDF, imgenes
Matemticas Matrices y matrices, lgebra lineal, otras reas de las matemticas
Anlisis de datos Datos bsicos de las operaciones, estadstica descriptiva, covarianza y correlacin, filtrado y convolucin, derivadas numricas e integrales, transformadas de Fourier, anlisis de series temporales
De Programacin y tipos de datos
La evaluacin de funciones de expresin, el control del programa, la funcin de asas, programacin orientada a objetos, el manejo de errores, los operadores, tipos de datos, fechas y horas, los temporizadores
Programacin orientada a objetos
Funciones para trabajar con clases y objetos
Grficos Parcelas de la Lnea, anotar los grficos, grficos especializados, las imgenes, la impresin, Handle Graphics
La visualizacin en 3-D De superficie y las parcelas de malla, control de vista, la iluminacin y la transparencia, la visualizacin de volmenes
GUI para el Desarrollo GUA, la programacin de interfaces grficas de usuario
Interfaces externos Interfaces con librariess compartida, Java,. NET, COM y ActiveX, servicios Web y los dispositivos de puerto serie, y C y Fortran rutinas
http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg01.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhiZeFfDAJqK2z9ZjDJO5-IdV1t2EA#f16-11063http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg01.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhiZeFfDAJqK2z9ZjDJO5-IdV1t2EA#f16-11063http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg02.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhhF1tQdmzr_dxa5yS6mhdkQjIImHQ#f16-5702http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg02.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhhF1tQdmzr_dxa5yS6mhdkQjIImHQ#f16-5702http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg03.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhhdbbiewaiGEWE8vJYgxjME5XEtsQ#f16-5872http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg04.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhhJQNZjH84mKYJnqrOkwTjC1IjhNQ#f16-48518http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg05.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhi1fRU885BXhkST02syBX9novln9A#f16-42340http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg05.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhi1fRU885BXhkST02syBX9novln9A#f16-42340http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg06.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhgKgPnVGEtd0tQ0omZDAlp9G3Frqg#brk7uzkhttp://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg06.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhgKgPnVGEtd0tQ0omZDAlp9G3Frqg#brk7uzkhttp://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg07.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhijKnxjxT0DeZK1ACc4oaWwYIAgqw#f16-8602http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg08.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhiyXMG9BK1fWhe-eqDw1EBIa5PlfA#f16-8867http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg09.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhhjX37hAWj8LDbmYB0y2kVV9BVvnA#f16-40727http://translate.googleusercontent.com/translate_c?hl=es&sl=en&u=http://www.mathworks.com/access/helpdesk/help/techdoc/ref/f16-6011seg10.html&prev=/search%3Fq%3Dmatlab%2Bwhat%2Bdir%2Bls%26hl%3Des&rurl=translate.google.com.pe&usg=ALkJrhjwP-B11Tz4_NmHkuFAcv2IJPKrBA#f16-35614
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 21
Operaciones con Matrices y Vectores Definicin de variables No es necesario definir variables.
Las variables mejoran la legibilidad del procedimiento y facilitan las correcciones y modificaciones.
En matlab las variables deben seguir las siguientes reglas:
No pueden comenzar con un nmero.
Puede tener nmeros en la estructura del nombre de la variable.
Las maysculas y minsculas se diferencian en el nombre de las variables.
Los nombres de las variables no pueden tener los siguientes smbolos: +, -, *, /, ., ,, ;, ^, ~, &, |, \.
Definicin de matrices desde teclado Para introducir una matriz:
Se separan los nmeros con espacio o comas.
Se separan las columnas con punto y coma.
Se agrupa toda la matriz entre corchetes.
Vectores fila
>> x=[1,2,3,5,7,11,13]; [ . . . . . . . ]
>> x2=[1 2 3 5 7 11 13]; [ . . . . . . . ]
>>a(5)=7; [ . . . . . ]
Vectores columna
>> X=[1;2;3;5;7;11;13]
X =
.
.
.
.
.
.
.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
22 | L. Maria Pimentel Herrera
.
Matriz
>> A = [16 3 2 13; 5 10 11 8; 9 6 7 12; 4 15 14 1]
A =
................
................
................
................
>> M = [1 2 3
4 5 6
7 8 9]
M=
............
............
............
Tipos de matrices Matriz diagonal
diag(v) genera una matriz diagonal con el vector v como diagonal.
diag(A) crea un vector con los elementos de la diagonal principal.
>> T=[17 25 29];
>> S=diag(T)
S=
............
............
............
>> diag(A)
ans =
....
....
....
....
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 23
Matriz simtrica diagonal constante
toeplitz(v) define una m atriz simtrica de diagonal constante con v como primera fila y primera columna.
>> v=[1 2 3]; toeplitz(v)
ans=
............
............
............
Matriz diagonal constante
toeplitz(w, v) define una matriz no simtrica de diagonal constante con w como primera columna y v como primera fila.
>> w=[1 2 3 4]; v=[1 5 7 8]; toeplitz(w,v)
ans=
................
................
................
................
Matriz de unos
ones(n) genera una matriz de n n con todos los valores iguales a uno.
ones(n,m) genera una matriz de n m con todos los valores iguales a uno.
>> unos=ones(3,4)
unos=
................................................
Matriz nula
zeros(n) genera una matriz de n n con todos los valores iguales a cero.
Zeros(m,n) crea una matriz mxn, todos ceros.
>> ceros=zeros(2,5)
ceros=
....................
....................
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
24 | L. Maria Pimentel Herrera
Matriz identidad
eye(n) genera una matriz identidad de n n.
eye(m,n) crea la matriz identidad de orden mxn con unos en la diagonal y ceros en los dems elementos.
>> ident=eyes(3)
ident=
............
............
............
Matriz aleatoria uniforme
rand(n) genera una matriz de n n con elementos de valor aleatorio entre 0 y 1 (distribucin uniforme).
>> aleatorio=rand(2,3);
aleatorio=
................................................................................................
Matriz aleatoria normal
randn(n) genera una matriz de n n cu yos elementos siguen una distribucin normal (media 0 y varianza 1).
>> normal=randn(2);
normal=
............................................................
>> n1=randn(3,2)
n1=
..............................
..............................
..............................
Nota
ones(m, n), zeros(m, n), rand(m, n) generan matrices de m n.
ones(size(A)), zeros(size(A)), eye(size(A)) generan matrices de la misma forma que A.
Muchos problemas lineales implican el manejo de matrices de gran tamao, en las que se da la circunstancia que la mayor parte de los elementos son nulos. Esto implica,
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 25
eventualmente, el almacenamiento de num erosos ceros y operaciones intiles. Para resolver este problema, MATLAB dispone de una f uncin (sparse) que permite prescindir de l as cantidades nulas, almacenando la posicin (fila y columna) de l os elementos no nulos. Un ejemplo elemental se construye con la matriz unidad:
>> clear
>> I=eye(1000);
>> whos
Name Size Bytes Class
I 1000x1000 8000000 double array
Grand total is 1000000 elements using 8000000 bytes
Observamos que el almacenamiento de I consume aproximadamente 8 Mb.
Determinamos una expresin sparse de I:
>> I=eye(1000);
>> S=sparse(I);
>> whos
Name Size Bytes Class
I 1000x1000 8000000 double array
S 1000x1000 16004 sparse array
Grand total is 1001000 elements using 8016004 bytes
Matrices especiales Matriz Transpuesta
A El caracter ' (apstrofe) denota la transpuesta de l a matriz. Transponer significa intercambiar filas por columnas. Si tenemos la matriz A y llamamos B = A', B es la transpuesta de la matriz A.
>> A=[5 9 11; 6 8 10]; C=A'
C=
Matriz Simtrica
Es una matriz cuadrada que cumple la condicin: A=AT.
>> M=[2 3 17; 3 -6 1; 17 1 7], Mt=M'
M=
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
26 | L. Maria Pimentel Herrera
............
............
............
Mt=
............
............
............
Matriz Antisimtrica
Es una matriz cuadrada que cumple la condicin: A=-AT, siendo ceros los elementos de la diagonal principal y sus dems elementos son de la forma aij=aji.
>> B=[0 1 -4; -1 0 -3; 4 3 0], Ba=-B'
B=
............
............
............
Ba=
............
............
............
Matriz de indices no nulos
Devuelve los ndices de los elementos de la matriz A, que no son ceros, en forma vertical.
>> B=[0 1 -4; -1 0 -3; 4 3 0]; ind=find(B)
ind=
2.00 3.00 4.00 6.00 7.00 8.00
>> ind2=find(B>0)
Ind2 =
3.00 4.00 6.00
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 27
>> [f,c,v]=find(B>0); fila=f',columna=c',verdad=v'
fila =
3.00 1.00 3.00
columna =
1.00 2.00 2.00
verdad =
1 1 1
Matriz Idempotente
Es una matriz cuadrada que cumple la condicin: A=A2.
>> V=[-1 2 4; 1 -2 -4; -1 2 4], Vi=V^2
V=
............
............
............
Vi=
............
............
............
Matriz Peridica
Es una matriz cuadrada que cumple la condicin: A=Ak+1 entonces A es peridica y kZ+. El periodo es igual a k.
>> P=[1/3 2; 1/9 2/3], Pp=P^2
P=
..............................
..............................
Pp=
..............................
..............................
Matriz Nilpotente
Es una matriz cuadrada que cumple la condicin: Ap=0, donde pZ. A es nilpotente para p.
>> N=[1 1 3; 5 2 6; -2 -1 -3], Ni=N^3
N=
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
28 | L. Maria Pimentel Herrera
............
............
............
Ni=
............
............
............
Matriz Hilbert
hilb(N) es una matriz de N por N con elementos 1/(i+j-1).
>> format rat ; hilb(3)
ans =
1 1/2 1/3
1/2 1/3 1/4
1/3 1/4 1/5
Matriz Hermitiana
Es una matriz cuadrada y compleja, que es igual a la transpuesta de su conjugada. Los elementos de su diagonal principal son nmeros reales..
>> H=[1 3+i i; 3-i 3 1-i; -i 1+i 2], Ht=conj(H), He=transpose(Ht)
H=
............
............
............
Ht=
............
............
............
He=
............
............
............
Matriz Inversa
inv(A) A-1 si A es cuadrada e invertible, se cumple A-1*A=A*A-1=I.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 29
>> A=[1 1 1; 1 2 3; 1 3 4]; format rat; S1=inv(A)
S1=
............
............
............
Matriz Seudoinversa
pinv(A) pseudoinversa de A , si X = pinv(A) produce una matriz X condimensiones de A, se cumple A*X*A = A, X*A*X = X y A*X y X*A son matrices hermitianas.
>> K=magic(4); K1=det(K), K2=pinv(K)
K1=
K2=
............................................................
............................................................
............................................................
............................................................
Matriz Ortogonal
Una matriz cuadrada es ortogonal si se cumple A-1= AT, es equivalente a A*AT =I .
>> a=sqrt(2); M=[1/a 0 -1/a; 1/a 0 1/a; 0 1 0], L=M*M'
M=
............
............
............
L=
............
............
............
Operaciones de Matrices Sumando y Restando Matrices
Las operaciones suma (+) y resta (-) son definidas para las matrices siempre y cuando stas tengan la misma dimensin. Es decir, si A y B son matrices 3 x 3, entonces A + B se puede calcular.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
30 | L. Maria Pimentel Herrera
Las operaciones suma y resta tambin estn definidas si uno de los operadores es un escalar, es decir, una matriz 1 x 1.
>> A= [4 2 0 1; -2 3 6 5; 2 1 8 1];
>> B= [3 17 8 5; -7 12 15 10; 23 19 0 -25];
>> A+B
ans=
................
................
................
................
>> f=[2 4 6]; c=[1; 2; 3];
>>S1= f+c'
S1=
................
................
................
................
>>S2=f+25
S2=
................
................
................
................
>>S3=c-17
S3=
................
................
................
................
Multiplicando Matrices
La operacin de m ultiplicacin de m atrices est definida siempre que el nmero de columnas de la primera matriz sea igual al nmero de filas de la segunda matriz.
A * B da la matriz resultante del producto AB
(si dicha operacin es posible).
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 31
A. * B da el producto elemento por elemento
(si size(A) = size(B), es decir, si tienen el mismo tamao)
>> A*B
ans=
................
................
................
................
Producto de una matriz por un vector
El producto de una matriz y un vector es un caso especial del producto matriz-matriz y naturalmente, un escalar como pi, puede multiplicar, ser multiplicado por, cualquier matriz.
>> 7*A
ans=
................
................
................
................
>> X*B
ans=
....
....
....
....
Dividiendo Matrices
En divisin de matrices, si A es una matriz cuadrada no-singular, entonces A\B y B/A corresponden a la multiplicacin izquierda y derecha de B por el inverso de A, esto es, inv(A) * B y B * inv(A) respectivamente. El resultado es obtenido directamente sin la computacin del inverso.
>> A=[0 1;1 0]; B=[1 2;3 4]; X=A\B
X =
3 4
1 2
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
32 | L. Maria Pimentel Herrera
>> Y=B/A
Y =
2 1
4 3
X = A\B es una solucin a A * X = B,
es igual a inv(A) * B si existe inv(A): la barra inversa es la divisin por la izquierda.
X = B/A es una solucin a X * A = B
A\B es definido cuando B tiene la misma cantidad de filas que A. Si A es cuadrada, el mtodo usado es la Eliminacin Gaussiana. El resultado es una matriz X con las mismas dimensiones que B.
Si A no es cuadrada, se factoriza utilizando la ortogonalizacin de Householder con pivoteo de columnas.
Los factores son usados para resolver sistemas de ecuaciones sub-determinados y sobre -determinados. El resultado es una matriz X m-por-n donde m es el nmero de columnas de A y n es el nmero de columnas de B. Cada columna de X tiene, al menos, k componentes diferentes de cero, donde k es el rango efectivo de A.
B/A esta definido en trminos de A\B por B/A = (A' \B') '.
Ejemplo: resolver el sistema:
2 a + 3 b + c = 6
4 a + b + 2 c = 7
6 a + b + 7 c = 4
>>A=[2 3 1; 4 1 2; 6 1 7];
>>B=[6; 7; 4]
>>X=A\B
x=
.....
.....
.....
Exponentes con Matrices
La expresin A^n eleva A a la n-sima potencia y est definido si A es una matriz cuadrada y n un escalar.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 33
>>p1=A^2
p1 =
............
............
............
Funciones Matriciales Trascendentales y Elementales
MATLAB considera expresiones como exp(A) y sqrt(A) como operaciones de arreglos, definidas en los elementos individuales de A. Tambin puede ca lcular funciones trascendentales de matrices, como la matriz exponencial y la matriz logartmica. Estas operaciones especiales estn definidas solamente para matrices cuadradas.
>>sqrt(A)
ans =
............
............
............
>>exp(A)
ans =
............
............
............
Funciones matemticas elementales
Trigonomtricas
sin(x) Funcin seno
cos(x) Funcin coseno
tan(x) Funcin tangente
asin(x) Inversa del seno
acos(x) Inversa del coseno
atan(x) Inversa de la tangente
atan2(x,y) Inversa de la tangente de los cuatro cuadrantes
sinh(x) Funcin seno hiperblico
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
34 | L. Maria Pimentel Herrera
cosh(x) Funcin coseno hiperblico
tanh(x) Funcin tangente hiperblica
asin(x) Inversa del seno hiperblico
acosh(x) Inversa del coseno hiperblico
atanh(x) Inversa de la tangente hiperblica
Tambin hay funciones para ngulos en sexagesimales.
Exponencial
exp(x) Exponencial ex.
log(x) Logaritmo natural.
log10(x) Logaritmo decimal.
log2(x) Logaritmo en base 2.
sqrt(x) Raz cuadrada.
Complejo
abs(x) Valor absoluto o magnitud de un nmero complejo.
real(x) Parte real de un nmero complejo.
imag(x) Parte imaginaria de un nmero complejo.
angle(x) ngulo de un nmero complejo.
conj(x) Conjugado complejo
Redondeo y resto
ceil(x) Redondeo hacia ms infinito.
fix(x) Redondeo hacia cero.
floor(x) Redondeo hacia menos infinito.
round(x) Redondea hacia el entero ms prximo.
rem(x,y) Resto despus de la divisin.
mod(x,y) Mdulo despus de la divisin.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 35
Operaciones elemento a elemento Matlab define algunas operaciones que sern realizadas elemento a elemento.
Ya vimos anteriormente que la operacin A*B realizaba el producto matricial.
Pero que pasa si queremos que cada elemento de A quede multiplicado por cada elemento de B (suponiendo que tienen las mismas dimensiones)?
Existe otro operador para tales fines:
Definimos las matrices a y b:
>>A = [1 2 3 ; 4 5 6; 7 8 9]; B= [9 8 7 ; 6 5 4; 3 2 1];
>>C1=A.*B
C1=
............
............
............
Como se ve, si anteponemos un punto al operador, la operacin se realiza elemento a elemento. La divisin (/) y la potencia (^) tambin permiten esta utilizacin.
>>C2=A./B
C2=
............
............
............
>>C3=A.\B
C3=
............
............
............
>>C4=B.^A
C4=
............
............
............
Otras operaciones pueden ser consultadas en la librera ops.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
36 | L. Maria Pimentel Herrera
Generador de vectores Para generar un v ector cuyos elementos sean nmeros crecientes o
decrecientes en un intervalo regular existe el operador : (dos puntos).
El uso ms sencillo de este operador sera:
:
>>d = 1:4
d =
1 2 3 4
Especificar un intervalo determinado, escribimos:
::
>>e = 1:0.5:4
d =
1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000
Si el rango es decreciente, slo escribimos el intervalo negativo.
>>f = 4:-0.5:1
f =
4.0000 3.5000 3.0000 2.5000 2.0000 1.5000 1.0000
Genera un vector con n valores entre x1 y x2 igualmente espaciados.
linspace(x1,x2,n)
>>g = linspace(0,20,7)
g=
0 10/3 20/3 10 40/3 50/3 20
>>h = linspace(0,pi,5)
h=
0 0.7854 1.5708 2.3562 3.1416
Comandos relacionados con tamao de datos
length() se aplica solo a vectores. Devuelve el largo del vector (es igual para filas y columnas).
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 37
>>length(g)
ans=
......
size() se aplica tanto a vectores como matrices. Devuelve un vector de dos elementos: cantidad de filas y cantidad de columnas.
>>size(A)
ans=
...... ......
Cambiar elementos en una matriz A dada
A(3, 2) = 7 coloca un 7 en el elemento (3, 2).
A(3,:) = v sustituye los valores de la tercera fila por los de v= [2; 4; 5].
A(:, 2) = w sustituye los valores de la segunda columna por los de w= 2:4.
El smbolo de los dos puntos : significa todo (todas las columnas o todas las filas).
A([2 3],:) = A([3 2],:) intercambia las filas 2 y 3 de A.
Crear submatrices de una matriz A de m n Matlab utiliza los parntesis para acceder a elementos de la matriz. Los subndices empiezan en 1, por lo tanto el primer elemento es a(1,1) Nota: Matlab no admite el cero como ndice de vectores ni matrices !!!!
A(i, j) muestra el elemento (i, j) de la matriz A (escalar = matriz de 1 1).
A(i, :) muestra la fila i-sima de A (como vector de fila).
A(:, j) muestra la columna j-sima de A (como vector de columna).
A(2: 4,3: 7) muestra las filas de la 2 a la 4 y las columnas de la 3 a la 7
(en forma de matriz de 3 5).
A([2 4],:) muestra las filas 2 y 4 y todas las columnas
(en forma de matriz de 2 n).
A(:) muestra una sola columna larga formada a partir de las columnas de A (matriz de mn 1).
triu(A) coloca ceros en todos los elementos por debajo de la diagonal (triangular superior).
tril(A) coloca ceros en todos lo elementos por encima de la diagonal (triangular inferior).
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
38 | L. Maria Pimentel Herrera
>> a=rand(4,5);
>> a(3,5)=56.8
Se pueden utilizar vectores para definir ndices
>> a(2:3,1:4)=zeros(2,4);
o bien: a(2:3,1:4)=0;
Se pueden utilizar vectores para definir ndices
>> a([2,3],[2,4])=ones(2,2);
o bien: a([2,3],[2,4])=1;
El operador ':' se utiliza para indicar "todos los elementos"
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 39
Triangular inferior
>>tril(M)
ans =
............
............
............
>>tril(M,-1)
ans =
............
............
............
Triangular superior
>>triu(M)
ans =
............
............
............
>>triu(M,1)
ans =
............
............
............
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
40 | L. Maria Pimentel Herrera
Concatenacin de vectores y matrices Supongamos que queremos formar una matriz con diferentes vectores y/o matrices de dimensiones compatibles, o queremos unir dos vectores para formar uno mas largo.
>>v = [1 2 3]; w=[4 5 6];
Podemos concatenarlos de varias maneras...
Si los unimos uno a continuacin del otro:
>>y = [v w]
y =
1 2 3 4 5 6
Notar que al separar v y w por un espacio estamos diciendo que ambos pertenecen a la misma fila, por lo que es entendible el resultado obtenido.
Si por el contrario, interponemos un punto y coma entre ambos:
>>y2 = [v ; w]
Y2 =
1 2 3
4 5 6
Matlab coloca cada vector en una fila diferente.
Funcin cat
cat(d,M,N) concatena matrices con una dimensin especificada.
>>M=[2 3; 1 9]; N=[8 4; 5 7]; MN=cat(1,M,N)
MN =
2 3
1 9
8 4
5 7
>>NM= cat(2,B,C,B)
NM =
35 1 14 16 35 1
3 32 18 11 3 32
>> A=magic(6)
A =
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 41
35 1 6 26 19 24
3 32 7 21 23 25
31 9 2 22 27 20
8 28 33 17 10 15
30 5 34 12 14 16
4 36 29 13 18 11
>> B=A(1:2,1:2); C=A(5:6,5:6); cat(1,B,C)
ans =
35 1
3 32
14 16
18 11
>> cat(2,B,C)
ans =
35 1 14 16
3 32 18 11
Funcin horzcat
concatena matrices horizontalmente.
>> A = magic(5); A(4:5,:) = []
A =
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
>> B = magic(3)*17
B =
136 17 102
51 85 119
68 153 34
>> C = horzcat(A, B)
C =
17 24 1 8 15 136 17 102
23 5 7 14 16 51 85 119
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
42 | L. Maria Pimentel Herrera
4 6 13 20 22 68 153 34
Funcin vertcat
concatena matrices verticalmente.
>> A = magic(5); A(:, 4:5) = []
A =
17 24 1
23 5 7
4 6 13
10 12 19
11 18 25
>> B = magic(3)*7
B =
56 7 42
21 35 49
28 63 14
>> C = vertcat(A,B)
C =
17 24 1
23 5 7
4 6 13
10 12 19
11 18 25
56 7 42
21 35 49
28 63 14
Funcin repmat
crea una nuev a matriz de co pias, m veces verticalmente y n v eces horizontalmente.
repmat(A,m,n)
>> K = repmat([1 2; 3 4],2,3)
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 43
K =
1 2 1 2 1 2
3 4 3 4 3 4
1 2 1 2 1 2
3 4 3 4 3 4
Funcin blkdiag
crea una nueva matriz de bloque diagonal a partir de matrices existentes, los dems elementos de la matriz son ceros.
>> A=7;B=[9 8; 2 5];C=pascal(3);D=[1 7 12]; blkdiag(A,B,C,D)
ans =
7 0 0 0 0 0 0 0 0
0 9 8 0 0 0 0 0 0
0 2 5 0 0 0 0 0 0
0 0 0 1 1 1 0 0 0
0 0 0 1 2 3 0 0 0
0 0 0 1 3 6 0 0 0
0 0 0 0 0 0 1 7 12
Nmeros y matrices asociados a A
Determinante
det(A) es el determinante (si A es una matriz cuadrada).
>> A=[3 2 1; 4 3 0; 7 9 15]; det(A)
ans =
30
>> det(magic(4))
ans =
0
Rango
rank(A) es el rango (nmero de pivotes = dimensin del espacio de filas y del espacio de columnas).
>> rank(A)
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
44 | L. Maria Pimentel Herrera
ans =
3
>> rank(magic(4))
ans =
3
Orden
size(A) es el par de nmeros [m n].
>> size(A)
ans =
3 3
Traza
trace(A) es la traza = suma de l os elementos de la diagonal = suma de autovalores.
>> trace(A)
ans =
21
>> trace(magic(4))
ans =
34
Producto Escalar
El producto interior (producto escalar pr oducto punto) se consigue de l a siguiente manera: x' * y asumiendo que x y y son vectores columnas. Note que y' * x produce el mismo resultado.
>> X=[1, 2, 3, 5];
>> Y=[4, 7, -8, 3];
>> X*Y
ans=
.
dot(u,v) calcula el producto escalar de dos vectores u y v.
>> u=[-7 17 5]; v=[-2 3 4]; dot(u,v)
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 45
ans =
85
Producto Vectorial
crossl(u,v) calcula el producto vectorial de dos vectores u y v en R3.
>> u=[1 7 5]; v=[2 3 8]; cross(u,v), cross(v,u)
ans =
41 2 -11
ans =
-41 -2 11
Norma
norm(A,p) calcula la norma de A, donde p puede ser 1, 2 o inf.
Si p=1, calcula la suma de l os valores absolutos de todos los elementos de A por columna su equivalentes es max(sum(abs(A))).
Si p=2, norma Euclidiana si es un vector, opcin por defecto, su equivalente es max(svd(A)).
Si p=inf, calcula el mximo valor absoluto de sus elementos, su equivalente es max(sum(abs(A'))).
>> v=[3 4 5]; norm(v,1)
ans =
12
>> norm(v,2), norm(v)
ans =
7.0711
ans =
7.0711
>> norm(v,inf)
ans =
5
>> norm(v,-inf)
ans =
3
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
46 | L. Maria Pimentel Herrera
>> N=[8 5; 1 7]; norm(N)
ans =
10.8035
Base ortogonal para el espacio nulo
null(A) es una matriz cuyas columnas n - r forman una base ortogonal para el espacio nulo de A.
>> A = [1 2 3
1 2 3
1 2 3];
Z = null(A);
A*Z
ans =
1.0e-015 *
0.2220 0.2220
0.2220 0.2220
0.2220 0.2220
>> Z'*Z
ans =
1.0000 -0.0000
-0.0000 1.0000
Base ortogonal para el espacio de columnas
B=orth(A) es una matriz cuyas columnas r forman una base ortogonal para el espacio de columnas de A, se cumple B*B=eye(rank(A)).
>> A=[3 2 1; 4 3 0; 7 9 15]; B=orth(A)
B =
-0.1548 -0.4852 -0.8606
-0.1622 -0.8468 0.5066
-0.9745 0.2180 0.0524
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 47
Operadores relacionales
< menor que
mayor que
>= mayor igual a
== igual a
~= distinto a
Operadores lgicos
& Conjuncin, devuelve el valor de 1 en ca da posicin, donde l os elementos de ambas matrices no son ceros, 0 para los otros casos.
| Disyuncin, devuelve el valor de 1 en cada posicin, donde los elementos de por lo menos de una m atriz no e s nula, 0 en o tros casos.
~ Negacin, complementa cada elemento de la matriz de entrada.
xor Disyuncin exclusiva, devuelve el valor de 1, en ca da posicin donde los elementos de ambas matrices son diferentes, 0 para los otros casos.
Inputs and or not xor
A B A & B A | B ~A xor(A,B)
0 0 0 0 1 0
0 1 0 1 1 1
1 0 0 1 0 1
1 1 1 1 0 0
>> A=3*ones(2), B=eye(2)*3
A =
3 3
3 3
B =
3 0
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
48 | L. Maria Pimentel Herrera
0 3
>> A&B
ans =
1 0
0 1
>> A|B
ans =
1 1
1 1
>> ~B
ans =
0 1
1 0
>> xor(A,B)
ans =
0 1
1 0
Variables lgicas
Tambin existen variables lgicas que toman los valores 0 (falso) o 1 (verdadero)
>>v=[-7 3 17];
>>abs(v)>=2 % Vector lgico cuyas coordenadas valen 1 si la % coordenada correspondiente de v es >= 2 y 0 si no lo es
ans =
0 1 1
>>vector=v(abs(v)>=2) % Vector formado por la coordenadas de v que % verifican la desigualdad
vector =
2 3
>>v2=[23 3 12]
v2 =
23 3 12
http://www.mat.ucm.es/~infante/matlab/notas.htm#lgicas#lgicas
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 49
>>logica=v==v2 % Asignacin de un valor lgico (el doble signo igual
% es el igual lgico)
logica =
0 1 0
>>logic2=v~=v2 % Distinto (~ es el operador de negacin)
logic2 =
1 0 1
Ejercicios 1. Crear los siguientes vectores:
a. X=[3 7 e12]
b. Y=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 ]
c. K=
1391642567
2. Crear un vector Z de cuatro nmeros complejos.
3. Listar el tercer elemento del vector.
4. Listar los 5 primeros elementos del vector Y.
5. Listar los 5 ltimos elementos del vector Y.
6. Listar los elementos de posiciones impares del vector Y.
7. Listar los elementos de posiciones 2, 4, 5 y 7 del vector Y.
8. Crear los vectores M=[5 4 3 2 1] y C=[17 9 8 25 12].
9. Fusionar los vectores M y C en un vector J.
10. Obtener la transpuesta del vector K.
11. Obtener la transpuesta del vector Z.
12. Crear las siguientes matrices:
a. R=
87651234
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
50 | L. Maria Pimentel Herrera
b. S=
44332211
13. Sumar las matrices R y S.
14. Multiplicar las matrices R y S.
15. Multiplicar R con la transpuesta de S.
16. Multiplicar R y S componente a componente.
17. Eleve 3 a cada elemento de R.
18. Obtener la inversa de cada elemento de R.
19. Hallar la matriz inversa de K.
20. E = eye(4); E(2, 1) = -3 crea una matriz de eliminacin elemental de 4 4. E*A
resta 3 veces la fila 1 de la fila 2 de A.
21. B = [A b] crea una matriz aumentada con b como columna adicional.
22. E = eye(3); P = E([2 1 3],:) genera una matriz de permutacin.
23. Ntese que triu(A) + tril(A) - diag(diag(A)) es igual a A.
24. Sea E=[-1 2 4; 1 -2 -4; -1 2 4], expresar la matriz F con las columnas de
izquierda a derecha de la matriz E.
............
............
............
25. Expresar E como la suma de una matriz simtrica y una antisimtrica
SIMTRICA ANTISIMTRICA
............
............
............
....................................
26. Determinar si E es idempotente.
27. Sea A=[-3 -6 2; 2 4 -1; 2 3 0], expresar la matriz A con las filas de arriba abajo.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 51
............
............
............
28. Determinar si la matriz A es involuta
............
............
............
29. Hallar el rango de la matriz A .
30. Hallar el determinante de la matriz A
31. Hallar la inversa de la matriz A
............
............
............
32. Sea C=[1 1 3; 5 2 6; -2 -1 -3], determinar si la matriz C es nilpotente.
............
............
............
33. Hallar el determinante de la matriz C
34. Hallar el rango de la matriz C .
35. Resolver el sistema:
2 a + 3 b + c = 6
4 a + b + 2 c = 7
6 a + b + 7 c = 4
Mediante la funcin inv.
36. Resolver el siguiente sistema de ecuaciones.
2x + 0y + 5z = 100
3x + 5y + 9z = 251
1x + 5y + 7z = 301
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
52 | L. Maria Pimentel Herrera
Polinomios
conv Producto de polinomios.
deconv Divisin de polinomios.
poly Polinomio con races especificado.
polyder Polinomio derivados
polyeig Polinomial problema de valores propios
polyfit Ajuste de la curva polinmica
polyint Integrar polinomio analticamente
polyval Evaluacin Polynomial.
polyvalm Matriz de evaluacin polinomio
residue Convertir entre la expansin de l as fracciones parciales y polinomio de coeficientes
roots Races del polinomio.
>>p=[1 0 2 0 3] % Polinomio x^4+2*x^2+3
p =
1 0 2 0 3
>>q=[2 1 0] % Polinomio 2*x^2+x
q =
2 1 0
>>polyval(p,-1) % Evaluacin del polinomio x^4+2x^2+3 en x=-1
ans =
6
>>pro=conv(p,q) % Producto de los polinomios p y q
pro =
2 1 4 2 6 3 0
>>deconv(pro,p) % Cociente entre pro y p; obviamente el resultado es q
ans =
2 1 0
>>roots(pro) % Races del polinomio pro
ans =
0
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 53
0.6050+1.1688i
0.6050-1.1688i
-0.6050+1.1688i
-0.6050-1.1688i
-0.5000
>>poly([i -i 1/2 pi]) % Polinomio mnico que tiene por races a los
% nmeros i, -i, 0.5 y pi
ans =
-3.6416 2.5708 -3.6416 1.5708
>> A=[4,2;3,3]; p=poly(A)
p =
1 -7 6
El resultado son los coeficientes del polinomio caracterstico ordenado de acuerdo a las potencias decrecientes de la variable, es decir:
P() = 2 7+6
Otra forma de calcular el polinomio caracterstico es usando el comando:
vpa(polynsym(p)), donde n indica el nmero de cifras decimales con que se quiere obtener los coeficientes del polinomio.
>> vpa(poly2sym(p))
ans =
x^2-7.*x+6
Expresa el polinomio caracterstico en la variable x.
Sean p(x)=3x2+6x+9 y q(x)=x2+2x)
>> px=[3 6 9], qx=[1 2 0]
> polyder(px,qx)
ans =
12 36 42 18
>> conv(px,qx)
ans =
3 12 21 18 0
>> prod=conv(px,qx)
prod =
3 12 21 18 0
>> deconv(prod,px)
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
54 | L. Maria Pimentel Herrera
ans =
1 2 0
>> polyder(prod)
ans =
12 36 42 18
>>polyint(qx)
ans =
1/3 1 0 0
Ejercicios
37. Definir los siguientes polinomios: p(x)=3*x^4+5*x^3+2*x^2+8*x+6 y
q(x)=6*x^4+2*x^3+x^2+7*x+8 y hallar
a. El valor del polinomio p(-1) y q(-3)
p(-1)=
q(-3)=
b. r(x) es el producto de los polinomios p(x) y q(x)
r(x)=
c. d(x) es el cociente de los polinomios p(x) y q(x)
d(x)=
d. las soluciones de los polinomios p(x) y q(x)
p(x)
x1= x2= x3= x4=
q(x)=
x1= x2= x3= x4=
38. Hallar las soluciones de los polinomios
a. f(s)=s4 + 3s3 15s2 2s + 9
s1= s2= s3= s4=
b. g(s)=s4 + 1
s1= s2= s3= s4=
c. h(x)=x3 + 5x2 2
x1= x2= x3=
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 55
lgebra Lineal
Matriz de logaritmos y exponenciales
expm Matriz exponencial.
logm Logaritmo Matrix.
sqrtm Raz matriz cuadrada.
>> A = [1 1 0; 0 0 2; 0 0 -1 ]
A =
1 1 0
0 0 2
0 0 -1
>> Y=expm(A)
Y =
2.7183 1.7183 1.0862
0 1.0000 1.2642
0 0 0.3679
>> logm(Y) % A=logm(Y)
ans =
1.0000 1.0000 0.0000
0 0 2.0000
0 0 -1.0000
>> X=[10 7; 15 22]
X =
10 7
15 22
>> Y=sqrtm(X)
Y =
2.8347 0.9575
2.0518 4.4761
>> Y*Y % X=Y*Y
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
56 | L. Maria Pimentel Herrera
ans =
10.0000 7.0000
15.0000 22.0000
Anlisis de matriz
cond Nmero de condicin con respecto a la inversin.
condeig Nmero de condicin con respecto a los valores propios.
det Determinante de la matriz.
norm Vector y matriz de las normas.
normest 2-estimacin de la norma.
null Espacio nulo.
orth Rango de espacio de la matriz.
rank Rango de la matriz.
rcond Matriz de estimacin de nmero de condicin de reciprocidad.
rref Reduccin de forma escalonada por fila.
subspace ngulo entre dos subespacios.
trace Suma de los elementos de la diagonal.
>> A=magic(3);
>> trace(A)
ans=
.
>>rank(A)
ans=
.
>>detA)
ans=
.
>>norm(A)
ans=
.
>> B=magic(4);
A = 8 1 6 3 5 7 4 9 2
B = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 57
>> rank(B)
ans =
3
>> rref(B)
ans =
1 0 0 1
0 1 0 3
0 0 1 -3
0 0 0 0
>> c=sqrt(2); C=[1/c 0 -1/c; 1/c 0 1/c; 0 1 0], OR=orth(C)
C =
0.7071 0 -0.7071
0.7071 0 0.7071
0 1.0000 0
OR =
-0.7071 0 0.7071
-0.7071 0 -0.7071
0 -1.0000 0
>> OR'*OR
ans =
1.0000 0 0.0000
0 1.0000 0
0.0000 0 1.0000
Valores propios y valores singulares
balance Diagonal de escala para mejorar la precisin de valores propios
cdf2rdf Cambio de forma diagonal complejo para bloquear real forma diagonal
eig Valores y vectores propios
eigs Valores propios ms grande y vectores propios de la matriz
gsvd Generalizada descomposicin de valor singular
hess Hessenberg forma de matriz de
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
58 | L. Maria Pimentel Herrera
ordeig Valores propios de las matrices de quasitriangular
ordgz Reordenar valores propios en QZ factorizacin
ordschur Reordenar valores propios en la factorizacin Schur
poly Polinomio con races especificado
polyeig Polynomial problema de valores propios
rsf2csf Cambio de forma real de Schur para formar complejos de Schur
schur Descomposicin de Schur
sqrtm Raz matriz cuadrada
ss2tf Cambio de estado de los parmetros de filtro de espacio para la transferencia de forma funcin de
svd Descomposicin del valor singular
vds Encuentre los valores singulares y de vectores
P()=det(A-I) donde A es una matriz cuadrada, I es la matriz identidad y es un parmetro. Al desarrollar det(A-I) obtenemos el polinomio caracterstico.
>>A=[5 4 2; 4 5 2; 2 2 2];
>> p=poly(A)
p =
1 -12 21 -10
Al resolver det(A-I)=0 obtenemos los valores propios o caractersticos del polinomio, y
(A-I)v=0 obtenemos los vectores propios asociados a los valores propios.
En Matlab tenemos eig, si [V,D]=eig(A) produce una matriz diagonal D de los valores propios y una matriz V de columnas que corresponden a los vectores propios. De tal forma que se cumpla: AV=VD.
>> [V,D]=eig(A)
V =
-601/1157 503/941 2/3
112/1181 -2107/2850 2/3
383/451 1514/3697 1/3
D =
1 0 0
0 1 0
0 0 10
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 59
Anlisis de datos
Operaciones bsicas
prod Producto de elementos de la matriz.
sum Suma de los elementos de la matriz.
cumprod Producto acumulado.
cumsum Suma acumulada.
sort Ordenar elementos de la matriz en or den ascendente o descendente.
sortrows Ordenar filas en orden ascendente.
>>x=[15 2 -7 5 6 3 4 1 8 9 21 17.5];
>>sum(x) % suma de los elementos de un vector.
ans =
84.50
>>cumsum(x) % devuelve el vector suma acumulativa de los elementos de un vector.
% producto de los elementos de un vector
.x =
Columns 1 through 7
15.00 17.00 10.00 15.00 21.00 24.00 28.00
Columns 8 through 12
29.00 37.00 46.00 67.00 84.50
>>x=[15 3 4 1 8 9];
>> cumprod(x) % devuelve el vector producto acumulativo de los elementos de un vector.
ans =
15.00 45.00 180.00 180.00 1440.00 12960.00
>> [y,l]=sort(x) % ordena de menor a mayor los elementos de un vector x.
y =
1.00 3.00 4.00 8.00 9.00 15.00
l =
4.00 2.00 3.00 5.00 6.00 1.00 5 6 8 9 15
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
60 | L. Maria Pimentel Herrera
Estadstica Descriptiva
Puede utilizar las funciones de MATLAB siguiente para calcular las estadsticas descriptivas para los datos.
Nota: Para la matriz de datos, estadsticas descriptivas para cada columna se calculan de forma independiente.
corrcoef los coeficientes de correlacin
cov covarianza.
max mayor valor.
mean promedio o valor medio.
median mediana.
min menor valor.
mode los valores ms frecuentes en la matriz.
std desviacin estndar.
var varianza.
>> x=[15 2 -7 5 6 3 4 1 8 9 21 17.5];
>>[xm,im]=max(x)
% mximo elemento de un vector. Devuelve el valor mximo xm y la posicin que ocupa im.
xm =
21
im =
11
>>[xmi,mi]=min(x)
% mnimo elemento de un vector. Devuelve el valor mnimo xmi y la posicin que ocupa mi.
xmi =
-7
mi =
3
>> mean(x) % valor medio de los elementos de un vector.
ans =
169/24
>> std(x) % desviacin tpica.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 61
ans =
644/83
Nota: en realidad estas funciones se pueden aplicar tambin a matrices, pero en este caso se aplican por separado a cada columna de la matriz, dando como valor de retorno un vector resultante de aplicar la funcin a cada columna de la matriz considerada como vector. Si estas funciones se quieren aplicar a las filas de la matriz basta aplicar dichas funciones a la matriz transpuesta.
Derivadas e integrales
cumtrapz Acumulativa de integracin numrica trapezoidal
del2 Integral laplaciano
diff Las diferencias y aproximar los derivados
gradient Gradiente numrico
int Integral
polyder Polinomio derivados
polyint Integrar polinomio analticamente
trapz La integracin numrica trapezoidal
>>f='sin(x)' % Funcin sin(x) definida mediante una cadena de caracteres
f =
sin(x)
, >>diff(f) % calcular derivadas
ans = cos(x)
>>diff(f,2) % Derivada segunda de f
ans = -sin(x)
>>int('log(x)') % Primitiva de la funcin logaritmo
ans = x*log(x)-x
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
62 | L. Maria Pimentel Herrera
>>diff('x*log(x)-x') % Comprobacin
ans = log(x)
Ejercicios
Las notas obtenidas por 10 alumnos en Fsica 1 y Fsica 2 son:
F 1 11 9 13 10 8 12 10 15 10 9 F 2 13 11 15 12 10 14 12 9 12 16
Completar el cuadro
F 1 F 2 Mxima nota Mnima nota Acumulado (suma) Media aritmtica Desviacin estndar Varianza Covarianza
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 63
Programacin MATLAB es una aplicacin que se puede programar muy fcilmente. Se comenzar viendo las bifurcaciones y bucles, y la lectura y escritura interactiva de variables, que son los elementos bsicos de cualquier programa de una cierta complejidad.
Es posible hacer una coleccin de comandos y agruparlos en un archivo de tipo texto y de extensin m (.m) llamado archivo-m. Estos archivos pueden ser scripts o funciones. El script es un archivo-m que contiene una se rie de comandos que se ejecutarn al ejecutar dicho archivo en MatLab. La funcin, es un archivo-m que permite la entrada y salida de argumentos adems de la ejecucin de comandos. Para crear un archivo-m se usa cualquier editor de textos, asegurndose de al macenar dicho archivo con la extensin (.m).
Para crear un archivo-M escogemos New del men File y seleccionamos M-file. Una vez guardado este archivo-M en el disco, Matlab ejecutar las rdenes en dicho archivo simplemente escribiendo su nombre (sin extensin) en l a ventana de comandos de Matlab.
Guardar el archivo con el nombre ejemplo.m
Defina o agregue la ruta donde esta guardando sus archivos:
>> path(path,'C:\Documents and Settings\Mary\Escritorio\matlab\CEPSUNI2010')
Ejecute el programa de la lnea de comando de la ventana Command Window:
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
64 | L. Maria Pimentel Herrera
>> ejemplo
Inicio clculos
La traza vale: 34.000000
Ordenes de gestin de archivos what devuelve un listado de todos los archivos-M del directorio actual.
dir lista todos los archivos en el directorio o carpeta actual.
ls contenido de la carpeta, igual a dir.
type test visualiza el archivo-M test.m en la ventana de comando.
delete test elimina el archivo-M test.m.
cd path cambia al directorio o carpeta dada por path.
chdir path lo mismo que cd path.
cd muestra el directorio o carpeta de trabajo presente.
chdir lo mismo que cd.
pwd lo mismo que cd.
which test visualiza el camino del directorio de test.m.
Bifurcaciones y bucles MATLAB posee un l enguaje de pr ogramacin que como cualquier otro lenguaje dispone de sentencias para realizar bifurcaciones y bucles. Las bifurcaciones permiten realizar una u otra operacin segn se cumpla o no una determinada condicin.
Ejemplos grficos de bifurcaciones.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 65
Los bucles permiten repetir las mismas o anlogas operaciones sobre datos distintos. Mientras que en C /C++/ Java el "cuerpo" de estas sentencias se determinaba mediante llaves {...}, en MATLAB se utiliza la palabra end con anloga finalidad. Existen tambin algunas otras diferencias de sintaxis.
Bucles con control al principio y al final.
La Figura muestra dos posibles formas de bucle, con el control situado al principio o al final del mismo. Si el control est situado al comienzo del bucle es posible que las sentencias no se ejecuten ninguna vez, por no haberse cumplido la condicin cuando se llega al bucle por primera vez. Sin embargo, si la condicin est al final del bucle las sentencias se ejecutarn por lo menos una vez, aunque la condicin no se cumpla. Muchos lenguajes de programacin disponen de bucles con control al principio (for y while en C /C++/Java) y al final (do while en C/C++/Java). En MATLAB no hay bucles con control al final del bucle, es decir, no existe construccin anloga a do ... while.
Sentencia if En su forma ms simple, la sentencia if se escribe en la forma siguiente:
if condicion sentencias end
Existe tambin la bifurcacin mltiple, en l a que pueden concatenarse tantas condiciones como se desee, y que tiene la forma:
if condicion1 bloque1 elseif condicion2 bloque2 elseif condicion3 bloque3
else % opcin por defecto para cuando no se cumplan las condiciones 1,2,3 bloque4
end donde la opcin por defecto else puede ser omitida: si no est presente no se hace nada en caso de que no se cumpla ninguna de las condiciones que se han chequeado.
Una observacin muy importante: la condicin del if puede ser una condicin matricial, del tipo A==B, donde A y B son matrices del mismo tamao. Para que se considere que la condicin se cumple, es necesario que sean iguales dos a dos todos los elementos de las matrices A y B (aij=bij, 1im, 1jn). Basta que haya dos elementos aij y bij diferentes para que las matrices ya no sean iguales, y por tanto las sentencias
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
66 | L. Maria Pimentel Herrera
del if no se ejecuten. Anlogamente, una condicin en la forma A=B exige que todos los elementos sean diferentes dos a dos (aijbij, 1im, 1jn). Bastara que hubiera dos elementos aij y bij iguales para que la condicin no se cumpliese. En resumen:
if A==B exige que todos los elementos sean iguales dos a dos if A~=B exige que todos los elementos sean diferentes dos a dos como se ha di cho, MATLAB dispone de funciones especiales para ayudar en el chequeo de co ndiciones matriciales. Por ejemplo, la funcin isequal(A, B) devuelve un uno si las dos matrices son idnticas y un cero en caso de que difieran en algo.
Ejemplo 1. Dados dos nmeros a y b si a es mayor que b entonces intercambiar los
valores.
if a > b
tmp=a; a=b; b=tmp;
end
>> a=16; b=-7;
>> prog01
>> a
a =
-7
>> b
b =
16 Ejemplo 2. Si el nmero dado es 7 entonces lo cambia por cero, en caso contrario
lo cambia a 1.
>>n=17;
if n==17 n=0 else n=1 end;
>>prog02
n =
0
prog02.m
Donde el 0 proviene de entrar al primer if, y el 1, de entrar al else del segundo if.
prog01.m
Si a es mayor que b entonces intercambia el valor de las variables.
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 67
>>n=7
>>prog02
n =
1
Ejemplo 3. Si el nmero dado es diferente de cero entonces lo cambia por uno, en
caso contrario lo cambia a 3.
if m~=0 m=1 else m=3 end; >>m=9
>>prog03
m =
1
>>m=0
>>prog03
m =
3
Ejemplo 4. Si nota es mayor o i gual a 13 i mprimir aprobado en ca so contrario
desaprobado.
if nota>=13 fprintf('aprobado\n') else fprintf('desaprobado\n') end; >>nota=11
>>prog04
desaprobado
>>nota=15
>>prog04
aprobado
Ejemplo 5. Imprimir una matriz magic de orden:
n si el resto de dividir entre 2 es diferente de cero,
prog03.m
Donde el 1proviene de entrar al primer if, y el 3, de entrar al else del segundo if.
prog04.m
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
68 | L. Maria Pimentel Herrera
n-2 si el resto de dividir entre 3 es diferente de cero,
n-1 en otro caso.
if rem(n,2) ~= 0
M = magic(n) elseif rem(n,3) ~= 0
M = magic(n-2) else
M = magic(n-1) end
>>n=7
>>prog05
>>n=8
>>prog05
>>n=6
>>prog05
Sentencia switch
La sentencia switch realiza una funcin anloga a un conjunto de if...elseif concatenados. Su forma general es la siguiente:
switch switch_expresion case case_expr1,bloque1 case {case_expr2, case_expr3, case_expr4,...} bloque2 ... otherwise, % opcin por defecto bloque3 end
prog05.m
UNIVERSIDAD NACIONAL DE INGENIERA Centro de Extensin y Proyeccin Social
L. Maria Pimentel Herrera | 69
Al principio se evala la switch_expresion, cuyo resultado debe ser un nmero escalar o una ca dena de ca racteres. Este resultado se compara con las case_expr, y se ejecuta el bloque de se ntencias que corresponda con ese resultado. Si ninguno es igual a sw itch_expresion se ejecutan las sentencias correspondientes a ot herwise. Segn puede v erse en el ejemplo anterior, es posible agrupar varias condiciones dentro de unas llaves (const