Решение

Исследование соединений

Код программы:

clear all;

close all;

k1=[5];

k2=[5];

TT1=[30];

TT2=[90];

TTT1=[0.2];

TTT2=[2];

KSI1=[0.4];

KSI2=[0.4];

kk1=[10];

kk2=[10];

t1=[TT1 1];

t2=[TT2 1];

i1=tf(k1,t1);

i2=tf(k2,t2);

ipo3=i1*i2;% последовательное соединение для инерционного звена

ipa3=i1+i2;% параллельное соединение для инерционного звена

iobr3=feedback(i1,i2);

% обратное соединение для инерционного звена

c1=[1 0];

ii1=tf(k1,c1);

c2=[1 0];

ii2=tf(k2,c2);

iipo3=ii1*ii2;

% последовательное соединение для интегрируещего звена

iipa3=ii1+ii2;

% параллельное соединение для интегрируещего звена

iiobr3=feedback(ii1,ii2);

% обратное соединение для интегрируещего звена

b1=[TTT1*TTT1 2*KSI1*TTT1 1];

iii1=tf(kk1,b1);

b2=[TTT2*TTT2 2*KSI2*TTT2 1];

iii2=tf(kk2,b2);

iiipo3=iii1*iii2;

% последовательное соединение для колебательного звена

iiipa3=iii1+iii2;

% параллельное соединение для колебательного звена

iiiobr3=feedback(iii1,iii2);

% обратное соединение для колебательного звена

% инерционное звено

figure;

subplot(2,2,1)

bode(ipo3,'r',ipa3,'b',iobr3,'g')

subplot(2,2,2)

nyquist(ipo3,'r',ipa3,'b',iobr3,'g')

subplot(2,2,3)

step(ipo3,'r',ipa3,'b',iobr3,'g')

subplot(2,2,4)

impulse(ipo3,'r',ipa3,'b',iobr3,'g')

% интегрирующее звено

figure;

subplot(2,2,1)

bode(iipo3,'r',iipa3,'b',iiobr3,'g')

subplot(2,2,2)

nyquist(iipo3,'r',iipa3,'b',iiobr3,'g')

subplot(2,2,3)

step(iipo3,'r',iipa3,'b',iiobr3,'g')

subplot(2,2,4)

impulse(iipo3,'r',iipa3,'b',iiobr3,'g')

% колебательное звено

figure;

subplot(2,2,1)

bode(iiipo3,'r',iiipa3,'b',iiiobr3,'g')

subplot(2,2,2)

nyquist(iiipo3,'r',iiipa3,'b',iiiobr3,'g')

subplot(2,2,3)

step(iiipo3,'r',iiipa3,'b',iiiobr3,'g')

subplot(2,2,4)

impulse(iiipo3,'r',iiipa3,'b',iiiobr3,'g')

Графики: