Using hooke-jeeves method with ode45

Hello, I got an ode45 function and I really need to reverse kinematics. I cant transform the formulas so natural way is to do it by optimalization. Ive heard that Hooke-Jeeves should be helpful but to be honest I cant really combain it ;/ Can anyone tell me if it is possible to do it and if yes help me with reasonable sources? I would be very thankful.

[Jan]

The question has only a weak connection to Matlab. Perhaps somebody has an idea, but the problem is off-topic.

[Jakub]

Why is that? Im coding it all in matlab and I guess ode45 is matlab only function?

Answer 1

[Star Strider]

See if the File Exchange contribution Unconstrained optimization using Hooke & Jeeves does what you want.

[Jakub]

It looks cool. Can anyone help me with it tho?

[t,x]=ode45(@discfltEOM,tspan,x0,options,CoefUsed);
&x0 are input parameters in a vector

%Hooke-Jeeves
%   Sx      : objective function
%   guess    : initial point
%   par      : parameters needed to function
%   ip      : (0): no plot (default), (>0) plot figure ip with pause, (<0) plot figure ip
%   Lb, Ub  : lower and upper bound vectors to plot (default = guess*(1+/-2))
%   problem  : (-1): minimum (default), (1): maximum
%   tol      : tolerance (default = 1e-4)
%   mxit    : maximum number of stages (default = 50*(1+4*~(ip>0)))
%   stp      : stepsize vector for the independent variables (default = max(0.01*abs(guess+~guess),0.1))
%   amp      : stepsize enlargement factor (1,oo) (default = 1.5)
%   red      : stepsize reduction factor (0,1) (default = 0.5)
%   xopt    : optimal point
%   Opt      : optimal value of Sx
%   Nav      : number of objective function evaluations
[xopt,Opt,Nav]=hkjeeves(Sx,guess,ip,Lb,Ub,problem,tol,mxit,stp,amp,red,par1, par2,...)

I dont really know who to combain those fuctions ;/ Yeah I realize Im newbie with MatLab.

[Star Strider]

Kinematics are not in my areas of expertise. I just found it on the File Exchange for you.