To explain what you are seeing (the envelope you see specifically is not a plotting artifact, but a fundamental DSP concept worth explaining) consider if you instead sampled a 500e6 Hz cosine wave at your sampling rate of 2000e6 Hz, which is exactly 4x the sampling rate. This would result in samples 1, 0, -1, 0, 1.... and since Matlab as you plotted would connect the adjacent samples, it would look like the bottom picture. (And if you used plot(x,'o') you would see just three lines).
Now consider your case where you have offset the frequency by 1 Hz, with this the sample location will "roll" along the sine wave, one cycle in the 1 us duration you are using.
This is exactly what you are seeing. The magnitudes 1 and -1 given in the first case are now moving slowly as the sample locations change slowly along the plot due to the slightly higher frequency, and if you use plot(x, 'o') for that case you will more clearly see the sample locations I am speaking of rather than the lines connecting the adjacent samples that are positive and negative (you will see what appears at 4 shapes varying sinusoidally and in quadrature. The envelope is just the artifact of the relatively low sampling rate (you have slightly less than 4 samples per cycle)- not how Matlab plots it.