lin_seed.m

用MATLAB实现遗传算法中的好多问题

% 用线性模型设置种子群体
U_x='[ones( num_data,1) u1 ';
for i=2:num_inp,
   U_x=[U_x 'u' num2str(i) ' '];
end
U_x=['U_data=' U_x '];'];
for k=num_inp:-1:1,
   U_x = strrep(U_x,['u' num2str(k)],['u(1:num_data,' num2str(k) ')'] );
end
eval(U_x);   
Y_data=[y(1:num_data)];
theta1=inv(U_data'*U_data)*U_data'*Y_data;
% 将线性模型符号化成串
pop_x=['(' '([' num2str(theta1(1))  '])' ];
for i=1:num_inp,
 pop_x=[pop_x '+'];
   if i < num_inp
      pop_x=[pop_x '('];
   end
   pop_x=[ pop_x '(([' num2str(theta1(i+1)) '])*(u' num2str(i) '))' ];
end

for j=1:num_inp,
   pop_x=[pop_x ')'];
end
seed_pop=pop_x;