Matlab - Ayuda iteracion matlab

   
Vista:

Ayuda iteracion matlab

Publicado por Alan reinoso.alan@hotmail.com (1 intervención) el 28/02/2018 14:11:11
Buenos días, necesito un programa en Matlab con el cual mediante iterciones pueda hallar el valor de x para el cual se cumpla la igualdad de los términos A y B, es decir x para que A=B.

1
2
3
4
5
6
7
8
9
10
11
12
13
T=500;
Fa0=2000
Fb0=1000
ya=0.16
thetab=Fb0/Fa0;
r=10;
k1=(1.79E4*(exp(-20800/(1.987*T)
k2=(2.57E16*(exp(-47400/(1.987*T)))
Pa=(ya*(1-x))/(1-2*ya*x)
Pb=(ya*(thetab-3*x))/(1-2*ya*x)
Pc=(2*ya*x)/(1-2*ya*x);
A=(k1)*((Pa*((pb)^1.5)*(200^1.5))/Pc));
B=r+(k2)*(Pc/((pb)^1.5)*(200^1.5)));
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder
Imágen de perfil de JOSE JEREMIAS CABALLERO

Ayuda iteracion matlab

Publicado por JOSE JEREMIAS CABALLERO jjcc94@hotmail.com (4632 intervenciones) el 28/02/2018 15:07:05
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clear all
T=500;
Fa0=2000;
Fb0=1000;
ya=0.16;
thetab=Fb0/Fa0;
syms x
r=10;
k1=(1.79E4*(exp(-20800/(1.987*T))));
k2=(2.57E16*(exp(-47400/(1.987*T))));
Pa=(ya*(1-x))/(1-2*ya*x);
Pb=(ya*(thetab-3*x))/(1-2*ya*x);
Pc=(2*ya*x)/(1-2*ya*x);
A=k1*((Pa*(Pb)^1.5*(200^1.5))/Pc);
B=r+k2*(Pc/((Pb)^1.5)*(200^1.5));
D=vpa(A-B,3);  %% aca puede aplicar un método iterativo para resolver por iteraciones.
x=vpa(solve(D,'x'),5)
%d=inline(D);
%Dx=d(x)


1
2
3
4
5
>> iteraciones
 
x =
 
0.000046279


Saludos
JOSE JEREMIAS CABALLERO
Asesor de Proyectos con Matlab
Servicios de programación matlab


http://matlabcaballero.blogspot.com
https://www.facebook.com/matlabcaballero
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar
Revisar política de publicidad