按键盘上方向键 ← 或 → 可快速上下翻页，按键盘上的 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（'带通滤波后的频谱图'）；