Matlab编程实现数字滤波器怎么设计

毕业设计

你好好查点资料，这个有太多东西可以写，我给你个简单的例子：
Wc=0.5*pi;
N=21;
window=hanning(N+1); %选择要用的窗函数
b=fir1(N,Wc/pi,window); %加窗的线性相位设计函数
[H,W]=freqz(b,1,512); %求幅频响应
subplot(3,1,1)
plot(W/pi,20*log10(abs(H))); %画图增益
xlabel('Frequency(piHz)');
ylabel('Hanning');
window1=hamming(N+1);
b1=fir1(N,Wc/pi,window1);
[H,W]=freqz(b1,1,512); %求幅频响应
subplot(3,1,2)
plot(W/pi,20*log10(abs(H))); %画图增益
xlabel('Frequency(piHz)');
ylabel('Hamming');
window2=blackman(N+1);
b2=fir1(N,Wc/pi,window2);
[H,W]=freqz(b2,1,512); %求幅频响应
subplot(3,1,3)
plot(W/pi,20*log10(abs(H))); %画图增益
xlabel('Frequency(piHz)');
ylabel('Blackman');

数字滤波器是具有一定传输选择特性的数字信号处理装置,其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。它的基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高通、带通、带阻等类型,与模拟滤波器相比,数字滤波器除了具有数字信号处理的固有优点外,还有滤波精度高(与系统字长有关)、稳定性好(仅运行在0与l两个电平状态)、灵活性强等优点。数字滤波器按单位脉冲响应的性质可分为无限长单位脉冲响应滤波器iir和有限长单位脉冲响应滤波器(fir)两种。本文介绍(iir)数字滤波器的设计与分析。 1 iir数字滤波器设计思路与步骤 iir 数字滤波器可用一个n阶差分方程 y(n)=σbrx(n-r)+σaky(n-k), 或用它的z域系统函数: 对照模拟滤波器的传递函数: 不难看出,数字滤波器与模拟滤波器的设计思路相仿,其设计实质也是寻找一组系数{b,a},去逼近所要求的频率响应,使其在性能上满足预定的技术要求;不同的是模拟滤波器的设计是在s平面上用数学逼近法去寻找近似的所需特性h(s),而数字滤波器则是在z平面寻找合适的h(z)。iir数字滤波器的单位响应是无限长的,而模拟滤波器一般都具有无限长的单位脉冲响应,因此与模拟滤波器相匹配。由于模拟滤波器的设计在理论上已十分成熟,因此数字滤波器设计的关键是将h(s)→h(z),即,利用复值映射将模拟滤波器离散化。已经证明,冲击响应不变法和双线性变换法能较好地担当此任,则在此基础上,数字滤波器的设计就可首先归结为模拟滤波器的设计了。