fiirandiir.txt

fiir和iir滤波器设计源码

%IIR滤波
t=0:0.001:0.5;
x=sin(50*2*pi*t);
z=randn(size(t));
y=x+z;
figure(1),subplot(221);
plot(t(1:50),x(1:50),'r')
title('原图');
xlabel('时间t');
ylabel('幅值');
figure(1),subplot(223);
plot(t(1:50),y(1:50),'b')
title('原图加噪声');
xlabel('时间t');
ylabel('幅值');
figure(1),subplot(222);
plot(t(1:50),z(1:50),'c')
title('噪声');
xlabel('时间t');
ylabel('幅值');
%***************************
Y = fft(y,512);
Pyy = Y.* conj(Y) / 512;
f = 1000*(0:256)/512;
figure(2),plot(f,Pyy(1:257),'m')%频谱图
title('加噪声后的频谱图');
xlabel('频率(Hz)')；
ylabel('幅值y')
%FIR滤波
t=0:0.001:0.5;%采样频率为1000
x=sin(50*2*pi*t);
subplot(221);
plot(t,x,'y');
z=randn(size(t));
y=x+z;
subplot(222);
plot(1000*t(1:50),y(1:50),'g');%只显示前面0.05秒
title('加随机噪声');
xlabel('时间t');
Y = fft(y,512);%快速傅立叶变换
Pyy = Y.* conj(Y) / 512;%平均幅值
f =1000*(0:256)/512;
subplot(223);
plot(f,Pyy(1:257),'r')%y的频谱图
title('加噪声后的频谱图');
xlabel('频率(Hz)');
[n,Wn] = buttord([40 60]/500,[20 80]/500,3,50);
b=fir1(80,[40 60]/500);%滤波器的阶数为80
figure(2)
freqz(b,1,512,1000);%图7,8，512个样本，采样频率1000H(z)
y = filter(b,1,y);%1为分母系数b为H(z)分子系数
figure(3),plot(t(1:500),y(1:500));