信号的采样与插值重建(包含matlab)

实验任务

双频信号的采样和插值重建:对双频信号进行采样(符合奈奎斯特低通采样定理),绘制波形和频谱;编制sinc函数插值函数,对采样进行恢复,绘制波形。

实验原理

image-20210518225034152

image-20210518225123223

image-20210518225141718

image-20210518225150904

image-20210518225203981

matlab代码

clear all;
close all;
%%parameter settings
%采样相关参数设置
N = 100;     %采样点数
Fs = 20;     %采样频率 20Hz
Ts = 1/Fs;   %采样间隔 
DS=1/(Ts*N);        %频率分辨率
T = N*Ts;    %采样时长
n = 0:N-1; 
ts = n*Ts;%采样时间序列
f_s = -Fs/2:DS:Fs/2-DS;%采样频率序列
%原始信号相关参数设置
f1 = 1;       %信号1频率1Hz
f2 = 2;        %信号2频率2Hz
Tn = 0.01;%原始信号时间间隔
Fn = 1/Tn;
Nn = N*Ts/Tn;
Dn = Fn/Nn;
f_n = -Fn/2:Dn:Fn/2-Dn;
tn = 0:Tn:T-Tn;%原始时间序列
x = sin(2*pi*f1*tn)+sin(2*pi*f2*tn);%原始信号
x1 = sin(2*pi*f1*ts)+sin(2*pi*f2*ts);  %采样信号
figure(1);
mag=abs(fft(x))/Fn;
subplot(211);
plot(tn,x);
title(['原始信号波形']);
subplot(212);
plot(f_n,fftshift(mag));
title(['原始信号频谱']);
figure(2);
mag1=abs(fft(x1))/Fs;
subplot(211);
stem(ts,x1);
title(['采样信号波形']);
subplot(212);
plot(f_s,fftshift(mag1));
title(['采样信号频谱']);
%%sinc函数内插恢复
step = Ts/Tn;
y=zeros(1,length(tn));%y表示内插恢复后的信号
t = -T+Tn:Tn:T-Tn;
mid = ceil(length(t)/2);
h = sinc(t/Ts);
figure(3);
plot(t,h);
axis([-T,T,0,2]);
title("sinc函数波形");
for index = 1:1:length(ts)t_shift = t(mid-(index-1)*step:end-(index-1)*step);%相当于t-kTsy = y + x1(index)*sinc(t_shift/Ts);
end
figure(4);
subplot(2,1,1);
plot(tn,y);
title("sinc函数内插恢复后的信号波形");
subplot(2,1,2);
magy = abs(fft(y))/Fn;
plot(f_n,fftshift(magy));
title("sinc函数内插恢复后的信号频谱");

运行结果

f1 = 1;    %信号1频率1Hzf2 = 2;     %信号2频率2Hz

image-20210518225346781

image-20210518225404366

image-20210518225422225

image-20210518225438087

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/261480.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

dc概论之IO约束

版权声明:转载时请以超链接形式标明文章原始出处和作者信息及本声明http://bb2hh.blogbus.com/logs/39654476.html 注明:如需转载,请注明作者出处,谢谢~,Author:pythonlong以下根据资料和个人体…

对寄存器ESP和EBP的一些理解

PS:EBP是当前函数的存取指针。即存储或者读取数时的指针基地址;ESP就是当前函数的栈顶指针。每一次发生函数的调用(主函数调用子函数)时,在被调用函数初始时,都会把当前函数(主函数)…

配置FCKeditor_2.6.3+fckeditor-java-2.4

下载 http://www.fckeditor.net/download FCKeditor.Javahttp://sourceforge.net/project/showfiles.php?group_id75348&package_id129511 -1- copy fckeditor文件夹到项目的webRoot目录下 -2- 添加Jar包 slf4j-simple-1.5.0.jar slf4j-api-1.5.0.jar java-core-2.4-beta-…

普通调幅(AM)与抑制载波双边带调幅(DSB)matlab编程实现

