cstr.m

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

clear 
x1(1,1)=100;
x2(1,1)=3;
% u1(1,1)=1;
% u2(1,1)=1.05;
 u1(1:1200,1)=[0.2*ones(600,1);0.6*ones(600,1)];
 u2(1:1200,1)=[1*ones(200,1);1.5*ones(200,1);0.5*ones(200,1);2*ones(200,1);zeros(200,1);0.7*ones(200,1)];%,0.95*ones(50,1),1.15*ones(50,1),0.85*ones(50,1)];
k1=0.2;
k2=1;
CB1=24.9;
CB2=0.1;
for k=2:1200
 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);
%u1(k,1)=2*rand;
%u2(k,1)=2*rand;
end