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

k = inp_param(1);

T = inp_param(2);

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

case 8, %неустойчивоезвено 2‑гопор., N_zv = 8, 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 9, %режекторныйфильтр, N_zv = 9, inp_param = [T, ksi1, ksi2]

T = inp_param(1);

ksi1 = inp_param(2);

ksi2 = inp_param(3);

W = (T^2*p^2+2*ksi1*T*p+1)/(T^2*p^2+2*ksi2*T*p+1);

end

%end of function TF_zv

Назначение всех использованных при создании этого файла функций Matlab приводится в приложении.

В программе звенья пронумерованы в соответствии с табл. 2.1, указание на номер моделируемого звена выполняется с помощью параметра N_zv рассматриваемой функции TF_zv. Второй параметр inp_param является векторным и предназначен для передачи значений параметров звена; в зависимости от вида звена он может содержать от одного до трех элементов. Для организации выбора описания звена служит оператор switch…case, при помощи которого производится анализ значения переменной N_zv и переход на вычисление соответствующей передаточной функции.

%Изучение типовых динамических звеньев: колебательное звено

%Очистка всех переменных в памЯти

clear all

%Очистка командного окна

clc

%Закрытие всех предыдущих рисунков

set (0,'ShowHiddenHandles', 'on')

delete (get(0,'Children'))

%Описание колебательного звена (N_zv = 5) через его передаточную функцию

%при различных значениЯх параметров. Параметры колебательного звена

%задаютсЯ вектором inp_param = [k, T, ksi] (см. текст файла TF_zv.m)

%изменЯем k