stem(n,angle(y)); title('phase');
subplot(4,1,4); stem(n,x1);
title('resumed signal');
运行结果：
15
实验五：
第一题（书上有参考代码，不再赘述） 运行结果：
16
第二题
源代码：
N=64; k=0:1:N-1; w=2*pi*k/N; z=exp(j*w); Xk=(z)./(z-0.8); xn=real(ifft(Xk,N)); xtilde=xn'*ones(1,2); xtilde=(xtilde(:))'; subplot(2,1,1); stem(0:127,xtilde); title('64μ?μ?IFFT·ù?è?×');xlabel('n');
ylabel('xtilde(n)'); N=16; k=0:1:N-1; w=2*pi*k/N; z=exp(j*w); Xk=(z)./(z-0.8); xn=real(ifft(Xk,N)); xtilde=xn'*ones(1,8); xtilde=(xtilde(:))'; subplot(2,1,2); stem(0:127 ,xtilde); title('16μ?μ?IFFT·ù?è?×');xlabel('n');
ylabel('xtilde(n)');
17
运行结果：
18
实验七：
第一题： 源代码： n=50;
Wn=0.45;
b1=fir1(n,Wn,boxcar(n+1)); [H,w]=freqz(b1,1,512,2); subplot(4,1,1); plot(w,abs(H));grid;
title('矩形窗低通FIR数字滤波器'); b2=fir1(n,Wn,hanning(n+1)); [H,w]=freqz(b2,1,512,2); subplot(4,1,2);
plot(w,abs(H));grid;
title('汉宁窗低通FIR数字滤波器'); b3=fir1(n,Wn,hamming(n+1)); [H,w]=freqz(b3,1,512,2); subplot(4,1,3);
plot(w,abs(H));grid;
title('海明窗低通FIR数字滤波器'); b4=fir1(n,Wn,blackman(n+1)); [H,w]=freqz(b4,1,512,2); subplot(4,1,4);
plot(w,abs(H));grid;
title('希拉克曼窗低通FIR数字滤波器');
运行结果：
矩形窗低通FIR数字滤波器21021021021000.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.7希拉克曼窗低通FIR数字滤波器0.80.9100.10.20.30.40.50.60.7海明窗低通FIR数字滤波器0.80.9100.10.20.30.40.50.60.7汉宁窗低通FIR数字滤波器0.80.91
19