Use the 'markerindices' argument, introduced in 2016b
plot(x, y, 'markerindices', [1:10:lenght(x)])
puts markers on every 10th point. In older released you need to use something like
plot(x, y, 'k-',...
Which is really annoying as you get 2 different line objects.