WYDZIAŁ INŻYNIERII METALI I INFORMATYKI PRZEMYSŁOWEJ
Nazwisko i imię: KRAJEWSKI PAWEŁ	AUTOMATYKA I PODSTAWY POMIARÓW
kierunek studiów III rok Metalurgii	Prowadzący dr inż. Dmytro Svyetlichnyy	Data 21.04.2006	Ocena:
specjalność PPM	Temat ćwiczenia: Charakterystyki podstawowych elementów dynamicznych

• Zbadać charakterystyki elementu oscylacyjnego o tramsmitancji:
  , zmieniając po kolei współcynnik wzmocnienia k, pulsacje oscylacji własnych ₀ oraz zredukowany współczynnik tłumienia . Wartości bazowe: k = 1,  ₀ = 1,   = 0,5.

clc

disp('Transmitancje ukladu oscylacyjnego I-go rzedu przy zmianie wartosci parametru k o t')

hold off

for k=[1 2 3 4]

G=tf([k],[1 1 1]

figure(1), step(G), hold on, grid on, legend('k=1','k=2','k=3','k=4'), title('Odpowiedz na wymuszenie skokowe elementu oscylacyjnego przy zmianie wartosci wzmocniemia k')

figure(2), impulse(G), hold on, grid on, legend('k=1','k=2','k=3','k=4'), title('Odpowiedz na wymuszenie impulsowe elementu oscylacyjnego przy zmianie wartosci wzmocniemia k')

figure(3), bode(G), hold on, grid on, legend('k=1','k=2','k=3','k=4'), title('Charakterystyka logarytmiczna elementu oscylacyjnego przy zmianie wartosci wzmocniemia k')

figure(4), nyquist(G), hold on, grid on, legend('k=1','k=2','k=3','k=4'), title('Charakterystyka amplitud-fazowa Nyquista elementu oscylacyjnego przy zmianie wartosci wzmocniemia k')

end

for o=[1 2 3 4]

G=tf([o]^2,[1 [o] [o]^2])

figure(5), step(G), hold on, grid on, legend('o=1','o=2','o=3','o=4'), title('Odpowiedz na wymuszenie skokowe elementu oscylacyjnego przy zmianie wartosci oscylacji wlasnych o')

figure(6), impulse(G), hold on, grid on, legend('o=1','o=2','o=3','o=4'), title('Odpowiedz na wymuszenie

impulsowe elementu oscylacyjnego przy zmianie wartosci oscylacji wlasnych o')

figure(7), bode(G), hold on, grid on, legend('o=1','o=2','o=3','o=4'), title('Charakterystyka logarytmiczna elementu oscylacyjnego przy zmianie wartosci oscylacji wlasnych o')

figure(8), nyquist(G), hold on, grid on, legend('o=1','o=2','o=3','o=4'), title('Charakterystyka amplitud-fazowa Nyquista elementu oscylacyjnego przy zmianie wartosci oscylacji wlasnych o')

end

for t=[0.5 1 1.5 2]

G=tf(1,[1 2*[t] 1])

figure(9), step(G), hold on, grid on, legend('t=0.5','t=1','t=1.5','t=2'), title('Odpowiedz na wymuszenie skokowe elementu oscylacyjnego przy zmianie wartosci tlumienia t')

figure(10), impulse(G), hold on, grid on, legend('t=0.5','t=1','t=1.5','t=2'), title('Odpowiedz na wymuszenie impulsowe elementu oscylacyjnego przy zmianie wartosci tlumienia t')

figure(11), bode(G), hold on, grid on, legend('t=0.5','t=1','t=1.5','t=2'), title('Charakterystyka logarytmiczna elementu oscylacyjnego przy zmianie wartosci tlumienia t')

figure(12), nyquist(G), hold on, grid on, legend('t=0.5','t=1','t=1.5','t=2'), title('Charakterystyka amplitud-fazowa Nyquista elementu oscylacyjnego przy zmianie wartosci tlumienia t')

end

2. Zbadać charakterystyki elementu całkującego o transmitancji
przy różnych wartościach T.

clc

disp('Charakterystyka ukladu calkujacego przy zmianie wartosci parametru T')

hold off

for T=[1 2 3 4]

G=tf(1,[[T] 0])

figure(1), step(G), hold on, grid on, legend('T=1','T=2','T=3','T=4'), title('Odpowiedz na wymuszenie skokowe elementu calkujacego przy zmianie wartosci T')

figure(2), impulse(G), hold on, grid on, legend('T=1','T=2','T=3','T=4'), title('Odpowiedz na wymuszenie impulsowe elementu calkujacego przy zmianie wartosci T')

figure(3), bode(G), hold on, grid on, legend('T=1','T=2','T=3','T=4'), title('Charakterystyka logarytmiczna elementu calkujacego przy zmianie wartosci T')

figure(4), nyquist(G), hold on, grid on, legend('T=1','T=2','T=3','T=4'), title('Charakterystyka amplitudowo-fazowa Nyquista elementu calkujacego przy zmianie wartosci T')

end

3. Zbadać charakterystyki elementu inercyjnego z opóźnienem o transmitancji
   przy różnych wartościach opóźnienia .

clc

disp('Charakterystyka ukladu inercyjnego przy zmianie wartosci opoznienia t')

hold off

for t=[1 2 3 4]

G=tf(1,[2 1],'outputdelay', [t])

figure(1), step(G), hold on, grid on, legend('t=1','t=2','t=3','t=4'), title('Odpowiedz na wymuszenie skokowe elementu inercyjnego przy zmianie wartosci opoznienia t')

figure(2), impulse(G), hold on, grid on, legend('t=1','t=2','t=3','t=4'), title('Odpowiedz na wymuszenie impulsowe elementu inercyjnego przy zmianie wartosci opoznienia t')

figure(3), bode(G), hold on, grid on, legend('t=1','t=2','t=3','t=4'), title('Charakterystyka logarytmiczna elementu inercyjnego przy zmianie wartosci opoznienia t')

figure(4), nyquist(G), hold on, grid on, legend('t=1','t=2','t=3','t=4'), title('Charakterystyka amplitud-fazowa Nyquista elementu inercyjnego przy zmianie wartosci opoznienia t')

end

---