%MATLAB-skript: mfd_lsim.m
tstopp=50;
t=[0:.1:tstopp]'; %Generer array for t
U=4;%Spranghoyden i u.
%Lager sprangtidsserie:
u=U*ones(length(t),1);%ones-funksjonen lager matrise av 1-ere
m=20; Kf=2; D=4;
A=[0,1;-Kf/m,-D/m];
B=[0;1/m];
C=[1 0];%Sørger for at utgangen er y=x1 (posisjonen).
D=[0];
%Definerer LTI-modell med navn mfd vha. ss-funksjonen:
mfd=ss(A,B,C,D);
x0=[0;0];
y=lsim(mfd,u,t);%Simulerer med lsim. Respons i y. 
%Plotting:
figure(1)
plot(t,y),ylabel('[m]'),title('Posisjon, y')
grid