CODIGOS - SCRIPT - FUNCIONES - OTROS

Esta sección está destinada para el almacenamiento y descarga de Script, Funciones, Códigos u otros elementos  de  Matlab, los cuales he ido creando con el tiempo, espero además, que algunos de mis amigos y amigas que quieran compartir sus conocimientos de programación de este lenguaje también lo hagan a través de sus comentarios.

En esta sección algunos de los archivos serán explicados con detalle para su uso eficiente y  apoyo al aprendizaje de los amigos y amigas que recién comienzan. Por lo demás, la mayoría de estos códigos no serán códigos con elevados detalles, como lo son la aparición de ventanas, mensajes de error y otras elementos de programación, eso lo dejaremos para otra sección, en donde analizaremos algunas forma de como fortalecer los códigos basicos que ud crearán por si solos. Pueden citar las líneas de programación que aquí encontrarán refiriéndose a este blog. Por otro lado, espero  posteriormente crear un sumidero o depositorio de todos los códigos que se han subido a este apartado,  además de su pequeña descripción.  Espero amigos y amigas que uds puedan aportar sus propias creaciones para fortalecer este espacio. Saludos a todos y espero les guste. (Dejen sus comentarios, de una u otra manera ayudan).

 ================================================================

 SCRIPTS

(1)   Ecuación de Van de Walls

Este script grafica las isotermas y las isobaras a las:

Temperaturas [ 300 400 500 600] grados celsius y 

Presiones      [45 55 65 75] atm

% Scrip Ecuación de Van der Walls

% (P +  (a / v^2)) * (V - b) = R*T

% Presión (atm); Volumnen (ltr); Temperatura (Kelvin); a y b son parametros

% que dependen del fluido, ejemplo para el benceno

% a = 18,78 (atm*ltr/mol^2)

% b = 0.128 (ltr/mol)

% 1 °Centigrado = 273.15 °Kelvin

 clear, clc, home 

% Definición de constantes para el benceno, puede ud cambiarlas para    %  otro gas (a y b)

a = 18.78;       % [Atm * Ltr^2 / Mol^2]

b = 0.1208;     % [Ltr / Mol]

R = 0.0821; % (atm*ltr/mol*°K)

% T define la o las temperaturas que se desean modelar (ver función lispace)

T = linspace(300,300,4)+ 273.15; %°KELVIN

V = 1:100; % rango del volumen

% Ecuación con Parametros de Temperatura Constante

ecu1 = R./(V-b); % primera parte de la ec una vez despejada la ecuación para P

ecu2 = a./V.^2;  % segunda parte de la ec una vez despejada la ecuación para P

P = []; % Matriz vacia para ir almacenanado los datos

for i =1:length(T)

 P(i,:) = (ecu1*T(i) - ecu2);

end

% ploteo de volumen y presión a las temperaturas antes descritas

subplot(1,2,1), plot(V,P(1,:), 'b-',V,P(2,:), 'g--',V,P(3,:), 'r:',V,P(4,:), 'k.-.')

axis([0 50 0 15]), title('Benceno (Isotermas)')

xlabel('Volumen [Ltr]','FontWeight','bold')

ylabel('Presión [Atm]','FontWeight','bold')

legend('T° = 300 °C','T° = 400 °C','T° = 500 °C','T° = 600 °C'), grid on

%% Ecuaciones con Parametros de Presión Constante.

P1 = linspace(45,75,4);

T1 = [];

for i =1:length(P1);

 T1(i,:) = ((P1(i)+(a*V.^-2)).*(V-b))./R;

end

% ploteo de volumen y temperatura a las presiones antes descritas

subplot(1,2,2), plot(V,T1(1,:), 'b-',V,T1(2,:), 'g--',V,T1(3,:), 'r:',V,T1(4,:), 'k.-.')

axis([0 100 0 7*10^4 ]), title('Benceno (Isobaras)'), 

xlabel('Volumen [Ltr]','FontWeight','bold')

ylabel('Temperatura [°K]','FontWeight','bold')

legend('P = 45 Atm','P = 55 Atm','P = 65 Atm','P = 75 Atm'), grid on

 

Descargar script en archivo .m desde Rapidshare