Fitting a function to data (fminsearch) with limits

165 次查看(过去 30 天)
Hey!
I try to fit a model to measured data. I already got a real close with the fminsearch option. However, I know the model will never fit completly to the measured data. I know the range wherein the unknown factors should be, how can I include this in the function.
This is a part of my code;
p= [0.31;114;3.5];
y2=fitfunc(TempZone1day,TempZone2day,TempZone3day,TempZone4day,TempZone5day,buitenTempday,FlowZone1day,FlowZone2day,FlowZone3day,FlowZone4day,KNMIwindday,SMA,p);
figure(1)
plot(tijd,WarmteCvday,'+',tijd,y2,'-')
a0=p;
aBest = fminsearch(@(a) SumErrfun1(a,TempZone1day,TempZone2day,TempZone3day,TempZone4day,TempZone5day,buitenTempday,FlowZone1day,FlowZone2day,FlowZone3day,FlowZone4day,KNMIwindday,SMA,WarmteCvday,tijd),a0);
disp(aBest)
For example, I know that p(1) should be within 0.1-0.5, p(2) within 100-200 and p(3) within 2-6. Because now aBest (the best combination) is [-0.0328; 61.8202; 0.4375], which is not even a possible option. How can I include these ranges?

采纳的回答

John D'Errico
John D'Errico 2016-8-30
编辑:John D'Errico 2016-8-30
fminsearch has no capability to take bounds on the search. If the objective is such that a better result lies outside of where you want it, too bad. :)
Having said that, you can use fminsearchbnd , a tool found on the file exchange. It does allow bounds on the variables. Just download and install that tool on your search path, then use it instead.

更多回答(3 个)

Jie Jian
Jie Jian 2020-1-9
Or you can use the function 'mapping_parameters.m' to transfer unbounded parameters to bounded ones

kursat cihan
kursat cihan 2020-8-30
John D'Errico....much love from Germany, helped me a lot!!!
Bachelor Thesis in Material Modelling, used it for a parameter optimization in bringing simulations together with experimental data...PEACE

Stefan Schuberth
Stefan Schuberth 2022-11-8
you can use q=f1*atan(p)+f2 to construct a limited parameter q from an unlimited parameter p :)

类别

Help CenterFile Exchange 中查找有关 Nonlinear Optimization 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by