Реферат Колона для перегону коньячного спирту

грамне забезпечення: програма, яка моделює поведінку системи. При тестуванні даної програми було показано регулятор працює адекватно. br/>

Література

1.Стопакевіч А.А.Теорія систем і системний аналіз.Учебнік для вузов.-Київ: ВІПОЛ, 1996.-200с.

2.Демченко В.А.Автоматізація і моделювання технологічних процесів АЕС та ТЕС.-Одеса: Астропріт, 2001.-308с.

3.Стопакевіч А.А.Матлаб. Методичні вказівки до лабораторних робіт, курсового та дипломного проектірованія.-Одеса, 2000.-18с.

Додаток 1

Текст програми

A = [-1/35 0 0 0 0 0 0 0 0 0 0 0 0 0;

0 -1/129 0 0 0 0 0 0 0 0 0 0 0 0;

0 4/48 -2/48 0 0 0 0 0 0 0 0 0 0 0;

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; p> 0 0;

0 -5.1/186; p> 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 =...