I have a system of 2 ODE's (ode45) and want to fit them to measured data. For that i have 4 variables i want to optimate.
This is my code:
% ode45
tspan=linspace(0,3785,3785);ic=[300 300];x=ones(1,4);f = @(t,T) [(x(1)*(T(2)-T(1))+x(2))/100; x(3)*(T(1)-T(2))+x(4)*(300-T(2))];[t,T]= ode45(f,tspan,ic);% ydata from excel
filename='Messung_1.xlsx';sheet =12;Range1='C6:C3790';Range2='D6:C3790';t_Kl=xlsread(filename,sheet,Range1)+273.15;t_Nl=xlsread(filename,sheet,Range2)+273.15;ydata=[t_Kl t_Nl];% objective function
objective= @(x) sum(((T(:,1)-ydata(:,1))./ydata(:,1)).^2)+sum(((T(:,2)-ydata(:,2))./ydata(:,2)).^2); x0=ones(1,4);pbest=fmincon(objective,x0);
When i try to solve this: Initial point is a local minimum that satisfies the constraints.
Can someone help me please?
Best Answer