当前位置:首页 >> 信息与通信 >>

数字信号处理实验五_图文

实验五:FIR数字滤波器设计及软件实现

3.实验结果及分析:

a.其程序如下: fp=120,fs=150,Ap=0.2,As=60,Fs=1000; wc=(fp+fs)/Fs; B=2*pi*(fs-fp)/Fs; N=ceil(11*pi/B); hn=fir1(N-1,wc,blackman(N)); Hw=abs(fft(hn,1024)); xt=xtg; ywt=fftfilt(hn,xt,N); Y=fft(ywt,1024); subplot(3,1,1),plot(abs(Y)); subplot(3,1,2),plot(abs(Hw)); subplot(3,1,3),plot(ywt); 运行结果如下:

b.等波纹的程序如下: xt=xtg; Fs=1000;T=1/Fs; fp=125;fs=150; Rp=0.1;As=60; wc=(fs+fp)/Fs; B=2*pi*(fs-fp)/Fs; N=ceil(11*pi/B); hn=fir1(N-1,wc,blackman(N)); Hw=abs(fft(hn,1024)); yt=fftfilt(hn,xt,N); figure; subplot(311); plot(Hw); subplot(312); plot(yt);

N1=1000;F1=fft(yt,N1); Tp1=N1*T; k=0:N1-1;f=k/Tp1; subplot(313);stem(f,abs(F1)/max(abs(F1)));grid;title('yt 信 号 的 频 谱 '); axis([0,300,0,1]); xlabel('f/Hz'); ylabel('幅度')

其运行结果如下: