2021-05-11 Matlab遗传算法工具箱的使用及实例(非线性规划)

悠悠 2023-01-20 08:08 19阅读 0赞

## Matlab遗传算法工具箱的使用及实例(非线性规划) ##

本文将介绍MATLAB遗传算法工具箱求解非线性规划问题。在阅读本文之前，建议读者阅读上一期“MATLAB遗传算法工具箱求解线性规划问题”。文章传送门：

[Matlab遗传算法工具箱的使用及实例(线性规划)][Matlab]

**一、引例**

上一期介绍了遗传算法求解线性规划的问题。我们来看看下面这个例子，能否用上次讲的方法解决。

![图片][8b0a2cc6f12ec5e0529d2d2049904854.png]

上述例子，第二个约束条件含有二次项，并不是线性的，用上次的方法好像无法直接解决。下面我们就来介绍一下非线性规划的遗传算法的实现。

**二、非线性规划的标准形式**

2.1 非线性规划的标准形式

和线性规划一样，在调用遗传算法工具箱之前，也得学习一下非线性规划的标准形式。因为，调用工具箱需要我们将非标准形式的模型转化为标准形式。

![图片][6fbf1b53c35124f1405a58b5e3ac9bb2.png]

其中，f(x)是目标函数，线性或者非线性都可以。式\[1\]、式\[2\]、式\[5\]同线性规划，为相应维数的矩阵和向量。式\[3\]表示非线性的不等式约束，式\[4\]是非线性等式约束。c(x)和ceq(x)为非线性向量函数（这里可能不好理解，下面结合一个例子解释）

**三、非线性规划模型的实例**

3.1 还是前面出现的模型

1) 编写适应度函数(即目标函数)

目标函数的编写就不多说了，和之前线性规划的编写方法相同

function f = fitnessfun(x)    f = x(1)^2 + x(2)^2 + x(3)^3 + 8;end

2) 编写非线性约束函数

根据式\[1\]、\[2\]，非线性不等式约束的函数为：，

根据式\[3\]、\[4\]，等式约束的函数为：

![20210511210634410.gif][]

非线性约束函数和需要在MATLAB中编写为如下格式（写在一个function里）

function [c, ceq] = 函数名(x)

参数和返回值解释：

x即为自变量的值，为行向量。c即不等式约束的函数值，是列向量（可以存在多个约束）。ceq即等式约束函数的函数值，同样是列向量（可以存在多个约束）

也就是本题返回值要写成这样的形式（这样看的更直观一些）：

![图片][fdbbeaa6b272313b61cc34eb8e3833d9.png]

因此，这题的function应定义为

function \[c,ceq\] = nonlconfun(x)

% 入口参数 x:为自变量的行向量    % 出口参数 c:非线性不等式约束的函数值。由于本题有两个不等式约束    %           因此c(1,1)为第一个不等式约束函数值，c(2,1)为第二个不等式约束函数值    %        ceq:非线性等式约束的函数值，由于本题有两个等式约束，同c一样    %            ceq(1,1)为第一个等式约束函数值，ceq(2,1)为第二个等式约束函数值    c(1,1) = -x(1)^2 + x(2) - x(3)^3;    c(2,1) = x(1) + x(2)^2 + x(3)^3 - 20;    ceq(1,1) = -x(1) - x(2)^2 + 2;    ceq(2,1) = x(2) + 2*x(3)^2 - 3;end

3) 设置遗传算法工具箱的相关参数

这里我们设置种群规模300，迭代次数为800，返回值为options.

options = gaoptimset('PopulationSize', 300, 'Generations', 800);

4) 调用ga工具箱

我们知道，ga工具箱的调用格式为：

[x_best,fval] = ga(fun, nvars, A, b, Aeq, beq, lb, ub, nonlcon, options);

由于本题没有线性约束，因此参数A、b、Aeq、beq都可以设置为空。

对自变量x的约束为：【0，+无穷】，因此，lb = \[0,0,0\]，而ub可以不设置

因此，这部分代码为

options = gaoptimset('PopulationSize', 300, 'Generations', 800); % 遗传算法相关配置fun = @fitnessfun; % 设置适应度函数句柄nonlcon = @nonlconfun; % 设置非线性约束函数句柄nvars = 3; % 自变量个数A = [];  b = [];Aeq = [];  beq = [];lb = [0;0;0];  ub = [];[x_best, fval] = ga(fun, nvars, A, b, Aeq, beq, lb, ub, nonlcon, options);

完整的代码为：

clear;clcoptions = gaoptimset('PopulationSize',200, 'Generations', 600); % 遗传算法相关配置fun = @fitnessfun; % 设置适应度函数句柄nonlcon = @nonlconfun; % 设置非线性约束函数句柄nvars = 3; % 自变量个数A = [];  b = [];Aeq = [];  beq = [];lb = [0;0;0];  ub = [];[x_best, fval] = ga(fun, nvars, A,b,Aeq,beq,lb,ub,nonlcon,options); function [c,ceq] = nonlconfun(x)    c(1,1) = -x(1)^2 + x(2) - x(3)^3;    c(2,1) = x(1) + x(2)^2 + x(3)^3 - 20;    ceq(1,1) = -x(1) - x(2)^2 + 2;    ceq(2,1) = x(2) + 2*x(3)^2 - 3;end function f = fitnessfun(x)    f = x(1)^2 + x(2)^2 + x(3)^3 + 8;end

