PDF de programación - MATLAB comandos básicos

MATLAB

comandos básicos
comandos básicos

Catedra Calculo de Reactores

A-10

Prof. Alberto Quintero

Definición de Variables

Para definir una variable en el programa solo necesitamos
escribir el comando

nombre de la variable=valor de la variable

Si deseamos que el resultado de una sentencia no aparezca en
pantalla colocamos “;” al final de la sentencia

Declaración de Matrices

Los elementos de cada fila se separan por “,” o espacio, mientras
que las columnas se diferencia por el “;”. Los elementos que
integran la matriz deben ir entre corchetes.

nombre de la matriz=[elementos de la fila 1;elementos de la fila 2;…..]

Declaración de Matrices

También se pueden crear vectores en forma rápida utilizando

nombre del vector=inicio:paso:fin

O vectores equiespaciados, con un número de elementos fijos

nombre del vector=linspace(primer elemento,elemento final,numero de elementos)
nombre del vector=linspace(primer elemento,elemento final,numero de elementos)

Indexación de Matrices

Si se define A como una matriz de mxn se puede acceder a sus
elementos a través de los siguientes comandos:

A(i,j) muestra el elemento de la fila “i” y la columna “j”.

A(i,:) muestra los elementos de la fila “i”.

A(:,j) muestra los elementos de la columna “j”.

A([1 3],j) muestra los elementos de la fila 1 y 3 de la columna “j”.
A([1 3],j) muestra los elementos de la fila 1 y 3 de la columna “j”.

A(i,[1 4]) muestra los elementos de las columnas 1 y 4 de la fila “i”.

También se pueden asignar valores específicos a elementos de
una matriz:

A(1,2)=17 guarda el número 17 en la fila 1, columna 2.

A(i,:)=17 guarda el número 17 en todos los elementos de la fila “i”.

A(:,j)=17 guarda el número 17 en todos los elementos de la columna “j”.

Matrices Especiales

En ciertas ocasiones se necesita definir matrices especiales,
como identidades, matrices nulas, unitarias, diagonales, etc.

A’ matriz transpuesta de A.

eye(n) crea una matriz identidad de nxn.

eye(n,m) crea una matriz identidad de mxn.

diag(V) crea una matriz diagonal con los elementos del vector V.
diag(V) crea una matriz diagonal con los elementos del vector V.

diag(V,n) crea una matriz diagonal desplazada en n con los elementos del
vector V. Si n>0 la diagonal se desplaza a la derecha, si n<0 la
diagonal se desplaza a la izquierda.

diag(A) crea un vector con la diagonal de la matriz A.

ones(m,n) crea una matriz de elementos de valor 1 de mxn.

zeros(m,n) crea una matriz de elementos de valor 0 de mxn.

rand(m,n) crea una matriz de valores aleatorios de mxn.

Operaciones con Matrices

En el caso de realizar diversas operaciones con matrices
previamente definidas se utilizan diversos operadores entre
ellos:

A*B multiplicación matricial de A por B.

A.*B multiplicación de los términos de A por los términos de B.

A/B multiplicación matricial de A por la inversa de B.
A/B multiplicación matricial de A por la inversa de B.

A\B multiplicación matricial de B por la inversa de A.

A./B división de los términos de A por los términos de B.

A+B suma de los términos de A con los de B.

A-B resta de los términos de A con los de B.

dot(U,V) producto punto de los vectores U y V.

Cross(U,V) producto cruz de los vectores U y V.

OJO: Para realizar cualquier operación elemento a elemento se antecede un
“.” al operador, ejemplo: A.^2 eleva cada elemento de la matriz A al cuadrado

Solución de Sistemas de

Ecuaciones Lineales

Supongamos que se tienen las siguientes ecuaciones lineales

y=x+3

y=-x+2

Recordando un poco de matemáticas, sabemos que el sistema
Recordando un poco de matemáticas, sabemos que el sistema
se puede escribir de la siguiente forma para resolverlo:

y-x=3

y+x=2

A*x=b → x=inv(A)*b

Problemas para practicar

1.-Se tiene una torre de destilación donde ingresa una corriente
de composición 30%p/p en Etanol y 70%p/p en Metanol, la
corriente de tope tiene una composición de 90%p/p en metanol
y la corriente de fondo tiene una composición de 5%p/p en
metanol.Si se alimentan 100Kg/min a la torre, cuanto seran los
flujos de salida de tope y fondo?

Resultados:

Tope=76,47Kg/min

Fondo=23,53Kg/min

Tienen 20 minutos para resolver el ejercicio.