多抽样率数字信号处理（摘抄）

最新推荐文章于 2023-08-21 11:22:10 发布

2021乐乐

最新推荐文章于 2023-08-21 11:22:10 发布

阅读量1.5k

点赞数

分类专栏：信号处理文章标签：抽样和取样

本文链接：https://blog.csdn.net/weixin_42528089/article/details/106473565

版权

信号处理专栏收录该内容

15 篇文章 8 订阅

订阅专栏

10.0：序言

在应用中需要将定抽样率的信号转化为具有不同抽样率的等价信号。例如音频信号常常需要在三种不同速率中进行抽样率的转换。多抽样率数字信号处理系统使用了下抽样和上抽样这两种基本的抽样率器件。以及各种常规期间，如加法器、乘法器及延时器。

离散时间系统如果在系统多个部分存在不同的抽样率，就称为多抽样率系统。在变抽样率中，基本的抽样速率转换器一定是和低通数字滤波器。

10.1：基本抽样率转换器件

10.1.1 时域特性，

一个上抽样因子为L的上抽样器，L是正整数，生成的输出序列 $x_{u}[n]$ 是输入序列x[n]的L倍。

$\huge x_{u}[n]=\left\{\begin{matrix} x[n/L], n=0,\pm L,\pm 2L,... \\ 0, others \end{matrix}\right.$

上抽样例子:

clc,clear all,close all
N =input('输入长度=');% 50
L=input('上抽样因子=');% 3
fo=input('输入信号频率='); % 0.12
n=0:N-1;
x=sin(2*pi*fo*n);
%%% y
y=zeros(1,L*length(x));
y([1:L:length(y)])=x; %%% 扩展的关键
figure;
subplot(311); stem(n,x);  title(' the input of  array');  xlabel('时间序号n'),ylabel('振幅');
subplot(312); stem(n,y(1:length(x)));  title('上抽样后的序列');  xlabel('时间序号n'),ylabel('振幅');
subplot(313);plot(y);  title(' the original of y');

生成的图片：

从图中可以看到，将原始信号x，有50扩展到150，L=3 。但是幅值未发生变化.因此，上抽样器，也称为抽样率扩展器，或者扩展器。

下抽样器： XI下抽样因子为M 的抽样器，M为正整数，生成的输出序列y[n]的抽样率为输入序列x[n]的1/M 倍。即是每一次抽样保留输入序列中的第M个样本，而除去M-1个样本。

$\large y[n]=x[nM]$

下抽样例子

clc,clear all,close all;
N =input('输入长度=');% 50
M=input('下抽样因子=');% 3
fo=input('输入信号频率='); % 0.042
n=0:N-1;
m=0:N*M-1;
x=sin(2*pi*fo*m);
%
y=x([1:M:length(x)]);%%% 下采样的关键
subplot(311); stem(n,x(1:N));  title(' the input of  array');  xlabel('时间序号n'),ylabel('振幅');
subplot(312); stem(n,y);  title('下抽样后的序列');  xlabel('时间序号n'),ylabel('振幅');
subplot(313);plot(y);  title(' the original of y');

生成的图片

在抽样率转换过程中的不同阶段有着各自对应的时间序号。

10.1.2 频域特性

上抽样后的频谱：有频谱镜像的效果, 需要使用低通滤波器处理。

clc,clear all,close all;
freq=[0 0.45 0.5 1];
mag=[0 1 0 0];
x=fir2(99,freq,mag);
%%%
figure(1);
[Xz,w]=freqz(x,1,512);
plot(w/pi,abs(Xz));grid 
xlabel('\omega/\pi');ylabel('幅度');
title(' 输入频谱');


L=input('上抽样因子=');% 3
y=zeros(1,L*length(x));
y([1:L:length(y)])=x;
[Yz,w]=freqz(y,1,512);
figure(2);
plot(w/pi,abs(Yz));grid on;
xlabel('\omega/\pi');ylabel('幅度');
title(' 输出频谱');

下抽样后的频谱，容易产生频谱混叠现象.

clc,clear all,close all;
freq=[0 0.42 0.48 1];
mag=[0 1 0 0];
x=fir2(101,freq,mag);
%%%
figure(1);
[Xz,w]=freqz(x,1,512);
plot(w/pi,abs(Xz));grid 
xlabel('\omega/\pi');ylabel('幅度');
title(' 输入频谱');

M=input('下抽样因子=');% 3
y=x([1:M:length(x)]);
[Yz,w]=freqz(y,1,512);
figure(2);
plot(w/pi,abs(Yz));grid on;
xlabel('\omega/\pi');ylabel('幅度');
title(' 输出频谱');

3 倍下抽样后的频谱, 占据了整个频谱范围.