44t.m

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

%此程序是正确的，他是直接从数字出发，同我们实验二里的第四个小问相同的思路。
wp1=0.3*pi;wp2=0.6*pi;ws1=0.4*pi;ws2=0.7*pi;Rp=0.5;Rs=40;
N=81;n=[0:1:N-1];
wc1=(wp1+ws1)/2;wc2=(wp2+ws2)/2;
hd=ideal_lp(wc1,N)+ideal_lp(pi-wc2,N).*exp(j*pi*n);
w_han=(hanning(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,-100,10])