about gatool

Question

Natalia 2012-1-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/25395-about-gatool

Hello! I have a fitness functiion:

 function z = my_fun(x)  
 Skorost8.input(1,1).mf(1,1).params=[x(1),x(2),x(3)];
 Skorost8.input(1,1).mf(1,2).params=[x(4),x(5),x(6)];
 Skorost8.input(1,1).mf(1,3).params=[x(7),x(8),x(9)];
 Skorost8.input(1,1).mf(1,4).params=[x(10),x(11),x(12)];
 Skorost8.input(1,2).mf(1,1).params=[x(13),x(14),x(15)];
 Skorost8.input(1,2).mf(1,2).params=[x(16),x(17),x(18)];
 Skorost8.input(1,2).mf(1,3).params=[x(19),x(20),x(21)];
 load_system('petina.mdl');
 sim('petina.mdl');
 close_system('petina.mdl');
 z=simout.signals.values;

And petina.mdl file, which it use. The mdl-file write to workspace a massive of values, which must be a minimal (z=simout.signals.values;)

When I try to use it with gatool I get:

 >> x=ga(@my_fun,21);
  ??? Subscripted assignment dimension mismatch.
Error in ==> fcnvectorizer at 14
            y(i,:) = feval(fun,(pop(i,:)));
Error in ==> makeState at 47
        Score =
        fcnvectorizer(state.Population(initScoreProvided+1:end,:),FitnessFcn,1,options.SerialUserFcn);
Error in ==> gaunc at 41
state = makeState(GenomeLength,FitnessFcn,Iterate,output.problemtype,options);
Error in ==> ga at 279
        [x,fval,exitFlag,output,population,scores] = gaunc(FitnessFcn,nvars, ...
Caused by:
    Failure in user-supplied fitness function evaluation. GA cannot continue.

How to use correctly gatool in this case? thank you very much!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sean de Wolski 2012-1-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/25395-about-gatool#answer_33312

在 MATLAB Online 中打开

dbstop if error

Which will stop deep inside of ga(). Now put a break point in first line of your objective function. Run the feval line that causes the error and step through your fitness function to find the error and analyze the inputs.

5 个评论
显示 3更早的评论隐藏 3更早的评论

Natalia 2012-1-5

Thank you very much for your answer!

I have done dbstop if error and have put a break point.

And what I have...

>> dbstop if error

>> x=ga(@my_fun,21);

6 Skorost8.input(1,1).mf(1,1).params=[x(1),x(2),x(3)];

7 Skorost8.input(1,1).mf(1,2).params=[x(4),x(5),x(6)];

8 Skorost8.input(1,1).mf(1,3).params=[x(7),x(8),x(9)];

9 Skorost8.input(1,1).mf(1,4).params=[x(10),x(11),x(12)];

11 Skorost8.input(1,2).mf(1,1).params=[x(13),x(14),x(15)];

12 Skorost8.input(1,2).mf(1,2).params=[x(16),x(17),x(18)];

13 Skorost8.input(1,2).mf(1,3).params=[x(19),x(20),x(21)];

16 load_system('petina.mdl');

17 sim('petina.mdl');

18 close_system('petina.mdl');

19 z=simout.signals.values;

End of function my_fun.

End of function globaloptim\private\validate>@(x)fitness(x,FitnessFcnArgs{:}).

??? Subscripted assignment dimension mismatch.

Error in ==> fcnvectorizer at 14

y(i,:) = feval(fun,(pop(i,:)));

Error in ==> makeState at 47

Score =

fcnvectorizer(state.Population(initScoreProvided+1:end,:),FitnessFcn,1,options.SerialUserFcn);

Error in ==> gaunc at 41

state = makeState(GenomeLength,FitnessFcn,Iterate,output.problemtype,options);

Error in ==> ga at 279

[x,fval,exitFlag,output,population,scores] = gaunc(FitnessFcn,nvars, ...

Caused by:

Failure in user-supplied fitness function evaluation. GA cannot continue.

>>

On the first step it generates input values and passes my function up to the end, substituting values into the fis-struct in workspace. Open and simulate mdl. Not clearly why there is an error. May be the resalt of function is a massiv of values?

Sean de Wolski 2012-1-5

Yes, a 1x450 would be 450 objectives (which sounds intuitively like a lot (way too many) to me!). Especially since you only want to optimize 21 variables.

If you look in the doc for ga() on line one it explains the output:

"x = ga(fitnessfcn,nvars) finds a local unconstrained minimum, x, to the objective function, fitnessfcn. nvars is the dimension (number of design variables) of fitnessfcn. The objective function, fitnessfcn, accepts a vector x of size 1-by-nvars, and returns a scalar evaluated at x."

Could you explain what you're trying to optimize?

Natalia 2012-1-5

I try to optimize the parametrs of fuzzy-logic controller, which describe into the fis file and use in my mdl. gamultiobj has worked. I don't have a erorrs :) Thank you :)

请先登录，再进行评论。