实验任务 实现单频信号的AM、DSB调制,绘制调制前后的波形和频谱 实验原理 matlab代码 AM调制 clear all; ts 0.0025; %信号抽样时间间隔 fs 1/ts; %抽样频率 t 0:ts:10-ts; %时间向量 df fs/length(t); %fft的频率分辨率 msg2 cos(2*pi*2*t); Pn fft(msg2…

electron 读取文件夹内容_如何使用Electron Framework选择,读取,保存,删除或创建文件...

本文概述为了处理文件(CRUD)的生命周期, 我们将使用对话框和文件系统组件。对话框模块提供了用于显示本机系统对话框(例如打开文件或警报)的API, 因此Web应用程序可以提供与本机应用程序和Node.js文件系统相同的用户体验。加载所需的依赖项我们需要加载以下依赖项, 以执行我们要…

MySQL中MySQL X.X Command Line Client一闪而过的问题

问题介绍:我安装完MySQL(我安装的是5.5)后,使用MySQL 5.5 Command Line Client,每次点击,总是一闪而过。解决方法如下:首先进入cmd 切入MySQL的安装目录,然后切入 bin 目录 ,输入mysqld-nt --sk…

4pam调制与解调(matlab实现)

实验原理 这是一种使用脉冲幅度调制技术的线路编码。PAM4信号有四个电压电平,每个幅度电平分别对应逻辑比特00、01、10和11。换言之,PAM4编码的每个符号由2个比特组成,它们对应一个电压电平,即幅度。 00对应-3,01对应-1,10对应1,…

WinForm开发,窗体显示和窗体传值相关知识总结

以前对WinForm窗体显示和窗体间传值了解不是很清楚最近做了一些WinForm开发,把用到的相关知识整理如下 A.WinForm中窗体显示显示窗体可以有以下2种方法:Form.ShowDialog方法 (窗体显示为模式窗体)Form.Show方法 (窗体显示为无模式…

阿尔法蛋机器人tf卡_阿尔法蛋机器人哪款好适合几岁孩子,超能蛋早教机真实效果评测(价格309元)...

阿尔法蛋机器人是安徽淘云科技旗下的儿童智能机器人产品,作为早教机还是非常有名的,那么阿尔法蛋机器人哪款好呢?不同机器人的功能不同,价位也不等,适合的年龄段不同,这里推荐一款比较全面的阿尔法蛋超能蛋…

memset,memcpy,strcpy 的区别

一.函数原型 strcpy extern char *strcpy(char *dest,char *src); #include <string.h> 功能&#xff1a;把src所指由NULL结束的字符串复制到dest所指的数组中 说明:src和dest所指内存区域不可以重叠且dest必须有足够的空间来容纳src的字符串。 返回指向dest的指针 memcp…

普通调幅(AM)与包络检波(matlab实现)

普通调幅&#xff08;AM&#xff09;与包络检波 1.实验任务 采用matlab实现普通调幅信号的包络检波&#xff0c;并绘制包络检波后的信号的波形与频谱。 2.实验原理&#xff08;请读者参见matlab通信仿真文件夹&#xff09; 3.matlab代码 time.33;%采样时长 Ts1/10000;%采样…

免费Google地图API使用说明(转)

转自&#xff1a;http://www.cnblogs.com/mypig/articles/1936154.html01事件监视 02 03GEvent.addListener用来注册事件监视器&#xff0c;在这个例子中&#xff0c;在用户移动或拖拽地图后&#xff0c;输出地图中心点的经/纬. 04 05 06var map new GMap(document.getElem…

双栈共享存储空间(C++代码实现)

理论部分 理论来源&#xff1a;b站up主:跟懒猫老师快乐数据结构 C代码 #include<iostream> using namespace std; const int STACKSIZE 6; //两栈共享存储空间编程 //使用类模板编程 template<class DataType> class BothStack { private:DataType *data;//属性…

J storm战队成员_DOTA2J.Storm战队介绍-DOTA2ESL孟买站预选赛J.Storm战队介绍_牛游戏网攻略...

《DOTA2》ESL孟买站预选赛J.Storm战队介绍&#xff01;由ESL主办ESL孟买站各大赛区预选赛热血开战&#xff01;J.Storm战队作为本次预选赛直邀战队在比赛中会有怎样的表现呢&#xff1f;下面一起来看看小编为大家带来的J.Storm战队参赛阵容介绍吧DOTA2 ESL孟买站信息一览近日&a…

【转】HTML - embed 与 object 之争

在 HTML 里嵌入文本和图片之外的事物&#xff0c;就会用到嵌入标签&#xff0c;而嵌入标签在各浏览器之间的不统一&#xff0c;一直是让开发人员很头痛的问题。一切都要从嵌入 SUN 公司的 Applet Java 小程序开始。 当时&#xff0c;对于 Applet 的嵌入&#xff0c;浏览器老大哥…

毕业三年

2013年6月&#xff0c;毕业整整3年了。曾经的同伴&#xff0c;有人创业了&#xff0c;有人升职加薪&#xff0c;也有人职场失利。而我&#xff0c;做了2年多的独立开发者。 在大学的时候我学的是.net&#xff0c;在毕业前还去考了微软的认证&#xff0c;考过了&#xff0c;但考…

Muduo 网络编程示例之五: 测量两台机器的网络延迟

Muduo 网络编程示例之五&#xff1a; 测量两台机器的网络延迟 陈硕 (giantchen_AT_gmail) Blog.csdn.net/Solstice t.sina.com.cn/giantchen 这是《Muduo 网络编程示例》系列的第五篇文章。 Muduo 全系列文章列表&#xff1a; http://blog.csdn.net/Solstice/category/779646.…

camerax 自动聚焦_Android Camera-CameraView和CameraX使用

前面两篇介绍了Camera1和Camera2的使用&#xff0c;发现Camera API从1到2的变化非常大&#xff0c;Camera2的复杂度提升了不少&#xff0c;官方为了让我们更容易使用Camera&#xff0c;出了个一个官方的库cameraview。不过这个库已经Deprecated&#xff0c;官方建议使用Jetpack…

ECMA5

<!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><title>Document</title> </head> <body><script>var arr [2,4,1,6,7,10,9];// 遍历数组// arr.forEach(function(value,index,arra…

数据结构表达式求值(C++)

实验原理 代码思路 使用两个顺序栈编程&#xff0c;一个顺序栈用来存储操作数&#xff0c;一个顺序栈用来存储操作符 关键为表达式求值函数的编写 对该函数&#xff0c;循环结束条件为栈顶元素为’#’,且读入的字符也为’#’ 循环体 &#xff08;1&#xff09;判断是否为操作…