How to generate eye diagram in Matlab ( without communication toolbox)

Question

I want to plot an eye dagram of output of modulator.

In Maltab I tried plot( I, Q), but it doesnt give me the diagram.

In Matlab there is a function eyediagram(x,n,period,offset,plotstring,h) in Communicatio Toolbox. Unfortunately I dont have the licence and can't use it.

Does anyone know another possobility how to generate the eye diagram?

EDIT 1 Eye diagram with my data

EDIT 2

I have added a dely and my new result is

a = 10;
eyediagram(bb(a;end), ovs)

EDIT 3

sps = 4;
data = randi([0 1],1000,1);
gmskMod = comm.GMSKModulator('BitInput',true,'PulseLength',5, 'BandwidthTimeProduct' , 0.25, ...
    'SamplesPerSymbol',sps);
modSigGMSK = gmskMod(data);
eyediagram(modSigGMSK,sps,1,sps/2)

Answer 1

Plotting I vs Q (which is real vs imaginary) is a polar plot, not the eye diagram. The eye diagram is plotted with I alone or Q alone over one or more symbol periods (typically just one or two) and repeatably synchronized with the symbol rate so that the result is multiple plots superimposed showing how the waveform behaves from symbol to symbol.

This can be done easily in MATLAB using the reshape function on either I or Q which will create a matrix that when plotted will be the eye diagram. However the waveform should be resampled to be an integer number of samples and many samples per symbol so that the trajectories between the actual samples are visible.

This link contains example code I have written in Python that can easily be ported to Python, as well as other open source code available for creating eye diagrams.

Here is an example of that done for a 16QAM waveform where the blue dots are the original waveform samples at 4 samples per symbol (with a time offset error that the eye diagram shows us), and the trajectories shown are after resampling to 64 samples per symbol and plotting the eye diagram to show two symbol durations:

Note for comparison that if we plotted the first 1200 samples of just I alone without the retrace over the two symbol periods as done above, the plot would look like this (the dashed lines show the two symbol period boundaries used for plotting the eye diagram above:

Please see this other post with further details on what the eye diagram is an how to use it.

This post shows how the eye diagram will appear if there is a phase offset rather than time (delay) offset as in the example above. (Note phase and delay are not the same thing! Phase is a rotation, not a delay but a delay will cause a rotation at any given frequency).

Below shows the expected eye diagram over two symbols for GMSK with BT=1.

Note that if we rotate the constellation 45° then the eye diagram would look like the following (this is just changing the reference for the constellation but doesn't change anything with the actual waveform).