数字信号处理课程设计---基于matlab的滤波器的设计

数字信号处理课程设计---基于matlab的滤波器的设计内容摘要：

率， wn 巴特沃兹（滤波器）的自然频率 [b,a]=butter(n,wn,39。 s39。 )。 %butter 函数是求 Butterworth 数字滤波器的系数，在求出系数后对信号进行滤波时用 filter 函数 [num,den]=bilinear(b,a,100000)。 %对模拟滤波器进行双线性变换，求得数字滤波器的传输函数系数 [h,w]=freqz(num,den)。 %数字滤波器的频率响应，返回量 H 则包含了离散系统频响 在 0—— pi 范围内 N 个频率等分点的值，向量 w 则包含 范围内 N 个频率等分点 f=w/pi*50000。 figure(5)。 plot(f,abs(h))。 axis([0,10000,0,])。 grid。 xlabel(39。 频率 /Hz39。 )。 ylabel(39。 幅度 /dB39。 )。 title(39。 LPF 的频率响应 39。 )。 ② 带 通滤波器的性能指标： fb1=5000。 fb2=30000。 fc1=3000。 fc2=32020。 5 程序如下： fb1=5000。 fb2=30000。 fc1=3000。 fc2=32020。 As=10。 %通带最大衰减 Ap=1。 %阻带最小衰减 wb1=2*100000*tan(2*pi*fb1/(2*100000))。 wb2=2*100000*tan(2*pi*fb2/(2*100000))。 wc1=2*100000*tan(2*pi*fc1/(2*100000))。 wc2=2*100000*tan(2*pi*fc2/(2*100000))。 [n,wn]=buttord([w1 w2],[ wr1 wr2],Ap,As,39。 s39。 )。 %Wb 和 Wc 分别是通带和阻带归一化的边缘频率， wn巴特沃兹（滤波器）的自然频率 [b,a]=butter(n,wn,39。 s39。 )。 %butter 函数是求 Butterworth 数字滤波器的系数，在求出系数后对信号进行滤波时用 filter 函数 [num,den]=bilinear(b,a,100000)。 %对模拟滤波器进行双线性变换，求得数字滤波器的传输函数系数 [h,w]=freqz(num,den)。 %数字滤波器的频率响应，返回量 H 则包含了离散系统频响 在 0—— pi 范围内 N 个频率等分点的值，向量 w 则包含 范围内 N 个频率等分点 f=w/pi*50000。 figure(5)。 plot(f,abs(h))。 axis([0,50000,0,])。 grid。 xlabel(39。 频率 /Hz39。 )。 ylabel(39。 幅度 39。 )。 title(39。 带通滤波器的频率响应 39。 )。 6 ③ 高 通滤波器的性能指标： fb=30000。 fc=32020。 程序如下： fb=30000。 %通带截止频率 fc=32020。 %阻带截止频率 As=10。 %通带最大衰减 Ap=1。 %阻带最小衰减 wb=2*100000*tan(2*pi*fb/(2*100000))。 %抽样频率为 1000hz wc=2*100000*tan(2*p。