运行结果为：

fval的值为10.9886即为所求，此时x\_best的值为\[0.9917 1.0041 0.9988\].

因此当时x1 = 0.9917;x2=1.0041;x3=0.9988，函数取到最小值10.9886

（遗传算法具有一定随机性，每次的运行结果有差别，建议多运行几遍程序，找一个最好的结果）

3.2 求下列问题的解

![图片][f1476d8646940a0524b75a1a6ea970cb.png]

1) 编写目标函数

function f = fitnessfun(x)    f = -2 * x(1) - 3 * x(1)^2 - 3 * x(2) - x(2)^2 - x(3);end

2）编写非线性约束函数

式\[1\]、\[2\]、\[3\]为非线性不等式约束，式\[4\]为非线性等式约束，因此

![20210511210634416.gif][]

使用MATLAB编写为：

function [c,ceq] = nonlconfun(x)    % 入口参数 x:为自变量的行向量    % 出口参数 c:非线性不等式约束的函数值。由于本题有多个不等式约束，以列向量的形式返回即可    %         ceq:非线性等式约束的函数值。    c(1,1) = x(1) + 2 * x(1)^2 + x(2) + 2 * x(2)^2 + x(3) - 10;    c(2,1) = x(1) + x(1)^2 + x(2) + x(2)^2 - x(3) - 50;    c(3,1) = 2 * x(1) + x(1)^2 + 2 * x(2) + x(3) - 40;    ceq = x(1)^2 + x(3) - 2;end

3）调用ga函数

（这里就不调用gaoptimset函数配置遗传算法的参数了，直接使用默认参数）

式\[5\]、\[6\]为线性不等式约束，标准形式为Ax<=b

![图片][326ee64ecf6866b2556b50727356474b.png]

可转化为：，因此：

![图片][b9707ae0aa453e5383d261063d391846.png]

本题没有线性等式约束，也没有自变量约束，因此Aeq、beq、lb、ub均设置为空矩阵。写好ga所需的入口参数，即可调用ga函数，调用方式如下：

fun = @fitnessfun; % 设置适应度函数句柄，(在函数名前加@即可)nonlcon = @nonlconfun; % 设置非线性约束函数句柄nvars = 3; % 自变量个数A = [-1,-2,0;-1,0,0];  b = [-1;0];  % 线性不等式约束Aeq = [];  beq = []; % 线性等式约束lb = [];  ub = [];   % 自变量定义域[x_best, fval] = ga(fun, nvars, A, b, Aeq, beq, lb, ub, nonlcon, []);

整体程序为：

clear;clcfun = @fitnessfun; % 设置适应度函数句柄，(在函数名前加@即可)nonlcon = @nonlconfun; % 设置非线性约束函数句柄nvars = 3; % 自变量个数A = [-1,-2,0;-1,0,0];  b = [-1;0];  % 线性不等式约束Aeq = [];  beq = [];   % 线性等式约束lb = [];  ub = [];     % 自变量定义域[x_best, fval] = ga(fun, nvars, A,b,Aeq,beq,lb,ub,nonlcon,[]); function [c,ceq] = nonlconfun(x)    c(1,1) = x(1) + 2 * x(1)^2 + x(2) + 2 * x(2)^2 + x(3) - 10;    c(2,1) = x(1) + x(1)^2 +

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FpbmdmZW5neGQx_size_16_color_FFFFFF_t_70][]

[Matlab]: http://mp.weixin.qq.com/s?__biz=MzU0NTkyMzk2MA==&mid=2247485478&idx=1&sn=32a1300da66bcb07802e37d44aa63191&chksm=fb643e31cc13b727053842b7cc9a16b5efa174b6e014c8c19d83899fd075b120d910305e9c40&scene=21#wechat_redirect
[8b0a2cc6f12ec5e0529d2d2049904854.png]: /images/20221021/164eac0939b74053b803340b49503359.png
[6fbf1b53c35124f1405a58b5e3ac9bb2.png]: /images/20221021/7b6a538fa54a4c59b49e22e345daa9a8.png
[20210511210634410.gif]: /images/20221021/7aa26e6cbb2e46519ed6ff460d490786.png
[fdbbeaa6b272313b61cc34eb8e3833d9.png]: /images/20221021/5448839ce86e45a5ad3255ef1178291c.png
[f1476d8646940a0524b75a1a6ea970cb.png]: /images/20221021/fc341f41538c44d38602b3c141c92d51.png
[20210511210634416.gif]: /images/20221021/5dbf5903fe794cfe924887fe8b2dd830.png
[326ee64ecf6866b2556b50727356474b.png]: /images/20221021/b0bdabb182154283bfd358a9597767a8.png
[b9707ae0aa453e5383d261063d391846.png]: /images/20221021/45739a170d9b49e6aca1df9512cd84a0.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FpbmdmZW5neGQx_size_16_color_FFFFFF_t_70]: /images/20221021/1c9e75ffb958412b9902abb55ab61bca.png