- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
. .
实验六、FIR数字滤波器设计及其网络结构
班级: 学号: 姓名: 成绩:
1实验目的
(1)熟悉线性相位FIR数字滤波器的时域特点、频域特点和零极点分布;
(2)掌握线性相位FIR数字滤波器的窗函数设计法和频率采样设计法;
(3)了解IIR数字滤波器和FIR数字滤波器的优缺点及其适用场合。
2 实验内容
(1)设计计算机程序.根据滤波器的主要技术指标设计线性相位FIR数字低通、高通、带通和带阻滤波器;
(2)绘制滤波器的幅频特性和相频特性曲线.验证滤波器的设计结果是否达到设计指标要求;
(3)画出线性相位FIR数字滤波器的网络结构信号流图。
3实验步骤
(1)设计相应的四种滤波器的MATLAB程序;
(2)画出幅频相频特性曲线;
(3)画出信号流图。
4 程序设计
%% FIR低通
f=[0.2,0.35];
m=[1,0];
Rp=1;Rs=40;
dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);
dat2=10^(-Rs/20);
rip=[dat1,dat2];
[M,f0,m0,w]=remezord(f,m,rip);
M=M+2;
hn=remez(M,f0,m0,w);
w=0:0.001:pi
xn=[0:length(hn)-1];
H=hn*exp(-j*xn*w);
figure
subplot(2,1,1)
plot(w/pi,20*log10(abs(H)));grid on;xlabel(\omega/\pi),ylabel(|H(e^j^w)|/dB)
subplot(2,1,2)
plot(w/pi,angle(H)/pi);xlabel(\omega/\pi),ylabel(\phi(\omega)/\pi)
%% FIR高通
f=[0.7,0.9];
m=[0,1];
Rp=1;Rs=60;
dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);
dat2=10^(-Rs/20);
rip=[dat2,dat1];
[M,f0,m0,w]=remezord(f,m,rip);
hn=remez(M,f0,m0,w);
w=0:0.001:pi
xn=[0:length(hn)-1];
H=hn*exp(-j*xn*w);
figure
subplot(2,1,1)
plot(w/pi,20*log10(abs(H)));grid on;xlabel(\omega/\pi),ylabel(|H(e^j^w)|/dB)
subplot(2,1,2)
plot(w/pi,angle(H)/pi);xlabel(\omega/\pi),ylabel(\phi(\omega)/\pi)
%% FIR带通
f=[0.2,0.35,0.65,0.8];
m=[0,1,0];
Rp=1;Rs=60;
dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);
dat2=10^(-Rs/20);
rip=[dat2,dat1,dat2];
[M,f0,m0,w]=remezord(f,m,rip);
M=M+3
hn=remez(M,f0,m0,w);
w=0:0.001:pi
xn=[0:length(hn)-1];
H=hn*exp(-j*xn*w);
figure
subplot(2,1,1)
plot(w/pi,20*log10(abs(H)));grid on;xlabel(\omega/\pi),ylabel(|H(e^j^w)|/dB)
subplot(2,1,2)
plot(w/pi,angle(H)/pi);xlabel(\omega/\pi),ylabel(\phi(\omega)/\pi)
%% FIR带阻
f=[0.2,0.35,0.65,0.8];
m=[1,0,1];
Rp=1;Rs=60;
dat1=(10^(Rp/20)-1)/(10^(Rp/20)+1);
dat2=10^(-Rs/20);
rip=[dat1,dat2,dat1];
[M,f0,m0,w]=remezord(f,m,rip);
hn=remez(M,f0,m0,w);
w=0:0.001:pi
xn=[0:length(hn)-1];
H=hn*exp(-j*xn*w);
figure
subplot(2,1,1)
plot(w/pi,20*log10(abs(H)));grid on;xlabel(\omega/\pi),ylabel(|H(e^j^w)|/dB)
subplot(2,1,2)
plot(w/pi,angle(H)/pi);xlabel(\omega/\pi),y
您可能关注的文档
- 唯美渐变时尚项目汇报PPT模板.ppt
- 无肝素透析护理.ppt
- 武汉剑桥六年级检测试题.doc
- 细胞周期与调控机制.ppt
- 消防安全公开课课件.ppt
- 小数的意义和读写集体备课.doc
- 小细胞肺癌2016的指南解读.ppt
- 小学数学骨干教师培训的讲座.ppt
- 校园无线局域网毕业论文_PPT设计说明书.ppt
- 新发展理念公务员读本98分试卷答案.doc
- 【可行性报告】2023年芯类产品项目可行性研究分析报告.docx
- 兽用抗菌药研究报告-中国兽用抗菌药市场深度全景调研及投资前景分析报告2023年.pptx
- 女装高领毛衣研究预测报告_中国女装高领毛衣行业发展趋势及投资风险研究预测报告(2024-2028年).pptx
- 【可行性报告】2023年水质检测相关项目可行性研究报告.docx
- 【可行性报告】2023年消音降噪设备行业项目可行性分析报告.docx
- 电力安防研究报告-中国电力安防行业发展分析与投资前景预测报告2024年.pptx
- 【可行性报告】2023年钯相关行业可行性分析报告.docx
- 养生茶研究预测报告-中国养生茶行业供需分析及发展前景研究预测报告(2024-2028年).pptx
- 【可行性报告】2023年玻璃器皿项目可行性研究分析报告.docx
- 企业活动策划方案.pdf
文档评论(0)