%TEL230 Reguleringsteknikk
%13.10.03. Finn Haugen
%Løsningsforslag til Matlab-øvingen i leksjon 4:

clear

q=0.2;
V=1;

teller=[1];
nevner=[V,q];
H_blandetank=tf(teller,nevner);

tslutt=60;
t=[0:0.1:tslutt]';
wA=0.8*ones(length(t),1);
cA=lsim(H_blandetank,wA,t);

figure(1)
subplot(2,1,1)
plot(t,wA)
grid
ylabel('wA [mol]')
xlabel('t [sek]')
subplot(2,1,2)
plot(t,cA)
grid
ylabel('cA [mol]')
xlabel('t [sek]')


