hnwplot.m

来自「MATLAB 基础编程 包含一些很基础的程序」· M 代码 · 共 12 行

function a=hnwplot(hn,k)
N=length(hn);
n=0:N-1;
hw=fft(hn,512);w=0:511;w=2*w/512;
subplot(1,2,k+1)
plot(w,20*log10(abs(hw)));grid;title('幅度特性(dB)');
axis([0,max(w),-100,5]);
xlabel('ω/π');ylabel('Magnitude(dB)')
subplot(1,2,k);stem(n,hn,'.');
axis([0,N-1,min(hn)-0.05,max(hn)+0.1]);xlabel('n');ylabel('h(n)')
line([0,N-1],[0,0])