Informatica Aplicada A La Ing. Quimica - Ucm

INFORMATICA APLICADA A LA INGENIERÍA QUIMICA: 2ª parte MATLAB
2.1 Resolución de un sistema de ecuaciones algebraicas
•Se plantea resolver el siguiente problema: sistema de n ecuaciones

algebraicas con n incógnitas, que puede incluir ecuaciones no lineales,
con la siguiente forma general.
f1(x1,…,xn)=0
f2(x1,…,xn)=0
.
.
.
fn((x1,,…,x
, n))=0
Encontrar el vector (x1,…….,xn) que hace que todas las ecuaciones sean igual a
cero Matlab tiene una sentencia para resolver esto: fsolve.m.
cero. fsolve m Utilizamos una
forma mejorada de este que es csolve.m

• Sintáxis:
[x,rc]=csolve(‘nombrefunción’,x0,’gradfun’,crit,itmax,p1,p2,…)
Nombre función = es un fichero función que calcula los valores de fi.
•Sintáxis de ‘nombrefunción’:
function y=nombrefuncion(x)
y(1,:)=f1(x(1,:), x(2,:),…,x(n,:));
y(n :)=fn(x(1,:),…,x(n,:));
y(n,:)=f (x(1 :) x(n :));
La variable x es una matriz con n filas (n = nº de incógnitas) y con un
número variable de columnas. Cada columna es una iteración que hace
ell programa.
X10 X11 X12 ------ X1n y10 (X10 X20‐‐‐‐‐ Xn0) y11 (X11 X21‐‐‐‐‐ Xn1) y1n (X1n X2n‐‐‐‐‐ Xnn)
3
‘ ‘ ‘ ‘ ‘ ‘
‘
‘ ‘ ‘ ‘ ‘ ‘
‘
‘ ‘ ‘ ‘ ‘ ‘
‘
‘ ‘ ‘ ‘ ‘ ‘
‘
Xn0 Xn1 Xn2 - - - - - - Xnn yn0 (X10 X20‐‐‐‐‐ Xn0) yn1 (X11 X21‐‐‐‐‐ Xn1) ynn (X1n X2n‐‐‐‐‐ Xnn)
Cuando se logre que sustituyendo la

columna de valores de x (desde x1
hasta xn) la columna de valores de las
funciones se aproximen todas a cero,
habré llegado a la solución.
• x0: es un vector columna con la primera estimación de x (x10,….,xn0)
• ‘gradfun’: función que calcula el gradiente de ‘nombrefunción’. Sirve

para acelerar la llegada
p g del comando a la solución del sistema de
ecuaciones. No es absolutamente necesario emplearlo ( en vez de
escribir grandfun, se escribe []).
i n
• crit: número pequeño positivo que cumple y

i1
n  crit en la solución
(y de esta manera se determina que se ha llegado a la solución).
• itmax: el número máximo de iteraciones permitidas al comando.
• p1,p2,…: parámetros adicionales que se introducen en esta posición

para calcular y1,…,yn. (Se pueden introducir con una sentencia global y
no sería necesario introducirlos aquí).
aquí)
[x rc]=csolve(‘nombrefunción’
[x,rc]=csolve( nombrefunción ,x0,
x0 ’gradfun’
gradfun ,crit,itmax,p1,p2,…)
crit itmax p1 p2 )
•rc: indicador del grado de acercamiento a una solución verdadera.

rc = 0, se ha llegado a una solución verdadera
rc = 4, se ha llegado al máximo número de interacciones.
•x: vector columna con la solución del sistema de ecuaciones si rc=0
El cálculo de la matriz [y(1,:); y(2,:);…….; y(n,:)] es muy complicado

cuando las ecuaciones f1,……….,ffn son complejas…….
complejas
DIFICULTAD EN LAS OPERACIONES ELEMENTO A ELEMENTO.

