数字信号处理信号抽样与插值MATLAB实现——实例
因为之前MATLAB没有好好学,而且DSP也没有好好学,做起作业来十分困难,为了好好写作业,我又重新研究了DSP和MATLAB,学了一下信号的抽样和插值
先看一看x(n)原本的图像
代码如下
1 2 3 4 5 6 7 8 9 10 11 12
| %定义基本参数 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倍的抽样,该怎么做
实现代码如下
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
| %定义基本参数 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()函数的意思就是抽样
1
| decimate(xn,M)%表示对xn进行M倍的抽样,并不能进行插值
|
注意看图像的横坐标,之前的T=16,现在的T=4
下面对信号进行插值
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
| %定义基本参数 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('幅度')
|
插值函数
1
| interp(xn,L);%表示对xn进行L的插值,L必须为正整数
|
同样,观察横坐标,既然L=3,那么T就变成了48
在画图时要注意stem()里面的参数,常常会因为抽样点和信号对不上画不出来
完整代码
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37
| 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('幅度')
|
![](https://mfjblog-pic.oss-cn-beijing.aliyuncs.com/img/b94f7fca5ad34b91a6521bf20cd336d4.png)