域上分别设计低通和高通滤波器程序-第1节

按键盘上方向键 ← 或 → 可快速上下翻页，按键盘上的 Enter 键可回到本书目录页，按键盘上方向键 ↑ 可回到本页顶部！
————未阅读完？加入书签已便下次继续阅读！


　悬赏分：200　…　解决时间：2008…11…29　23：14　
写四个程序　时域上的低通与高通
　　　　　　　　　　　频域上的低通与高通滤波器程序　
问题补充：急于交试验报告　　还请哪位高手帮帮忙　　参数没有要求可以自己设置

提问者：　匿名最佳答案％设计低通滤波器：

＇N；Wc＇=buttord（）

％估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc

＇a；b＇=butter（N；Wc）；　％设计Butterworth低通滤波器

＇h；f＇=freqz（）；　％求数字低通滤波器的频率响应

figure（2）；　％　打开窗口2

subplot（221）；　％图形显示分割窗口

plot（f；abs（h））；　％绘制Butterworth低通滤波器的幅频响应图

title（巴氏低通滤波器''）；

grid；　％绘制带网格的图像

sf=filter（a；b；s）；　％叠加函数S经过低通滤波器以后的新函数

subplot（222）；

plot（t；sf）；　％绘制叠加函数S经过低通滤波器以后的时域图形

xlabel（'时间　（seconds）'）；

ylabel（'时间按幅度'）；

SF=fft（sf；256）；　％对叠加函数S经过低通滤波器以后的新函数进行256点的基—2快速傅立叶变换

w=　％新信号角频率

subplot（223）；

plot（））；　％绘制叠加函数S经过低通滤波器以后的频谱图

title（'低通滤波后的频谱图'）；

％设计高通滤波器

＇N；Wc＇=buttord（）

％估算得到Butterworth高通滤波器的最小阶数N和3dB截止频率Wc

＇a；b＇=butter（N；Wc；'high'）；　％设计Butterworth高通滤波器

＇h；f＇=freqz（）；　％求数字高通滤波器的频率响应

figure（3）；

subplot（221）；

plot（））；　％绘制Butterworth高通滤波器的幅频响应图

title（'巴氏高通滤波器'）；

grid；　％绘制带网格的图像

sf=filter（）；　％叠加函数S经过高通滤波器以后的新函数

subplot（222）；

plot（t；sf）；　；％绘制叠加函数S经过高通滤波器以后的时域图形

xlabel（'Time（seconds）'）；

ylabel（'Time　waveform'）；

w；　％新信号角频率

subplot（223）；

plot（））；　％绘制叠加函数S经过高通滤波器以后的频谱图

title（'高通滤波后的频谱图'）；

％设计带通滤波器

＇N；Wc＇=buttord（＇）

％估算得到Butterworth带通滤波器的最小阶数N和3dB截止频率Wc

＇a；b＇=butter（N；Wc）；　％设计Butterworth带通滤波器

＇h；f＇=freqz（）；　％求数字带通滤波器的频率响应

figure（4）；

subplot（221）；

plot（f；abs（h））；　％绘制Butterworth带通滤波器的幅频响应图

title（'butter　bandpass　filter'）；

grid；　％绘制带网格的图像

sf=filter（a；b；s）；　％叠加函数S经过带通滤波器以后的新函数

subplot（222）；

plot（t；sf）；　％绘制叠加函数S经过带通滤波器以后的时域图形

xlabel（'Time（seconds）'）；

ylabel（'Time　waveform'）；

SF=fft（）；　％对叠加函数S经过带通滤波器以后的新函数进行256点的基—2快速傅立叶变换

w=（　％新信号角频率

subplot（223）；

plot（'））；　％绘制叠加函数S经过带通滤波器以后的频谱图

title（'带通滤波后的频谱图'）；