- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
基于Matlab的FIR滤波器设计与实现-sunev-博客园
sunev
博客园
闪存
首页
新随笔
联系
管理
订阅
随笔- 43? 文章- 0? 评论- 162?
前面一篇文章介绍了通过FDATool工具箱实现滤波器的设计,见“基于Matlab中FDATool工具箱的滤波器设计及相关文件的生成”,这里通过几个例子说明采用Matlab语言设计FIR滤波器的过程。
Matlab7.1
以低通滤波器为例,其常用的设计指标有:
通带边缘频率fp(数字频率为Ωp)
阻带边缘频率fst?(数字频率为Ωst)
通带内最大纹波衰减δp=-20log10(1-αp),单位为 dB
阻带最小衰减αs=-20log10(αs),单位为 dB
阻带起伏αs
通带峰值起伏αp
其中,以1、2、3、4条最为常用。5、6条在程序中估算滤波器阶数等参数时会用到。
数字频率 = 模拟频率/采样频率
第一步:通过过渡带宽度和阻带衰减,计算滤波器的阶数B和β值。
第二步:通过程序设计滤波器。
程序如下:
b?=?fir1(29,0.4,kaiser(30,4.55));
[h1,w1]=freqz(b,1);??
plot(w1/pi,20*log10(abs(h1)));?
axis([0,1,-80,10]);?
grid;
xlabel('归一化频率/p')?;
ylabel('幅度/dB')?;
?
波形如下:
?
采用[n,Wn,beta,ftype]?=?kaiserord(f,a,dev)函数来估计滤波器阶数等,得到凯塞窗滤波器。
这里的函数kaiserord(f,a,dev)或者kaiserord(f,a,dev,fs):
f为对应的频率,fs为采样频率;当f用数字频率表示时,fs则不需要写。
a=[1?0]为由f指定的各个频带上的幅值向量,一般只有0和1表示;a和f长度关系为(2*a的长度)- 2=(f的长度)
devs=[0.05?10^(-2.5)]用于指定各个频带输出滤波器的频率响应与其期望幅值之间的最大输出误差或偏差,长度与a相等,计算公式:
阻带衰减误差=αs,通带衰减误差=αp,可有滤波器指标中的3、4条得到。
fs缺省为2Hz。
?
程序如下:
fcuts?=?[0.3??0.5];?%归一化频率omega/pi,这里指通带截止频率、阻带起始频率
mags?=?[1?0];
devs?=?[0.05?10^(-2.5)];
[n,Wn,beta,ftype]?=?kaiserord(fcuts,mags,devs);??%计算出凯塞窗N,beta的值
hh?=?fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');?
freqz(hh);
?
波形如下:
?
实际中,一般调用MATLAB信号处理工具箱函数remezord来计算等波纹滤波器阶数N和加权函数W(ω),调用函数remez可进行等波纹滤波器的设计,直接求出滤波器系数。函数remezord中的数组fedge为通带和阻带边界频率,数组mval是两个边界处的幅值,而数组dev是通带和阻带的波动,fs是采样频率单位为Hz。
?
?
解:在MATLAB中可以用remezord?和remez两个函数设计
?
程序如下:?
?
fedge=[800?1000];?
?
mval=[1?0];?
?
dev=[0.0559?0.01];
?
fs=4000;?
?
[N,fpts,mag,wt]=remezord(fedge,mval,dev,fs);?
?
b=remez(N,fpts,mag,wt);?
?
[h,w]=freqz(b,1,256);?
?
plot(w*2000/pi,20*log10(abs(h)));?
?
grid;?
?
xlabel('频率/Hz')?;
?
ylabel('幅度/dB');
?
波形如下:
?
程序如下:
[n,wn,bta,ftype]=kaiserord([0.3?0.45?0.65?0.8],[0?1?0],[0.01?0.1087?0.01]);%用kaiserord函数估计出滤波器阶数n和beta参数
h1=fir1(n,wn,ftype,kaiser(n+1,bta),'noscale');
[hh1,w1]=freqz(h1,1,256);
figure(1)
subplot(2,1,1)
plot(w1/pi,20*log10(abs(hh1)))
grid
xlabel('归一化频率w');ylabel('幅度/db');
subplot(2,1,2)
plot(w1/pi,angle(hh1))
grid
xlabel('归一化频率w');ylabel('相位/rad');
?
您可能关注的文档
- 城市防洪规划编制大纲.doc
- 培养创新观念培育创造能力铸造创新型人才.doc
- 培优提高十一.doc
- 培训导师协议.doc
- 培训需求分析.doc
- 基于ARCSDE的版本管理与历史回溯的设计与实现.doc
- 基于AR模型谱估计方法及其应用仿真.doc
- 基于CAN总线的纯电动汽车电池管理系统研究【终稿】.doc
- 基于C语言的RS232串行接口通信实现.doc
- 基于dSPACE的三容液位控制系统仿真研究-迟长婷.doc
- 给荔红简约中国风员工述职报告.pptx
- 2023门窗行业招商加盟数据分析报告.pdf
- 《基于基因表达数据的胰腺癌预后研究现状综述》3900字.docx
- 《连云港美致医疗器械公司员工激励问题及对策分析》9300字.docx
- 《惠州名道时尚武夷岩茶股份有限公司企业文化建设研究》9000字.docx
- 2024年02月中国共产党崇左市委员会党史研究室(广西)2024年招考1名编外工作人员笔试参考题库后附答案详解.docx
- 2024年02月中共清远市委办公室(广东)2024年公开选调2名公务员笔试参考题库后附答案详解.docx
- 2024年02月中共雅安市雨城区委组织部2024年公开考核招考4名高层次人才笔试参考题库后附答案详解.docx
- 2024年02月中国基本建设优化研究会2024年招考笔试参考题库后附答案详解.docx
- 2024年02月中国伊斯兰教协会2024年度面向社会公开招考2名工作人员笔试参考题库后附答案详解.docx
文档评论(0)