资源描述：
《FIR数字滤波器设计与软件实现.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、湖南师范大学职业技术学院（工学院）实验数据报告单实验课程：数字信号处理实验题目：FIR数字滤波器设计与软件实现实验日期：2019年5月20日专业：电信年级：2016班：一班姓名：陈哲瀚学号：5实验目的：《数字信号处理》P293实验内容：《数字信号处理》P293-P295程序清单或流程图：（1）Xtg:functionxt=xtgN=1000;Fs=1000;T=1/Fs;Tp=N*T;t=0:T:(N-1)*T;fc=Fs/10;f0=fc/10;mt=cos(2*pi*f0*t);ct=cos(2*pi*fc*t);xt=mt.*ct;nt=2*ran
2、d(1,N)-1;fp=150;fs=200;Rp=0.1;As=70;fb=[fp,fs];m=[0,1];dev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1)];[n,fo,mo,W]=remezord(fb,m,dev,Fs);hn=remez(n,fo,mo,W);yt=filter(hn,1,10*nt);xt=xt+yt;fst=fft(xt,N);k=0:N-1;f=k/Tp;subplot(2,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(t)');axi
3、s([0,Tp/5,min(xt),max(xt)]);title('(a)信号加噪声波形')subplot(2,1,2);plot(f,abs(fst)/max(abs(fst)));grid;title('(b)信号加噪声的频谱')axis([0,Fs/2,0,1.2]);xlabel('f/Hz');ylabel('幅度')（2）clc;clearall;closeall;NN=1000;fs1=1000;T=1/fs1;xt=xtg;t=0:T:(NN-1)*T;Tp=NN*T;k=0:NN-1;f=k/Tp;%rs=60dB，因此采用布莱克曼窗
4、%转换为数字频率fp=120;fs=150;wp=2*pi*fp/fs1;ws=2*pi*fs/fs1;bt=ws-wp;%计算滤波器的阶数N0=ceil(11*pi/bt);N1=N0+mod(N0+1,2);wc=(wp+ws)/2/pi;hn=fir1(N1-1,wc,blackman(N1));figure;[H,W]=freqz(hn,1,512,2);%滤波器的幅频曲线plot(W,20*log10(abs(H)));gridtitle('FIR滤波器幅频曲线')xlabel('omega/pi');ylabel('幅度(dB)')y1=
5、filter(hn,1,xt);y1_fft=fft(y1,NN);figure;subplot(2,1,1);plot(t,y1);gridtitle('信号')xlabel('t/s');ylabel('x(t)');axis([0,0.5*Tp,min(y1),max(y1)]);%axis([0.2*Tp,0.7*Tp,min(y1),max(y1)]);subplot(2,1,2);stem(f,abs(y1_fft)/max(abs(y1_fft)),'.');title('信号频谱')xlabel('f/Hz');ylabel('幅度')a
6、xis([0,fs1/5,0,1.2]);grid运行结果：实验总结：用窗函数法设计的滤波器，如果在阻带截止频率附近刚好满足，则离开阻带截止频率越远，阻带衰减富余量越大，即存在资源浪费。几种常用的典型窗函数的通带最大衰减和阻带最小衰减固定，且差别较大，又不能分别控制。所以设计的滤波器的同代最大衰减和阻带最小衰减通常都存在较大富余。Blackman窗函数阻带最小衰减位74Db，而指标位60dB。指导老师签名:
