信号的采样与插值重建(包含matlab)
实验任务
双频信号的采样和插值重建:对双频信号进行采样(符合奈奎斯特低通采样定理),绘制波形和频谱;编制sinc函数插值函数,对采样进行恢复,绘制波形。
实验原理
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
信号的采样与插值重建(包含matlab)相关推荐
- matlab信号内插,【 MATLAB 】MATLAB 实现模拟信号采样后的重建(三)应用三次样条函数spline实现内插...
前三篇博文讲了三种方法进行内插重建信号: 这篇文章使用三次样条函数spline来实现内插重建,并分析重建误差. 采用的案例依然是上篇博文中的案例: 模拟信号: 对该信号使用两种不同的采样频率采样. a ...
- matlab 信号插零,【 MATLAB 】MATLAB 实现模拟信号采样后的重建(二)零阶保持(ZOH)...
这篇博文我们使用零阶保持器(ZOH)来重建信号,采用的案例依然是上篇博文中的案例: 模拟信号: 对该信号使用两种不同的采样频率采样. a. 在 fs = 5000 对信号进行采样 b. 在 fs = ...
- 【 MATLAB 】MATLAB 实现模拟信号采样后的重建(三)应用三次样条函数spline实现内插
前三篇博文讲了三种方法进行内插重建信号: sinc函数内插重建 零阶保持(ZOH)内插重建 一阶保持内插(FOH)重建 这篇文章使用三次样条函数spline来实现内插重建,并分析重建误差. 采用的案例 ...
- 【 MATLAB 】MATLAB 实现模拟信号采样后的重建(三)一阶保持(FOH)内插
上篇博文采用了零阶保持(ZOH)的方式进行了重构:[ MATLAB ]MATLAB 实现模拟信号采样后的重建(二)零阶保持(ZOH) 这篇博文我们使用一阶保持(FOH)内插来重建信号,采用的案例依然是 ...
- 【 MATLAB 】MATLAB 实现模拟信号采样后的重建(二)零阶保持(ZOH)
上篇博文采样sinc函数内插的方式实现了模拟信号的重建:[ MATLAB ]MATLAB 实现模拟信号采样后的重建(一) 这篇博文我们使用零阶保持器(ZOH)来重建信号,采用的案例依然是上篇博文中的案 ...
- 【 MATLAB 】MATLAB 实现模拟信号采样后的重建(一)
为了让MATLAB数字信号处理的相关博文能够得到一个梳理,我开通了一个专栏:数字信号处理的MATLAB实现 模拟信号经过采样后得到x(n),从x(n)中重建模拟信号在数学上可用公式来描述: 式中, 是 ...
- 信号的采样与恢复matlab实验报告,实验七 连续信号的采样与恢复
一.实验目的 通过MATLAB仿真验证抽样定理,进一步加深对抽样定理的理解. 二.实验原理 1. 连续信号的采样 对某一连续时间信号f(t)的采样原理图为: 由图可知, ,其中,单位冲激采样信号 的表 ...
- matlab计算信号得频谱,用MATLAB分析离散信号的频谱与信号的采样
<用MATLAB分析离散信号的频谱与信号的采样>由会员分享,可在线阅读,更多相关<用MATLAB分析离散信号的频谱与信号的采样(7页珍藏版)>请在人人文库网上搜索. 1.实验六 ...
- 连续时间信号的采样与恢复 Matlab仿真
目录 前言 采样 为什么需要采样 采样是什么 奈奎斯特采样定理 采样恢复 采样定理Matlab仿真(含代码) 前言 记录在学数字信号处理过程中一些个人理解与感悟,本文主要讲的是采样定理在Matlab实 ...
最新文章
- SmartRules让MindManager的交互图变得更加智能
- Spring AOP原理及拦截器
- 【通信原理】【实验】实验二:数字基带传输实验(GZHU)
- 【蓝桥杯】基础练习 十六进制转十进制
- spring 自定义日志_Spring和Hibernate的自定义审核日志
- 浏览器渲染页面的过程
- antd 日期选择框如何提交_基于Ant Design的Modal组件来实现一个可拖拽的React模态框...
- hdu1789----贪心+回溯
- boost::test
- 在Access数据库中使用SQL查询分析器
- 摩尔庄园不同服务器账号互通吗,摩尔庄园手游服务器互通吗?摩尔庄园手游不同区可以一起玩吗?...
- 用JavaScript实现,用户输入自己的汽油编号,然后输入自己加多少升,就弹出价格
- 2021哔哩哔哩1024程序员节日第一弹:算法与安全
- SSM开发笔记-尚硅谷-佟刚-Spring4.0.0
- 惹事之后,亚马逊AI助手Alexa将增新功能,真的能比自己更懂你?
- 《流浪地球》中的逃生气囊球和马斯克有什么关系?
- 餐中餐(5)Lucene--存储文件加载(Part 1: Segments加载)
- G4560成功上Bigsur系统
- H5 适配 iPhone X
- 【Metashape精品教程5】影像质量评估