clear;num=[0,0,1];den=[1,1,10];
w=logspace(-1,2,100);
Gjw=polyval(num,j*w)./polyval(den,j*w);
Gain=20*log10(abs(Gjw));
Phase=rad2deg(angle(Gjw));
subplot(2,1,1);
semilogx(w,Gain);grid('on');
xlabel('Frequency[rad/s]');ylabel('Gain[dB]');
subplot(2,1,2);
semilogx(w,Phase);grid('on');
xlabel('Frequency[rad/s]');ylabel('Phase[deg]');
