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

0 0 0 0 0 0 0 -1/13.5 0 0 0 0 0 0;

0 0 0 0 0 0 0 0 -1/98 0 0 0 0 0;

0 0 0 0 0 0 0 0 4/133 -2/133 0 0 0 0;

0 0 0 0 0 0 0 0 0 0 -1/50 0 0 0;

0 0 0 0 0 0 0 0 0 0 4/12 -2/12 0 0;

0 0 0 0 0 0 0 0 0 0 0 0 -1/186 0;

0 0 0 0 0 0 0 0 0 0 0 0 4/150 -2/150];

B=[-42.5/35 0;

-72.5/129 0;

0 0;

1720/38 0;

0 0;

730/110 0;

0 0;

0 0.994/13.5;

0 0.459/98;

0 0;

0 -6.9/50;

0 0;

0 -5.1/186;

0 0];

C=[1 -1 1 0 0 0 0 1 -1 1 0 0 0 0;

0 0 0 -1 1 -1 1 0 0 0 -1 1 -1 1];

D=[0 0;0 0];

dt=0.1/max(abs(eig(A)));

t=0:dt:999;

G=length(t);

[Ad Bd]=c2d(A,B,dt);

y=dstep(Ad,Bd,C,D,1,G);