El cálculo
ál l ded la
l matriz
t i [y(1,:);
[ (1 ) y(2,:);…….;
(2 ) y(n,:)]
( )] es muy complicado
li d
cuando las ecuaciones f1,……….,fn son complejas…….
DIFICULTAD EN LAS OPERACIONES ELEMENTO A ELEMENTO.
‘ ‘ ‘ ‘ ‘ ‘ ‘
‘ ‘ ‘ ‘ ‘ ‘ ‘
‘ ‘ ‘ ‘ ‘ ‘ ‘
‘ ‘ ‘ ‘ ‘ ‘ ‘
Xn0 Xn1 Xn2 - - - - - - Xnn yn0 (X10 X20‐‐‐‐‐ Xn0) yn1 (X11 X21‐‐‐‐‐ Xn1) ynn (X1n X2n‐‐‐‐‐ Xnn)
• Para facilitar el cálculo de las funciones f1,…,fn a partir de la

matriz x(1,:),…,x(n,:),
x(1 :) x(n :) se extrae cada columna x1,…x xn por
separado, y se calcula:
y1=f1(x1,…,xn)
y2=ff2(x
( 1,…,xn)
.
.
.
yn=fn(x1,…,xn)
• Construimos la función ‘nombrefunción’ :

filas=size(x 1) %número de filas de la matriz x
filas=size(x,1)
columnas=size(x,2) %número de columnas de la matriz x
El comandod size(matriz,
i ( t i número)
ú ) calcula
l l ell número
ú d filas
de fil de
d lal
matriz cuando el número es 1, y el número de columnas cuando es 2.
• Se busca calcular una matriz y a partir de una matriz x con la siguiente

forma:
 x1,0 ꞏ ꞏ ꞏ x1,columnas   y1,0 ꞏ ꞏ ꞏ y1,columnas 
 ꞏ ꞏ   ꞏ 
  ꞏ
 
 ꞏ ꞏ   ꞏ ꞏ 
   
 ꞏ ꞏ   ꞏ ꞏ 
 x filas,0 ꞏ ꞏ ꞏ x filas,columnas   y filas,0 ꞏ ꞏ ꞏ y filas,columnas 
• El siguiente paso es extraer una por una cada columna de la matriz x:
ffor jj=1:columnas
1 l
for i=1:filas
v(i) =x(i,j);
end
2 1 Resolución
2.1 R l ió de d un sistema
i de
d ecuaciones
i algebraicas
l b i
•(v
( 1,,…,v
, n) contiene ((x1,,..,x
, n) para
p calcular (y1,,…,y
,yn))’
Con este ciclo, el vector (v1,…,vn) contiene el conjunto de valores con el

que hay que calcular (y1,…,yn).
y(1,j)=f
y(1 j)=f1(v(1),…,v(n));
(v(1) v(n));
.
.
.
y(n,j)=fn(v(1),…,v(n));
%Se ha calculado la columna j de la matriz deseada
end %vuelve a calcular la siguiente columna de y

principal.m
n=lenght(x); %número de ecuaciones (o de incógnitas)

x0=ones(1,n) %Suponemos x(1)=1, x(2)=1, …….x(n)=1
x0=x0’ %Necesario
i para csolve
l
[
[x,rc]=csolve(‘nombrefuncion’,x0,[
] ( [ ],1e-7,10000)
] ) funcion.m
x %vector solución
rc
( , );
filas=size(x,1);
columnas=size(x,2);
for j=1:columnas
for i=1:filas
v(i)=x(i,j);
end
y(1,:)=f1(x(1,:), x(2,:),…,x(n,:));
LAS FUNCIONES
.
DE f(x1, x2….xn)
.
y(n,:)=fn(x(1,:),…,x(n,:));
end
• Ejemplo: resolver el siguiente sistema de ecuaciones:
3  x12 x 2  x 2  7x 3  0 funcion.m
x1  x 33  0 function y=funcion(x)
x2  x3  0 filas=size(x,1);
filas size(x,1);
columnas=size(x,2);
for j=1:columnas
f ii=1:filas
for 1 fil
v(i)=x(i,j);
end
y(1,j)=3+v(1)^2*v(2)+v(2)-7*v(3);
y(2,j)=v(1)-v(3)^3;
y(3 j)=v(2) v(3);
y(3,j)=v(2)-v(3);
end
principal.m
n 3; %número de ecuaciones (o de incógnitas)

n=3;
x0=ones(1,n) %Suponemos x(1)=1, x(2)=1, x(3)=1
x0=x0’ %Necesario para csolve
[x,rc]=csolve(‘funcion’,x0,[ ],1e-7,10000)
x %vector solución
rc
PROGRAMA
n=3
x0=ones(1,n)
x0=x0'
0 0'
[x,rc]=csolve('funcionparacsolve',x0,[],1e-7,10000)
FUNCION
csolve
FUNCION
function y=funcionparacsolve(x)
filas=size(x,1)
columnas=size(x 2)
columnas=size(x,2)
for j=1:columnas
for i=1:filas
v(i)=x(i,j)
end
y(1,j)=3+v(1)^2*v(2)+v(2)-7*v(3);
y(2,j)=v(1)-v(3)^3;
y(3 j)=v(2)-v(3);
y(3,j) v(2) v(3);
end
PROGRAMA
n=3
x0=ones(1,n)
x0=x0'
0 0'
[x,rc]=csolve('funcionparacsolve',x0,[],1e-7,10000)
FUNCION
csolve
FUNCION
function y=funcionparacsolve(x)
filas=size(x,1)
columnas=size(x 2)
columnas=size(x,2)
for j=1:columnas
for i=1:filas
v(i)=x(i,j)
end
y(1,j)=3+v(1)^2*v(2)+v(2)-7*v(3);
y(2,j)=v(1)-v(3)^3;
y(3 j)=v(2)-v(3);
y(3,j) v(2) v(3);
end

Informatica Aplicada A La Ing. Quimica - Ucm

Cargado por

Copyright:

Formatos disponibles

Informatica Aplicada A La Ing. Quimica - Ucm

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Informatica Aplicada A La Ing. Quimica - Ucm

Cargado por

Copyright:

Formatos disponibles

INFORMATICA APLICADA A LA INGENIERÍA QUIMICA: 2ª parte MATLAB

2.1 Resolución de un sistema de ecuaciones algebraicas

•Se plantea resolver el siguiente problema: sistema de n ecuaciones

2.1 Resolución de un sistema de ecuaciones algebraicas

Nombre función = es un fichero función que calcula los valores de fi.

Cuando se logre que sustituyendo la

2.1 Resolución de un sistema de ecuaciones algebraicas

• x0: es un vector columna con la primera estimación de x (x10,….,xn0)

• ‘gradfun’: función que calcula el gradiente de ‘nombrefunción’. Sirve

• crit: número pequeño positivo que cumple y

(y de esta manera se determina que se ha llegado a la solución).

• itmax: el número máximo de iteraciones permitidas al comando.

• p1,p2,…: parámetros adicionales que se introducen en esta posición

2.1 Resolución de un sistema de ecuaciones algebraicas

•rc: indicador del grado de acercamiento a una solución verdadera.

•x: vector columna con la solución del sistema de ecuaciones si rc=0

El cálculo de la matriz [y(1,:); y(2,:);…….; y(n,:)] es muy complicado

DIFICULTAD EN LAS OPERACIONES ELEMENTO A ELEMENTO.

DIFICULTAD EN LAS OPERACIONES ELEMENTO A ELEMENTO.

• Para facilitar el cálculo de las funciones f1,…,fn a partir de la

• Construimos la función ‘nombrefunción’ :

• Se busca calcular una matriz y a partir de una matriz x con la siguiente

• El siguiente paso es extraer una por una cada columna de la matriz x:

Con este ciclo, el vector (v1,…,vn) contiene el conjunto de valores con el

%Se ha calculado la columna j de la matriz deseada

end %vuelve a calcular la siguiente columna de y

n=lenght(x); %número de ecuaciones (o de incógnitas)

2.1 Resolución de un sistema de ecuaciones algebraicas

• Ejemplo: resolver el siguiente sistema de ecuaciones:

2.1 Resolución de un sistema de ecuaciones algebraicas

n 3; %número de ecuaciones (o de incógnitas)

También podría gustarte