matlab抽取与插值

什么是抽取?

我们假设一个数字信号 x ( n ) , n = 1 , 2 , . . . , N x(n),n=1,2,...,N x(n),n=1,2,...,N共有 N N N个点,抽取就是每个几个点抽1个点,比如2倍抽取,那么抽取后的信号为 y ( n ) , y ( 1 ) = x ( 1 ) , y ( 2 ) = x ( 3 ) , y ( 3 ) = x ( 5 ) , . . . , y ( N / 2 ) = x ( N − 1 ) y(n),y(1)=x(1),y(2)=x(3),y(3)=x(5),...,y(N/2)=x(N-1) y(n),y(1)=x(1),y(2)=x(3),y(3)=x(5),...,y(N/2)=x(N−1),可以看到抽取之后信号的长度变短了。

什么是插值?

插值和抽取相反,插值之后信号会变得越来越密集,注意绝对时间长度是不变的。例如数字信号 x ( n ) , n = 1 , 2 , . . . , N x(n),n=1,2,...,N x(n),n=1,2,...,N共有 N N N个点,2倍插值之后,就会变成 2 N 2N 2N个点,但是假设 N N N个点时间向量对应的是 t 0 < t < t 1 t_0<t<t_1 t0<t<t1,那么插值之后对应的时间还是它,只不过采样率提高了。

抽取仿真

假设, x ( n ) = c o s ( 2 π f n ) , f = 1 H z x(n)=cos(2\pi fn),f=1Hz x(n)=cos(2πfn),f=1Hz,信号的采样率 f s = 16 H z f_s=16Hz fs=16Hz,每个周期内可以采到16个点,我们进行 M = 2 M=2 M=2倍的抽取。

抽取可以用matlab的decimate函数。

matlab 复制代码
%% 信号的抽取和插值
f = 1;  %信号频率
fs = 16;%采样频率
N = 33; %采样点数
n = 0:1:N-1;
xn = cos(2*pi*f*n*(1/fs));
figure; 
stem(n,xn); %绘制原信号
title('x(n)');
xlabel('n');
ylabel('amplitude');

M = 2; %抽取的倍数
xn1 = decimate(xn,M); %抽取之后的信号
stem(xn1)%注意对上m和xn1
title('进行M=4倍抽取后的x(n)图像')
xlabel('n')
ylabel('幅度')

绘制的信号如下:

图1 原始信号


图2 抽取后的信号

插值可以用interp函数

matlab 复制代码
L = 3; %L倍插值
xn2 = interp(xn,L);  %以L倍采样率进行插值
stem(xn2)
title('进行L=3倍插值后的x(n)图像')
xlabel('n')
ylabel('幅度')

3倍插值后的信号

图3 3倍插值后的信号

指定插值方法的interp1函数使用方法

上面的插值,只需要指定倍数就可以了,matlab还有一种方法能够能够做插值,并且可以指定插值方法

interp1函数是Matlab中用于数据插值的工具,它能够根据给定的一组数据点(x、y),以及所要得到的插值点(xq),计算出插值点对应的函数值(vq)。interp1函数支持的插值方式包括:线性、样条、分段常数、分段线性和分段三次等。其中线性插值是最简单的插值方法,其他插值方法可以满足更高的精度要求。

复制代码
vq = interp1(x, y, xq) 
vq = interp1(x, y, xq, method) 
vq = interp1(x, y, xq, method, extrapval)

interp1函数的语法格式包含3~5个输入参数和一个输出参数,具体参数说明如下:

  • x:一维向量,表示数据点的自变量;
  • y:一维向量,表示数据点的因变量;
  • xq:一维向量,表示插值点的自变量;
  • method: 字符串,表示插值方式,可选值为'linear'、'nearest'、'next'、'previous'、'pchip'、'cubic'和'spline';linter是线性插值,spline是样条插值,pchip是分段多项式插值;
  • extrapval:数字或字符串,表示在插值点xq超出x的范围时用于外推的值,默认为NaN
  • 一维向量,表示插值点对应的函数值
相关推荐
看海天一色听风起雨落4 分钟前
Python学习之装饰器
开发语言·python·学习
Want5959 分钟前
C/C++圣诞树①
c语言·开发语言·c++
老赵的博客21 分钟前
c++ 杂记
开发语言·c++
jimmy.hua23 分钟前
[C++刷怪笼]:set/map--优质且易操作的容器
开发语言·c++
w2sfot1 小时前
Passing Arguments as an Object in JavaScript
开发语言·javascript·ecmascript
郝学胜-神的一滴2 小时前
避免使用非const全局变量:C++中的最佳实践 (C++ Core Guidelines)
开发语言·c++·程序人生
搞一搞汽车电子2 小时前
S32K3平台eMIOS 应用说明
开发语言·驱动开发·笔记·单片机·嵌入式硬件·汽车
星马梦缘2 小时前
Matlab机器人工具箱使用2 DH建模与加载模型
人工智能·matlab·机器人·仿真·dh参数法·改进dh参数法
总有刁民想爱朕ha3 小时前
车牌模拟生成器:Python3.8+Opencv代码实现与商业应用前景(C#、python 开发包SDK)
开发语言·python·数据挖掘
小菜全3 小时前
uniapp新增页面及跳转配置方法
开发语言·前端·javascript·vue.js·前端框架