There are multiple possible solutions to your question (with different results), but all of them will fail the "without any distortion" requirement.
- you can just take every 1000th element (bad if there is high frequency noise)
new_var = sim_var(1:1000:end);
- do the same, but with a filter first (f.ex. moving average) - better for noise
sim_var2 = movmean(sim_var,1000);
new_var2 = sim_var2(1:1000:end);
attached is a little example to show your the differences between the first two approaches
- Interpolate the values in your postprocessing calculation using "interp1"
- choose a different Stepsize for Simulation (100001 sounds like fixed step) to get to 1001 points