用matlab進行《電路理論》課程設計，急求matlab編程

%頻率響應函數：

xlabel('角頻率') %圖形設置

ylabel('增益abs(H)')

num=[0 2 0 ]

den=[1 2 10001]

[h,w]=freqs(num,den);

amp=abs(h);%求幅值

subplot(2,1,1); %分圖繪圖

semilogx(w,amp);%頻率和幅值

title('H(w)的幅頻特性');

xlabel('角頻率（w）');

ylabel('幅度');

ang=angle(h);%相角

subplot(2,1,2);

semilogx(w,ang*180/pi);%頻率和相角

title('H(w)的相頻特性');

xlabel('角頻率');

ylabel('相位（度）');

hold off