How to set up fmincon

Question

sam plant 2018-12-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/435557-how-to-set-up-fmincon

重新打开： madhan ravi 2018-12-13

I have an objective function which takes input values for a projectile with initial speed v0 and angle theta. Given a target it then calculates the minimum distance of the trajectory to the target and outputs the minimum distance d.

Of course i need to use fmincon to get the minimum distance of the objective function less than 0.1 so that the projectile hits the target but I'm having trouble understanding how to set fmincon up with inputs.

Any Help would be appreciated!

4 个评论
显示 2更早的评论隐藏 2更早的评论

madhan ravi 2018-12-13

solve_ode45??

sam plant 2018-12-13

that's another function i had to write using ode45, i have to write 3 functions inside of each other

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2018-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/435557-how-to-set-up-fmincon#answer_352175

编辑：Torsten 2018-12-13

function main
v0=...;   % start value for initial velocity
theta0=...; % start value for angle
target=[... , ...]; % target position
x0 = [v0, theta0];
x = fminsearch(@(x)objective(x,target),x0)
end
function M=objective(initial,target)
v0=initial(1);
theta=initial(2);
h=0;
dt=0.01;
maxtime=60;
range=solve_ode45(v0,theta,h,dt,maxtime); %calculates projectiles path
compare = [target(1)-range(:,3),target(2)-range(:,4)]; %compares displacement of target to line at all points
d=sqrt(compare(:,1).^2+compare(:,2).^2); %magnitude of displacement to give distance at all points
M=min(d);   %finds the minimum distance
end 

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

sam plant 2018-12-13

thank you so much! it's worked!

请先登录，再进行评论。

Answer 2

Alan Weiss 2018-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/435557-how-to-set-up-fmincon#answer_352172

Take a look at Optimize an ODE in Parallel. You probably need to specify a time at which the projectile hits the closest to the target. In the example, we use fzero to find when the projectile hits the ground, and you will probably have to define a similar stopping criterion and use fzero also. Or maybe use an interpolation of the trajectory to find the minimum distance. You can easily interpolate an ODE45 solution by using deval.

Alan Weiss

MATLAB mathematical toolbox documentation