46t.m

包含50个数字信号处理方面的matlab程序

wp=0.25*pi;ws=0.3*pi;Rp=0.25;Rs=50;
N=190;n=[0:1:N-1];
wc=(wp+ws)/2;
hd=ideal_lp(wc,N);
w_han=(blackman(N))';
h=hd.*w_han;
[db,mag,pha,grd,w]=freqz_m(h,1);
subplot(211);stem(n,h,'.');xlabel('n');ylabel('h(n)');
subplot(212);plot(w/pi,db);xlabel('w/pi');ylabel('pha');axis([0,1,-150,10])