24小时热门版块排行榜    

查看: 3442  |  回复: 0

F先森

新虫 (初入文坛)

[求助] 带约束的多目标最优计算问题,gamultiobj 函数的参数设定?

目标函数的M文件:
function f=my_first_multi
f(1)=3*x(1)+9.7*x(2)+7.4*x(3)+6.8*x(4)+11*x(5)+3.5*x(6)+43.58*x(7)+0.26*x(8)+0.76*x(9)+0.18*x(10)+0.55*x(11);
f(2)=24.38*x(1)+44.55*x(2)+61.84*x(3)+28.84*x(4)+164.7*x(5)+64.77*x(6)+134.68*x(7)+55.68*x(8)+290.15*x(9)+43.96*x(10)+61.55*x(11);

调用gamultiobj 的基本程序和参数设定:
>> clear
>> fitnessfcn=@my_first_multi;
>> nvars=11;
>> lb=[24,0,0,0,0,0,0,0,0,0,0];
>> ub=[36,inf,inf,inf,inf,inf,inf,inf,inf,inf,inf];
>> A=[0.0323,0.004,0.06,0.0011,0.0001,0.035,0.0001,0,0,0,0;-0.0323,-0.004,-0.06,-0.0011,-0.0001,-0.035,-0.0001,0,0,0,0;0.006,0.0036,0.03,0.0022,0,0.02,0,0.0093,0,0.014,0;-0.006,-0.0036,-0.03,-0.0022,0,-0.02,0,-0.0093,-0,-0.014,0;0.006,0.0085,0,0.0052,0,0.005,0,0,0,0,0;-0.006,-0.0085,0,-0.0052,0,-0.005,0,0,0,0,0;0.00007,0.0003,0.0003,0.0002,0,0.0004,0,0,0,0,0;-0.00007,-0.0003,-0.0003,-0.0002,0,-0.0004,0,0,0,0,0;0.000035,0.001,0.0004,0.0003,0,0.0006,0,0.003,0,0.001,0;-0.000035,-0.001,-0.0004,-0.0003,0,-0.0006,0,-0.0003,0,-0.001,0;0,0.18,0.6575,0.002,0,0.02,0,0,0,0,0;0,-0.18,-0.6575,-0.002,0,-0.02,0,0,0,0,0;0,0.08,0,0.001,0.9998,0.041,0,0,0,0,0;0,-0.08,0,-0.001,-0.9998,-0.041,0,0,0,0,0;0.0325,0.01925,0.16,0.01175,0,0.10715,0,-0.94,0,-0.61,0;-0.039,-0.0231,-0.192,-0.0141,0,-0.12858,0,0.83,0,0.58,0;0.00035,0.001,0.0004,0.0003,0,0.0006,0.79,0.0003,0,0.001,0;-0.00035,-0.001,-0.0004,-0.0003,0,-0.00006,-0.79,-0.0003,0,-0.001,0;];
>> b=[3.537;-1.895;0.306;0;1.6;0;0.04;0;0.387;0;24.783;-22.174;7.347;-6.122;0;0;0.36;0;];
>> Aeq=[0.455,0.303,0.237,0.407,0.0016,0.08,0.0016,0.00985,0,0.0155,-1];
>> beq=[4.04];
>> options=gaoptimset('ParetoFraction',0.3,'PopulationSize',100,'Generations',200,'StallGenLimit',200,'TolFun',1e-100,'PlotFcns',@gaplotpareto);
>> [x,fval]=gamultiobj(fitnessfcn,nvars,A,b,Aeq,beq,lb,ub,options);


以上是一个双目标的约束优化计算问题,运行结果是  Could not find a feasible initial point. 在度娘查了查说有可能是参数设定有问题,本人也做了好多次尝试 还是无果!求大神显灵,谢谢!感激不尽。。。
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 F先森 的主题更新
信息提示
请填处理意见