Matlab - 1D inverse modelling

   
Vista:

1D inverse modelling

Publicado por LauPra (6 intervenciones) el 21/01/2013 12:01:25
Hola,

he leído un estudio en el que realizan 1D inverse modelling mediante Levenberg-Marquardt y quisiera hacer lo mismo con mis datos. Para ello tengo una ecuación con dos parámetros y una incógnita, y una condición ¿puedo resolver la incógnita y los parámetros usando matlab?

Gracias,
LauPra.
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

1D inverse modelling

Publicado por JOSE JEREMIAS CABALLERO jjcc94@hotmail.com (3431 intervenciones) el 21/01/2013 18:03:14
De que se trata el tema que comentas "estudio en el que realizan 1D inverse modelling mediante Levenberg-Marquardt ", pues no conozco ese tema, por eso no te había contestado ni dado sugerencia.

Quizás necesitas de una asesoría mas personalizada, respecto a tu pregunta, pero para ellos debes conocer bien tu fundamento teórico de tu tema y luego buscar un programador en matlab, para que ayude a hacer el código.

Saludos.
JOSE JEREMIAS CABALLERO
Asesor de Proyectos con Matlab
programador en matlab
Servicios de programación matlab
jjcc94@hotmail.com

http://matlabcaballero.blogspot.com
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

1D inverse modelling

Publicado por LauPra (6 intervenciones) el 22/01/2013 10:20:32
Hola José Jeremías,

en primer lugar darte las gracias por tu respuesta y por tu interés en ayudarme a mí y a otros como yo, según veo por tu elevado número de intervenciones.
En las últimas semanas he preparado las ecuaciones necesarias para hacer el modelado inverso con mis datos, sin embargo una vez en este punto no sé resolverlas en matlab.

Planteo el problema a continuación:

taprp_s11=[24 29 3];
taprp_s21=[12 12 4];
tahcp_s1=[22 21 7];
tahcp_s2=[23 22 8];

Rprp_s11=(taprp_s11-ttop)/(tsub-ttop);
Rprp_s21=(taprp_s21-ttop)/(tsub-ttop);
Rhcp_s1=(tahcp_s1-ttop)/(tsub-ttop);
Rhcp_s2=(tahcp_s2-ttop)/(tsub-ttop);


% PARAMETROS tsub, ttop

% INCOGNITA Z*sub

s11=1.1;
s21=2.1;
s1=1;
s2=2;

Za=0.075;

Z*sub = s11 * (Rprp_s11/(2*(1-(Rprp_s11)^2)^0.5) - Za;
Z*sub = s21 * (Rprp_s21/(2*(1-(Rprp_s21)^2)^0.5) - Za;
Z*sub = s1 * (((1/(4*Rhcp_s1)^2) - 0.25)^0.5) - Za;
Z*sub = s2 * (((1/(4*Rhcp_s2)^2) - 0.25)^0.5) - Za;

% CONDICION minimizar un sumatorio
Zsub=[0.15 0.10 0.20]; % medido
sumatorio desde i=1 hasta i=n de (Zsub(i)-Z*sub(i))^2 = mínimo; % no se plantearlo en matlab


¿Se entiende bien?

Un saludo,
LauPra.
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