Mark specific frequencies in smithchart
There is an easy workaround to plot corresponding frequencies from S-parameter vector:   

hold on
plot(S22(1), 'r*')
plot(S22(159),'r*')
plot(S22(end),'r*')  

The result looks like this: