Курсовая работа: Колонна для перегона коньячного спирта
0
0 0 -1/38 0 0 0 0 0 0 0 0 0 0;
0
0 0 4/9 -2/9 0 0 0 0 0 0 0 0 0;
0
0 0 0 0 -1/110 0 0 0 0 0 0 0 0;
0
0 0 0 0 4/134 -2/134 0 0 0 0 0 0 0;
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);
figure(1)
subplot(2,2,1);plot(y(:,1));grid;ylabel('u1,МПа');title('Razgon
u1,1%');
subplot(2,2,3);plot(y(:,2));grid;ylabel('u2,');
%u2
y=dstep(Ad,Bd,C,D,2,G);
subplot(2,2,2);plot(y(:,1));grid;ylabel('y1,М3/с');title('Razgon
u2,1%');
subplot(2,2,4);plot(y(:,2));grid;ylabel('y2,кг/с');
A1=[Ad
zeros(14,2);C eye(2)];
B1=[Bd;
zeros(2)];
C1=[zeros(2,14)
eye(2)];
Q2=1e6*[1
0;0 1];
V=C'*Q2*C;
Q=[V
zeros(14,2);zeros(2,14) eye(2)];
R=eye(2);
Q1=eye(16);
R1=eye(2);
K=dlqr(A1,B1,Q,R);
L=dlqr(A1',C1',Q1,R1)';
K1=K(:,1:14);
K2=K(:,15:16);
L1=L(1:14,:);
L2=L(15:16,:);
Ar=[Ad-Bd*K1
-Bd*K2-L1 L1; C eye(2)-L2 L2; zeros(2,14) zeros(2) eye(2)];
Br=[zeros(14,2);
zeros(2); eye(2)];
Cr=[-K
zeros(2)];
Az=[Ad
Bd*Cr; Br*C Ar];
Bf=[Bd;
zeros(18,2)];
Bz=[zeros(14,2);
Br];
Cz=[C
zeros(2,18)];
x=zeros(14,1);xr=zeros(18,1);
u=zeros(2,1);
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');
|