遗传算法

uaidenp · 发表于 2020-7-21 09:34

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

/*arguments参数*/
c1=[16000;500;350;500;2500;1000;0;0;0;0;0;0;0;0]
c2=[3250;1000;350;500;2500;1000;28;27;7;44;0;0;0;0]
mypi=[1600;425;213;213;2600;750;11;11;3;18;40800;1200;3000;8000]
e=[0;0;0;0;0;0;0;0;0;0;10200;300;750;2000]
r=[5;4;3;3;3;3;1;1;1;1;2;0;1;0]

/*fun1*/
function f=Fun1(x,E,mypi)
f=x*(E-mypi)-85*x(3)*x(4)-3400*x(12)*x(13)*x(14);
end

/*main函数*/
ObjectFunction = @Fun1;
nvars = 14;%变量个数
LB = [0]%定义域下限
UB = [1]%定义域上限
A=[-1,0,0,0,0,0,0,0,0,0,0,0,0,0;-1,1,0,0,0,0,0,0,0,0,0,0,0,0;0,-1,1,0,0,0,0,0,0,0,0,0,0,0;0,-1,0,1,0,0,0,0,0,0,0,0,0,0;0,0,0,0,-1,0,1,0,0,0,0,0,0,0;0,0,0,0,-1,0,0,1,0,0,0,0,0,0;0,0,0,0,-1,0,0,0,1,0,0,0,0,0;0,0,0,0,-1,0,0,0,0,1,0,0,0,0;-1,0,0,0,0,0,0,0,0,0,1,0,0,0;0,0,0,0,0,0,0,0,0,0,-1,1,0,0;0,0,0,0,0,0,0,0,0,0,-1,0,1,0;0,0,0,0,0,0,0,0,0,0,-1,0,0,1]
b=[1;0;0;0;0;0;0;0;0;0;0;0]
[x,fval] = ga(ObjectiveFunction,nvars,A,b,[],[],LB,UB,NonCon)%调用ga函数

/*非线性约束NonCon*/
function [c,ceq] = NonCon(x)
%UNTITLED2 此处提供此函数的摘要
% 此处提供详细说明
c=[[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8),x(9),x(10),x(11),x(12),x(13),x(14)]*[16000;500;350;500;2500;1000;0;0;0;0;0;0;0;0]-268*x(3)*x(4)-350*x(4)*x(5)-250*x(4)*x(6)-600*x(4)*x(5)*x(6)-20000;[x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8),x(9),x(10),x(11),x(12),x(13),x(14)]*[3250;1000;350;500;2500;1000;28;27;7;44;0;0;0;0]-155*x(2)*x(3)-255*x(2)*x(4)-188*x(4)*x(3)-200*x(4)*x(5)-175*x(4)*x(6)-125*x(5)*x(6)-375*x(4)*x(5)*x(6)-6000];
ceq=[];
end

这个一直报错输入参数不足，错误在Noncon.m的line4，请教大神