Reacciones de autocatalisis Grupo 9A - Historial de revisiones

David Carmona Rodriguez en 10:30 6 mar 2015

2015-03-06T10:30:28Z

David Carmona Rodriguez en 10:22 6 mar 2015

2015-03-06T10:22:00Z

David Carmona Rodriguez en 10:18 6 mar 2015

2015-03-06T10:18:05Z

David Carmona Rodriguez en 22:57 5 mar 2015

2015-03-05T22:57:32Z

David Carmona Rodriguez en 22:41 5 mar 2015

2015-03-05T22:41:52Z

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434

David Carmona Rodriguez en 22:34 5 mar 2015

2015-03-05T22:34:31Z

David Carmona Rodriguez en 22:33 5 mar 2015

2015-03-05T22:33:43Z

David Carmona Rodriguez en 22:30 5 mar 2015

2015-03-05T22:30:45Z

David Carmona Rodriguez en 22:24 5 mar 2015

2015-03-05T22:24:04Z

David Carmona Rodriguez en 22:11 5 mar 2015

2015-03-05T22:11:34Z

@@ Línea 100: / Línea 100: @@
 %DATOS DEL PROBLEMA
 t0=0; tN=10; y0=0.01; K1=1; K2=1.01;h=0.1;
 %GENERACION DEL VECTOR TIEMPO t EN FUNCION DE h
 t=t0:h:tN;
@@ Línea 149: / Línea 148: @@
 t=t0:h:tN;
 %calculo del numero de intervalos
-N=(tN-t0)/h; %calculo del numero de intervalos
+N=(tN-t0)/h;
 %preparacion del vector solucion aproximada
 y=zeros(1,length(t));

@@ Línea 442: / Línea 442: @@
 Mediante este procedimiento se obtienen las siguientes graficas:
 [[Archivo:Definita.jpg|800px|sinmarco|centro]]
 ===Planteamiento, resolución por Heun e interpretación del PVI (Apartado 7)===

@@ Línea 442: / Línea 442: @@
 Mediante este procedimiento se obtienen las siguientes graficas:
 [[Archivo:Definita.jpg|800px|sinmarco|centro]]
 ===Planteamiento, resolución por Heun e interpretación del PVI (Apartado 7)===

← Revisión anterior		Revisión del 22:57 5 mar 2015
Línea 227:		Línea 227:
	====Interpretacion de las graficas y soluciones====		====Interpretacion de las graficas y soluciones====

−	En una reacción autocatalítica si comenzamos con una cantidad pequeña de B, la velocidad de reacción aumentará a medida que se vaya formando más B. En el otro extremo, cuando haya desaparecido prácticamente todo el componente A, la velocidad ha de tender a cero. ~~'''~~Este comportamiento se puede apreciar en la gráfica anterior, en la que la velocidad varía a lo largo de una función cuyo máximo corresponde a concentraciones iguales de A y de B~~'''?~~.	+	En una reacción autocatalítica si comenzamos con una cantidad pequeña de B, la velocidad de reacción aumentará a medida que se vaya formando más B. En el otro extremo, cuando haya desaparecido prácticamente todo el componente A, la velocidad ha de tender a cero. Este comportamiento se puede apreciar en la gráfica anterior, en la que la velocidad varía a lo largo de una función cuyo máximo corresponde a concentraciones iguales de A y de B.

@@ Línea 259: / Línea 259: @@
 {{matlab|codigo=
-%Metodo de Euler para sistemas
+%METODO DE EULER PARA SISTEMAS
-%Datos y condiciones iniciales:
 k=1;t0=0; tN=10;h=0.1; y0=0.01; x0=1;
-%Definimos el vector del tiempo, t pertenece a [0,10] recorrido con paso h.
+%generacion del vector tiempo, conocido h
 t=t0:h:tN;
-%PREPARACION DEL VECTOR SOLUCION APROXIMADA
+%calculo del numero de intervalos
 N=(tN-t0)/h;
-%Creamos los vectores x e y, solución.
+%preparacion del vector solucion aproximada
-%concentración A y B, respectivamente.
 y=zeros(1,length(t));
 x=y;
-%Valor de arranque de los vectores.
+%inicializacion
 y(1)=y0;
 x(1)=x0;
+for i=1:N
 for i=1:N
-    %Euler:
      y(i+1)=y(i)+h*(x(i)*y(i));
      x(i+1)=x(i)+h*(-x(i)*y(i));
 end
-%Dibujamos ambas gráficas.
+%Graficas
-plot(t,y,'r');
 hold on
-plot(t,x);
+plot(t,y,'b');
 hold off
-title('Gráfica Concentración-Tiempo');
+xlabel('segundos (s)');
-xlabel('tiempo (segundos)');
+ylabel('Concentracion (mol/L)');
-ylabel('concentración (mol/l)');
+legend('Concentracion de B','Concentracion de A','Location','Best');
-legend('Concentración de B','Concentración de A');
 }}

