%
% �¶‘ÌƒVƒXƒeƒ€�HŠw‚ÌŠî‘b
% �}4.19
%
clear all;
close all;
Gs = tf([0.16], [0.001 0.03 0.3 1], 'ioDelay', 0.18);
f=0.1:0.1:4;
w=2*pi()*f;
[mag, pha] = bode(Gs,w);
mag1 = 20 * log10(squeeze(mag));
pha1 = squeeze(pha);
subplot(2, 1, 1);
semilogx(f, mag1);
xlabel('Frequency (Hz)');
ylabel('Gain (dB)');
subplot(2, 1, 2);
semilogx(f, pha1);
xlabel('Frequency (Hz)');
ylabel('Phase (deg)');