FIR滤波器设计与软件实现

Fir滤波器设计和软件实现 加实验报告，比较简单的一种，低通

clear all;close all;
xt=xtg;N=1000;
fp=120;fs=150;Rp=0.1;As=60;Fs=1000;%指标
%======================窗函数========================
B=2*pi*（fs-fp）/Fs;
M=ceil（11*pi/B）;%blackman窗过渡带宽度Bt精确值为11pi/N
wc=（fp+fs）/Fs;                     %书7.2.20
hn=fir1（Mwcblackman（M+1））;
HW=fft（hn1024）;
Hw=abs（HW）;
ywt=fftfilt（hnxtN）;
%画图
figure（2）
subplot（211）
n=0:M;stem（nhn‘.‘）;grid on;
axis（[0M+11.2*min（hn）1.2*max（hn）]）;
xlabel（‘n‘）;ylabel（‘h（n）‘）;title（‘（a）单位脉冲响应‘）
subplot（212）;
f=[0:1023]*Fs/1024;
plot（f20*log10（Hw））;grid on;
axis（[0Fs/2-12010]）;
xlabel（‘f/Hz‘）;ylabel（‘幅度（dB）‘）;title（‘（b）低通滤波器幅频特性‘）
figure（3）
subplot（211）
t=[0:N-1]/Fs;Tp=N/Fs;
plot（tywt）;grid;
axis（[0Tp/2-11]）;xlabel（‘t/s‘）;ylabel（‘y_w（t）‘）;
title（‘（c）滤除噪声后的信号波形‘）
subplot（212）
YWn=abs（fft（ywt1024））;
f=[0:1023]*Fs/1024;
stem（fYWn‘.‘）;grid
axis（[0Fs/201.2*max（YWn）]）;
xlabel（‘f‘）;ylabel（‘|Y_w（n）|‘）;title（‘（d）输出信号频谱‘）
%======================等波纹最佳逼近========================
fb=[fpfs];m=[10]; 		% 确定remezord函数所需参数fmdev
dev=[（10^（Rp/20）-1）/（10^（Rp/20）+1）10^（-As/20）];
[MefomoW]=remezord（fbmdevFs）;	% 确定remez函数所需参数
hn=remez（MefomoW）;		% 调用remez函数进行设计
Hw=abs（fft（hn1024））;		% 求设计的滤波器频率特性
yet=fftfilt（hnxtN）;       % 调用函数fftfilt对xt滤波
YEn=fft（yet1024）;
%画图
figure（4）;
subplot（211）
n=0:Me;stem（nhn‘.‘）;grid on;
axis（[0M+11.2*min（hn）1.2*max（hn）]）;
xlabel（‘n‘）;ylabel（‘h（n）‘）;title（‘（e）单位脉冲响应‘）
subplot（212）
f=[0:1023]*Fs/1024;
plot（f20*log10（Hw））;grid;title（‘（f）低通滤波器幅频特性‘）
axis（[0Fs/2-8010]）;
xlabel（‘f/Hz‘）;ylabel（‘幅度‘）
figure（5）
subplot（211）;
plot（tyet）;grid;
axis（[0Tp/2-11]）;xlabel（‘t/s‘）;ylabel（‘y_e（t）‘）;
title（‘（g）滤除噪声后的信号波形‘）
subplot（212）
stem（fYEn‘.‘）;grid;
xlabel（‘f‘）;ylabel（‘|Y_e（n）|‘）;title（‘（h）输出信号频谱‘）
axis（[0Fs/201.2*max（YEn）]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1948  2014-12-12 21:01  实验5\exp5.m

     文件       1269  2014-12-08 09:45  实验5\xtg.m

     文件     222087  2014-12-22 20:28  实验5\实验五：FIR熟悉滤波器设计与软件实现.docx

     目录          0  2014-12-22 20:29  实验5

----------- ---------  ---------- -----  ----

               225304                    4