数字信号处理信号抽样与插值MATLAB实现——实例

因为之前MATLAB没有好好学，而且DSP也没有好好学，做起作业来十分困难，为了好好写作业，我又重新研究了DSP和MATLAB，学了一下信号的抽样和插值
先看一看x(n)原本的图像

代码如下

%定义基本参数
f=1;%原信号频率
fs=16;%抽样信号频率
N=52;%采样点数
n=0:N-1;

%画原序列图像
xn=cos(2*pi*n*(f/fs));
stem(n,xn)
title('原x(n)图像')
xlabel('n')
ylabel('幅度')

题目中说了f/fs=16，也就是说T=16，在图像中也吻合
那么我们要将其进行4倍的抽样，该怎么做

实现代码如下

%定义基本参数
f=1;%原信号频率
fs=16;%抽样信号频率
N=52;%采样点数
n=0:N-1;

xn=cos(2*pi*n*(f/fs));

%画M=4倍抽样的x(n)
M=4;
m=0:N/M-1;%52/4=13 13-1=12 0:12有13个点
xn1=decimate(xn,M);%对序列进行抽样
subplot(3,1,2)
stem(m,xn1)%注意对上m和xn1
title('进行M=4倍抽取后的x(n)图像')
xlabel('n')
ylabel('幅度')

抽样函数

其中decimate()函数的意思就是抽样

decimate(xn,M)%表示对xn进行M倍的抽样，并不能进行插值

注意看图像的横坐标，之前的T=16，现在的T=4

下面对信号进行插值

%定义基本参数
f=1;%原信号频率
fs=16;%抽样信号频率
N=52;%采样点数
n=0:N-1;

xn=cos(2*pi*n*(f/fs));

%画L=3倍插值的x(n)
L=3;
xn2=interp(xn,L);
p=0:N*L-1;
stem(p,xn2)
title('进行L=3倍插值后的x(n)图像')
xlabel('n')
ylabel('幅度')

插值函数

interp(xn,L);%表示对xn进行L的插值，L必须为正整数

同样，观察横坐标，既然L=3，那么T就变成了48

在画图时要注意stem()里面的参数，常常会因为抽样点和信号对不上画不出来

完整代码

clear
close all
clc

%定义基本参数
f=1;%原信号频率
fs=16;%抽样信号频率
N=52;%采样点数
n=0:N-1;

%画原序列图像
xn=cos(2*pi*n*(f/fs));
%subplot(3,1,1)
stem(n,xn)
title('原x(n)图像')
xlabel('n')
ylabel('幅度')

%画M=4倍抽样的x(n)
M=4;
m=0:N/M-1;%52/4=13 13-1=12 0:12有13个点
xn1=decimate(xn,M);%对序列进行抽样
subplot(3,1,2)
stem(m,xn1)%注意对上m和xn1
title('进行M=4倍抽取后的x(n)图像')
xlabel('n')
ylabel('幅度')

%画L=3倍插值的x(n)
L=3;
xn2=interp(xn,L);
p=0:N*L-1;
subplot(3,1,3)
stem(p,xn2)
title('进行L=3倍插值后的x(n)图像')
xlabel('n')
ylabel('幅度')