可以用fmincon。
t=[1 6:6:72];
c=[0.093 0.064 0.058 0.055 0.053 0.05 0.044 0.043 0.042 0.04 0.039 0.039 0.038];
function [c,ceq] = mycon1104(x)
c = 1-x(2)-x(3)-x(2)*x(3);
ceq = [];
end %约束函数,作为m文件
x = fmincon(@(x) fun(x),rand(1,3),[],[],[],[],[],[],@(x) mycon1104(x))
得到x,可能是指数函数的原因,比较难找到最优值,或者可能在很多点都可以达到。
x =0.0754 0.0121 7.6007
fun(x)=5.8530e-04
x = -0.0754 65.6695 0.0122
fun(x)= 5.8390e-04
-------