MATLAB INGENIERIA QUIMICA-UNICA
PRACTICA: FUNCIONES
Una función tiene nombre, valor de retorno y argumentos.
Los argumentos de cada función van a continuación del nombre entre paréntesis (y
separados por comas si hay más de uno). Los valores de retorno son el resultado de la
función y sustituyen a ésta en la expresión donde la función aparece.
Los nombres de funciones deben estar en minúsculas.
Variable_resultado = función (argumento)
x = linspace(-2,2)
y = sin (x)
z = sin (x*pi/180)
v = -10
w = log (abs (v))
MATLAB nos ofrecerá ayuda en cualquier momento escribiendo help, podemos escribir
help elfun (elfun: elementary functions) para que nos informe de mas funciones que tiene
disponibles.
>> help elfun
Elementary math functions.
Trigonometric.
Exponential.
Complex.
Rounding and remainder
>> help linspace
linspace Linearly spaced vector.
linspace(X1, X2) generates a row vector of 100 linearly equally spaced points between X1 and
X2.
linspace(X1, X2, N) generates N points between X1 and X2.
For N = 1, linspace returns X2.
Class support for inputs X1,X2: float: double, single
Numeros Complejos.- Cuando representamos un número complejo con una parte real y
una parte imaginaria, como por ejemplo, en 2 + i3, estamos usando notación
rectangular.
Un número complejo podría describirse también con un ángulo θ y un radio r relativo al
origen. Esta forma se denomina notación polar.
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
ARCHIVOS SCRIPT
Un archivo.m consiste de una secuencia de instrucciones normales de MATLAB. Se
llaman archivos.m. Esto es así porque siempre tienen una extensión de .m como la última
parte de su nombre de archivo.
Hay dos tipos de archivos.m: archivos de comandos (Script) y archivos de funciones.
Archivo Script o Archivo de Comandos.- Es una secuencias de comandos MATLAB.
Cuando se ejecuta un archivo de comandos o también llamado archivo script, MATLAB
ejecuta los comandos en el orden en que éstos han sido escritos, igual que si se
ejecutaran uno a uno en la ventana de comandos.
La función disp de MATLAB se utiliza para mostrar un escalar, un mensaje (string), un
vector o una matriz.
Recuerda que el empleo del signo ; (punto y coma) suprime automáticamente la
exhibición del contenido de la variable.
>> a=42;
>> b=6.78;
>> x=5;
>> y=x*x;
>> z=2*x+cos(x)-x/y+log(x/4);
>> disp('Buenas Tardes')
Buenas Tardes
>> disp(z)
10.3068
>> disp(y)
25
>> d=1:0.5:3;
>> disp(d)
1.0000 1.5000 2.0000 2.5000 3.0000
>> A=[1 2 3; 4 5 6; 7 8 9];
>> disp(A)
1 2 3
4 5 6
7 8 9
>> disp(pi)
3.1416
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
La función fprintf puede escribir en pantalla una combinación de datos y/o texto.
Para la escritura sólo de texto se utiliza la función de igual manera que disp pero
finalizando la cadena en el carácter \n (salto de línea).
fprintf(‘texto\n’)
Para escritura de texto y datos en pantalla, o sólo datos, se utiliza la función con la
siguiente sintaxis
fprintf(‘formato’, variables)
>> clc
>> a=42;
>> b=6.78;
>> fprintf ('El valor es %d\n',a);
El valor es 42
>> clc
>> fprintf ('Una es %d y la otra es %f\n',a,b);
Una es 42 y la otra es 6.780000
La función input de MATLAB permite aceptar datos de entrada a través del teclado.
Una alternativa más cómoda de asignar valores a una variable es utilizando el
comando input.
>> clear all
>> clc
>> r = input ('Teclee el radio: ')
Teclee el radio:
>> Nombre=input('Introduce tu nombre y apellidos: ','s')
Introduce tu nombre y apellidos:
>> nombre=input('Introduce tu nombre y apellidos: ')
Introduce tu nombre y apellidos:
Por ejemplo: Sea la ecuacion de segundo grado x2 +x -6=0
a=input('primer coeficiente, a: ');
b=input('segundo coeficiente, b: ');
c=input('término independiente, c: ');
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
EJERCICIOS_FUNCIONES
1. Evalúe las siguientes expresiones, y luego verifique sus respuestas
introduciendo las expresiones en MATLAB.
A. round (-2.6)
B. fix (-2.6)
C. sign (-2.6)
D. rem (15,2)
E. sqrt (3^2 + 4^2)
F. log(100)
G. log10(100)
H. log10 (100) + log10(0.001)
2. Calcula el seno de x :
x en radianes,
x = sin (30)
X en grados.
X = sind (30)
3. Escriba expresiones MATLAB para calcular lo siguiente:
>> cos(0:pi/2:2*pi)
>> exp(log(2))
>> (exp(5)+exp(-5))/2
4. Escriba los siguientes ejemplos de números complejos:
c1 = 1 – 2i
c2 = 3* (2 - sqrt (-1)*3)
c3 = sqrt (-2)
c4 = (c1 + c2) / c3
5. Convierta los valores complejos a la forma polar. Verifique sus respuestas
usando funciones MATLAB.
a) 3 - i2 b) 0.5 + i
Escriba los resultados en notación polar.
6. Convierta los valores exponenciales complejos a la forma rectangular. Verifique
sus respuestas usando funciones MATLAB.
a) ei b) 0.5ei2.3
Escriba los resultados en notación rectangular
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
7. Si A=2 1
4 3
Halle el rango y la inversa de A.
8. Halle la determinante de las siguientes matrices:
A= 1 3 0
-1 5 2
1 2 1
B= 9 2 1 -1
-6 8 -8 7
-9 -1 -5 2
8 -6 -1 -6
9. Sea A la matriz siguiente: 2 1 0
-2 5 -1
3 4 6
Gire la matriz A 90° en dirección contraria a las agujas del reloj.
10. Sea A la matriz siguiente: 1 2
4 8
-2 0
a) Invierta la matriz A de izquierda a derecha
b) Invierta la matriz A de arriba hacia abajo.
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
EJERCICIOS_ARCHIVOS SCRIPT
1.- Escriba un archivo script que calcule el área de un círculo a partir de su radio
Solución:
- En la ventana del editor/depurador, escribir las siguientes instrucciones:
% Calculo del area de un circulo a partir de su radio
% Notese que el signo de porcentaje sirve para agregar comentarios
% Nuestro primer archivo script
%
% Entrada de datos
r = input ('Ingrese el radio del circulo: ');
% Proceso mediante el empleo de una formula que calcule el area
A = pi*r^2;
% Mostrar la salida de datos
fprintf('El area del circulo es %f\n',A);
Guardar como:
Ejecutar Para ejecutar el archivo guardado, solo debe digitar el nombre y pulsar la
tecla ENTER
2.- Escriba un archivo script que evalúe con 15 dígitos la función sin(x)/x para valores de
x muy próximos a cero. En base a estos resultados, infiera una respuesta para el cálculo
del limx 0 sin(x)/x.
Solución:
- En la ventana del editor/depurador, escribir las siguientes instrucciones:
% Solucion del Ejemplo 2
format long
% formato de 15 digitos
% Entrada de datos
% Definimos los valores de la variable independiente
x=[0.1, 0.01, 0.001, 0.0001];
% Formula a emplear en el proceso
% Calculamos los valores deseados
y= sin(x)./x;
% Salida de datos
disp('y = ');
disp (y)
Guardar como:
Ejecutar
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
3.- Cree el siguiente archivo script :
% El programa muestra el uso de algunas funciones matematicas.
% Punto y coma suprime la exhibición de las variables
clc
x = [2 4 6]
y = sin(x)
t = 5;
s = t+2;
w = s-4;
% Las variables s, t, w no son exhibidas.
v=s^t+sin(s)/cos(t)+log(w);
disp('v = ')
disp(v)
z=v+exp(t/2)+s*sqrt(w);
disp('z =: ')
disp(z)
fprintf ('En un formato diferente z= %f\n',z)
Guardar el archivo
Ejecute el archivo
4. Cree un archivo script que permita calcular la Presion de un gas usando la ecuacion
de los gases ideales.
Solucion: En el editor de texto , escribimos las siguientes instrucciones:
% Calcula la P de un gas usando la ecuacion de los gases ideales
T=input (' TEMPERATURA EN K = ');
V=input (' VOLUMEN EN L = ');
n=input (' n en mol = ');
R=input (' Constante R EN atm-l/molK = ');
Vm= V/n;
P=R*T/Vm;
disp('P =');
disp(P);
Guardar el archivo
Ejecute el archivo con los siguientes datos: 0.652 mol de oxígeno gaseoso que ocupa un
cilindro de 10.0 litros a 30° C, R = 0.0821 atm-L/mol °K
Escriba el resultado
5. Cree un archivo script que permita calcular las constantes a y b
27 R 2Tc 2 RTc
a b
64 Pc 8 Pc
Guardar el archivo
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
Ejecute el archivo con los siguientes datos:
R = 0.082 Pc = 45.8 Tc = 190.5
Trascriba las líneas de instrucciones y el resultado obtenido
6. Cree un archivo script que permita calcular el volumen de un gas usando la ecuacion
de los gases ideales.
Guardar el archivo
Ejecute el archivo con los siguientes datos: 2.15 mol de oxigeno gaseoso a 27°C y 1.25
atm. R = 0.0821 atm-L/mol °K
Trascriba las líneas de instrucciones y el resultado obtenido
7. Cree un archivo script que permita calcular las moles de un gas usando la ecuacion
de los gases ideales.
Guardar el archivo
Ejecute el archivo con los siguientes datos: cilindro de 5.25 litros a 27° C y 1.30 atm.
R = 0.0821 atm-L/mol °K
Trascriba las líneas de instrucciones y el resultado obtenido
8.- Cree un archivo script que permita calcular la Presion de un gas usando la ecuacion
de los gases ideales.
Guardar el archivo
Ejecute el archivo
Trascriba las líneas de instrucciones y el resultado obtenido
9.- Cree un archivo-m que permita calcular la Presion de un gas usando la ecuacion de
Van Der Waals
Guardar el archivo
Ejecute el archivo
Trascriba las líneas de instrucciones y el resultado obtenido
10.- Cree un archivo-m que permita calcular la Presion de un gas usando la ecuacion de
Redlich-Kwong.
Guardar el archivo
Ejecute el archivo
Trascriba las líneas de instrucciones y el resultado obtenido
Docente: Ing. Felipe Estuardo Yarasca Arcos
MATLAB INGENIERIA QUIMICA-UNICA
Datos para los ejercicios 8, 9 y 10:
Para el Gas metano:
T = 298, V = 5, n = 50, Tc = 190.5, Pc = 45.8, R= 0.0821
Vm= V/n
De tablas termodinamicas: a = 2.25, b = 0.0428 Para la ecuacion de VDW
ECUACION DE LOS GASES IDEALES
RT
P
Vm
P=R*T/Vm
ECUACION DE VAN DER WAALS (VDW)
La constante a se halla en las tablas termodinámicas
La constante b se halla en las tablas termodinámicas
RT a
P
Vm b Vm 2
P=((R*T)/(Vm-b)-(a/Vm^2))
ECUACION DE REDLICH-KWONG (R-K)
Los valores de a y b se hallan por formula
0.42748 R 2Tc 2.5 0.08664 RTc
a b
Pc Pc
RT a
P 1/ 2
Vm b T Vm(Vm b)
Docente: Ing. Felipe Estuardo Yarasca Arcos
