% Soll- und Ist-Verlauf
soll1 = [0 0 0 0 0 0 0 50 50 50 50 50 50 50];
t1 = 1:1:14;
ist1 = [0 0 0 0 0 0 0 0 0 50 50 50 50 50];
t = 1:0.01:14;
soll = interp1(t1, soll1, t, 'spline');
ist = interp1(t1, ist1, t, 'spline');

% Störungsüberlagerung
for i=1:length(ist)
   ist(i) = ist(i) + 5*rand; 
end

figure(1)
plot(t, soll,'r', t, ist, 'g', t, soll-ist,'b')
legend ('Soll-Verlauf','Ist-Verlauf','Differenz')