untitled2.m

CSTR产生数据的程序。由于CSTR具有广泛的应用性

options = odeset('RelTol',1e-6,'AbsTol',[1e-2 1e-2 1e-2 1e-2]);
%U(1:1200,:)=[[1.6*ones(600,1);0.6*ones(600,1)],[2*ones(200,1);0.1*ones(200,1);1.1*ones(200,1);0.1*ones(200,1);1.1*ones(200,1);2*ones(200,1);]];
%%%%%%
X(1,:)=[100,6.273];
U(1,:)=[1,1.05];
ps=0.05;
for i=2:1200
   [t,x] = ode45(@cstr_d,[0 2000],[X(i-1,:),U(i-1,:)],options );
    X(i,1)=x(2,1);
    X(i,2)=x(2,2);
   if rand<ps
       U(i,1)=1.9*rand+0.1;     
   else
      U(i,1)=U(i-1,1) ;
   end
   if rand<ps     
       U(i,2)=1.9*rand+0.1;
   else   
      U(i,2)=U(i-1,2) ;
   end
end
%%%%%%%
%[t,x] = ode45(@cstr_d,[0 2000],[X(1,:),0,0],options );
% X(1,1:2)=[100,2.72];
% % u1=1;
% % u2=1.05;
%  for i=2:1200
%     [t,x] = ode45(@cstr_d,[0 2000],[X(i-1,:),U(i-1,:)],options );
%      X(i,1)=x(2,1);
%      X(i,2)=x(2,2);
%  end
%[t,z] = ode45(@cstr_d,[0 2000],[X(1,:),u1,u2],options );
% for i=1:200
%     y(i)=-0.2*sqrt(X(i,1))+X(i,1);
%     X(i+1,1)=y(i);
% end


%  x1(k,1)=u1(k-1)+u2(k-1)-k1*sqrt(abs(x1(k-1)))+x1(k-1);
%  x2(k,1)=(CB1-x2(k-1))*u1(k-1)/x1(k-1)+(CB2-x2(k-1))*u2(k-1)/x1(k-1)-k2*x2(k-1)/(1+x2(k-1))^2+x2(k-1);