fir数字滤波器设计实现.txt

本程序是关于FIR和IIR滤波器的设计现

%FIR数字滤波器设计实现
clear all;
pi=3.1415926;

%读取数字信号的数据，打开名为"noise.txt"的信号数据文件
f_id=fopen('noise.txt','r');
[noise,count]=fscanf(f_id,'%f');

%提取时间量
j=1;
count1=count;
for i=1:2:count1
    time(j)=nois(i);
    j=j+1;
end

%提取幅度量
j=1;
for i=2:2:count1
    nois_m(j)=100*nois(i);
    j=j+1;
end

%坐标时间变换
count2=length(time);
dd=(max(time)-min(time))/2;
for k=1:count2
    time(k)=time(k)+dd;
end

%保存时域波形
nois1=[time;nois_m];
save noise.dat nois1 -ascii
X=nois1;

%画出滤波前信号波形图
plot(time,noise_m);
grid;
xlabel('时间/us');
ylabel('幅度');

n=35;    %FIR滤波器阶数
f=[0,0.01,0.1,1.0];    %FIR带通滤波器边界频率参数，0-0.001阻带 0.001-0.1通带  0.1-1.0阻带
a=[0,1,0];    %描述滤波器的理想特性参数，0-阻带 1-通带
up=[0.01,1.1,0.01];    %通带、阻带最大波动值
lo=[0,0.9,0];    %通带、阻带最小波动值
b=fircls(n,f,a,up,lo);    %计算滤波器单位脉冲响应h(n)
save b-firfilter.dat b -ascii;    %保存h(n)

sig_filter1=1*fftfilt(b,nois_m);    %实现FIR滤波，返回滤波后的信号
figure;
plot(time,sig_filter1);
sig_filter=[time,sig_filter1];
save sig_firfilter.dat sig_filter -ascii;    %保存滤波后的信号

grid;
xlabel('时间/us');
ylabel('响应幅度值');