Лабораторная работа: Исследование линейных систем

%ФункциЯ описаниЯ звена в виде передаточной функции

%

function W = TF_zv (N_zv, inp_param)

p = tf('p');

switch N_zv

case 1, %интегрирующеезвено, N_zv = 1, inp_param = [k]

k = inp_param(1);

W = k/p;

case 2, %инерционноезвено, N_zv = 2, inp_param = [k, T]

k = inp_param(1);

T = inp_param(2);

W = k/(T*p+1);

case 3, %реальноедифференцирующеезвено, N_zv = 3, inp_param = [k, T]

k = inp_param(1);

T = inp_param(2);

W = (k*p)/(T*p+1);

case 4, %интегро-дифференцирующеезвено, N_zv = 4, inp_param = [T1, T2]

T1 = inp_param(1);

T2 = inp_param(2);

W = (T1*p+1)/(T2*p+1);

case 5, %колебательноезвено, N_zv = 5, inp_param = [k, T, ksi]

k = inp_param(1);

T = inp_param(2);

ksi = inp_param(3);

W = k/(T^2*p^2+2*ksi*T*p+1);

case 6, %неминимально-фазовое звено 1‑го пор., N_zv = 6, inp_param = [T1, T2]

T1 = inp_param(1);

T2 = inp_param(2);

W = (1‑T1*p)/(1+T2*p);