Курсовая работа: Колонна для перегона коньячного спирта

yy=[]; uu=[];f=[.0010;.0010];z=[0;0];

for i=1:2000,

y=C*x; e=-z+y;

u=Cr*xr; xr=Ar*xr+Br*e;

y=C*x; x=Ad*x+Bd*(u+f);

yy=[yy; y']; uu=[uu; u'];

end

x1=x; xr1=xr; u1=u;

figure(2)

subplot(2,2,1); plot(yy(:,1));grid;ylabel('y1,MPa');title('Perehod proces Braga ');

subplot(2,2,3); plot(yy(:,2));grid;ylabel('y2,');

subplot(2,2,2); plot(uu(:,1));grid;ylabel('u1,M3/c');title('Perehod proces Par ');

subplot(2,2,4); plot(uu(:,2));grid;ylabel('u2,kg/c');