%頻率響應函數:
Xlabel('角頻率')%圖形設置
Ylabel('增益絕對值(H)')
num=[0 2 0 ]
den =[1 2 10001]
[h，w]=freqs(num，den)；
amp = ABS(h)；%振幅
支線劇情(2，1，1)；%細分繪圖
semilogx(w，amp)；%頻率和振幅
標題(' H(w ' H(w)'的幅頻特性)；
Xlabel('角頻率(w)')；
Ylabel('振幅')；
ang =角度(h)；%相位角
支線劇情(2，1，2)；
semilogx(w，ang * 180/pi)；%頻率和相位角
標題(' H(w ' H(w)'的相頻特性)；
Xlabel(“角頻率”)；
Ylabel('相位(度)')；
拖延