% skrypt symulacyjny regulacji dwustawnej
% Karol Kozlowski, 2008
%

model = 'dwustawna';

t0 = 5;    %czas za³¹czenia
SP = 0.3;   %nastawa
T  = 60;    %czas symulacji

sim(model);

n  = size(simout,1);

t  = simout(:,1);
y0 = simout(:,2);
u  = simout(:,3);
yt = simout(:,4);
yr = (max(yt) + min(yt(n/2:n))) * ones(1, size(t,1))/2;

h1 = subplot(2,1,1);
p1 = plot(t, y0, t, yt, t, yr, 'r');
xlabel('czas [s]');
l1 = legend('y_0', 'y(t)', 'y_r');
grid on;

h2 = subplot(2,1,2);
p2 = plot(t, u, 'r');
xlabel('czas [s]');
legend('u');
grid on;

axis(h1,[min(t) max(t) min(yt) max(yt)*1.2])
axis(h2,[min(t) max(t) min(u) max(u)*1.5])
set(p1, 'LineWidth',2);
set(p2, 'LineWidth',2);
set(l1, 'Orientation', 'horizontal');
