- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
山东建筑大学信电学院课程设计说明书
PAGE
15 -
山东建筑大学
课 程 设 计 说 明 书
题 目: 信号的采样与恢复、采样定理的仿真
课 程: 数字信号处理课程设计
院 (部): 信息与电气工程学院
专 业: 通信工程
班 级:
学生姓名:
学 号:
指导教师: 魏莉
完成日期: 2017年1月
目 录
TOC \o 1-3 \h \z \u 目 录 - 1 -
摘 要 - 3 -
1、设计目的与要求 - 4 -
2、设计原理 - 4 -
3、设计内容与步骤 - 5 -
3.1用MATLAB产生连续信号及其对应的频谱 - 5 -
3.2对连续信号进行抽样并产生其频谱 - 6 -
3.3通过低通滤波恢复原连续信号 - 9 -
4、总结 - 14 -
5、致谢 - 15 -
6、参考文献 - 16 -
摘 要
运用数字信号处理知识实现对信号的采样、恢复以及采样定理的仿真,可借助于MATLAB强大的运算和图形显示功能,首先生成一个连续时间信号并生成其频谱,然后对该连续信号抽样,并对采样后的频谱进行分析,最后通过设计低通滤波器滤出抽样所得频谱中多个周期中的一个周期频谱,并显示恢复后的时域连续信号。
在MATLAB中连续信号的频谱由于无法实现真正的连续,所以通过对信号扩大采样点数的方法实现近似连续,取样过程中,通过控制采样频率实现大于或小于二倍连续信号的最高频率,依此来验证采样定理。信号恢复时,滤波器的参数需要很好的设置,以实现将抽样后的信号进行滤波恢复原连续信号。
1、设计目的与要求
1.掌握利用MATLAB在数字信号处理中的基本应用,并会对结果用所学知识进行分析。
2.对连续信号进行采样,在满足采样定理和不满足采用定理两种情况下对连续信号和采样信号进行FFT频谱分析。
3.从采样信号中恢复原信号,对不同采样频率下的恢复信号进行比较分析。
4.基本要求:每组一台电脑, 电脑安装MATLAB6.5版本以上软件。
2、设计原理
将连续信号采样并进行频谱分析,其频谱产生了周期延拓,每隔一个采样频率 fs,重复出现一次。为保证采样后信号的频谱形状不失真,采样频率必须大于信号中最高频率成分的两倍,这称之为采样定理。时域采样定理从采样信号恢复原信号必需满足两个条件:
(1) 必须是带限信号,其频谱函数在 > 各处为零;(对信号的要求,即只有带限信号才能适用采样定理。)
(2) 取样频率不能过低,必须满足>2。(对取样频率的要求,即取样频率要足够大,采得的样值要足够多,才能恢复原信号。)如果采样频率大于或等于,即(为连续信号的有限频谱),则采样离散信号能无失真地恢复到原来的连续信号 。
综合以上,得采样定理:
(1)对连续信号进行等间隔采样形成采样信号,采样信号的频谱是原连续信号的频谱一采样频率为周期进行周期延拓得到的;
(2)设连续信号是带限信号,如果采样角频率大于等于2倍的最高截止频率,则采样信号通过理想低通滤波器可唯一恢复出原连续信号,否则会造成采样信号中的频谱混叠现象,不可无失真的恢复原连续信号。
以采样定理为依据设计本次实验,产生一个连续时间信号,并进行频谱分析,根据采样定理要求对所产生的连续时间信号进行采样和频谱分析,并将此频谱与连续信号的频谱进行比较。验证采样定理,改变采样频率,重复以上过程。设计低通滤波器,恢复原连续信号。
3、设计内容与步骤
3.1用MATLAB产生连续信号及其对应的频谱
生成一个频率为200的连续信号。MATLAB程序源代码如下:
f1=200;t=(1:100)/3000;
w=linspace(0,100,length(t));
x1=sin(2*pi*f1*t);
n=0:(length(x1)-1);
figure(1);subplot(2,1,1);plot(x1,k);
grid on;
axis([0 100 -1.2 1.2]); %横纵坐标的范围
line([0,100],[0,0]); %画横坐标轴线
xlabel(t);ylabel(x1(t));title(连续信号);
subplot(2,1,2);
plot(w,abs(fft1(t,x1,n)),k); %其对应频域信号
xlabel(w);ylabel(x(w));
title(其对应频域信号);
grid on;
结果如下:
图3.1.1
在对其进行频谱分析的时候,需要调用以下子函数fft1:
function result=fft1(w,ha
文档评论(0)