脚本如下:
fun = @(t)-(183545*t-367490003);
lb = 2004;
ub = 2030;
t0 = 2015;
x = fmincon(fun,t0,[],[],[],[],lb,ub,@c)
起始值t0随便写一个。
另外在c.m中定义非线性约束函数:
function [c,ceq] = c(t)
c = 462.2*t.^2-2e6*t+2e9; % nonlinear inequalities at x.
ceq = 0; % nonlinear equalities at x.
运行结果2030。
追问还想再请问下,我发现算出来的结果和我设的ub值有关,如果我把ub值改成2050年,那算出来最大值就是2050。这个结果是不是不对呢??
追答其实你可以验证一下你的约束函数,先看看其在ub处是不是<0。我验证过2030是的。如果这个条件成立的话,一直取ub就很显然了。