Тексты программ:

Получение графиков для передаточной функции W(p)=10/(2p+1):

sys=tf([10],[2 1]);

step(sys); grid;

impulse(sys); grid;

bode(sys); grid;

Получение частотной характеристики моделируемой системы1:

R=1.35; L=0.135; Ce=0.0535; Cm=0.0535; J=0.0000235; k=1; u=27; M=0.0135; dt=0.0001; T=3;

A=[-R/L -Ce/L;Cm/J 0];

B=[k/L 0;0 -1/J];

C=[1 0;0 1];

D=[0 0;0 0];

sys=ss(A,B,C,D);

bode(sys); grid;

Получение частотной характеристики моделируемой системы2:

R=1.35; L=0.135; Ce=0.0535; Cm=0.0535; J=0.0000235; koc=1; k=1; kred=10; u=1; M=0.0135;

dt=0.0001; T=3;

A=[-R/L -Ce/L -(koc*k)/L;Cm/J 0 0;0 1/kred 0];

B=[k/L 0;0 -1/(kred*J);0 0];

C=[1 0 0;0 1 0;0 0 1];

D=[0 0;0 0;0 0];

sys=ss(A,B,C,D);

bode(sys); grid;

Графики:

Передаточная функция W(p)=10/(2p+1):

Реакция на ступенчатое воздействие:

Реакция на импульсное воздействие:

Частотная характеристика:

Моделируемая система 1(с регулированием по разомкнутому контуру):

Частотные характеристики:

Моделируемая система 2(с регулированием по замкнутому контуру):

Ч астотные характеристики