数字信号处理MATLAB上机实验2014年浅析.ppt

* kwp=fix(wp*M /(2*pi))+1 kws=fix(ws*M /(2*pi))+1 k=[0:1:M-1]; wk = (2*pi/M)*k; Hrs=[ones(1,kwp),zeros(1,M-2*kwp+1),ones(1,kwp-1)] % Hrs = [1,1,1,1,0.5,zeros(1,11),0.5,1,1,1] wd = [0,0.3,0.6,1]; Hrd = [1,1,0,0]; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; angH = [-alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2)]; * H = Hrs.*exp(j*angH); h = real(ifft(H,M)); [Hm,wm] = freqz(h,1,512); mag=abs(Hm); pha=angle(Hm); db=20*log10((mag+eps)/max(mag)) n=[0:255]; fs=2000; x1=sin(2*pi/fs*100*n); x2=sin(2*pi/fs*800*n); x=x1+x2; * y=filter(h,1,x); figure(4); subplot(2,2,1); plot(wk(1:11)/pi,Hrs(1:11),o,