opt_plot.m

王小平《遗传算法——理论、应用与软件实现》随书光盘

% opt_plot : For generating a plot of optimum T1,T2 for 
%            varying values of CAF1 (for 2  CSTR)

clear all

T10=341;
T20=341;

options=foptions;
%options(1)=1;
%options(14)=

dQF1=0;

%for k=1:15,
%   dCAF1(k)=-100 + k*10;
%  x0=[T10 T20]';
%   
%   x=fmins('cstr_ss8',x0,options,[],dCAF1(k),dQF1);
%   T1(k)=x(1);
%   T2(k)=x(2);
%   
%   xp=fmins('gp_mod',x0,options,[],dCAF1(k));
%   T1p(k)=xp(1);
%   T2p(k)=xp(2);
%
%end

dCAF1=-100;

for k1=1:20,
   T1(k1)=310+k1*2.5;
   for k2=1:20,
      T2(k2)=320+k2*1.75;
      
      f=gp_mod([T1(k1) T2(k2)],dCAF1);
      CB2(k1,k2)=-f;
   end  
end

dCAF1=150;

for k1=1:20,
   T1(k1)=310+k1*2.5;
   for k2=1:20,
      T2(k2)=320+k2*1.75;
      
      f=gp_mod([T1(k1) T2(k2)],dCAF1);
      CB2p(k1,k2)=-f;
   end  
end

surf(T2,T1,CB2-CB2p)