MATLAB抽样定理实验

这篇具有很好参考价值的文章主要介绍了MATLAB抽样定理实验。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

目录

一、实验目的

二、实验原理

三、实验要求

四、实验内容

1、连续时间信号时域波形及其幅度谱

2、信号进行抽样

3、频谱分析

4、由各抽样信号恢复出连续时间信号,计算并画出误差函数


一、实验目的

1、掌握抽样定理工作原理

2、练习使用Matlab编程进行抽样定理验证与分析

3、分析并通过实验观察过采样、临界采样和欠采样3种不同条件下恢复信号误差,并由此总结抽样频率对信号恢复产生误差影响,加深对时域低通抽样定理的理解。

二、实验原理

抽样定理:

设时间连续信号f(t),其最高截止频率为fm,如果用时间间隔为T≤12fm的开关信号对f(t)进行抽样时,则f(t)就可被样值信号唯一地表示。

在一个频带限制在(0,fh)内的时间连续信号f(t),如果以小于等于的时间间隔对它进行抽样,那么根据这些抽样值就能完全恢复原信号。或者说,如果一个连续信号f(t)的频谱中最高频率不超过fh,这种信号必定是个周期性的信号,当抽样频率fs≥2fh时,抽样后的信号就包含原连续信号的全部信息,而不会有信息丢失,当需要时,可以根据这些抽样信号的样本来还原原来的连续信号。

时域低通抽样定理:

奈奎斯特抽样定理:设有一个频带限制在(0,fh)Hz内的时间连续信号f (t),如果以不低于2fh次/秒的频率对它进行抽样,那么所得的抽样值将包含f (t)的全部信息,并且可以用低通滤波器从这些样值中重建f (t)。假设f (t)的频谱为F(w),我们抽样所用的信号是单位冲击序列:

其中:Ts为抽样时间间隔,那么抽样后的信号fs(t)为:

其信号频谱为:

抽样后信号f (t)的频谱 由无限多个以ωs的各次谐波为中心点所组成,当然幅度只有原来的1/Ts

显然为了要使相邻的边带不发生混叠,必须满足如下条件ωs≥2ωh,或fs≥2fh

当抽样满足抽样定理要求,频谱不发生混叠时,在接收端只要用理想低通滤波器就可以从抽样信号中无失真地恢复原信号。

带通信号抽样定理:

f(t)频带为,仍按fs=2fh抽样。

经观察可发现带通信号的最高频率fh如果是其带宽的整数倍的话,例如fh=2B,当抽样频率fs=2(fsfl )=2B时,其频谱并不发生混叠。

如果最高频率fh不是信号带宽B的整数倍,即:fh=KB

其中K的整数部分为n,小数部分为k,即:K=n+k

我们可以假想一个比B宽的带宽B',使正好是它的整数倍。fh=nB'

只要我们以2B'抽样频率fhf (t)进行抽样必然不会出现频谱混叠。因此

MATLAB抽样定理实验

从式可见,随着n的增大,趋向于2B,当n比较大时,式可简化为:fk=2B

三、实验要求

1、画出某连续时间信号时域波形及其幅度谱。

如:MATLAB抽样定理实验

或者自己编写一信号,或者找某一语音信号

2、对信号进行抽样(可以理想抽样,可以周期窄矩形脉冲抽样),分别画出抽样频率为200Hz、300Hz、500Hz时抽样后信号时域波形。

3、对上述3种抽样频率下得到的抽样信号进行频谱分析,绘制幅度谱曲线,对比各抽样频率下抽样信号时域波形和幅度谱有无差异。

4、由各抽样信号恢复出连续时间信号,分别画出时域波形,计算并画出误差函数,并与原连续信号波形进行对比,总结不同抽样频率对信号恢复产生的影响。

四、实验内容

1、连续时间信号时域波形及其幅度谱

t=0:1/2000:0.1;
t1=0:1/200:2;
t2=0:1/300:2;
t3=0:1/500:2;

f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
f1=cos(100*pi*t1)+2*sin(200*pi*t1)+3*cos(300*pi*t1);
f2=cos(100*pi*t2)+2*sin(200*pi*t2)+3*cos(300*pi*t2);
f3=cos(100*pi*t3)+2*sin(200*pi*t3)+3*cos(300*pi*t3);

subplot(121)
plot(t,f);
xlabel('t');
title('f(t)的时域波形');