@@ Línea 290: / Línea 290: @@
 [[Archivo:GRAFICA_BUENA2.png|800px|centro|thumb|Método Numérico de Euler en Sistema]]
 {{matlab|codigo=
@@ Línea 334: / Línea 336: @@
 [[Archivo:GRAFICA_BUENA2.png|800px|centro|thumb|Método Numérico de Runge-Kutta en Sistema]]
 == Segunda reacción propuesta: Reacción consecutiva propuesta por Lotka ==

← Revisión anterior		Revisión del 22:24 5 mar 2015
Línea 290:		Línea 290:

	[[Archivo:GRAFICA_BUENA2.png\|800px\|centro\|thumb\|Método Numérico de Euler en Sistema]]		[[Archivo:GRAFICA_BUENA2.png\|800px\|centro\|thumb\|Método Numérico de Euler en Sistema]]
		+
		+	{{matlab\|codigo=
		+	%METODO DE RUNGE-KUTTA PARA SISTEMA
		+	t0=0; tN=10; y0=0.01; x0=1; h=0.1;
		+	%generacion del vector tiempo, conocido h
		+	t=t0:h:tN;
		+	%calculo del numero de intervalos
		+	N=(tN-t0)/h; %calculo del numero de intervalos
		+	%preparacion del vector solucion aproximada
		+	y=zeros(1,length(t));
		+	x=y;
		+	%inicializacion
		+	x(1)=x0;
		+	y(1)=y0;
		+	%preparamos las funciones a usar en el bucle
		+	U=inline('x*y','t','y','x');
		+	V=inline('-x*y','t','y','x');
		+	%Recorremos, con un bucle, los vectores, y en ellos almacenamos los
		+	%resultados obtenidos con RK4.
		+	for k=1:N
		+	%Runge-Kutta de orden 4 (RK4):
		+	k1_y=U(t(k),y(k),x(k));
		+	k1_x=V(t(k),y(k),x(k));
		+	k2_y=U(t(k)+(h/2),y(k)+hk1_y/2,x(k)+hk1_x/2);
		+	k2_x=V(t(k)+(h/2),y(k)+hk1_y/2,x(k)+hk1_x/2);
		+	k3_y=U(t(k)+(h/2),y(k)+hk2_y/2,x(k)+hk2_x/2);
		+	k3_x=V(t(k)+(h/2),y(k)+hk2_y/2,x(k)+hk2_x/2);
		+	k4_y=U(t(k)+h,y(k)+k3_yh,x(k)+k3_xh);
		+	k4_x=V(t(k)+h,y(k)+k3_yh,x(k)+k3_xh);
		+
		+	y(k+1)=y(k)+(h/6)(k1_y+2k2_y+2*k3_y+k4_y);
		+	x(k+1)=x(k)+(h/6)(k1_x+2k2_x+2*k3_x+k4_x);
		+	end
		+	%Graficas
		+	hold on
		+	plot(t,y,'b');
		+	plot(t,x,'r');
		+	hold off
		+	xlabel('segundos (s)');
		+	ylabel('Concentracion (mol/L)');
		+	legend('Concentracion de B','Concentracion de A','Location','Best');
		+	}}
		+
		+	[[Archivo:GRAFICA_BUENA2.png\|800px\|centro\|thumb\|Método Numérico de Runge-Kutta en Sistema]]

	== Segunda reacción propuesta: Reacción consecutiva propuesta por Lotka ==		== Segunda reacción propuesta: Reacción consecutiva propuesta por Lotka ==