Моделирование систем. Елизаров И.А - 94 стр.

UptoLike

ВУЗ: 

ТГТУ | Тамбов

Составители: 

Рубрика: 

92
Kp=20;
Ti=20;
p_z=9.9951e+005; % заданное значение давления пара
M_v_z=2500; % номинальное (заданное) значение массы воды
% характеристики датчика и ИМ
Td=10; %постоянная времени датчика
Tim=20; %постоянная времени ИМ
% цикл решения уравнений модели методом Эйлера
for i=2:n
% поддержание номинальной массы воды в котле
eps=M_v(i-1)-M_v_z; % ошибка регулирования
sum=sum+eps;
G_v=G_v_n-Kp*(eps+sum*dt/Ti); % уравнение регулятора
% ограничение на расход воды
if G_v<=0
G_v=0;
end
if G_v>=5;
G_v=5;
end
% определение давление пара по температуре воды
p(i-1)=interp1(yy,xx,t_v(i-1),'linear');
% определение измеренного значения давления
pi(i)=pi(i-1)+dt*((p(i-1)-pi(i-1))/Td);
eps2=-pi(i)+p_z; % ошибка регулирования
sum2=sum2+eps2;
mu2_=mu2_n+Kp2*(eps2+sum2*dt/Ti2); % уравнение регулятора
% ограничение на проводимость клапана
if mu2_ < 0
mu2_=0;
end
if i==2
mu2_old=mu2_;
end
% определение реальной проводимости клапана
mu2=mu2_old+dt*(mu2_-mu2_old)/Tim;
mu2_old=mu2;
% решение уравнений модели
V_p=V_0-M_v(i-1)/ro_v;
M_p(i)=p(i-1)*V_p*M/(R*(t_v(i-1)+273));
if p_pot > p(i-1)

Страницы