dt=0.0005;
N=10000;
k=-N:N;
W=pi*k/(N*dt);
subplot(122)
t = -1:dt:1;f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
F=dt*f*exp(-j*t'*W);
F=abs(F);
plot(W,F),grid on,
title('f(t)的幅度谱');

MATLAB抽样定理实验

我们可以看到时域f(t)的时域波形图是一个周期函数。观察频谱,可以发现频谱结构是关于原点对称的谱线,因为是信号中有三种不同频率的正弦信号,其傅里叶变换是三种大小不一样的冲激函数。 

2、信号进行抽样

2.1理想抽样时域波形

t=0:1/20000:0.1;
t1=0:1/200:0.1;
t2=0:1/300:0.1;
t3=0:1/500:0.1;

f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
f1=cos(100*pi*t1)+2*sin(200*pi*t1)+3*cos(300*pi*t1);
f2=cos(100*pi*t2)+2*sin(200*pi*t2)+3*cos(300*pi*t2);
f3=cos(100*pi*t3)+2*sin(200*pi*t3)+3*cos(300*pi*t3);

MATLAB抽样定理实验

 理想抽样可以得到一系列周期的抽样信号

3、频谱分析

3.1理想抽样频谱分析

设有一个最高频率小于的信号,将这个信号和周期性单位冲激脉冲相乘,得到的就是理想抽样信号。即:

 

从而可以得到

figure(2),

N=10000;
k=-N:N;
W=pi*k/(N*dt);

dt=1/10000;t = -1:dt:1;f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
F=dt*f*exp(-j*t'*W);
F=abs(F);
subplot(221),plot(W,F),grid on,axis([-350*pi 350*pi 0 1]);
title('f(t)的幅度谱');

dt=1/200;t = -1:dt:1;f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
F=dt*f*exp(-j*t'*W);
F=abs(F);
subplot(222),plot(W,F),grid on,axis([-350*pi 350*pi 0 1]);
title('f(t)200Hz理想采样的幅度谱');

dt=1/300;t = -1:dt:1;f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
F=dt*f*exp(-j*t'*W);
F=abs(F);
subplot(223),plot(W,F),grid on,axis([-350*pi 350*pi 0 1]);
title('f(t)300Hz理想采样的幅度谱');

dt=1/500;t = -1:dt:1;f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
F=dt*f*exp(-j*t'*W);
F=abs(F);
subplot(224),plot(W,F),grid on,axis([-350*pi 350*pi 0 1]);
title('f(t)500Hz理想采样的幅度谱');

MATLAB抽样定理实验

当抽样频率为200Hz、300Hz的时候发生了频谱混叠,因为抽样会造成信号的频域频谱搬移,但是因为采样信号频率,就会造成频谱的混叠,经过低通滤波器无法还原原本的信号。但是观察500Hz的采样信号,发现频谱搬移后没有发生混叠,经过低通滤波器可以复原信号。 

4、由各抽样信号恢复出连续时间信号,计算并画出误差函数

连续信号可以用抽样函数的无穷级数来表示,该级数的系数就是抽样值。

t=0:1/20000:2;
t1=0:1/200:2;
t2=0:1/300:2;
t3=0:1/500:2;

f=cos(100*pi*t)+2*sin(200*pi*t)+3*cos(300*pi*t);
f1=cos(100*pi*t1)+2*sin(200*pi*t1)+3*cos(300*pi*t1);
f2=cos(100*pi*t2)+2*sin(200*pi*t2)+3*cos(300*pi*t2);
f3=cos(100*pi*t3)+2*sin(200*pi*t3)+3*cos(300*pi*t3);

%信号重建

subplot(311)
dt=1/200;
ft=f1*dt*(310*pi)/pi*sinc(310*(ones(length(t1),1)*t-t1'*ones(1,length(t))));
error=abs(ft-f);
plot(t,error),grid on;
xlabel('\itt'),ylabel('\ite\rm(\itt\rm)');
title('\fontname{宋体}200Hz抽样重建信号与原信号的绝对误差');
set(gca,'Fontname', 'Times New Roman', 'FontSize',12);

%信号重建

subplot(312)
dt=1/300;
ft=f2*dt*(310*pi)/pi*sinc(310*(ones(length(t2),1)*t-t2'*ones(1,length(t))));
error=abs(ft-f);
plot(t,error),grid on;
xlabel('\itt'),ylabel('\ite\rm(\itt\rm)');
title('\fontname{宋体}300Hz抽样重建信号与原信号的绝对误差');
set(gca,'Fontname', 'Times New Roman', 'FontSize',12);

%信号重建
subplot(313)
dt=1/500;
ft=f3*dt*(310*pi)/pi*sinc(310*(ones(length(t3),1)*t-t3'*ones(1,length(t))));
error=abs(ft-f);
plot(t,error),grid on;
xlabel('\itt'),ylabel('\ite\rm(\itt\rm)');
title('\fontname{宋体}500Hz抽样重建信号与原信号的绝对误差');
set(gca,'Fontname', 'Times New Roman', 'FontSize',12);

MATLAB抽样定理实验

MATLAB抽样定理实验

观察四个时域波形,可以发现在200Hz和300Hz抽样时和原信号差距比较大,500Hz抽样时和原信号相差很小。我们绘制误差函数可以得到如下结果,因为我们只是截取了部分信号,所以不可避免会有频谱泄露现象,在信号开始和结束附近误差较大。总结:只有当采样频率时,才能保证信号正常恢复,否则会造成频谱混叠。文章来源地址https://www.toymoban.com/news/detail-447481.html

到了这里,关于MATLAB抽样定理实验的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 【通信原理(含matlab程序)】实验二:FM的调制和解调

    💥💥💞💞欢迎来到本博客❤️❤️💥💥 本人持续分享更多关于电子通信专业内容以及嵌入式和单片机的知识,如果大家喜欢,别忘点个赞加个关注哦,让我们一起共同进步~ 理解FM基本原理、FM信号时域和频域主要特点; 理解FM信号调制和解调的实现方法; 掌握matlab程序

    2024年02月08日
    浏览(45)
  • Matlab 离散小波变换函数 dwt2() 原理介绍与实验

    离散小波变换(DWT)的原理介绍和说明请参考文章: 【DWT笔记】傅里叶变换与小波变换 这篇文章写的通俗易懂,小白也能看懂。 1.1 小波变换简介 离散小波变换(DWT)的原理介绍和说明请参考文章:【DWT笔记】傅里叶变换与小波变换 这篇文章写的通俗易懂,小白也能看懂。

    2024年02月04日
    浏览(44)
  • 物联网控制原理与技术--基于Matlab/利用MATLAB进行频域分析(伯德图)的应用(超详细/设计/实验/作业/练习)

    (1)熟练掌握运用MATLAB命令绘制控制系统伯德图的方法; (2)了解系统伯德图的一般规律及其频域指标的获取方法; (3)熟练掌握运用伯德图分析控制系统稳定性的方法; 1、Windows 10 2、Matlab 2012a 1. 用MATLAB作伯德图 控制系统工具箱里提供的bode()函数可以直接求取、绘制给

    2024年02月09日
    浏览(45)
  • 【Spring Boot】缓存预热与缓存清除--原理、作用和目的

    在本文中,我们将详细讲解两个与缓存相关的重要概念: 缓存预热 和 缓存清除 。 我们将了解缓存的基本概念,然后分别介绍缓存预热和缓存清除的原理、作用和目的。 缓存是一种数据存储技术,用于存储经常访问的数据,以便在需要时快速获取 。通过缓存数据,可以 减

    2023年04月21日
    浏览(56)
  • 【数字信号处理】带通采样定理及其MATLAB仿真

    按照奈奎斯特采样定理(低通采样),采样频率 f s f_{s} f s ​ 要大于等于信号中最高频率 f m a x f_{max} f ma x ​ 的2倍,才可以保证采样后的数字信号通过DAC转换后,可以无失真的恢复为原信号。然而,如果信号的频率分布在某一有限频带上,并且信号的最高频率 f m a x f_{max} f

    2024年02月16日
    浏览(49)
  • 南京邮电大学电工电子基础B实验四(戴维南与诺顿定理)

    一、 实验目的 1、学习几种常用的等效电源的测量方法 2、比较几种测量方法所适用的情况 3、分析各种方法的误差大小及其产生的原因 二、 主要仪器设备及软件 硬件:交流电源、电容、电感、电阻、波特图仪。 软件:Multisim14.0 三、 75页实验表格 四、 仿真电路 五、 测量方

    2023年04月15日
    浏览(58)
  • 【AHP层次分析法】原理+应用步骤+旅游目的地选择实例应用

    实现标准之间相对重要程度,并给出决策方案中每个标准的权重,利用权数求出各个方案的优劣次序。 1.建立层次结构模型 该层主要有三个方面: 目标层 准则层 领域层(各种解决问题的措施和方案) 这里选择了一个旅游问题的层次分析模型来直观的展示三个层的关系: 如

    2024年02月04日
    浏览(44)
  • 实验6-cp –r系统命令的实现--源路径(目录)中的所有文件和子目录,以及子目录中的所有内容,全部拷贝到目标路径(目录)中--操作系统实验

    掌握Linux目录操作方法,包括打开目录、关闭目录、读取目录文件 掌握Linux文件属性获取方法,包括三个获取Linux文件属性的函数、文件属性解析相关的宏 掌握POSIX与ANSI C文件I/O操作方法,包括打开文件、关闭文件、创建文件、读写文件、定位文件 利用POSIX API(文件操作也可

    2024年02月08日
    浏览(54)
  • 概率统计·样本及抽样分布【随机样本、抽样分布】

    总体:试验的全部可能的观察值称为总体 个体:总体中每个可能的观察值称为个体 总体期望=样本平均期望 总体方差/n=样本平均方差 X 1 ,X 2 ……X n 相互独立(x 1 ,x 2 ……x n 是观察值),称为总体X的一个简单随机变量(样本) 联合=(全部)边缘相乘 函数表示化(不含未

    2024年02月10日
    浏览(57)
  • 【burpsuite安全练兵场-服务端3】目录遍历漏洞-6个实验(全)

       博主:网络安全领域狂热爱好者(承诺在CSDN永久无偿分享文章)。 殊荣:CSDN网络安全领域优质创作者,2022年双十一业务安全保卫战-某厂第一名,某厂特邀数字业务安全研究员,edusrc高白帽,vulfocus、攻防世界等平台排名100+、高校漏洞证书、cnvd原创漏洞证书等。 擅长:

    2024年02月03日
    浏览(42)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包