ayuda con funcion objetivo
Publicado por principiante (3 intervenciones) el 21/02/2011 15:33:43
Hola a todos,
Necesito realizar un ajuste utilizando Matlab. Mis datos experimentales deben de ajustarse a una ecuación del tipo:
Y=1-exp(-(k*t)^c1)
donde;
k=c2*exp(-c3/T)
Y es la variable dependiente.
T y t son las variables independientes.
c1,c2 y c3 son los coeficientes que deseo calcular.
Anteriormente había escrito y me proporcionaron la siguiente alternativa para resolver mi problema:
Primero, definir la función objetivo, que es que la diferencia entre las Y original y encontrada sea cero.
Si esta función se llama OF_exp quedaría así:
------------------------------------
function U = OF_exp(C)
global Y T t
k = C(2)* exp(-C(3)./ T);
Y2 = 1 - exp(-(k .* t) .^ C(1));
U = norm(Y - Y2, 1);
-----------------------------------
Luego construyo la función que utiliza a fminsearch para minimizar la función objetivo.
Mi consulta es la siguiente:
Es posible dejar el ajuste en base al error absoluto?
Como debería escribir la última línea de la función OF_exp(C) de modo que:
U= sumatoria ((Ycalc-Yexp)/Yexp)^2
Agradezco de antemano cualquier sugerencia
Necesito realizar un ajuste utilizando Matlab. Mis datos experimentales deben de ajustarse a una ecuación del tipo:
Y=1-exp(-(k*t)^c1)
donde;
k=c2*exp(-c3/T)
Y es la variable dependiente.
T y t son las variables independientes.
c1,c2 y c3 son los coeficientes que deseo calcular.
Anteriormente había escrito y me proporcionaron la siguiente alternativa para resolver mi problema:
Primero, definir la función objetivo, que es que la diferencia entre las Y original y encontrada sea cero.
Si esta función se llama OF_exp quedaría así:
------------------------------------
function U = OF_exp(C)
global Y T t
k = C(2)* exp(-C(3)./ T);
Y2 = 1 - exp(-(k .* t) .^ C(1));
U = norm(Y - Y2, 1);
-----------------------------------
Luego construyo la función que utiliza a fminsearch para minimizar la función objetivo.
Mi consulta es la siguiente:
Es posible dejar el ajuste en base al error absoluto?
Como debería escribir la última línea de la función OF_exp(C) de modo que:
U= sumatoria ((Ycalc-Yexp)/Yexp)^2
Agradezco de antemano cualquier sugerencia
Valora esta pregunta
0
