zongy.m

在MATLAB平台中设计一有限冲击响应滤波器

load Data_7;
%使用哈明窗设计滤波器
fp1=3;fp2=6;fs=20;
N=35;   %滤波器的阶数
w1=2*pi*fp1/fs;w2=2*pi*fp2/fs;%将模拟滤波器的指标转换为数字滤波器的指标
wind1=hamming(N);       %使用哈明窗
b1=fir1(N-1,[w1/pi,w2/pi],wind1);%使用标准频率的加窗设计函数fir1
[h1,f1]=freqz(b1,1,512,fs);      %在0到fs/2频率范围内选取512个频率点(记录在f中),计算相应的频率响应
subplot(211);plot(f1,20*log10(abs(h1))); 
xlabel('频率(KHz)');ylabel('滤波器的幅频响应曲线');grid;
subplot(212);plot(f1,unwrap(angle(h1)));
xlabel('频率(KHz)');ylabel('滤波器的相频响应曲线');grid;

figure;
[gd1 f11]=grpdelay(b1,1,512,20);
tt1=3.4;nn1=tt1*f1;
plot(nn1,gd1);grid;
xlabel('Samples(n)');ylabel('滤波器的群延迟');

figure;
n=-200:200;                
t=n/fs;
x=s+i;      %叠加之后的混合信号
sf1=filter(b1,1,x);  %利用滤波函数对混合信号进行滤波
i1=filter(b1,1,i);  %利用滤波函数对噪声信号进行滤波
subplot(5,1,1);plot(t,s);axis([0,8,-10,10]);ylabel('原始信号');grid;
subplot(5,1,2);plot(t,i);axis([0,8,-10,10]);ylabel('噪声信号');grid;
subplot(5,1,3);plot(t,x);axis([0,8,-10,10]);ylabel('混合信号');grid;
subplot(5,1,4);plot(t,sf1);axis([0,8,-10,10]);ylabel('滤波信号');grid;
subplot(5,1,5);plot(t,i1);axis([0,8,-10,10]);ylabel('滤波后的噪声部分');xlabel('单位(ms)');grid;


%使用矩形窗设计滤波器
figure;
wind2=boxcar(N);       %使用矩形窗
b2=fir1(N-1,[w1/pi,w2/pi],wind2);%使用标准频率的加窗设计函数fir1
[h2,f2]=freqz(b2,1,512,fs);      %在0到fs/2频率范围内选取512个频率点(记录在f中),计算相应的频率响应
subplot(211);plot(f2,20*log10(abs(h2))); 
xlabel('频率(KHz)');ylabel('滤波器的幅频响应曲线');grid;
subplot(212);plot(f2,unwrap(angle(h2)));
xlabel('频率(KHz)');ylabel('滤波器的相频响应曲线');grid;

figure;
[gd2 f22]=grpdelay(b2,1,512,20);
tt2=3.4;nn2=tt2*f22;
plot(nn2,gd2);grid;
xlabel('Samples(n)');ylabel('滤波器的群延迟');

figure;
n=-200:200;                
t=n/fs;
x=s+i;      %叠加之后的混合信号
sf2=filter(b2,1,x);  %利用滤波函数对混合信号进行滤波
i2=filter(b2,1,i);  %利用滤波函数对噪声信号进行滤波
subplot(5,1,1);plot(t,s);axis([0,8,-10,10]);ylabel('原始信号');grid;
subplot(5,1,2);plot(t,i);axis([0,8,-10,10]);ylabel('噪声信号');grid;
subplot(5,1,3);plot(t,x);axis([0,8,-10,10]);ylabel('混合信号');grid;
subplot(5,1,4);plot(t,sf2);axis([0,8,-10,10]);ylabel('滤波信号');grid;
subplot(5,1,5);plot(t,i2);axis([0,8,-10,10]);ylabel('滤波后的噪声部分');xlabel('单位(ms)');grid;