Transferencia de Masa I
Grupo #3 de Laboratorio.
Tema: Efecto de la temperatura en la estimación de la difusividad másica mediante la
ley de Fick.
Integrantes:
- Andrade María
- Baque Adrián
- Chiluisa Alex
- Garzón Mishell
Código de Matlab ejecutado durante la práctica de laboratorio
clear, clc, close all % Limpiar la memoria
fprintf('Universidad de las Fuerzas Armadas ESPEL\n')
fprintf('Procesamiento de Imagenes\n')
fprintf('Quinto Petroquimica\n')
fprintf('Grupo 3\n')
%Introducir el nombre de la imagen y el formato en
%('nombre,formato'), ejemplo ('c3.jpeg)%
img1 = imread('c3.jpeg');
%imshow(img);
AB1 = imcomplement(img1);
%figure(1)
AB2 = rgb2gray(AB1);
imshow(AB2)
intensidad1=improfile;
%Seleccionar n para el radio inicial de la difusión para cada temperatura
los
%cuales son: para 22°C es 0.000885 , para 30 °C es 0.00115 y para 60°C es
%0.000975%
n=0.000975;
a=0;
%El código dentro de for-end sirve para mostrar los datos en una columna
%los cuales son tomados para realizar el resto de cálculos en el Excel%
for i=1:length(intensidad1)
longitud1(i)=a;
a=a+n;
end
longitud1=longitud1';
a=0;
%close all
%La figura (2), se ubicará acorde a la posición que el usuario
%determine en el comando subplot jugando con los nnúmeros (#,#,#). Luego,
el comando plot, grafica la
%longitud vs la intensidad y los puntos se pueden modificar gracias a 'g+'
%en el cual, 'g' es la inicial del nombre del color en inglés y '+' es el
sibolo que
%aparecerá en cada punto, ejmplo 'b*' %
pág. 1
figure(2)
subplot (1,1,1), plot(longitud1,intensidad1,'g+');
%Para poner nombre a los ejes, se utiliza los comandos xlabel & ylabel que
%son los que proporcionan el nombre al eje X y Y respectivamente%
xlabel('Longitud (cm)')
ylabel('Intensidad')
%Finalmente, para colocar un nombre a la gráfica, se utiliza el comando
%tittle tal como se muestra a continuación%
title('Gráfica de puntos a utilizar')
%Para la figura (1), el comando inshow grafica la intensidad de la imagen%
figure(1)
imshow(img1),title('Difusión de Masa');
subplot (1,1,1), plot(intensidad1), grid on ,grid minor;
xlabel('Longitud (pixeles)')
ylabel('Intensidad')
title('Gráfica de difusión');
[f1,c1]=size(img1);
%El comando rgb2gray cambia el color de una imagen de colores a blanco y
%negro%
a=rgb2gray(img1);
figure(3)
subplot(1,3,1)
imshow(a)
xlabel('Grupo #3')
title('Imagen Blanco y Negro')
t=(img1);
figure(3)
subplot(1,3,3)
imshow(t)
xlabel('Grupo #3')
title('Imagen a color')
subplot(1,3,2)
% Muestra una gráfica del histograma. Si la imagen de entrada es una
% imagen indexada, el histogrma muestra la distribucion de los valores de
% pixel por encima de un colobar del mapa de colores%
imhist(a)
title('Perfil de Intensidad de toda la imagen')
pág. 2