%
% �¶‘ÌƒVƒXƒeƒ€�HŠw‚ÌŠî‘b
% �}4.67
%
close all
clear all
r = 0.039;
R = 0.9;
C= 1.333;
num = [1/r 1/(r*R*C)];
den = [1 1/C*(1/r + 1/R)];
Gs = tf(num, den);
[mag pha w] = bode(Gs);
mag1 = 20 * log10(squeeze(mag));
pha1 = squeeze(pha);
subplot(2, 1, 1);
semilogx(w, mag1);
xlabel('Angular frequency (rad/s)');
ylabel('Gain (dB)');
subplot(2, 1, 2);
semilogx(w, pha1);
xlabel('Angular frequency (rad/s)');
ylabel('Phase (deg)');
