ga(遗传算法)优化计算中，计算结果显示“返回非实数”(return real value.)

bigfeatheryu

我用遗传算法，写一个非线性等式约束，运行时提示错误：
GA running.
Error using ==> gads\private\constrValidate at 25
Constraint function must return real value.


源文件：function [c,ceq]=constraint(A)
c=1;
ceq=-0.6547*(-0.5773)/A(1)+1.327*0.28867/A(2)+(-0.5)*0.5773/A(3)+1.077*(-0.5773)/A(4)+0.5*0.5773/A(5)+0.8273*0.2887/A(6)+(-1.655)*(-0.5773)/A(7)-1;
end


希望高手解答，谢谢

messenger

你这个问题是因为约束函数返回了非实数值。可能是因为表达式计算出现了复数，也可能是你的命令代码哪个地方写错了，因为你的代码没有贴全，现在也搞不清楚。

你可以试试，将计算式用real命令强制取实部，看看是否还出现错误提示，即
ceq=real(-0.6547*(-0.5773)/A(1)+1.327*0.28867/A(2)+(-0.5)*0.5773/A(3)+1.077*(-0.5773)/A(4)+0.5*0.5773/A(5)+0.8273*0.2887/A(6)+(-1.655)*(-0.5773)/A(7)-1);

zhaochunyu347

如果这样强取实部的话会不会造成最后计算不精确呢？或者说取了实部，还能起到预先想要的约束效果吗？
第二个问题，不做任何修改，使用fmincon函数优化时可以得到结果，也没什么提示，二者有什么区别，哪个更好些 2# messenger

messenger

3# zhaochunyu347

取实部不是说一定可以得到正确的结果，而是要发现问题所在。

ga优化要求约束函数为实数，这在help中说得很清楚，而且是以note的形式特别强调。

bigfeatheryu

加了real还是有相同的错误，

我的目标函数文件是：
function V=func_GA(A)
V=(A(1)+A(2)+A(3)+A(4)+A(5)+A(6)+A(7));
end

约束函数文件是：
function [c,ceq]=constraint(A)
c=1;
ceq=(-0.6547*(-0.5773)/A(1)+1.327*0.28867/A(2)+(-0.5)*0.5773/A(3)+1.077*(-0.5773)/A(4)+0.5*0.5773/A(5)+0.8273*0.2887/A(6)+(-1.655)*(-0.5773)/A(7)-1);
end

帮助看看有什么问题，谢谢高手帮助啦

2# messenger

zhaochunyu347

如果在约束条件为复数的情况下用fmincon计算出来对吗？ 4# messenger

messenger

5# bigfeatheryu


最好把ga优化命令也给出，比如初值什么的，这样别人可以重复问题，以利发现问题所在。

加了real还是有相同的错误，

我的目标函数文件是：
function V=func_GA(A)
V=(A(1)+A(2)+A(3)+A(4)+A(5)+A(6)+A(7));
end

约束函数文件是：
function [c,ceq]=constraint(A)
c=1;
ceq=(-0.6547*(-0.5773 ...
bigfeatheryu 发表于 2009-10-17 13:10

messenger

6# zhaochunyu347


没试过，你自己查一下帮助吧。不过，如果你本身表达式应该为实数，计算出来却是复数，应该找出问题所在，否则稀里糊涂地按复数计算也不会正确。另外，如果用real命令也会出来那个错误提示，显然不是出现复数的问题。

如果在约束条件为复数的情况下用fmincon计算出来对吗？ 4# messenger
zhaochunyu347 发表于 2009-10-17 16:46

bigfeatheryu

我用的GUI界面，没用命令格式，貌似GA不用给初值把？
7# messenger

zhaochunyu347

还有个问题就是约束条件可以互相约束吗？就是说第一个约束条件确定了一个变量a的范围，第二个约束条件用到变量a，此时的a是题目本身的范围还是约束条件一确定的范围
谢谢！！ 8# messenger

messenger

9# bigfeatheryu

报歉，命令搞混了。另外，你的约束函数里的，c=1;是什么？

我用的GUI界面，没用命令格式，貌似GA不用给初值把？
7# messenger
bigfeatheryu 发表于 2009-10-17 22:14

messenger

10# zhaochunyu347

应该可以吧。比如优化中经常遇到，a<0，b<a的情况，这样b也<0呀。

还有个问题就是约束条件可以互相约束吗？就是说第一个约束条件确定了一个变量a的范围，第二个约束条件用到变量a，此时的a是题目本身的范围还是约束条件一确定的范围
谢谢！！ 8# messenger
zhaochunyu347 发表于 2009-10-18 08:55

bigfeatheryu

约束函数要求返回c，和ceq，并且将c的表达式大于零，ceq的表达式等于零作为GA的优化的约束条件。
11# messenger

messenger

13# bigfeatheryu

c=1;，你这么写是错的吧，如果没记错的话，这不是要求约束函数1<0吗？

约束函数要求返回c，和ceq，并且将c的表达式大于零，ceq的表达式等于零作为GA的优化的约束条件。
11# messenger
bigfeatheryu 发表于 2009-10-20 21:56

bigfeatheryu

本文题已经解决，自变量下限不能是0

hpctju

同样问题，也解决，自变量下限不能为0