How to set up fmincon

11 次查看（过去 30 天）
sam plant
sam plant 2018-12-13
重新打开： 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更早的评论
sam plant
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

请先登录，再进行评论。

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

采纳的回答

Torsten
Torsten 2018-12-13
编辑：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 个）

Alan Weiss
Alan Weiss 2018-12-13
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

类别

在 Help Center 和 File Exchange 中查找有关 Nonlinear Optimization 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by