例えば、下記のドキュメントページを参考に、記述できるかと思います。
上記ページで紹介したコードを使った、サンプルコードを参考までに紹介します。
sys = tf(10,[1 10]);
[mag,phase,w]=bode(sys);
mag = squeeze(mag)
phase = squeeze(phase);
f = w/2/pi;
yyaxis left
semilogx(f,20*log10(mag))
xlabel('周波数 f [Hz]')
ylabel('振幅 [dB]')
yyaxis right
semilogx(f,phase)
ylabel('位相 [deg]')
