%Samp9_9
clf;h=fir1(30,0.2,boxcar(31));  %設計30階截止頻率為0.2的FIR濾波器
r=randn(16384,1);       %產生白噪聲隨機向量
x=filter(h,1,r);        %對隨機向量濾波
psd(x,1024,10000,kaiser(512,5),0,0.95);   %采用Kaiser窗進行功率譜估計
xlabel('頻率/Hz');ylabel('功率譜/dB')