>> w=logspace(0,3,1000); >> s=j*w; >> H=225./(s.*s+3*s+225); >> AdB=20*log10(abs(H)); >> phase=angle(H)*(180/pi); >> subplot(2,1,1),semilogx(w,AdB),grid >> xlabel('w , rad/s'),ylabel('Amplitude , dB') >> subplot(2,1,2),semilogx(w,phase),grid >> xlabel('w , rad/s'),ylabel('Phase , degre')