Исследование типовых динамических звеньев, страница 2

hp1=ilaplace(k/(s*(T*s+1)))       %Обратное преобразование Лапласа

hp2=ilaplace(2*k/(s*(T*s+1)))     %изображений переходных функций

hp3=ilaplace(4*k/(s*(T*s+1)))    

hi1=ilaplace(k/(T*s+1))           %Обратное преобразование Лапласа

hi2=ilaplace(2*k/(T*s+1))         %изображений весовых функций

hi3=ilaplace(4*k/(T*s+1))        

ltiview(h1,h2,h3)                 %Просмотр характеристик.

а)                                б)                                            в)                                     г)

д)                                е)                                            ж)                                    з)

и)                                               к)

Рис. 1.3. Результаты моделирования весовых функций типовых звеньев

На рисунках 1.4 – 1.7 представлены результаты исследования апериодического звена в зависимости от коэффициента усиления k.

Программа исследования поведения апериодического звена в зависимости от его постоянной времени представлена ниже.

Рис.1.4. Переходные функции                                             Рис. 1.5. Весовые функции

Рис. 1.6. АЧХ и ФЧХ                          Рис. 1.7. Амплитудно-фазовые характеристики

h1=tf([k],[T,1]);                 %Передаточная функция при Т=2

h2=tf([k],[2*T,1]);               %Передаточная функция при Т=4

h3=tf([k],[4*T,1]);               %Передаточная функция при Т=8

figure(1),step(h1,h2,h3),grid     %Переходныефункции

figure(2),impulse(h1,h2,h3),grid  %Весовыефункции

figure(3),bode(h1,h2,h3),grid     %АЧХиФЧХ

figure(4),nyquist(h1,h2,h3),grid  %Амплитудно-фазовыехарактеристики

symss                            %Ввод символьных переменных

hp1=ilaplace(k/(s*(T*s+1)))       %Обратное преобразование Лапласа

hp2=ilaplace(k/(s*(2*T*s+1)))     %изображений переходных функций

hp3=ilaplace(k/(s*(4*T*s+1)))    

hi1=ilaplace(k/(s+1))             %Обратное преобразование Лапласа

hi2=ilaplace(k/(2*s+1))           %изображений весовых функций

hi3=ilaplace(k/(4*s+1))          

ltiview(h1,h2,h3)                 %Просмотр характеристик.

Результаты исследования представлены на рис. 1.8 – 1.11.

Рис.1.8. Переходные функции                                           Рис. 1.9. Весовые функции

Рис. 1.10. АЧХ и ФЧХ                      Рис. 1.11. Амплитудно-фазовые характеристики

Определение временных и частотных характеристик инерционно-форсирующего звена  в зависимости от коэффициента k:

k=3;                              %Коэффициент усиления

T0=0.2;                           %Постоянная времени Т0

T=3;                              %Постоянная времени Т

h1=tf(k*[T0,1],[T,1]);            %Передаточная функция при k=3

h2=tf(2*k*[T0,1],[T,1]);          %Передаточная функция при k=6

h3=tf(4*k*[T0,1],[T,1]);          %Передаточная функция при k=12

figure(1),step(h1,h2,h3),grid          %Переходныефункции

figure(2),impulse(h1,h2,h3),grid       %Весовыефункции

figure(3),bode(h1,h2,h3),grid          %АЧХиФЧХ

figure(4),nyquist(h1,h2,h3),grid       %АФХ

symss                                 %Ввод символьных переменных

hp1=ilaplace(k*(T0*s+1)/(s*(T*s+1)))   %Обратное преобразование

hp2=ilaplace(2*k*(T0*s+1)/(s*(T*s+1))) %Лапласаизображений

hp3=ilaplace(4*k*(T0*s+1)/(s*(T*s+1))) %переходныхфункций

hi1=ilaplace(k*(T0*s+1)/(T*s+1))       %Обратное преобразование

hi2=ilaplace(2*k*(T0*s+1)/(T*s+1))     %Лапласа изображений

hi3=ilaplace(4*k*(T0*s+1)/(T*s+1))     %весовых функций

ltiview(h1,h2,h3)                      %Просмотр характеристик.

Результаты исследования представлены на рис. 1.12 – 1.15.