因为之前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('幅度')