Matlab有GA优化工具箱。
例1
在这里,就用2017年重庆高考数学最后一道题来开刀了。
当然,只是算一个大概,大概……
理论证明嘛,无非就是各种不等式去试……
1. 默认求 min f(x,y)。
2. 非线性约束条件按规范写。
不等式 g(x,y) >= 0;
等式 h(x,y) = 0;
e.g. x^2 + y^2 = 2 要写为 x^2 + y^2 - 2
2元数量值函数。
NonCon
function [c, ceq] = NonCon( x )
c = []; % 非线性不等式约束条件为空
ceq = x(1)^3 + x(2)^3 - 2; % 非线性等式约束条件
end
FitFun1
function y = FitFun1(x)
% x可以是一个向量,FitFun()是一个一(多)元数量值函数
y = (x(1)+x(2))*(x(1)^5+x(2)^5); % 第1问
% y = -(x(1)+x(2)); % 第2问
end
结果如下
----------------------------
Optimization running.
Objective function value: 4.018996707929087
Optimization terminated: average change in the fitness value less than options.FunctionTolerance and constraint violation is less than options.ConstraintTolerance.
当然,这个结果并不是试卷上要让我们证明的答案。一分也拿不到。真是惨。