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

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

%IRR数字滤波器设计实现
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(X);
xlabel('时间');
ylabel('幅度');

%滤波器特性指标，根据需要可以调整
wp=10;    %通带边界频率
ws=30;    %阻带边界频率
rp=0.5;   %通带衰减
rs=85;    %阻带衰减
Fs=200;   %取样频率

%选择滤波器阶数
[N,Wn]=butter(wp/(Fs/2),ws/(FS/2),rp,rs);

%利用滤波器阶数N和等效截止频率Wn求得滤波器的系数a和b
[b,a]=butter(N,Wn);

%保存数字滤波器的系数
save b-parem.dat -ascii;
save a parem.dat -ascii;

%画滤波器幅-频特性曲线
figure;
[H,W]=freqz(b,a);
plot(W*Fs/(2*pi),abs(H));
grid;
xlabel('频率/Hz');
ylabel('幅值');

%利用系数为a和b的滤波器对信号进行滤波
Y=filter(b,a,X);

%画出滤波后的波形
figure;
plot(Y);
xlabel('时间');
ylabel('幅值');
save filter-noise.dat Y -ascii;

%画出滤波前信号X的功率谱密度
figure;
psd(X(:,2),[],Fs);
[px,fx]=psd(X(:,2),[],Fs);
pxa=[px,fx];
save x-psd.dat pxa -ascii;

%画出滤波后信号Y的功率谱密度
figure;
psd(Y(:,2),[],Fs);
[py,fy]=psd(Y(:,2),[],Fs);
pya=[py,fy];
save y-psd.dat pya -ascii;