non linear minimization

1 次查看(过去 30 天)
gianluca
gianluca 2012-3-9
Hi, I've this function which describes the temperature T vs depth z
T(i) = T0 - (H(i).*z(i).^2)./k(i) + ((Q.*z(i))./k(i))
where H(i) and k(i) are known vectors (n x 1) of the same size of z(i) which describe the heat source and thermal conductivity vs depth, respectively. I know temperature measurements at some depth
T = [0 18;1000 52; 1154 55]
I would solve this function for Q (surface terrestrial heat flow). Probably the commands lsqnonlin or fsolve could help me but I do not know how to do. Thanks for the help.
Gianluca

回答(1 个)

Andrew Newell
Andrew Newell 2012-3-9
The simplest approach would be to invert the equation to get
Q = ((T-T0).*k +(H.*z.^2))./z;
This will give you three estimates of Q. Take the mean of them.

类别

Help CenterFile Exchange 中查找有关 Thermal Analysis 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by