Post on 06-Dec-2015
description
Operador punto
Este operador le indica a Matlab que se desea que se hagan todas las operaciones elemento a elemento de dos vectores.
Ejemplo:Ki=A.*Psat;
Comando interpolación cubica
Ejemplo:%Formula interpolaciónra=[0.010 0.005 0.002 0.001]x=[0 0.2 0.4 0.6]%Para una x de 50Y2=interp1(x,ra,0.5,’cubic’)
*Nota: si le quitamos el ’cubic’ se hará una interpolación lineal.
Comando ode15s
function main_principal
global poner parámetros que se usan en la función principal y la función a utilizar%ParametrosPoner parametros %Condiciones de simulaciónPoner condiciones de simulación (condición(es) inicial(es) e intervalo de simulación)
*Nota: si son varias condiciones iniciales se ponen como vectores y el intervalo siempre es en forma de vector.
%Aplicación del ode15s[variable independiente, variable dependiente]=ode15s(@función_a_utilizar, intervalo_de_simulación, condición_inicial);
%Si quiero obtener una variable del vector respuesta (dependiente)na= variable dependiente (:,#columnda);
Ejemplo:Temperatura= xs(:,2);
%Si quiero obtener el ultimo valor de un vector uso el comando “end” Ejemplo:temp=Temperatura(end)
%Si quiero graficar puedo poner solo una gráfica (comando plot) y si quiero poner varias (adiciono comando subplot)
Ejemplo:figure (1)subplot(largo, ancho, cuadrante alias posición por ejemplo 1)plot(ts,na)xlabel('Tiempo (hrs)')
ylabel('Concentración (lbmol/m3)')title 'Tiempo vs 'Concentración'legend('C_A','C_B', 'C_C',0)grid onsubplot(largo, ancho, cuadrante alias posición por ejemplo 2)plot(ts,Temperatura,'red')xlabel('Tiempo (hrs)')ylabel('Temperatura (°R)')title 'Tiempo vs Temperatura'grid on
function dcidt= función_a_utilizar (variable independiente ejemplo “t”, variable dependiente ejemplo “x”)global poner parámetros que se usan en la función principal y la función a utilizar%Renombrando variables
Ejemplona= x(1);T= x(2);
%Ecuaciones auxiliaresSi es que tiene. %Ecuaciones diferenciales
Ejemplodnadt = -k*na;dTbdt = (DHr*-k*na-U*A*(T-Tm))/Cv;dcidt=[dnadt;dTbdt];
Comando fsolve
%Valores de Exitflag
1 Function converged to a solution x.
2 Change in x was smaller than the specified tolerance.
3 Change in the residual was smaller than the specified tolerance.
4 Magnitude of search direction was smaller than the specified tolerance.
0 Number of iterations exceeded options.MaxIter or number of function evaluations exceededoptions.MaxFunEvals.
-1 Output function terminated the algorithm.
-2 Algorithm appears to be converging to a point that is not a root.
-3 Trust region radius became too small (trust-region-dogleg algorithm).
-4 Line search cannot sufficiently decrease the residual along the current search direction.
function main_programax0=supuesto(s) inicial(es)options=optimset('Display','iter');[x,fval,exitflag]=fsolve(@funcion_a_usar,x0,options);
function f= funcion_a_usar (poner variable a encontrar raiz)
%ParametrosSi es que tiene
%Ecuaciones auxiliaresSi es que tiene %Ecuación(es) (Poner ecuaciones)f(1)=…..;
f2)=…..;
f(3)=…..;
f=f’;
Comando elección de incisos
%Eleccion del incisoeleccion = input('Inciso a = 1, Inciso b = 2, Inciso n = # ==> '); %Parametrosif eleccion==1%Lo que quieres que pase en este incisoelseif eleccion==2%Lo que quieres que pase en este incisoend
Graficar en 3D
figure (1)surf(eje X, eje Y, eje Z)xlabel('eje X')ylabel('eje Y')zlabel('eje Z')title 'Poner título'
Comando linspace
ci=linspace(valor inicial, valor final o mismo valor inicial, tamaño deseado del vector a generar)
Este comando crea un vector con las mismas distancias entre cada elemento.