模糊规则优化matlab,遗传算法优化模糊pid控制规则

用遗传算法优化模糊控制规则，有部分代码，不懂计算适应度那块调用子函数，子函数要怎么写？以及M程序如何和simulink相互调用？

clear all;

close all;

b=newfis('yichuan');

b=addvar(b,'input','e',[-3,3]); %Parametere

b=addmf(b,'input',1,'NB','zmf',[-3,-1]);

b=addmf(b,'input',1,'NM','trimf',[-3,-2,0]);

b=addmf(b,'input',1,'NS','trimf',[-3,-1,1]);

b=addmf(b,'input',1,'Z','trimf',[-2,0,2]);

b=addmf(b,'input',1,'PS','trimf',[-1,1,3]);

b=addmf(b,'input',1,'PM','trimf',[0,2,3]);

b=addmf(b,'input',1,'PB','smf',[1,3]);

b=addvar(b,'input','ec',[-3,3]); %Parameterec

b=addmf(b,'input',2,'NB','zmf',[-3,-1]);

b=addmf(b,'input',2,'NM','trimf',[-3,-2,0]);

b=addmf(b,'input',2,'NS','trimf',[-3,-1,1]);

b=addmf(b,'input',2,'Z','trimf',[-2,0,2]);

b=addmf(b,'input',2,'PS','trimf',[-1,1,3]);

b=addmf(b,'input',2,'PM','trimf',[0,2,3]);

b=addmf(b,'input',2,'PB','smf',[1,3]);

b=addvar(b,'output','kp',[-0.3,0.3]);%Parameter kp

b=addmf(b,'output',1,'NB',&#