- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
山东建筑大学信电学院课程设计说明书
PAGE
PAGE I
山东建筑大学
课 程 设 计 说 明 书
题 目: 基于MATLAB的信号分析与处理
课 程: 数字信号处理课程设计
院 (部): 信息与电气工程学院
专 业: 通信工程
班 级: 通信111班
学生姓名:
学 号:
指导教师:
完成日期: 2014年1月
目 录 TOC \o 1-3 \h \z \u
摘 要 Ⅰ
1 设计目的和要求 1
2 设计原理 2
3 设计内容 3
3.1 程序源代码 4
3.2 调试分析与过程描述 7
3.3 结果分析 12
总 结 13
致 谢 14
参考文献 15
摘 要
这次是基于MATLAB的信号分析与处理。所谓数字滤波器,就是输入、输出都是数字信号的,通过数值计算处理改变输入信号所含频率成分的相对比例,或者滤除某些频率成分的数字器件或程序。常用的经典滤波器有低通、高通、带通、带阻。
首先产生一个连续信号,包含低频、中频、高频分量;对其进行采样,得到数字信号;对数字信号进行FFT频谱分析,绘制其频谱图;根据信号频谱分析的结果,分别设计高通、低通、带通滤波器,绘制滤波器的幅频及相频特性;用所设计的滤波器对信号滤波,并绘制出滤波后的频谱图。
关键词:MATLAB; FFT;滤波器;信号产生;频谱分析
1设计目的和要求
产生一个连续信号,包含低频,中频,高频分量,对其进行采样,进行频谱分析,分别设计三种高通,低通,带通滤波器对信号进行滤波处理,观察滤波后信号的频谱。
2设计原理
信号的采样要符合奈奎斯特采样定律,一般为被采信号最高频率的2倍,只有这样,才能保证频域不混叠,也就是采样出来数字信号中包含了被采信号的所有信息,而且没有引入干扰。这就是信号的时域采样。
频谱分析是指对信号进行频域谱的分析,观察其频域的各个分量的功率大小,其理论基础是傅立叶变换,现在一般采用数字的方法,也就是将时域信号数字化后做FFT,可以得到频域的波形。
数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。可以设计系统的频率响应,让它满足一定的要求,从而对通过该系统的信号的某些特定的频率成分进行过滤,这就是滤波器的基本原理。
IIR滤波器的设计原理:
IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法来进行设计,通常采用模拟滤波器原型有butterworth函数、chebyshev函数、bessel函数、椭圆滤波器函数等。
IIR数字滤波器的设计步骤:
(1)按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标;
(2)根据模拟滤波器技术指标设计为响应的模拟低通滤波器;
(3)很据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器;
(4)如果要设计的滤波器是高通、带通或带阻滤波器,则首先把它们的技术指标转化为模拟低通滤波器的技术指标,设计为数字低通滤波器,最后通过频率转换的方法来得到所要的滤波器。
本课程设计设计思想:首先利用MATLAB分别产生低频、中频、高频信号,然后进行叠加得到连续时间信号;对所产生的连续时间信号进行采样,得到数字信号;对信号进行FFT频谱分析,绘制其频谱图;根据信号频谱分析的结果,分别设计高通,低通,带通滤波器,得到滤波器的幅频及相频特性。
设计内容
3.1程序源代码
f1=40;
f2=300;
f3=2000;
t=(1:500)/8000;%每个周期采8000个点,Fs1=8000
x1=sin(2*pi*t*f1);
figure(1);
subplot(2,1,1);
plot(x1);
xlabel(t);
ylabel(x1(t));
title(连续信号);
grid;
x2=sin(2*pi*t*f2);
subplot(2,1,2);
plot(x2);
xlabel(t);
ylabel(x2(t));
title(连续信号);
grid;
x3=sin(2*pi*t*f3);
figure(2);
subplot(2,1,1);
plot(x3);
xlabel(t);
ylabel(x3(t));
title(连续信号);
grid;
subplot(2,1,2); %只放大显示x3(t)的0~50的部分
plot(x3);
axis([0 50 -1 1.2]);
xlabel(t);
ylabel(x3(t));
title(连续信号);
grid;
x=sin(2*pi*t*f1)+sin(2*pi*t*f2)+sin(2*pi*
文档评论(0)