I have no idea what the other values are that you are plotting on the bodeplot axes. Apparently, the only legend entries bodeplot permits are the systems it plots.
It might be easier to get the data from bodeplot and plot everything on ‘ordinary’ axes.
This outlines how to get those data from bodeplot —
sys = tf([1 0 0],[1 1 1 1 1])
h = bodeplot(sys)
Frq = h.Responses.Data.Frequency
Mag = h.Responses.Data.Magnitude
Phs = h.Responses.Data.Phase
.


