bode 関数のオプションなどでラップを有効にしたデータを、出力で得ることはできないようでした。 freqresp 関数で複素周波数応答を出し、angle 関数で位相を出すことで、ラップを有効にした時と同じような結果が得られるかと思います。
sys = tf(1,[2 3 5 1]);
[H,w] = freqresp(sys); % 複素周波数応答
H = squeeze(H);
mag = abs(H);
phase = angle(H)/pi*180;
% ボード線図をプロット
subplot(2,1,1)
semilogx(w,20*log10(mag))
subplot(2,1,2)
semilogx(w,phase)
