Курсовая работа: Автоматическая система регулирования температуры

U=zeros(1,numel(time));

Up=zeros(1,numel(time));

alpha=zeros(1,numel(time));

Un=zeros(1,numel(time));

Pn=zeros(1,numel(time));

err1=zeros(1,numel(time));

T=zeros(1,numel(time));

Ud=zeros(1,numel(time));

%Начальные условия

%--------------------------

err(1)=U0;

d_err(1)=0;

U(1)=Kp*err(1);

if (U(1)>=10)

Up(1)=10;

end;

alpha(1)=pi/(1+K0*Up(1));

Un(1)=(((pi-alpha(1)+(sin(2*alpha(1))/2))/pi)^0.5)*Uc;

Pn(1)=(Un(1))^2/Rn;

err1(1)=Pn(1)-((T(1)-Tc)*gamma);

T(1)=0; %начальная температура объекта

Ud(1)=0;

%-------------------------

for k=2:1:numel(time)

err(k)=U0-Ud(k-1);

d_err(k)=(err(k)-err(k-1))/dt;

U(k)=U(k-1)+((dt/Tp)*(Kp*Tp*d_err(k-1)+Kp*err(k-1)));

if (U(1)>10)

Up(1)=10;