



免费预览已结束,剩余1页可下载查看
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
课程名称 数字信号与处理 实验成绩 指导教师 实 验 报 告院 姓名 日期 实验六:IIR数字滤波器设计与软件实现一、实验目的1熟悉用双线性变换法设计IIR数字滤波器的原理与方法。2学会调用MATLAB信号处理工具箱中滤波器设计函数(或滤波器设计分析工具fdatool)设计各种IIR数字滤波器,学会根据滤波需求确定滤波器指标参数。3掌握IIR数字滤波器的MATLAB实现方法。4通过观察滤波器输入输出信号的时域波形及其频谱,建立数字滤波的概念二、实验内容及步骤1实验程序%计算时域离散系统损耗函数并绘图function myplot(B,A) H,W=freqz(B,A,1000);m=abs(H);plot(W/pi,20*log10(m/max(m);grid on;xlabel(omega/pi);ylabel(幅度(db);axis(0,1,-80,5);title(损耗函数曲线);%时域序列连续曲线绘图function tplot(xn,T,yn) %xn:信号数据序列;yn:绘图信号的纵坐标名称n=0:length(xn)-1;t=n*T;plot(t,xn);xlabel(t/s);ylabel(yn);axis(0,t(end),min(xn),1.2*max(xn);%产生信号序列向量st,并显示st的时域波形和频谱function st=mstg %st=mstg 返回三路调幅信号相加形成的混合信号,长度N=1600N=1600 %N为信号st的长度。Fs=10000;T=1/Fs;Tp=N*T; %采样频率Fs=10kHz,Tp为采样时间t=0:T:(N-1)*T;k=0:N-1;f=k/Tp;fc1=Fs/10; %第1路调幅信号的载波频率fc1=1000Hz,fm1=fc1/10; %第1路调幅信号的调制信号频率fm1=100Hzfc2=Fs/20; %第2路调幅信号的载波频率fc2=500Hzfm2=fc2/10; %第2路调幅信号的调制信号频率fm2=50Hzfc3=Fs/40; %第3路调幅信号的载波频率fc3=250Hz,fm3=fc3/10; %第3路调幅信号的调制信号频率fm3=25Hzxt1=cos(2*pi*fm1*t).*cos(2*pi*fc1*t); %产生第1路调幅信号xt2=cos(2*pi*fm2*t).*cos(2*pi*fc2*t); %产生第2路调幅信号xt3=cos(2*pi*fm3*t).*cos(2*pi*fc3*t); %产生第3路调幅信号st=xt1+xt2+xt3; %三路调幅信号相加fxt=fft(st,N); %计算信号st的频谱%=以下为绘图部分,绘制st的时域波形和幅频特性曲线=subplot(2,1,1)plot(t,st);grid;xlabel(t/s);ylabel(s(t);axis(0,Tp/8,min(st),max(st);title(a) s(t)的波形)subplot(2,1,2)stem(f,abs(fxt)/max(abs(fxt),.);grid;title(b) s(t)的频谱)axis(0,Fs/5,0,1.2);xlabel(f/Hz);ylabel(幅度)% IIR数字滤波器设计及软件实现clear all;close allFs=10000;T=1/Fs; %采样频率%调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st st=mstg;%低通滤波器设计与实现=fp=280;fs=450;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标(低通滤波器的通、阻带边界频)N,wp=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wpB,A=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和Ay1t=filter(B,A,st); %滤波器软件实现% 低通滤波器设计与实现绘图部分figure(2);subplot(2,1,1);myplot(B,A); %调用绘图函数myplot绘制损耗函数曲线yt=y_1(t);subplot(2,1,2);tplot(y1t,T,yt); %调用绘图函数tplot绘制滤波器输出波形%带通滤波器设计与实现=fpl=440;fpu=560;fsl=275;fsu=900;wp=2*fpl/Fs,2*fpu/Fs;ws=2*fsl/Fs,2*fsu/Fs;rp=0.1;rs=60; N,wp=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wpB,A=ellip(N,rp,rs,wp); %调用ellip计算椭圆带通DF系统函数系数向量B和Ay2t=filter(B,A,st); %滤波器软件实现figure(3);subplot(2,1,1);myplot(B,A);subplot(2,1,2);yt=y_2(t);tplot(y2t,T,yt);%高通滤波器设计与实现=fp=890;fs=600;wp=2*fp/Fs;ws=2*fs/Fs;rp=0.1;rs=60; %DF指标(低通滤波器的通、阻带边界频)N,wp=ellipord(wp,ws,rp,rs); %调用ellipord计算椭圆DF阶数N和通带截止频率wpB,A=ellip(N,rp,rs,wp,high); %调用ellip计算椭圆带通DF系统函数系数向量B和Ay3t=filter(B,A,st); %滤波器软件实现figure(4);subplot(2,1,1);myplot(B,A);subplot(2,1,2);yt=y_3(t);tplot(y3t,T,yt);2实验结果三路调幅信号st的时域波形和幅频特性曲线(a) st的时域波形(b) 低通滤波器损耗函数及其分离出的调幅信号y1(t)(c) 带通滤波器损耗函数及其分离出的调幅信号y2(t)(d)高通滤波器损耗函数及其分离出的调幅信号y3(t)3思考题(1)请阅读信号产生函数mstg,确定三路调幅信号的载波频率和调制信号频率。答:三路调幅信号的载波频率分别为250Hz、500Hz、1000Hz.调制信号频率分别为100Hz、50Hz、25Hz(2)信号产生函数mstg中采样点数N=800,对st进行N点FFT可以得到6根理想谱线。如果取N=1000,可否得到6根理想谱线?为什么?N=2000呢?请改变函数mstg中采样点数N的值,观察频谱图验证您的判断是否正确。答:分析发现,st的每个频率成分都是25Hz的整数倍。采样频率Fs=10kHz=25400Hz,即在25Hz的正弦波的1个周期中采样400点。所以,当N为400的整数
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 水产店牛蛙管理办法
- 柑橘容器苗管理办法
- 沈阳市疫情管理办法
- 望都县风控管理办法
- 武汉四类药管理办法
- 新时代农民管理办法
- 武汉充电桩管理办法
- 司机班人员管理办法
- 地铁对制度管理办法
- 县人防工程管理办法
- 2025年综合类-专业综合第十二章其他-中国邮政储蓄银行理财考试历年真题摘选带答案(5卷单选100题合辑)
- 2025年中国大唐集团有限公司应届毕业生招聘笔试历年参考题库附带答案详解
- JJG 693-2011可燃气体检测报警器
- 广东省高速公路优质工程质量管理规定宣贯
- 信息化武器装备知到章节答案智慧树2023年中北大学
- qgis制作基站扇区图层简明教程
- (重要)医疗器械经营质量管理制度及目录、工作程序
- 七年级道德与法治试卷分析
- 山东省一级保护古树名木名录
- 2023年常德市澧县中医院医护人员招聘笔试题库及答案解析
- 消化道穿孔【肠外科】课件
评论
0/150
提交评论