Source: https://es.scribd.com/doc/34906816/Tutorial-Matlab
Timestamp: 2016-02-13 20:08:25+00:00

Document:
SubirSign inJoinBooksAudiobooksComicsSheet MusicScribd Selects BooksHand-picked favorites from our editorsScribd Selects AudiobooksHand-picked favorites from our editorsScribd Selects ComicsHand-picked favorites from our editorsScribd Selects Sheet MusicHand-picked favorites from our editorsTop BooksWhat's trending, bestsellers, award-winners & moreTop AudiobooksWhat's trending, bestsellers, award-winners & moreTop ComicsWhat's trending, bestsellers, award-winners & moreTop Sheet MusicWhat's trending, bestsellers, award-winners & moreCategoriesArts & IdeasBiography & MemoirBusiness & LeadershipChildren'sComputers & TechnologyCooking & FoodCrafts & HobbiesFantasyFiction & LiteratureHappiness & Self-HelpHealth & WellnessHistoryHome & GardenHumorLGBTMystery, Thriller & CrimePolitics & EconomyReferenceReligionRomanceScience & NatureScience FictionSociety & CultureSports & AdventureTravelYoung AdultCategoriesArts & IdeasBiography & MemoirBusiness & LeadershipChildren'sComputers & TechnologyCooking & FoodFantasyFiction & LiteratureHappiness & Self-HelpHealth & WellnessHistoryHome & GardenHumorLGBTMystery, Thriller & CrimePolitics & EconomyReferenceReligionRomanceScience & NatureScience FictionSociety & CultureSports & AdventureTravelYoung AdultCategoriesAdaptationsChildren’sCrime & MysteryFictionHumorMangaNonfictionRomanceSciFi, Fantasy & HorrorSuperheroesYoung AdultPublishersArcanaArchie ComicsBOOM! StudiosDynamiteIDW PublishingKingstone ComicsMarvel ComicsSpace Goat ProductionsTop Cow ComicsTop Shelf ProductionsValiant Comics ZenescopeDifficultyBeginnerIntermediateAdvancedMixedInstrumentBrassDrums & PercussionGuitar, Bass, and FrettedPianoStringsVocalWoodwindsGenreClassicalCountryFolkJazz & BluesMovies & MusicalsPop & RockReligious & HolidayStandardsP. 1Tutorial MatlabTutorial Matlab|Views: 987|Likes: 58Publicado pormirypauMore info:Published by: mirypau on Jul 27, 2010Copyright:Attribution Non-commercialAvailability:Read on Scribd mobile: iPhone, iPad and Android.download as PDF, TXT or read online from ScribdFlag for inappropriate content|Agregar a la colecciónSee moreSee lesshttps://es.scribd.com/doc/34906816/Tutorial-Matlab07/07/2013pdftextoriginalSectionsFuncionales especiales y elementalesMétodos numéricos no linealesEstadística y análisis de FourierOperaciones algebráicas y lógicas2.8 NON LINEAR CONTROL DESIGN TOOLBOX2.9 NAG FOUNDATION TOOLBOX3. INICIANDO MATLAB4. USO DE COMANDOS4.2 Instrucciones de MATLAB y Variables4.3 Obteniendo Información del Espacio de Trabajo4.4 Variables Permanentes4.6 Saliendo y Guardando el Espacio de Trabajo4.7 Manipulación de Vectores y Matrices4.8 Operaciones de Matrices4.9 Operaciones de Arreglos4.10 Ejemplos: Operaciones Aritméticas5. PROGRAMANDO CON MATLAB5.1 Generalidades5.1.1 Archivos-M: Comandos y Funciones5.1.2 Otras funciones5.1.3 Declaración function5.2 Operadores relacionales5.3 Operadores lógicos5.4 Caracteres especiales5.5 Control de flujo5.5.1 Declaración FOR simple5.5.2 Declaración FOR anidada5.5.3 Declaración WHILE5.5.4 Declaraciones IF, ELSE, ELSEIF y BREAK5.6.1 Creación de una matriz5.6.2 Cambio del orden de una matriz: reshape5.6.3 Modificación individual de elementos5.6.4 Modificaciones adicionales de una matriz5.7.1 Declaración fopenEjemplo5.7.2 Declaración fclose5.7.3 Declaración fread5.7.4 Declaración fwrite5.7.5 Declaración fprintf5.8 Variables globales5.9 Vectorización de algoritmos y estructuras (for, while)5.10 Gráficas en Dos DimensionesComando PlotSímbolo ColorSímbolo Estilo de línea5.10.6 Comandos gráficos5.11 Gráficos en 3 dimensiones5.12 Archivos de disco5.12.1 Manipulación de Archivos de Disco5.12.2 Ejecutando Programas Externos5.12.3 Importando y Exportando Datos5.13 INDICE ALFABETICO6. SIMULINK6.1 Acelerador de Simulink6.2 Generador de código-C en Simulink7. COMANDOS DE MATLAB7.1 General purpose commands:Control System Toolbox Commands:8. APLICANDO MATLAB AL CONTROL DE PROCESOS8.1 Respuesta en el dominio del tiempo8.2 Respuesta en el dominio de la frecuencia8.3 Lugar de las raíces8.4 Controladores PID9. TRUCOS EN MATLAB®Paper semilogarítmico gratis: papelbod.mIng Sergio Casillas de la Torre Tutorial de MatlabTUTORIAL DE MATLAB
TUTORIAL DE MATLAB 1. ¿QUÉ ES MATLAB? 1.1 Uso de Matrices 1.2 Origen de MatLab 1.3 Plataformas 1.4 Productos 2. LIBRERÍA DE APLIC ACIONES DE MATLAB 2.1 SIGNAL PROCESSING TOOLBOX 2.2 THE MATLAB C MATH LIBRARY 2.2.1 Desarrollo de aplicaciones utilizando la MATLAB C Math Library 2.2.2 Utilización de MATLAB y de su compilador 2.2.3 Velocidad y Precisió n 2.2.4 Lista parcial de funciones Funciones matemáticas Funcionales especiales y elementales Algebra lineal numérica Polinomios e interpolación Métodos numéricos no lineales Estadística y análisis de Fourier Operaciones algebráicas y lógicas 2.2.5 Utilidades 2.2.6 Requerimientos 2.3 THE MATLAB COMPILER TOOLBOX 2.3.1 Generación Automática de ficheros MEX. 2.3.2 Rendimiento del compilador 2.3.3 Opciones de ajuste del rendimiento 2.3.4 Requerimientos del sistema 2.3.5 Limitaciones del código compilado 2.4 SYMBOLIC MATH TOOLBOX 2.5 OPTIMIZATION TOO LBOX 2.6 IMAGE PROCESSING TOOLBOX 2.7 Neural Network Toolbox
1 4 5 5 5 5 7 7 7 8 8 9 9 9 9 9 9 10 10 10 10 10 11 11 12 12 12 13 13 14 15 16
Ing Sergio Casillas de la Torre Tutorial de Matlab
2.8 NON LINEAR CONTROL DESIGN TOOLBOX 2.9 NAG FOUNDATION TOOLBOX 3. INICIANDO MATLAB 4. USO DE COMANDOS 4.2 Instrucciones de MATLAB y Variables 4.3 Obteniendo Información del Espacio de Trabajo 4.4 Variables Permanentes 4.6 Saliendo y Guardando el Espacio de Trabajo 4.7 Manipulación de Vectores y Matrices 4.8 Operaciones de Matrices 4.9 Operaciones de Arreglos 4.10 Ejemplos: Operaciones Aritméticas 5. PROGRAMANDO CON M ATLAB 5.1 Generalidades 5.1.1 Archivos -M: Comandos y Funciones 5.1.2 Otras funciones 5.1.3 Declaración function 5.2 Operadores relacionales 5.3 Operadores lógicos 5.4 Caracteres especiales 5.5 Control de flujo 5.5.1 Declaración FOR simple 5.5.2 Declaración FOR anidada . 5.5.3 Declaración WHILE 5.5.4 Declaraciones IF, ELSE, ELSEIF y BREAK 5.6.1 Creación de una matriz 5.6.2 Cambio del orden de una matriz: reshape 5.6.3 Modificación individual de elementos 5.6.4 Modificaciones adicionales de una matriz 5.7.1 Declaración fopen Ejemplo 5.7.2 Declaración fclose 5.7.3 Declaración fread 5.7.4 Declaración fwrite 5.7.5 Declaración fprintf 5.8 Variables globales 5.9 Vectorizació n de algoritmos y estructuras (for, while) 5.10 Gráficas en Dos Dimensiones COMANDO PLOT Símbolo Color Símbolo Estilo de línea 5.10.6 Comandos gráficos
17 18 20 20 22 23 23 23 24 25 28 29 33 33 33 37 41 41 42 43 44 44 45 46 47 50 50 50 51 57 57 57 57 58 58 58 59 60 60 60 61 63
5.11 Gráficos en 3 dimensiones 5.12 Archivos de disco 5.12.1 Manipulación de Archivos de Disco 5.12.2 Ejecutando Programas Externos 5.12.3 Importando y Exportando Datos 5.13 INDICE ALFABETICO 6. SIMULINK 6.1 Acelerador de Simulink 6.2 Generador de código -C en Simulink 7. COMANDOS DE MATLAB 7.1 General purpose commands: Control System Toolbox Commands : 8. APLICANDO MATLAB AL CONTROL DE PROCES OS 8.1 Respuesta en el dominio del tiempo 8.2 Respuesta en el dominio de la frecuencia 8.3 Lugar de las raíces 8.4 Controladores PID 9. TRUCOS EN MATLAB® Paper semilogarítmico gratis: papelbod.m
66 73 73 73 73 74 75 77 77 78 78 81 86 86 91 95 97 99 99
MatLab es un programa interactivo para computación numérica y visualización de datos. Es ampliamente usado por Ingenieros de Control en el análisis y diseño, posee además una extraordinaria versatilidad y capacidad para resolver problemas en matemática aplicada, física, química, ingeniería, finanzas y muchas otras aplicaciones. Está basado en un sofisticado software de matrices para el análisis de sistemas de ecuaciones. Permite resolver complicados problemas numéricos sin necesidad de escribir un programa. MATLAB es un entorno de computación y desarrollo de aplicaciones totalmente integrado orientado para llevar a cabo proyectos en donde se encuentren implicados elevados cálculos matemáticos y la visualización gráfica de los mismos. MATLAB integra análisis numérico, cálculo matricial, proceso de señal y visualización gráfica en un entorno completo donde los problemas y sus soluciones son expresados del mismo modo en que se escribirían tradicionalme nte, sin necesidad de hacer uso de la programación tradicional. El nombre de MATLAB proviene de la contracción de los términos MATrix LABoratory y fue inicialmente concebido para proporcionar fácil acceso a las librerías LINPACK y EISPACK, las cuales repre sentan hoy en dia dos de las librerías más importantes en computación y cálculo matricial. MATLAB es un sistema de trabajo interactivo cuyo elemento básico de trabajo son las matrices. El programa permite realizar de un modo rápido la resolución numérica de problemas en un tiempo mucho menor que si se quisiesen resolver estos mismos problemas con lenguajes de programación tradicionales como pueden ser los lenguajes Fortran, Basic o C. MATLAB goza en la actualidad de un alto nivel de implantación en escuelas y centros universitarios, así como en departamentos de investigación y desarrollo de muchas compañías industriales nacionales e internacionales. En entornos universitarios, por ejemplo, MATLAB se ha convertido en una herramienta básica, tanto para los pro fesionales e investigadores de centros docentes, como una importante herramienta para la impartición de cursos universitarios, tales como sistemas e ingenieria de control, álgebra lineal, proceso digital de imagen, señal, etc. En el mundo industrial, MATLAB está siendo utilizado como herramienta de investigación para la resolución de complejos problemas planteados en la realización y aplicación de modelos matemáticos en ingeniería. Los usos más característicos de la herramienta los encontramos en áreas de c omputación y cálculo numérico tradicional, prototipaje algorítmico, teoría de control automático, estadística, análisis de series temporales para el proceso digital de señal.
Ing Sergio Casillas de la Torre Tutorial de Matlab MATLAB dispone también en la actualidad de un amplio abanico de programas de apoyo especializados, denominados Toolboxes, que extienden significativamente el número de funciones incorporadas en el programa principal. Estos Toolboxes cubren en la actualidad prácticamente casi todas las áreas principales en el mundo de la ingeniería y la simulación, destacando entre ellos el 'toolbox' de proceso de imágenes, señal, control robusto, estadística, análisis financiero, matemáticas simbólicas, redes neurales, lógica difusa, identificación de sistemas, simulación de sistemas dinámicos, etc. Ad emás también se dispone del programa Simulink que es un entorno gráfico interactivo con el que se puede analizar, modelizar y simular la dinámica de sistemas no lineales.
MatLab emplea matrices porque con ellas se puede describir infinidad de cosas de una forma altamente flexible y matemáticamente eficiente. Una matriz de pixeles puede ser una imagen o una película. Una matriz de fluctuaciones de una señal puede ser un sonido o una voz humana. Y tal vez más significativamente, una matri z puede describir una relación lineal entre los componentes de un modelo matemático. En este último sentido, una matriz puede describir el comportamiento de un sistema extremadamente complejo. Por ejemplo una matriz puede representar el vuelo de una avión a 40.000 pies de altura, o un filtro digital de procesamiento de señales.
MatLab fue originalmente desarrollado en lenguaje FORTRAN para ser usado en computadoras mainframe. Fue el resultado de los proyectos Linpack y Eispack desarrol lados en el Argonne National Laboratory. Su nombre proviene de MATrix LABoratory. Al pasar de los años fue complementado y reimplementado en lenguaje C. Actualmente la licencia de MatLab es propiedad de MathWorks Inc .
La empresa MathWorks ofrece MatLab como su principal producto para computación numérica, análisis y visualización de datos. También ofrece Simulink 5
Ing Sergio Casillas de la Torre Tutorial de Matlab como un anexo a MatLab y que interactua con él en lenguaje de MatLab y lenguaje de bajo nivel C. Simulink es usado para simulación modelado no lineal avanzado. Se ofrecen además numerosas herramientas especiales en "Toolboxes" para resolver problemas de aplicaciones específicas, por ejemplo control, procesamiento de señales, redes neurales, etc. Estas herramientas son c o l ecciones de rutinas escritas en MatLab.
2.1 SIGNAL PROCESSING TOOLBOX
Diseño de filtros FIR mediante el algorítmo óptimo de Parks-McClellan. Procesamiento de la transformada rápida de Fourier FFT, incluyendo la transformación para potencias de dos y su inversa, y transformada para no potencias de dos.
La MATLAB C Math Library proporciona al usuario la capacidad computacional de MATLAB en una libreria en formato objeto enlazable. El objetivo principal de la C Math Library es soportar el desarrol lo de aplicaciones 'stand alone' utilizando MATLAB y su compilador. Puede ser utilizada independientemente de MATLAB por programadores avezados en lenguaje C que necesiten prestaciones computacionales robustas y de alto rendimiento. Junto con el compilador de MATLAB, la C Math Library permitirá a los programadores de aplicaciones utilizar MATLAB para la creación de aplicaciones 'stand alone'. Para los usuarios clásicos de MATLAB, se elimina así cualquier necesidad de volver a reescribir algoritmos en lenguaje C para ser utilizada por programas externos. Para aquellos usuarios que sean nuevos en la tecnología MATLAB, esta tecnología ofrece una nueva vía para la reducción del tiempo de desarrollo y puesta a punto de aplicaciones. La MATLAB C Math Library propo rciona una amplia gama de funciones clásicas del programa MATLAB, proporcionadas como librerias objeto, incluyendo básicamente las siguientes categorías de funciones presentes en MATLAB y ficheros M compilados:
Ing Sergio Casillas de la Torre Tutorial de Matlab • • • • • • • Matrices elementales y manipulación de vectores. Matrices especiales. Estadística básica y análisis de datos. Polinomios e interpolación. Gestión de cadenas de caracteres. Entradas y Salida s. Gestión de memoria y errores.
La construcción y desarrollo de aplicaciones utilizando esta l ibrería es un proceso de amplias perspectivas una vez se tiene un dominio adecuado de su operativa. El producto está dividido en dos categorías (como librerías objeto): la librería (built-in library) contiene versiones de las funciones de MATLAB en lenguaje C del tipo numérico, lógico y utilidades. Por otra parte la librería de toolboxes (toolbox library) contiene versiones compiladas de la mayoría de ficheros M de MATLAB para cálculo numérico, análisis de datos y funciones de acceso a ficheros y matrices. En equipos UNIX estas librerias pueden ser igualmente obtenidas como librerías de tipo estático (static libraries) o bien como librerías compartidas (shared libraries). Respecto al mundo PC, estas librerías pueden obtenerse como DLL's en el entorno Microso ft Windows o como librerias compartidas en equipos Apple MacIntosh.
Para construir una aplicación del tipo 'stand alone' que incorpore código originalmente desarrollado como ficheros M de MATLAB , deberán seguirse los pasos siguientes: 1. Utilizar el compilador de MATLAB para convertir ficheros M en C mediante la utilización de la instrucción mcc -e (la cual es externa a MATLAB). 2. Compilar el código C fuente en código objeto utilizando un compilador ANSI C. 3. Enla zar el código resultante con la MATLAB C Math Library y con cualquier tipo de ficheros y programas específicos que hayan sido previamente definidos por el usuario. 8
Los algoritmos utilizados en la MATLAB C Math Library han sido desarrollados por un grupo de renombrados expertos en programación algorítmica de funciones de tipo matemático (algebra lineal y cálculo numérico). Las funciones de álgebra lineal han sido obtenidas de las librerias mundialmente reconocidas LINPACK y EISP ACK. La MATLAB C Math Library contiene más de 300 funciones numéricas, lógicas y de utilidad. Todas estas funciones le permitirán operar en datos de tipo escalar, vectorial o matricial con la misma facilidad sintáctica.
Funciones gamma, beta y elípticas. Transformación de sistemas de coordenadas. Matriz identidad y otras matrices elementales. Matrices de Hilbert, Toeplitz, Vandermonde, Hada mard, etc. Partes reales, imaginarias y complejas conjugadas. Funciones trigonométricas y de potencias.
Valores propios y descomposición de matrices. Funciones generales de evaluación de matrices. Dete rminantes, normas, rangos, etc. Matrices inversas y factorización de matrices. Matriz exponencial, logarítmica y raíces cuadradas.
Interpolación 1 - D y 2 -D. Construcción polinomial. Interpolación por splines cúbicos. Diferenciación de polinomios. Evaluación de polinomios. Multiplicación y división de polinomios. Residuos de polinomios y residuos.
Convolución 1 - D y 2 - D. Filtros digitales 1 - D y 2 -D . Transformadas de Fourier 1 - D y 2 -D y su inversa. Coeficientes de correlación y matrices de covarianza. Deconvolución. Magnitudes y ángulos de fase. Funciones max, min, sum, mean y otras funciones de estadística básica.
Suma, resta, multiplicación, división y potencias de matrices. Matrix traspuesta. Operadores lógicos AND, OR, NOT y XOR.
Gestión y mantenimiento de errores. Conversión de tipos de datos Fortran. Funciones de fecha y hora. Clasificación de matrices. Conversión de números a cadenas y viceversa.
El nuevo compilador de MATLAB - The MATLAB Compiler- permite crear código C optimizado procedente de ficheros M - M files- de MATLAB. Este compilador puede ser utilizado de dos modos: 1. Como un generador MEX automático. Pueden convertirse ficheros M en funciones C ejecutables que se ejecutaran desde dentro de MATLAB. Como un generador de código C fuente. 2. Pueden construirse aplicaciones que se ejecutaran independientemente de MATLAB. Estas aplicaciones externas requieren de la MATLAB C Math Library, que está disponible separadamente.
Mediante la conversión automática de ficheros M en código C fuente, el compilador MAT LAB elimina consumo de tiempo y la conversión manual de código. Todo el proceso de conversión, compilación y enlazado se inicia a través de una simple instrucción de MATLAB.
El compilador de MATLAB automatiza l a creación de ficheros MEX de C (MATLAB Ejecutables). Los ficheros MEX contienen código objeto que es dinámicamente enlazado como 'runtime' en el entorno MATLAB por el intérprete del programa. El proceso en cuestión se realiza en tres pasos: 1. El compilador de MATLAB traduce las funciones MATLAB en sus funciones equivalente en lenguaje C. 2. La instrucción MATLAB cmex llama al compilador y al enlazador del sistema para construir un fichero MEX objeto. 3. El intérprete de MATLAB enlaza automáticamente la función de MATLAB como 'runtime'. Mientras se efectúa una conversión de los ficheros M en ficheros MEX, el compilador realiza llamadas a las rutinas de la libreria C para muchas de las instrucciones contenidas en el propio núcleo de MATLAB. Existen algunas funciones, incluyendo las rutinas 'Handle Graphics', para las cuales se generan de nuevo llamadas 'callbacks' a MATLAB. Pueden convertirse convenientemente ficheros M en código fuente C para incorporarlos posteriormente en los ficheros externos desarrollados en lenguaje C, si ese es el caso. Esta opción es ideal para usuarios que quieren sacar la máxima ventaja de MATLAB desde cualquier otra aplicación o producir código C eficiente a partir de los algoritmos desarrollados con MATLAB. Los desarrollos 11
Ing Sergio Casillas de la Torre Tutorial de Matlab del tipo 'stand -alone' requieren para ello de la MATLAB C Math Library. Obsérvese que las funciones gráficas de MATLAB no están incluidas. Para construir aplicaciones 'stand-alone' se debería seguir los siguientes pasos: 1. Utilizar el compilador de MATLAB para convertir ficheros M en C con la instrucción externa mcc - e. 2. Compilar el código C fuente en código objeto utilizando un compilador C. 3. Enlazar el código resultante con las librerías matemáticas C de MATLAB y los ficheros específicos que dispongamos.
Mediante la compilación de los ficheros M se puede obtener un rendimiento significativo. La velocidad de mejora de este rendimiento, depende fuertemente de cada aplicación. En algunos casos el rendimiento puede mejorar hasta en 200 veces la ejecución si la comparamos con el modo de trabajo interpretado del programa. Las operaciones matriciales y vectoriales ejecutadas desde MATLAB ya están fuertemente optimizadas en su diseño. Sin embargo, mediante la utilización del compilador se obte ndrán significativas mejoras.
El compilador de MATLAB ofrece varias opciones que permiten generar el programa final de la forma más eficiente. Por ejemplo, Ud. puede directamente: • • • Tratar todas las variables en fic heros como datos enteros y/o reales. Utilizar una variable concreta como variable escalar, vectorial, entera, real o una combinación de estas. Desactivar el control de parámetros de entrada y el redimensionamiento dinámico de vectores.
Para utilizar el compilador de MATLAB para crear ficheros MEX se necesita la versión de MATLAB 4.2c y tener instalado uno de los siguientes compiladores de lenguaje C: PC/Microsoft Windows Metaware High C/C++ V.3.0 o sup erior. Watcom C V.10.0 o superior Power MacIntosh MetroWerks CodeWarrior C V.7
Ing Sergio Casillas de la Torre Tutorial de Matlab MPW MrC V.1.0b2 o PPCC version 1.0.5 680x0 MacIntosh MPW C Versión 3.4 UNIX y VMS
Ciertas instrucciones, como load y eval, no están soportadas por el compilador de MATLAB . Este no puede generar código de los diagramas de bloques de SIMULINK. Los toolboxes de MATLAB puede n incluir ficheros MEX y otros componentes que no son compilables.
El Toolbox de Matemática Simbólica, añade a MATLAB la capacidad de realizar cálculos simbólicos basados en MAPLE V © soportando además (The Extended Symbolic Ma th Toolbox) las librerías especializadas, y los programas realizados para este último. Entre otros, los principales tipos de operaciones soportados son los siguientes: • • • • • Algebra simbólica: Derivación, integración y simplificación de expresiones matemáticas. Algebra lineal exacta: Inversas, canónicas de matrices simbólicas. determinantes, autovalores y formas
Aritmética de precisión variable: Evaluación de expresiones matemáticas con diversos grados de precisión. Resolución de ecuaciones: Resolución numéri ca y simbólica de ecuaciones algebraicas y diferenciales. Funciones matemáticas especiales: Evaluación de la mayoría de las funciones utilizadas en matemáticas aplicadas.
Ing Sergio Casillas de la Torre Tutorial de Matlab kernel de MAPLE utilizando la Sintaxis y el estilo del lenguaje MATLAB. The Extended Symbolic Math Toolbox aumenta esta funcionalidad incluyendo todas las características de programación de MAPLE, y el acceso a los paquetes de funciones de más de veinte campos de las matemáticas especiales aplicadas. Es posible utilizar este Toolbox sin conocimiento previos de MAPLE, ya que los ficheros contenidos en él son totalmente autónomos. Sin embargo, si lo que se desea es obtener toda la potencia de cálculo del entorno, será necesario un amplio conocimiento del manejo y la programación de MAPLE
2.5 OPTIMIZATION TOO LBOX
El toolbox de optimización consta de un conjunto de funciones que resuelven problemas de extremos, con o sin condiciones, de funciones reales las cuales son generalmente multivariables y no lineales. Asimismo, posee funciones para la resolución de algunos tipos de problemas matriciales en extremos. Resulta conveniente para una comp rensión y mejor manejo de la toolbox poseer conocimientos básicos previos de análisis de funciones reales, matrices y teoría de extremos. Algunas de las áreas básicas que cubre este toolbox para MATLAB son las siguientes: • Cálculo de un extremo local (máxi mo o mínimo) de una función real f(x), en general multivariable y no lineal, sin imponer ninguna restricción o condición a la solución. Como caso particular, se incluye una rutina especial para problemas de mínimos cuadrados no lineales. Cálculo de un extremo local (máximo o mínimo) de una función real f(x), en general multivariable y no lineal, condicionado a que la solución satisfaga ciertas condiciones de desigualdad (g(x)<=0) y/o igualdad (g(x)=0). Problemas de aproximación a un conjunto de objetivos. Cálculo de soluciones de un sistema de ecuaciones continuas y, en general, no lineales. Solución de problemas minimax. Programación lineal. Programación cuadrática. Problemas de mínimos cuadrados no negativos.
Este Too lbox proporciona a MATLAB de un conjunto de funciones que amplia las capacidades del producto para realizar desarrollo de aplicaciones y de nuevos algoritmos en el campo del proceso y análisis de imagenes. El entorno matemático y de creación de MATLAB es ideal para el procesado de imágenes, ya que estas imágenes son, al fin y al cabo, matrices. Este toolbox incorpora funciones para: • • • • • Diseño de filtros. Mejora y retocado de imágenes. Análisis y estadística de imágenes. Operaciones morfológicas, geométricas y de color. Transformaciones 2D.
El proceso de imágenes es un campo de trabajo absolutamente crucial para aquellos colectivos e industrias que estén trabajando en áreas como diagnóstico médico, astronomía, geofísica, ciencia medioambientales, análisis de datos en laboratorios, inspección industrial, etc. Los programas actuales de procesado y análisis de imágenes se clasifican actualmente en dos categorías: librerías de bajo nivel para programadores profesionales y paquetes de aplicación con capacidades li mitadas de personalización. Ambos tipos de aplicaciones están, generalmente, pensados para tareas básicas de visualización de datos y 'rendering'. Sin embargo, muchos de ellos adolecen de la posibilidad de efectuar análisis numéricos de los mismos. El Image Processing Toolbox entra dentro de la categoría de familias de funciones que, desde el entorno de trabajo de MATLAB , permitirá al profesional efectuar una exploración exhaustiva y desde un punto de vista matemático de las imágenes y gráficos que se dese en tratar o analizar. Algunas de las funciones más importantes incluidas dentro de este toolbox son las siguientes: • • • • • • • • Análisis de imágenes y estadística. Diseño de filtros y recuperación de imágenes. Mejora de imágenes. Operaciones morfológicas. Definición de mapas de colores y modificación gráfica. Operaciones geométricas. Transformación de imágenes. Proceso de bloques
Este toolbox proporciona funciones para el diseño, inicialización, simulación y entrenamiento de los model os neuronales de uso más extendido en la actualidad: Perceptrón, redes lineales, redes de retropropagación, redes de base radial, aprendizaje asociativo y competitivo, aplicaciones autoorganizativas, aprendizaje de cuantización vectorial, redes de Elman y redes de Hopfield. Mediante la inclusión de un amplio abanico de funciones y procedimientos escritos para MATLAB, el usuario puede mediante el Neural Network Toolbox efectuar el diseño de arquitecturas complejas, combinando los modelos que ya estan proporc ionados por defecto en el toolbox. Asimismo, el usuario puede definir sus propias funciones de transferencia e inicialización, reglas de aprendizaje, funciones de entrenamiento y estimación de error para usarlas posteriormente con las funciones básicas. El toolbox, aporta las facilidades y prestaciones gráficas de MATLAB para el estudio del comportamiento de las redes: visualización gráfica de la matriz de pesos y vector de desplazamiento mediante diagramas de Hinton, representación de errores a lo largo del entrenamiento, mapas de superficie de error en función de pesos y vector de desplazamiento, etc. Estos gráficos resultan muy útiles en el estudio de la convergencia y estabilidad de los algoritmos de aprendizaje. Este toolbox incluye un manual de introducción al campo de las redes neuronales junto con una colección de demostraciones y aplicaciones muy didácticas, útiles para el estudio y la profundización en las cuestiones fundamentales de los paradigmas de redes neuronales básicos. Asimismo, se proporcionan las referencias bibliográficas más significativas referidas a los distintos modelos que aparecen en la aplicación. A pesar de que el estudio de las redes neuronales se inició ya hace algunas decadas, las primeras aplicaciones sólidas dentro de este campo no han tenido lugar hasta hace unos doce años y aun ahora constituyen un área de investigación en rápido desarrollo. Este toolbox tiene por tanto una orientación diferente a aquellos destinados a campos como el de sistemas de control u optimización dond e la terminología, fundamentos matemáticos y procedimientos de diseño estan ya firmemente establecidos y se han aplicado durante años. Este toolbox pretende que sea utilizado para la valoración y diseño de diseños neuronales en la industria y sobre todo en educación e investigación. Esta herramienta tiene el soporte de MATLAB 4.2c y SIMULINK. La librería de SIMULINK contiene modelos de capas de redes neuronales de cada tipo de neurona implementada en el toolbox de redes neuronales. Es posible por tanto dise ñar sistemas SIMULINK para simular redes neuronales creadas usando esta herramienta. Simplemente, las capas se conectan de acuerdo con la arquitectura de la red y se proporcionan como entrada a la caja de diálogo de cada capa la matriz de pesos apropiada y el vector de desplazamiento. Usando el generador de código C de SIMULINK es posible generar automáticamente el código correspondiente a un diseño neuronal. 16
Ing Sergio Casillas de la Torre Tutorial de Matlab Dentro de las aplicaciones básicas de este toolbox, cabe destacar aquellas que están orientadas a aquellas que se enmarcan dentro del campo de la industria aeroespacial y automoción (simulación, sistemas de control, autopilotaje), banca, defensa (reconocimiento de patrones, procesamiento de señales, identificación de imágenes, extracción de característi cas, compresión de datos), electrónica (control de procesos, análisis de errores, modelado no lineal, síntesis de voz, visión por ordenador), economía (análisis financiero, análisis predictivo), industria (control de procesos, identificación en tiempo real, sistemas de inspección), medicina, robótica (control de trayectorias, sistemas de visión), reconocimiento y síntesis del habla, telecomunicaciones (control de datos e imágenes, servicios de información automatizada, traducción del lenguaje hablado en tiempo real, diagnosis, sistemas de enrutamiento), etc. El toolbox contiene muchos ejemplos de algunas de estas aplicaciones.
Se trata del primer producto comercialmente disponible en la actualidad para el diseño de controladores automáticos en entornos de sistemas no lineales. Este nuevo toolbox está pensado para ser utilizado exhaustivamente por ingenieros que diseñan controladores para industrias avanzadas, destacando el sector del automóvil, ingenieria aeroespacial, control de procesos y empresas petroquímicas. Según indica Jim Tung, Vicepresidente del área de desarrollo de The MathWorks Group, Inc. "El proceso de aproximación tradicional en el diseño de controladores en sistemas no lineales ha sido hasta la fecha linealizarlos de algún modo para aplicar posteriomente un método de diseño lineal que requiere de importantes ajustes manuales. El toolbox NCD permite por primera vez a los ingenieros de control diseñar directamente sus controladores en un ambiente no lineal, obviando la aproximación lineal y otros procedimientos auxiliares que antes se necesitaban de modo imperativo. Los resultados ahora son de elevada calidad, controladores más robustos y un ciclo de diseño mucho más rápido. El toolbox NCD extiende, además, las prestaciones que incorpora SIMULINK, el entorno de desarrollo de diagramas de bloques para la modelación y análisis de sistemas dinámicos de The MathWorks, Inc. El usuario puede incluir uno o más bloques NCD en el sistema y describir posteriormente de modo totalmente gráfico las restricciones, tolerancias y límites de permisividad de cada uno de estos bloques. Los métodos avanzados de optimización y la simulación del proceso son posteriormente analizados y ajustados mediante la inclusión de unas c i e r t a s variables de contorno para poder obtener los tiempos de respuesta deseados. Este toolbox puede ser utilizado para ajustar una amplia variedad de controladores que se utilicen en un sistema, destacando los controladores PID, LQR, LQG y estructuras H infini to. El diseñador de sistemas puede utilizar el método de Montecarlo para el diseño y análisis de controladores robustos,
Ing Sergio Casillas de la Torre Tutorial de Matlab siempre que se detecten determinadas variaciones en los componentes del sistema. El toolbox NCD es un componente avanzado del entorno integrado de desarrollo que ofrecen a los especialistas los programas MATLAB y SIMULINK. Por ello, los diseñadores podrán beneficiarse de muchos de los toolboxes desarrollados para este entorno en materia de diseño de sistemas lineales. Por ejemplo, podrán utilizarse toolboxes para el análisis de sistemas lineales para el diseño inicial; posteriormente, podrán utilizarse modelos no lineales más sofisticados utilizando SIMULINK. Además, puede invocarse NCD para un mejor ajuste paramétrico y para la optimización de los controladores. Este toolbox se encuentra actualmente disponible para una amplia variedad de plataformas informáticas, destacando ordenadores personales tipo PC o Apple MacIntosh, numerosas estaciones UNIX y ordenadores Digital VAX VMS.
Este toolbox proporciona un acceso interactivo, desde dentro de MATLAB, a un amplio conjunto de funciones matemáticas y estadísticas contenidas en las clásicas NAG Fortran Libraries de la empresa The Numerical Algorithms Group Incorpora más de 200 ficheros M, los cuales cubren un amplio espectro de áreas de interés, entre las que cabe destacar optimización, ecuaciones diferenciales ordinarias y en derivadas parciales, cuadratura, estadística, etc. La NAG Foundation Toolbox añade también rutinas concretas para campos específicos tales como la resolución de problemas con condiciones de contorno, problemas de cuadratura adaptativa multidimensional, ajuste de curvas y superficies y el acceso a los algoritmos LAPACK para la resolución de ecuaciones lineales. Los nombre de las funciones han sido directamente tomados de las especificaciones de función clásica que añade The Numerical Algorithms Group para sus librerías. Como resultado de esto, aquellos usuarios de las librerías Fortran de NAG que a la vez sean usuarios de MATLAB, encontraran bastante cómodo acceder a las rutinas NAG utilizando la nomenclatura original. La NAG Foundation Toolbox es resultado de la colaboración corporativa que actualmente están llevando a cabo The MathWorks Group y T he Numerical Algoriths Group para proporcionar un rápido acceso desde MATLAB a un importante de rutinas matemáticas contenidas en la NAG Foundation Library. Actualmente, este toolbox incorpora 250 rutinas matemáticas. Algunas de las áreas de cobertura de l a NAG Foundation Toolbox son las siguientes: • • Ceros de polinomios Raíces de una o más ecuaciones de tipo trascendental.
Ing Sergio Casillas de la Torre Tutorial de Matlab • • • • • • • • • • • • • • • • • • Suma de series. Cuadraturas. Ecuaciones diferenciales ordinarias. Ecuaciones diferenciales en derivadas parciales. Estadística no paramétrica. Análisis de series temporales. Rutinas de clasificación. Aproximación de funciones especiales. Aproximación de curvas y superficies. Maximización y minimización de funciones. Factorización de matrices. Valores y vectores propios. Resolución de ecuaciones lineales simultáneas. Ecuaciones lineales (LAPACK). Estadística básica. Análisis de correlación y regresiones. Métodos multivariantes. Generación de números aleatorios.
Después de ejecutar el programa MatLab desde el sistema operativo empleado, por ejemplo haciendo doble click sobre el icono de MatLab en ambientes Windows, aparece el indicador de comandos el cual está listo para recibir instrucciones en lenguaje MatLab. Este indicador es de la siguiente forma: >> Al iniciar el uso de MatLab están disponibles dos comandos de ayuda y demostración. Para ejecutarlos se escribe el comando en la línea de comandos después del símbolo >> y se presiona la tecla Enter. Por ejemplo: >>help permite obtener una ayuda sobre los diferentes comandos de MatLab. >>demo hace una demostración de las diferentes aplicaciones de MatLab. Para cerrar o finalizar el uso de MatLab se usa el comando quit . >>quit
La primera forma de interactuar con MatLab es a través d e la línea de comandos. Puede ejecutarse un comando si este está escrito después del símbolo >> y se presiona la tecla Enter. MATLAB trabaja esencialmente con matrices numéricas rectangulares. La manera más fácil de entrar matrices pequeñas es enumerando los elementos de ésta de tal manera que: • • • los elementos estén separados por blancos ó comas. los elementos estén cerrados entre corchetes, [ ]. muestre el final de cada fila con ; (punto y coma).
Ejemplo: A = [ 1 2 3; 4 5 6; 7 8 9 ] resultaría en la matri z A= 123 456 789 MATLAB guarda esta matriz para utilizarla luego bajo el nombre de A. Si la matriz a introducir es muy grande se puede utilizar el siguiente formato:
A = [1 2 3 456 7 8 9] El comando load y la función fread pueden leer matrices generadas en sesiones anteriores ó generadas por otros programas. Ya que MatLab se basa en el álgebra de matrices como ejemplo crearemos una matriz. Estas pueden estar formadas por un sólo elementos (escalar), por una fila o una columna (vector) o por una serie de filas y columnas (matriz propiamente dicha). >>A=1 define A como un escalar de valor 1. Al definir A automáticamente MatLab presenta en pantalla su valor. A= 1 Para no presentar el valor de la variable creada, debe agregarse punto y coma (;) al final del comando. Después de crear una variable, puede presentarse escribiendo la variable después del prompt (>>). >>A Se pueden redefinir variables, por ejemplo: >>A=[1 2 3] define A como un vector de tres elementos, A(1)=1, A(2)=2 y A(3)=3. Estos elementos deben separase con espacios en blanco o comas (,). Para definir una matriz se deben separar las filas con punto y coma (;) o con retorno (Enter). >>A=[1 2 3; 4 5 6] o >>A=[1 2 3 4 5 6] ambos comandos producen el mismo efecto: A= 123 456 su valor en pantalla
Los elementos de una matriz pueden ser cualquier expresión de MATLAB. Ejemplo: x = [-1.3,sqrt(3),(1+2+3) *4/5] resultaría en x = -1.3000 1.7321 4.8000
Para añadir otra fila a la matriz A de arriba podemos hacer lo siguiente: r = [10 11 12]; A = [A; r] y resultaría A= 1 2 3 456 789 10 11 12
Si omites el nombre de la variable y el signo "=", MATLAB automáticamente crea la variable a n s para guardar el resultado. También distingue las letras mayúsculas de las minúsculas. Todos los nombres de funciones deben ser en letras minúsculas.
Los ejemplos que hemos dado se han guardado en variables que están en el espacio de trabajo de MATLAB. Para listar las variables en el espacio de trabajo se utiliza el comando who . Para ver información adicional acerca de estas variables se utiliza el comando whos.
Las variables permanentes son aquellas con significado especial, y que no se pueden eliminar. Estas son por ejemplo las variables a n s y e p s. La variable e p s es una tolerancia para determinar. Por ejemplo la singularidad y el rango. Su valor inicial es la distancia de 1.0 al próximo número de punto flotante mayor.
Las funciones que utiliza MATLAB son intrínsecas al procesador de éste. Otras funciones están disponibles en la librería externa de archivos -M. Además de éstas funciones todo usuario también puede crear otras funciones. Puedes combinar las funciones de acuerdo a tu necesidad. Ejemplo: x = sqrt(log(z) )
Para salir de MATLAB se escribe quit ó exit . Al terminar una sesión de MATLAB, las variables en el espacio de trabajo se borran. Si deseas guardar tu espacio de trabajo escribes save . save guarda todas las variables en un archivo llamado matlab.mat. Se puede utilizar save y load con otros nombres de archivos, ó para guardar solo variables seleccionadas Ejemplo: save temp X Y Z Este ejemplo guarda las variables X, Y, Z en el archivo temp.mat. Usando el comando load temp las obtienes nuevamente del archivo temp.mat. load y save también pueden importar y exportar información de archivos ASCII.
Índices Podemos referirnos a elementos individuales de matrices encerrando sus índices en paréntesis. Ejemplo: A= 123 456 789
A(3, 3) = A(1, 3) + A(3, 1) resultaría A= 123 456 7 8 10
Un índice puede ser un vector. Si x y v son vectores, entonces x(v) es [x(v(1)), x(v(2)), ...,x(v( n))]. Para matrices, los índices de vectores permiten acceso a submatrices contiguas y no -contiguas.
Ing Sergio Casillas de la Torre Tutorial de Matlab Por ejemplo, suponga que A es una matriz 10 por 10. Entonces A(1:5, 3) especifica la submatriz 5 x 1, ó vector columna, que consiste de los primeros cinco elementos en la tercera columna de A. También A(1:5, 7:10) es la submatriz 5 x 4 de las primeras cinco filas y las últimas cuatro columnas. Utilizando solo los dos puntos denota todo lo correspondiente a la fila ó columna. Podríamos tener una instrucci ón como: A(:, [3 5 10]) = B(:, 1:3) que reemplaza la tercera, quinta y décima columna de A con las primeras tres columnas de B.
Manipulación de Matrices d i a g - extrae ó crea una diagonal t r i l - parte inferior triangular triu - parte superio r triangular ' - transposición
Las operaciones suma (+) y resta (-) son definidas para las matrices siempre y cuando éstas tengan la misma dimensión. Es decir, si A y B son matrices 3 x 3, entonces A + B se puede calcular. Las operaciones suma y resta también está definidas si uno de los operandos es un escalar, es decir, una matriz 1 x 1. Ejemplo: x= -1 0 2
y=x - 1 resultaría y = -2 -1 1
Ejemplo: >>A=[1 2 3;4 5 6]; B=[6 5 4; 3 2 1]; define las matrices A y B. Para sumarlas se escribe la operación: >>A+B El resultado de la operación es por defecto almacenado en la variable ans e inmediatamente presentado en pantalla: ans = 777 777 Para almacenar la suma de A y B en la variable C: >>C=A+B C= 777 777
El producto de una matriz y un vector es un caso especial del producto matrizmatriz y naturalmente, un escalar como pi, puede multiplicar, ó ser mu ltiplicado por, cualquier matriz.
En división de matrices, si A es una matriz cuadrada no-singular, entonces A\B y B/A corresponden a la multiplicación izquierda y derecha de B por el inverso de A, esto es, inv(A) * B y B * inv(A) res pectivamente. El resultado es obtenido directamente sin la computación del inverso.
A \B es definido cuando B tiene la misma cantidad de filas que A. Si A es cuadrada, el método usa do es la Eliminación Gaussiana. El resultado es una matriz X con las mismas dimensiones que B. Si A no es cuadrada, se factoriza utilizando la ortogonalización 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 número de columnas de A y n es el número 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 términos de A\B por B/A = (A'\B ' ) ' .
La expresión A^n eleva A a la n- ésima potencia y esta definido si A es una matriz cuadrada y n un escalar.
M ATLAB considera expresiones como exp(A) y sqrt(A) como operaciones de arreglos, definidas en los elementos individuales de A. También puede calcular funciones trascendentales de matrices, como la matriz exponencial y la matriz
Ing Sergio Casillas de la Torre Tutorial de Matlab logarítmica. Estas operaciones matrices cuadradas. especiales están definidas solamente para
poly - polinomio característico d e t - determinante trace - traza kron - producto tensorial de Kronecker eig - calcula los valores propios de la ma triz
El término operaciones de arreglo se refiere a las operaciones de aritmética elemento por elemento. Un punto (.) antes de un operador indica una operación de arreglos elemento por elemento.
El símbolo .* denota multiplicación de arreglos elemento por elemento. Ejemplo: x = [1 2 3]; y = [4 5 6]; z = x. *y resulta z= 4 10 18 Las expresiones A./B y A. \B dan los cocientes de los elementos individuales. Ejemplo: z = x.\ y resulta z= 4.0000 2.5000 2.0000
>> a' ans = 212
>> b' ans = 123 29
>> a.*b ans = 2 2 6 >> a*b' ans = 2 46 123 246
Ing Sergio Casillas de la Torre Tutorial de Matlab >> a/3 ans = 0.6667 0.3333 0.6667
>> a./3 ans = 0.6667 0.3333 0.6667 >> a^b ??? Error using ==> ^ Matrix dimensions must agree.
Ing Sergio Casillas de la Torre Tutorial de Matlab >> 2.^a ans = 4 2 4
Aproximadamente 16 dígitos significativos computadoras utilizando aritmética flotante IEEE. El rango aproximado es: 10^-308 a 10^308.
4/3 a) format short 1.3333 b) format short e 1.3333e+00 c) format long 1.33333333333333 d) format long e 1.33333333333333e00 e) format bank 1.33 f) format hex 3ff 5555555555555
Programar en MatLab es usar una serie de comandos que permitan realizar una tarea o función específica. Estos pueden ser escritos uno por uno a través de la línea de comandos: >>A=[1 2 3;4 5 6;7 8 9] A= 123 456 789 >>A' ans = 147 258 369 El primer comando A=[1 2 3;4 5 6;7 8 9] define la matriz A y el siguiente comando A' calcula y presenta en pantalla la transpuesta de A.
Los archivos de disco que contienen instrucciones de MATLAB se llaman archivos- M. Esto es así porque siempre tienen una extención de ".m" como la última parte de su nombre de archivo. Un archivo - M consiste de una secuencia de instrucciones normales de MATLAB, que probablemente incluyen referencias a otros archivos - M. Un archivo- M se puede llamar a sí mismo recursivamente. Puedes crear archivos-M utilizando un editor de texto ó procesador de palabras. Hay dos tipos de archivos - M: los de comandos y las funciones. Los archivos de co mandos, automatizan secuencias largas de comandos. Los archivos de funciones, permiten añadir a MATLAB funciones adicionales expandiendo asi la capacidad de este programa. Ambos, comandos y funciones, son archivosordinarios de texto ASCII.
Archivos de Comandos Cuando un archivo de comandos es invocado, MATLAB simplemente ejecuta los comandos encontrados en dicho archivo. Las instrucciones en un archivo de comando operan globalmente en los datos en el espacio de trabajo. Los comandos son utilizados para hac er análisis, resolver problemas, ó diseñar secuencias
Ing Sergio Casillas de la Torre Tutorial de Matlab largas de comandos que se conviertan en interactivas. Por ejemplo, suponga que el archivo fibo.m contiene los siguientes comandos de MATLAB:
% Un archivo-M para calcular los elementos de la serie de Fibonacci f = [1 1]; i = 1; while f(i) + f(i+1) < 1000 f(i+2) = f(i) + f(i+1); i = i + 1; end plot(f)
Si escribimos fibo en una ventana de MATLAB seguido de "enter" vemos que MATLAB calcula los primeros 16 números de Fibonacci, y luego grafica estos. Luego que la ejecución del archivo es completada, las variables f y i permanecen en el espacio de trabajo. Los programas de demostraciones incluidos en MATLAB son ejemplos de como usar comandos para hacer tareas más complicadas. Para utilizar estos escriba demo s en el "prompt" de MATLAB.
Archivos de Funciones Un archivo - M que contiene la palabra function al principio de la primera línea, es un archivo de función. En una función, a diferencia de un comando, se deben de pasar los argumentos. Las variables definid as y manipuladas dentro de la función son locales a esta y no operan globalmente en el espacio de trabajo. Los archivos de funciones se utilizan para extender a MATLAB, i.e., crear nuevas funciones para MATLAB utilizando el lenguaje propio de MATLAB. E l a rchivo mean.m contiene las instrucciones: function y = mean(x) % Valor medio. % Para vectores, mean(x) retorna el valor medio de los elementos del vector x. % Para matrices, mean(x) es un vector fila conteniendo el valor medio de cada columna. [m, n] = size (x); if m == 1 m = n; end y = sum(x)/m;
Ing Sergio Casillas de la Torre Tutorial de Matlab (Las lineas que comienzan con "%" son interpretadas como comentarios por MATLAB). La existencia de este archivo en el disco duro define una nueva función en MATLAB llamada mean. Si z es un vector de los enteros desd e 1 a 99, por ejemplo, z = 1:99; entonces, el valor promedio es encontrado escribiendo mean(z) que resultaría ans = 50
Veamos algunos detalles de mean.m: La primera línea declara el nombre de la función, los argumentos de entrada, y los argumentos de salida. Sin esta línea sería un archivo de comando. % indica que el resto de la línea es un comentario. Las primeras líneas documentan el archivo - M y aparecen en la pantalla cuando escribimos help mean. Las variables m, n, e y son locales a mean y n o existen en el espacio de trabajo. (O si existen, permanecen sin cambios.) No es necesario asignar los enteros de 1 al 99 en la variable x. Utilizamos mean con una variable llamada z. Este vector que contenía los enteros de 1 a 99 fue pasado ó copiado a mean donde se convirtió en una variable local llamada x. Ejemplo % Ejemplo de un archivo- m % Creación del vector x usando el comando for n=5; for i=1:n x(i)=i^2; end x % Fin del archivo-m Este ejemplo es un archivo -m tipo comando. Para ejecutarlo, en la línea de comandos se debe escribir el nombre del archivo: >>ejemplo x = 1 4 9 16 25
Ing Sergio Casillas de la Torre Tutorial de Matlab Ejemplo % Calcula el promedio de los elementos de un vector y dibuja dicho vector % Sintaxis : promedio(x) donde x es el vector a promediar function p = promedio(x) n=length(x); p=0; for i=1:n p=p+x(i); end p=p/n; plot(x);
Para ejecutar la función, se hace la llamada en la línea de comandos incluyendo el parámetro. La función promedio usa por parámetro un vector. Este vector debe ser definido previamente. >>A=[1 2 4 3 7 5 6 1 2 0 8 5]; >>promedio(A) ans = 3.6667
Ing Sergio Casillas de la Torre Tutorial de Matlab Esta imagen es el res ultado del comando plot(x) al ejecutar la función promedio. MatLab posee un conjunto de archivos- m incorporados (built-in). Puede agregársele archivos - m definidos por el usuario almacenando los mismos en el directorio principal de MatLab. Los comentarios incluidos en estos scripts y funciones se visualizan al usar el comando help seguido del nombre del archivo. >>help promedio Calcula el promedio de los elementos de un vector y dibuja dicho vector Sintaxis : promedio(x) donde x es el vector a promediar
Para ver el contenido de un archivo- m se usa el comando type seguido del nombre del archivo.
Funciones Matemáticas Algunas funciones trigonométricas utilizadas por MATLAB son: sin - seno c o s - coseno t a n - tangente asin - seno inverso ac o s - coseno inverso atan - tangente inversa
Algunas funciones elementales son: real(a) Parte real imag(a) Parte imaginaria conj(a) Conjugado de a fft(x) Transformada discreta de Fourier del vector x fft(x,n) FFT de n puntos muestrales ifft(x) Transfor mada inversa rápida de Fourier del vector x ifft(x,n) FFT inversa de n puntos muestrados zeros Inicializa a ceros zeros(n) Matriz de nxn de ceros zeros(m,n) Matriz de mxn de ceros y=zeros(size(A) Matriz del tamaño de A, todos ceros
Ejemplo size Regresa el número de filas y columnas A= 0 7 -6 100 010
>> [m n]=size(A) m = 3 n= 3 Funciones matriciales tril(A) Matriz triangular inferior triu(A) Matriz triangular superior pascal Triangulo de Pascal tocplitz Tocplitz
Ejemplos >> A = 0 7 -6 100 010 >> t o e p l i t z ( A ) ans = 0 1 0 7 0 1 -6 0 0 1 0 1 0 7 0 1 -6 0 0 1 0 1 0 7 0 1 -6 Producto de dos matrices triangulares. Esta factorización se utiliza para obtener el inverso y el determinante. También es la base para la solución de sistemas lineales. Para ob tener la factorización LU de A escribimos, [L, U] = lu(A).
Factorización Ortogonal ó Factorización QR. Se utiliza para matrices cuadradas ó rectangulares. Esta factorización se utiliza para resolver sistemas lineales con más ecuaciones que desconocidas. Esta factorización también es la base para las funciones null y orth , que generan bases ortonormales para el espacio nulo y rango de una matriz rectangular dada.
Descomposición de Valores Singulares La descomposición de Valores Singulares es importante problemas que envuelvan matrices. para el análisis de
La asignación triple [U, S, V] = svd(A) produce los tres factores en la descomposición de valores singulares A = U*S*V'. Las matrices U y V son ortogonales y la matriz S es diagonal. La función s v d ( A ) devuelve solamente los elementos de la diagonal de S, que son los valores singulares de A.
Descomposición de Valores Propios La Descomposición de Valores Propios se utiliza para obtener los valores y vectores propios de una matriz cuadrada A. La función e i g ( A ) devuelve los valores propios de A en un vector columna. La asignación [X,D]=eig(A) produce una matriz diagonal D cuyos elementos diagonales son los valores propios de A y las columnas de X son los vectores propios correspondientes.
Las Funciones de norma, rango y acondicionamiento asociadas son: cond - número de condición en la norma 2 norm - norma 1, norma 2, norma F, norma rank - rango rcond - estimado del número de condición
Funciones de Funciones MATLAB representa funciones matemáticas mediante archivos- M función. Un ejemplo de una función es el archivo -M llamado h u m p s . m. de tipo
Ejemplo: El archivo- M llamado h u m p s . m contiene las siguientes instrucciones: function y = humps(x)
Ing Sergio Casillas de la Torre Tutorial de Matlab y = 1./((x- .3).^2 +.01) + 1./((x- .9).^2 +.04) - 6; % y para la gráfica de la función escribimos x = - 1:.01:2; plot(x, humps(x))
Integración Numérica (Cuadratura) El área bajo la gráfica de la función f(x) se puede aproximar integrando f(x) numéricamente mediante una regla de cuadratura. Para integrar la función definida por h u m p s . m desde 0 hasta 1 escribimos: q = quad('humps', 0, 1) q= 29.8583
Note que el argumento de quad contiene un nombre de una función. Por esto quad se llama una función de función, i.e., es una función que opera en otras funciones.
Ecuaciones No-lineales y Funciones de Optimización Las funciones de funciones para ecuaciones no -lineales y optimización incluyen: fmin - mínimo de una función de una variable fmins - mínimo restricciones) de una función multi - variable (minimización no-lineal sin
fzero - cero de una función de una variable c o n s t r - minimización con restricciones fsolve - solución de ecuación no-lineal leastsq - cuadrados mínimos no-lineales
Funciones para Ecuaciones Diferenciales Las funciones de MATLAB para resolver problemas de valo r inicial para ecuaciones diferenciales ordinarias son: ode23 - método Runge- Kutta de largo de paso variable que combina un método de orden dos con uno de orden tres. ode45 - método Runge-Kutta -Fehlberg de largo de paso variable que combina un método de orden cuatro con uno de orden cinco.
Ing Sergio Casillas de la Torre Tutorial de Matlab Ejemplo to=0; tf=10; [t,x]=ode23(`edif',to,tf,xo);
[t,x]=ode23(`deriv',to,tf,xo); ode45
[t,x]=ode23(`deriv',to,tf,xo,to1,trace); ode45 trace => 0 - no resuntados intermedios 1 - resultados intermedios
default tol: o de23 -> 1.0e - 0 3 ode45 - > 1.oe- 0 6
Sintaxis: function nombre_1=nombre_2(parametro_1, ..., parametro_n) Ejemplos: function y=promedio(x) function i=inodal(t,v) function xpunto=vdpol(t,x) xpunto=zeros(2,1); xpunto(1)=x(1).*(1 -x ( 2 ) . ^2 )-x(2); xpunto(2)=x(1);
Ing Sergio Casillas de la Torre Tutorial de Matlab Ejemplo: if n< maxn ... if n>=0, break, end
Los operadores &, | y ~ son los operadores de lógica "y", "ó" y "no" respectivamente. El resultado de C = A & B es una matriz cuyos elementos son unos donde A y B sean ambos distintos de cero, y ceros donde A ó B sean cero. A y B deben de ser matrices con las mismas dimensiones, a menos que una de ellas sea un escalar. El resultado de C = A | B es una matriz cuyos elementos son unos donde A ó B tienen un elemento diferentede cero, y ceros donde ambas tienen elementos cero. A y B deben de se r matrices con las mismas dimensiones, a menos que una sea un escalar. El resultado de B = ~A es una matriz cuyos elementos son uno donde A tiene un elemento cero, y ceros donde A tiene elementos diferentes de cero.
Funciónes any, all La función any(x) de vuelve 1 si cualquiera de los elementos de x es diferente de cero, de lo contrario devuelve 0. La función all(x) devuelve 1 solamente si todos los elementos de x son diferentes de cero. Estas funciones se usan en cláusulas i f. Por ejemplo: if all(A <.5) . . . end Para argumentos matriciales, a n y y a l l trabajan por columnas para devolver un vector fila con el resultado para cada columna. Aplicando la función dos veces, any(any(A)), siempre reduce la matriz a una condición escalar. Las funciones relacionales y lógicas en MATLAB son: a n y - condiciones lógicas a l l - condiciones lógicas find - halla índices de arreglos de valores lógicos exist - verifica si existen variables isinf - detecta infinitos
Ing Sergio Casillas de la Torre Tutorial de Matlab finite - verifica para los valores finitos
Los caracteres especiales de MatLab son: [ ] Se utilizan para formar vectores y matrices ( ) Define precedencia en expresiones aritméticas. Encierra argumentos de funciones en forma usual , Separador de elementos de una matriz, argumentos declaraciones en líneas con declaraciones múltiples ; Termina filas de una matriz, separador de declaraciones % Comentario Ejemplos: [6.0 9.0 3.4 ] sqrt(2) for i=1:n, a(i)=0, end for i=1:n; a(i)=0; end % inicia vector a en 0 de funciones y
Sintaxis for variable=incio:paso:final declaración 1; ... declaración n; end for variable=inicio:final declaración 1; ... declaración n; end
El ciclo FOR permite que una instrucción, ó grupo de instrucciones, pueda repetirse un número determinado de veces. Por ejemplo, for i = 1:n, x(i) = 0, end asigna 0 a los primeros n elementos de x. Si n es menor de 1, el ciclo sigue siendo válido pero MATLAB no ejecuta la instrucción intermedia. Si x no esta definido, ó si tiene menos de n elementos, entonces un espacio adicional es localizado automáticamente a x cada vez que sea necesario.
5.5.2 Declaración FOR anidada.
Sintaxis for variable 1 = inicio1:paso1:fin1 for variable2 = inicio2:paso2:fin2 declaración 1; ... declaración n; end end
Ejemplo y=1 for t1=0:0.1:1 for t2=1: - 0.1:0 y(1)=sin(t1*t2) end i=i+1; end
Ejemplo for i = 1:m for j = 1:n A(i, j) = 1/(i+j- 1); end end A La "A" al terminar el ciclo muestra en la pantalla el resultado final. Es importante que para cada for halla un end.
Sintaxis: while expresion proposición 1; ... proposición 2; end Ejemplos e=1.0; while (1.0+e)>1.0001 e=e/2.0; end
it=1; t=0; wo =2.0*pi*60.0; while it<=npts, ut=sin(wo*t);t=t+dt;end
El ciclo WHILE permite a una instrucción, ó grupo de instrucciones, repetirse un número indefinido de veces, bajo el control de una condición lógica. El siguiente ciclo while halla el primer entero n para el cual n! es un número de 100 digitos:
Un cálculo más práctico ilustrando el ciclo while es en el cómputo del exponencial de una matriz, llamado expm(A) en MATLAB. Una posible definición de la función exponencial es mediante la serie:
La idea es sumar todos los términos necesarios hasta producir un resultado que, en la precisión finita la de computadora, no cambie aunque más términos sean añadidos. Para esto pr ocedemos de la forma siguiente:
Ing Sergio Casillas de la Torre Tutorial de Matlab E = zeros(size(A)); F = eye(size(A)); k = 1; while norm(E+F- E, 1) > 0 E = E + F; F = A*F/k k = k+1; end Aqui A es la matriz dada, E representa la suma parcial de la serie, F es un término individual en la serie, y k es el índice de este término.
Sintaxis a) if expresión proposición 1; ... proposición n; end
Ing Sergio Casillas de la Torre Tutorial de Matlab c) if expresión p r oposición 1; ... proposición n; elseif proposición 1; ... proposición m; else proposición 1; ... proposición r; end
d) if expresión, break, end Ejemplos if dv(i) > maxer maxer=dv(i); nmaxe=i; end
sum=0.0; y=1; while i<=so n=input(`Introduzca n, interrump e con valor negativo `); if n<0, break, end; if n==0 sum=sum+n; elseif n<=10 sum=sum+n/2; else sum=sum+n/10; end end
Ing Sergio Casillas de la Torre Tutorial de Matlab A continuación se muestra como un cálculo se puede dividir en tres casos, dependiendo del signo ó paridad de un entero n:
En el segundo, partiendo de un entero positivo n, si este es par, se divide entre dos; si es impar, se multiplica por tres y se le suma uno. ¿Habrá algún entero para el cual el proceso nunca termine? Aquí se ilustran los enunciados while y i f , también se muestra la función input (en este caso es una entrada del teclado), y el enunciado break, que provee salidas abruptas de los ciclos. Veamos:
% Problema "3n+1" clásico de la teoria de número s . while 1 n = input('Entre n, negativo termina. '); if n <= 0, break, end while n > 1 if rem(n, 2) == 0 n = n/2 else n = 3*n+1 end end end
Ejemplo >> A=[1 2 3; 4 5 6; 7 8 9] A= 123 456 789
Sintaxis: matriz_modificada = reshape(matriz_original, filas, columnas) Ejemplo >> A=[1 4 7 10; 2 5 8 11; 3 6 9 12] A= 1 4 7 10 2 5 8 11 3 6 9 12
5.6.3 Modific ación individual de elementos
>> A(1,1)=A(1,2)+A(2,1) A= 5 2 3 4
Ing Sergio Casillas de la Torre Tutorial de Matlab Modificación de los elementos >> A(:)=10:15 A= 10 13 11 14 12 15
>> x=5:-1 : 1 x = 5 4 3 2 1
Acceso a submatrices contiguas y no contiguas Ejemplos Si la matriz original A es de 10*10, entonces: A(1:3,5) de A matriz de 3x1 que tiene los tres primeros elementos de la columna 5
Ing Sergio Casillas de la Torre Tutorial de Matlab Matrices vacias La declaración x = [ ] asigna una matriz de dimensión 0x0 a x
Generación de tablas >> x=(0.0:0.2:3.0); >> y=exp(- x).*sin(x); >> [x;y] ans = Columns 1 through 7 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.200 0 0 0.1627 0.2610 0.3099 0.3223 0.3096 0.2807 Columns 8 through 14 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 0.2430 0.2018 0.1610 0.1231 0.0896 0.0613 0.0383 Columns 15 through 16 2.8000 3.0000 0.0204 0.0070
Ing Sergio Casillas de la Torre Tutorial de Matlab Determinante de A: det(A) >> A=[1 2 3;4 5 6; 7 8 9] A= 123 456 789 >> det(A) ans =0
Diagonal de A: diag(A) >> diag(A) ans = 1 5 9 Valores y vectores característicos : e i g ( A ) >> A=[0 7 - 6; 1 0 0;0 1 0] A= 0 7 -6 100 010
>> eig(A) ans = -3.0000 2.0000 1.0000
v - Vectores caracterís ticos d - valores característicos
Ing Sergio Casillas de la Torre Tutorial de Matlab >> [v d]=eig(A) v= 0.9435 -0.8729 0.5774 -0.3145 -0.4364 0.5774 0.1048 -0.2182 0.5774 d= -3.0000 0 0 0 2.0000 0 0 0 1.0000
- Exponencial de una matriz : expm(A) >> A=[0 7 - 6; 1 0 0;0 1 0] A= 0 7 -6 100 010
>> e x p m ( A ) ans = 5.2541 11.0757 - 13.6115 2.2686 5.2541 - 4.8044 0.8007 2.2686 - 0.3510
- Factorización LU de A: lu(A) >> [L U]=lu(A) L= 0 1.0000 0 1.0000 0 0 0 0.1429 1.0000 U = 1.0000 0 0 0 7.0000 - 6.0000 0 0 0.8571
- Inversa de A: inv(A) >> inv(A) ans = 0 1.0000 0 0 0 1.0000 -0.1667 0 1.1667 - Ecuación característica de la matriz A: poly(A) >> p=poly(A) p= 1.0000 0.0000 - 7.0000 6.0000 - Raices de la ecuación característica : roots(p) >> r=roots(p) r = -3.0000 2.0000 1.0000
Sintaxis id = fopen(`nombre.dat', `permiso') donde permiso puede ser: `r' Abre archivo para lectura `r+ Abre archivo para lectura y escritura `w' Borra el contenido del archivo existente o crea un nuevo archivo y lo abre para escritura `w+' Idem que `w' únicamente que el archivo se abre para lectura y escritura `a' Crea y abre un nuevo archivo o abre un archivo `a+' Idem que `a' únicamente que el archivo es abierto para lectura y escritura
Ejemplo fid = fopen(`archiv o.dat','r') fid = - 1 , e r r o r 0, lectura/escritura normal [fid, mensaje = fopen(`archivo.dat','r')
Ing Sergio Casillas de la Torre Tutorial de Matlab Ejemplo: A = fread(fid,10,'float')
Salida con formato Ejemplos: fprintf(fid,'titulo \ n'); fprintf(fid,'%f %12.7f\ n', y); Formato %s - cadena decimal %d - número decimal %f - punto flotante % g - formato g
Son variables, de las cuales una sola copia es compartida por el programa principal y sus funciones. Sintaxis: global variable1, ..., variable_N Ejemplo function x=ccdifs(t,x) global ka,kb xp=[x(1)- ka*x(1)*x(2); - x(2)+kb*x(1)*x(2)]; ... global ka,kb ka=0.01 kb=0.02 [t,x]=ode23('ccdifs',0,10,[1:1]);
Para que los programas en MATLAB ejecuten más rápido, debemos vectorizar estos siempre que sea posible. Esto es, debemos convertir los ciclos for y while a operaciones de vectores ó de matrices. Por ejemplo, un modo de calcular la función "sin" para 1001 números entre 1 y 10 es: i = 0; for t = 0:.01:10 i = i + 1; y(i) = sin(t); end Una versión vectorizada del mismo código es t = 0:.01:10; y = sin(t); En una computadora lenta, el primer ejemplo tomó 15 segundos, mientras que el segundo tomó 0.6 segundos.
Vectores Pre- Asignados Si no podemos vectorizar un pedazo de código, podemos hacer que los ciclos for vayan más rápido pre-asignando cualquier vector en el cual el resultado de salida sea guardado. Veamos un ejemplo: y = zeros (1,100); for i = 1:100 y(i) = det(X^i); end Si no pre -asignamos el vector "y", el interpretador de MATLAB irá aumentando el tamaño de "y" por uno cada vez que se itera en el ciclo. Permite incrementar la velocidad de proceso de MATLAB Sintaxi s variable=inicio:incremento:final Ejemplo i=1, wo=2*pi*fo; for t=0:dt:per f(i)=sin(wo*t); i=i+1; end ó t=0:dt:per; fi=sin(wo*t);
plot - crea una gráfica de vectores ó columnas de matrices. loglog - crea una gráfica utilizando una escala logarítmica para ambos ejes. semilogx - crea una gráfica utilizando una escala logarítmica para el eje-x y una escala lineal para el eje-y . semilogy - crea una gráfica utilizando una escala log arítmica para el eje -y y una escala lineal para el eje-x . Puedes añadir títulos, encabezamientos de ejes, líneas entre cortadas y texto a tus gráficas utilizando: tittle - añade título a la gráfica xlabel - añade encabezamiento al eje-x ylabel - añade encabezamiento al eje-y text - añade una cadena de texto en una localización específica gtext - añade texto a la gráfica utilizando el ratón grid - crea líneas entrecortadas
Sintaxis: a) plot(y) b) plot(x,y) c) plot(x,y,'tipo_línea') d) plot(x1,y1,'tipo_línea_1',x2,y2,'tipo_línea_2', ... , xn,yn,'tipo_línea_n') Si y es un vector, p l o t ( y ) produce una gráfica lineal de los elementos de y versus el índice de estos. Si especifica dos vectores como argumentos, plot(x, y) produce una gráfica de y versus x.
Ing Sergio Casillas de la Torre Tutorial de Matlab b azul w blanco k negro Símbolo Estilo de línea . punto o circulo x marca + mas * asterisco - sólido : punteado -. segmento punto -- s e g m e n t o
plot(Y) - dibuja una línea para cada columna de Y. El eje-x es encabezado por el vector índice de fila, 1:m, donde m es el número de filas en Y. Si plot es usado con dos argumentos y si X ó Y tienen más de una fila ó columna, entonces: si Y es una matriz, y x es un vector, plot(x,Y) grafica las filas ó columnas de Y versus el vecto r x ; si X es una matriz y y es un vector, plot(X,y) grafica cada fila ó columna de X versus el vector y; si X y Y son ambas matrices del mismo tamaño, plot(X, Y) grafica las columnas de X versus las columnas de Y.
Puede importar y graficar datos generados fuera de MATLAB utilizando el comando load .
Hay diferentes formas de graficar funciones y = f(x). Una de estas formas es evaluar la función en miles de puntos en el intervalo de interés. La siguiente función oscila infinitamente rápido en el intervalo, 0 x 1. Podemos gráficarla como sigue: x = (0:1/2000:1)'; plot(x, cos(tan(pi*x))) Para hacer esto más eficiente podemos usar la función fplot la cual concentra su evaluación sobre las regiones donde la rapidez de cambio de la función es más grande.
Para evaluar una función, se crea un archivo de esta función y se le pasa el nombre del archivo a fplot. El siguiente archivo -M de tipo función define la función anterior como fofx. function y = fofx(x) y = cos(tan(pi*x)); % Este archivo se guarda con el nombre de fofx.m. Ahora la instrucción fplot('fofx', [0 1]) produce la gráfica:
hold Permite añadir líneas al dibujo previo o n Activa hold off Desactiva hold Ejemplo plot(x); hold on; plot(y',':');plot(yz,'-.')
loglog Sintaxis a) loglog(x,y) b) loglog(x,y,'tipo_línea') c)loglog(x1',y1','tipo_línea_1',..., xn,yn,'tipo_línea_n') Ejemplo x=logspace(- 1,3); loglog(x,exp(x)) donde logspace tiene las formas: logspace(a,b) logspace(a,b,n) a,b exponentes de los límites. Es decir, 10^a y 10^b
Ing Sergio Casillas de la Torre Tutorial de Matlab fill Dibuja el area interior de una curva en determinado color Sintaxis: a) fill(x,y,'c') b) fill(x1,y1,'c1',...,xn,yn,cn) Ejemplo t=0:0.5:2*pi; x=sin(t); fill(t,x,'b') t=0:0.5:2*pi; x=sin(t); y=cos(t); fill(y,x,'r')
subplot Di buja la pantalla en mxn subdivisioens, numeradas por el parámetro p, de izquierda a derecha, iniciando por la fila superior Sintaxis: subplot(m,n,p) Ejemplo: vt=smvars(:,1); it=smvars(:,2); rang=smvars(:,3); ikd=smvars(:,4); subplot(2,2,1); plot(vt) subplo t(2,2,2) plot(it) suplot(2,2,3) plot(rang) subplot(2,2,4) plot(ikd)
bar Crea una gráfica de barras Sintaxis: a) bar(y); b) bar(x,y); c) [xb,yb]=bar(y); => plot(xb,yb) d) [xb,yb]=bar(x,y); => plot(xb,yb)
Ing Sergio Casillas de la Torre Tutorial de Matlab Ejemplo x= - 2.8:0.2:2.8 bar(x,exp(- x.*x) Nota: Los valores de x deben estar igualmente espaciados
stairs Igual que b a r, únicamente sin líneas internas
fplot Dibuja la gráfica de una función Sintaxis: a) fplot(`función', [inicio,final]) b) fplot(`función', [inicio,final],n) c) fplot(`función', [inicio,final],n,ángulo) d) [x,y]=fplot(`función', [inicio,final]) => plot(x,y) n - número de puntos ángulo - ángulo entre segmentos sucesivos de la función Ejemplo fplot(`sin',[0,pi]) fplot(`tanh',[ -2 2]) function y=func(x) y=200*sin(x(:))./x(:); fplot(`func',[-3 0 30],60,2)
colormap Colorea con sombreado el interior de una curva o polígono Sintaxis colormap(colorbase) donde colorbase es: gray hot cool copper pink Ejemplo t=0:0.05:2*pi; x=sin(t); y=cos(t); colormap(hot(130)); ... % Nota: 130 es opcional el rango 0-255 fill(y,x,x) => sombreado horizontal fill(y,x,y) => sombreado vertical
contour, contour3 Genera dibujos compuestos de líneas de valores de datos constantes obtenidos de una matriz de entrada Sintaxis: a) contour(z) b) contour(z,n)
Ing Sergio Casillas de la Torre Tutorial de Matlab c) contour(x,y,z) d) contour(x,y,z,n)
Ejemplo contour3(peaks,30)
Ejemplo. Evalue y dibuje la funcion z=x*exp( -x^2 -y^2) sobre el rango -2<=x<=2, -2 < = y < = 2 [X,Y]=meshgrid(-2:2:2); z=x.*exp(- x ^ 2-y^2); ... mesh(Z)
meshz Añade un plano de referencia o cortina al dibujo
surf, surfc Crean superficies sombreadas en 3 dimensiones Sintaxis: a) surf(x,y,z,c) b) surf(x,y,z) c) surf(x,y,z,c) d) surf(x,y,z) e) surf(z,c) f) surf(z) g) surfc(...) => misma Sintaxis que surf
Ejemplo [x,y]=meshgrid(- 3:.2:3); z=peaks(x,y); surf(x,y,z) k=5; n=2^k-1; [x,y,z]=sphere(n); c=hadamard(2^k);... surf(x,y,z,c); colormap(hot)
hadamard Matriz hadamard compuesta de 1's y -1's, empleada en procesamiento de señales y análisis numérico Ejemplo (matriz de 4*4) 1 1 1 1 1 -1 1 -1 1 1 -1 -1 1 -1 -1 1
surfl Superficie sombreada tridimensioanl con efecto de reflexión de luz Sintaxis: a) surfl(z) b) surfl(z,s ) c) surfl(x,y,z) d) surfl(x,y,z,s) s - dirección de la luz
shading flat - cada segmento de la superficie tiene un valor constante determinado por el color de los puntos extremos del segmento o sus esquinas shading interp - el color en cada segmento varia linealmente e interpolo los valores extremos o esquinas shading faceted superpuestas utiliza sombreado "flat" con líneas de malla negras
axis Escala y apariencia de los ejes Sintaxis: a) axis([xmin, xmax, ymin, ymax]) b) axis([xmin, xmax, ymin, ymax, zmin, zmax]) c) axis(`auto') d) axis(`ij') e) axis(`xy') f) axis(`square') g) axis(`equal') h) axis(`off') 70
Ing Sergio Casillas de la Torre Tutorial de Matlab i) axis(`on') donde: axis(`auto') realiza el escalamiento de ejes a su modo de autoescalamiento por defecto. axis(`ij') dibuja nuevamente la gráfica. El eje y es vertical y es numerado de arriba hacia abajo. El eje j es horizontal y es numerado de izquierda a derecha. axis(`xy') regresa la forma de ejes cartesianos que existe por defecto . El eje x es horizontal y se numera de izquierda a derec ha. El eje y es vertical y se numera de abajo hacia arriba axis(`square') determina que la región de los ejes es cuadrada axis(`equal') indica que los factores de escalamiento y marcas incrementales a lo largo de los ejes x y y son iguales. axis(`off') des activa las etiquetas de los ejes y las marcas axis(`on') activa las etiquetas de los ejes y las marcas Para el ejemplo último: ... axis([ - 3 3 - 3 3 - 8 8])
rand matrices y números aleatorios distribuidos uniformemente Sintaxis: a) rand(n) - matriz de nxn b) rnad(m,n) - matriz de mxn
Ing Sergio Casillas de la Torre Tutorial de Matlab load carga en el area de trabajo un archivo (imagen, sonido, datos, etc) Sintaxis a) load archivo b) load archivo.ext donde: ext - extensión
image crea un objeto imagen y lo presenta Sintaxis: a) image(x) b) image(x,y,x) c) presenta la matriz c como una imagen d) especifica los límites de los datos de la imagen en los ejes x e y. En b), x e y son vectores Ejemplo load clown colormap(map) image(x)
brighten hace más brillante o más obscura la imagen Sintaxis: a) brighten(alfa) b) brighten(map, a l f a ) donde: 0<alfa<1 más brillante -1<alfa<0 más obscuro Del ejemplo anterior: ... brighten(0.6) ó brighten(- 0.6)
sound convierte un vector en sonido (en computadoras sparc y macintosh) Sintaxis a) sound(y) b) sound(y,Fs) donde: Fs frecuencia especificada en Hz
Ejemplo load train sound(y,Fs) t=(0:length(y)- 1)/Fs; plot(t,y)
Algunos comandos utilizados para la manipulación de archivos de disco son dir, type, delete y c d . Si la extención no se especifica, MATLAB utiliza .m automáticamente. El comando diary crea un diario de tu sesión de MATLAB en un archivo de disco. Para más información utiliza la Guía de Referencia de MATLAB ó el comando help.
5.12.2 Ejecutando Prog ramas Externos
El simbolo "!" le indica a MATLAB que el resto de la línea de entrada es un comando para el sistema operativo. Por ejemplo, ! edt darwin.m invoca un editor llamado e d t en un archivo llamado darwin.m. Luego que este programa sea completado, el sistema operativo devuelve el control a MATLAB.
Puedes introducir datos de otros programas a MATLAB por varios métodos. Similarmente, puedes exportar datos de MATLAB a otros programas. También puedes hacer que tus p rogramas manipulen datos directamente en archivos - MAT, 73
Ing Sergio Casillas de la Torre Tutorial de Matlab el cúal es el formato de archivo utilizado por MATLAB. Para información acerca de las técnicas utilizadas para importar y exportar datos consulte la sección de Importando y Exportando Datos de la guía de MATLAB ó utilice al comando help de MATLAB.
axis, all, any, bar, break, brighten, clf, for, colormap, conj, contour, contour3, det, diag, else, elseif, eig, expm, fclose, fill, fill3, fft, fftn, function, fopen, fplot, fread, fwrite, grid, gtext, hadamard, hold, if, ifft, ifftn, imag, image, inv, loglog, load, lu, mesh, meshc, meshz, meshgrid, ode23, ode45, peaks, plot, plot3, poly, reshape, rot90, sphere, tril, triu, pascal, polar, real, toplitz, rand, reshape, semilog, semilog y, shading (flat interp faceted), size, sound, stairs, subplot, surf, surfc, surfl, text, title, xlabel, ylabel, while, zeros.
Simulink es una herramienta para el modelaje, análisis y simulación de una amplia variedad de sistemas físicos y matemáticos, inclusive aquellos con elementos no lineales y aquellos que hacen uso de tiempos continuos y discretos. Como una extensión de MatLab, Simulink adiciona muchas características específicas a los sistemas dinámicos, mientras conserva toda la funcionalidad de propósito general de MatLab. Así Simulink no es completamente un programa separado de MatLab, sino un anexo a él. El ambiente de MatLab está siempre disponible mientras se ejecuta una simulación en Simulink. Simulink tiene dos fases de uso: la definición del modelo y el análisis del modelo. La definición del modelo significa construir el modelo a partir de elementos básicos construidos previamente, tal como, integradores, bloques de ganancia o servomotores. El análisis del modelo significa realizar la simulación, linealización y determinar el punto de equilibrio de un modelo previamente definido. Para simplificar la definición del modelo Simulink usa diferentes clases de ventanas llamadas ventanas de diagramas de bloques. En estas ventanas se pue de crear y editar un modelo gráficamente usando el ratón. Simulink usa un ambiente gráfico lo que hace sencillo la creación de los modelos de sistemas. Después de definir un modelo este puede ser analizado seleccionando una opción desde los menús de Simulink o entrando comandos desde la línea de comandos de MatLab. Simulink puede simular cualquier sistema que pueda ser definido por ecuaciones diferenciales continuas y ecuaciones diferenciales discretas. Esto significa que se puede modelar sistemas continuos en el tiempo, discretos en el tiempo o sistemas híbridos. Simulink usa diagramas de bloques para representar sistemas dinámicos. Mediante una interface gráfica con el usuario se pueden arrastrar los componentes desde una librería de bloques existentes y luego interconectarlos mediante conectores y alambre. La ventana principal de Simulink se activa escribiendo simulink en la línea de comandos de MatLab, y se muestra a continuación: Haciendo doble click en cualquiera de las librerías presentes en esta ventana se abrirá otra ventana conteniendo una cantidad de bloques relativos a dicha librería. Para realizar un sistema debe abrirse una nueva ventana de diagrama de bloques seleccionando la opción file del menú principal del Simulink y allí la opción new. En e sta nueva ventana se colocarán todos los bloques interconectados que formarán el sistema deseado.
Ing Sergio Casillas de la Torre Tutorial de Matlab Como ejemplo se ha tomado un generador de ondas seno de la librería de fuentes "sources" y un osciloscopio de la librería "sinks", ambos se unieron mediante un conector usando el ratón. Este sistema se almacena como un archivo - m. Haciendo doble click sobre cada elemento del sistema se pueden ver y modificar sus características. Por ejemplo, al generador seno se le puede modificar su amplitud, frecuencia y fase . Al osciloscopio se le definen las escalas horizontal y vertical. Para ejecutar el programa se usa la opción simulation en el menú de la ventana del archivo - m creado. En este submenú está la opción start que permite ejecutar el programa. También está la opción parameters que activa el panel de control de Simulink en donde se definen los métodos y parámetros usados para la simulación, tal como se muestra a continuación: Al ejecutar el programa seno.m creado mediante simulink, se puede observar la respuesta al hacer doble click en el osciloscopio. Existen numerosos bloques y funciones incorporados en las librerías de simulink que pueden ser empleados para simular cualquier sistema. Por ejemplo, para implementar un sistema que emplea un controlador PID tenemo s:
Para incrementar la velocidad de Simulink se debe instalar el acelerador "Accelerator". Este permite automáticamente generar una versión mejorada de los modelos los cuales correrán diez veces más rápido que el original. El acelerador puede ser usado sobre modelos continuos, discretos en el tiempo y híbridos. El acelerador trabaja generando y compilando un código - C para un modelo dado. Una vez se completa la compilación, la simulación es ejecutada en la ventada de modelos de Simulink exactamente igual que antes sólo que más rápidamente. El propósito del acelerador es aumentar la velocidad de simulación. Si el programa MatLab posee instalado el "Accelerator" podrá iniciarse la acción aceleradora seleccionando la opción simulation en el menú principal del Simulink y dentro de esta seleccionando la opción Accelerate. Esta acción es totalmente transparente en el sentido de que el incremento de la velocidad se presenta sin ningún otro requerimiento por parte del usuario.
Una vez se ha creado un modelo dinámico en Simulink, se puede invocar el generador de código- C que permite convertir el diagrama de bloques implementado en un código C. Este puede ser útil para varios propósitos: puede ser usado para control en tiempo real, simulación en tiempo real o simulación acelerada en tiempo no real. Sus aplicaciones pueden ser control de movimiento, control de procesos, sistemas automotores, equipos médicos, robótica, etc. El código-C es diseñado tal que puede ser ejecutado en tiempo real. No requiere ser escrito manualmente por un programador pues es creado a nivel de diagramas de bloques en Simulink. El código generado puede correr sobre un amplio rango de hardware ubicado en estaciones de trabajo, PC o microprocesadores. Este código es la forma en la que puede usare el Simulink para adquisición de datos.
Managing commands and functions: help - On -line documentation. what - Directory listing of M-, MAT - and MEX -files. type - List M-f i l e . lookfor - Keyword search through the HELP entries. which - Locate functions and files. demo - Run demos. path - Control MATLAB's search path.
Managing variables and the workspace: who - List current variables. whos - List current variables, long f orm. load - Retrieve variables from disk. save - Save workspace variables to disk. clear - Clear variables and functions from memory. pack - Consolidate workspace memory. size - Size of matrix. length - Length of vector. disp - Display matrix or text.
Wor king with files and the operating system: cd - Change current working directory. dir - Directory listing. delete - Delete file. getenv - Get environment value. ! - Execute operating system command. unix - Execute operating system command & return result. d iary - Save text of MATLAB session.
Controlling the command window: cedit - Set command line edit/recall facility parameters. clc - Clear command window. home - Send cursor home. format - Set output format. echo - Echo commands inside script files. more - Control paged output in command window. Starting and quitting from MATLAB: quit - Terminate MATLAB. startup - M -file executed when MATLAB is invoked. matlabrc - Master startup M-f i l e .
General information: info - Information about MATLAB and The MathWorks, Inc. subscribe - Become subscribing user of MATLAB. hostid - MATLAB server host identification number. whatsnew - Information about new features not yet documented. ver - MATLAB, SIMULINK, and TOOLBOX version information. Operators and special charact ers: Char Name HELP topic + Plus arith - Minus arith * Matrix multiplication arith .* Array multiplication arith ^ Matrix power arith .^ Array power arith \ Backslash or left division slash / Slash or right division slash ./ Array division slash
Ing Sergio Casillas de la Torre Tutorial de Matlab kron Kronecker tensor product kron : Colon colon ( ) Parentheses paren [ ] Brackets paren . Decimal point punct .. Parent directory punct ... Continuation punct , Comma punct ; Semicolon punct % Comment punct ! Exclamation point punct ' Transpose and quote punct = Assignment punct == Equality relop < > Relational operators relop & Logical AND relop | Logical OR relop ~ Logical NOT relop xor Logical EXCLUSIVE OR xor Logical characteristics: exist - Check if variables or functions are defined. any - True if any elemen t of vector is true. all - True if all elements of vector are true. find - Find indices of non- zero elements. isnan - True for Not-A- Number. isinf - True for infinite elements. finite - True for finite elements. isempty - True for empty matrix. issparse - True for sparse matrix. isstr - True for text string. isglobal - True for global variables.
Control System Toolbox Commands: Model building: append - Append system dynamics. augstate - Augment states as outputs. blkbuild - Build state-space system from block diagram. cloop - Close loops of system. connect - Block diagram modeling. conv - Convolution of two polynomials. destim - Form discrete state estimator from gain matrix. dreg - Form discrete controller/estimator from gain matrices. drmodel - Generate random discrete model. estim - Form continuous state estimator from gain matrix. feedback - Feedback system connection. ord2 - Generate A,B,C,D for a second-order system. pade - Pade approximation to time delay. parallel - Parallel system connection. reg - Form continuous controller/estimator from gain matrices. rmodel - Generate random continuous model. series - Series system connection. ssdelete - Delete inputs, outputs, or states from model. ssselect - Select subsystem from larger system.
Model conversi ons>: c2d - Continuous to discrete - time conversion. c2dm - Continuous to discrete- time conversion with method. c2dt - Continuous to discrete conversion with delay. d2c - Discrete to continuous- time conversion. d2cm - Discrete to continuous- time conversion with method. poly - Roots to polynomial conversion. residue - Partial fraction expansion.
Ing Sergio Casillas de la Torre Tutorial de Matlab ss2tf - State - space to transfer function conversion. ss2zp - State- space to zero -pole conversion. tf2ss - Transfer function to state- space conversion. tf2zp - Transfer function to zero -pole conversion. zp2tf - Zero- pole to transfer function conversion. zp2ss - Z e r o- pole to state - space conversion.
Model reduction: balreal - Balanced realization. dbalreal - Discrete balanced realization. dmodred - Discrete model order reduction. minreal - Minimal realization and pole-zero cancellation. modred - Model order reduction.
Model realizations: canon - Canonical form. ctrbf - Controllability staircase form. obsvf - Observability staircase form. ss2ss - Apply similarity transfor m. Model properties : covar - Continuous covariance response to white noise. ctrb - Controllability matrix. damp - Damping factors and natural frequencies. dcgain - Continuous steady state (D.C.) gain. dcovar - Discrete covariance response to white noise. ddamp - Discrete damping factors and natural frequencies. ddcgain - Discrete steady state (D.C.) gain. dgram - Discrete controllability and observability gramians. dsort - Sort discrete eigenvalues by magnitude. eig - Eigenvalues and eigenvectors. esort - Sort continuous eigenvalues by real part.
Ing Sergio Casillas de la Torre Tutorial de Matlab gram - Controllability and observability gramians. obsv - Observability matrix. printsys - Display system in formatted form. roots - Polynomial roots. tzero - Transmission zeros. tzero2 - Transmission zeros using random perturbation method.
Time response: dimpulse - Discrete unit sample response. dinitial - Discrete initial condition response. dlsim - Discrete simulation to arbitrary inputs. dstep - Discrete step response. filter - SISO z - transform simulation. impu lse - Impulse response. initial - Continuous initial condition response. lsim - Continuous simulation to arbitrary inputs. ltitr - Low level time response function. step - Step response. stepfun - Step function. Frequency response: bode - Bode plot (frequency response). dbode - Discrete Bode plot (frequency response). dnichols - Discrete Nichols plot. dnyquist - Discrete Nyquist plot. dsigma - Discrete singular value frequency plot. fbode - Fast Bode plot for continuous systems. freqs - Laplace - transform frequency response. freqz - Z- transform frequency response. ltifr - Low level frequency response function. margin - Gain and phase margins. nichols - Nichols plot.
Ing Sergio Casillas de la Torre Tutorial de Matlab ngrid - Draw grid lines for Nichols plot. nyquist - Nyquist plot. sigma - Singular value frequency plot.
Root locus: pzmap - P o l e- zero map. rlocfind - Interactive root locus gain determination. rlocus - Evans root-locus. sgrid - Draw continuous root locus wn,z grid. zgrid - Draw discrete root locus wn,z grid.
Gain selection: acker - SISO pole placement. dlqe - Discrete linear-quadratic estimator design. dlqew - General discrete linear quadratic estimator design. dlqr - Discrete linear-quadratic regulator design. dlqry - Discrete regulator design with weighting on outputs. lqe - L i n e a r-quadratic estimator design. lqed - Discrete estimator design from continuous cost function. lqe2 - Linear quadratic estimator design using Schur method. lqew - General linear-quadratic estimator design. lqr - L i n e a r-quadratic regulator design. lqrd - Discrete regulator design from continuous cost function. lqry - Regulator design with weighting on outputs. lqr2 - Linear quadratic regulator design using Schur method. place - Pole placement. Equation solution : are - Algebraic Riccati equation solution. dlyap - Discrete Lyapunov equation solution. lyap - Continuous Lyapunov equation solution. lyap2 - Lyapunov equation solution using diagonalization.
Demonstrations: ctrldemo - Introduction to the Control Toolbox. boildemo - LQG design of boiler system. jetdemo - Classica l design of jet transport yaw damper. diskdemo - Digital control design of hard disk controller. kalmdemo - Kalman filter design and simulation.
8. APLICANDO MATLAB AL CONTROL DE PROCES OS 8.1 Respuesta en el dominio del tiempo
Para obtener la respuesta de un sistema en el tiempo ante una entrada estándar, debe primero definirse el sistema. Para ello puede definirse en MatLab la función de transferencia propia del sistema o las ecuaciones de estado. La función de transferencia de un sistema es una relación formada por un numerador y un denominador:
En MatLab debe definirse el numerador Y(s) y el denominador U(s) como vectores, cuyos elementos son los coeficientes de los polinomios del numerador y del denominador en potencias decrecientes de S. Por ejemplo, para definir la función de transferencia:
Para determinar la respuesta en el tiempo para una entrada escalón unitario de este sistema se usa el comandos step indicando el vector del numerador y del denominador entre paréntesis. step(num,den) >>step(y,u)
Para el ejemplo anterior tenemos: >>[a,b,c,d]=tf2ss(y,u) a = -0.2500 -1.0000 1.0000 0 b = 1 0 c = 0 1 d = 0
Ing Sergio Casillas de la Torre Tutorial de Matlab Se puede hacer la conve rsión de una ecuación de estado a su equivalente función de transferencia, mediante el comando ss2tf. Se deben especificar los vectores para almacenar los coeficientes del polinomio numerador y del denominador. Su Sintaxis es: [num,den]=ss2tf(a,b,c,d)
Para obtener la respuesta escalón de un sistema a partir de las ecuaciones de estado se usa el comando step con la Sintaxis: step(A,B,C,D)
Para obtene r la respuesta en el tiempo para una entrada impulso unitario se usa el comando impulse, con Sintaxis idéntica a la utilizada con el comando step: Si se define el sistema en MatLab por los polinomios del numerador y denominador de la función de transferencia tenemos: » y=[1 5 4]; » u=[1 6 11 6]; » impulse(y,u)
Si por el contrario el sistema se define en MatLab por las ecuaciones de estado: » [A,B,C,D]=tf2ss(y,u) A= -6 - 1 1 - 6 100 010
Ing Sergio Casillas de la Torre Tutorial de Matlab B= 1 0 0 C= 154 D = 0 » impulse(A,B,C,D)
Ing Sergio Casillas de la Torre Tutorial de Matlab MatLab permite, además de obtener la respuesta en el tiempo para una entrada escalón o impulso, también obtener respuesta para otras entradas tal como rampas o sinusoides. El comando lsim permite obtener la respuesta en el tiempo para un sistema con una entrada u, donde u se define como una función del tiempo. La Sintaxis de este comando es: lsim(A,B,C,D,U,T) usando las matrices de estado o lsim(NUM,DEN,U,T) usando la función de transferencia. Para obtener la respuesta en el tiempo para una función rampa, se define U de la siguiente forma: >>T=0:0.1:10 >>U=T; >>NUM=[1]; >>DEN=[1 0.25 1]; >>[Y,X]=lsim(NUM,DEN,U,T); >>PLOT(T,Y,T,U)
Al hacer U=T se está definiendo la función rampa. T es el vector de tiempo variando desde 0 hasta 10 seg. NUM y DEN son los vectores de los coeficientes decrecientes en potencia de S de los polinomios del numerador y del denominador respectivamente. En la variable Y se almacena la salida del sistema en función del tiempo T. El comando plot permite presentar en la ventana de figuras la variable Y (salida) y la entrada U (rampa) en función del tiempo, obteniéndose:
Para el estudio de un sistema en el dominio de la frecuencia existen tres herramientas disponibles en MatLab como son: los diagramas de Bode, de Nyquist y de Nichols. Para obtener el diagrama de Bode de una función de transferencia, se definen dos vectores cuyos elementos son los coeficientes de los polinomios del nu merador y del denominador en potencias decrecientes de S. Estos vectores son usados en el comando bode con la siguiente Sintaxis: bode(num,den). Se define la función de transferencia:
Otro formato mediante el cual el comando bode presenta el diagrama de bode, es a través de las ecuaciones de estado representadas por las matrices de estado (A,B,C,D). Su Sintaxis es: bode(A,B,C,D). Para especificar un rango deseado de frecuencias en las cuales se desea obtener el diagrama de Bode, se emplea un vector de frecuencias en el que se especifica la frecuencia inicial, el incremento y la frecuencia final. Por ejemplo: >>W=0:0.1:100; >>bode(y,u,W) Este comando muestra el diagrama de Bode entre 0 y 100 rad/s. Otra herramienta de análisis en el dominio en la frecuencia que ofrece MatLab es el diagrama de Nichols. Para obtener el diagrama de Nichols se utiliza el comando nichols, cuya Sintaxis es idénti ca a la del comando bode: nichols(A,B,C,D,W) si se emplean las matrices de estado o nichols(num,den,W) si se emplea la función de transferencia.
Si se define y como el vector de los coeficientes del polinomio del numerador y u como el deldenominador: >>y =[0 0 100]; >>u=[0.04 1 0]; >>nichols(y,u)
Ing Sergio Casillas de la Torre Tutorial de Matlab Otra herramienta de análisis en el dominio en la frecuencia que ofrece MatLab es el diagrama de Nyquist. Para obtenerlo se utiliza el comando ny quist, cuya Sintaxis es idéntica a la del comando bode y nichols: nyquist(A,B,C,D,W) si se emplean las matrices de estado o nyquist(num,den,W) si se emplea la función de transferencia. Si se define y como el vector de los coeficientes del polinomio del numerador y u como el del denominador: >>y=[1]; >>u=[1 6 5]; >>nyquist(y,u) MatLab presenta en la ventana de figuras el diagrama de Nyquist:
Ing Sergio Casillas de la Torre Tutorial de Matlab Para obtener el margen de ganancia, el margen de fase, la frecuencia de cruce de ganancia y la frecuencia de cruce de fase MatLab dispone del comando margin. Las diferentes formas de utilizar este comando son: [Gm,Pm,Wcg,Wcp] = MARGIN(A,B,C,D) retorna los valores de margen de ganancia (Gm), margen de fase (Pm), frecuencia de cruce de ganancia (Wcg) y la frecuencia de cruce de fase (Wcp) cuando se trabaja con las matrices de estado (A,B,C,D). [Gm,Pm,Wcg,Wcp] = MARGIN(NUM,DEN) cuando se trabaja con la función de transferencia. [Gm,Pm,Wcg,Wcp] = MARGIN(MAG,PHASE,W) toma los vectores de magnitud, fase y frecuencia del diag rama de Bode. MARGIN(A,B,C,D) dibuja el diagrama de Bode y muestra con líneas verticales los márgenes deganancia y de fase.
Para obtener el lugar de las raíces, MatLab dispone del comando rlocus. Las diferentes Sintaxis para utilizar este comando son: rlocus(NUM,DEN) calcula y dibuja el lugar de las raíces cuando se trabaja con la función de transferencia donde NUM y DEN son los vectores de los coeficientes en potencia descendiente de S de los polinomios del numerador y denominador de la función de transferencia G(S). MatLab generará automáticamente un conjunto de valores de la ganancia K. rlocus(NUM,DEN,K): calcula y dibuja el lugar de las raíces cuando se trabaja con la función de transferencia y ha sido previamente definido el rango de valores de K. Por ejemplo de 0 a 100 con incrementos de 10: k=0:10:100 R = rlocus(NUM,DEN,K) o [R,K] = rlocus(NUM,DEN) no dibuja el lugar de las raíces pero almacena en la matriz R, de longitud igual al número de elementos de K, la localización de las raíces. R tendrá tantas columnas como raíces existan, estas pueden además ser complejas. rlocus(A,B,C,D), R=rlocus(A,B,C,D,K), o [R,K]=rlocus(A,B,C,D) son equivalentes a las Sintaxis anteriores pero empleando las matrices de estado para hallar el lugar de las raíces.
Ing Sergio Casillas de la Torre Tutorial de Matlab Para la siguiente forma modificada de la ecuación característica de un sistema se desea hallar el lugar de las raíces mediante MatLab: >>num=[0,0,0,1]; >>den=[1,3,2,0]; >>rlocus(num,den) MatLab dispone del comando rlocfind que permite determinar los polos del sistema para una valor determinado de k. Su Sintaxis es: [K,POLES] = rlocfind(num,den) permite determinar los polos para un valor determinado de k, cuando se trabaja con la función de transferencia. Por medio del curso en el lugar de las raíces se selecciona una localización, MatLab retorna el valor de k para esta localización y los polos asociados a esta ganancia.
Cuando se trabaja con las matrices de estado, las Sintaxis para el comando rlocfind es: [K,POLES] = rlocfind(A,B,C,D). Al ejecutar el comando rlocfind con la función de transferencia anterior, MatLab activa la ventana de figuras en espera de que el usuario seleccione un punto del lugar de las raíces mediante el cursor. En este caso el punto seleccionado fue 2.4623 en la parte real y - 0.0132 en la parte imaginaria. » [k,poles]=rlocfind(num,den) Select a point in the graphics window selected_point = -2.4623 - 0.0132i k = 1.6655 poles = -2.4625 -0.2688 + 0.7773i -0.2688 - 0.7773i Para seleccionar el punto en el cual calcular los polos del lugar de las raíces sin usar el cursor se agrega un parámetro al comando rlocfind. Este debe ser el punto o los puntos en donde se desea tomar el valor de k. La nueva Sintaxis es: [K,POLES] = rlocfind(A,B,C,D,P) o [K,POLES] = rlocfind(num,den,P) P debe definirse previamente indicando la parte real e imaginaria del mismo. Por ejemplo: P=3+0i o P=1 -0.555i.
Para implementar los diferentes tipos de controladores (P, PD, PI, PID) en MatLab se hace uso de la función de transferencia propia del sistema a objeto de estudio. Si dicho sistema es de la forma:
donde G(S) es la función de transferencia de la planta o proceso; mientras que C(S) es la función de transferencia del controlador. Para el caso del controlador proporcional, C(S)=Kp, que es una constante o valor escalar. El controlador PI es C(S)=Kp + Ki/S que puede representarse como una relación ente dos polinomios. El controlador PID es C(S)=Kp + Ki/S + Kd S que se representa como:
que es de nuevo una relación entre dos polinomios. Los coeficientes decrecientes en potencias de S de estos polinomio pueden ser almacenados en vectores en MatLab. Si se multiplica el controlador C(S) por la función de transferencia del proceso o planta G(S) se formará la función de transferencia de lazo abierto. Por ejemplo un G(S) puede ser:
Ing Sergio Casillas de la Torre Tutorial de Matlab Para obtener la respuesta en lazo abierto ante una entrada escalón unitario tenemos: >>Kp=50; >>Ki=1; >>Kd=10; >>num=[Kd Kp Ki]; >>den=[1 10 20 0 0]; >>step(num,den)
Para obtener la respuesta de lazo cerrado en el tiempo para una entrada escalón unitario se emplea el comando cloop, el cual genera los polinomios del numerador (numc) y denominador (denc) de la función de transferencia de lazo cerrado con realimentación unitaria a partir de los polinomios de la función de transferencia de lazo abierto (num y den). Su Sintaxis es: [numc,denc]=cloop(num,den,sign) El signo de la realimentación viene dado por sign. Para el ejemplo anterior, tenemos: >>Kp=500; >>Ki=1; >>Kd=100; >>num1=[Kd Kp Ki]; >>den1=[1 0]; >>num2=1; >>den2=[1 10 20 0]; >>[numc,numd]=cloop(conv(num1,num2),conv(den1,den2), - 1); >>step(numc,denc)
Se usa el comando conv para obtener la convolución y multiplicación polinomial de dos vectores. La salida obtenida mediante el comando step se muestra a continuación:
Para cotejar sus diagramas de Bode: >>bode(num,den) donde n u m y d e n son vectores que contienen los coeficien tes del numerador y denominador de H(s) en orden de potencias descendentes de s.
Nota: Esto da las curvas exactas, no las aproximaciones asintóticas con líneas rectas.
Precaución: El punto "." puede significar operación elemento -por -elemento o punto decimal. Cuando escribimos un dígito pegado al punto como "2.", el interpretador cree que es el número "2.0". Entonces si queremos calcular A2B2, donde A y B son arreglos y no matrices (o sea, queremos operación elemento -por- elemento), debemos escribir >>A.^2 .*B.^2 (notar el espacio después del primer 2) y no >>A.^2.*B.^2
Ing Sergio Casillas de la Torre Tutorial de Matlab Para cambiar el c olor de trasfondo de la gráfica: >> whitebg('c') donde c es el código del color descrito en help plot.
Para establecer propiedades de la gráfica, es más fácil hacerlo al crearla que después. Por ejemplo, para graficar con una línea gruesa, >>plot(x,y,'linewidth',3) (En el momento de creación) >>set(get(gca,'children'),'linewidth',3) (Después de creada)
Más de este usuario16619641 Libro Serway Electric Id Ad y Magnetismo5ta Edicionespanol16619641 Libro Serway Electric Id Ad y Magnetismo5ta EdicionespanolrelacionesrelacionadasrelacionesrelacionadasRAZONRAZONRazon de Cambio BicicletaRazon de Cambio Bicicleta5052034-CALCULO-DIFERENCIAL-E-INTEGRAL-I-FAS1-RAZON-DE-CAMBIO-5052034-CALCULO-DIFERENCIAL-E-INTEGRAL-I-FAS1-RAZON-DE-CAMBIO-Trabajo de Dibujo TecnicoTrabajo de Dibujo TecnicoSISTREPRESISTREPREPresentaci Online AsPresentaci Online Asmatlab70primeromatlab70primeroINTRODUCC..INTRODUCC..

References: resolución 
 resolución 
 Resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 Resolución