Расчет параметров настройки ПИД-регулятора

Передаточная функция ПИД-регулятора определяется по формуле

,

где .

Т. к. для ПИД-регулятора необходимо определить три коэффициента, то построим кривую равной степени затухания с учетом времени дифференцирования _.

Script 7:

>> w=0.15:0.001:0.26;

>> Wex=(0.9*((j-m).*w).^2+7*(j-m).*w+2.2)./...

(336*((j-m).*w).^3+146*((j-m).*w).^2+21*(j-m).*w+1);

>> Win=1./Wex;

>> R=real(Win);

>> I=imag(Win);

>> Ki=w*(m^2+1).*(I+w*2.04);

>> Kp=m.*I-R+2*m.*w*2.04;

>> plot(Kp,Ki);xlabel('Axis Kp');ylabel('Axis Ki');grid

Рисунок 6 – Кривая равной степени затухания

Коэффициенты k_p и k_i определяем по кривой равной степени затухания (рисунок 6). k_p = 2.05 и k_i = 0.18. Тогда передаточная функция ПИД-регулятора будет иметь вид

.

Построим переходную характеристику САР с ПИД-регулятором.

Script 8:

>> Wop=tf([0.9 7 2.2],[336 146 21 1]);

>> Wap3=tf([2.04 2.05 0.18],[1 0]);

>> W3=series(Wap3,Wop)

Transfer function:

1.836 s^4 + 16.13 s^3 + 19 s^2 + 5.77 s + 0.396

-----------------------------------------------

336 s^4 + 146 s^3 + 21 s^2 + s

>> Fi3=feedback(W3,1)

Transfer function:

1.836 s^4 + 16.13 s^3 + 19 s^2 + 5.77 s + 0.396

-----------------------------------------------

337.8 s^4 + 162.1 s^3 + 40 s^2 + 6.77 s + 0.396

>> step(Fi3);grid

Рисунок 7 – Переходная характеристика САР с ПИД-регулятором

По формуле (3) определяем степень затухания ψ

Script 9:

>> ((1.33-1)-(1.05-1))/(1.33-1)

ans =0.8485

Полученная степень затухания примерно совпадает с заданной, значит параметры настройки регулятора выбраны верно.