- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
信号抽样及抽样定理分析
实验五 信号抽样及抽样定理
一、实验目的
学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析
学会运用MATLAB改变抽样时间间隔,观察抽样后信号的频谱变化
学会运用MATLAB对抽样后的信号进行重建
实验原理
(一)信号抽样
信号抽样是利用抽样脉冲序列从连续信号中抽取一系列的离散值,通过抽样过程得到的离散值信号称为抽样信号,记为。从数学上讲,抽样过程就是信号相乘的过程,即
因此,可以使用傅里叶变换的频域卷积性质来求抽样信号的频谱。常用的抽样脉冲序列有周期矩形脉冲序列和周期冲激脉冲序列。
上式表明,信号在时域被抽样后,它的频谱是原连续信号频谱以抽样角频率为间隔周期的延拓,即信号在时域抽样或离散化,相当于频域周期化。在频谱的周期重复过程中,其频谱幅度受抽样脉冲序列的傅里叶系数加权,即被加权。
可以看出,是以为周期等幅地重复。
(二)抽样定理
如果是带限信号,带宽为,则信号可以用等间隔的抽样值来唯一表示。经过抽样后的频谱就是将的频谱在频率轴上以抽样频率为间隔进行周期延拓。因此,当时,周期延拓后频谱不会产生频率混叠;当时,周期延拓后频谱将产生频率混叠。通常把满足抽样定理要求的最低抽样频率称为奈奎斯特频率,把最大允许的抽样间隔称为奈奎斯特间隔。
(二)抽样定理
如果是带限信号,带宽为,则信号可以用等间隔的抽样值来唯一表示。经过抽样后的频谱就是将的频谱在频率轴上以抽样频率为间隔进行周期延拓。因此,当时,周期延拓后频谱不会产生频率混叠;当时,周期延拓后频谱将产生频率混叠。通常把满足抽样定理要求的最低抽样频率称为奈奎斯特频率,把最大允许的抽样间隔称为奈奎斯特间隔。
(三)信号重建
抽样定理表明,当抽样定理小于奈奎斯特间隔时,可以使用抽样信号唯一表示原信号,即信号的重建。为了从频谱中无失真的恢复原信号,可以采用截止频率为的理想低通滤波器。
上式表明连续信号可展开为抽样函数的无穷级数,该级数的系数为抽样值。
利用MATLAB中的函数来表示,所以可获得由重建的表达式,即
三、实验环境
Matlab7.0
四、实验内容
1、设有三个不同频率的正弦信号,频率分别为,,;现在使用抽样频率对这三个信号进行抽样,使用MATLAB命令画出各抽样信号的波形和频谱,并分析其频率混叠现象。
程序如下:
clear all;close all;clc;
ts=2.5*10^(-4);
dt=0.000001;
t1=0:dt:0.01;
ft=sin(2*pi*100*t1).*(t1=0);
subplot(221);
plot(t1,ft),grid on
%axis([-4 4 -0.1 1.1])
xlabel(time(sec)),ylabel(f(t))
title(100Hz正弦信号)
n=500;
k=-n:n;
w=pi*k/(n*dt);
fw=dt*ft*exp(-i*t1*w);
subplot(222)
plot(w,abs(fw)),grid on
%axis([-10 10 -0.2 1.1*pi]);
xlabel(/comega),ylabel(f(w))
title(100Hz正弦信号的频谱)
t2=0:ts:0.01;
fst=sin(2*pi*100*t2).*(t2=0);
subplot(223)
plot(t1,ft,:),hold on
stem(t2,fst),grid on
%axis([-4 4 -0.1 1.1])
title(100Hz抽样后的信号),hold off
fsw=ts*fst*exp(-i*t2*w);
subplot(224)
plot(w,abs(fsw)),grid on
%axis([-10 10 -0.2 1.1*pi])
title(100Hz抽样信号的频谱)
clear all;close all;clc;
ts=2.5*10^(-4);
dt=0.000001;
t1=0:dt:0.01;
ft=sin(2*pi*200*t1).*(t1=0);
subplot(221);
plot(t1,ft),grid on
%axis([-4 4 -0.1 1.1])
xlabel(time(sec)),ylabel(f(t))
title(200Hz正弦信号)
n=500;
k=-n:n;
w=pi*k/(n*dt);
fw=dt*ft*exp(-i*t1*w);
subplot(222)
plot(w,abs(fw)),grid on
%axis([-10 10 -0.2 1.1*pi]);
xlabel(/comega),ylabel(f(w))
title(200Hz正弦信号的频谱)
t2=0:ts:0.01;
fst=sin(2*pi*200*t2
您可能关注的文档
- 住宅楼样板引路施工方案分析.doc
- 住宅质量通病防治细则分析.doc
- 住宅楼模板施工技术交底分析.doc
- 伦敦标志性建筑分析.ppt
- 低血糖识别和预防分析.ppt
- 佐敦手机一键启动汽车系统的安装和功能分析.ppt
- 体力活动、体育与健康分析.ppt
- 体力活动测量方法综述分析.ppt
- 体外循环基础分析.ppt
- 体系结构-41309050219分析.doc
- 分析let s单元56ago2卷纸zheng unit56.pdf
- 塑胶材料其它分类原料pa9t 12.pdf
- md16x16数字媒体切换器设备.pdf
- 者参考项目发起人学科类型单位序列承包商修订页代码顺序典型.pdf
- 届世界天然气大会阿姆斯特丹2006add10288.pdf
- 期测试记录表每周weekly g1g6 journeys tests level 6 lesson26.pdf
- modernize-whitepaper现代化您应用程序白皮书.pdf
- anybackup产品典型案例分析.pdf
- 约克金融工程课程tfeslide32.pdf
- 广州市妇女儿童医疗中心历份教学药历01tjy.pdf
文档评论(0)